Indian Journal of Biotechnology Vol 15, April 2016, pp 201-209

Bacterial diversity in sediments of river Mahananda (Siliguri) as determined by 16S rRNA gene analysis

Shriparna Mukherjee1,2, Dhananjay Kumar1,3 and Ranadhir Chakraborty1* 1Omics Laboratory, Department of Biotechnology, University of North Bengal (NBU), PO NBU, Siliguri 734 013, India 2Department of Botany, P D Women’s College, Jalpaiguri 735 101, India 3Department of Biotechnology, Birla Institute of Technology, Mesra, Raanchi 835 215 India

Received 4 November 2014; revised 16 April 2015; accepted 10 May 2015

Analyses of 16S ribosomal (r) DNA clone library have enabled to assess the bacterial diversity in the sediment of an anthropogenically stressed river Mahananda at Siliguri, India. A total of 178 clones were examined by amplified rDNA restriction analysis (ARDRA) using HaeIII (four base specific restriction ). According to the restriction patterns generated from ARDRA, 71 clones were selected as operational taxonomic units (OTUs) and sequenced. No clones were identical to any of the known 16S rRNA sequences or to each other. Sequencing analysis revealed eight distinct major lineages of [α-, β-, γ-, ε-, Bacteroidetes, Firmicutes, Chloroflexi & Planctomycetes]. About 69.66% of the sequenced clones were grouped into Proteobacteria [comprising α (12.90%), β (45.96%), γ (2.41%) & ε (38.70%)] and 29.21% were included amongst the divisions other than Proteobacteria represented by Bacteroidetes (34.61%), Firmicutes (57.69%), Chloroflexi (3.84%) and Planctomycetes (3.84%). 16S rRNA gene phylotypes only distantly related to any of the previously identified sequences (non-affiliated rRNA genes) represent a significant fraction of the total cloned sequences. This is the first report describing microbial diversity of Mahananda river sediment by using metagenomic approach. The results suggest that the sediment of river Mahananda is the reservoir of unexplored microbial resources.

Keywords: ARDRA, 16S rDNA, microbial diversity, sediment

Introduction In the past, studies have been carried out to explore A special dominion in aquatic ecosystems is the microbial diversity in freshwater sediments of sediment. Compared to the corresponding water lakes, estuary and marine environments8-10. Despite bodies, the biomass and taxon richness of microbes such studies, data on the bacterial composition of are much higher in sediment part1. Due to their high the surface sediments of organically polluted/ abundance relative to the water column, these bacteria anthropogenically stressed river is very limited11. In play a key role in the decomposition, nutrient cycling one of our earlier reports, bacterial diversity of and carbon flux2. A complete knowledge regarding surface water of the river Mahananda of Siliguri, the resident microbial populations will, therefore, India was studied by analysis of 16S rRNA genes be required to understand how the chemical derived from environmental DNA12. In the present environment of sediment is generated and stabilized3. study, diversity of the benthic bacterial assemblages For centuries, rivers have been used as the dumping from sediments of the anthropogenically stressed river grounds for urban effluents, agricultural and industrial Mahananda was determined using 16S rRNA gene wastes that contain substances varying from simple clone library method. nutrients to highly toxic chemicals including heavy metals4-6. Several studies have been focused on the Materials and Methods adverse effects of different pollutants, fecal pathogens, organic and inorganic fertilizers on the Sediment Sampling, Analyses and Site Climate 7 The river Mahananda flows through the middle of river ecosystem . the Siliguri city. The maximum river gorge below the ——————— main Mahananda bridge is about 50 m wide (varies *Author for correspondence: from 10-50 m in different seasons). Sediment samples Tel: +91+353_2776354; Fax: +91-353-2699001 Mobile: +91-9434872273 were collected from both the banks and mid river. [email protected] An ethanol flamed spatula was used to collect soil 202 INDIAN J BIOTECHNOL, APRIL 2016

