Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Distinctive microbial community and genome structure in coastal seawater from a human-made port and nearby offshore island in northern Taiwan facing the Northwestern Pacific Ocean

  • Chi-Yu Shih,

    Roles Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Bachelor Degree Program in Marine Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Taiwan Ocean Genome Center, National Taiwan Ocean University, Keelung, Taiwan

  • Shiow-Yi Chen,

    Roles Data curation, Investigation, Methodology

    Affiliations Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan

  • Chun-Ru Hsu,

    Roles Conceptualization, Data curation

    Affiliations Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan

  • Ching-Hsiang Chin,

    Roles Data curation, Software, Visualization

    Affiliation Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan

  • Wei-Chih Chiu,

    Roles Methodology, Visualization

    Affiliation Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan

  • Mei-Hung Chang,

    Roles Methodology, Visualization

    Affiliation Biotools Co., Ltd., New Taipei City, Taiwan

  • Lee-Kuo Kang,

    Roles Methodology, Resources, Supervision

    Affiliation Bachelor Degree Program in Marine Biotechnology, National Taiwan Ocean University, Keelung, Taiwan

  • Cing-Han Yang,

    Roles Software, Visualization

    Affiliations Department of Computer Science and Information Engineering, National Taipei University of Technology, Taipei, Taiwan, Department of Computer Science and Engineering, National Taiwan Ocean University, Keelung, Taiwan

  • Tun-Wen Pai,

    Roles Formal analysis, Methodology, Software

    Affiliations Department of Computer Science and Information Engineering, National Taipei University of Technology, Taipei, Taiwan, Department of Computer Science and Engineering, National Taiwan Ocean University, Keelung, Taiwan

  • Chin-Hwa Hu,

    Roles Conceptualization, Resources, Supervision

    Affiliations Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan

  • Pang-Hung Hsu ,

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

    phsu@ntou.edu.tw (PHH); wstzou@ntou.edu.tw (WST)

    Affiliations Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan

  • Wen-Shyong Tzou

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing

    phsu@ntou.edu.tw (PHH); wstzou@ntou.edu.tw (WST)

    Affiliations Taiwan Ocean Genome Center, National Taiwan Ocean University, Keelung, Taiwan, Departent of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung, Taiwan, Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung, Taiwan

Abstract

Pollution in human-made fishing ports caused by petroleum from boats, dead fish, toxic chemicals, and effluent poses a challenge to the organisms in seawater. To decipher the impact of pollution on the microbiome, we collected surface water from a fishing port and a nearby offshore island in northern Taiwan facing the Northwestern Pacific Ocean. By employing 16S rRNA gene amplicon sequencing and whole-genome shotgun sequencing, we discovered that Rhodobacteraceae, Vibrionaceae, and Oceanospirillaceae emerged as the dominant species in the fishing port, where we found many genes harboring the functions of antibiotic resistance (ansamycin, nitroimidazole, and aminocoumarin), metal tolerance (copper, chromium, iron and multimetal), virulence factors (chemotaxis, flagella, T3SS1), carbohydrate metabolism (biofilm formation and remodeling of bacterial cell walls), nitrogen metabolism (denitrification, N2 fixation, and ammonium assimilation), and ABC transporters (phosphate, lipopolysaccharide, and branched-chain amino acids). The dominant bacteria at the nearby offshore island (Alteromonadaceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Rhodobacteraceae) were partly similar to those in the South China Sea and the East China Sea. Furthermore, we inferred that the microbial community network of the cooccurrence of dominant bacteria on the offshore island was connected to dominant bacteria in the fishing port by mutual exclusion. By examining the assembled microbial genomes collected from the coastal seawater of the fishing port, we revealed four genomic islands containing large gene-containing sequences, including phage integrase, DNA invertase, restriction enzyme, DNA gyrase inhibitor, and antitoxin HigA-1. In this study, we provided clues for the possibility of genomic islands as the units of horizontal transfer and as the tools of microbes for facilitating adaptation in a human-made port environment.

Introduction

As interference and pollution by humans have continued to deteriorate our environment, many studies have focused on the pollution of rivers and seawater by investigating the microbiome to explore the composition of microorganisms in the water as a means for inferring their effects on public health and the environment. For example, researchers discovered the presence of multidrug resistance genes in bacteria in the Ganges River in India [1]. Metagenomics analysis demonstrated the disease-causing bacterial pathogens Burkholderia, Shigella, and Salmonella in the downstream metropolitan sewage discharge area of São Pedro, Brazil [2]. Researchers in China also found metal tolerance genes and drug resistance genes in bacteria collected from mangrove sediments on Hainan Island [3] and the Pearl River Delta Dawan District [4], China. Some research has also focused on changes in microbiota and potential bioremediation in coastal areas and commercial harbors [59].

Fishing ports are unique in that the sources of pollution are numerous. The pollution could be derived from domestic sewage, waste oil and sewage discharged from ships, and wastewater discharged from the surrounding fish market and restaurants. Some fishing ports in Taiwan have both economic and recreational functions, and the environment of the fishing ports recently raised concern due to their pollution [10]. The Badouzi Fishing Port in Keelung, Taiwan, is located at the northern tip of Taiwan, facing the East China Sea and the Northwestern Pacific Ocean. In contrast to the Badouzi Fishing Port, the offshore Heping Island is a recreational park without any fishing or economic activity. This research aims to investigate the microbiomes in the seawater of these two adjacent locations. By employing both 16S rRNA gene amplicon and whole-genome shotgun sequencing, we discovered the differences in the microbiomes in these two locations and identified the potential disease-causing bacteria with antibiotic resistance genes and metal tolerance genes. Network analysis also uncovered the coexistence and mutual-exclusion relationship in the microbial community as the response of the microbial consortia to pollution. Microbial genomes from the fishing port also revealed genomic islands potentially capable of horizontal transfer. This study is the first one for revealing microbial community and its associated genome structure in coastal seawater in northern Taiwan facing the Northwestern Pacific Ocean.

Materials and methods

Sample collection and DNA extraction

Heping Island and Badouzi Fishing Port are located at the northeastern point of Taiwan, facing the East China Sea in the Northwestern Pacific Ocean (Fig 1A, S1 and S2 Figs in S1 File). The distance between Heping Island and Badouzi Fishing Port is approximately 4 kilometers. Coastal surface seawater was collected on Heping Island (25.16194 N, 121.76296 E) and Badouzi Fishing Port (25.14120 N, 121.79287 E) in triplicate at different locations on the mornings of July 8 and July 9, 2021, respectively (N = 3 for sampling in Heping Island and Badouzi Fishing Port). Seawater was passed through a 1.2 μm filter first to remove debris and grainy particles and then further passed through a 0.47 μm filter to collect the prokaryotic cells. The filters were stored at -20°C until further processing. Total genomic DNA from samples was extracted using the column-based method (QIAamp PowerFecal DNA Kit, Qiagen). All the samples were processed separately.

thumbnail
Fig 1. Map of the seawater sampling sites.

