Bacterial diversity in the sediments of a temperate artificial lake, Rapel reservoir

Rapel reservoir is an eutrophic system in Chile that has undergone a series of anthropogenic impacts in recent decades. To provide a better understanding of the processes occurring in the reservoir, we examined bacterial composition from surface sediments using traditional microbiology techniques and molecular biology tools. Our results showed significant temporal variation in the physical and chemical composition of the water column, but no depth-related differences during the study period. To detect temporal changes in bacterial composition, cultivable heterotrophic bacteria, heterotrophic iron oxidizing bacteria, and sulfate-reducing bacteria were extracted from the surface sediments and their concentration measured. Microbial diversity in sediments was represented by closest relatives of eight different bacterial phyla. The most frequently recovered phylotypes in the clone library of 16S rDNA were related to sulfate-reducing bacteria belong to the Deltaproteobacteria group.


INTRODUCTION
Rapel reservoir is located in Central Chile (34º S, 71° W) at an altitude of 130 m.This mesotrophic temperate lake was created during 1968 for hydroelectric power generation.The lake has three sub-basins with different physical, chemical and hydrodynamic characteristics (Contreras et al. 1994).The shallowest zone corresponds to Alhué basin (Fig. 1), which contrasts with the other subbasins as it receives a low discharge (< 1 m 3 s -1 ) but sulfate rich (1,500 mg L -1 ) inflowing stream that carries mine-water effluents.Rapel reservoir is shallow (mean depth = 6-8 m), and wind-effects ensure that the water column is permanently mixed and oxygenated at alkaline pH values (Vila et al. 1997(Vila et al. , 2000)).Total mixing of the water column during most of the year is similar to many Chilean lakes (Vila et al. 2000, Soto 2002).The sediments in this system have a high concentration of reduced chemicals, and gas production is high with CH 4 and N 2 supersaturation in the surface layers (Adams et al. 2000).
A series of anthropogenic activities have produced significant changes in composition and abundance of phytoplankton communities and in the water quality in the lake (Vila & Pardo 2003).Several hypotheses could explain the change of the phytoplankton composition in this reservoir.The most accepted are high sulfate levels in the water column (Norici et al. 2005) and possible metal solubilization from the sediments (Vila et al. 2000).As it has been cited for reservoirs with different trophic state, the bacterial activity in the sediments has been linked to the trophic conditions in the water body due to sedimentation of particulate organic matter (Wobus et al. 2003).In Chilean freshwater ecosystems, microbial diversity has been examined as a means of assessing water quality (Miranda & Castillo 1998).The role of microorganisms in sediment gas production under anoxia or low oxygen conditions is well known.Most work on the microbiology of sediments has been undergone through measurement of decomposition processes and the isolation of microorganisms (Spring et al. 2000).This approach is limited because only between 0.001 to 10 % of the total number of visible cells can be retrieved by isolation (Amann et al. 1995).The introduction of molecular methods based on PCR and sequencing of small-subunit ribosomal RNA genes made it possible to estimate the phylogenetic diversity without culture.A combination of different techniques would allow a better understanding of the complexity of bacterial interactions.
Previous research in Rapel reservoir has focused on limnological aspects (Montecino & Cabrera 1982, Reynolds et al. 1985, Vila et al. 1987, 1997), and until now, the bacterial diversity in this reservoir has not been examined.To identify and elucidate the possible processes that occur in the interphase sedimentswater of the reservoir, bacterial diversity were determinate in the surface sediments of the Alhué basin by using traditional microbial techniques and molecular approaches, namely analysis of 16S rDNA by PCR-DGGE and clone library.At the same time, we measured temporal variation in a series of physico-chemical parameters.We compared the results of clone libraries from the Alhué site which received mine-water discharges with El Manzano, an unimpacted site.

Sample analysis
Water and sediment samples for physical and chemical analyses were obtained monthly between January and June 2002 from two sites (Alhué and El Manzano) at the Alhué basin of Rapel reservoir (Fig. 1).These two stations were selected because Alhué directly receives mining wastes, and El Manzano does not.Water samples were collected with a Van Dorn 3 L bottle at depths of 0.3, 6 and 8 m (maximum depth).Temperature was recorded with a digital Hanna HI thermometer, pH with a Hanna HI 8314 pH meter, and conductivity with a VWR conductivity meter.Dissolved oxygen was measured with the Winkler method, BOD 5 , P-PO 4 3-, N-NO 3 -, S-SO 4 2-and chlorophyll a were analyzed according to Standard Methods (APHA 1998).Data were analyzed with two-factor ANOVA and Kendall Tau non parametric correlation (P < 0.05) (Zar 1996).
Samples were obtained from the upper 15 cm of the lake sediments using a Phleger core sampler and the sediment was homogenized for subsequent DNA extraction.For DNA analyses, we used sediment samples taken from Alhué (Rap1) and El Manzano sites (Rap2) in July 2003 (Fig. 1).

