Bacterial Community Structure and Diversity of Closely Located Coastal Areas

Bacterial community structure and diversity of two closely located stations are usually considered similar which can be verified by more intensive investigations using relatively large amount of datasets from the next generation sequencer. This study was conducted to assess the bacterial community structure and diversity between two closely located coastal stations, the port side and the sea side of the Oarai, Ibaraki, Japan from March 2013 to July 2014 using 454 GS Junior sequencer. Two stations underwent similar changes in physicochemical properties but the community structure and diversity was different. The Proteobacteria (the class Alphaproteobacteria, followed by the Gammaproteobacteria) and the Bacteroidetes (the class Flavobacteriia) were two abundant phyla in both the stations. But, the Flavobacteriia was more abundant in the port side, contributed about 26% to 48%, compared to the sea side (about 12% to 39%). Conversely, the relative abundance of the Gammaproteobacteria was higher on the sea side, about 10% to 17%, compared to the port side (about 4% to 12%). Among others, the phyla Cyanobacteria, Deferribacteres, Verrucomicrobia and the class Betaproteobacteria were also relatively abundant at the sea side. Because of their dominancy, the class Flavobacteriia and Alphaproteobacteria were further analysed at a lower phylogenetic level and marked differences were observed between the stations. Bacterial biodiversity in terms of the species richness (Chao index) and evenness (inverse Simpson) indicated higher patterns of diversity in the sea side area compared to the port side. Non-metric Multidimensional Scaling fitting with the environmental features (metaMDS), redundancy analysis (RDA) and Bray-Curtis clustering analysis also showed marked differences in bacterial community structure and diversity between the stations. However, some OTUs were commonly found in both the stations in all the sampling periods. So, the bacterial community structure and diversity of the coastal areas are distinguishable even between two closely located sampling points. Corresponding author.


Introduction
The community structures of bacteria are considered as one of the most fundamental information in microbial ecology as it provides basic information regarding the environment.Bacterial communities are usually modified by many environmental conditions [1]- [3] and their spatiotemporal changes and biogeographical distributions are of special attention [4].However, in common practice bacterial community structures are considered similar between two closely located areas of an aquatic environment and treated them as replicate to one another.To verify this concept, more intensive investigations at a relatively finer scale are required which had been hampered due to methodological limitations, mainly the difficulties in culturing prokaryotic cells.Recent developments of molecular techniques, however, considerably overcome this problem by directly obtaining the genetic information without cultivation [5]- [8].Furthermore, the introduction of the next generation sequencer [NGS] made it possible to obtain by a far large amount of sequencing data within a short period of time and showed the presence of numerous previously unknown sequences or operational taxonomic units (OTUs).These facilities allowed us to assess any similarity or dissimilarity even between two closely located areas at different times of the year.
Microbial habitats are fluctuating widely in coastal environments because of the influences of terrestrial, freshwater and oceanic conditions.Some areas are also affected by anthropogenic activities.Organic matters, nutrients, pollutants and microorganism may be brought into coastal environments depending on the geographical characteristics, season, local weather, currents and so on.Coastal microbial communities consisting of highlyactive and diversified microbes have an important role in alleviating pollution and environmental damage due to nutritional supply from terrestrial sources [9].Also, the higher levels of bacterial diversity in the coastal estuarine habitats are considered to be causally related to the mixing of bacterial communities from different environments through the act of river influx and tidal exchange [10].A number of physicochemical parameters significantly influence the bacterial diversity of this brackish water habitat such as temperature, salinity, and dissolved nutrients [11].Although, the large populations of bacteria are well documented in coastal water research, their variations in terms of community structure and diversity between closer points were not considered well.
The purpose of this study is to assess the similarities or dissimilarities in bacterial community structure and diversity between two coastal areas of Oarai, Ibaraki, Japan at different time scale.The studied stations are located closely, only about one kilometer far from one another.One of them is the Oarai port area (port side), partly bounded by a sea bank, and the other is the Oarai beach area on the open seashore (sea side).We assumed that although the study areas are located closely, their community structure and biodiversity will be different.Because of more influences by both freshwater and marine water, bacterial biodiversity will be higher on the sea side station compare to the port side.As the stations are located closely, sea waters seem to be exchanged between them and thus, the basic physicochemical properties will be similar.

