Seasonal Variations in the Bacterial Community of Gwangyang Bay Seawater
Seasonal Variations in the Bacterial Community of Gwangyang Bay Seawater
Journal of Life Science. 2014. May, 24(5): 522-531
Copyright © 2014, Korean Society of Life Science
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : February 26, 2014
  • Accepted : May 21, 2014
  • Published : May 30, 2014
Export by style
Cited by
About the Authors
성찬, 박
Woogene B&G Rearch Institute, Seoul 150-972, Korea
지희, 이
Department of Biology, Sunchon National University, Suncheon 540-950, Korea
주원, 강
Department of Biology, Sunchon National University, Suncheon 540-950, Korea
근식, 백
Division of Life Sciences, Korea Basic Science Institute, Daejeon 305-806, Korea
치남, 성
Department of Biology, Sunchon National University, Suncheon 540-950, Korea

Seasonal variations in the bacterial community of Gwangyang Bay seawater were analyzed using both isolation and cultivation-independent methods. Amplified rDNA restriction analysis was applied to 200 bacterial isolates. Bacterial isolates were composed of four phyla: Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes. Pyrosequencing was conducted, in addition to denaturing gradient gel electrophoresis (DGGE) of genomic DNA extracted directly from the water samples. The bacterial sequences obtained by pyrosequencing of 16S rRNA genes consisted of 24 phyla in the spring and summer, 39 in the fall, and 32 in the winter. The diversity index was high in the fall, whereas the dominancy index was high in the spring. In the spring, phylum Firmicutes was dominant, whereas phylum Proteobacteria dominated in the other three seasons. The second most dominant phyla were Proteobacteria in the spring, Firmicutes in the summer, and Bacteroidetes both in the fall and winter. Bacilliaceae was the most predominant family in the spring. Rhodobacteraceae and Bacilliaceae dominated in the summer, and Rhodobacteraceae dominated in the winter. Neither was dominant in the fall Twenty-seven bands purified from DGGE profiles were cloned and analyzed phylogenetically. In the spring, phylum Firmicutes dominated, followed by Proteobacteria. Proteobacteria dominated in all other seasons. Thus, two cultivation-independent methods for determination of seasonal variation patterns at the phylum level were in accordance with each other.
Gwangyang bay is located in the southern coastal region surrounded by Yeosu Peninsula and Namhae Island, Republic of Korea (34° 51′ 16″~34° 56′ 55″ N, 127° 37′ 23″~127° 50′ 86″ E). Gwangyang bay, which receives freshwater via Seomjin River, is semi-closed and linked to the South Sea via two channels, Yeosu and Noryang, to the south and the east, respectively. The surface area of the bay is 230 km 2 and contains 20 islets. A large-scale steel works, container port, and heavy-chemical estate were constructed in the coastal region of the bay. In recent years, this area has been designated as a special management area for coastal pollution, receiving great attention from the public due to its increasing environmental problems. There have been several investigations into the functions of environmental factors, water quality, and planktons and benthos [ 2 , 6 , 22 ], whereas studies on microorganisms are few. Therefore, we carried out systematic cataloguing of the microbial community of Gwangyang bay seawater.
The objective of this study was to examine seasonal variations in the bacterial community of Gwangyang bay seawater. To achieve this goal, cultivation–dependent and –independent analyses, such as pyrosequencing and denaturing gradient gel electrophoresis (DGGE) followed by phylogenetic analysis, were performed.
Materials and Methods
- Sampling and environmental factors
Seawater samples were taken four times (April, July, October 2008, and February 2009) at a surface depth of 1 m in the mouth of Yeosu channel in Gwangyang bay, Republic of Korea (34° 51′ 20.3″ N, 128° 47′ 29.7″ E). Samplings were carried out at ebb tide. Temperature, pH, salinity, dissolved oxygen, heterotrophic bacterial number, and chlorophyll a content of the samples were checked according to Standard Methods [1] .
- Isolation of bacteria
Isolation of bacteria was achieved by the standard dilution plating technique using marine agar 2216 (MA, Becton Dickinson) at 25℃ for 7 days. Purified isolates were routinely cultured on MA and maintained as a glycerol suspension (20%, w/v) at –80℃.
- Amplified rDNA restriction analysis (ARDRA) of isolates
Bacterial DNA preparation and PCR amplification were carried out as previously described [7] . Universal primers such as 27F ( E. coli numbering 8~27; 5′-AGA GTT TGA TCM TGG CTC AG-3′) and 1492R ( E. coli numbering 1492~1510; 5’-GGY TAC CTT GTT ACG ACT T-3’) were used for amplification of the 16S rRNA gene [21] . PCR amplification of the nearly full-length 16S rRNA gene was performed in 50 μl reaction mixtures. PCR products were detected after 1% agarose gel electrophoresis.
Amplified 16S rRNA genes were digested with Hae III (GG ↓CC; TaKaRa) as previously described [3] . The resulting fragments were electrophoresed on a 4% NuSieve 3:1 Agarose (Cambrex, USA) gel, followed by EtBr staining and photography by UV transillumination. Similarities were estimated using the Gelcompar II program 4.6 software (Applied Maths).
16S rRNA genes obtained from representative clones of all phylotypes were partially sequenced using automated DNA sequencer (ABI3700, Applied Biosystems) with the sequencing primer 27F [21] . DNA sequences were edited using PHYDIT. All sequences were analyzed for the presence of chimera using the CHIMERA CHECK program ( [35] ). Partial 16S rRNA gene sequences were identified using EzTaxon-e database ( [19] ).
- DNA extraction from water sample
Total microbial DNA was directly extracted from the seawater samples. Suspended solids of water samples were removed with filter paper (Whatman no. 2), followed by concentration using a filter unit (0.22 μm) (Sterivex-GP, Millipore). DNA extraction was achieved by following previously described method [34] . The quality and concentration of DNA were examined by standard electrophoresis, and spectrophometry using UV/vis spectrophotometer spectrophotometer (Ultraspec2100 pro, GE Healthcare Life Sciences). DNA concentration was adjusted to 300 ng/μl for the next studies.
- Pyrosequencing
The four genomic DNA samples were diluted to 1:5 and 1:100. These eight diluted genomic DNA samples were then amplified separately using the bacterial primer pair 27F ( E. coli numbering 8~27; 5′- xxxAC GAG TTT GAT CMT GGC TCA G-3′) and 518R ( E. coli numbering 518~534; 5’- xxxAC WTT ACC GCG GCT GCT GG-3’) [9] , generating the PCR V3 region of the 16S rRNA fragments of c. 400 bp, where xxx was designed for the sample identification barcoding key and AC represents linker. Each 50 μl PCR reaction included the following: dNTP (2.5 mM), 5 μl; Ex Taq DNA polymerase (5 U), 0.25 μl; template DNA (100 ng/μl), 1 μl; BSA (4 μg/μl), 5 μl; reaction buffer [10 mM Tris-HCl, 40 mM KCl, 2.5 mM MgCl 2 ], 5 μl; and each primer (20 μM), 1 μl. PCR was conducted on a thermal cycler (TP600, TaKaRa) under the following conditions: 94℃ for 5 min; 30 cycles consisting of 94℃ for 30 s, 50℃ for 30 s, and 72℃ for 40 s; and a final extension at 72℃ for 10 min. PCR products were detected after 1% agarose gel electrophoresis, purified using the gel extraction kit (QIAquick, QIAGEN), and pooled to obtain eight amplicon libraries. Amplicon length and concentration were estimated, and an equimolar mix of all eight amplicon libraries was used for pyrosequencing. Pyrosequencing of the eight amplicon libraries on the Roche 454 GS FLX Titanium system at Macrogen (Korea) resulted in 150,264 reads satisfying the sequence quality criteria employed [11] .
- Pyrosequencing data analysis
Initially, the raw sequences obtained from pyrosequencing were processed using an adequate software (GL FLX, Roche) for sorting according to the key (i.e., sequences from the suspended and attached samples), with low quality sequences discarded, and the primer sequences trimmed. In addition, sequences less than 300 bp long were discarded. The processed sequences were subjected to taxonomical classification using the EzTaxon-e database (http://eztaxon- /[19] ). Clustering of tags into operational taxonomic units (OTUs) and generation of rarefaction curves were carried out using the pyrosequencing pipeline of the Ribosomal RNA database project (http://pyro. [10] ). Classification criteria based on similarity for the species, genus, family, order, class, and phylum are ≥97%, ≥94%, ≥90%, 85%, ≥80%, and ≥75%, respectively. The program DOTUR [32] was employed to calculate diversity estimates, namely the abundance-based coverage estimator ACE [5] and species diversity estimator Chao1 [4] , by using sampling without replacement. Rarefaction curves were calculated by using the analytical approximation algorithm described by Hurlbert [16] , and 95% confidence intervals were estimated as described by Heck et al . [14] . Assignment of a sequence to the same group (phylotype) was determined with similarities ≥97%.
- Denaturing gradient gel electrophoresis
Amplification of partial 16S rRNA genes for DGGE was performed the same as in pyrosequencing, except that 803R ( E. coli numbering 803~826; 5’-GAA AAA TCC CCG CAG CAA CCA TAG-3’) was used as reverse primer. Purified DNAs were amplified again for the V3 region of 16S rRNA genes. Primer pairs used were 341F-GC (5′-CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GCC TAC GGG AGG CAG CAG-3′) and 518R (5′-WTT ACC GCG GCT GCT GG-3′). Preparation of the reaction mixture was the same as above, except for the primers. PCR conditions were as follows: initial denaturation at 94℃ for 5 min, followed by 20 cycles of 30 s at 94℃, annealing at 65-55℃ (in the first cycle, annealing was conducted at 65℃, after which the temperature was reduced by 0.5℃ in each following cycle) for 30 s, and extension at 72℃ for 30 s. This procedure was followed by 10 cycles of 30 s at 94℃, 30 s at 55℃, and 30 s at 72℃. Finally, a primer extension at 72℃ for 7 min was conducted. PCR products were detected after 1% agarose gel electrophoresis and purified using the QIAquick gel extraction kit.
DGGE was conducted on the D-Code System (Bio-Rad) using 8% (w/v) polyacrylamide gel (acrylamide: bisacrylamide ratio of 29:1) in 1× TAE buffer (20 mM Tris acetate, 10 mM acetic acid, 0.5 mM EDTA, pH 8.0). PCR products containing approximately equal amounts of DNA of similar sizes were separated on a gel containing an increasing linear gradient of denaturants (urea and formamide) from 40% at the top of the gel to 60% at the bottom. Electrophoresis was conducted at 60℃ by applying 60 V to the gel for 15.5 hr. After electrophoresis, gels were stained with 1× TAE buffer containing 1× SYBR Green (Sigma) and then processed using Gel Doc XR system (Bio-Rad).
- DGGE profile analysis
Twenty-seven bands representing dominant taxa in each sample were selected for further analyses. DNAs were eluted, purified, and cloned using pGEM-T Easy Vector, followed by transformation into E. coli DH5α competent cells. Clones were selected by white colony selection on LB agar containing X-Gal (5-bromo-4-chloro-3-indoly- β -D-galactopyranoside) (Promega, USA), IPTG (isopropyl- β -D-thiogalacto-pyranoside) (Promega), and ampicillin (50 mg/ml). Selected colonies (three colonies per band) were reamplified directly using pGTf (5′- TAC GAC TCA CTA TAG GGC GA -3′) and pGTr (5′- ACT CAA GCT ATG CAT CCA ACG -3′) primers (8). PCR conditions were the same as in the amplification of partial 16S rRNA genes. Sequencing was performed with an automated DNA sequencer (ABI3700, Applied Biosystems) according to the manufacturer’s instructions. Sequence analysis was carried out using the EzTaxon-e database (http://eztaxon-e.ezbiocloud. net/ [19] ) and PHYDIT.
- Cluster analysis
Cluster analysis of DGGE profiles was performed by using GelCompar II version 4.6 software (Applied Maths). Quantitative information of each OTU derived from pyrosequencing was used to calculate the similarity of the bacterial composition between seasons using UniFrac distance metric [24] . The unweighted pair group method with averages (UPGMA [33] ) was used to express the clustering result.
- Environmental factors
Seasonal variations in environmental parameters are described in Table 1 . Heterotrophic bacterial number and chlorophyll a content were higher in autumn and winter, respectively. The pH value and salinity ranged from 8.0 to 8.3% and from 29.6 to 31.9%, respectively.
Physicochemical properties of seawater samples collected in Gwangyang bay
PPT Slide
Lager Image
※ Tem., temperature; HB, heterotrophic bacteria; Chl. a, chlorophyll a
- Bacterial community discovered by ARDRA with isolates
From the isolates, 80 phylotypes, including 39 singlemember phylotypes, were produced by ARDRA using Hae III endonuclease. Distribution patterns of the phylotypes produced from bacterial isolates are listed in Table 2 . The community of bacterial isolates was composed of four phyla: Firmicutes (101 strains; 50.5%), Proteobacteria (66; 33%), Actinobacteria (22; 11%), and Bacteroidetes (11; 5.5%). Firmicutes was dominant in spring and summer, whereas Proteobacteria was in autumn and winter (Table 2 ). Isolates were affiliated as 46 genera composing 27 families. Members of the genus Bacillus dominated in all seasons, whereas those of the genus Erythrobacter did in autumn and winter.
Bacterial community of Gwangyang bay seawater discovered by ARDRA with isolates
PPT Slide
Lager Image
Bacterial community of Gwangyang bay seawater discovered by ARDRA with isolates
- Bacterial community discovered by pyrosequencing
A total of 140,987 bacterial sequences longer than 300 bp were obtained by pyrosequencing, and the average length was 439 bp. Taxonomic assignment of each sequence was carried out in an automated manner, and the results were a list of the taxonomic assignment of each amplicon as well as the number of amplicons in each major taxonomic rank. Fig. 1 shows the rarefaction analysis based on best matches between each tag and sequences in the EzTaxon-e database database [19] , along with their frequency of recovery. These rarefaction curves describe unprecedented levels of bacterial complexity for marine samples, yet none has reached the curvilinear or plateau phase. The statistical richness estimates of ACE and Chao1 determined by the DOTUR program are given in Table 3 . Shannon-Weaver diversity, ACE, and Chao1 values were higher in autumn, whereas Simpson dominance index was higher in spring. Taxonomic/phylogenetic assignment of bacterial sequences indicated that the sample contained 19 divisions, 21 candidate divisions, and five new divisions. Seasonal variations in the number of divisions were as follows: 24, 24, 39, and 32 in spring, summer, autumn, and winter, respectively. Phylum Proteobacteria was the most abundant (48.85% of total bacteria), followed by Firmicutes (36.72%), Bacteroidetes (6.23%), Planctomycetes (3.12%), and Actinobacteria (1.96%). The other 38 divisions presented only in minor components (data not shown). In the phylum Proteobacteria, class Alphaproteobacteria (80.71% of Proteobacteria) was the most abundant, followed by the class Gammaproteobacteria. A small portion of the reads (1,268 reads) was assigned to 21 candidate divisions for which no cultural representatives are available (data not shown). Candidate division OD1 was the most abundant group (0.24% of total bacteria), followed by candidate divisions OP11 (0.22%), SAR202 (0.12%), and GN02 (0.11%).
PPT Slide
Lager Image
Rarefaction curves for sequence reads obtained from the pyrosequencing analysis using DNA extracts obtained from Gwangyang bay seawter. Rarefaction curves were calculated from operational taxonomic unit (OTU) analysis of 97% identical sequence type.
Summary of pyrosequencing data of DNA extracts from Gwangyang bay seawater
PPT Slide
Lager Image
Summary of pyrosequencing data of DNA extracts from Gwangyang bay seawater
In spring, phylum Firmicutes dominated (85.24% of the total reads) and consisted of 12 genera with the genus Bacillus predominating (81.5%) (Fig. 2A ). In summer, phylum Proteobacteria (52.18%) dominated, followed by phyla Firmicutes (39.68%), Bacteroidetes (2.91%), and Actinobacteria (2.64%) (Fig. 2B ). Phylum Proteobacteria consisted of 171 genera with the genus Tropicibacter (family Rhodobacteraceae ) predominating (7,536 reads; 15.41%). Phylum Firmicutes consisted of nine genera with the genus Bacillus predominating (35.04%). In autumn, phylum Proteobacteria (69.85%) dominated, followed by phyla Bacteroidetes (9.84%), Planctomycetes (6.94%), Acidobacteria (2.94%), Actinobacteria (2.31%), and Firmicutes (1.89%) (Fig. 2C ). Phyla Deinococcus- Thermus and Tenericutes appeared exclusively in this season. In winter, phylum Proteobacteria (77.74%) dominated, followed by phyla Bacteroidetes (15.54%), Actinobacteria (2.27%), and Verrucomicrobia (1.34%) (Fig. 2D ). The family Rhodobacteraceae (class Alphaproteobacteria) predominating (50.71%). Phyla Fibrobacteres and Spirochaetes appeared exclusively in this season.
PPT Slide
Lager Image
Relative abundance of bacterial pyrosequences on the phylum level obtained from Gwangyang bay seawater. Percentage of each phylum was calculated and shown if higher than 1%. (A), spring; (B), summer; (C), autumn; (D), winter.
Variations between seasons in the SAR11 clade, which is prevalent in seawater DNA (Morris et al. , 2002), slightly occurred in spring (0.22%) and summer (0.3%), whereas they were more abundant in autumn (1.47%) and winter (8.5%).
- Dominant bacterial community discovered by DGGE
Overall bacterial compositions of DGGE profiles are summarized in Fig. 3 and Table 4 . Six families in two phyla (Firmicutes and Proteobacteria) were detected in spring. Clones from intensive bands were affiliated with the genera Bacillus, Sphingobium, Sphingomonas, Acinetobacter, Tateyamaria , and Nevskia . In summer, three phyla (Proteobacteria, Firmicutes, and Actinobacteria) composed of six families were detected. The intensive bands represent the genera Seohaeicola, Tateyamaria, Nautella, Sulfitobacter, Thiohalocapsa, Mesorhizobium, Bacillus, Oceanobacillus, Microlunatus , and Propionibacterium . Members of phyla Proteobacteria and Actinobacteria increased in amount, whereas that of Firmicutes decreased. In autumn, the most dominant bacteria were the members of phylum Proteobacteria. Phyla Actinobacteria and Bacteroidetes were detected in only one band each. In winter, four phyla composed of seven families were detected, and the members of families Rhodobacteraceae and Chromatiaceae of phylum Proteobacteria were abundant.
PPT Slide
Lager Image
DGGE profiles of bacterial DNA extracted from Gwangyang bay seawater. Lanes: sp, spring; su, summer; au, autumn; wi, winter.
Sequence information for the bands obtained from DGGE using DNA extracts from Gwangyang bay seawater
PPT Slide
Lager Image
+: Detected sequences in the clone obtained from DGGE band.
- Cluster analysis
Bacterial composition in spring derived from DGGE was different from that of the other seasons (Fig. 4 ). This relationship was also confirmed by cluster analysis using the pyrosequencing data, in which the diversity index of spring was lower than that of the other seasons.
PPT Slide
Lager Image
Relatedness of seasonal bacterial communities of Gwangyang bay seawater. Data are obtained from pyrosequencing (A) and DGGE (B).
Recently, most studies on the bacterial diversity have been carried out based on cultivation-independent 16S rRNA gene clones. For sequencing of the 16S rRNA gene, a modified Sanger method or next-generation techniques [ 12 , 20 , 25 ] are typically used. However, there have been few comprehensive studies using isolation techniques. A previous comprehensive study [25] compared the retrieving capacity of the cultivation-dependent method with that of pyrosequencing, a cultivation-independent method, still having weakness not obtaining pure isolates. Instead, it only used DNA information from viable cells colonized on plates, enriched culture, or grown in microtiter plates. Meanwhile, this study used cultivation-dependent and -independent methods to investigate the bacterial community structure of Gwangyang bay seawater. This study is the first to comprehensively study the bacterial diversity of Gwangyang bay seawater.
The bacterial community of Gwangyang bay seawater based on pyrosequencing data was composed of 45 phyla, whereas that of the Yellow Sea contains 53 phyla, including candidate or new divisions [25] . This can be attributed to the different environmental conditions and sampling methods between the two experiments; i.e. two depth samples were used in the latter while four seasonal samples at a single site were used in this study. Bacterial communities of seawater mainly consisted of Proteobacteria (α-subdivision, β-subdivision, γ-subdivision), Actinobacteria, and Bacteroidetes [ 15 , 20 , 23 , 25 ]. In addition to the common predominant members, dominance of Firmicutes [ 15 , 18 , 23 ] and Cyanobacteria [23] was reported. In this study, dominance of Firmicutes was a special case, especially in spring. This was due to the different environmental conditions, such as water temperature and content of chrolophyll a. In addition, geographical properties might have altered bacterial composition, as our sampling site was located in a bay influenced by the Seumjin river. In addition to common constituents, Planctomycetes was considerably detected in spring, summer, and autumn in amounts of 3.84, 2.05, and 6.96%, respectively. Abundance of Planctomycetes has been reported in marine sediments of the Eastern Mediterranean Sea [ 30 , 31 ]. However, members of phylum Cyanobacteria were found only in small portions (at most 1.07%) in summer. The presence of only a small number of environmental uncultured bacteria that were distantly related to cultivated ones in this study was also reported in the Yellow Sea [25] .
Bacterial communities derived from the cultivation-dependent analysis and DGGE profiles revealed only four phyla, Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes. While pyrosequencing data showed the phyla Planctomycetes, Acidobacteria, Cyanobacteria, and Verrucomicrobia exceeded 1% at least in any one season, the four phyla were not detected in both of the cultivation-dependent and DGGE analyses. Nevertheless, four common dominant phyla were retrieved in the pyrosequencing analysis
Although only a small number of environmental parameters were obtained, the content of chlorophyll a and heterotrophic bacterial abundance were related to bacterial diversity. More information on environmental factors such as nutrients and metal ions should be obtained and compared with the bacterial community.
ARDRA with cultivated bacterial strains revealed only a small portion of the bacterial community in this study. Moreover, members of phylum Firmicutes were the most predominant. However, this analysis showed the distribution tendency of the bacterial composition; the number of Proteobacteria increased while that of Firmicutes decreased from spring to winter. Estimation of bacterial composition based on the isolation technique had a limitation in terms of the number of isolates. In addition, there were several biases in the cultivation and selection stages; most of the bacteria were uncultivable under single culture medium and incubation conditions, and slow-growing bacteria could be rejected in the selection process [ 13 , 17 ]. Despite the limitations of the culture-dependent analysis, several strains were found to be novel species, which were supported by the comparative sequence similarity and phylogenetic analyses. Four species have already been validly published: Marinivirga aestuarii (strain KYW371 [28] ), Altererythrobacter namhicola (strain KYW48 [26] ), Pseudidiomarina aestuarii (strain KYW314 [29] ), and Shewanella marina (strain C4 [27] ).
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education, Science and Technology (No. 2010-0006861).
1989 Standard Methods for the Examination of Water and Wastewater
Baek S. H. , Kim D. S. , Hyun B. G. , Choi H. W. , Kim Y. O. 2011 Characteristics of horizontal community distribution and nutrient limitation on growth rate of phytoplankton during a winter in Gwangyang bay, Korea Ocean Polar Res 33 99 - 111    DOI : 10.4217/OPR.2011.33.2.099
Baik K. S. 2005 Prokaryotic diversity of Woopo wetland estimated by phylogenetic analysis
Chao A. 1984 Nonparametric-estimation of the number of classes in a population Scand J Stat 11 265 - 270
Chao A. , Lee S. M. 1992 Estimating the number of classes via sample coverage J Am Stat Assoc 87 210 - 217    DOI : 10.1080/01621459.1992.10475194
Choi C. G. , Huh S. H. 2008 Composition of marine algal community at the intertidal zone in Gwangyang bay, South Sea, Korea Korran J Fisher Aqua Sci 41 201 - 207    DOI : 10.5657/kfas.2008.41.3.201
Chun J. , Goodfellow M. 1995 A phylogenetic analysis of the genus Nocardia with 16S rRNA gene sequences Int J Syst Bacteriol 45 240 - 245    DOI : 10.1099/00207713-45-2-240
Chun J. , Huq A. , Colwell R. R. 1999 Analysis of 16S-23S rRNA intergenic spacer regions of Vibrio cholerae and Vibrio mimicus Appl Environ Microbiol 65 2202 - 2208
Chun J. , Kim K. Y. , Lee J.-H. , Choi Y. 2010 The analysis of oral microbial communities of wild-type and Toll-like receptor 2-deficient mice using a 454 GS FLX Titanium Pyrosequencer BMC Microbiol 10 101 -    DOI : 10.1186/1471-2180-10-101
Cole J. R. , Wang Q. , Cardenas E. , Fish J. , Chai B. , Farris R. J. , Kulam-Syed-Mohideen A. S. , McGarrell D. M. , Marsh T. , Garrity G. M. , Tiedje J. M. 2009 The ribosomal database project: improved alignments and new tools for rRNA analysis Nucleic Acids Res 37 141 - 145
Droege M. , Hil B. 2008 The Genome Sequencer FLX System-longer reads, more applications, straight forward bioinformatics and more complete data sets J Biotechnol 31 3 - 10
Galand P. E. , Bourrain M. , De Maistre E. , Catala P. , Desdevises Y. , Elifantz H. , Kirchman D. L. , Lebaron P. 2012 Phylogenetic and functional diversity of bacteria and archaea in a unique stratified lagoon, the clipperton atoll (N Pacific) FEMS Microbiol Ecol 82 616 - 628    DOI : 10.1111/j.1574-6941.2012.01432.x
Ghosh A. , Maity B. , Chakrabarti K. , Chattopadhyay D. 2007 Bacterial diversity of East Calcutta wet land area: possible identification of potential bacterial population for different biotechnological uses Microb Ecol 54 452 - 459    DOI : 10.1007/s00248-007-9244-z
Heck K. L. , Van Belle G. , Simberloff D. 1975 Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size Ecol 56 1459 - 1461    DOI : 10.2307/1934716
Henriques I. S. , Almeida A. , Cunha A. , Correia A. 2004 Molecular sequence analysis of prokaryotic diversity in the middle and outer sections of the Portuguese estuary Ria de Aveiro FEMS Microbiol Ecol 49 269 - 279    DOI : 10.1016/j.femsec.2004.04.003
Hurlbert S. H. 1971 The nonconcept of species diversity: a critique and alternative parameters Ecol 52 577 - 586    DOI : 10.2307/1934145
Joshi A. A. , Kanekar P. P. , Kelkar A. S. , Shouche Y. S. , Vani A. A. , Borgave S. B. , Sarnaik S. S. 2008 Cultivable bacterial diversity of alkaline Lonar lake, India Microb Ecol 55 163 - 172    DOI : 10.1007/s00248-007-9264-8
Kim B.-S. , Kim B. K. , Lee J.-H. , Kim M. , Lim Y. W. , Chun J. 2008 Rapid phylogenetic dissection of prokaryotic structure in tidal flat using pyrosequencing J Microbiol 46 357 - 363    DOI : 10.1007/s12275-008-0071-9
Kim O. S. , Cho Y. J. , Lee K. , Yoon S. H. , Kim M. , Na H. , Park S. C. , Jeon Y. S. , Lee J. H. , Yi H. , Won S. , Chun J. 2012 Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species Int J Syst Evol Microbiol 62 716 - 721    DOI : 10.1099/ijs.0.038075-0
Kirchman D. L. , Cottrell M. T. , Lovejoy C. 2010 The structure of bacterial communities in the western Arctic Ocean as revealed by pyrosequencing of 16S rRNA genes Environ Microbiol 12 1132 - 1143    DOI : 10.1111/j.1462-2920.2010.02154.x
Lane D. J. , Stackebrandt E. , Goodfellow M. 1991 16S/23S rRNA sequencing;Nucleic acid techniques in bacterial systematics Wiley Chichester, U.K.
Lee Y. S. , Kang C.-K. , Choi Y.-K. , Leem S. Y. 2007 Origin and spatial distribution of organic matter at Gwangyang bay in the fall Ocean Science J 12 1 - 8
Ling J. , Dong J. D. , Wang Y. S. , Zhang Y. Y. , Deng C. , Lin L. , Wu M. L. , Sun F. L. 2012 Spatial variation of bacterial community structure of the Northern South China Sea in relation to water chemistry Ecotoxicol 21 1669 - 1679    DOI : 10.1007/s10646-012-0941-0
Lozupone C. , Knight R. 2005 UniFrac: a new phylogenetic method for comparing microbial communities Appl Environ Microbiol 71 8228 - 8235    DOI : 10.1128/AEM.71.12.8228-8235.2005
Na H. , Kim O.-S. , Yoon S.-H. , Kim Y. , Chun J. 2011 Comparative approach to capture bacterial diversity of coastal waters J Microbiol 49 729 - 740    DOI : 10.1007/s12275-011-1205-z
Park S. C. , Baik K. S. , Choe H. N. , Lim C. H. , Kim H. J. , Ka J. O. , Seong C. N. 2011 Altererythrobacter namhicola sp. nov., and Altererythrobacter aestuarii sp. nov., isolated from seawater Int J Syst Evol Microbiol 61 709 - 715    DOI : 10.1099/ijs.0.021196-0
Park S. C , Baik K. S. , Kim M. S. , Kim D. , Seong C. N. 2009 Shewanella marina sp. nov., isolated from seawater Int J Syst Evol Microbiol 59 1888 - 1894    DOI : 10.1099/ijs.0.005470-0
Park S. C. , Choe H. N. , Hwang Y. M. , Baik K. S. , Kim S. N. , Lee Y. S. , Jung J. S. , Seong C. N. 2013 Marinivirga aestuarii gen. nov., sp. nov., a novel member of the family Flavobacteriaceae, isolated from marine environment, and emended descriptions of the genera Hyunsoonleella, Jejuia and Pontirhabdus and Hyunsoonleella jejuensis, Jejuia pallidilutea and Pontirhabdus pectinivorans Int J Syst Evol Microbiol 63 1524 - 1531    DOI : 10.1099/ijs.0.042341-0
Park S. C. , Lim C. H. , Baik K. S. , Lee K. H. , Lee J. S. , Seong C. N. 2010 Pseudidiomarina aestuarii sp. nov., a marine bacterium isolated from shallow coastal seawater Int J Syst Evol Microbiol 60 2071 - 2075    DOI : 10.1099/ijs.0.018051-0
Polymenakou P. N. , Bertilsson S. , Tselepides A. , Stephanou E. G. 2005 Bacterial community composition in different sediments from the Eastern Mediterranean Sea: a comparison of four 16S ribosomal DNA clone libraries Microb Ecol 50 447 - 462    DOI : 10.1007/s00248-005-0005-6
Polymenakou P. N. , Lampadariou N. , Mandalakis M. , Tselepides A. 2009 Phylogenetic diversity of sediment bacteria from the southern Cretan margin, Eastern Mediterranean Sea Syst Appl Microbiol 32 17 - 26    DOI : 10.1016/j.syapm.2008.09.006
Schloss P. D. , Handelsman J. 2005 Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness Appl Environ Microbiol 71 1501 - 1506    DOI : 10.1128/AEM.71.3.1501-1506.2005
Sokal D. M. , Michener C. D. 1958 A statistical method for evaluating systematic relationship. Univ. Kansas Sci Bull 38 1409 - 1438
Somerville C. C. , Knight I. T. , Straube W. L. , Colwell R. R. 1989 Simple, rapid method for direct isolation of nucleic acids from aquatic environments Appl Environ Microbiol 55 548 - 554
2012 DECIPHER, a search-based approach to chimera identification for 16S rRNA sequences Appl Environ Microbiol 78 717 - 725    DOI : 10.1128/AEM.06516-11