samples from the 2 to 10 cm of soil from various sequencing was performed with the ABI PRISMTM points. Collected soil samples were brought to the Dye Terminator Cycle Sequencing Ready Reaction laboratory in properly labeled, autoclaved and sealed kit and the reactions were analyzed on an ABI PRISM polyethylene bags on ice. Prior to the total DNA 377 DNA sequencer. The recombinant plasmids isolation, all the soil samples collected from different were sequenced with T7 and SP6 primers from both spots were mixed to homogeneity to make a single ends and the middle portion was sequenced with composite sample. Textural composition of the 530F primer. Only the sequences with overlapping sediment samples (determination of sand, silt & clay) reads were included in further analysis. was not determined. Phylogeny and Statistical Analysis Isolation of Total Soil DNA A CHIMERA_CHECK program of the Ribosomal Total soil (sediment) DNA was isolated using a Database Project II (RDP II) was used to detect the modified CTAB-SDS based DNA extraction 14 13 presence of possible chimeric sequences . In order to technique . RNase (10 mg/mL) treatment, followed find closely related sequence(s) and to obtain a by phenol extraction was carried out prior to PCR preliminary phylogenetic affiliation of the clones, amplification of 16S rRNA gene sequences. submissions of unaligned sequences were made to PCR Amplification and 16S rRNA gene Library Construction the sequence match program of the RDP and to the PCR was performed in 50 L reaction mixtures advanced BLAST search program of NCBI. The most (40 M deoxynucleoside triphosphates, 12 pMol of similar reference sequences were retrieved and each primer, 1 U of Taq DNA polymerase [Promega], aligned with clone sequences using Clustal_X15. 1× PCR buffer [Promega] for each DNA sample, Phylogenetic trees were constructed using MEGA using universal bacterial primer 27f (5′- version 4.0 by the neighbor-joining algorithm and the AGAGTTTGATCCTGGCTCAG-3′) and 1492r (5′- Jukes-Cantor distance estimation method with TACGGTTACCTTGTTACGACTT-3′). Maintenance bootstrap analyses for 1000 replicates16. of PCR condition and construction of 16S rRNA gene clone library were made according to the methods Nucleotide Sequence Accession Numbers described earlier12. The accession numbers for sequences determined from sediment clone library are as follows: FR690816 ARDRA Method to FR690827, FR691000 to FR691009, FR691051, About 178 clones from sediment clone library were FR691485 to FR691531 and FR729922. chosen for ARDRA analysis. The plasmid inserts were PCR amplified by using vector-specific primers Results SP6 (5′-GATTTAGGTGACACTATAG-3′) and T7

(5′-TAATACGACTCACTATAGGG-3′). All the DNA Extraction, Library Construction and Sequencing PCR assays were carried out in a total volume of analysis 50 L mixtures ((40 M deoxynucleoside triphosphates, The total DNA was extracted from the sediment of 12 pMol of each primer, 1 U of Taq DNA polymerase river Mahananda using modified CTAB-SDS based [Promega], 1× PCR buffer [Promega] for each DNA DNA isolation technique. One clone library was sample). One cycle of PCR was run at 94°C for 5 min, established. The recombinant clones were selected followed by 25 cycles each comprised of 94°C for based on α-complementation (blue-white screening) 30 sec, 47°C for 30 sec, 72°C for 2 min. A final technique and also confirmed by the re-PCR analysis extension was performed at 72°C for 7 min. PCR and restriction enzyme digestion. Sequencing analyses products were digested with HaeIII in accordance included 71 clones. One nonbacterial chloroplastic with the manufacturer’s instructions, electrophoresed DNA was found amongst the recombinant clones. In on a 2% (w/v, solution of agarose in TAE buffer) this study, no primer pair was chosen to amplify the agarose gel, and stained with ethidium bromide. archaeal 16S rRNA gene sequences and, therefore, The band sizes were determined by using 100 bp did not get any archaeal sequences. Our phylogenetic DNA ladder as size standard. analysis revealed that 70 clones from the sediment

Sequencing and clone library fell into 5 major lineages of the The clones corresponding to the unique banding domain Bacteria: Proteobacteria, Bacteroidetes, pattern were used for sequencing. Nucleotide Firmicutes, Chloroflexi and Planctomycetes. MUKHERJEE et al: BACTERIAL DIVERSITY IN MAHANANDA RIVER SEDIMENT 203

Proteobacteria α-Proteobacteria Phylum Proteobacteria was represented by 69.66% The subclass ∝- Proteobacteria was represented by of the total sequences. Proteobacterial sequences fell 9 OTUs (16 clones). The cloned sequences belong to within 4 different classes: α-Proteobacteria (12.90%), orders Sphingomonadales, Rhizobiales, Rhodobacterales β-Proteobacteria (45.96%), γ-Proteobacteria (2.41%) and Caulobacterales. The majority of the sequences and ε-Proteobacteria (38.70%). belonged to Sphingomonadales (62.5%). The sequences

β-Proteobacteria belonging to Rhodobacterales were 18.75%, while It represented the most dominant subclass of Rhizobiales and Caulobacterales accounted for Proteobacteria comprising of 20 OTUs (consisting of 12.5 and 6.25%, respectively. OTU MS004 was 57 clones) forming four clusters in the phylogenetic related (98%) to uncultured bacterial clone tree. The clusters included four orders, Burkholderiales (DQ337059) isolated from the subsurface water (represented by 52.63% of the clones), Rhodocyclales of the Kalahari shield, Africa [GenBank description]. (24.56%), Neisseriales (19.29%) and Methylophilales In the phylogenetic tree MS004 branched deeply (3.50%). Sequences affiliated to cluster I were related with Novosphingobium aromaticivorans (U20756). to the cultured representatives of Burkholderiales MS003 and MS009, most predominating sequence (Hydrogenophaga, , Acidovorax, types represented by five clones, were related to Comamonas, Mitsuaria & Herbaspirillum). A sub (98%) uncultured bacterium clone S0022 (FJ820395) cluster of four OTUs (MS014, MS017, MS018 & recovered from a hypertrophic freshwater lake MS028) did not branch with any of the cultured [GenBank description]. MS005 was most similar to representatives. They shared 99% similarity with an Pleomorphomonas oryzae (99% similarity) and uncultured bacterium clone (AY212659) isolated MS006 was similar to Phenylobacterium falsum (98% from water 20 m downstream of manure (GenBank similarity). These relationships were also supported description). MS022 related to (99% similarity) by a very high bootstrap confidence (100%) value on uncultured Hydrogenophaga sp. clone DS062 the phylogenetic tree (Fig. 1). MS002 and MS008 (GenBank description) but in the phylogenetic tree clustered with Rhodobacter changlensis. clustered with Macromonas bipunctata (strict anaerobic Gram-negative heterotroph), which is a closest γ-Proteobacteria representative of the genus Hydrogenophaga. MS026 Only two OTUs (MS031 & MS032) belong to this branched with Zoogloea ramigera (97% similarity) by subclass and exhibited 97-98% similarity to 100% bootstrap confidence value. It was found that Thermomonas brevis (AB355702) previously isolated the two OTUs, MS029 (FR691489) and MS030 from activated sludge18. (FR691490) clustered with the terrestrial strain Methylobacillus pratensis (an obligately methylotrophic Other Bacteria non-methane utilizing bacterium)17 (Fig. 1). Only 29.21% of the sequenced clones were

ε-Proteobacteria included in the divisions other than Proteobacteria. The second most abundant fraction within Proteobacteria was represented by the members of Bacteroidetes Epsilonproteobacteria. About 38.70% of the cloned This fermentative group was represented by 12 sequence types have shown similarity to the OTUs with 18 clones. Sequences formed three epsilonproteobacterial sequences in the database. clusters in the phylogenetic tree. Sequences affiliated Sequences were phylogenetically closely related to the to cluster I branched with Cloacibacterium normanense, originally isolated from treated clones retrieved from waste water (GQ844389), heavy 19 metal contaminated estuarine sediment (HQ132428), municipal waste water . Sequence affiliated to uncultured Arcobacter sp. clone (AY692045) isolated Cluster II branched with Leadbetterella byssophila 20 from anaerobic waste water, and Sulfuricurvum sp. originally isolated from cotton waste compost . enrichment culture clone (EU498374) isolated from Sequences affiliated to cluster III were closely related chloroethane degrading culture [GenBank description]. to the uncultured bacterium clones isolated from lake It is important to note that all the sequenced clones water (FJ612126), microbial fuel cells (FJ375426), exhibited close resemblance with the uncultured wastewater (HM445933), tar-oil contaminated aquifer representatives as also evident in the phylogenetic tree sediments (EU266902) and other environmental (Fig. 1). sources (AJ306737) [GenBank description] (Fig. 2). 204 INDIAN J BIOTECHNOL, APRIL 2016

Fig 1—Neighbor-joining tree showing positions of proteobacterial sequences in sediment clone library including reference sequences obtained from GenBank. [Methanothermococcteus okinawensis is used as an outgroup. GenBank accession numbers are in parentheses.] MUKHERJEE et al: BACTERIAL DIVERSITY IN MAHANANDA RIVER SEDIMENT 205

Fig 1 (Contd.)—Neighbor-joining tree showing positions of proteobacterial sequences in sediment clone library including reference sequences obtained from GenBank. [Methanothermococcteus okinawensis is used as an outgroup. GenBank accession numbers are in parentheses.]

206 INDIAN J BIOTECHNOL, APRIL 2016

Fig 2—Neighbor-joining tree showing positions of non-proteobacterial sequences in sediment clone library including reference sequences obtained from GenBank. [Methanothermococcteus okinawensis is used as an outgroup. GenBank accession numbers are in parentheses.]

Firmicutes (FJ462045), bacterial consortia (AJ488074) etc. Firmicutes were second only to β-Proteobacteria [GenBank description] and in the phylogenetic tree in abundance in this clone library comprising of they formed a separate cluster (Fig. 2). 17 OTUs (32 clones). Sequences fell within this phylum were grouped with both uncultured clones Planctomycetes and cultured representatives. Most of the OTUs Two OTUs were group within the Planctomycetes. branched with Clostridium spp. MS033, MS042, A FASTA search revealed a close relationship (96% MS043, MSO44, MS049, MS068 were closely related similarity) of MS062 to an uncultured bacterium to uncultured bacterial clones retrieved from clone (FJ849362) isolated from arctic stream wastewater sludge (FJ534955), continental crust sediment [GenBank description] and MS063 to (DQ088761), effluent from the chemical industry an uncultured bacterium clone (HM129983) MUKHERJEE et al: BACTERIAL DIVERSITY IN MAHANANDA RIVER SEDIMENT 207

Fig 2 (Contd.)—Neighbor-joining tree showing positions of non-proteobacterial sequences in sediment clone library including reference sequences obtained from GenBank. [Methanothermococcteus okinawensis is used as an outgroup. GenBank accession numbers are in parentheses.].

(97% similarity) isolated from a Tibetan lake Discussion [GenBank description]. Phylogenetic analysis showed This is the first report describing the bacterial that these two OTUs branched deeply with Pirellula diversity in part using 16S rRNA gene clone library staleyi and Isosphaera pallida with very high method in sediments of river Mahananda. The present bootstrap value (Fig. 2). study has obtained a primary insight into the major bacterial populations in Mahananda river sediment Chloroflexi that are only qualitative and not quantitative. The Only one OTU, MS071, with two clones in analysis revealed the presence of phylogenetically sediment library, showed sequence similarity with an diverse bacterial populations in this ecological niche. uncultured Chloroflexi bacterium (EU016438) [GenBank Five major phyla of the bacterial domain were description] and branched deeply with Dehalococcoides detected. All these are the phyla commonly detected with a 99% bootstrap confidence value. in various freshwater sediment environments10,21. 208 INDIAN J BIOTECHNOL, APRIL 2016

Newton et al22 suggested that generally most 2 Bai Y, Shi Q, Wen D, Li Z, Jefferson W A et al, Bacterial members of the phyla Proteobacteria (α & β) and communities in the sediments of Dianchi lake, a partitioned eutrophic waterbody in China, PLoS One, 7 (2012) e37796. Bacteroidetes are restricted to only freshwater 3 Nealson K H, Sediment bacteria: Who’s there, what are they ecosystems and thus represent native organisms. doing, and what’s new? Annu Rev Earth Planet Sci, Several studies have shown the dominance of 25 (1997) 403-434. Proteobacteria in most of the sediments2,9,10. 4 Weber R, Gaus C, Tysklind M, Johnston P, Forter M et al, Comparable results were achieved by analyzing Dioxin- and POP-contaminated sites—Contemporary and future relevance and challenges: Overview on background, sediment clone library of Mahananda River. In the aims and scope of the series, Environ Sci Pollut Res Int, present study, majority of bacterial sequences 15 (2008) 363-393. obtained belonged to the phylum Proteobacteria 5 Schulz-Zunkel C & Krueger F, Trace metal dynamics in which represented 69.66% of the total sequences. As floodplain soils of the river Elbe: A review, J Environ Qual, most Proteobacteria are known to easily metabolize 38 (2009) 1349-1362. 21 6 Skoulikidis N T, The environmental state of rivers in the degradable organic substrates , their dominance in Balkans—A review within the DPSIR framework, Sci Total the sediment of this eutrophic and heavily polluted Environ, 407 (2009) 2501-2516. river is not surprising. The γ-proteobacterial members 7 van der Zaan B, Smidt H, de Vos W M, Rijnaarts H & are less abundant in freshwater environments. Among Gerritse J, Stability of the total and functional microbial the non-proteobacterial members, Bacteroidetes communities in river Sediment mesocosms exposed to anthropogenic disturbances, FEMS Microbiol Ecol, 74 (34.61%) and Firmicutes (57.69%) are found (2010) 72-82. more abundant compared to Planctomycetes (3.84%) 8 Ghosh A, Dey N, Bera A, Tiwari A, Sathyaniranjan K B et and Chloroflexi (3.84%). Members of the phylum al, Culture independent molecular analysis of bacterial Bacteroidetes and Firmicutes are apparently communities in the mangrove sediment of Sundarban, India, ubiquitous in freshwater environments. Frequent Saline System, 6 (2010) 1-11. 9 Haller L, Tonolla M, Zopfi J, Peduzzi R, Wildi W et al, occurrence of Chloroflexi, or green nonsulfur bacteria Composition of bacterial and archaeal communities in 23 are recorded in deep marine sediments , sulfide and freshwater sediments with different contamination levels (Lake methane-rich cold seep sediments24, vermicompost 25 Geneva, Switzerland), Water Res, 45 (2011) 1213-1228. and fresh water sponge26. In this study, only one 10 Tamaki H, Sekiguchi Y, Hanada S, Nakamura K, Nomura N et al, Comparative analysis of bacterial diversity in OTU, MS071, branched deeply with Dehalococcoides freshwater sediment of a shallow eutrophic lake by molecular with a 99% bootstrap confidence value. Studies reveal and improved cultivation-based techniques, Appl Environ the future promise for Dehalococcoides like species Microbiol, 71 (2005) 2162-2169. in bioremediation applications27. Genomic and 11 Umamaheswari S & Saravanan N A, Genetic composition of proteomic studies made on Planctomycetes isolates benthic bacterial assemblages in the river Cauvery, Tamil Nadu, India, Int J Lakes Rivers, 3 (2010) 111-118. revealed their ability to degrade phytoplankton- 28 12 Mukherjee S, Kumar D, Nanda A K & Chakraborty R, 16S derived carbohydrates . Recent comparative genomic rRNA gene sequence analyses of the metagenome derived and metagenomic analyses of a few representatives of from waters of river Mahananda at Siliguri: An approach to Planctomycetes offered new insights with major understand bacterial diversity, Indian J Biotechnol, 12 (2013) medical, ecological and biotechnological relevance 29. 80-87. 13 Zhou J, Bruns M A & Tiedje J M, DNA recovery from soils Although analysis of only 71 clones is not apt to of diverse composition, Appl Environ Microbiol, 62 (1996) infer bacterial diversity in totality, yet it has revealed 316-322. that the sediment of the river Mahananda harbored a 14 Maidak B L, Cole J R, Parker C T Jr, Garrity G M, Larsen N phylogenetically diverse population of organisms et al, A new version of the RDP (Ribosomal Database from bacterial domain. Since the recovery of gene Project), Nucl Acids Res, 27 (1999) 171-173. sequence information from DNA templates does not 15 Thompson J D, Gibson T J, Plewniak F, Jeanmougin F & Higgins D G, The Clustal_X windows interface: Flexible necessarily indicate in-situ cellular activity, application strategies for multiple sequence alignment aided by quality of metatranscriptomics in future endeavors shall lead analysis tools, Nucleic Acids Res, 25 (1997) 4876-4882. to more in-depth knowledge about the structure- 16 Kumar S, Tamura K & Nei M, MEGA 3: Integrated software function component of the bacterial community. for molecular evolutionary genetics analysis and sequence alignment, Brief Bioinform, 5 (2004) 150-163. References 17 Doronina N V, Trotsenko Y A, Kolganova T V, Tourova T P 1 Zinger L, Amaral-Zettler L A, Fuhrman J A, Horner-Devine & Salkinoja-Salonen M S, Methylobacillus pratensis sp. M C, Huse S M et al, Global patterns of bacterial beta- nov., a novel non-pigmented, aerobic, obligately diversity in seafloor and seawater ecosystems, PLoS One, methylotrophic bacterium isolated from meadow grass, Int J 6 (2011) e24570. Syst Evol Microbiol, 54 (2004) 1453-1457. MUKHERJEE et al: BACTERIAL DIVERSITY IN MAHANANDA RIVER SEDIMENT 209

18 Phuong K, Kakii K & Nikata T, Intergeneric coaggregation 24 Reed A J, Lutz R A & Vetriani C, Vertical distribution and of non-flocculating Acinetobacter spp. isolates with other diversity of bacteria and archaea in sulfide and methane-rich sludge-constituting bacteria, J Biosci Bioeng, 107 (2009) cold seep sediments located at the base of the Florida 394-400. Escarpment, Extremophiles, 10 (2006) 199-211. 19 Allen T D, Lawson P A, Collins M D, Falsen E & Tanner R 25 Fracchia L, Dohrmann A B, Martinotti M G & Tebbe C C, S, Cloacibacterium normanense gen. nov., sp. nov., a novel Bacterial diversity in a finished compost and vermicompost: bacterium in the family Flavobacteriaceae isolated from Differences revealed by cultivation-independent analyses of municipal waste water, Int J Syst Evol Microbiol, 56 (2006) PCR-amplified 16S rRNA genes, Appl Microbiol Biotechnol, 1311-1316. 71 (2006) 942-952. 26 Schmitt S, Denies P, Behnam F, Wagner M & Taylor M W, 20 Weon H Y, Kim B Y, Kwon S W, Park I C, Cha I B et al, Leadbetterella byssophila gen. nov., sp. nov., isolated from Chloroflexi bacteria are more diverse, abundant and similar cotton-waste composts for the cultivation of oyster in high than in low microbial abundance sponges, FEMS mushroom, Int J Syst Evol Microbiol, 55 (2005) 2297-2302. Microbiol Ecol, 78 (2011) 497-510. 27 Fennell D E, Nijenhuis I, Wilson S F, Zinder S H & 21 Wilms R, Köpke B, Sass H, Chang T S, Cypionka H et al, Häggblom M M, Dehalococcoides ethenogenes strain 195 Deep biosphere-related bacteria within the subsurface of tidal reductively dechlorinates diverse chlorinated aromatic flat sediments, Appl Environ Microbiol, 8 (2006) 709-719. pollutants, Environ Sci Technol, 38 (2004) 2075-2081. 22 Newton R J, Jones S E, Eiler A, McMohan K D & Bertilsson 28 Rabus R, Gade D, Helbig R, Bauer M, Glöckner F O et al, S, A guide to the natural history of freshwater lake bacteria, Analysis of N-acetylglucosamine metabolism in the marine Microbiol Mol Biol Rev, 75 (2012) 14-49. bacterium Pirellula sp. strain 1 by a proteomic approach, 23 Inagaki F, Nunoura T, Nakagawa S, Teske A, Lever M et al, Proteomics, 2 (2002) 649-655. Biogeographical distribution and diversity of microbes in 29 Wagner M & Horn M, The Planctomycetes, Verrumicrobia, methane hydrate-bearing deep marine sediments on the Chlamydiae and sister phyla comprise a superphylum with Pacific Ocean Margin, Proc Natl Acad Sci USA, 103 (2006) biotechnological and medical relevance, Curr Opin 2815-2820. Biotechnol, 17 (2006) 241-249.