INTRODUCTION
Salar de Lagunilla is located in the Chilean altiplano, at 3800 m of altitude, is considered an athalassohaline system, because the composition of water is different from the seawater (Oren 2002). The waters of the lagoon have a variable salinity and pH neutral where the principal ionic composition corresponding to Na+, K+ and Cl-. Abiotic condition in the altiplano including, their geographical isolation and remoteness, high UV radiation and low dissolved oxygen concentration, strong variation in atmospheric pressure and, low temperature in the night and high temperature in the daily (Rodríguez-Valera 1988; Ventosa 2006). In addition, the Chilean altiplano is characterized by two very different seasons, a wet period (characterized by the input of storm-water during the months of January and February) and dry period with minimal or low rainwater in the rest of year (Karzulovic & Garda 1979). All these factors are crucial for biota inhabiting these environments, to be adapted to this extreme condition.
The study of microbial diversity and the distribution of microorganisms that inhabit in salt lakes, are crucial since these types of microorganism play a fundamental role in biogeochemical cycles such as carbon and nutrients. As example, the study of the microbial composition in saline environments might present potentially interesting in biotechnological processes as the bioremediation, UV- blockade and new antibiotics production. Also the study of their participation in biogeochemical processes as degradation and remineralization of the organic matter could provide information about he evolution of the Earth (Dundas 1998; Hoehler et al. 2001).
In the present study, we describe the composition and abundance of prokaritic comunity in the sediment samples at, five contrasting sites from Salar de Lagunilla using Fluorescent in Situ Hybridization (FISH) technique has been widely used in studies of microbial ecology (Fazi et al. 2007; Pernthaler & Amann 2005; Glöckner et al. 2000; Llobet-Brossa et al. 1998; Snair et al. 1997; Amann et al. 1991), providing an important tool to study the composition of microbial communities. In addition, we also studied the changes in temporary space induced by environmental factors, whit sediment samples taken during three years.
MATERIAL AND METHODS
SITE DESCRIPTION AND SAMPLING
The study was developed in the Salar de Lagunilla, located in the Andes mountains to the coast 155 km north of Chile (19°56`01``S 68º50`54``W). There were established five sites of permanent sampling in the lagoon including from north to south, realizing the compilation of the samples between March 2006 and February 2008. (Figure 1)
The samples were extracted from the superficial sediments of the lagoon by core (10-cm of length and 3 cm width), saved immediately to hermetic bags and fixed directly in 1% formalin-phosphate-buffered saline (PBS) (composed of 0,05 M Na2HPO4, 0,85% NaCl [pH 7 in water]). The samples were added in each section of filtere. The hibridization solution contained 20 to 35% (p/v) formamide (depending were transported at 4ºC and processed inmediately upon return to the laboratorythen, withing 5 h of sampling. From each sample was performed in surface sediment mix which took 0.5 grams, the formalin-fixed samples wer washed in PBS and centrifuged at 10.000 rpm for 5 minutes and then he supernatant is stored under following conditions in ethanol- PBS (1:1) at 20ºC. Sample was diluted and treated by mild sonication for 30 pulses each performed a second and then continues with a series of 5 min of centrifugation at 4,000 rpm. Washed sample was taken 30 μl for be filtered in 0.2 μm pore-size polycarbonate filter, 22 mm in diameter with 1 ml of sterile distilled water, yielding two filtered samples, which were subsequently analyzed by FISH.
CY3-LABELED OLIGONUCLEOTIDE PROBES
Thespecific probes for recognition of the Bacteria domain (ALF1b, BET42a, GAM42a, SBR385, CF319a) and Archaea (ARCH915), were used for the specific detection of the object cells. Probes with in order to optimize the hybridization conditions, non-fluorescent competitors probes were used under the same conditions of BET42a and GAM42a probes (Manz et al. 1992). The specific oligonucleotide sequences of CY3-labeled probes and the hybridization conditions are summarized in the Table 1.
STAINING AND CELL COUNTS OF SAMPLES
To perform FISH, Filters were quartered with a razor blade, ande the section (6 to 7) were placed face up on the glass slide. Aproximately 20 μl of preheated hibridization solution was placed on each section filter, and 50 ng specific probes on the experiment; see table 1), and a mixture (5M NaCL, 1M Tris-HCl [pH 7,2], 10% sodium dodecyl sulfate [SDS]) and incubated at 46 ° C for 2 hours.
Afterwards, the filters sections were incubated in 10 ml of prewarmed washing buffer at 48 º C for 15 min. The washing solution consisted of (20 mM Tris-HCl (pH 8), 5 mM EDTA, 0.01% SDS and variable concentration of NaCl). A formamide concentration of 35% was used 80 mM and 20% was used 0,225 mM of NaCl. The filter sections were dried on blotting paper, placed back on a glass slide, and covered whith 50 μl of DAPI solution, final concentration, 2 to 5 μg/ml, for 5 min at room temperature in the dark. Then, the filters were gently washed in 2ml of 0,2 μm-filtered destiled water and ethanol, dried on blotting paper, and mounted on glass slides.
The filter sections were inspected whith a fluorescence microscope, equipped whith a 50-W high-pressure mercury bulb and specific filter sets (DAPI [Zeiss 01]), Cy3 (Chroma HQ 41007). Each microscopic field was first viewed with CY3 filter, before switching to the DAPI filter, to avoid bleaching of CY3 during the DAPI count. For each samples and probe, more tan 500 cells were counts; for the DAPI examination, more then 1500 cells were counted per samples. All probe- specific cell counts were presented as the percentage of cell counts by DAPI. The mean abundances and standard deviations were calculated from the counts of 10 to 20 randomly chosen fields on each filter section. The account of the abundance of prokaryotic cells per gram was performed using the protocol and algorithms proposed by (Fry 1986).
STATISTICAL ANALYSIS
Of the data matrix was constructed Bray-Curtis similarity. A non-metric multidimensional scaling analysis (NMDS) was performed to study the changes in the distribution patterns and variability of the prokaryotic community. Them similarity analysis comparing the sampling sites were performed by ANOSIM (Analysis of Similarity). In this sense, ANOSIM shows an R-value statistic, which determinate the distance ranges within and between groups or samples. In general this value can be divided into 3 categories; R> 0,75: indicates that there are large differences and treatments or groups are well separated, R> 0.5: indicates (groups/samples) as separated, but overlapping, and R <0, 25: indicates little or no differences, and treatments or groups are barely separable (Clarke 1993). To avoid differences in the distribution of data by the sample size, proportions were transformed to log (x +1) (Ramette 2007). All tests were carried out using PRIMER software version 6 (Clarke & Gorley, 2005).
To study if the composition of biological community was influenced by environmental parameters, a canonical correspondence analysis (CCA) was performed using the Community Ecological Analysis software (Ecom II version 2.01), (Seaby & Henderson, 2007), was used to compared the prokaryotic abundance (FISH) matrices with the enviromental variables (physicochemical matrices). Boyh FISH data and enviromental variables were log(x+1) transformed. Forward selection was used to rank the environmental variables in importance for determining the species data (Multicollinearity test). The environmental (explanatory) variables used were nitrogen concentration measured as NO percentage of limestone measured as (CaCO3), percentage of easily oxidized carbon (CO) and electrical conductivity (EC [mS/cm]).
a Escherichia coli, numbering (Brosius . 1981).
b Percent formamide (FA) in situ hybridization buffer.
RESULTS
ESTIMATING THE ABUNDANCE OF PROKARYOTES
Table 2 shows the variation of cell counts per gram (cell g-1) and percentages of FISH counts during the months of the study. The maximum value of DAPI counts (6.98 * 10 7cells g-1) was detected during August 2006 and the minimum (2.39 *10 7cells g-1 ) was detected in April 2007; being the average DAPI counta total of 3.26*10 7cells g-1. In addtion, the percentage of FISH detectable cells to the Bacteria and Archaea domains was 9.05E +05 and 1.08e +06 (cells/g) respectively.
These sediments have a hight diversity of cell morphologies as, filamentous bacillary forms, bow-shaped rod, coccoid, flagellum cell, bacilli in chains, and among them were large consortia bacteria. Most of the cell identifica in the sediment showed a homogeneous morphology of thin short road. (Figure 2).
The groups with major variability were the domain Archaea and Alpha-proteobacteria with relative percentages of hybridized cells (FISH) were from 11.8±5.2 to 34.7±10.7 % and from 7.9±2.0 to 31±8.8% respectively, whereas that the unalterable groups were Gamma-proteobacteria and Cytofaga-flavobacteria with relative percentages of hybridized cells from 15.7±4 to 22.5±8.5% and from 14.6±10 to 28.7±2.3 % respectively. On the other hand, the most abundant group corresponded to the Archaea domain with relative abundances of 34.7 % in April, and 29% during the months of March and August (2006), and from November to February (2007). In the same way, the less abundant detected group was Beta-proteobacteria, showing a relative abundance of 5.5% during March 2006. Members of sulfate-reducing bacteria group from the subclass Delta- proteobacteria, showed an increase in their relatives abundances during the months of November, 2007 and proteobacteria group were the group with less abundances in the superficial sediments system of the salar de Lagunilla.
PHYSICOCHEMICAL RELATIONSHIP OF THE PROKARYOTIC COMMUNITY.
To describe the relations between the abundances of the prokaryotic community and the physicochemical variables (Table 3), there was made a canonical correspondence analysis (CCA). The results for the humid period (March, 2006, January, 2007 and February, 2008), makes clear 48.7 % of the accumulative variance, on the other hand the dry period (August, 2006, November, 2006, April, 2007, August, 2007 and November, 2007), the variable significant for the first two axes explain 45.2 % of the accumulative variance.
The table 4 shows the environmental variables, which contribute, significantly in the Dry period, from bigger to minor contribution in model they are, NO3 mg/K, Limy (CaCO3 g), % p/p CO (carbon easily oxidabl) and CE mS/ cm. For the humid period the significant concentrations were: % p/p CO, NO3 mg/k, limestone (CaCO3 g), CE mS/ cm and PO4 +3 mg/k expressed in the axis 1.
Biplot diagram ordination of CCA analysis shows the patterns of changes in prokaryotic community, which can be explained by environmental variables.
a Detection rate compared with DAPI. Mean and standard deviation were calculated by counting 20 fields in the section of the filter. b Unparse
Month of Study | pH | CE (mS/cm) | %CO | % Sulfato | % Caliza (CaCO3) | NO3 (mg/K) | PO4+3 (mg/K) |
March-2006 | 8.1 | 9.0 | 2.1 | 1.5 | 16.4 | NDa | NDa |
August-2006 | 7.5 | 6.3 | 1.6 | 1.1 | 13.7 | NDa | NDa |
November-2006 | 7.5 | 8.0 | 2.1 | 0.9 | 14.0 | NDa | NDa |
January -2007 | 7.5 | 8.5 | 1.7 | 0.9 | 13.8 | 24.4 | 7.1 |
April -2007 | 7.9 | 4.8 | 1.8 | 0.9 | 13.0 | 28.7 | 9.3 |
August-2006 | 8.3 | 4.1 | 1.4 | 0.8 | 15.0 | 25.5 | 7.9 |
November-2006 | 7.1 | 1.7 | 1.4 | 0.6 | 15.1 | 22.7 | 9.6 |
February -2008 | 7.8 | 8.5 | 1.4 | 0.8 | ND | 18.9 | 11.2 |
NDa, a undetermined
The environmental variables are represented by the vectors, which it is constructed by the points of the variables on the maximum inertia projected in the space. The arrows represent the vector direction. The vector length indicates the importance of the variable in the model, and the position represented by the pictures and indicate the relation with the environmental variables (vector) and his arrangement with environmental variables in the model.
The Figure 3 shows the analysis corresponding to the dry of these indicates that so correlated the environmental period, for this period the available NO3 and the Limestone variable is with the axes. Also the angles between the vectors indicate the correlation that exists between them (a small angle indicates high correlation). The samples are (CaCO3) they are the most important variables for the model, since 1 is narrowly tied to the axis.
The group Cytofaga-flavobacteria shows a good correlation of negative trend with the limy variable (CaCO3) and Gamma-protobacteria correlates negatively with the electrical conductivity (CE), explained in the axis 2, this one can represent important quantities of you work out soluble.
For the humid period (Figure 3), the most influential variables in the model are the organic carbon easily to oxidize (CO) and the available NO3. The members of the Sulfato-reductors it is correlated by the NO3 positively, Beta- proteobacterias is correlated negatively with CaCO3 , which explains to itself in the axis 2, Gamma-protobacteria with organic carbon easily oxidable (CO), Alpha-proteobacteria and the Domain Archaea show correlations positive with electrical conductivity, which explains to itself in the axis 2 and Cytofaga-flavobacteria is not correlated by any variable, in general the model shows a correlation with the samples and the variable electrical conductivity (CE).
SPATIO-TEMPORAL VARIABILITY
In term of the relative abundances of the prokaryotic community, these exhibited a great variability during the study period. In the figure 1 NMDS, shown that March and August, 2006 are different from the rest of the study, these are characterized by presenting high abundances of microorganism of Archaea domain (34.4 % and 33.5 % respectively), whereas the sum of the relative abundances between the groups of the domain Bacteria represent 47.8 % and 65,1 % respectively, with regard to the rest of the groups analyzed in comparison to other studies periods. On the other hand the figure 4 shows the analysis spreads of NMDS, in this analysis one does not show significant differences between the stations of sampling. For this study there were identified 5 sites of study (L1-L2-L3-L4-L5).
The abundance of the prokaryotic community was compared to be able to observe differences intra and interannual. The ANOSIM showed significant differences with overlapping during the period of study (R = 0.527, p <0.001), also an ANOSIM was realized to observe the spatial differences, where significant differences were not observed (R = 0.102, p <0.005).
DISCUSSION
The Highland Lakes as the Salar de lagunilla represent unique and extreme habitats clearly dominated by various forms of microbial life. However, microbiological studies have been performed only occasionally (Demergaso et al. 2004; Demergaso et al. 2008; Dorador et al. 2008a, 2008b, 2009, 2010). Salar de Lagunilla is considered as a moderately athalassohaline wetland (Table 3), comparable in its phylo-genetic groups with other high-altitude cold aquatic systems (Tibetan Lakes: Dong et al. 2006; Jiang et al. 2006. Atacama Lakes: Demergasso et al. 2004).
In this work, using FISH technique was possible to establish a first approximation of the microbial composition and the most abundant specific groups in the sediments of the Salar de Lagunilla. However, the use of this technique has its limitations, such as the weak signal and lack of specificity of the oligonucleotides. These limitations have been described and solutions have been proposed by Penthaler & Amann, (2005), thus been possible to be more objective with our results. Moreover it is important to note that the contribution achieved through these studies will help to generate other future investigations.
Bacterial assemblages were dominated by Cytofaga-Flavobacteria-Bacteriodes (CFB), (28.7% abundances of hybridized cells). This percentage is high when considering the percentage of coverage with the probe, which is about 38% (Amann et al. 2008). This group is characteristic in the oceans where it is abundant between the groups (Kirchman 2002), just as this group has been reported in sediments and in environments with high concentrations of salinity (Antón et al. 1999). These results are consistent with others studies, where has been observed in athalossahalino lake of the Atacama by PCR-denaturing gradient gel electrophoresis (DGGE) and sequencing of 16S rRNA gene fragments where the tendency for increasing contribution of (CFB) for higher salinities and altituded (Demergasso et al. 2004). In addition, in the sediment samples from three evaporitic basins located in the highlands of northern Chile, has observed the high diversity of bacteroides through of 16S rRNA gene clone libraries created with targeted Bacteroidetes-specific primers and separation of specifically amplified gene fragments by denaturing gradient gel electrophoresis (DGGE) (Dorador et al. 2008a). Similar results have also been found in high-altitude lakes from the Tibetan Plateau (Wu et al. 2006).
Cytofaga-flavobacteria have been described for their role in the degrader of complex molecules and high molecular weight cellulose and, has also been seen in the presence of abundant phytoplankton booms and association with some algae that may result in enhanced algal growth (Grossant 1999).
The increase of abundant CF in the wet period, where organic matter and salinity increases, due to the death of organisms from other trophic levels and increased salinity, and your afinity for CaCO3 apparently show that their role is similar to their role in the ocean. Therefore, this group could be the responsible for the degradation of high molecular weight compounds in altiplanic lakes.
The Proteobacteria group is characteristic as a cosmopolitan and grouped according to their habitat, creating divisions in the sea, fresh water or soil (Glockner et al. 2000; Bowman et al. 2003 & Crump et al. 2004). Alpha and Beta-proteobacteria are the group with the lowest rates of hybridization; this is consistent with studies in the ocean (Glöckner et al. 1999; Cottrell et al. 2000). These groups are not abundant in saline environments, but it includes mostly the halotolerant bacteria. However this group has been associated with salt deposits, but participating in processes of nitrification (Nold & Zwart 1998).
Gammaproteobacteria group showed variability in relative abundance in the order of 19%, which has been observed in the work of (Cottrell et al. 2000). This group is characteristic as oxides in aerobic and anaerobic environments (Ollivier et al. 1994), and it has also been shown that members of this group are facultative able to adapt to different environmental conditions. In studies recently, has been observed the present of anoxygenic phototrophic bacteria in the Atacama salar, across application of the functional pufLM genes, was determined a new phylogenetic lineage of phototrophic Gammaproteobacteria (Thiel et al. 2011). Possibly, this group have a role important as primary producers in the salar of Lagunilla
Members of the sulfate-reducing bacteria, showed variability in the abundance in the order of a 13 - 27% hybridization of cells in the period of study. This group is important for the degradation of organic matter in sediments, using sulfate as electron acceptor (Ollivier et al. 1994). In the salar de Huasco, has been observed the present of deltaproteobacteria, with a study of 16S rRNA gene sequences related to sulfate-reducing bacteria of the genera Desulfobacterium (Dorador et al.2009), this group can have a significant fraction of carbon mineralization in sediments hypersaline.
The salar de Lagunilla is considered as an athalassohaline enviroments, also this system is highly variable depending on the groundwater inputs in different zones of the lake basin and stormwater inputs, especially in the months of January and February, and bacteria have to adapt locally to the new conditions. In this studie the abundance of archaeas comunities exhibit hight contribution in this system above bacterium. Howeber has been described that dominance of Archaea over Bacteria in enviroments where NaCl concentration close to saturation, here the halobacterias are the dominant (Oren 2002). However, has been observed the predominance of methanogemic archaea and ammonia-oxidizing Archaea in studies in Salar de Huasco (Dorador et al., 2010). Probably, the Salar de Lagunilla presents dominance by methanogemic archaea, because physical-chemical characteristics of the Lagunilla´s sediments. This sediment has high carbonates saturation and anoxic sediments forming minerals, being a favorable condition for the adaptability of these organisms.
The environmental variables affected the diversity and composition of bacterial communities strongly. CCA analysis indicated that the salinity gradients, CaCO3 and carbon organic are the environmental variables more significantly close 50% of the variance, in specie composition bacterial. The enrichment of carbon and CaCO3, and the salinity gradient clearly increased the abundance of community bacterial in the sediment; however, the concentration of nitrogen (measured as nitrate) and phosphorus (measured as fosfate) do not significantly affect the abundance and bacterial diversity.
The temporal variability in bacterioplankton communities was significantly highly was highly correlated with salinity, %CO (organic carbon) and CaCO3; In contrast, spatial variability was relatively small in communities across environments altiplanicos, confirming the observations of microbial diversity and abundance based on FISH. Few studies address temporal and spatial variability using different habitats, Kirchman et al. (2010) identified a pattern where differences into winter and summer in the Arctic Ocean of the bacterioplankton communities was minimal compared with spatial variability across their sampling range. Similarly, Fortunato et al. 2012, compares the bacterioplankton community composition from river to surface ocean, by depth from surface to deep ocean, and through time seasonally over an year, by the use of 16S ARN, where they observed changes in the spatial variability and correlations with salinity and Depth.
Our study, show temporary changes that may be associated with increased salinity in dry periods and increased rainfall (wet period), which brings an eutrophication. Although, it seems that the salinity gradient is the best indicator of biogeographic distribution in these environments. We also identified significant overlap in communities of Archaea Domain in the dry period.
This study reveals the first analyses of the composition of a prokariotic community, temporary differences are demonstrated and as the environmental factors they are incidental in the structure, distribution and microbial variability of the populations in the Salar de Lagunilla. These results can serve for future investigations, so that the community analyses could be combined by the evolution of the substrates of the sediments.