Bacterial diversity in relatively pristine and anthropogenically-influenced mangrove ecosystems (Goa, India).

To appreciate differences in benthic bacterial community composition at the relatively pristine Tuvem and the anthropogenically-influenced Divar mangrove ecosystems in Goa, India, parallel tag sequencing of the V6 region of 16S rDNA was carried out. We hypothesize that availability of extraneously-derived anthropogenic substrates could act as a stimulatant but not a deterrent to promote higher bacterial diversity at Divar. Our observations revealed that the phylum Proteobacteria was dominant at both locations comprising 43-46% of total tags. The Tuvem ecosystem was characterized by an abundance of members belonging to the class Deltaproteobacteria (21%), ~ 2100 phylotypes and 1561 operational taxonomic units (OTUs) sharing > 97% similarity. At Divar, the Gammaproteobacteria were ~ 2× higher (17%) than at Tuvem. A more diverse bacterial community with > 3300 phylotypes and > 2000 OTUs mostly belonging to Gammaproteobacteria and a significantly higher DNT (n = 9, p < 0.001, df = 1) were recorded at Divar. These findings suggest that the quantity and quality of pollutants at Divar are perhaps still at a level to maintain high diversity. Using this technique we could show higher diversity at Divar with the possibility of Gammaproteobacteria contributing to modulating excess nitrate.

Introduction

Mangroves constitute nearly 75% of tidal vegetation in tropical regions (Alongi ). The proximity of these coastal wetlands to human inhabitation, aquaculture farms, waste discharge from industrial units and domestic sewage release make them vulnerable to high nutrient inputs especially inorganic nitrogen (N). Microbial communities in mangrove sediments play an important role in the biogeochemical cycles of coastal ecosystems (Thorsten and José, 2001). Anaerobic conditions and substrate availability in organic C-rich mangrove sediments (Krishnan and LokaBharathi, 2009) favour alternate respiratory pathways like denitrification and sulfate reduction. Recently, it has been shown that bacterially-mediated benthic denitrification effectively reduces N load in an anthropogenically-influenced tropical mangrove ecosystem of Goa, India (Fernandes and LokaBharathi, 2011). However, little is known about the microbial communities inhabiting these regions.

In-depth analysis of prokaryotic communities is crucial in understanding ecosystem functioning. So far, few culture-independent studies have been carried out to assess the taxonomic diversity of bacteria in mangrove sediments (Liang ; Marcial Gomes ; Sjöling ). The sequences recovered in these studies revealed the dominance of Proteobacteria and the presence of all five sub-divisions within this class. As these phylogenetic studies used traditional clone library and sequencing approaches, only a small fraction of the mangrove bacterial communities could be sampled. New technology, such as pyrosequencing of hypervariable regions of the 16S rRNA gene, is a cost-effective and a better alternative to examine the phylogenetic diversity of microbial populations (Huse ; Schwarz ). This technology has been applied to the water column of the oceans (Andersson ; Huber ; Kirchman ; Sogin ), soils (Roesch ) and more recently to marine sediments (Chariton ; Kim ). To our knowledge, very few studies have used the pyrosequencing approach to examine the prokaryotic diversity in mangrove sediments.

Though physiology cannot be assumed from phylogenetic analysis (Kim ), bacterial diversity data could be useful for understanding the ecological roles that microbes may play in the environment. Hence, the present study was carried out to gain an insight into the taxonomic diversity of bacteria in two tropical mangrove systems of Goa, India using the pyrosequencing approach. The relatively pristine site Tuvem was compared to Divar which is influenced by extraneous nutrient input. The Divar mangrove ecosystem lies along the Mandovi estuary. This estuary receives a considerable supply of inorganic nitrogenous nutrients from various sources such as mining wastes (De Souza, 1999), land runoff (Sardessai and Sundar, 2007) and sewage effluents (Ansari ). We hypothesize that these extraneously-derived substrates at Divar could support higher bacterial diversity. Our observations revealed higher bacterial diversity at Divar suggesting that the levels of pollutants have not reached a limiting level. Higher sequence clusters of class Gammaproteobacteria and denitrification activity (DNT) were recorded at Divar which is suggestive of their important ecological role in controlling the levels of excess nitrate in the ecosystem through denitrification.