(A) The sampling sites (in green pentagon) are in northern Taiwan facing the Northwestern Pacific Ocean. Sampling points in a fishing port (Badouzi Fishing Port in orange circles) and nearby offshore island (Heping Island in blue boxes) are shown (map was made from open source https://github.com/Wesely/Taiwan-Python-Map). (B) Principal component analysis (PCA) based on OTU level for the microbiota composition of all samples (PC1 = 76.9%, PC2 = 16.8%). Heping Island (HPisland) and Badouzi Fishing Port (BDZfishP) are indicated in blue boxes and orange circles, respectively.

https://doi.org/10.1371/journal.pone.0284022.g001

16S rRNA gene amplicon sequencing and analysis

For 16S rRNA gene sequencing, the V3-V4 region was amplified by a specific primer set (341F: 5’-CCTACGGGNGGCWGCAG-3’, 806R: 5’-GACTACHVGGGTAT CTAATCC-3’) [11], and the sequencing library was prepared according to the 16S Metagenomic Sequencing Library Preparation procedure (Illumina). The library was sequenced on an Illumina MiSeq platform, and paired 300‐bp reads were generated. Demultiplexing was carried out based on barcode identification. Raw amplicon paired-end reads were assembled using FLASH (v1.2.11) [12]. Low-quality reads (Q <20) were discarded in the QIIME (v1.9.1) pipeline [13]. A read was truncated if the quality score of three consecutive bases was < Q20; the resulting read was retained in the dataset only if it was at least 75% of the original length (using split_libraries_fastq.py script in QIIME [14]). UCHIMES was employed to filter sequence chimeras [15, 16]. The effective tags were clustered at 97% sequence identity to generate operational taxonomic units (OTUs) by using the UPARSE [17] function in the USEARCH (v7.0.1090) pipeline [18]. For each representative sequence, the RDP classifier (v2.11) algorithm [19] was employed to annotate taxonomy classification based on the information retrieved from SILVA (an on-line resource for quality checked and aligned ribosomal RNA sequence data, https://www.arb-silva.de/) (138.1) [20, 21] with an 80% minimum confidence threshold. To analyze the sequence similarities among different OTUs, multiple sequence alignment was conducted by using PyNAST software (v1.2) [22] against the core-set dataset in the Silva database.

To normalize the variations in sequence depth across samples, OTU abundance information was rarefied to the minimum sequence depth using the QIIME script (single_rarefaction.py). Subsequent analysis of alpha and beta diversities was performed using the normalized data. The relative abundance and evenness accounting for diversity were evaluated by the Shannon and Simpson indices using the QIIME pipeline. A rarefaction curve was constructed by a random selection of a certain amount of sequencing data of each sample to represent the number of observed species [23]. For beta diversity, T-SNE was performed using the R package Rtsne [24].

Whole-genome shotgun sequencing and analysis

Total genomic DNA from samples was extracted using the column-based method (QIAamp PowerFecal DNA Kit, Qiagen). DNA degradation and potential contamination were monitored on 1% agarose gels. DNA quantification was checked using a Qubit® dsDNA Assay Kit in a Qubit® 4.0 Fluorometer (Life Technologies, CA, USA). A total amount of 1 μg DNA per sample was used as input material for the library preparations. Sequencing libraries were generated using Illumina Nextera DNA Flex Library Prep (Illumina, USA) following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample. The library quality was assessed on the Qubit 4.0 Fluorometer and a Qsep100TM system. Subsequently, the library was sequenced on an Illumina NovaSeq platform, and paired 150 bp reads were generated.

The original data obtained by high-throughput sequencing (Illumina NovaSeq 6000 platform) were transformed into raw sequenced reads by CASAVA base calling and stored in FASTQ format. The obtained raw paired-end reads were filtered by Trimmomatic (v0.38) [25] to discard low-quality reads and trim adaptor sequences and to eliminate poor-quality bases. The obtained high-quality data (clean reads) were used for subsequent analyses. Bowtie2 (v2.3.4.1) [26] was used to filter out the contaminating sequences. The filtered reads were then assembled using MEGAHIT (v1.1.3) [27]. Prodigal (v2.6.3) [28] was used to predict the open reading frames (ORFs) from the assembled contigs with lengths > 500 bps. ORFs were filtered to remove those shorter than 100 bps. A nonredundant gene catalog was constructed using CD-HIT (v4.7) with 95% identity [29, 30]. Filtered reads from each sample were mapped to the initial gene catalog using BWA (v0.7.17-r1188). SAMtools (v1.8) was used to create BAM files [31], and the jgi_summarize_bam_contig_depths script from the MetaBAT2 pipeline [32] was run on all BAM files to calculate the coverage for each sample. MetaBAT2 (v2.12.1) was used to carry out metagenomic binning with a minimum contig length threshold of 2500 bps for generating metagenome-assembled genomes (MAG). The genome percentage completeness and contamination of all bins were assessed using CheckM (v1.1.3) [33]. The obtained unigenes were blasted against the NCBI Refseq [34] database using DIAMOND (v0.9.22.123) [35]. The taxonomy assignment was determined by using the lowest common ancestor (LCA) algorithm. The abundances of each taxonomic group were calculated by summing the abundance of genes annotated to a feature. The constructed nonredundant gene catalog was annotated by several functional databases to assign their function. Gene annotation was conducted by aligning sequences against the Refseq (a database that provides comprehensive, integrated, non-redundant, well-annotated set of sequences, including genomic DNA, transcripts, and proteins, https://www.ncbi.nlm.nih.gov/refseq/) [34], eggNOG (a database of hierarchical, functionally and phylogenetically annotated orthology resource, http://eggnog5.embl.de/) [36], KOfam (a database and web server to assign KEGG Orthologs to protein sequences by homology search against a database of profile hidden Markov models, https://www.genome.jp/tools/kofamkoala/) [37], VFDB (a database of virulence factors of bacterial pathogens, http://www.mgc.ac.cn/VFs/) [38], CARD (a database of resistance genes, their products and associated phenotypes, https://card.mcmaster.ca/) [39], NCyc (a database for metagenomic profiling of nitrogen cycling genes, https://github.com/qichao1984/NCyc) [40], dbCAN2 (a database for carbohydrate-active enzymes, https://bcb.unl.edu/dbCAN2/) [41], and BacMet (a database of bacterial genes that are experimentally confirmed to confer resistance to metals and/or antibacterial biocides, http://bacmet.biomedicine.gu.se/) [42] databases using DIAMOND (v0.9.22.123), HMMER (v3.2.1) and other database-specific annotators. Several methods were developed to discover genomic islands in microbial genomes [43], and we employed IslandViewer 4 [44] to search for genomic islands in three assembled genomes. The phylogenomics of MAGs were based on anvi’o workflow [45] and the tree-drawing tool iTOL [46].

Coding sequences from genomic sequences of 8 strains of Phaeobacter italicus (RefSeq assembly accession numbers: GCF_019801785.1, GCF_001258055.1, GCF_001404195.1, GCF_900113345.1, GCF_020171225.1, GCF_017743945.1, GCF_019801155.1, and GCF_001251095.1) were downloaded from NCBI. The DNA sequences of genomic islands from the assembled genome (bin) were searched by BLASTN against coding sequences of the 8 strains stated above with the e-value set to 1e-5.

We analyzed the community network structure by analyzing OTU results (abundance value) using CoNet [47]. CoNet deduced the Spearman correlation among all pairs and ran bootstrapping 1000 times. Correlation coefficients larger than 0.8 and p values smaller than 0.05 were retained and shown by Cytoscape [48].

Results

Taxonomic analysis using 16S rRNA gene amplicon sequencing reads

We investigated the microbiome composition by sequencing the 16S rRNA gene V3-V4 regions of bacteria. There were 707 and 712 OTUs collected from 224,703 and 202,954 tags for seawater from Heping Island and Badouzi Fishing Port, respectively (S1 Table in S2 File). The alpha diversity and evenness index indicated that the microbial community in coastal seawater in the fishing port (5.11 of Shannon index and 0.58 of Pielou’s evenness) had higher phylogenetic diversity and evenness than that in the offshore island (4.16 of Shannon index and 0.47 of Pielou’s evenness) (S3-S11 Figs in S1 File, S2 Table in S2 File). PCA showed that the microbial communities in fishing ports and offshore islands were well clustered separately; however, OTUs in fishing ports were more diverse in the second PCA dimension (Fig 1B). A Venn diagram showed that 43% of OTUs were shared between the microbiota in fishing ports and offshore islands (Fig 2A).

thumbnail
Fig 2. Composition of OTU and LDA analyses.

(A) Venn diagrams showing the shared and unique OTUs of the microbiota in the costal seawater of Heping Island and Badouzi Fishing Port. (B) Linear discriminant analysis effect size (LEfSe) cladogram of the microbial composition in seawater of Heping Island and Badouzi Fishing Port. Green and red nodes/shades specify taxa significantly higher in relative abundance in Heping Island and Badouzi Fishing Port, respectively. The taxa abundance is proportional to the circle size. The threshold of the linear discriminant analysis (LDA) score was 4.0. (C) Histogram of the LDA scores showing differential abundances of microbial taxa between Heping Island and Badouzi Fishing Port.

https://doi.org/10.1371/journal.pone.0284022.g002

We employed SILVA to conduct the taxonomic analysis of OTUs. We found significant differences in the abundance of characteristic microbiota in both sampling sites as determined via LEfSe (Linear discriminant analysis Effect Size, an analysis that can identify genomic biomarkers characterizing statistical differences among biological groups analysis) (Fig 2B and 2C). At the phylum level, Proteobacteria and Bacteroidetes comprised the highest percentages of OTUs in the seawater from Heping Island (50.5% and 43.5%, respectively) and Badouzi Fishing Port (65.3% and 28.7%, respectively) at the phylum level. The sums of the OTU percentages of Proteobacteria and Bacteroidetes from the two sampling locations were the same (94.0%), indicating that the increase or decrease in these two phyla could account for the major difference in the bacterial community in seawater from Heping Island and Badouzi Fishing Port. Compared with the community structure of Heping Island, the decrease in Bacteroidetes was on par with the increase in Proteobacteria at Badouzi Fishing Port (for Proteobacteria, 65.3%-50.5% = 14.8%; for Bacteroidetes, 28.7–43.5% = -14.8%).

At the family level (S3 Table in S2 File), at Heping Island, Cryomorphaceae (34.9%), Litoricolaceae (14.2%), Rhodobacteraceae (11.5%), Alteromonadaceae (10.0%), Flavobacteriaceae (7.3%), Clade_I (2.2%), and Actinomarinaceae (2.1%) were the top families of the bacterial community; at Badouzi Fishing Port, Vibrionaceae (22.6%), Rhodobacteraceae (21.4%), Flavobacteriaceae (17.2%), Nitrincolaceae (8.2%), Cryomorphaceae (8.0%), and Pseudoalteromonadaceae (2.8%) were the top families of the bacterial community (Fig 3A).

thumbnail
Fig 3. Microbial composition (in relative abundance) in Badouzi Fishing Port and Heping Island.

The distribution of major contributing families (top ten) based on the analysis of (A) 16S rRNA gene amplicon sequencing and (B) whole-genome shotgun sequencing.

https://doi.org/10.1371/journal.pone.0284022.g003

These data indicated that in comparison with Heping Island in terms of the percentage of community structure, the changes in the populations of Cryomorphaceae (Bacteroidetes, 26.9% down), Litoricolaceae (Proteobacteria, 14.1% down), and Alteromonadaceae (Proteobacteria, 7.3% down) were commensurate with the changes in the populations of Vibrionaceae (Proteobacteria, 22.3% up), Rhodobacteraceae (Proteobacteria, 9.9% up), Flavobacteriaceae (Bacteroidetes, 9.9% up) and Nitrincolaceae (Proteobacteria, 8.1% up) at Badouzi Fishing Port.

Taxonomic analysis using whole-genome shotgun sequencing

The numbers of paired-end raw reads sequenced were 173 and 198 million for seawater from Heping Island and Badouzi Fishing Port, respectively (S4 Table in S2 File). A total of 2,352,957 contigs (with 947,562 contigs longer than 500 nt), 1,359,362 ORFs, and 409,712 unique genes were generated. Finally, the numbers of unique genes annotated by Refseq and eggNOG was 225,025 and 301,223, respectively, for seawater from Heping Island. For seawater from Badouzi Fishing Port, the numbers of unique genes annotated by Refseq and eggNOG was 230,953 and 305,820, respectively.

At the phylum level, Proteobacteria and Bacteroidetes were still the two most abundant phyla. For Heping Island, 55.0% and 38.3% of microbes detected belong to Proteobacteria and Bacteroidetes, respectively. For Badouzi Fishing Port, 80.3% and 18.0% of microbes detected belonged to Proteobacteria and Bacteroidetes, respectively.

At the family level, at Heping Island, the unclassified family was the most abundant (16.1%), followed by Flavobacteriaceae (13.8%), Rhodobacteraceae (13.3%), Litoricolaceae (12.0%), Cryomorphaceae (11.9%), Alteromonadaceae (8.3%), Halieaceae (3.6%), and Rhodospirillaceae (2.2%) (Fig 3B, S5 Table in S2 File).

At Badouzi Fishing Port, Rhodobacteraceae was the most abundant (31.8%), followed by Vibrionaceae (15.2%), Flavobacteriaceae (12.3%), Oceanospirillaceae (12.2%), unclassified (7.6%), Alteromonadaceae (3.8%), Halieaceae (3.5%), Rhodospirillaceae (2.7%) and Pseudoalteromonadaceae (2.5%) (Fig 3B, S5 Table in S2 File).

If we compared the percentages of the dominant taxa at Heping Island with those at Badouzi Fishing Port, we observed the decreases (in percentage) of Litoricolaceae (11.9% decrease), Cryomorphaceae (10.1% decrease), and Alteromonadaceae (4.6% decrease) and the increases of Rhodobacteraceae (18.5% increase), Vibrionaceae (14.6% increase), and Oceanospirillaceae (11.4% increase).

General functional analysis using Kyoto Encyclopedia of Genes and Genomes (KEGG)

We investigated the functional categories by employing KEGG analysis using ORFs as the input sequences and found 334 pathways associated with ORFs. Pathways common to both collections with high abundance included (top 20, Heping Island and Badouzi Fishing Port, respectively) ko01100 Metabolic pathways (29.5%, 23.2%), ko01110 Biosynthesis of secondary metabolites (13.3%, 9.8%), ko01120 Microbial metabolism in diverse environments (7.5%, 6.3%), ko01240 Biosynthesis of cofactors (6.4%, 4.6%), ko01230 Biosynthesis of amino acids (5%, 3.7%), ko01200 Carbon metabolism (4.5%, 3.3%), ko02010 ABC transporters (2.5%, 2.8%), ko02020 Two-component system (2%, 2.6%), ko03010 Ribosome (3.7%, 2%), ko00230 Purine metabolism (2.7%, 2%), ko00240 Pyrimidine metabolism (2.4%, 1.7%), ko02024 Quorum sensing (1.7%, 1.7%), ko00190 Oxidative phosphorylation (1.9%, 1.4%), ko00630 Glyoxylate and dicarboxylate metabolism (1.5%, 1.3%), ko00260 Glycine, serine and threonine metabolism (1.6%, 1.3%), ko00620 Pyruvate metabolism (1.5%, 1.2%), ko01212 Fatty acid metabolism (1.6%, 1.2%), ko00280 Valine, leucine and isoleucine degradation (1.2%, 1.1%), ko01250 Biosynthesis of nucleotide sugars (1.7%, 1%), and ko00860 Porphyrin and chlorophyll metabolism (1.1%, 1%).

We also found several pathways unique to sampling sites in Badouzi Fishing Port: ko01053 Biosynthesis of siderophore group nonribosomal peptides (9.6 folds compared with that in Heping Island), ko05110 Vibrio cholerae infection (33.3 folds), ko05150 Staphylococcus aureus infection (112.3 folds) and ko05135 Yersinia infection (210.2 folds). Furthermore, we also found several pathways unique to sampling sites in Heping Island: ko00944 Flavone and flavonol biosynthesis (10.4 folds compared with that in Badouzi Fishing Port) and ko04144 Endocytosis (6.3 folds).

Nitrogen metabolism

We employed the NCyc database to understand nitrogen metabolism. nirK, nirS, nosZ, and nmo were more abundant at both Heping Island and Badouzi Fishing Port than other genes (Fig 4A, S6 Table in S2 File). (t-test p = 2.2e-16, 1.4e-10, 1.4e-10, and 3.2e-07 for nirK, nirS, nosZ, and nmo, respectively).

thumbnail
Fig 4. Heatmaps showing the relative abundances of genes of metabolisms.

(A) nitrogen cycling metabolism, (B) carbohydrate metabolism, and (C) ABC transporters. Samples were collected from three sites in Badouzi Fishing Port (BDZfishP1, BDZfishP2, BDZfishP3) and Heping Island (HPisland1, Hpisland2, Hpisland3).

https://doi.org/10.1371/journal.pone.0284022.g004

nirK, nirS, nosZ, and nmo were all responsible for denitrification. Denitrification consists of four enzymatic steps, starting from nitrate and producing the intermediates nitrite, nitric oxide, and nitrous oxide. nirK and nirS catalyzed the reduction of nitrite. nosZ catalyzed the reduction of nitrous oxide. nmo is an enzyme (nitronate monooxygenase) that catalyzes the oxidative denitrification of alkyl nitronates by using oxygen (O2).

We found that most of the nirK gene sequences were affiliated with Flavobacteriaceae, Litoricolaceae, and Alteromonadaceae at Heping Island and with Vibrionaceae, Rhodobacteraceae, and Oceanospirillaceae at Badouzi Fishing Port (S12 Fig in S1 File, S6 Table in S2 File). Most of the nirS gene sequences were affiliated with Cryomorphaceae, Litoricolaceae, Alteromonadaceae at Heping Island and with Rhodobacteraceae at Badouzi Fishing Port (S13 Fig in S1 File). Most of the nosZ and nmo gene sequences were affiliated with Litoricolaceae at Heping Island and Rhodobacteraceae at Badouzi Fishing Port (S14 and S15 Figs in S1 File, S6 Table in S2 File).

We also found genes responsible for ammonium assimilation (glnA, gdh) at both Heping Island (Cryomorphaceae) and Badouzi Fishing Port (Rhodobacteraceae) (S16-S18 Figs in S1 File, S6 Table in S2 File). A gene responsible for N2 fixation (nifH) was also found in bacteria collected from both Heping Island (Flavobacteriaceae) and Badouzi Fishing Port (Oceanospirillaceae) (S19 Fig in S1 File, S6 Table in S2 File). Other genes responsible for N2 fixation were also found in bacteria collected from both Heping Island (nifW from Alteromonadaceae, Litoricolaceae and Rhodobacteraceae) and Badouzi Fishing Port (nifKD from Vibrionaceae).

Moreover, narJ at Heping Island (Flavobacteriaceae), a chaperone protein for the assembly of nitrate reductase (narG, a nitrate reductase that can convert nitrate to nitrite) and molybdenum cofactor, was also observed.

Carbohydrate metabolism

We employed the dbCAN2 database to search for genes responsible for carbohydrate metabolism. In 452 families of bacteria found in the sampling sites, 183 families of them were demonstrated to have genes responsible for degrading, modifying, or creating glycosidic bonds (40.5%). Both GT2 and GT4 were highly abundant at Heping Island and Badouzi Fishing Port (Fig 4B, S7 Table in S2 File). GT stands for the family of glycosyltransferases that catalyze the transfer of sugar moieties from activated donor molecules to specific acceptor molecules in the biosynthesis of disaccharides, oligosaccharides, and polysaccharides. GT2 includes cellulose synthase, chitin synthase, and hyaluronan synthase, which are involved in biofilm and capsule components in bacteria. GT4 includes sucrose synthase and sucrose-phosphate synthase, which are involved in sucrose and sucrose 6-phosphate, respectively [49].

We also found that GH13 gene sequences were abundant at Heping Island. GH stands for a family of glycoside hydrolases that can hydrolyze the glycosidic bonds between two or more carbohydrates or between a carbohydrate and a noncarbohydrate moiety. GH13 includes alpha-amylase, pullulanase, sucrose phosphorylase, and glucosidase, which can degrade phytoplankton-derived glucan. GH23 was abundant at Badouzi fishing port. GH23 is a peptidoglycan lyase important for the remodeling of bacterial peptidoglycan (PG) in the cell wall of bacteria and its pathogenicity.

We found that most of the GT2 gene sequences were affiliated with Rhodobacteraceae at Heping Island and Badouzi Fishing Port (S20 Fig in S1 File, S7 Table in S2 File). GT4 gene sequences were affiliated with Cryomorphaceae and Alteromonadaceae at Heping Island and with Rhodobacteraceae and Vibrionaceae at Badouzi Fishing Port (S21 Fig in S1 File, S7 Table in S2 File). At Heping Island, GH13 gene sequences were affiliated with Cryomorphaceae and Flavobacteriaceae (S22 Fig in S1 File, S7 Table in S2 File). At Badouzi Fishing Port, GH23 gene sequences were affiliated with Rhodobacteraceae and Vibrionaceae (S23 Fig in S1 File, S7 Table in S2 File).

ABC transporters and phosphotransferase

To understand ABC transporter metabolism in the collected bacteria, we employed the ggnog database and found 3,129 genes. Gene sequences for pstB (ATP-binding protein for phosphate transport), ccmC (export of heme to the periplasm for the biogenesis of c-type cytochromes), pstC (phosphate transport system permease protein), and lptB (lipopolysaccharide export) were abundant at both Heping Island and Badouzi Fishing Port (Fig 4C, S8 Table in S2 File). Gene sequences for bmpA (sugar transport system substrate-binding protein), ybiT (ATP-binding protein), ttg2C (phospholipid or cholesterol transport), ccmB (permease involved in cytochrome c biogenesis), troA (manganese/zinc/iron transport system substrate-binding protein), and devA (ABC exporter, ATP-binding subunit) were abundant at Heping Island. msbA (ATP-dependent lipid A-core flippase), livM (permeases for branched-chain amino acid transport) and livG (permeases for branched-chain amino acid transport) were abundant at Badouzi Fishing Port.

At Badouzi Fishing Port and Heping Island, we found that most of the pstB sequences were affiliated with Rhodobacteraceae. pstC sequences were affiliated with Litoricolaceae at Heping Island and Vibrionaceae at Badouzi Fishing Port. ccmC sequences were affiliated with Litoricolaceae and Alteromonadaceae at Heping Island and Oceanospirillaceae at Badouzi Fishing Port. lptB sequences were affiliated with Litoricolaceae and Cryomorphaceae at Heping Island and with Rhodobacteraceae and Vibrionaceae at Badouzi Fishing Port.

At Heping Island, bmpA, ybiT, and troA were affiliated with Litoricolaceae, Cryomorphaceae, and Thermoactinomycetaceae, respectively. At Badouzi Fishing Port, msbA and livM and livG were affiliated with Rhodobacteraceae, and livG was affiliated with Gammaproteobacteria.

Phosphotransferase system (PTS) was utilized by bacteria to uptake sugar as the source energy. We also employed eggNOG database and found 53 genes encoding phosphotransferase. Gene sequences for pstN (Phosphotransferase system mannitol fructose-specific IIA domain) were abundant at both Heping Island and Badouzi Fishing Port. Gene sequences for thrB (Phosphotransferase enzyme family) was abundant at Heping Island. cmtB, frwC (Phosphotransferase system mannitol fructose-specific IIA domain), celC (Phosphotransferase system cellobiose-specific component IIA), crr (Phosphotransferase system IIA components), chpT (Histidine phosphotransferase C-terminal domain), and fruA (Phosphotransferase system fructose-specific component IIB) were abundant at Badouzi Fishing Port.

At Badouzi Fishing Port, we found that most of the celC, cmtB, crr, fruA, frwC, and pstN sequences were affiliated with Vibrionaceae. chpT and pstN sequences were affiliated with Rhodobacteraceae. pstN sequences was also affiliated with Oceanospirillaceae.

At Heping Island, pstN and thrB were affiliated with Litoricolaceae. pstN sequences was also affiliated with Rhodobacteraceae.

Antibiotic resistance genes

We found 68 genes conferring antibiotic resistance in the collected bacterial samples by employing the CARD database. We found 452 families of bacteria in the sampling sites and 41 families of them have AR genes of presence (9.1%). While 33.8% of 68 AR genes were found in the sea water of Heping Island, 97.1% of AR genes were found in the sea water of Badouzi Fishing Port.

Regarding resistance mechanisms, “antibiotic efflux” was the most prevalent (35 genes, 51.5%). Thirteen genes (19.1%) functioned through “antibiotic inactivation”, 12 genes (17.6%) through antibiotic target alteration, 7 genes (10.3%) through antibiotic target replacement, and one gene through antibiotic target protection. We also categorized potential bacteria-resistant antibiotics into 20 antibiotic categories. Genes conferring resistance to diaminopyrimidines and macrolides were abundant in the seawater of Heping Island (Fig 5A, S9 Table in S2 File). Genes responsible for the resistance of diaminopyrimidines include dfrA26, dfrA17, dfrA15, dfrA3, and dfrA20 (dihydrofolate reductase, antibiotic target replacement). Genes responsible for the resistance to macrolides include macB (ABC transporter), oleC (ABC transporter), mexW (RND-type membrane protein of the efflux complex MexVW-OprM), CpxR (activation of expression of the RND efflux pump MexAB-OprM), LpeB (subunit of the LpeAB efflux pump), and abeS (efflux pump of the SMR family). Conversely, genes conferring resistance to ansamycin, nitroimidazole, and aminocoumarin were abundant in the seawater of Badouzi Fishing Port (Fig 5A, S9 Table in S2 File). The gene responsible for the resistance of ansamycin was rpoB2 (beta-subunit of RNA polymerase, antibiotic target alteration, or replacement). The gene involved in the resistance of nitroimidazole was msbA (multidrug resistance transporter, antibiotic efflux). The gene responsible for the resistance of aminocoumarin was novA (type III ABC transporter, antibiotic efflux).

thumbnail
Fig 5. Heatmaps showing the relative abundances of genes conferring resistance to antibiotics and metals and encoding virulence factors.

(A) Resistance to 20 categories of antibiotics and (B) Metal tolerance to 13 categories of metals and (C) Virulence factors. Samples were collected from three sites in Badouzi Fishing Port (BDZfishP1, BDZfishP2, BDZfishP3) and Heping Island (HPisland1, Hpisland2, Hpisland3).

https://doi.org/10.1371/journal.pone.0284022.g005

At Badouzi Fishing Port, we found that most of the gene sequences conferring ansamycin resistance were affiliated with Vibrionaceae (S31 Fig in S1 File, S9 Table in S2 File). Gene sequences conferring nitroimidazole resistance were affiliated with Rhodobacteraceae, Shewanellaceae, and Oceanospirillaceae (S32 Fig in S1 File, S9 Table in S2 File). Gene sequences conferring aminocoumarin resistance were affiliated with Rhodobacteraceae.

At Heping Island, we found that most of the gene sequences conferring diaminopyrimidine or macrolide resistance were affiliated with Litoricolaceae and Alteromonadaceae (S33 Fig in S1 File, S9 Table in S2 File).

More specifically (S9 Table in S2 File), at Heping Island, Litoricola lipolytica (in the family of Litoricolaceae within Gammaproteobacteria) harnessed five antibiotic resistance genes (CRP, Pseudomonas, dfrA3, Nocardia, msbA). At Badouzi Fishing Port, Vibrio within Gammaproteobacteria contained 13 antibiotic resistance genes (ugd, MexF, abeS, CARB-18, catB8, CARB-19, CRP, qacH, hmrM, FosC2, MCR-3, mexW and catB9). Rhodobacteraceae within Alphaproteobacteria harbored four antibiotic resistance genes (msbA, novA, smeR, and ugd). Oleibacter marinus and Marinobacterium sp. AK27 (in the family Oceanospirillaceae within Gammaproteobacteria) contained five (msbA, CRP, dfrA26, acrB, Streptomyces) and three antibiotic resistance genes (Nocardia, MexF, and ugd), respectively.

Metal tolerance genes

We found 85 genes in the BacMet database involved in metal resistance. In 452 families of bacteria found in the sampling sites, 44 families of them have metal tolerance genes of presence (9.7%). While 51.8% of 85 metal tolerance genes were found in the sea water of Heping Island, 98.8% of metal tolerance genes were found in the sea water of Badouzi Fishing Port.

There were thirteen types of metal resistance (gene number): multimetal (35 genes, channel, enzyme, efflux, transcription regulator, membrane transporter, binding protein), arsenic (9 genes, transcription regulator, chaperone, enzyme, binding protein), aluminum (one gene, enzyme), chromium (3 genes, efflux and enzyme), copper (11 genes, enzyme and transcription regulator), iron (8 genes, enzyme, binding protein, transcription regulator, and membrane protein), lead (one gene, enzyme), mercury (4 genes, membrane transporter and transcription regulator), nickel (3 genes), selenium (2 genes, enzyme), silver (2 genes, efflux and enzyme), tellurium (3 genes), and zinc (5 genes, enzyme, efflux, membrane transporter, binding protein, and transcription regulator) (S10 Table in S2 File).

At Heping Island, we found abundant genes resistant to copper, arsenic, and multimetal (Fig 5B, S10 Table in S2 File). Most of the gene sequences conferring multimetal and arsenic resistance were affiliated with Alteromonadaceae (S34, S35 Figs in S1 File). Gene sequences conferring copper resistance were affiliated with Cryomorphaceae (S36 Fig in S1 File, S10 Table in S2 File).

In the seawater of Badouzi Fishing Port, we found abundant genes conferring resistance to copper, chromium, iron, and multimetal compounds (Fig 5B, S10 Table in S2 File). Most of the gene sequences conferring multimetal resistance were affiliated with Rhodobacteraceae and Vibrionaceae (S34 Fig in S1 File, S10 Table in S2 File). Gene sequences conferring copper resistance were affiliated with Rhodobacteraceae (S36 Fig in S1 File, S10 Table in S2 File). Gene sequences conferring iron resistance were affiliated with Vibrionaceae and Oceanospirillaceae (S37 Fig in S1 File, S10 Table in S2 File). Gene sequences conferring chromium resistance were affiliated with Rhodobacteraceae, Vibrionaceae, and Oceanospirillaceae (S38 Fig in S1 File, S10 Table in S2 File).

More specifically (S10 Table in S2 File), at Heping Island, Litoricola lipolytica (in the family Litoricolaceae within Gammaproteobacteria) harbored six metal-tolerance genes. Both Phaeocystidibacter luteus (in the family Cryomorphaceae within Flavobacteriales) and Alteromonadaceae within Gammaproteobacteria contained one metal tolerance gene. At Badouzi Fishing Port, Vibrio within Gammaproteobacteria contained 24 metal-tolerance genes. Phaeobacter italicus (in the family Rhodobacteraceae within Alphaproteobacteria) harbored four metal-tolerance genes. Marinobacterium stanieri and Oleibacter marinus (in the family Oceanospirillaceae species within Gammaproteobacteria) harbored nine and two metal-tolerance genes, respectively.

Virulence factors

We searched the virulence factor database (VFDB) for virulence factors of the bacteria under investigation. We found 369 genes belonging to 122 groups. In 452 families of bacteria found in the sampling sites, 102 families of them have virulence genes of presence (22.6%). While 45.3% of 369 virulence genes were found in the sea water of Heping Island, 98.1% of virulence genes were found in the sea water of Badouzi Fishing Port.

At the group level, we discovered that LOS (CVF494, 367 genes) was abundant at both Heping Island and Badouzi Fishing Port; ND (AI144, 278 genes), VF0408 (53 genes) and VF0273 (244 genes) were abundant at Badouzi Fishing Port (Fig 5C, S11 Table in S2 File).

Genes of the LOS (lipooligosaccharide) (CVF494) group bear similarity to those from Haemophilus influenzae Rd KW20 and encode proteins involved in the structuring and biosynthesis of lipid A-containing complex glycolipids in the outer membranes and associated pathogenicity [5052]. Genes of the ND (AI144) group bear similarity to those from Aeromonas hydrophila subsp. hydrophila ATCC 7966 and encode proteins responsible for chemotaxis and flagellar synthesis and function. They were demonstrated to be involved in cirrhosis in humans and the invasion and survival of Aeromonas hydrophila in Anguilla japonica macrophages [53, 54]. Genes of the VF0273 group bear similarity to those from Pseudomonas aeruginosa PAO1 and encode proteins responsible for swimming motility of bacteria relying on flagellar activity and involved in biofilm formation and pathogenic adaptations [5557]. Genes of the VF0408 group bear similarity to those from Vibrio parahaemolyticus RIMD 2210633 and encode proteins for the T3SS1 system. T3SS1 is a type three secretion system (a needle-like structure) and causes the cytotoxicity of host cells, involving the induction of autophagy, cell rounding, and cell lysis [58].

We also found that most of the gene sequences of LOS (CVF494) were affiliated with Bacteroidetes and Alteromonadaceae at Heping Island and with Vibrionaceae and Oceanospirillaceae at Badouzi Fishing Port (S39 Fig in S1 File, S11 Table in S2 File). At Badouzi Fishing Port, ND (AI144) and T3SS1 (VF0408) were affiliated with Vibrionaceae, and Flagella (VF0273) was affiliated with Oceanospirillaceae (S40 and S41 Figs in S1 File, S11 Table in S2 File).

MAGs (metagenome-assembled genomes) in the coastal seawater and the characterization of metabolic pathways

We assembled microbial genomes of high quality and abundance collected from coastal surface seawater of Heping Island and Badouzi Fishing Port (completeness larger than 80%, contamination less than 5%), yielding 19 bins. These genomes were affiliated with two phyla of bacteria (Bacteroidota and Proteobacteria) and three major classes (Bacteroidia, Alphaproteobacteria and Gammaproteobacteria) (Fig 6A). The largest genome of Bins in size contains 4.5MBases (genus Vibrio). Genus UBA3478 (Bin.013, Flavobacteriaceae), UBA10364 (Bin.150, Schleiferiaceae) and UBA8309 (Bin.146, Puniceispirillaceae) were abundant in Heping Island. Genus Glaciecola (Bin.022, Alteromonadaceae), HIMB11 (Bin.007, Rhodobacteraceae), Thalassobius (Bin.053, Rhodobacteraceae), Phaeobacter (Bin.095, Rhodobacteraceae), Vibrio (Bin.040, Vibrionaceae) and OM182 (Bin.086, Pseudohongiellaceae) were abundant in seawater of Badouzi Fishing Port (Fig 6B).

thumbnail
Fig 6. 19 metagenome assembled genomes (MAGs) discovered in the coast seawater of Heping Island and Badouzi Fishing Port.

(A) Phylogenetic analysis of MAGs. Color codes indicate the taxonomic assignment of “class”. (B) Heatmap showing the abundance of 19 MAGs. Samples were collected from three sites in Badouzi Fishing Port (BDZfishP1, BDZfishP2, BDZfishP3) and Heping Island (HPisland1, HPisland2, HPisland3). Accompanied table includes statistics of MAG including total length (bp), N50 (bp) and GC (%).

https://doi.org/10.1371/journal.pone.0284022.g006

We also mapped the metabolic pathway inferred from MAGs (Figs 7 and 8). Overrepresented gene functions in seawater of both sampling sites were metabolic pathways, biosynthesis of secondary metabolites, microbial metabolism in diverse environments, biosynthesis of cofactors, carbon metabolism, biosynthesis of amino acids, ABC transporters, fatty acid metabolism, purine metabolism and two-component system. On the other hand, through clustering of functions, we also found several groups of abundant gene functions associated with specific MAGs. For example, in Heping Island (Fig 7), we discovered that quorum sensing, glycine, serine/threonine metabolism, oxidative phosphorylation, porphyrin/chlorophyll metabolism, and glyoxylate/dicarboxylate metabolism were associated with genus HIMB11 (Bin.007, Rhodobacteraceae), Thalassobius (Bin.053, Rhodobacteraceae), Phaeobacter (Bin.095, Rhodobacteraceae) and UBA8309 (Bin.146, Puniceispirillaceae). Benzoate degradation, tryptophan metabolism, butanoate metabolism and fatty acid degradation were associated with genus SZUA_381 (Bin.115, Crocinitomicaceae), UBA974 (Bin.151, UA16) and Saprospiraceae (Bin.130).

thumbnail
Fig 7. Heatmaps showing the relative abundance of genes with metabolic functions associated with 19 MAGS discovered in the coast seawater of Heping Island.

(A) Top ten metabolic functions of high abundance. (B) Other metabolic functions. Color codes indicate the taxonomic assignment of “class”.

https://doi.org/10.1371/journal.pone.0284022.g007

thumbnail
Fig 8. Heatmaps showing the relative abundance of genes with metabolic functions associated with 19 MAGS discovered in the coast seawater of Badouzi Fishing Port.

(A) Top ten metabolic functions of high abundance. (B) Other metabolic functions. Color codes indicate the taxonomic assignment of “class”.

https://doi.org/10.1371/journal.pone.0284022.g008

Clustering of functions associated with specific MAGs were also found in Badouzi Fishing Port (Fig 8). For example, oxidative phosphorylation, glycine, serine/threonine metabolism, glyoxylate/dicarboxylate metabolism and porphyrin/chlorophyll metabolism were associated with genus HIMB11 (Bin.007, Rhodobacteraceae), Thalassobius (Bin.053, Rhodobacteraceae), Phaeobacter (Bin.095, Rhodobacteraceae), Roseicyclus (Bin.106, Rhodobacteraceae), Marinobacterium (Bin.078, Nitrincolaceae) and UBA8309 (Bin.146, Puniceispirillaceae). Pyrimidine metabolism, pyruvate metabolism, carbon fixation pathways, aminoacyl-tRNA biosynthesis, biosynthesis of nucleotide sugars and homologous recombination were associated with genus Vibrio (Bin.040, Vibrionaceae), Polaribacter (Bin.062, Flavobacteriaceae), Flavobacterium (Bin.076, Flavobacteriaceae), MAG_120531 (Bin.077, Flavobacteriaceae), MAG_120531 (Bin.127, Flavobacteriaceae) and Glaciecola (Bin.022, Alteromonadaceae).

Genomic islands (GIs) of the abundant bacteria Rhodobacteraceae in the coastal seawater of the fishing port

We assembled microbial genomes of high quality and abundance from the bacteria collected from Badouzi Fishing Port and found one affiliated with Rhodobacteraceae that resolved into the species Phaeobacter italicus, which fit the criteria (completeness 98.54%, contamination 0%, S12 Table in S2 File). Interestingly, while we compared the genome sequence of P. italicus from Baoudozi fishing port with that of P. italicus in NCBI, we found four genomic islands (GIs) of P. italicus at Badouzi Fishing Port that were not detected in the publicly available genomes in NCBI. Lower GC% contents than the neighboring regions of the genome, tRNA genes located on one end of the genomic island, and two direct repeats in the border of the genomic islands all mark the features of genomic islands. Within these four genomic islands, we discovered genes encoding phage tyrosine integrase, prophage DNA invertase, DNA gyrase inhibitor, restriction enzyme, antitoxin HigA-1, and many proteins of unknown function (S13 Table in S2 File). Of these four GIs, two can map to a single contig of publicly available genomes of P. italicus in NCBI (Figs 9 and 10). Within these two genomic islands, all genes but one have the same transcription direction, suggesting an operon structure. The protein-coding genes (e.g. antitoxin HigA-1, metallopeptidase, SNARE associated Golgi protein, Diguanylate cyclase, and CbiX protein for Cobalamin (Vitamin B12) biosynthesis, S13 Table in S2 File) within GIs would confer the advantage for the bacteria to adapt in the environment. This discovery implies that the horizontal transfer of genomic islands could facilitate the adaptation of P. italicus in the seawater environment of the fishing port.

thumbnail
Fig 9. GC contents and direct repeat sequences on the tRNA structures of two genomic islands in the genome of Phaeobacter italicus discovered in the costal seawater of Badouzi Fishing Port.

(A, B) GC content as a percentage, employing a scanning window of 3000 bps. (C, D) tRNA secondary structures and the direct repeat sequences on one side of the genomic islands (perfect repeat nucleotides are colored in red circles, and imperfect repeat nucleotides are colored in orange). The tRNAs in C and D belong to the genomic islands shown in A and B, respectively. (A) GI_2 (B) GI_4 (see S13 Table in S2 File).

https://doi.org/10.1371/journal.pone.0284022.g009

thumbnail
Fig 10. Genomic islands in the genome of Phaeobacter italicus discovered in the costal seawater of Badouzi Fishing Port.

Gene clusters in the genomic island are shown, and sequence comparisons of the genomic island (upper) with the corresponding genome of Phaeobacter italicus available in NCBI (lower) are shown in (A) GI_2 (B) GI_4 (see S13 Table in S2 File). tRNA genes are colored in gray. “DR” (in red circles) indicates the direct repeat on the two flanking sides of the genomic islands.

https://doi.org/10.1371/journal.pone.0284022.g010

The community network structure of the microbiome of the coastal seawater at the fishing port and the nearby island

We generated the community network structure by analyzing OTU results from 16S rRNA gene V3V4 amplicon sequencing (Fig 11). We found that the community network was composed of one main cluster associated with a small cluster. The main cluster of the network can be roughly divided into two parts: one with mainly coexistence and the other with mutual exclusion. The network with mainly coexistence was composed of the bacteria AEGEAN_169_marine_group, Clade_I, Clade_II, Cyanobiaceae, Puniceicoccaceae, Porticoccaceae, Actinomarinaceae, SAR116_clade, and Kiritimatiellaceae; the other network with mutual exclusion consisted of Rhodobacteraceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Marinilabiliaceae. The small cluster consisted of Pseudohongiellaceae, Nitrincolaceae, Saccharospirillaceae, Crocinitomicaceae, Arcobacteraceae, Rhodocyclaceae, and Saccharimonadaceae, all with coexistence relationships among the constituent nodes while connecting with the main cluster (through Litoricolaceae) by mutual exclusion. Interestingly, Alteromonadaceae formed its coexistence relationship with Litoricolaceae and NS9_marine_group, the latter with Microtrichaceae, then Kiritimatiellaceae, and then Cryomorphaceae. However, Alteromonadaceae maintained its mutual exclusion relationships with Vibrionaceae, Marinilabiliaceae, and Staphylococcaceae.

thumbnail
Fig 11. Microbial community network in the costal seawater of Heping Island and Badouzi Fishing Port.

The color of the network edges denotes positive (blue) and negative (red) correlation, and the edge width is proportional to the correlation coefficient. Only network edges with an absolute value of the correlation coefficient larger than 0.8 and bootstrap p value smaller than 0.01 were retained. Family names of bacteria are designated.

https://doi.org/10.1371/journal.pone.0284022.g011

Discussion

The bacterial community structures in coastal seawater from a fishing port and a nearby island were different

We employed next-generation sequencing to decipher the microbiome compositions and infer the biochemical properties of bacteria in seawater at two nearby locations facing the Northwestern Pacifica Ocean and the East China Sea. After investigating the derived biochemical functions from gene-coding protein sequences based on the genome sequences, some features emerged (Fig 11, S14 Table in S2 File). At Heping Island, Alteromonadaceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Rhodobacteraceae were dominant. They harbored genes with the functions of denitrification, nitrate assimilation, phytoplankton-derived glucan degradation, phosphate transport, lipopolysaccharide transport, sugar transport, antibiotic resistance (diaminopyrimidines, macrolide), and metal tolerance (copper and multimetal). Conversely, at Badouzi Fishing Port, we found the emergence of much higher proportions of Oceanospirillaceae, Rhodobacteraceae, and Vibrionaceae. More gene sequences with the functions of nitrogen fixation, denitrification, nitrate assimilation, biofilm formation, cell wall remodeling, transport of branched-chain amino acids, transport of phosphate, transport of lipopolysaccharide, flagellar motility, antibiotic resistance (nitroimidazole, aminocoumarin, ansamycin), metal tolerance (multimetal, copper, chromium, iron), chemotaxis and flagella and T3SS1 were discovered.

thumbnail
Fig 12. A conceptual figure for major discoveries of the microbiome from the costal seawater of Heping Island and nearby Badouzi Fishing Port.

Dominant species in two sampling sites are demonstrated. Specific functions of genes encoding for the nitrogen metabolism, carbohydrate metabolism, antibiotics resistance, and virulence factors are summarized.

https://doi.org/10.1371/journal.pone.0284022.g012

Alphaproteobacteria, Gammaproteobacteria, and Flavobacteriia are heterotrophic bacteria and have been discovered to respond to phytoplankton blooms by actively utilizing the nutrients (such as dissolved organic matter) released by phytoplankton [59, 60]. Cryomorphaceae (within Flavobacteriia) can degrade the high molecular weight substrates of coastal phytoplankton and have been found in coastal surface seawater worldwide [61]. Cryomorphaceae was also found to play important roles in the nutrient dynamics in inshore reefs in the Great Barrier Reef [62]. Litoricolaceae (within Gammaproteobacteria) was discovered in surface seawater in the Yellow Sea, Korea [63]. A metabolomic study showed that Litoricola marina (in the family Litoricolaceae) feeds on nucleotides and nucleosides in the phytoplankton bloom area [64]. A phycosphere study in the seawater of east-south Asia (Xiamen, Hong Kong, and the South China Sea) also showed that Rhodobacteraceae (within Alphaproteobacteria) was abundantly associated with marine Synechococcus at all three sampling locations; however, Alteromonadaceae (within Gammaproteobacteria) was abundant in coastal water around Xiamen and Hong Kong, and Flavobacteriaceae (within Flavobacteriia) were abundant only in the South China Sea [65]. In our study, both 16S rRNA gene amplicon and whole-genome shotgun sequencing results indicated that Alteromonadaceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Rhodobacteraceae were the five most abundant families in the seawater of Heping Island. This implies that in July of the sampling time, phytoplankton existed in the seawater of Heping Island.

The community structure inferred from our study suggested that in the seawater of Heping Island, Alteromonadaceae, Litoricolaceae, and Cryomorphaceae formed a coexistence relationship with a group of bacteria, including AEGEAN_169_marine_group, Clade_I, Clade_II, Cyanobiaceae, Puniceicoccaceae, Porticoccaceae, Actinomarinaceae, SAR116_clade, and Kiritimatiellaceae, most of which were reported to be abundant in phytoplankton blooms [66, 67]. Conversely, Rhodobacteraceae and Vibrionaceae showed higher abundances at Badouzi Fishing Port and formed mutual-exclusion relationships with the abovementioned bacteria, which maintained a relationship of coexistence. Moreover, the isolated group comprising Pseudohongiellaceae, Nitrincolaceae, Saccharospirillaceae, Crocinitomicaceae, Arcobacteraceae, Rhodocyclaceae, and Saccharimonadaceae formed mutual-exclusion relationships with Litoricolaceae and showed a higher abundance at Badouzi Fishing Port than at Heping Island. Arcobacteraceae was associated with human feces in sewage [68] and was also found in sewer sediment [69]. Rhodocyclaceae was demonstrated to function in single-chamber microbial fuel cells as a denitrifier [70]. Ethylbenzene dehydrogenase, which conducts the anaerobic bacterial degradation of ethylbenzene and propylbenzene, was also from Rhodocyclaceae [71]. Saccharimonadaceae was found to be more abundant in colorectal cancer (CRC) model mice when chromium was administered [72]. Saccharimonadaceae is also an intestinal metabolizer and microbial degrader in earthworms [73, 74].

Potential impacts posed by the bacterial community in the coastal seawater of fishing ports

Rhodobacteraceae was shown to be able to utilize a different form of carbon source [75] and to have large amounts of ABC transporters [65]. Rhodobacteraceae was also the most abundant bacterium involved in biofilm formation in eastern Mediterranean coastal seawater [76] and the western Pacific Ocean [77]. Rhodobacteraceae was the abundant taxon in the coastal waters of the reefs, inlets, and wastewater outfalls of southeast Florida, USA [78], where Rhodobacteraceae showed significant correlations with oxygen level, salinity, pH, and nitrate concentration. In the southern coastal seawater of India, both Rhodobacteraceae and Vibrionaceae were involved in biofilm formation [79]. Investigation of the microbial community in the Persian Gulf also showed that after chronic exposure and oil spill events, Oceanospirillales (containing Oceanospirillaceae) was enriched in the high aliphatic content of the pollution, and Rhodobacterales (containing Rhodobacteraceae) was enriched in polyaromatic pollution, leading to the hypothesis of the division of labor for bioremediation [80]. Climate, environment, and human activity have been shown to correlate with the abundance of the Vibrio genus [81, 82]. The spatiotemporal distribution and abundance of Vibrio [83] in the marine aquaculture environment of Dongshan Bay in the southwest Taiwan Strait were investigated, and 28 species were found to vary in abundance between seasons. Vibrionaceae was reported to be associated with cage-cultured marine fishes, plants, algae, zooplankton, and animals (fish, oysters, etc.) [84, 85]. Vibrionaceae contains many pathogens and is found in the aquatic environment [86]. At Badouzi Fishing Port, where leaked oil and dead fish body pollution pose serious problems in the sea environment, we observed the abundant presence of Rhodobacteraceae, Vibrionaceae, and Oceanospirillaceae, with genes conferring biochemical functions that lead to adaptation to the environment and allow them to thrive. More virulence genes also appeared at Badouzi Fishing Port, which poses public health problems.

Antibiotic resistance seen in the microbiome poses a serious threat to the environment and humans. A recent study using shotgun metagenomic sequencing of the intestinal microbiome of deep-sea fish from the Atlantic Ocean revealed that fish in the deep sea almost do not have antibiotic resistance genes [87]. In contrast, metagenomic samples from the TARA Oceans project revealed a wide distribution of antibiotic resistance genes in microbes in the ocean [88]. A study in the coastal area of northwestern Sicily, Italy, showed that seawater with polyethylene pollution led to the accumulation of antibiotic resistance genes in the microbiome [89]. Another study in collected fish from the mainland and marine environments in China as well as from Chile and Nigeria showed that antibiotic resistance in marine fish shared high similarity and high abundance and was distinct from that observed in terrestrial animals [90]. In this study, we found more antibiotic resistance genes (nitroimidazole, aminocoumarin, and ansamycin) in the microbiome in the seawater of the Badouzi Fishing Port. The pattern of antibiotic resistance at Badouzi Fishing Port was different from that at Heping Island. This result suggested that the pollution at Badouzi Fishing Port posed an additional potential threat to the community in the fishing port. Metal tolerance has recently gained more attention due to the impact of metal pollution on microbial communities and ecological functions [91]. We also identified abundant metal tolerance genes at Badouzi Fishing Port. For example, Vibrio contains several multimetal-tolerant genes, of which the mechanisms of tolerance include channel (glpF), efflux (corC), enzyme (corB, dsbB, modB, modC, recG, ruvB), and membrane transporter (fecE, mntH/yfeP). Moreover, genes that render chromium resistance were apparent in Vibrionaceae collected in the sea of the Badouzi Fishing Port. Tolerance to chromium by Vibrio parahaemolyticus in China was reported before [92], and the tolerance concentration of Vibrio can be as high as 3200 μg/mL [93]. Antimicrobial resistance can be propagated through the food chain from antimicrobial-resistant bacteria to the contaminated food consumed by animals and humans [94]. Even the fishmeal containing antimicrobial-resistant genes were detected in mariculture sediment, posing another round of antimicrobial resistance propagation [95].

We found an increase in genes affiliated with denitrification and N2 fixation in the microbiome in seawater of the Badouzi Fishing Port. This result suggested that the input of a nitrogen source provided some advantage in the environment to certain bacteria. A laboratory study showed that abrupt increases in nitrate or ammonium in surface water from the central North Pacific Ocean stimulated the growth of Oceanospirillaceae and Rhodobacteraceae in incubations in the dark [96]. Interestingly, a study on the biocrust showed that one of the most abundant transcriptionally active N-transforming microorganisms in the investigated biocrusts was affiliated with Rhodobacteraceae [97].

Scrutiny of the contigs from the collection of whole-genome shotgun sequencing revealed abundant genes coding for the T3SS1 apparatus at the Badouzi Fishing Port. A protein sequence similarity search against the NCBI database indicated that four T3SS1s belonged to Vibrio harveyi (three contigs), Vibrio alginolyticus (one contig), and Vibrio diabolicus (one contig). Recently, a study on the whole genomes of 110 Vibrio species from coastal areas in China demonstrated that T3SS1 genes were universally detected in Vibrio parahaemolyticus, Vibrio alginolyticus, Vibrio harveyi, and Vibrio campbellii [98]. Our effort in binning the reads from whole-genome shotgun sequencing yielded a Vibrio sequence of sufficient completeness to the genus level. We are currently employing long-read sequencing technology (such as Oxford Nanopore sequencing or PacBio single-molecule real-time sequencing) [99101], aiming to assemble a complete metagenome-assembled bacterial genome and reach a thorough investigation of whole TSS1 clusters.

Through investigating the community network structure (Fig 11), both coexistence and mutual-exclusion networks were discovered. Dominated bacteria in the seawater of Heping Island were mainly composed of a coexistence network; however, the dominated bacteria in the seawater of Badouzi Fishing Port posed an “antagonistic” relation with the coexistence network, forming mutual-exclusion interaction. Ecological niche overlap and resource competition [102104] were demonstrated for causing the community structure change. In our study, we suggest that the dominant bacteria in the seawater of Heping island shared resources and, possibly, formed a metabolite exchange scenario in the natural ecological environment; nonetheless, in the seawater of Badouzi Fishing Port, the pollution and human interference caused large scale change of ecological environment, leading to the prospering opportunist bacteria without the need to cooperate with other bacteria.

Genomic islands as the unit of horizontal transfer among microbes in the seawater of fishing port

Genomic islands could serve as a means for horizontal transfer between bacteria, facilitating their adaptation and evolution [105, 106]. We found candidate genomic islands which have features such as lower GC% contents, tRNA genes located on one end of the genomic island, and two direct repeats in the border of the genomic islands. We also found genes encoding phage integrase and recombinase that could trigger the excision and insertion of genomic elements in an appropriate environment. Genomic islands could harbor pathogenesis genes of medical importance. We found several genomic islands in the assembled genomes from whole-genome shotgun sequencing, many of which contained genes of high interest. For example, in a genomic island of P. italicus (in the family Rhodobacteraceae), we found the antitoxin HigA-1. HigA-1 belongs to the HigBA system, a toxin-antitoxin system of type II toxin-antitoxin [107109]. HigB protein can endonucleotically cleave RNA in the ribosome 30S-RNA complex [110, 111]. The HigBA toxin-antitoxin system was suggested to play an important role in stress conditions and its adaptation to environmental challenges. HigA can bind HigB to inactivate its molecular function by binding the promoter of the higBA operon, which contains two palindromes [112, 113], requiring three water molecule-mediated amino acid-nucleotide interactions [112]. In our study, the presence of HigA in a genomic island could stem from horizontal transfer and be essential for the survival of P. italicus in the fishing port environment.

Conclusions

We discovered, in the northwestern Pacific Ocean, the dominant bacteria Oceanospirillaceae, Rhodobacteraceae, and Vibrionaceae at a fishing port and Alteromonadaceae, Cryomorphaceae, Flavobacteriaceae, Litoricolaceae, and Rhodobacteraceae at a nearby offshore island. Bacteria present at fishing ports develop multiple strategies to adapt to the heavily human-interfered environment, while bacteria at offshore islands are related to phytoplankton. Community structure reflected the change in consortia of bacteria and their relative relationships of both coexistence and mutual exclusion in the sea environment. Genomic islands deduced from the bacterial genomes of the fishing port indicate that horizontal transfer could take place in the seawater of the fishing port, and, using gene exchange, the challenge to and adaptation of bacteria in the polluted sea environment could pose threat to humans and associated activities.

Acknowledgments

We would like to express our gratitude to the Heping Island Park authority for approving the sampling of seawater during the study period.

References

  1. 1. Mittal P, Prasoodanan Pk V, Dhakan DB, Kumar S, Sharma VK. Metagenome of a polluted river reveals a reservoir of metabolic and antibiotic resistance genes. Environmental Microbiome. 2019;14(1):5. pmid:33902720
  2. 2. Medeiros JD, Cantao ME, Cesar DE, Nicolas MF, Diniz CG, Silva VL, et al. Comparative metagenome of a stream impacted by the urbanization phenomenon. Braz J Microbiol. 2016;47(4):835–45. Epub 2016/08/16. pmid:27522532; PubMed Central PMCID: PMC5052392.
  3. 3. Li Y, Zheng L, Zhang Y, Liu H, Jing H. Comparative metagenomics study reveals pollution induced changes of microbial genes in mangrove sediments. Sci Rep. 2019;9(1):5739. Epub 2019/04/07. pmid:30952929; PubMed Central PMCID: PMC6450915.
  4. 4. Chen J, McIlroy SE, Archana A, Baker DM, Panagiotou G. A pollution gradient contributes to the taxonomic, functional, and resistome diversity of microbial communities in marine sediments. Microbiome. 2019;7(1):104. Epub 2019/07/17. pmid:31307536; PubMed Central PMCID: PMC6632204.
  5. 5. Miettinen H, Bomberg M, Nyyssonen M, Reunamo A, Jorgensen KS, Vikman M. Oil degradation potential of microbial communities in water and sediment of Baltic Sea coastal area. PLoS One. 2019;14(7):e0218834. Epub 2019/07/03. pmid:31265451; PubMed Central PMCID: PMC6605675.
  6. 6. Sauret C, Tedetti M, Guigue C, Dumas C, Lami R, Pujo-Pay M, et al. Influence of PAHs among other coastal environmental variables on total and PAH-degrading bacterial communities. Environ Sci Pollut Res Int. 2016;23(5):4242–56. Epub 2015/07/01. pmid:26122564.
  7. 7. Jeanbille M, Gury J, Duran R, Tronczynski J, Agogue H, Ben Said O, et al. Response of Core Microbial Consortia to Chronic Hydrocarbon Contaminations in Coastal Sediment Habitats. Front Microbiol. 2016;7:1637. Epub 2016/10/30. pmid:27790213; PubMed Central PMCID: PMC5061854.
  8. 8. Neethu CS, Saravanakumar C, Purvaja R, Robin RS, Ramesh R. Oil-Spill Triggered Shift in Indigenous Microbial Structure and Functional Dynamics in Different Marine Environmental Matrices. Scientific Reports. 2019;9(1):1354. pmid:30718727
  9. 9. Nogales B, Bosch R. Microbial Communities in Hydrocarbon-Polluted Harbors and Marinas. In: McGenity T (eds) Microbial Communities Utilizing Hydrocarbons and Lipids: Members, Metagenomics and Ecophysiology Handbook of Hydrocarbon and Lipid Microbiology 2018:1–15.
  10. 10. Chen G-Y. A Study on the Integrated Strategies of Marine Pollution Control in the Fishing Ports: Taking Badouzih Fishing Port as the Case. MSc Thesis, Institute of Marine Resource Management, National Taiwan Ocean University. 2007. Available from: http://ethesys.lib.ntou.edu.tw/cdrfb3/record/#G0M94370002.
  11. 11. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1. Epub 20120828. pmid:22933715; PubMed Central PMCID: PMC3592464.
  12. 12. Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63. pmid:21903629
  13. 13. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 2010;7(5):335–6. pmid:20383131
  14. 14. Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10(1):57–9. Epub 2012/12/04. pmid:23202435; PubMed Central PMCID: PMC3531572.
  15. 15. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21(3):494–504. Epub 2011/01/08. pmid:21212162; PubMed Central PMCID: PMC3044863.
  16. 16. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. Epub 2011/06/28. pmid:21700674; PubMed Central PMCID: PMC3150044.
  17. 17. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8. Epub 2013/08/21. pmid:23955772.
  18. 18. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1. Epub 2010/08/17. pmid:20709691.
  19. 19. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7. Epub 2007/06/26. pmid:17586664; PubMed Central PMCID: PMC1950982.
  20. 20. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Research. 2012;41(D1):D590–D6. pmid:23193283
  21. 21. Yilmaz P, Parfrey LW, Yarza P, Gerken J, Pruesse E, Quast C, et al. The SILVA and "All-species Living Tree Project (LTP)" taxonomic frameworks. Nucleic Acids Res. 2014;42(Database issue):D643–8. Epub 2013/12/03. pmid:24293649; PubMed Central PMCID: PMC3965112.
  22. 22. Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26(2):266–7. Epub 2009/11/17. pmid:19914921; PubMed Central PMCID: PMC2804299.
  23. 23. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41. Epub 2009/10/06. pmid:19801464; PubMed Central PMCID: PMC2786419.
  24. 24. Van der Maaten L. Accelerating t-SNE using Tree-Based Algorithms. Journal of Machine Learning Research. 2014;15:25.
  25. 25. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. Epub 2014/04/04. pmid:24695404; PubMed Central PMCID: PMC4103590.
  26. 26. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. Epub 2012/03/06. pmid:22388286; PubMed Central PMCID: PMC3322381.
  27. 27. Li D, Luo R, Liu CM, Leung CM, Ting HF, Sadakane K, et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11. Epub 2016/03/26. pmid:27012178.
  28. 28. Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11(1):119. pmid:20211023
  29. 29. Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2. Epub 2012/10/13. pmid:23060610; PubMed Central PMCID: PMC3516142.
  30. 30. Li W, Fu L, Niu B, Wu S, Wooley J. Ultrafast clustering algorithms for metagenomic sequence analysis. Briefings in bioinformatics. 2012;13(6):656–68. Epub 2012/07/10. pmid:22772836; PubMed Central PMCID: PMC3504929.
  31. 31. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. Epub 2009/06/10. pmid:19505943; PubMed Central PMCID: PMC2723002.
  32. 32. Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165. Epub 2015/09/04. pmid:26336640; PubMed Central PMCID: PMC4556158.
  33. 33. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55. Epub 2015/05/16. pmid:25977477; PubMed Central PMCID: PMC4484387.
  34. 34. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Database issue):D61–5. Epub 2006/11/30. pmid:17130148; PubMed Central PMCID: PMC1716718.
  35. 35. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. Epub 2014/11/18. pmid:25402007.
  36. 36. Huerta-Cepas J, Szklarczyk D, Heller D, Hernandez-Plaza A, Forslund SK, Cook H, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47(D1):D309–D14. Epub 2018/11/13. pmid:30418610; PubMed Central PMCID: PMC6324079.
  37. 37. Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36(7):2251–2. Epub 2019/11/20. pmid:31742321; PubMed Central PMCID: PMC7141845.
  38. 38. Liu B, Zheng D, Jin Q, Chen L, Yang J. VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019;47(D1):D687–D92. Epub 2018/11/06. pmid:30395255; PubMed Central PMCID: PMC6324032.
  39. 39. Alcock BP, Raphenya AR, Lau TTY, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020;48(D1):D517–D25. Epub 2019/10/31. pmid:31665441; PubMed Central PMCID: PMC7145624.
  40. 40. Tu Q, Lin L, Cheng L, Deng Y, He Z. NCycDB: a curated integrative database for fast and accurate metagenomic profiling of nitrogen cycling genes. Bioinformatics. 2019;35(6):1040–8. Epub 2018/08/31. pmid:30165481.
  41. 41. Zhang H, Yohe T, Huang L, Entwistle S, Wu P, Yang Z, et al. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018;46(W1):W95–W101. Epub 2018/05/18. pmid:29771380; PubMed Central PMCID: PMC6031026.
  42. 42. Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DG. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014;42(Database issue):D737-43. Epub 2013/12/07. pmid:24304895; PubMed Central PMCID: PMC3965030.
  43. 43. Bertelli C, Tilley KE, Brinkman FSL. Microbial genomic island discovery, visualization and analysis. Briefings in bioinformatics. 2019;20(5):1685–98. Epub 2018/06/06. pmid:29868902; PubMed Central PMCID: PMC6917214.
  44. 44. Bertelli C, Laird MR, Williams KP, Simon Fraser University Research Computing G, Lau BY, Hoad G, et al. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017;45(W1):W30–W5. Epub 2017/05/05. pmid:28472413; PubMed Central PMCID: PMC5570257.
  45. 45. Eren AM, Esen OC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: an advanced analysis and visualization platform for ’omics data. PeerJ. 2015;3:e1319. Epub 20151008. pmid:26500826; PubMed Central PMCID: PMC4614810.
  46. 46. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–W6. pmid:33885785; PubMed Central PMCID: PMC8265157.
  47. 47. Faust K, Raes J. CoNet app: inference of biological association networks using Cytoscape. F1000Res. 2016;5:1519. Epub 2016/11/22. pmid:27853510; PubMed Central PMCID: PMC5089131.
  48. 48. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. pmid:14597658; PubMed Central PMCID: PMC403769.
  49. 49. Bi Y, Hubbard C, Purushotham P, Zimmer J. Insights into the structure and function of membrane-integrated processive glycosyltransferases. Curr Opin Struct Biol. 2015;34:78–86. Epub 2015/09/06. pmid:26342143; PubMed Central PMCID: PMC4684724.
  50. 50. Swords WE. Nontypeable Haemophilus influenzae biofilms: role in chronic airway infections. Front Cell Infect Microbiol. 2012;2:97. Epub 2012/08/25. pmid:22919686; PubMed Central PMCID: PMC3417564.
  51. 51. Apicella MA. Nontypeable Haemophilus influenzae: the role of N-acetyl-5-neuraminic acid in biology. Front Cell Infect Microbiol. 2012;2:19. Epub 2012/08/25. pmid:22919611; PubMed Central PMCID: PMC3417534.
  52. 52. Schweda EK, Richards JC, Hood DW, Moxon ER. Expression and structural diversity of the lipopolysaccharide of Haemophilus influenzae: implication in virulence. Int J Med Microbiol. 2007;297(5):297–306. Epub 2007/04/25. pmid:17452015.
  53. 53. Qin Y, Lin G, Chen W, Huang B, Huang W, Yan Q. Flagellar motility contributes to the invasion and survival of Aeromonas hydrophila in Anguilla japonica macrophages. Fish Shellfish Immunol. 2014;39(2):273–9. Epub 2014/05/27. pmid:24859591.
  54. 54. Bajaj JS, Shamsaddini A, Acharya C, Fagan A, Sikaroodi M, Gavis E, et al. Multiple bacterial virulence factors focused on adherence and biofilm formation associate with outcomes in cirrhosis. Gut Microbes. 2021;13(1):1993584. Epub 2021/11/09. pmid:34743650; PubMed Central PMCID: PMC8582993.
  55. 55. Doyle TB, Hawkins AC, McCarter LL. The complex flagellar torque generator of Pseudomonas aeruginosa. J Bacteriol. 2004;186(19):6341–50. Epub 2004/09/18. pmid:15375113; PubMed Central PMCID: PMC516612.
  56. 56. Lo YL, Shen L, Chang CH, Bhuwan M, Chiu CH, Chang HY. Regulation of Motility and Phenazine Pigment Production by FliA Is Cyclic-di-GMP Dependent in Pseudomonas aeruginosa PAO1. PLoS One. 2016;11(5):e0155397. Epub 2016/05/14. pmid:27175902; PubMed Central PMCID: PMC4866697.
  57. 57. Reimmann C, Ginet N, Michel L, Keel C, Michaux P, Krishnapillai V, et al. Genetically programmed autoinducer destruction reduces virulence gene expression and swarming motility in Pseudomonas aeruginosa PAO1. Microbiology (Reading). 2002;148(Pt 4):923–32. Epub 2002/04/05. pmid:11932439.
  58. 58. Calder T, de Souza Santos M, Attah V, Klimko J, Fernandez J, Salomon D, et al. Structural and regulatory mutations in Vibrio parahaemolyticus type III secretion systems display variable effects on virulence. FEMS Microbiol Lett. 2014;361(2):107–14. Epub 2014/10/08. pmid:25288215; PubMed Central PMCID: PMC4675326.
  59. 59. Teeling H, Fuchs BM, Becher D, Klockow C, Gardebrecht A, Bennke CM, et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science. 2012;336(6081):608–11. Epub 2012/05/05. pmid:22556258.
  60. 60. Buchan A, LeCleir GR, Gulvik CA, Gonzalez JM. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nature reviews Microbiology. 2014;12(10):686–98. Epub 2014/08/20. pmid:25134618.
  61. 61. Grieb A, Francis TB, Kruger K, Orellana LH, Amann R, Fuchs BM. Candidatus Abditibacter, a novel genus within the Cryomorphaceae, thriving in the North Sea. Syst Appl Microbiol. 2020;43(4):126088. Epub 2020/07/22. pmid:32690198.
  62. 62. Frade PR, Glasl B, Matthews SA, Mellin C, Serrao EA, Wolfe K, et al. Spatial patterns of microbial communities across surface waters of the Great Barrier Reef. Commun Biol. 2020;3(1):442. Epub 2020/08/17. pmid:32796904; PubMed Central PMCID: PMC7428009.
  63. 63. Choi A, Lee K, Oh HM, Feng J, Cho JC. Litoricola marina sp. nov. Int J Syst Evol Microbiol. 2010;60(Pt 6):1303–6. Epub 2009/08/12. pmid:19667367.
  64. 64. Han Y, Jiao N, Zhang Y, Zhang F, He C, Liang X, et al. Opportunistic bacteria with reduced genomes are effective competitors for organic nitrogen compounds in coastal dinoflagellate blooms. Microbiome. 2021;9(1):71. Epub 2021/03/26. pmid:33762013; PubMed Central PMCID: PMC7992965.
  65. 65. Xia X, Zheng Q, Leung SK, Wang Y, Lee PY, Jing H, et al. Distinct metabolic strategies of the dominant heterotrophic bacterial groups associated with marine Synechococcus. Sci Total Environ. 2021;798:149208. Epub 2021/08/11. pmid:34375229.
  66. 66. Nowinski B, Smith CB, Thomas CM, Esson K, Marin R, Preston CM, et al. Microbial metagenomes and metatranscriptomes during a coastal phytoplankton bloom. Scientific Data. 2019;6(1):129. pmid:31332186
  67. 67. Yang C, Li Y, Zhou B, Zhou Y, Zheng W, Tian Y, et al. Illumina sequencing-based analysis of free-living bacterial community dynamics during an Akashiwo sanguine bloom in Xiamen sea, China. Scientific Reports. 2015;5(1):8476. pmid:25684124
  68. 68. Collado L, Inza I, Guarro J, Figueras MJ. Presence of Arcobacter spp. in environmental waters correlates with high levels of fecal pollution. Environ Microbiol. 2008;10(6):1635–40. Epub 2008/01/25. pmid:18215159.
  69. 69. Morales Medina WR, Eramo A, Fahrenfeld NL. Metabolically Active Prokaryotes and Actively Transcribed Antibiotic Resistance Genes in Sewer Systems: Implications for Public Health and Microbially Induced Corrosion. Microbial Ecology. 2021. pmid:34117524
  70. 70. Park Y, Yu J, Nguyen VK, Park S, Kim J, Lee T. Understanding complete ammonium removal mechanism in single-chamber microbial fuel cells based on microbial ecology. Science of The Total Environment. 2021;764:144231. pmid:33385649
  71. 71. Heider J, Szaleniec M, Sunwoldt K, Boll M. Ethylbenzene Dehydrogenase and Related Molybdenum Enzymes Involved in Oxygen-Independent Alkyl Chain Hydroxylation. J Mol Microbiol Biotechnol. 2016;26(1–3):45–62. Epub 2016/03/10. pmid:26960184.
  72. 72. Zhang Z, Cao H, Song N, Zhang L, Cao Y, Tai J. Long-term hexavalent chromium exposure facilitates colorectal cancer in mice associated with changes in gut microbiota composition. Food Chem Toxicol. 2020;138:111237. Epub 20200304. pmid:32145354.
  73. 73. Gómez-Brandón M, Aira M, Santana N, Pérez-Losada M, Domínguez J. Temporal Dynamics of Bacterial Communities in a Pilot-Scale Vermireactor Fed with Distilled Grape Marc. 2020;8(5):642. Microorganisms. pmid:32354197
  74. 74. Luo S, Ren L, Wu W, Chen Y, Li G, Zhang W, et al. Impacts of earthworm casts on atrazine catabolism and bacterial community structure in laterite soil. Journal of Hazardous Materials. 2022;425:127778. pmid:34823960
  75. 75. Rambo IM, Dombrowski N, Constant L, Erdner D, Baker BJ. Metabolic relationships of uncultured bacteria associated with the microalgae Gambierdiscus. Environ Microbiol. 2020;22(5):1764–83. Epub 2019/11/28. pmid:31775181.
  76. 76. Elifantz H, Horn G, Ayon M, Cohen Y, Minz D. Rhodobacteraceae are the key members of the microbial community of the initial biofilm formed in Eastern Mediterranean coastal seawater. FEMS Microbiology Ecology. 2013;85(2):348–57. pmid:23551015
  77. 77. Dang H, Li T, Chen M, Huang G. Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl Environ Microbiol. 2008;74(1):52–60. Epub 2007/10/30. pmid:17965206; PubMed Central PMCID: PMC2223210.
  78. 78. Campbell AM, Fleisher J, Sinigalliano C, White JR, Lopez JV. Dynamics of marine bacterial community diversity of the coastal waters of the reefs, inlets, and wastewater outfalls of southeast Florida. Microbiologyopen. 2015;4(3):390–408. Epub 2015/03/06. pmid:25740409; PubMed Central PMCID: PMC4475383.
  79. 79. Sushmitha TJ, Rajeev M, Sriyutha Murthy P, Ganesh S, Toleti SR, Karutha Pandian S. Bacterial community structure of early-stage biofilms is dictated by temporal succession rather than substrate types in the southern coastal seawater of India. PLoS One. 2021;16(9):e0257961. Epub 2021/09/28. pmid:34570809; PubMed Central PMCID: PMC8476003.
  80. 80. Rezaei Somee M, Dastgheib SMM, Shavandi M, Ghanbari Maman L, Kavousi K, Amoozegar MA, et al. Distinct microbial community along the chronic oil pollution continuum of the Persian Gulf converge with oil spill accidents. Scientific Reports. 2021;11(1):11316. pmid:34059729
  81. 81. Vezzulli L, Grande C, Reid PC, Helaouet P, Edwards M, Hofle MG, et al. Climate influence on Vibrio and associated human diseases during the past half-century in the coastal North Atlantic. Proc Natl Acad Sci U S A. 2016;113(34):E5062–71. Epub 2016/08/10. pmid:27503882; PubMed Central PMCID: PMC5003230.
  82. 82. Liang J, Liu J, Wang X, Lin H, Liu J, Zhou S, et al. Spatiotemporal Dynamics of Free-Living and Particle-Associated Vibrio Communities in the Northern Chinese Marginal Seas. Appl Environ Microbiol. 2019;85(9). Epub 2019/03/03. pmid:30824453; PubMed Central PMCID: PMC6495765.
  83. 83. Xu W, Gong L, Yang S, Gao Y, Ma X, Xu L, et al. Spatiotemporal Dynamics of Vibrio Communities and Abundance in Dongshan Bay, South of China. Front Microbiol. 2020;11:575287. Epub 2020/12/17. pmid:33324364; PubMed Central PMCID: PMC7726330.
  84. 84. Mohamad N, Mustafa M, Amal MNA, Saad MZ, Md Yasin IS, Al-Saari N. Environmental Factors Associated with the Presence of Vibrionaceae in Tropical Cage-Cultured Marine Fishes. J Aquat Anim Health. 2019;31(2):154–67. Epub 2019/01/18. pmid:30653742.
  85. 85. Takemura AF, Chien DM, Polz MF. Associations and dynamics of Vibrionaceae in the environment, from the genus to the population level. Front Microbiol. 2014;5:38. Epub 2014/02/28. pmid:24575082; PubMed Central PMCID: PMC3920100.
  86. 86. Kopprio GA, Neogi SB, Rashid H, Alonso C, Yamasaki S, Koch BP, et al. Vibrio and Bacterial Communities Across a Pollution Gradient in the Bay of Bengal: Unraveling Their Biogeochemical Drivers. Front Microbiol. 2020;11:594. Epub 2020/05/01. pmid:32351470; PubMed Central PMCID: PMC7174592.
  87. 87. Collins FWJ, Walsh CJ, Gomez-Sala B, Guijarro-Garcia E, Stokes D, Jakobsdottir KB, et al. The microbiome of deep-sea fish reveals new microbial species and a sparsity of antibiotic resistance genes. Gut Microbes. 2021;13(1):1–13. Epub 2021/05/11. pmid:33970781; PubMed Central PMCID: PMC8115496.
  88. 88. Cuadrat RRC, Sorokina M, Andrade BG, Goris T, Davila AMR. Global ocean resistome revealed: Exploring antibiotic resistance gene abundance and distribution in TARA Oceans samples. Gigascience. 2020;9(5). Epub 2020/05/12. pmid:32391909; PubMed Central PMCID: PMC7213576.
  89. 89. Sucato A, Vecchioni L, Savoca D, Presentato A, Arculeo M, Alduina R. A Comparative Analysis of Aquatic and Polyethylene-Associated Antibiotic-Resistant Microbiota in the Mediterranean Sea. Biology (Basel). 2021;10(3). Epub 2021/04/04. pmid:33800749; PubMed Central PMCID: PMC8001005.
  90. 90. Chen YM, Holmes EC, Chen X, Tian JH, Lin XD, Qin XC, et al. Diverse and abundant resistome in terrestrial and aquatic vertebrates revealed by transcriptional analysis. Sci Rep. 2020;10(1):18870. Epub 2020/11/04. pmid:33139761; PubMed Central PMCID: PMC7608656.
  91. 91. Hao X, Zhu J, Rensing C, Liu Y, Gao S, Chen W, et al. Recent advances in exploring the heavy metal(loid) resistant microbiome. Comput Struct Biotechnol J. 2021;19:94–109. Epub 2021/01/12. pmid:33425244; PubMed Central PMCID: PMC7771044.
  92. 92. Jiang H, Yu T, Yang Y, Yu S, Wu J, Lin R, et al. Co-occurrence of Antibiotic and Heavy Metal Resistance and Sequence Type Diversity of Vibrio parahaemolyticus Isolated From Penaeus vannamei at Freshwater Farms, Seawater Farms, and Markets in Zhejiang Province, China. Frontiers in Microbiology. 2020;11. pmid:32676056
  93. 93. Adekanmbi AO, Olawuni OG, Olaposi AV. Metal contamination and coexistence of metal and antibiotic resistance in Vibrio species recovered from aquaculture ponds with and without history of antibiotic usage in Southwest Nigeria. Bulletin of the National Research Centre. 2021;45(1):122.
  94. 94. Samtiya M, Matthews KR, Dhewa T, Puniya AK. Antimicrobial Resistance in the Food Chain: Trends, Mechanisms, Pathways, and Possible Regulation Strategies. Foods. 2022;11(19):2966. Epub 2022/10/15. pmid:36230040; PubMed Central PMCID: PMC9614604.
  95. 95. Han Y, Wang J, Zhao Z, Chen J, Lu H, Liu G. Fishmeal Application Induces Antibiotic Resistance Gene Propagation in Mariculture Sediment. Environ Sci Technol. 2017;51(18):10850–60. Epub 2017/08/31. pmid:28853291.
  96. 96. Goldberg SJ, Nelson CE, Viviani DA, Shulse CN, Church MJ. Cascading influence of inorganic nitrogen sources on DOM production, composition, lability and microbial community structure in the open ocean. Environ Microbiol. 2017;19(9):3450–64. Epub 20170814. pmid:28618153.
  97. 97. Maier S, Kratz AM, Weber J, Prass M, Liu F, Clark AT, et al. Water-driven microbial nitrogen transformations in biological soil crusts causing atmospheric nitrous acid and nitric oxide emissions. ISME J. 2021. Epub 20211111. pmid:34764454.
  98. 98. Wu C, Zhao Z, Liu Y, Zhu X, Liu M, Luo P, et al. Type III Secretion 1 Effector Gene Diversity Among Vibrio Isolates From Coastal Areas in China. Front Cell Infect Microbiol. 2020;10:301. Epub 20200618. pmid:32637366; PubMed Central PMCID: PMC7318850.
  99. 99. Chen Y-H, Chiang P-W, Rogozin DY, Degermendzhy AG, Chiu H-H, Tang S-L. Salvaging high-quality genomes of microbial species from a meromictic lake using a hybrid sequencing approach. Communications Biology. 2021;4(1):996. pmid:34426638
  100. 100. Ciuffreda L, Rodriguez-Perez H, Flores C. Nanopore sequencing and its application to the study of microbial communities. Comput Struct Biotechnol J. 2021;19:1497–511. Epub 20210307. pmid:33815688; PubMed Central PMCID: PMC7985215.
  101. 101. Latorre-Perez A, Villalba-Bermell P, Pascual J, Vilanova C. Assembly methods for nanopore-based metagenomic sequencing: a comparative study. Sci Rep. 2020;10(1):13588. Epub 20200812. pmid:32788623; PubMed Central PMCID: PMC7423617.
  102. 102. Michalska-Smith M, Song Z, Spawn-Lee SA, Hansen ZA, Johnson M, May G, et al. Network structure of resource use and niche overlap within the endophytic microbiome. ISME J. 2022;16(2):435–46. Epub 2021/08/21. pmid:34413476; PubMed Central PMCID: PMC8776778.
  103. 103. Poudel R, Jumpponen A, Schlatter DC, Paulitz TC, Gardener BB, Kinkel LL, et al. Microbiome Networks: A Systems Framework for Identifying Candidate Microbial Assemblages for Disease Management. Phytopathology. 2016;106(10):1083–96. Epub 2016/08/03. pmid:27482625.
  104. 104. Goodman AZ, Papudeshi B, Doane MP, Mora M, Kerr E, Torres M, et al. Epidermal Microbiomes of Leopard Sharks (Triakis semifasciata) Are Consistent across Captive and Wild Environments. Microorganisms. 2022;10(10). Epub 20221021. pmid:36296361; PubMed Central PMCID: PMC9610875.
  105. 105. Langille MGI, Hsiao WWL, Brinkman FSL. Detecting genomic islands using bioinformatics approaches. Nature Reviews Microbiology. 2010;8(5):373–82. pmid:20395967
  106. 106. Juhas M, van der Meer JR, Gaillard M, Harding RM, Hood DW, Crook DW. Genomic islands: tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol Rev. 2009;33(2):376–93. Epub 2009/01/31. pmid:19178566; PubMed Central PMCID: PMC2704930.
  107. 107. Zhang S-P, Wang Q, Quan S-W, Yu X-Q, Wang Y, Guo D-D, et al. Type II toxin–antitoxin system in bacteria: activation, function, and mode of action. Biophysics Reports. 2020;6(2):68–79.
  108. 108. Kamruzzaman M, Wu AY, Iredell JR. Biological Functions of Type II Toxin-Antitoxin Systems in Bacteria. Microorganisms. 2021;9(6). Epub 2021/07/03. pmid:34208120; PubMed Central PMCID: PMC8230891.
  109. 109. Fraikin N, Goormaghtigh F, Van Melderen L. Type II Toxin-Antitoxin Systems: Evolution and Revolutions. J Bacteriol. 2020;202(7). Epub 2020/01/15. pmid:31932311; PubMed Central PMCID: PMC7167474.
  110. 110. Schureck MA, Maehigashi T, Miles SJ, Marquez J, Dunham CM. mRNA bound to the 30S subunit is a HigB toxin substrate. RNA. 2016;22(8):1261–70. Epub 2016/06/17. pmid:27307497; PubMed Central PMCID: PMC4931118.
  111. 111. Schureck MA, Dunkle JA, Maehigashi T, Miles SJ, Dunham CM. Defining the mRNA recognition signature of a bacterial toxin protein. Proceedings of the National Academy of Sciences. 2015;112(45):13862–7. pmid:26508639
  112. 112. Liu Y, Gao Z, Liu G, Geng Z, Dong Y, Zhang H. Structural Insights Into the Transcriptional Regulation of HigBA Toxin–Antitoxin System by Antitoxin HigA in Pseudomonas aeruginosa. Frontiers in Microbiology. 2020;10(3158). pmid:32038588
  113. 113. Schureck MA, Meisner J, Hoffer ED, Wang D, Onuoha N, Ei Cho S, et al. Structural basis of transcriptional regulation by the HigA antitoxin. Molecular Microbiology. 2019;111(6):1449–62. pmid:30793388