Insights on aquatic microbiome of the Indian Sundarbans mangrove areas

Abstract

Background

Anthropogenic perturbations have strong impact on water quality and ecological health of mangrove areas of Indian Sundarbans. Diversity in microbial community composition is important causes for maintaining the health of the mangrove ecosystem. However, microbial communities of estuarine water in Indian Sundarbans mangrove areas and environmental determinants that contribute to those communities were seldom studied.

Methods

Nevertheless, this study attempted first to report bacterial and archaeal communities simultaneously in the water from Matla River and Thakuran River of Maipith coastal areas more accurately using 16S rRNA gene-based amplicon approaches. Attempt also been made to assess the capability of the environmental parameters for explaining the variation in microbial community composition.

Results

Our investigation indicates the dominancy of halophilic marine bacteria from families Flavobacteriaceae and OM1 clade in the water with lower nutrient load collected from costal regions of a small Island of Sundarban Mangroves (ISM). At higher eutrophic conditions, changes in bacterial communities in Open Marine Water (OMW) were detected, where some of the marine hydrocarbons degrading bacteria under families Oceanospirillaceae and Spongiibacteraceae were dominated. While most abundant bacterial family Rhodobacteracea almost equally (18% of the total community) dominated in both sites. Minor variation in the composition of archaeal community was also observed between OMW and ISM. Redundancy analysis indicates a combination of total nitrogen and dissolved inorganic nutrients for OMW and for ISM, salinity and total nitrogen was responsible for explaining the changes in their respective microbial community composition.

Conclusions

Our study contributes the first conclusive overview on how do multiple environmental/anthropogenic stressors (salinity, pollution, eutrophication, land-use) affect the Sundarban estuary water and consequently the microbial communities in concert. However, systematic approaches with more samples for evaluating the effect of environmental pollutions on mangrove microbial communities are recommended.

About the Authors:

Paltu Kumar Dhal

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Supervision, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Life Science and Biotechnology, Jadavpur University, Kolkata, India

ORCID logo http://orcid.org/0000-0001-5790-3268

Germán A. Kopprio

Roles Writing – review & editing

Affiliation: Department of Biogeochemistry and Geology, Tropical Marine Microbiology, Leibniz Center for Tropical Marine Research, Bremen, Germany

Astrid Gärdes

Roles Funding acquisition, Project administration, Supervision, Writing – review & editing

Affiliation: Department of Biogeochemistry and Geology, Tropical Marine Microbiology, Leibniz Center for Tropical Marine Research, Bremen, Germany

Abstract

Background

Anthropogenic perturbations have strong impact on water quality and ecological health of mangrove areas of Indian Sundarbans. Diversity in microbial community composition is important causes for maintaining the health of the mangrove ecosystem. However, microbial communities of estuarine water in Indian Sundarbans mangrove areas and environmental determinants that contribute to those communities were seldom studied.

Methods

Nevertheless, this study attempted first to report bacterial and archaeal communities simultaneously in the water from Matla River and Thakuran River of Maipith coastal areas more accurately using 16S rRNA gene-based amplicon approaches. Attempt also been made to assess the capability of the environmental parameters for explaining the variation in microbial community composition.

Results

Our investigation indicates the dominancy of halophilic marine bacteria from families Flavobacteriaceae and OM1 clade in the water with lower nutrient load collected from costal regions of a small Island of Sundarban Mangroves (ISM). At higher eutrophic conditions, changes in bacterial communities in Open Marine Water (OMW) were detected, where some of the marine hydrocarbons degrading bacteria under families Oceanospirillaceae and Spongiibacteraceae were dominated. While most abundant bacterial family Rhodobacteracea almost equally (18% of the total community) dominated in both sites. Minor variation in the composition of archaeal community was also observed between OMW and ISM. Redundancy analysis indicates a combination of total nitrogen and dissolved inorganic nutrients for OMW and for ISM, salinity and total nitrogen was responsible for explaining the changes in their respective microbial community composition.

Conclusions

Our study contributes the first conclusive overview on how do multiple environmental/anthropogenic stressors (salinity, pollution, eutrophication, land-use) affect the Sundarban estuary water and consequently the microbial communities in concert. However, systematic approaches with more samples for evaluating the effect of environmental pollutions on mangrove microbial communities are recommended.

Figures

Fig 6

Table 2

Fig 1

Fig 2

Table 1

Fig 3

Fig 4

Fig 5

Fig 6

Table 2

Fig 1

Fig 2

Table 1

Citation: Dhal PK, Kopprio GA, Gärdes A (2020) Insights on aquatic microbiome of the Indian Sundarbans mangrove areas. PLoS ONE 15(2): e0221543. https://doi.org/10.1371/journal.pone.0221543

Editor: Jiang-Shiou Hwang, National Taiwan Ocean University, TAIWAN

Received: August 7, 2019; Accepted: February 4, 2020; Published: February 25, 2020

Copyright: © 2020 Dhal et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files

Funding: This work was part of the Joint NAM S&T Centre – ZMT Bremen Fellowship in Tropical Coastal Marine Research that was supported by the Leibniz Zentrumfür Marine Tropenforschung (ZMT), Bremen, Germany

Competing interests: The authors have declared that no competing interests exist.

1. Introduction

Sundarbans, the largest mangrove forest of the world, is situated in the joint delta of Ganges, Brahmaputra and Meghna rivers at Bay of Bengal [1,2]. This UNESCO World Heritage site comprises the Indian state of West Bengal and southwest Bangladesh [3]. Livelihood and well being of millions of people live in and around of Sundarbans, depends on its status and ecological services. Despite its high ecological and economical values, Sundarbans is seriously threatened by different anthropogenic activities. Since the early 19th century, landscapes of Sundarbans have also been changing due to saline and freshwater imbalances. Water quality of this ecosystem is largely affected by sewage pollutant originated from industries located upstream and urban areas of West Bengal. Sewage entering into coastal water contains diverse chemical and microbiological pollutants and a wide variety of organic and inorganic wastes [4,5], driving changes on its ecological and physiological health.