Materials and Methods

Study area and sampling

Investigations were carried out at mangrove forests located along the Mandovi and Chapora rivers in Goa, west coast of India (Fernandes ). One site was at Divar island (15°30′35″ N and 73°52′63″ E) which is separated from the mainland by the Mandovi estuary. The adjoining Mandovi estuary is used heavily for transportation of iron ore from mines located upstream. These iron ore beneficiation plants situated on the riverbank, discharge effluents directly into the estuary. This discharge contains high quantities (~ 8–10 Mg month−1) of NH4NO3 used in ferromanganese mining operations (De Souza, 1999). The second site was located at Tuvem (15°39′94″ N and 73°47′65″ E) along the river Chapora. This ecosystem is comparatively less influenced by anthropogenic activities and considered as relatively pristine (Krishnan ).

The Mandovi estuary receives high extraneous nitrate input during the monsoon (Divya ). Sampling was carried out at low tide during this season (September, 2008). A total of 12 cores were collected from each site which were used as follows:

Three cores were collected using PVC push cores (7.5 cm inner diameter, 20 cm length) and capped at both ends. They were transferred to a plastic tray for transportation to the laboratory under ambient temperature. The variation in physical parameters was measured on arrival to the laboratory i.e. within 1 h of sampling.

Six cores of dimensions described above were collected for chemical and microbiological analyses. These were transported to the laboratory in an ice box. On arrival, the cores were maintained at 4 °C until analysis.

Three sediment cores were collected for molecular biology analysis using pre-cut (open-ended), pre-sterilized 50 mL syringe cores. The contents were transferred to clean plastic bags, each of which was sealed with a rubber band. The cores were immediately preserved in a liquid nitrogen containing cryocan for transportation to the laboratory. On arrival, the samples were stored at −20 °C until analysis.

Earlier observations at Tuvem and Divar have shown have shown highest DNT to occur at the surface (0–2 cm) of mangrove sediments (Fernandes ; Fernandes and LokaBharathi, 2011). Hence, analyses of all parameters in the present manuscript have been restricted to the surficial sediments.

Physico-chemical analysis

Hydrogen ion concentration (pH) in the surficial sediments was measured using an Orion 4-Star Plus benchtop pH/ISE meter (Thermo Fisher Scientific Inc., Waltham, MA). Sediment oxidation-reduction potential (Eh) was measured using an Orion platinum redox in combination with an Ag/AgCl2 reference electrode (Thermo Fisher Scientific Inc.).

For nutrient analyses, the top two centimeters of each core was sectioned and transferred to 100 mL of sterile saline (8.5 g L−1 NaCl). The contents were gently homogenized using a glass rod. The resultant slurry was centrifuged at 4 °C for 10 min at 5000 rpm. A low spin speed was maintained during centrifugation to ensure minimal change in nutrient concentrations due to lysis of benthic infauna. The supernatant was filtered through a 0.2 μm filter. Nitrite, ammonium and nitrate were measured colorimetrically (± 0.01 μmol L−1) according to standard methods (Bendschneider and Robinson, 1952; Koroleff, 1969; Wood ). Weight of the sediment used in the extraction was estimated by drying the sediment in a hot air oven at 60 °C for 48 h.

Surficial sediment samples for analysis of total organic carbon (TOC) were dried in a hot air oven at 60 °C for 48 h. The TOC in the sediment was determined by the wet oxidation method (El Wakeel and Riley, 1957) with a precision of 0.01%.

Total bacterial counts

Sediment cores were sectioned from surface to 2 cm. Sub-samples of ≈ 5 g wet weight sediment were extruded using syringe cores. The sub-samples were transferred to 45 mL of filter sterilized full strength seawater (10−1 dilution). Tween 80 (50 μL) was added and the mixture was sonicated at 40 mHz for 15 s. The sub-samples were diluted with filter sterilized seawater and fixed with buffered formalin (2% final concentration).