Metals
Total iron and copper in the sediments were analyzed by atomic absorption spectroscopy.A sediments microscopic mineralogical composition analysis was conducted following Siegel (1999).

DNA extraction
Aliquots of 0.25 g of sediments were used for DNA extraction using Ultra Clean Soil DNA Isolation Kits (MoBio Lab., Inc., Solana Beach, California, USA).DNA was dissolved in 10 mM Tris-HCl buffer and stored at -18 ºC.DNA quality was checked by 0.8 % agarose gel electrophoresis and staining with ethidium bromide.

PCR amplification
Oligonucleotide primers Eub9-27F and Eub1542R were used to amplify eubacterial 16S rDNA (Brosius et al. 1981).Each PCR reaction contained 10X PCR-buffer with 2 mM MgCl 2 (Roche), 200 mM dNTP mixture (Gibco), 1 pmol of each primer, 2.5 U Taq DNA polymerase (Roche), 10-100 ng template DNA and sterile HPLC grade water to a final volumen of 50 μl.The PCR was performed using the following conditions: initial denaturing step of 5 min at 94 ºC followed by 34 cycles of denaturing at 94 ºC for 30 sec, annealing at 40 ºC for 45 sec and elongation at 72 ºC for 1.5 min.An annealing temperature at 40 °C was used in order to maximise the number and diversity of bacteria recovered.

DGGE analysis
DGGE was performed according to Muyzer et al (1993).The fragments were amplified with the primer combination P2-P3.PCR products were applied onto 7.5 % polyacrylamide gels containing a linear gradient of 30-60 % denaturant where 100 % denaturant was defined as 7 M urea and 40 % deionized formamide.Electrophoresis was carried out in a BioRad D Gene System (BioRad, Munich, Germany) at 60 °C, 200 V for 6 h.The gels were stained with silver nitrate.

Cloning and 16S rRNA gene sequence analysis
Purified amplicons prepared with proof reading DNA polymerase (pfu) were cloned into Escherichia coli using the pCR-Blunt vector (Invitrogen, San Diego, California, USA) according to the manufacturer instructions.Clones with correct rDNA insertion were initially identified using whole-cell PCR with vector-specific M13F/R primers and checked for appropriate length by agarose gel electrophoresis.For sequencing, the PCR products of correct size were purified in multiscreen cleanup plates (Millipore) and quantified after agarose electrophoresis and ethidium bromide staining via the Molecular Analyst package (BioRad).To screen clones, cycle sequencing was performed with primer 341F (Muyzer et al 1993) using the BigDye Terminator Cycle Sequencing Kit v3.1 and analyzed on an automatic capillary sequencer (model 3100 Gene Analyzer, Applied Biosystems).