Microbial communities of mangroves are responsible for nutrient cycling and play a vital role in productivity, conservation and rehabilitation of mangrove ecosystems [6]. Therefore, understanding their responses to environmental changes is essential to predict changes in service-provisioning [7]. Several recent studies described the microbial community compositions of surface sediments and water of Indian Sundarban mangrove areas. Surface sediments this area dominanted with Deltaproteobacteria followed by Gammaproteobacteria, Alphaproteobacteria, Betaproteobacteria, and Epsilonproteobacteria under phylum Proteobacteria. Abundant bacterial orders are Desulfobacterales, Desulfuromonadales, Myxococcales, and Bdellovibrionales. [8–10]. While bacterioplankton communities in the water of this region were found to be abundant with Gammaproteobacteria and Alphaproteobacteria. At the family level dominancy of Hyphomicrobiaceae, Rhodobacteraceae, Pseudomonadaceae, Erythrobacteraceae, Kordiimonadaceae, Hyphomonadaceae, and Ruminococcaceae were observed [11–13]. However, sampling locations of those studies on microbial communities in the Indian Sundarban mangrove water mainly restricted near to an island (Sagar Island) and other estuary of Mooriganga, Thakuran, Matla and Harinbhanga, therefore the major conclusions of these studies were made based on a limited number of samples. Moreover, not much effort has been made to investigate the archaeal community of this region except single report by [3] on surface and subsurface sediments of Indian Sundarban mangrove forest. Moreover, the above studies have rarely analyzed the bacterial and archaeal community structures of the same samples at the same time. Therefore, our knowledge on those communities as well as information on how they are controlled by environmental parameters is limited. In order to assess the microbial communities of marine ecosystem via high-throughput sequencing of amplified 16S rRNA genes with high resolution and fidelity, it is extremely important to select the proper primer set that can’t underestimated or overestimated any common marine taxa [14]. However, this will be the first attempt to visualize the accurate and well-resolved picture of bacterial and archaeal communities simultaneously of marine water in Sundarbans mangrove using next-generation amplicon sequencing of the 16S rRNA gene using recently developed 515F-Y/926R primers that target V4-V5 region of 16S rRNA gene. We also tried to explore the environmental determinants that contribute to the variation of their microbial communities. This study will provide baseline knowledge on microbial ecology of the World Heritage site and serve as a baseline for monitoring programs and predicting changes at impacted sites.

2. Material and methods

2.1 Ethics statement

The water samples were brought to ZMT, Germany for its environmental parameters and nutrients measurements. The extracted mutagenic DNA were also carried to ZMT, Germany for sequencing and further analysis. In this reference, The National Biodiversity Authority (NBA), Government of India is well aware of this activity (NBA/Tech Appl/9/Form B/52/17/17-18/2985) that the above mentioned samples are being used only for research purpose. Further no permits were prerequisite for the defined field studies, which complied with all relevant regulation neither the studied locations are not privately owned.

2.2 Study sites and sample collection

In the present study sampling was conducted in the Sundarbans mangrove ecosystem that shared between India and Bangladesh and lies in the Ganga-Brahmaputra-Meghna (GBM) delta. This mangrove ecosystem contains over 102 islands with a network complex of many rivers, rivulets, and creeks [15]. Sampling was carried out at two different locations on Thakuran River—Matla River estuarine complex of Maipith coastal areas in the Indian Sundarbans mangroves during March 2017 during low tide (Fig 1). They are designated as Island of Sundarban Mangroves (ISM) and Open Marine Water (OMW). ISM is an uninhabited small island with lesser anthropogenic disturbance situated in Thakuran and Matla river complex in low-lying coastal plain. This river has no reports for perennial fresh water source [15]. Water from three different coastal regions of this island is selected for sampling. OMW is an open marine site around same regions, which is supposedly continually influenced by the wastewaters from upstream regions of Matla River. Three independent replicated water samples (1 L) from each of three different sampling sites of both ISM and OMW were collected in sterile containers and immediately stored at a chilled box until further laboratory analysis.

[Figure omitted. See PDF.]

Fig 1. Map of the sampling area: Water samples were collected from three stations (KL 1, KP 1 and TH1) of a small Island are named as ISM and open marine water samples named as OWM (TH 2, KL 2 and BL 1).