Sampling Collection
The seawater samples were collected from two sampling stations, the "port side" (36˚30'99"N and 140˚58'46"E) and "sea side" (36˚31'74"N and 140˚59'20"E), of the Oarai coastal area of the North Pacific Ocean, Ibaraki prefecture, Japan (Figure 1).Sampling stations are located closely, about 1 km far from one to another.One of them, the port side, is close to the Oarai port area which is semi-enclosed by a sea bank.The sea side sampling point, one the other hand, is thought to be subjected to the river inflow to some extent because, the Naka River, one of the class-one rivers of Japan, and the Hinuma River flows into the Pacific at the north end of the coastline.
In order to observe the similarities and dissimilarities in bacterial community structure and diversity at different

Sample Filtrations and Preparation
About two litres of the seawater sample were filtered through 0.22 µm pore sized Sterivex-GP pressure filter unit (Millipore, Billerica, MA, USA) using a peristaltic pump to have the community compositions.The sterivex cartridge filters were immediately kept in sterile bags and stored at −80˚C until further processing.The filters in the sterivex units were cut aseptically and placed inside screw tubes just before DNA extraction.

Environmental Parameters
The water temperature was measured at the time of sampling by using a mercurial thermometer.The seawater salinity was determined by a handy refractometer (IS/Mill-E, As One, ATAGO, Japan).The chlorophyll-a dataset with 4-km resolution was obtained from Level-3 MODerate resolution Imaging Spectroradiometer (MODIS) aqua standard-mapped image distributed by the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (oceancolor.gsfc.nasa.gov/cms).Averages of eight consecutive days used to prepare chlorophyll-a illustrated maps for each sampling time using Ocean Data View (ODV) software [12].

DNA Extraction and Cleaning
DNA was extracted from the sterivex filters by combined use of ChargeSwitch Forensic DNA Purification Kits (Invitrogen TM , Carlsbad, USA) and ZircoPrep Mini (FastGene TM , Nippon Genetics Co. Ltd., Bunkyo-ku, Tokyo, Japan) beads beating with a slight modification of the manufacturer's protocol.A MicroSmash (MS-100R, Tomy Medico., Ltd., Tokyo, Japan) was used for beads beating at 5000 rpm and 4˚C for 30 seconds for each filter under sterile conditions with a great care to avoid contamination.The extracted DNA was also cleaned using NucleoSpin (MACHEREY-NAGEL GmbH & Co. KG, Neumann-Neander-Str., Düren, Germany) gDNA cleanup kit according to the manufacturer protocol and stored −30˚C until amplification.

16S rDNA Amplification and Pyrosequencing
The V1-V3 hyper variable regions of 16S rDNA gene were amplified from the extracted DNA templates by polymerase chain reaction (PCR).The 27F with multiplex identifier (MID): 5'-CCATCTCATCCCTGCG-TGTCTCCGACTCAGXXXXXXXXXXAGAGTTTGATCMTGGCTCAG-3', where X's represents the sample-specific multiplex identifier-MID [13] was used as the forward primer and the 519R with adaptor: 5'-CCTATCCCCTGTGTG-CCTTGGCAGTCTCAG(GWATTACCGCGGCKGCTG)-3' was used as the reverse primer.Each PCR reactions were carried out in a volume of 20 µL in triplicates while the mixture consisted of 2 µL DNA template, 13.1 µL molecular biological grade double distilled water, 0.6 µL (5 µM) each primer, 2 µL 10X TaKaRa Ex Taq Buffer, 1.6 µL TaKaRa dNTP mixture, and 0.1 µL TaKaRa Ex Taq HS Polymerase (TaKaRa, Japan).Thermal cycling was carried out for a total of 25 cycles as per the following conditions: initial denaturation at 94˚C for 4 mins, denaturation at 98˚C for 10 sec, annealing at 55˚C for 30 sec, elongation at 72˚C for 1 min and final elongation at 72˚C for 10 mins.After amplification, the desired length of the 16S rDNA gene was confirmed by agarose gel electrophoresis and any sort of contamination was carefully verified by observing the bands of the triplicates of the same samples.After confirming the desired length, amplified DNA products were purified and normalized using Agencourt AMPure XP (Beckman Coulter INC., USA) according to the guidance of the 454 Sequencing Amplicon Library Preparation Method Manual (GS Junior Titanium Series 2012, Roche, USA).The purified DNA amplicon was then quantified using Quant-iT Picogreen dsDNA Kit (Invitrogen, Carlsbad, USA).The bacterial 16S rDNA gene amplicons were then sequenced using the 454 GS Junior sequencer (Roche, USA) at Atmosphere and Ocean Research Institute (AORI), the University of Tokyo (Kashiwa, Chiba, Japan) according to the manufacturer's protocol for 454 GS Junior Titanium Series.

Sequence Analyses
The open-sourced MOTHUR program [14] was used for subsequent analysis, quality checking and arrangement of the obtained sequences following the guidelines available to the operation manual for the 454 [http://www.mothur.org/wiki/454_SOP].Initially, the unique sequences were selected and then the similar sequences were clustered and aligned against the SILVA bacterial databases [15].Then the pre-cluster method [16] was applied to reduce the sequencing errors by screening, filtering, and de-noising.The chimera.uchime command was used for checking and removing the chimeras.The sequences were subsequently classified against the ribosomal 160 database project (RDP) database and the inactive components such as chloroplast, mitochondria etc. organelles affiliated "former" bacterial sequences were removed from our dataset to improve the data quality.The qualified high-quality sequences were then used to generate distance matrix and clustered assigning to operational taxonomic units (OTUs) at 97% identity level [17].A representative sequence from every OTU was used for classification by running the MOTHUR program based on the SILVA bacterial databases.To standardize the number of sequences between samples, they were randomly re-sampled to the sample with the fewest reads (2674 reads) using the MOTHUR program based on the OTU files clustered at 0.03 cut-off levels.
The species richness and diversity indices were considered to evaluate the biodiversity and analyse the rarefaction.For the species richness, the Chao1 index [18], and for diversity, the inverse Simpson (Invsimpson) index [19] was calculated using the MOTHUR software at OTU definition at a distance of 0.03.

Statistical Analyses for Community Structure
To check the correlations between bacterial communities and environmental factors, nonmetric multidimensional scaling fitting with the environmental features (metaMDS) was carried out based on the relative abundance data of each OTUs.The permutation test was used following the "MASS" [20] and "Vegan" package [21] from R software (R Development Core Team 2012).Redundancy analysis (RDA) was also carried out using the R software with "Vegan" package [21] based on the relative abundance data of each OTUs and environmental information.Toassess the similarities or dissimilarities between the bacterial groups of the two sampling stations, the clustering analysis (Bray-Curtis) test [22] was also performed using the R software with "Vegan" package [21].

Environmental Parameters
Changes in the water temperature, salinity, and chlorophyll-a values at the two sampling stations are shown in Table 1.Salinity showed similar fluctuation patterns at both stations, with the maximum values, around 35, were obtained in February 2014 and the minimum, around 24, in July 2014 after the rainy season for both the stations.Other features, such as water temperature and water depth were also similar at both the stations.
The chlorophyll-a values obtained from the satellite data were used to prepare chlorophyll-a illustrated maps (Figure 2) for each date of sampling using Ocean Data View (ODV) software.Then the chlorophyll-a values of the sampling locations were obtained from this illustrated maps considering the location (latitude and longitude) of the sampling stations.However, as the sampling stations are closely located, it was not possible to obtain data separately for the two stations (Table 1).The chlorophyll-a values varied from about 2.5 µg•m −3 in February to about 5.0µg•m −3 in March.No chlorophyll-a data was obtained during July due to cloudy weather (Table 1 and Figure 2).

Bacterial Community Structure Analysis
After sequencing all the samples a total of 66,609 sequences were obtained which consists of 5249 different types of OTUs.The obtained sequences were analysed for community composition.The composition of the Table 1.Description of the environmental characteristics of the sampling sites throughout the study periods.The water temperatures and salinity were determined by using a mercurial thermometer and a handy refractometer respectively.Chlorophyll-a data were obtained from the satellite data of NASA's Ocean Color website.major groups mostly at phylum/class level is presented to Figure 3.

Samplings
The Phylum Proteobacteria was the most dominant one followed by the Bacteroidetes at almost all the sampling periods regardless of the sampling stations.Among different subgroups of Proteobacteria, the class Alphaproteobacteria was the most abundant one and the class Gammaproteobacteria was the second.In the case of the phylum Bacteroidetes, the class Flavobacteriia was mostly abundant.The relative abundance of the Flavobacteriia was higher in port side while the Gammaproteobacteria in the sea side at almost all the sampling periods.The class Alphaproteobacteria was almost equal or slightly higher in the port side station except April 2014.Among other groups the phylum Cyanobacteria, Verrucomicrobia and the class Betaproteobacteria were higher in abundance to the sea side as compared to the port side.Moreover, the unclassified members were also higher on the sea side.The class Flavobacteriia was most abundant in March, followed by the February at the port side, while the class Alphaproteobacteria in July at both the sampling points.The abundance of the phylum Cyano- bacteria was also relatively higher in July and lower or almost absent in March and February; while the phylum Actinobacteria was more abundant in February.The phylum Deferribacteres was found abundant in October at sea side station (Figure 3).
As the class Flavobacteriia and Alphaproteobacteria were 2 most dominant classes in both the stations throughout study periods, community composition within these classes at family or genus level was also evaluated.In the case of Flavobacteriia, the order Flavobacteriales was the only contributing group (Figure 4).Although, there were differences between the sampling periods, but in general, genera NS3a marine group, Polaribacter, and Winogradskyella of the family Flavobacteriaceae comprised most of the bacterial fractions at port side while the genera Flavobacterium, NS4 marine groups and NS5marine groups of the same family at the sea side station (Figure 4).Analysis of the members within Alphaproteobacteria also showed marked differences between the stations (Figure 5).The maximum contribution to the bacterial community was made by the members of the order Rhodobacterales.The genera Lentibacter, Nereida, Sulfitobacter and unclassified members of the family Rhodobacteraceae was relatively abundant in the port side station while the genus Roseobacter clade, order Rickettsiales and SAR11 contributed significantly in the sea side station.However, there were seasonal variations in their abundance.For example, the genus Sulfitobacter contributed mostly to the samples of March, February, and April while the order SAR11 to the samples of October, April and July (Figure 5).

Diversity of Bacterial Communities
Bacterial biodiversity was evaluated in terms of the species richness and richness-evenness considering the Chao index and inverse Simpson (invSimpson) index, respectively, as was shown in Figure 6.The Chao index values were higher at the sea side in March 2013, October 2013 and July 2014 while at the port side in February 2014 and April 2014.Compared to the port side in almost all the sampling periods except April 2014, Simpson index values were higher at the sea side.It seems reasonable, from these observations, to suppose that more diversified communities were made up in the sea side area.For both the indices, the highest values were observed in April and the lowest in October on the port side.Relatively higher Chao index values were observed in October and July on the sea side, while the lowest was shown in February (Figure 6).The rarefaction curves are showing the relationship between the numbers of obtained sequences and observed OTUs (Supplementary Figure 1A).

Seasonal Environmental Changes and Bacterial Community Structures of the Locations
Non-metric Multidimensional Scaling fitting with the environmental features (metaMDS) based on the relative abundance data of the samples and seasonal environmental data was used to categorize bacterial community composition of the studied stations at different sampling periods.The samples were separated according to the sampling periods (r 2 = 0.63 and P = 0.1, based on 1000 permutations) as well as according to the stations (r 2 = 0.22 and P = 0.17, based on 1000 permutations); the community composition was similar between the stations at February (Figure 7).The association between the bacterial community structure and environmental factors was examined by RDA.In Figure 8, the samples are plotted on a sample-to-sample basis, with respect to stations (r 2 = 0.3 and P = 0.02, based on 1000 permutations) and season (r 2 = 0.5 and P = 0.4, based on 1000 permutations).The water temperature (r 2 = 0.08 and P = 0.7, based on 1000 permutations) and salinity (r 2 = 0.1 and P = 0.6, based on 1000 permutations) had no significant influence in clustering (Figure 8).
We performed a Bray-Curtis clustering analysis based on the bacterial relative abundance data at phylogenetic level.Samples were aligned according to degrees of similarity in community composition on a sample-tosample basis (Figure 9).However, the April sample of the port side station was aligned next to the sea side clade, whilst the February sample of the sea side was aligned next to the port side clade.This indicates that there were similarities in bacterial community composition among the samples of the February and April (Figure 9).

Analyses of the Common OTUs across Sampling Stations and Seasons
About 20.92% to 30.72% of the OTUs were common in the port side station and about 7.25% to 31.53% in the  sea side station between different sampling periods (Table 2).The percentages of common OTUs shared between March and October, October and February, February and April, and April and July were 6.95, 6.88, 10.20, and 5.79 respectively at the port side station while 7.63, 5.12, 8.25, and 7.90 respectively for the sea side station (Table 3).The overall observations showed that bacterial community composition was fluctuated highly between two consecutive sampling periods and the community composition was also dissimilar between the stations at most of the sampling periods (Table 2 and Table 3).Table 2. Percentage of common OTUs between the port side and sea side station in different periods of sampling.The percentage was calculated after counting the total number of OTUs appeared at each station as well as among the stations for every sampling period.

Mar-2013
Oct Table 3. Percentage of common OTUs between two consecutive sampling periods at port side and sea side station in order to evaluate the degree of fluctuations between two sampling periods at OTUs level.The percentage was calculated after counting the total number of OTUs appeared at each sampling periods as well as among two consecutive sampling periods.Only 14 OTUs out of 5249 (0.27%) was found common throughout the study regardless of the sampling periods and stations.However, their contribution to the total abundance was 10.7% in July to 27.36% in April on the port side station while 24.46% in October to 48.31% in March on the sea side station (Figure 10).Bacterial groups of these common OTUs were Candidatus Actinomarina (Acidimicrobiales) of phylum Actinobacteria; Fluviicola, NS4 marine group, NS5 marine group (2 OTUs), Owenweeksia and Polaribacter (2 OTUs) (Flavobacteriales) of phylum Bacteroidetes; Roseobacter clade (2 OTUs) (Rhodobacteriales) and unclassified SAR11 of Alphaproteobacteria; unclassified Alteromonadales and SAR86 (Oceanospirillales) of Gammaproteobacteria; and 12up (Rhodocyclales) of Betaproteobacteria (Figure 10).

Discussions
The introduction and use of the next generation sequencer [NGS] made it possible to obtain numerous previously unknown sequences or operational taxonomic units (OTUs) that can be used to verify and assess the bacterial community structure at a relatively finer and intensive scale.This study was conducted to evaluate the similarities or dissimilarities in bacterial community structure and diversity between two closely located coastal areas of Oarai, Ibaraki, Japan at different time scale.Bacterial community structure was retrieved by obtaining high-throughput sequencing data using Roche 454 sequencer.The results indicated that two sampling stations underwent a similar change in physicochemical properties but the community structure and diversity was dissimilar between the stations.The class Alphaproteobacteria followed by the class Gammaproteobacteria of the phylum Proteobacteria and the class Flavobacteriia of the phylum Bacteroidetes were mostly abundant but the relative abundance of Flavobacteriia was higher atthe port side and Gammaproteobacteria at the sea side throughout the study period.Among others, the phyla Cyanobacteria, Deferribacteres, Verrucomicrobia and the class Betaproteobacteria were also relatively abundant at the sea side.It was found that the relative abundance of different bacterial groups was fluctuated markedly across time due to seasonal influences and there also were marked differences between the stations at almost all the sampling periods.Bacterial biodiversity in terms of the species richness (Chao index) and evenness (inverse Simpson) indicated high levels and patterns of diversity in the sea side area compared to those in the port side.Non-metric Multidimensional Scaling fitting with the environmental features (metaMDS), RDA and Bray-Curtis clustering analysis also showed marked differences in the bacterial community structure and biodiversity between the sampling stations.However, some OTUs were commonly found in both the stations in all the sampling periods.
Previous studies suggested that the phylum Proteobacteria was found everywhere as an abundant one.Among different classes of the Proteobacteria, the classes Alphaproteobacteria (mostly SAR11), followed by the Gammaproteobacteria are abundant in marine waters [23]- [26]; while the class Betaproteobacteria is abundant in freshwater habitats [8] [27]- [29].The phylum Bacteroidetes is also dominant in some freshwater and marine water habitats [27] [29]- [31].Our results are consistent with those published findings.The higher relative abundance of the phylum Verrucomicrobia at the sea side point may be because they predominate in shallow brackish wasters compared to open oceans [32].The higher abundance of the phylum Cyanobacteria of July's sampling agreed with the concept that the growth rate of Cyanobacteria is usually higher at high water temperature during the summer season [33]- [35].The relative abundances of the class Betaproteobacteria were higher at the sea side station, may be because this station is influenced by riverine communities (and so terrestrial communities as well) of the Naka River [36] [37] (Figure 1).
Further analyses of the community composition within the classes Flavobacteriia and Alphaproteobacteria showed marked differences between the studied stations.Within the class Flavobacteriia, the genera NS3a marine group, Polaribacter, and Winogradskyella of the family Flavobacteriaceae comprised most of the bacterial fractions on port side while the genera Flavobacterium, NS4 marine groups and NS5marine groups of the same family at the sea side station (Figure 4).Korlević et al. 2015 [35] reported that pyrotags related to the order Flavobacteriales were abundant, with high frequencies of clades NS2b, NS4, and NS5, which is consistent with our results.The genus Polaribacter was first isolated from a polar marine environment [38], however, it was also isolated from coastal areas of Japan [39] [40].The genus Winogradskyella was isolated from an alga collected from the Sea of Japan [41] [42] while the members of genus Flavobacterium are widely distributed [43].Among different members within Alphaproteobacteria, the genera Lentibacter, Nereida, Sulfitobacter and unclassified members of the family Rhodobacteraceae was relatively abundant in the port side station while the Roseobacter clade, order Rickettsiales and SAR11 in the sea side station (Figure 5).Previous reports showed that the marine Rhodobacterales is widespread the members of marine Roseobacter clade formed the most common and dominant surface-colonizing bacterial group [44] [45].The members of the genus Sulfitobacter are widely distributed in coastal and open ocean environments [46] [47], where they may play an important role in organic sulfur cycling.The type species of the genus Lentibacter was isolated from seawater samples in the coastal region of Qingdao, China, during a massive green algae bloom [48], indicating that this genus has an affinity to eutrophic waters that are present at port side, observed during samplings.Information on the genus Nereidain marine environments is very limited [49], but theSAR11 is a typical dominant group in the oceanic surface environment among other orders of Alphaproteobacteria [25] [29] [50]- [52].Previous reports indicated that these available members of the class Flavobacteriia and Alphaproteobacteria are common in the coastal marine environments, consistent with our findings.
The bacterial diversity at the sea side station was higher than that of the port side station; and within port side station, diversity was higher in April while within sea side station in October (Figure 6).Cury et al. 2011 [53] and Rodrigues et al. 2013 [54] reported that anthropogenic activities negatively influence the bacterial diversity in forest soil and coastal environments respectively, explaining the reason why the bacterial diversity at the port side was less than that of the sea side.Less pollution and better water quality i.e. better ecological condition of the coastal environment also supports higher bacterial diversity [55].Moreover, introduction and mixing of an exogenous bacterial group from river input may also affect the bacterial diversity in the sea side samples.A number of physical (especially temperature), chemical (salinity, nutrients, oxygen concentration, pollution, etc.), and biological (predation, competition, plankton bloom) factors influence the bacterial diversity in coastal areas [11] [34] [56] [57].Further consideration of specific environmental factors and investigation of their seasonal changes are required to explain the higher diversity in April at port side and in October at the sea side.

Conclusion
In conclusion, bacterial community structures and diversity were investigated by seasonal samplings at two closely located coastal stations.Although the community was mostly dominated by the phyla Proteobacteria and Bacteroidetes, there were variations in their relative abundance between the sampling stations and the periods of samplings.The overall observations also indicated that the bacterial communities in the sea side area grow in diversity compared to that in the port side area.So, bacterial community structure and diversity of the coastal

Figure 1 .
Figure 1.Sampling stations, the port side and the sea side; the inset map area filled with red color is the Ibaraki prefecture and the arrow indicating the position of Oarai.time scale, water samples were collected in March 2013 (early spring), October 2013 (autumn), February 2014 (winter), April 2014 (spring), and July 2014 (summer).At every sampling, about 5 liters of seawaters were collected in a previously sterilized screw-capped plastic bag and carried back to the laboratory in ice boxes within 2 hours after sampling.

Figure 2 .
Figure 2. Chlorophyll-a data from satellite at different dates of sampling (8 days average, prepared by using Ocean Data View software).Some areas in the figures kept blank (white) because of unavailability of data due to clouds.

Figure 3 .
Figure 3. Bacterial community structure of the studied stations at different sampling periods deduced from the 16S rRNA pyrosequencing analysis.The groups "Others" referred to the sum of those phyla did not individually contributed 1% of the relative abundance in at least one sample and "Unclassified" are the unidentified/unknown members.

Figure 4 .
Figure 4. Relative abundance of different members within the class Flavobacteriia showing differences in abundance and composition between the port side and sea side stations at different sampling periods.

Figure 5 .
Figure 5. Relative abundance of different members within the class Alphaproteobacteria showing differences in abundance and composition between the port side and sea side stations at different sampling periods.

Figure 6 .
Figure 6.Changes in (a): Chao index and (b): inverse-Simpson index of the studied stations indicating the biodiversity in terms of richness and richness-evenness at various sampling periods.

Figure 7 .
Figure 7. Non-metric Multidimensional Scaling fitting with the environmental features (metaMDS) showing clustering according to the sampling periods and stations.The first and the second part of the sample IDs' are expressing the sampling periods (Ma = March, Oc = October, Fe = February, Ap = April and Ju = July), and the sampling stations (PS = port side, SS = sea side) respectively.Clustering of two samples closely meaning they are relatively similar in composition.

Figure 8 .
Figure 8. Redundancy analysis (RDA) of the samples showed the clustering of the samples.The abbreviations are same as Figure 7.

Figure 9 .
Figure 9. Clustering analysis (Bray-Curtis) showing similarities or dissimilarities in bacterial community composition between the stations.The abbreviations are same as Figure 7.Samples are arranged into two major clades according to the stations.However, the April sample of the port side arranged with the April sample of the sea side in sea side clade while February sample of the sea side arranged with the February sample of the port side in port side clade indicating their compositional similarities.

Figure 10 .
Figure 10.Relative abundance of the bacterial groups of the common OTUs, those were found common at the entire study period regardless of the stations and sampling periods.

Sampling Sites Water Temp. (˚C) Salinity * Water depth Chlorophyll-a (µg•m −3 ) ** Closest location with data Obtained values
* PSU, practical salinity unit.** As the sampling stations are located closely, just about 1 km far from one another, it was not possible to get chlorophyll-a data separately for port side and sea side station rather a single value from the closest available location was noted here for both the stations.