Enumeration of total bacteria was carried out by epifluorescence microscopy (Hobbie ). A sub-sample was stained with acridine orange (final concentration 0.01% w/v), incubated for 2 min and then filtered through 0.2 μm black Isopore polycarbonate filters (Millipore). Bacterial cells were counted using a Nikon 50i epifluorescence microscope equipped with a 100X oil immersion objective. Cells were counted from 10–100 microscopic fields. Total bacterial counts were expressed as cells g−1 (dry sediment).

Denitrification activity

Denitrification rates were measured using sediment slurries by the acetylene inhibition technique based on the inhibition of the conversion of N2O to N2 (Sørensen, 1978). About 1 cm3 sediment from 0–2 cm was extruded using a syringe core and transferred aseptically to sterile 20 mL headspace vials. Three mL of sterilized ambient sea water from the sampling site containing 4.3 μmol NO3−-N L−1 was added to the sediment. No additional carbon or nitrate was amended to the seawater. The vials were capped with butyl stoppers and the slurry was briefly vortexed for 5 sec. The vials were purged with high purity N2 for 10 min to induce anaerobic conditions. The headspace over these slurries was amended with acetylene at 20 kPa (Bonin ) and the tubes were briefly vortexed. Triplicates were maintained and the vials were then incubated in the dark for 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 h. At the end of the incubation period, bacterial activity in all incubations was terminated using 0.1 mL of 1 mol L−1 HgCl2 solution.

Nitrous oxide in the headspace was analyzed using a Shimadzu 2010 gas chromatograph fitted with a electron capture detector and Porapak Q column (1/8″ SS column, 3.05 m length, 80/100 mesh). The oven and detector temperatures were 40 °C and 300 °C, respectively. High purity nitrogen at a flow rate of 35 mL min−1 was used as a carrier gas. The gas chromatograph was calibrated using a secondary standard 44 ± 0.38 nmol L−1 N2O in nitrogen (National Physical Laboratory, New Delhi). Denitrification activity was calculated based on the linear accumulation of N2O over time (Tiedje, 1982). The solubility coefficients of N2O were used to correct for dissolved N2O (Weiss and Price, 1980). Sediment used in each microcosm was filtered through a laboratory grade filter paper and dried at 60 °C for approximately 48 h and weighed. Denitrification rate was expressed as nmol N2O g−1 h−1 of dry sediment. Analysis of variance (ANOVA) in Analysis tool pack (Microsoft Excel) was used to assess significant difference in denitrification activity between both the locations.

DNA extraction, high-throughput pyrosequencing and data analyses

Genomic DNA was extracted from 0.25 g (dry weight) sediment from each of the 3 cores collected. The Mo Bio PowerSoil DNA isolation kit (Mo Bio Laboratories, California) was used to extract the DNA according to the manufacturer’s protocol. Though analyses in triplicate would have been ideal, analysis was restricted to one in this study as it was a preliminary assessment to draw a first order comparison. The DNA of single sample was concentrated and used for pyrosequencing analysis while the rest were used for other down-stream applications. The hypervariable region of 16S rRNA gene (BV6-rRNA tags) were amplified, subjected to high-throughput pyrosequencing and the sequence reads were then trimmed (Sogin ; Huber ). Pyrosequencing error rates were minimized (Huse ).