Three biological replicates from each of the six stations; therefore, total eighteen (18) samples (nine from ISM and another nine from OMW) were collected for this investigation. The GPS data of sampling sites were collecting during sampling and compiled together in an ArcGIS10.3 software environments and finally map has been prepared using open source database of GADM (https://gadm.org/).

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

2.3 Environmental parameters and nutrients analysis

Physiochemical parameters (salinity and pH) of all collected samples were measured using Eureka 2 Manta multiprobe (Eureka Environmental Engineering, Texas, USA). Total 50 mL of each sample was filtered through a 0.7 μm syringe filter and poisoned with 200 μL of 3.5 g/100 mL HgCl2 solution for nutrient analysis. Now each of the treated samples was filtered through a 0.7 μm pore size GF/F filter (GE Healthcare Bio-Sciences, Pittsburgh, PA, USA) for DOC and nutrients measurements. The total nitrogen (TN) content of each the samples were measured using the EuroVector EA 3000 elemental analyzer. For measurements of DOC, the filtered samples were acidified with concentrated HCl (pH <2) and analyzed by high-temperature catalytic oxidation method method using a TOC-VCPN analyzer (Shimadzu, Mandel, Canada). Seawater standards (Hansell laboratory, RSMAS University Miami, USA) were used for calibration and quality control, and ultrapure water as a blank. The Dissolved inorganic nutrients that includes combined nitrate and nitrite (NOx), phosphate (PO43-P), and silicate (Si(OH)4-) were analyzed using a continuous flow analyzer (Flowsys by Unity Scientific, Brookfield, USA) and detected spectrophotometrically as a colored complex [16] (https://doi.pangaea.de/10.1594/PANGAEA.889699).

2.4 DNA extraction, PCR amplification, and Illumina MiSeq sequencing

From each site, water (1 L) was filtered (0.2 μm) and DNA was extracted using the PowerWater® DNA Isolation Kit according to the manufacturer’s instructions (MoBio Laboratories Inc., Carlsbad, CA, USA). DNA concentrations and purity were measured spectrophotometrically. Presence of bacterial and archebacterial 16S rRNA gene was in the extracted metagenome was verified following previous method [17]. In order to classify taxonomically both bacterial and archaeal community structure simultaneously, sequencing of V4–V5 hypervariable regions of 16S rRNA gene were generated using primers 515F-Y (5′-GTGYCAGCMGCCGCGGTAA-3′) and 926R (5′-CCGYCAATTYMTTTRAGTTT-3′) [14] on the Illumina MiSeq platform (CeBiTec Bielefeld, Germany), in a 2 × 300 bp paired-end run.

2.5 High throughput sequencing data processing

Primer sequences were removed using cutadapt from the raw paired-end reads [18]. The primer-trimmed sequences are available on Sequence Read Archive (SRA) (accession no. SRP144285). Sequences were quality trimmed with trimmomatic v0.32 [19] using a sliding window of 4 bases and a minimum average quality of 15, and merged with PEAR v0.9.5 [20]. Quality-filtered sequences were clustered into OTUs with swarm algorithm using default parameters [21]. One single representative sequence per OTU was taxonomically classified with SINA (SILVA Incremental Aligner; v1.2.11; Silva reference database release 132) at a minimum alignment similarity of 0.9, and a last common ancestor consensus of 0.7 [22]. OTUs that were unclassified on the domain level and those matching to chloroplast and mitochondrial sequences were excluded from the analysis using well-standardized r script [23]. The final OTU tables are accessible at (https://doi.pangaea.de/10.1594/PANGAEA.890757).

2.6 Statistical analysis

Principal component analysis (PCA) was performed to cluster the sampling sites based on their environmental parameters. Differences in environmental parameters among ISM and OMW were assessed using general linear mixed models (GLMM) with sampling station as a random factor [24].

Alpha-Diversity indices were calculated to assess richness and evenness of the microbial communities [25] in the studied samples, based on repeated random subsampling of the amplicon data sets after randomly rarefying the data set to the minimum library size (50517 sequences). Significant differences in alpha-diversity indices between the studied stations were determined by using the non-parametric Kruskal test followed by p-value adjusted Wilcoxon tests [26].

To assess the differences in community structure between two sampling sites (beta-diversity), Bray–Curtis dissimilarities were calculated using the relative OTU abundances and also non-metric multidimensional scaling (NMDS) plot was produced. Analysis of similarity (ANOSIM) was calculated to assess the separation of bacterial communities between the two sites. P-values of all multiple pairwise comparisons were adjusted using the false discovery rate (fdr) correction method by [27]. In order to evaluate the environmental parameters as drivers of the variations in community compositions, redundancy analysis (RDA) was used with centered log ratio (clr)-transformed sequence counts using the R function aldex.clr of the ALDEx2 package via median values of 128 Monte-Carlo instances [28]. To compare the explanatory power of all measured environmental parameters, additional RDA models were constructed with environmental parameters as predictors. Forward model selection was used after checking for variance inflation to determine which of parameters would be included in the RDA models. When more than one parameter was included, pure effects were also tested accounting for the variation explained by the other factors in the model. Collinearity among predictors was determined via Variance inflation factors (VIFs) of the individual parameters. All of the parameters in any of the RDA models displayed VIFs less than 10. The adjusted R2 is provided as goodness-fit-stat. All statistical analyses were conducted in R using the core distribution, version 3.3.2 and R-Studio, version 1.0.153,with following packages: vegan [29], lmerTest for the GLMM [24], ALDEx2 [29] and multcomp [28].

3. Result

3.1 Environmental characterization

Environmental parameters (pH, salinity) and nutrients (including DOC, TN, NOx, nitrate, DIN, phosphate: PO43-P, and silicate: Si(OH)4-) concentrations for all samples were measured (Table 1). Samples were slightly alkaline (pH 8.0 to 8.7) in nature. The GLMMs analysis indicated that the measured water nutrients that differed significantly among the two sampling station were mainly TN, DIN and PO43-P (Table 1; S1 Table). The PCA ordination (Fig 2) showed that first two principal components (PC1 and PC2) represented 74.4% of data variation among sites. PC1 alone represents 60.1% of total variation and influenced by most of the measured parameters, while pH showed a strong correlation with PC2 (S1 Fig). Noteworthy to mention, the samples were separated into two clusters by PC1. One cluster is mainly composed with the samples from ISM (except TH2.3 of OMW) and other cluster accommodating samples collected from OMW. This ordination probably indicates elevated eutrophication in samples from the OMW compared to ISM.

[Figure omitted. See PDF.]

Fig 2. Principal component analysis (PCA) to ordinate the eighteen collected water samples collected samples from ISM and OMW based on their environmental parameters.

The arrows show the direction of the environmental parameters. DIN, dissolved inorganic nitrogen; TN, total nitrogen; DOC, dissolved organic carbon.

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

[Figure omitted. See PDF.]

Table 1. Environmental conditions in Indian Sundarban mangroves and details of the sampling sites.

https://doi.org/10.1371/journal.pone.0221543.t001

3.2. Microbial communities

Total numbers of reads generated per sample ranged between 50517 to 90468 (after merged) corresponding to 3,644 to 6,470 swarmed, non-singleton 16S OTUs (S2 Fig). After rarefaction, numbers of bacterial and archaeal OTUs ranged between 3390 to 5415 and 37 to 91, respectively (S2 Fig). None of the measured diversity indices (Average Shannon diversity and inverse Simpson indices) were found significance differences in between OMW and ISM (Fig 3; S2 Table) indicated by Kruskal—Wilcoxon test (p > 0.5), although values varied considerably.

[Figure omitted. See PDF.]

Fig 3. Alpha diversity of the water microbial community at two different sites (ISM and OMW) of Sundarban mangrove forest areas.

Values are calculated based on repeated random subsampling to the lowest number of sequences per sample. The median per group presented in black line.

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

The microbial community of marine estuary water from Sundarbans was dominated with bacteria occupying more than 96% of total community and archaea represented only 4%. Bacterial assemblage of two different sites this area showed to have a distinct community. In class level, among the dominant bacterial groups, Flavobacteria (ISM: 15.9% vs OMW: 8.6%), Alphaproteobacteria (ISM: 29.5% vs OMW: 28%), and Acidimicrobiia (ISM: 6.6% vs OMW: 5.0%) were dominant in ISM while OMW was dominated with mainly with Gammaproteobacteria (ISM: 22.6% vs OMW: 35.3%) (S3 Fig). At higher taxonomic resolution levels (Fig 4), bacterial communities were composed with a total of 474 and 915 different bacterial family and genus, respectively. The most dominant bacterial family was Rhodobacteraceae (18.6%), almost equally distributed between studied two sites. Other dominant bacterial families of ISM were Flavobacteriaceae (14.8%) and OM1 clade (5.2%) whereas in OMW, Oceanospirillaceae (16%) and Spongiibacteraceae (4%) were the most abundant.

[Figure omitted. See PDF.]

Fig 4. Taxonomic composition of dominant bacterial taxa on family level across eighteen samples under sites ISM and OWM (nine samples each).

Ten (10) most abundant bacterial families for each of the samples were reported here and rests less dominant members are label as “other”.

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

We observed dominancy of Marine Group (MG I) (currently known as Thaumarchaeota) and Euryarchaea MG II in archaeal community assemblages with 78.5% and 16.9% of relative abundance, respectively while the presence of Woesearchaeota (2.1%) was also evident (Fig 5). MG I was found in relatively higher abundant at ISM constituting on average 82.1% of sequences as opposed to 75% at OMW. Euryarchaeota MG II comprised about 20.3% at OMW compared to 13.7% at ISM. Among total twenty-nine (29) archaeal genera, Candidatus Nitrosopumilus and Candidatus Nitrosopelagicus accounted for the 40.7% and 21.4% of total relative abundance, respectively (S4 Fig).

[Figure omitted. See PDF.]

Fig 5. Taxonomic compositions of dominate archaeal phyla across eighteen samples represents two sites ISM and OWM (nine samples each).

Ten (10) most abundant phyla for each of the samples were reported here and rests less dominant members are label as “other”.

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

3.3 Environmental drivers of bacterial communities

At OTU resolution level also, distinct microbial communities were observed between OMW and ISM based on changes in community structure (beta diversity) which is quantified by non-metric multidimensional scaling (NMDS) plot by calculating Bray–Curtis dissimilarly (Fig 6). This pattern is confirmed by the ANOSIM test that indicated a significant difference in microbial community structure between ISM and OMW (ANOSIM, R = 0.24, p < 0.001). Redundancy analyses attempted to identify the water quality parameters that had strong explanatory power for microbial communities. We observed that total nitrogen (TN) and dissolved inorganic nutrients (DIN) accounted for almost 10% of the variability in microbial community of OMW where TN alone explain 6% variation of microbial community (RDA, R2 = 0.06, F(1,7) = 1.34, p < 0.05). In contrast, salinity and TN explained approximately 9% of the variability in community composition of ISM (RDA, R2 = 0.06, F(1,7) = 1.34, p < 0.05) and alone salinity responsible for explaining 7% microbial variation of this site (Table 2).

[Figure omitted. See PDF.]

Fig 6. Non-metric multidimensional scaling (NMDS) plot of bacterial community composition of the bacterial communities of each sampled at the inhabited island (ISM) and open marine areas (OWM).

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

[Figure omitted. See PDF.]

Table 2. Contribution of environmental parameters including nutrient content of six sampling sites to explaining the variation in microbial community composition based on redundancy analysis (RDA).

https://doi.org/10.1371/journal.pone.0221543.t002

4. Discussion

The pH values (8.0–8.7) indicates the water of ISM and OMW slightly alkaline in nature which supports the previous findings in similar samples from Sundarbans Mangrove forest areas [30]. Such ranges of pH may be attributed by the buffering capacity of water that support high biological activity [31,30]. The water of Sundarbans is characterized by elevated salinity values in line with previous reports [31]. The long-term changes in water properties in the eastern part of Sundarbans, sampling regions of our study, indicating increased trends on salinity and pH [3]. Our result shows differences in measured environmental parameters between two sites represented by three sampling stations and nine samples each leading to their segregation into two clusters (in PCA analysis) along with their sampling sites. This ordination as a result of different nutrient loads and this is reflected by potential eutrophication in water from OWM. Influences from the Thakuran and Matla rivers reported to have a strong impact on the estuary water quality represented by OWM [4,5], that also reflected in our study. The perennial discharge in Thakuran and Matla Estuary from upstream regions brings in a high suspended matter load throughout the year [30]. Those estuaries severely contaminated with huge organic load and sediment flux originated from upstream domestic sewage, aquaculture, intensive trawling activities, agricultural runoff as well as soil erosion [32,33] may supports the prospective cause for eutrophication in the OWM sites.

Because of the relevance of microbial community of Indian Sundarbans, several investigators attempted on surface sediments samples [34,8,3,10,35] as well as recently on water column of this regions [12,13] using 16S rRNA gene metagenomic approaches. Unlike previous studies, this investigation attempted to asses both the bacterial and archaeal community at a same time of water from relative less anthropogenic disturbance sites using an efficient primer set to target V4-V5 variable region of 16S rRNA gene in order to avoid the problems of underestimated or overestimated common marine taxa [14], therefore our investigation gives more accurate and well-resolved picture of microbial communities of these sites.

Although insignificant differences, elevated trends of α-diversity of the marine estuary water samples (OWM) might be an indications of relatively rich bacterial community compared to ISM of Sundarbans might be attributed toward their elevated eutrophication. This observation was supported by previous reports that indicate a higher diversity and equitability in the human impacted estuary because of proliferation of several different microorganisms [36,37].

The bacterial assemblages of the studied samples (specially OMW) showed similarity with the previously reported bacterial community of marine sediments and water samples of Sundarban Mangrove areas [8, 10–12]. The most dominant bacterial family presents both the station with almost equal proportion is Rhodobacteraceae. Dominancy of members of this family in marine water microbial community previously reported and known to be associated with marine phytoplankton blooms where it plays important role in transforming phytoplankton-derived organic matter [34,38,39]. The abundant OTUs of Rhodobacteraceae are classified as anoxygenic phototrophs Nautella, reported to serve an indicator of marine eutrophication, predominantly found in higher eutrophic OMW samples. The other dominant one under the same family is marine heterotrophs Ruegerias (almost equally distributed among both ISM and OMW) serve as the model of marine sulfur and carbon cycle [38,40,41].

Interestingly enrichment of several OTUs from the Flavobacteria is observed in the oligotrophic ISM compared to the eutrophic OMW. They are specialized in utilization of biopolymers and organic substances in oligotrophic environment i.e., when organic substances present even at very low concentrations. Higher abundance of these polymers degrading bacteria responsible for releasing the nutrients heterotrophs/copiotrophic organisms thus helps in the growth of heterotrophic bacterial community at oligotrophic marine environments [42–44]. Therefore, they might plays central role for balancing the heterotrophs/copiotrophic microbial ecology of ISM. The dominant OTUs under family Flavobacteriaceae were mainly classified as Aureimarina and NS5 marine group genus. Roles of Aureimarina in marine biogeochemistry has not been investigated much although few studies reported on their presences in coastal seawater and saline estuarine [45,46]. However, this is the first report of their abundance in marine water of Sundarbans. The NS5 marine group which are equally dominated in both the studied sites are reported to be ubiquitous in the seawater-related samples and known for phytoplankton-derived macromolecules [47–49].

This investigation identified Actinobacteria constituted a predominant fraction both in OMW and ISM but elevated amount in the later samples. Bacteria under this group are consisted of both copiotrophic and oligotrophic members with higher abundance in oligotrophic marine environments [50,51]. Their presence in mangrove estuary regions is well documented [52–54]. As marine Actinobacteria are the richest sources of secondary metabolites thus, have been well reported as potential sources of bioactive compounds [55]. Therefore, their abundance in our studied sites (specially ISM) would be potential hotspot for isolating bioactive molecules from Indian Sundarban mangrove forest. The OM1 clade (dominating in ISM), an uncultured Actinobacterial clade, frequently recovered from various marine environments with higher abundance at near coastal sites than open marine areas however supports our reports [56–58]. The dominant OTUs of this family were classified as Candidatus Actinomarina. Those photoheterotrophs are one of smallest free-living prokaryotes are reported to be ubiquitous in marine systems. Not many reports are described their role in the marine biogeochemical cycle but Actinobacteria are reported to be associated in carbon cycling to decompose the plant biomass via degrading the cellulose and hemicellulose materials, a dominate resources materials in mangrove, however may supports their dominancy in ISM. They are also known as chitin, hydrocarbons and organic contaminants degrader [59,60].

The OTUs affiliated to families Oceanospirillaceae and Spongiibacteraceae of Gammaproteobacteria showed increase abundance in the impacted site OMW. Bacteria from these families are known to be present in eutrophic marine environments. They are known as polymer degraders and can utilize polyhydroxy alkanoate compounds and proteorhodopsin, for harvesting an additional energy, supports their living in eutrophic water samples [61–63]. The dominant OTUs of Oceanospirillaceae are affiliated to chemoheterotrophic genus Marinobacterium. Their presence in mangroves as well as surface seawater have already been described in previous studies and known to be associated with hydrocarbon biodegradation [64,65]. The other dominated bacterial family in the samples from OMW is Spongiibacteraceae. They comprise mainly marine bacteria known as Oligotrophic Marine Gammaproteobacteria (OMG) group [66,67]. We recorded the dominant OTUs of this family are affiliated with BD1-7 which is a cosmopolitan group of Gammaproteobacteria is mostly autochthonous, reported to inhabits at diverse marine habitats [68,65,69]. In line with previous reports this investigation, therefore, indicates proliferation of bacterial groups under Gammaproteobacteria with respond to increased nutrient concentrations in estuary [36].

However, in contrast to sediments reported in previous investigations, an archaeal community of marine waters in the Sundarban mangroves is dominated with Thaumarchaeota Marine Group (MG I) and Euryarchaea MG II. The chemolithoautotrophic MG I which are in higher in number on ISM are responsible for oxidation of ammonia and showed ability in inorganic carbon fixation [70] thus important players in global Carbon (C) and Nitrogen (N) biogeochemical cycles. While the heterotrophic MG II, dominated in OMW, previously reported to be abundant in the marine aquatic environment [71]. Their abilities to acquired energy in presences of light through organic carbon degradation in the photic zone is also documented [71]. Dominant OTUs of MGI group are affiliated with ammonia-oxidizing archaeal, Candidatus Nitrosopumilus and Candidatus Nitrosopelagicus, play important roles in nitrogen and carbon cycling of marine ecosystem [3]. However, this investigation reports first on their presence in this areas. Therefore, the biological and geochemical processes at estuary water habitats in the Indian Sundarban Mangrove areas have likely influenced the archaeal community structure.

Overall, our study indicates along with the elevated level of average pH and salinity, the open marine water (OMW) showed eutrophication probably leads to an observed bacterial shift toward more copiotrophic and photoheterotrophic bacterial (Oceanospirillaceae and Spongiibacteraceae) and archaeal community (Euryarchaea MG II) and compared to the more oligotrophic microbial community (Aureimarina, NS5 marine group, OM1 clade and Thaumarchaeota MG I) of costal water of a small Island of Sundarban Mangroves (ISM). These microbial assembles thus might represent key players in biogeochemical cycle of this mangrove and the studied areas represent a hotspot for bacterial having potential to produce the commercially important secondary metabolites. This investigation also reports that total nitrogen and dissolved inorganic nitrogen are the major environmental contributors on determining the microbial communities for OMW and for ISM it is combination of total nitrogen and salinity.

5. Conclusion

This investigation provides the first details description of bacterial and archeal communities concurrently of Thakuran and Matla river complex of Maipith coastal areas in the Indian Sundarbans mangroves areas. This study indicates the eutrophication in open marine water (OMW) dominated with more copiotrophic and photoheterotrophic bacterial and archaeal community (Euryarchaea MG II) and while oligotrophic microbial community abundant in costal water of the Island of Sundarban Mangroves (ISM).

However, given the rising burden on Indian mangrove coastal ecosystems, this study suggests that sewages from urban areas lacking proper treatment can alter microbial communities that may play vital role in biogeochemical cycle (nitrogen cycle) of mangrove ecosystem and consequently may impact on the climate in the tropical country.

6. Acknowledgments

Foremost PKD would like to thank NAM S&T Centre New Delhi, India and ZMT, Bremen for selection. PKD also acknowledge to Jadavpur University for granting the leave to avail this fellowship. We sincerely thank the faculty and technical staff of ZMT Bremen for their generous help to carry this research work. Many thanks to Dr Halina Tegetmeye of CeBiTec Bielefeld, Germany for sequences. The assistance of conducting analyses at the laboratories at the ZMT in Bremen by Matthias Birkicht, Sonja Peters and Achim Meyer also acknowledge. Dr. Debajit Datta, Assistant professor, Jadavpur University is also acknowledged for helping in designing the map of the sampling location.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Heatmap of pairwise correlations between the different environmental parameters and three principal components.

Levels of correlations are indicated with different color bar.

https://doi.org/10.1371/journal.pone.0221543.s001

(TIFF)

S2 Fig. The rarefaction curve of the eighteen (18) samples, indicated by the number of OTUs as a function of the number of reads.

The curve approaching plateau indicates that the number of reads are enough to describe the OTUs representing the community.

https://doi.org/10.1371/journal.pone.0221543.s002

(TIFF)

S3 Fig. Taxonomic composition of the ten (10) most abundant bacterial phyla in the studied two sites ISM and OWM represented by eighteen samples (nine samples each).

https://doi.org/10.1371/journal.pone.0221543.s003

(TIF)

S4 Fig. Taxonomic composition of the ten (10) most abundant archaeal genus in the studied two sites ISM and OWM represented by eighteen samples (nine samples each).

https://doi.org/10.1371/journal.pone.0221543.s004

(TIF)

S1 Table. Kruskal-Wallis test for Environmental parameters at the six sampling sites.

https://doi.org/10.1371/journal.pone.0221543.s005

(DOCX)

S2 Table. Kruskal-Wallis test for alpha diversity of the six sampling sites.

https://doi.org/10.1371/journal.pone.0221543.s006

(DOCX)

Citation: Dhal PK, Kopprio GA, Gärdes A (2020) Insights on aquatic microbiome of the Indian Sundarbans mangrove areas. PLoS ONE 15(2): e0221543. https://doi.org/10.1371/journal.pone.0221543

References

1. Spalding MD, Blasco F, Field CD. 1997. World Mangrove Atlas. Okinawa (Japan): International Society for Mangrove Ecosystems.

2. Gopal B., Chauhan M. Biodiversity and its conservation in the Sundarban Mangrove Ecosystem. Aquat. Sci. 68, 338–354 (2006)

3. Bhattacharyya A, Majumder NS, Basak P et al. Diversity and Distribution of Archaea in the Mangrove Sediment of Sundarbans. Archaea 2015;2015: 968582. pmid:26346219

4. Mitra A, Gangopadhyay A, Dube A et al. Observed changes in water mass properties in the Indian Sundarbans (northwestern Bay of Bengal) during 1980–2007. Current Science 2009: 1445–52.

5. Banerjee K, Gatti RC, Mitra A. Climate change-induced salinity variation impacts on a stenoecious mangrove species in the Indian Sundarbans. Ambio 2017;46: 492–9. pmid:27804094

6. Holguin G, Gonzalez-Zamorano P, de-Bashan LE et al. Mangrove health in an arid environment encroached by urban development—a case study. Sci Total Environ 2006;363: 260–74. pmid:16112715

7. Mishra RR, Swain MR, Dangar TK et al. Diversity and seasonal fluctuation of predominant microbial communities in Bhitarkanika, a tropical mangrove ecosystem in India. Rev Biol Trop 2012;60: 909–24. pmid:23894955

8. Basak P, Majumder NS, Nag Set al. Spatiotemporal analysis of bacterial diversity in sediments of Sundarbans using parallel 16S rRNA gene tag sequencing. Microb Ecol 2015;69: 500–11. pmid:25256302

9. Chakraborty A, Bera A, Mukherjee Aet al. Changing bacterial profile of Sundarbans, the world heritage mangrove:impact of anthropogenic interventions. World Journal of Microbiology and Biotechnology 2015;31: 593–610. pmid:25655378

10. Basak P, Pramanik A, Sengupta Set al. Bacterial diversity assessment of pristine mangrove microbial community from Dhulibhashani, Sundarbans using 16S rRNA gene tag sequencing. Genom Data 2016;7: 76–8. pmid:26981367

11. Ghosh A, Bhadury P. Insights into bacterioplankton community structure from Sundarbans mangrove ecoregion using Sanger and Illumina MiSeq sequencing approaches: A comparative analysis. Genom Data 2017;11: 39–42. pmid:27942457

12. Ghosh A, Bhadury P. Investigating monsoon and post-monsoon variabilities of bacterioplankton communities in a mangrove ecosystem. Environmental Science and Pollution Research 2018;25: 5722–39. pmid:29230649

13. Ghosh A, Bhadury P. Exploring biogeographic patterns of bacterioplankton communities across global estuaries. MicrobiologyOpen. 2019;8:e741.

14. Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol 2016;18: 1403–14. pmid:26271760

15. Das Gautam. (2016). Geomorphic Environments of the Thakuran River of the Sunderbans, India. Earth Science India. 9. 974–8350.

16. Grasshoff K, Kremling K, Ehrhardt M. Methods of seawater analysis. 3rd edn. Weinheim, Germany: Wiley-VCH, 1999.

17. Dhal P. K.; Sar P. Microbial communities in uranium mine tailings and mine water sediment from Jaduguda U mine, India: A culture independent analysis J. Environ. Sci. Health, Part A: Toxic/Hazard. Subst. Environ. Eng.2014, 49, 694– 709 pmid:24521415

18. Martin M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal, 17(1), pp. 10–12. https://doi.org/10.14806/ej.17.1.200

19. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014;30: 2114–20. pmid:24695404

20. Zhang J, Kobert K, Flouri T et al. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 2014;30: 614–20. pmid:24142950

21. Mahe F, Rognes T, Quince Cet al. Swarm v2: highly-scalable and high-resolution amplicon clustering. PeerJ 2015;3: e1420. pmid:26713226

22. Pruesse E, Peplies J, Glockner FO. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 2012;28: 1823–9. pmid:22556368

23. Kunda P, Dhal PK, Mukherjee A. Endophytic bacterial community of rice (Oryza sativa L.) from coastal saline zone of West Bengal: 16S rRNA gene based metagenomics approach. Meta Gene 2018;18: 79–86.

24. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. Journal of Statistical Software 2017;82.

25. Hill MO. Diversity and evenness: a unifying notation and its consequences. Ecology 1973;54: 427–32.

26. Hassenruck C, Fink A, Lichtschlag A et al. Quantification of the effects of ocean acidification on sediment microbial communities in the environment: the importance of ecosystem approaches. FEMS Microbiol Ecol 2016;92: fiw027. pmid:26887661

27. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the royal statistical society Series B (Methodological) 1995: 289–300.

28. Fernandes AD, Reid JN, Macklaim JM et al. Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2014;2: 15. pmid:24910773

29. Oksanen J, Blanchet FG, Kindt R et al. 2016. Vegan: Community Ecology Package.

30. Sarkar KS and Bhattacharya DB (2010) Water Quality analysis of the Coastal Region of Sunderban Mangrove Wetland, India, and Using Multivariate Statistical Techniques. Sciyo, Croatia, pp 258.

31. Balasubramanian R, Kannan L. Physico-chemical characteristics of the coral reef environs of the Gulf of Mannar Biosphere Reserve, India. Int J Ecol Environ Sci 2005;31: 273–8.

32. Sarkar SK, Frančišković-Bilinski S, Bhattacharya Aet al. Levels of elements in the surficial estuarine sediments of the Hugli River, northeast India and their environmental implications. Environment international 2004;30: 1089–98. pmid:15337354

33. Mukherjee D, Mukherjee A, Kumar B. Chemical fractionation of metals in freshly deposited marine estuarine sediments of sundarban ecosystem, India. Environmental geology 2009;58: 1757–67.

34. Ghosh A, Dey N, Bera A et al. Culture independent molecular analysis of bacterial communities in the mangrove sediment of Sundarban, India. Saline Systems 2010;6: 1. pmid:20163727

35. Mallick I, Bhattacharyya C, Mukherji S et al. Effective rhizoinoculation and biofilm formation by arsenic immobilizing halophilic plant growth promoting bacteria (PGPB) isolated from mangrove rhizosphere: A step towards arsenic rhizoremediation. Sci Total Environ 2018;610–611: 1239–50.

36. Nogales B., Lanfranconi M.P., Piña-Villalonga J.M., and Bosch R. Anthropogenic perturbations in marinemicrobial communities.FEMS Microbiol Rev2011;35:275–298. pmid:20738403

37. Borin S, Brusetti L, Daffonchio D, Delaney E & Baldi F. Biodiversity of prokaryotic communities in sediments of different sub-basins of the Venice lagoon. Res Microbiol2009;160:307–314 pmid:19416751

38. Buchan A, González JM, Moran MA. Overview of the marine Roseobacter lineage. Applied and environmental microbiology 2005;71: 5665–77. pmid:16204474

39. Simon M, Scheuner C, Meier-Kolthoff JP et al. Phylogenomics of Rhodobacteraceae reveals evolutionary adaptation to marine and non-marine habitats. The ISME journal 2017;11: 1483. pmid:28106881

40. Dang H, Zhang X, Sun J et al. Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea. Microbiology 2008;154: 2084–95. pmid:18599836

41. Doberva M, Stien D, Sorres J et al. Large Diversity and Original Structures of Acyl-Homoserine Lactones in Strain MOLA 401, a Marine Rhodobacteraceae Bacterium. Frontiers in microbiology 2017;8: 1152. pmid:28690598

42. Lauro F. M. et al. The genomic basis of trophic strategy in marine bacteria. Proc. Natl Acad. Sci. USA 106, 15527–15533 (2009) pmid:19805210

43. Jessen C, Lizcano JFV, Bayer T et al. In-situ effects of eutrophication and overfishing on physiology and bacterial diversity of the Red Sea coral Acropora hemprichii. PLoS One 2013;8: e62091. pmid:23630625

44. Kegler HF, Hassenrück C, Kegler P et al. Small tropical islands with dense human population: differences in water quality of near-shore waters are associated with distinct bacterial communities. PeerJ 2018;6: e4555. pmid:29761035

45. Teeling H, Fuchs BM, Becher D et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 2012;336: 608–11. pmid:22556258

46. Campbell AM, Fleisher J, Sinigalliano Cet al. Dynamics of marine bacterial community diversity of the coastal waters of the reefs, inlets, and wastewater outfalls of southeast F lorida. MicrobiologyOpen 2015;4: 390–408. pmid:25740409

47. Tanaka D, Tanaka S, Yamashiro Y et al. Distribution of oil‐degrading bacteria in coastal seawater, Toyama Bay, Japan. Environmental Toxicology: An International Journal 2008;23: 563–9.

48. Dupont CL, Rusch DB, Yooseph S et al. Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage. The ISME journal 2012;6: 1186. pmid:22170421

49. Seo J-H, Kang I, Yang S-J, Cho J-C. 2017. Characterization of spatial distribution of the bacterial community in the South Sea of Korea. PLoS One 12:e0174159. pmid:28306743

50. Adrian Ho, Paolo Di Lonardo D., Bodelier Paul L. E. Revisiting life strategy concepts in environmental microbial ecology, FEMS Microbiology Ecology, Volume 93, Issue 3, March 2017.

51. Nouioui I., Carro L., Garcia-Lopez M., Meier-Kolthoff J. P., Woyke T., Kyrpides N. C., et al. (2018). Genome-based taxonomic classification of the phylum Actinobacteria. Front. Microbiol. 9:2007. pmid:30186281

52. Azman A.-S., Othman I., Velu S. S., Chan K.-G., and Lee L.-H. (2015). Mangrove rare actinobacteria: taxonomy, natural compound, and discovery of bioactivity. Front. Microbiol. 6:856. pmid:26347734

53. Sengupta S, Pramanik A, Ghosh A, Bhattacharyya M (2015) Antimicrobial activities of actinomycetes isolated from unexplored regions of Sundarbans mangrove ecosystem. BMC Microbiol 15:1–16

54. Gong B., Cao H., Peng C. et al. High-throughput sequencing and analysis of microbial communities in the mangrove swamps along the coast of Beibu Gulf in Guangxi, China. Sci Rep 9, 9377 (2019) pmid:31253826

55. Manivasagan P, Kang KH, Sivakumar K, Li-Chan EC, Oh HM, Kim SK. Marine actinebacteria: an important source of bioactive natural products. Environ Toxicol Pharmacol 2014; 38(1): 172–88. pmid:24959957

56. Giovannoni SJ, Stingl U. Molecular diversity and ecology of microbial plankton. Nature 2005;437: 343. pmid:16163344

57. Morris RM, Frazar CD, Carlson CA. Basin‐scale patterns in the abundance of SAR11 subclades, marine Actinobacteria (OM1), members of the Roseobacter clade and OCS116 in the South Atlantic. Environmental microbiology 2012;14: 1133–44. pmid:22225975

58. Ngugi DK, Stingl U. High-Quality Draft Single-Cell Genome Sequence of the NS5 Marine Group from the Coastal Red Sea. Genome announcements 2018;6: e00565–18. pmid:29930069

59. Scott JJ, Budsberg KJ, Suen G, Wixon DL, Balser TC, Currie CR. Microbial community structure of leaf-cutter ant fungus gardens and refuse dumps. PLOS ONE. 2010;5(3):e9922 pmid:20360970

60. Lewin GR, Carlos C, Chevrette MG, Horn HA, McDonald BR, Stankey RJ, Fox BG, Currie CR. Evolution and ecology of actinobacteria and their bioenergy applications. Annu Rev Microbiol. 2016;70:235–54 pmid:27607553

61. Mizuno CM, Rodriguez-Valera F, Ghai R. Genomes of planktonic acidimicrobiales: widening horizons for marine actinobacteria by metagenomics. MBio 2015;6: e02083–14. pmid:25670777

62. Hoffmann K, Hassenruck C, Salman-Carvalho Vet al. Response of Bacterial Communities to Different Detritus Compositions in Arctic Deep-Sea Sediments. Front Microbiol 2017;8: 266. pmid:28286496

63. Ribicic D, Netzer R, Hazen TCet al. Microbial community and metagenome dynamics during biodegradation of dispersed oil reveals potential key-players in cold Norwegian seawater. Marine pollution bulletin 2018;129: 370–8. pmid:29680562

64. Dos Santos HF, Cury JC, Do Carmo FLet al. Mangrove bacterial diversity and the impact of oil contamination revealed by pyrosequencing: bacterial proxies for oil pollution. PLoS One 2011;6: e16943. pmid:21399677

65. Spring S, Scheuner C, Göker M et al. A taxonomic framework for emerging groups of ecologically important marine gammaproteobacteria based on the reconstruction of evolutionary relationships using genome-scale data. Frontiers in microbiology 2015;6: 281. pmid:25914684

66. Cho J-C, Giovannoni SJ. Cultivation and growth characteristics of a diverse group of oligotrophic marine Gammaproteobacteria. Applied and Environmental Microbiology 2004;70: 432–40. pmid:14711672

67. Yilmaz P, Yarza P, Rapp JZ et al. Expanding the world of marine bacterial and archaeal clades. Frontiers in microbiology 2016;6: 1524. pmid:26779174

68. Huggett MJ, Rappé MS. Genome sequence of strain HIMB30, a novel member of the marine Gammaproteobacteria. Journal of bacteriology 2012;194: 732–3. pmid:22247531

69. Zhou J, Richlen ML, Sehein TR et al. Microbial community structure and associations during a marine dinoflagellate bloom. Frontiers in microbiology 2018;9.

70. Haro-Moreno JM, Rodriguez-Valera F, López-García P, Moreira D, Martin-Cuadrado A-B. (2017). New insights into marine group III Euryarchaeota, from dark to light. ISME J 11: 1102–1117. pmid:28085158

71. Liu H, Zhang CL, Yang C, Chen S, Cao Z, Zhang Zet al. Marine Group II dominates planktonic Archaea in water column of the Northeastern South China Sea. Front Microbiol. 2017;8:1098 pmid:28663746