Phylogenetic analysis
The 16S rDNA sequences were analyzed in BLAST search (http://www.ncbi.nlm.nih.gov/blast) in order to determine the closest relatives in the database.Sequences were aligned using CLUSTALX v1.81 (Thompson et al. 1997).Phylogenetic relationships of the sequences were calculated using distance matrix methods (Neighbor-joining) in MEGA 3.1 (Kumar et al. 2004), with 1000 bootstrap permutations.The reference strains and sequences were obtained via BLAST searches in GenBank.Sequences with similarities ≥ 99 % were considered to represent the same phylotype.Rarefaction curves (Simberloff 1972) were determined by RARFAC program (http://www.icbm.de/pmbio/downlist.htm).

Nucleotide sequence accession numbers
The nucleotide sequences of bacterial 16S rDNA sequences have been deposited in GenBank with accession numbers: EF192877-EF192915.

Water column characteristics
Table 1 details the physic and chemical characteristics of Alhué basin water column.No significant differences were found between samples collected from surface and deeper waters (two-way ANOVA without replication).Therefore, the results shown in Table 1 are mean values, with standard deviations from three depths.Temporal variations show the typical decrease of water temperature between summer and winter.Only N-NO 3 -concentrations showed significant differences among the sampling periods (F 1 , = 168.07,P < 0.0001), and there were no significant temporal variation in any of the other chemical variables (two-way ANOVA without replication).Both pH and conductivity increased with time.Oxygen saturation was below 100 % during January, March and June at all depths (Table 1).The change in chemical characteristics in the water column during June coincided with heavy rain and inflow versus outflow and manipulation of lake levels.
The concentration ranges of the nutrients N-NO 3 -, P-PO 4 3-and chlorophyll a corresponded to a meso-eutrophic lake (Ryding & Rast 1992) (Table 1).There was no difference between ammonium values recorded in surface and deeper waters.Sulfate values were homogenous in the water column.

Bacteria of the sediment water interphase
Fig. 2 shows the heterotrophic, iron oxidizing and sulfate-reducing bacteria concentration found during the study period.The level of heterotrophic bacteria is similar compared with other meso-eutrophic systems (Wobus et al. 2003) and with low seasonal variation (6.9x10 6 -9.7x10 6 CFU g dry weight -1 ).
Generally, heterotrophic iron oxidizing bacteria were more abundant during autumn, and showed a general decrease between January and April, but subsequently increased by three orders of magnitude in May (3.1x10 3 -3.7x10 6CFU g dw -1 ).It is important to note the presence of sulfate-reducing bacteria in the upper sections of the sediment.MPN values that fluctuated between 10 2 and 10 4 MPN g dw -1 , and showed a positive significant correlation between sulfate content in the water column and sulfate reducing bacteria found in sediments (Kendall Tau, n = 5, Z = 2.45, P < 0.01).3).The sequences shared similarity values ranging between 87 % and 95 % with cultured species, the higher values corresponding to those taxa most closely related to Deltaproteobacteria.

Clone library and bacterial diversity
Ninety-six clones per sample were arbitrarily selected for analysis and fifty-five sequences with the correct insert were recovered.Thirtyfive sequences were used for phylogenetic tree construction.Each sequence was used as representative clones and reference organisms were selected as the closest relatives to samples sequences from the database.Rarefaction curves were used to determine if a sufficient number of clones were screened to estimate total diversity in each clone library.A combined rarefaction curve for both Alhué and El Manzano sites reached an asymptotic maximum (Fig. 4).We found clones related to eight different phylogenetic groups.The majority of sequences (67 %) were grouped to the Proteobacteria phylum dominated by the Deltaproteobacterium group (24 %), (Table 2), which in this case corresponds to sulfatereducing bacteria.Other Proteobacteria represented in the clone library were Alphaproteobacteria (12 %), Gammaproteobacteria (19 %) and Betaproteobacteria (11 %).Further, Cytophaga-Flavobacteria-Bacteroidetes (CFB) (12 %), Firmicutes (5 %) Verrumicrobia (4 %) were present in the clone library (Fig. 5).Twelfth percent of the sequences did not have a clear classification inside the phylum Bacteria.Only 27 % percent of sequences shared more than 97 % sequence identity with the first entry in the BLAST search, the rest of the sequences had more than 90 % similarity.Sixty-two percent of the sequences were classified in the BLAST search as uncultured bacteria, which means that most of the sequences does not have similarity with already cultured bacteria.