The 454 tags (reads) served as query to identify its closest match in a reference database (V6RefDB) containing ≈ 40,000 unique V6 sequences (Sogin ). Singletons were removed from the dataset before further analyses. Sequence characteristics like average length and tag aggregates were estimated using R package (R Development Core Team, 2007). Multiple sequence alignment was done using ClustalX 1.83. Distance matrices were calculated using DNAdist from PHYLIP 3.69 (Felsenstein, 2005). DOTUR (Schloss and Handelsman, 2005) was used for clustering tags that ranged from unique sequences showing no variation to 5% dissimilarity. These clusters served as operational taxonomic units (OTUs) for generating rarefaction curves and for calculating two indices of diversity- Chao1 (Chao, 1987) and ACE (abundance-based coverage estimator (Chao and Lee, 1992). The OTU overlap between samples was plotted using the Venn Diagram Plotter (Pacific Northwest National Laboratory, Department of Energy). Information on primers is available at the VAMPS site (http://vamps.mbl.edu/). Tag sequences used in this study have been archived under ICM_ICR_Bv6-ICoMM Cooperative Run and can be accessed using the link http://vampsarchive.mbl.edu/vamps/combined/combined.php. The dataset ID no. for sample from Divar is ICR_0003_2008_09_05 while for Tuvem is ICR_0004_2008_09_05.

Results

Sediment characterization

At the time of sampling, precipitation due to the south-west monsoon resulted in fresh water conditions at the sampling locations as observed from salinity values of ≤ 1 psu (Table 1). Though pH was marginally higher at Tuvem, reducing conditions were more intense at this location with up to −202 ± 42 mV as compared to a redox potential of 55.5 ± 2.3 mV at Divar (Table 1). However, inorganic nitrogenous compounds (NH4+, NO3− and NO2−) and organic carbon content were higher at Divar than at the relatively pristine Tuvem. Denitrification activity at Divar occurred at a rate of up to 3.17 (± 0.99) nmol N2O g−1 h−1 which was significantly higher (n = 9, p < 0.001, df = 1) as compared to Tuvem. The mangrove sediments harboured up to 1010 bacterial cells g−1; the Divar sediments had marginally higher abundance (2.67 (± 1.76) × 1010 cells g−1) than at Tuvem (1.45 (± 0.20) × 1010 cells g−1).

Table 1

Environmental parameters at sampling locations.

	Tuvem	Divar
Temp (°C)	27.7 (± 0.1)	28.0 (± 0.1)
Salinity (psu)	0	1
pH	6.65 (± 0.01)	6.29 (± 0.03)
Eh (mV)	−202.2 (± 41.86)	55.5 (± 2.26)
NH4+-N (μmol L−1)	2.58 (± 0.69)	11.78 (± 0.44)
NO2−-N (μmol L−1)	3.37 (± 0.20)	5.85 (± 0.28)
NO3−-N (μmol L−1)	3.01 (± 0.11)	5.52 (± 0.10)
TOC (%)	1.05 (± 0.44)	1.98 (± 0.96)
Total bacterial cells (x1010 cells g−1)	1.45 (± 0.20)	2.67 (± 1.76)
Denitrification activity (nmol N2O g−1 h−1)	1.95 (± 0.85)	3.17 (± 0.99)

Bacterial diversity in mangrove sediments

More than 34,000 bacterial V6 amplicons were sequenced from the two sites (Table 2). The average read length of the sequences from both sites was 60.7 ± 2.6 bp. Taxonomic analysis revealed differences in community composition between the two sites. The phylum Proteobacteria was the most dominant, contributing 43 and 46% of total V6 tags at Tuvem and Divar respectively (Figure 1). About 22 and 18% of the bacterial tag sequences at Tuvem and Divar respectively have not been identified to the class level (Figure 2). Deltaproteobacteria was the next most abundant class at Tuvem comprising ~ 21% of the total tag sequences while at Divar their abundance was 15% (Figure 2). Other important classes recorded at the locations were Actinobacteria, Gammaproteobacteria and Alphaproteobacteria. Within class Deltaproteobacteria, up to 53% tags comprised of order Desulfobacterales (Figure 3). Some of the other deltaproteobacterial sequences found at these locations belonged to the orders Desulfuromonadales, Myxococcales and Synthrophobacterales. At Divar, the Gammaproteobacteria comprised 36% of the Proteobacteria followed by Deltaproteobacteria (33%) and Alphaproteobacteria (22%). Among the Gammaproteobacteria, > 60% of the tags remained unidentified to the order level (Figure 4). The next most abundant order was Chromatiales. About 76% of the tags at both the sites could not been identified to the genus level.

Table 2

Sequencing information and diversity estimates for all bacteria at Tuvem and Divar.

	Tuvem	Divar
DNA recovered from 0.25 g sediment (ng/μL)	10.8	15.2
Total number of bacterial V6 tag sequences	12954	21433
Total number of trimmed bacterial V6 tag sequences*	8030	15326
Total number of unique sequences	2661	3707
Total OTUs at 3% difference	1561	2198
Chao1 estimator of richness at 3% difference (min, max)	1988 (1914, 2076)	3092 (2968, 3237)
ACE estimator of richness at 3% difference (min, max)	2166 (2136, 2196)	3342 (3261, 3429)

After removal of singletons.

Figure 1

Phylum level distribution of tag sequences at Tuvem (n = 8030 tags) and Divar (n = 15326 tags).

Figure 2

Class level phylogenetic classification of bacterial communities at Tuvem (n = 8030 tags) and Divar (n = 15326 tags). Tags of known classes that occurred at an abundance of < 0.1% have been grouped in the “Other known” category.

Figure 3

Order level phylogenetic classification of Deltaproteobacteria at Tuvem (n = 1712 deltaproteobacterial tags) and Divar (n = 2357).

Figure 4

Order level phylogenetic classification of Gammaproteobacteria at Tuvem (n = 762 gammaproteobacterial tags) and Divar (n = 2542).

Using the furthest neighbor assignment algorithm implemented in DOTUR, more than 2600 unique V6 tag sequences were placed in 1561 operational taxonomic units (OTUs) that were < 3% dissimilar at Tuvem. Non-parametric statistical analysis (ACE) predicted the presence of 2166 OTUs at Tuvem. At Divar, 24% of the total trimmed sequences were unique with 2198 OTUs while the predicted OTUs were estimated to be ~ 3300. Similarly, species richness estimated by the Chao1 estimator was higher at Divar than Tuvem (Table 1). Rarefaction curves did not reach a maximum and continued to increase with increasing sample size (Figure 5), indicating that the complete diversity at both sites had not been sampled. Phylogenetic classification of the twenty most abundant clusters (≥ 97% similarity) at Tuvem showed that they belonged to the class Deltaproteobacteria while at Divar, sequences from taxonomic groups within the Gammaproteobacteria were dominant. Large portion of the bacterial community were site specific (Figure 6). Only about 550 of the OTUs were common between both sites with 35% representation from Tuvem and 25% from Divar.

Figure 5

Rarefaction curves based on 95 and 97% sequence similarity of sediment bacterial community at Tuvem (T) and Divar (D).

Figure 6

Venn diagram showing overlap of operational taxonomic units (OTUs) at Tuvem (n = 1561) and Divar (n = 2198). The OTUs represent clusters with 97% sequence similarity. The numbers in the figure indicates the OTUs.

Discussion

Earlier investigations at Tuvem and Divar have focussed on examining the occurrence of denitrification, associated processes and functional diversity of denitrifiers (Fernandes ; Fernandes ). The major denitrifying communities at these locations belong to uncultured microorganisms clustering within Proteobacteria (Fernandes ). As denitrification is widespread in many genera, it would be interesting to obtain a census of bacteria in the area. High throughput pyrosequencing permits culture-independent analysis of bacteria in soil DNA extracts (Roesch ). Under the aegis of the International Census of marine microbes (ICoMM), this technique was used to obtain an initial appreciation of the bacterial diversity in relatively pristine and anthropogenically-influenced mangrove ecosystems of Goa, India.

Previous studies have shown predominant bacterial phylotypes in mangrove sediments to cluster within Proteobacteria, Bacteroidetes, Gemmatimonadetes, Actinobacteria and Firmicutes (Zhang ). Application of the pyrosequencing technology revealed that sediments at Tuvem and Divar harbour all major sediment bacterial groups affiliated with the phyla Proteobacteria, Bacteroidetes, Firmicutes, Chloroflexi, Planctomycetes, and Actinobacteria. This study also revealed the existence of bacteria belonging to other phyla like Acidobacteria, Gemmatimonadetes and members of the candidate divisions OD1 and WS3 in mangrove sediments. The phylum Gemmatimonadetes has been found earlier in systems with high nutrient input (Li ). The candidate divisions OD1 and WS3 have been reported to occur in anoxic sediments and their occurrence has been attributed to sulfur cycling (Kirkpatrick ). Thus, their occurrence in mangrove sediments is suggestive of their participation in sulfate reduction.

Consistent with previous studies of marine sediments (Bowman ; Zhao and Zeng, 2008), our results too have shown Proteobacteria to be the most abundant in mangroves. As in Divar, the Proteobacteria at Tuvem constituted > 40% of the total tags. As determined by tag abundance, the proteobacterial community at Tuvem was dominated by members of the class Deltaproteobacteria, mainly members of the order Desulfobacterales. Dos Santos have also reported higher occurrence of Deltaproteobacteria in pristine mangrove sediments. Earlier studies by LokaBharathi have reported their occurrence in mangrove swamps of the Zuari estuarine system in Goa. At Tuvem, the redox-potential of the sediments is ~ −200 mV implying that the benthic environment is largely anaerobic, more so than at Divar. Anoxic environments are known to be dominated by Deltaproteobacteria (Schwarz ). This class of bacteria have been reported to occur in coastal (Paisse ; Zhang ), continental shelf (Hunter ) as well as cold-seep sediments (Reed ). Bacteria belonging to the order Desulfobacterales have been implicated to be involved in sulfur cycling (Vrionis ) in particular sulfate reduction (Reed ). At Tuvem and Divar, sulfate-reduction has been reported to occur by Attri . It is possible that Deltaproteobacteria in this anoxic habitat not only participate in sulfur cycling but also prevent accumulation of metals (Attri ) and inorganic nitrogenous compounds by obtaining energy from the reduction of Fe(III), Mn(IV) and nitrate (Greene ).

Members of the phylum Actinobacteria are ubiquitous in estuarine and oceanic environments. Nearly 15% of the sequences recorded at Tuvem and Divar belonged to representatives of the phylum Actinobacteria making it the next most abundant phylum following Proteobacteria. Many Actinobacteria are of economic importance (Ward and Bora, 2006) as they are a source of antibiotics (Adinarayana ; Ellaiah and Zeeck, 2006; Kim ; Manivasagan ). Like Deltaproteobacteria, they play multiple roles in the environment that include degradation of cellulose (Pankratov ), hydrocarbons (Harwati ), metal oxidation (Bryan and Johnson, 2008; Johnson ), and nitrate reduction (Van Keulen ). The presence and activity of Actinobacteria in these sediments could thus be vital in altering the benthic chemistry.

The bacterial community at the anthropogenically-influenced Divar was relatively more diverse than at Tuvem. The Divar community had ~ 3300 phylotypes which is an order of magnitude higher than previously reported in marine sediments (Zhang ). Bacterial diversity can be lower in stressed environments, for example in areas prone to heavy metals (Hu ) or hydrocarbon contamination (Greer, 2010). The Divar ecosystem receives high inputs of anthropogenically-derived organic and inorganic material. Consequently, availability of a wide variety of substrates could result in high taxonomic and metabolic diversity. Deltaproteobacteria were the predominant OTU clusters at Tuvem whereas at Divar, OTUs belonging to the Gammaproteobacteria were dominant. Molecular investigations in a Chinese mangrove ecosystem (Liang ) showed that the Gammaproteobacteria-affiliated sequences constituted the largest portion of their clone library. The Gammaproteobacteria are active mediators of the N, S and C cycles. Some of the Gammaproteobacteria recorded at Divar are classified under the order Chromatiales, Alteromonadales and Oceanospirillales. Species belonging to the genus Marinobacter, Shewanella, Alteromonas of order Alteromonadales have been commonly found in the marine environment (Zhuang ). These genera have also been recorded in the present investigation. Marinobacter spp. use a variety of hydrocarbons as the sole source of carbon and energy (Gauthier ). These halophilic bacteria also have a high capacity for denitrification (Yoshie ). Similarly, Shewanella (Brettar ; Zhao ) and Alteromonas(Haijun, 2002) are known to mediate the reductive phase of the N cycle. In the present study, denitrification, a key process in the nitrogen cycle, was examined. At the time of sampling, near fresh water conditions prevailed in the mangrove systems as a result of rainfall over the region which was evident from low salinity values of ≤ 1 psu. The Divar ecosystem is prone to high external nutrient input due to land runoff during the monsoon (Divya ). As expected, pore water analysis in the present study revealed higher concentration of inorganic nitrogenous compounds and organic carbon content at Divar. Denitrification activity was also found to be significantly higher at Divar as compared to Tuvem providing evidence for the active occurrence of alternate respiratory pathways which consequently help to reduce the nutrient load in the aquatic system. Further, degradative pathways like denitrification also result in the oxidation of organic matter. The Divar sediments are organic carbon rich and the existence of metabolically versatile bacteria could indicate their contribution to the C and N cycle. Fast growing Gammaproteobacteria dominate estuaries (Greer, 2010) and they have a preference for elevated concentrations of nutrients (Pinhassi and Berman, 2003). The Gammaproteobacteria form the most abundant denitrifying communities in marine sediments (Bhatt ). About 96% of cultured denitrifiers belong to the Gammaproteobacteria (Brettar ). The culturable and non-culturable (probed based on the phylogenetic diversity of nosZ genes) denitrifiers from Tuvem and Divar mostly belong to the class Gammaproteobacteria (Fernandes, 2010; Fernandes ). Nitrification has also been reported to occur at the same sampling locations and has been found to be mediated by bacteria belonging to the Beta and Gammaproteobacteria (Krishnan, 2010). The pyrosequencing approach has revealed a low abundance of nitrifying Betaproteobacteria-related phylotypes. An abundance of sulfide-oxidizing denitrifiers (class Gammaproteobacteria, order Chromatiales) has also been observed (Mori ). It is possible that the oxidative pathway of nitrogen in mangrove sediments is also more likely to be mediated by members belonging to the class Gammaproteobacteria.

Rarefaction curves at both locations did not reach an asymptotic stage indicating that the bacterial richness in the present study was not fully covered. Statistical richness estimates of Chao1 and ACE indicated that species richness in mangrove sediments was less than tidal mud flats of Dongmak, Korea (Kim ). Contrastingly, when compared to the water columns of ocean,(Kirchman ) the species richness in mangrove sediments is far greater. A considerable fraction of the low-abundance OTUs of the so called “rare biosphere” (Sogin ) were responsible for the high diversity observed in mangrove sediments and indicate that they have the potential to become dominant when favorable environmental conditions arise. Like in most molecular based taxonomic surveys of bacterial communities from environmental samples, a significant fraction of the bacterial population at both the locations are new as they have not been identified up to the class level. About 18% of the total tag sequences recovered at Divar and 22% at Tuvem belonging to various microbes await to be cultured and identified. These sequences have not been reported in earlier studies carried out in mangroves. It is evident that previous reports may have underestimated bacterial diversity since they incorporated tedious cloning procedures and a limited number of sequences. The Venn diagram illustrated that majority of the OTUs in the present study were specific to the sampling sites with < 35% of the sequence clusters being common between both the sites. Fewer samples in this study did not permit the use of statistical methods to assess the inter-relationships between environmental characteristics and bacterial diversity. However, the largely unique bacterial community, higher inorganic-N content, DNT and pre-dominance of clusters belonging to Gammaproteobacteria at Divar indicate that the population structure is influenced by the prevailing environmental conditions.

Though our analysis was limited to single sediment samples from two different locations, the pyrosequencing technique was useful in revealing sustenance of higher bacterial diversity and richness at the anthropogenically-influenced Divar as compared to Tuvem. This finding suggests that the levels of pollutants are not a deterrent for affecting bacterial diversity at Divar. On the contrary, they could play an important role in buffering the anthropogenic influences. Higher occurrence of tag sequences associated to class Gammaproteobacteria at Divar was also observed which is suggestive of their role in lowering nitrate levels through DNT. Thus, the prevailing environmental conditions could be crucial in influencing the composition of the autochthonous benthic bacterial communities which play a vital role in ecosystem functioning. More interestingly, a significant fraction of the novel bacterial population at both the locations await to be cultured and identified for various applications. Future studies would be directed to understand the dynamics of the microbial populations at a higher spatial and temporal resolution.