DISCUSSION
The significant physic and chemical seasonal differences reported here were similar to those previously reported for Rapel reservoir (Vila et al. 2000).The concentrations of heterotrophic bacteria we describe correspond to those previously described in the literature for reservoirs in the process of eutrophication (Wobus et al. 2003).This is indicative of important biological activity at the sediments of the Alhué basin; a gradual increase in bacterial abundance was observed between January and June and, this coincides with higher chlorophyll a values during the same period.
The presence of iron oxidizing bacteria in our samples may be indicative of the presence of Fe (II) associated with organic matter, and which could be oxidized to Fe (III).In spite of the low availability of degradable organic matter found (deduced from BOD 5 ), we have shown that there are high concentrations of iron in the sediments compared to other metals like Mn, Cu and Zn.According to metals analysis, iron is associated to metals and other elements (Pizarro et al. 2003).Higher values of iron oxidizing bacteria during May and June coincided with the reservoir volume increase and with the opening of the floodgates, due to heavy rain, which would imply particulate material incorporation and consequently an increase in the availability of organic matter.
Recently, it has been demonstrated that sulfate-reducing bacteria are present in oxygenated and anoxygenic environments, especially those close to surface sediments (Cypionka 2000, Dolla et al. 2006).This might explain the results of this study, e.g. the presence of SRB in sediments and water column.Cypionka (2000) noted that it is possible to detect two equivalent peaks of SRB at different depths, one in anoxic zones and the other at the water sediment interphase.DGGE analyses demonstrated a greater number of bands at the El Manzano site, but few qualitative differences with the Alhué sites.The clone library confirmed this result, as the number of phylotypes was higher in El Manzano site than at Alhué.In total, eight different Phyla of the Bacteria domain were found, with the sulfate-reducing bacteria group the most abundant.This group has an important role in heavy metal precipitation in sediments, which is dependent on pH and redox potential.They participate in other biogeochemical cycles including phosphite oxidation (Schink & Friederich 2000, Schink et al. 2002), reduction of Fe(III) (Coleman et al. 1993) and the precipitation of arsenic (Rittle et al. 1995, Kirk et al. 2004).Mining activities in the catchment area of Rapel reservoir have increased the amounts of sulfate and heavy metals entering the lake due to inflow of clear water tailings (Vila et al. 2000).Kirk et al. (2004) emphasize the role of sulfate-reducing bacteria in freshwater systems because in the absence of sulfate reduction, methanogenesis is the dominant type of microbial metabolism, and arsenic accumulates to high levels.In Rapel reservoir the flux of CH 4 from the sediments is high (Adams et al. 2000) but oxygen in the water column prevents bubble formation, and CO 2 is produced.Currently, nothing is known regarding the flux or production rates of H 2 S in the Rapel Reservoir.It may be that sulfate reduction might limit heavy metal solubilization, which would explain the low amounts of Cu, Zn and other metals in the water column (Vila et al. 2000) in spite of the constant inflow of mine tailings.
Desulfotignum phosphitoxidans, a sulfate reducing bacterium isolated from marine sediments that utilises phosphite as its sole electron donor and CO 2 its only carbon source (Schink et al. 2002), was the closest relative for most of the clones in the library.Some clones were related to Geobacter and Desulfobacterium.It has been described that this kind of bacteria are able to reduce Fe (III) under anaerobic conditions.This process is one of the most important geochemical changes in aquatic sediments (Lovley 2004).Some sequences related to the group of methanotrophic bacteria (Methylocaldum, Methylophilus) might play an important role in the methane oxidation, due to the high methane production in the Rapel Reservoir (as described by Adams et al. 2000).It is important to determine the role and dynamics of both processes (methanogenesis, sulfate reduction) because these would be responsible for rapid chemical changes in the lake, which are likely to have detrimental consequences for the biota (Vila et al. 2000).Árbol filogenético de las secuencias del ADNr 16S obtenidas desde los sedimentos de dos sitios en el embalse Rapel, determinadas por el análisis de distancia de Jukes-Cantor.Los números entre paréntesis (*) representan el número de clones con secuencia similar.Se muestran los porcentajes de "boostrap" mayores del 50 %.La escala representa el número de sustituciones por nucleótido.The fact that elevated concentrations of sulfate have entered the lake for several decades might explain the bacterial diversity we report here.For example, there is sufficient available substrate for the development of sulfate-depending organisms, which is confirmed by the high SRB-MPN.A study in four reservoirs of different trophic state (Wobus et al. 2003) showed that under mesotrophic conditions, the clone library was dominated by CFB but in a highly eutrophic reservoir the more abundant sequences correponded to Firmicutes, Cyanobacteria and Alphaproteobacteria.In comparison with Rapel reservoir, only 4 % of the sequences belong to Fimicutes, and Cyanobacteria have not been recorded for at least a decade (Vila et al. 2000), an unexpected situation, when one considers the trophic state of the lake (Van der Gucht et al. 2005).We emphasize the presence of SRB members considering the totally oxygenated water column but we need more chemical and physical studies in the sediment to understand the bacterial processes and the role of the sulfates in the sediments.

Fig. 4 :
Fig. 4: Rarefaction analysis of clone libraries from Alhué and El Manzano sites in Rapel reservoir, showing asymptote at n = 42.

Fig. 5 :
Fig. 5: Phylogenetic tree from 16S rDNA sequences determined by distance Jukes-Cantor analysis from sediments of two sites of Rapel reservoir.GenBank accession numbers are in parentheses.Number in brackets (*) represents the number of clones with similar sequence.Boostrap percentages greater than 50 % were reported.Scale bar represents substitution per nucleotide position.

TABLE 1
Variation in mean (± 1SD) physico-chemical and biological parameters recorded from the water column of Alhué basin, Rapel reservoir,(January-June, 2002)

TABLE 2
Taxonomy of 16S rDNA clone library sequencesTaxonomía de las secuencias de 16S rDNA en la genoteca