Original Research ARTICLE
Enrichment experiment changes microbial interactions in an ultra-oligotrophic environment
- 1Laboratorio de Ecología Molecular y Experimental, Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Coyoacán, México
- 2Red de Estudios Moleculares Avanzados, Instituto de Ecología A.C., Xalapa, México
- 3School of Life Sciences, Arizona State University, Tempe, AZ, USA
The increase of nutrients in water bodies, in particular nitrogen (N) and phosphorus (P) due to the recent expansion of agricultural and other human activities is accelerating environmental degradation of these water bodies, elevating the risk of eutrophication and reducing biodiversity. To evaluate the ecological effects of the influx of nutrients in an oligotrophic and stoichiometrically imbalanced environment, we performed a replicated in situ mesocosm experiment. We analyzed the effects of a N- and P-enrichment on the bacterial interspecific interactions in an experiment conducted in the Cuatro Cienegas Basin (CCB) in Mexico. This is a desert ecosystem comprised of several aquatic systems with a large number of microbial endemic species. The abundance of key nutrients in this basin exhibits strong stoichiometric imbalance (high N:P ratios), suggesting that species diversity is maintained mostly by competition for resources. We focused on the biofilm formation and antibiotic resistance of 960 strains of cultivated bacteria in two habitats, water and sediment, before and after 3 weeks of fertilization. The water habitat was dominated by Pseudomonas, while Halomonas dominated the sediment. Strong antibiotic resistance was found among the isolates at time zero in the nutrient-poor bacterial communities, but resistance declined in the bacteria isolated in the nutrient-rich environments, suggesting that in the nutrient-poor original environment, negative inter-specific interactions were important, while in the nutrient-rich environments, competitive interactions are not so important. In water, a significant increase in the percentage of biofilm-forming strains was observed for all treatments involving nutrient addition.
A central goal of ecology is the understanding of the driving principles underpinning biodiversity (Gaston, 2000). Several lines of research have tried to explain the differences in diversity among local communities, focusing on approaches ranging from network interactions in food webs to random assemblages resulting from the dispersion capacity of organisms (Hutchinson, 1959; Dykhuizen, 1998; Kassen et al., 2000; Bennie et al., 2006; Marini et al., 2007). These theories, as well as other more complex explanations of biodiversity (such as chaotic interactions between competing strains; Huisman and Weissing, 1999) are partially successful, depending on the system studied and are frequently complementary. The patterns of biodiversity are determined by the combined impacts of interactions of several biotic and abiotic environmental factors, and cost-benefit strategies followed by each species can change according to the complexity of the community and its nutrient availability (Marini et al., 2007; Werner et al., 2014).
For instance, it is known that the abundances and proportions of nitrogen and phosphorus in the environment have a major influence in the composition of species at macro and microscopic scales (Makino and Cotner, 2005; Jansson et al., 2006; Marini et al., 2007; Christofoli et al., 2010; Nelson and Carlson, 2011). Classic trade-off competition theory places resource availability at the center of the processes that influence community structure (Tilman et al., 1981; Smith, 1993; Brauer et al., 2012). MacArthur and Wilson (1967) coined the terms “r-selection” and “K-selection.” This theory takes into account biotic and abiotic factors such as climate, mortality, survivorship, population size, intra and interspecific competition, relative abundance and length of life. Under this theory, r-strategists are adapted to abundant nutrients, which are rapidly exploited. On the other hand, K-strategists are adapted to a long-term survival on limited resources (Pianka, 1970; Fuchs et al., 2000; Singer et al., 2011). Experimental evidence from different microcosm experiments shows that nitrogen and phosphorus enrichments result in significant changes in community structure in terms of uniformity and species richness (Schäfer et al., 2001; Nelson and Carlson, 2011).
The changes in community structure are probably driven mainly by the effects of nutrients on the growth rates of individual strains (Smith, 1993), but other factors could also play a role such as seasonal changes (Rodríguez-Verdugo et al., 2012; Bevivino et al., 2014), as well as changes in mutualistic interactions due to fluctuations in the supply and demand of “public goods” by cause of the entrance of nutrients to the ecosystem (Morris et al., 2012; Sachs and Hollowell, 2012; Werner et al., 2014).
The release of chemical compounds into the environment, which are toxic or inhibitory to the competitors, is one of the commonly observed antagonistic interactions (Riley and Gordon, 1999; Lenski and Riley, 2002; Riley and Wertz, 2002; Kirkup and Riley, 2004; Hibbing et al., 2010; Kohanski et al., 2010; Majeed et al., 2011, 2013; Pérez-Gutiérrez et al., 2013; Aguirre-von-Wobeser et al., 2014). These secretions highly influence community structure and maintain cohesion of bacterial populations by leading to the extinction of sensitive strains in liquid media, but to coexistence in a structured media (Validov et al., 2005; Greig and Travisano, 2008; Rypien et al., 2010; Cordero et al., 2012). In the non-transitive model rock-paper-scissors (RPS), one antagonist, one sensitive and one resistant strain coexist in structured media (Czarán et al., 2002; Kerr et al., 2002; Kirkup and Riley, 2004). In a natural environment, this non-transitive relation may occur if toxic production is costly, both sensitive and resistant strains exist and costs associated with resistance are less than those of toxin production. However, the relative magnitude of each of these features is critical for coexistence (Kerr et al., 2002).
The effects of nutrient addition in microbial community structure can also depend on the extant biodiversity prior to the increase of nutrient availability, as shown by an experiment performed on the bacterial community of Owasso Lake (Minnesota, USA), which has been reported as one of the least diverse bacterial communities known (Makino and Cotner, 2005). When this lake was enriched with nitrogen, phosphorus and carbon, the community showed a response that would be expected from a single strain, rather than from a community (Jürgens and Güde, 1990), suggesting that if bacterial diversity in a given environment is low, this diversity will remain low, because of its low potential to respond at a community level to changes in its environment. In contrast, in more diverse environments, nutrient enrichment usually homogenizes the composition of the community assemblage (Donohue et al., 2009).
The Cuatro Cienegas Basin (CCB) located at the Chihuahuan desert of north central Mexico is a good system to study the effects of the addition of nutrients on bacterial communities. CCB harbors a number of highly oligotrophic aquatic ecosystems that have very low available phosphorus levels, as well as a stoichiometric disequilibrium with nitrogen and thus are strongly limited by phosphorus (Elser et al., 2005, 2006). We have suggested that the high microbial species diversity in CCB is strongly shaped by the stress of low nutrient supplies and the impacts of interspecific competition (Souza et al., 2008). Previous studies indicate that CCB bacterial communities display strategies to cope with this lack of nutrients, including a high representation of genes involved in phosphorus assimilation (e.g., pho and pst in Mesquites river; Breitbart et al., 2009) and the presence of a large number of genes related to the production and resistance to antibiotics (in Pozas Rojas; Bonilla-Rosso et al., 2012; Peimbert et al., 2012).
Pseudomonas and other cultivable proteobacteria are abundant in different aquatic systems in CCB, including the Churince system (Escalante et al., 2009) and in the Los Hundidos region (Bonilla-Rosso et al., 2012; Peimbert et al., 2012). In a different study in the flow system Churince, it was found that the Pseudomonas genus exhibited a seasonal variation across summer and winter (Rodríguez-Verdugo et al., 2012).
In this paper, we analyze the seasonal response of bacteria cultivable in Pseudomonas Isolation Agar (PIA; Difco1, Detroit, MI) in response to the nutrient amendment. Our study seeks to understand the effects of nitrogen and phosphorus fertilization on the interaction potential of the cultivable gamma-proteobacteria community in two different habitats, water and sediment, by determining how nutrient enrichment affected various features related to microbial species interactions, such as the tendency to form biofilms and resistance to antibiotic compounds. We are clearly aware of the limitations of culture media in order to capture the diversity of a site, as we have described elsewhere (see Souza et al., 2006, 2012), nevertheless, only in culture we can study the physiology and sociology of particular strains, which is the particular aim of this study. Our findings shed light on how the supplies of key nutrients, such as N and P modulate community structure, as well as the nature and intensity of interspecific interactions.
Materials and Methods
The nutrient enrichment experiment was conducted in situ in a small shallow evaporitic pond, Lagunita (26.84810° N, 102.14160° W), lateral to the main Churince flow system. The Churince flow system, located at the western region of CCB is dominated by gypsum sediments and has a strong longitudinal gradient of salinity, temperature, pH and dissolved oxygen (Cerritos et al., 2011). Lagunita is characterized by low phosphorus concentrations (PO4 as low as 0.1 μM and often below detection) (Elser et al., 2005), but relatively high concentrations of inorganic N and thus high N:P ratios (>200:1 for total nutrients; Lee et al., in press). Lagunita is subjected to strong evaporation, the greatest water depth at the beginning of the experiment was 32 cm and decreased to as low as 12 cm by the end of the experiment.
Experimental Design of Mesocosm
The mesocosm experiment was conducted from May to June 2011. Each mesocosm consisted of a round clear plastic tube with a diameter of 40 cm. The tube had a depth of 20 cm into the sediment and approximately 20 cm above the water surface. The mesocosms were arranged in a randomized complete block with a total of 5 blocks, separated about 2-4 m from each other. Each block consisted of four treatments. A non enriched control treatment, a phosphorus enrichment treatment (P), amended with KH2PO4 and maintained at 2 day intervals at a final concentration of 1 μM; a nitrogen and phosphorus enrichment (NP), amended as above with KH2PO4 to 1 μM but also with NH4NO3, to achieve N:P of 16:1; and a nitrogen and phosphorus enrichment with extra nitrogen (NNP), amended as above with 1 μM KH2PO4 but also with NH4NO3 at an N:P ratio = 75:1.
Methodology for Chemical Analyses
For the chemical analyses, water was collected in acid-washed 2 L cubitainers. Water samples were filtered through pre-combusted (24 h at 450°C) GF/F filters (Whatman, Piscataway, NJ) for seston elementary analysis and stored at −20°C. To measure total dissolved nutrients, water samples were filtered through 0.2 μm polyethersulfone membrane filters. Samples for dissolved organic carbon (DOC), and total dissolved nitrogen (TDN) were acidified with 12 N HCL to pH < 2 and stored in the dark at room temperature, while the remaining filtrate was frozen for total dissolved phosphorus (TDP) and soluble reactive phosphorus (SRP) analyses.
GF/F with seston were thawed, dried at 60°C and packed into tin disks (Elemental Microanalysis, U.K.) for N analyses with a Perkin Elmer™ PE 2400 CHN Analyzer at the Arizona State University Goldwater Environmental Laboratory (ASU GEL). Another set of dried GF/F filters from the same water samples was used in order to estimate seston P content. These filters were digested in persulfate followed by a colorimetric analysis to determine PO3−4 (APHA, 2005). TDP concentrations were determined using the colorimetric assay after persulfate digestion as previously described; SRP was measured without the persulfate digestion. DOC and TDN were analyzed using the Shimadzu TOC-VC/TN analyzer at the ASU GEL. Total Phosphorus (TP) concentration was calculated as the sum of the seston and total dissolved pools.
Sample Collection and Processing
Water and sediment samples were taken for each one of the four treatments for each experimental block (20 samples total). Initial samples were obtained from surface water and top of sediment prior to the application of the treatments (T0; 14 May 2011) and after 21 days of enrichment (4 June 2011) using sterile BD Falcon vials (BD Biosciences, San Jose, CA). Cultivable bacteria were obtained by plating 100 μl of each water sample or 100 μl of a 1:10 sediment dilution, prepared with 0.9% NaCl solution. Strains were isolated from water and sediment using PIA medium. Strains were incubated on agar plates at room temperature for 2 days at the field, and then were kept at 4°C until isolation in the laboratory. Individual colonies were transferred to new PIA plates and then were grown at 30°C. A total of 960 isolates were obtained.
DNA Extraction and PCR Amplification
Phylogenetic identification of the strains was performed using the 16S rRNA gene. For the isolates, DNA extraction was very complicated, and thus several methods were tested (Chen and Kuo, 1993; Aljanabi and Martinez, 1997; Reischl et al., 2000, and DNeasy Blood and Tissue kits, Qiagen, Hilden, Germany). Ultimately, the DNeasy Blood and Tissue kit was the method used, and from the 960 isolates, good quality DNA was obtained for 152 strains. 16S rRNA genes were amplified using universal primers 27F (5′-AGA GTT TGA TCC TGG CTC AG-3′) and 1492R (5′-GGT TAC CTT GTT ACG ACT T-3′) (Lane, 1991) and high fidelity Phusion hot start DNA polymerase (Finnzymes, Espoo, Finland). All reactions were carried out in a Techne TC-3000 thermal cycler (Barloworld Scientific, Staffordshire, UK) with the following program: 94°C for 5 min, followed by 30 cycles consisting of 94°C for 1 min, 50°C for 30 s, 72°C for 1 min and 72°C for 5 min. Polymerase chain reaction (PCR) amplification products were electrophoresed on 1% agarose gels. Sanger sequencing was performed at the University of Washington High Throughput Genomics Center. The sequences have been uploaded to GenBank with accession numbers (KF317734-KF317770, KM352505-KM352636).
The 16S rRNA sequences were aligned with CLUSTALW (Larkin et al., 2007) and MUSCLE (Edgar, 2004), and the alignments were manually revised. For the reconstruction of the phylogenetic tree, a maximum likelihood analysis was done with PhyML version 3.0 (Guindon et al., 2010) with the TrN+G model. The substitution model was calculated with jModelTest 2.1.3 (Darriba et al., 2012). The degree of statistical support for the branches was determined with 1000 bootstrap replicates. Genera level identification of the strains was made using the classifier tool (Wang et al., 2007) from the Ribosomal Database Project (RDP) Release 10, update 30 (Cole et al., 2009; Table 1). We performed a local BLAST search (Altschul et al., 1990) to find the nearest neighbors using the 16S ribosomal RNA. These analyses were performed with 700 bp from the 5′ end of all the sequences.
Antibiotic Resistance Assays
The 960 isolates were grown on LB plates supplemented with the following antibiotic concentrations: Carbenicillin 500 μg/ml, Kanamycin, 200 μg/ml, Tetracycline, 150 μg/ml, Streptomycin, 200 μg/ml and Gentamicin 150 μg/ml. Isolates were labeled as resistant to an antibiotic if we observed bacterial growth; otherwise, they were considered as nonresistant.
For statistical inference, the antibiotics data were organized as contingency tables, and Barnard's exact test (Barnard, 1945) was performed. Independent tests were conducted for each treatment against the initial time, for water and sediment data. Since several treatments were compared to the initial condition, a Bonferroni correction was applied for multiple tests.
Biofilm Formation Assay
A microtiter dish assay was performed for the same 960 isolates as described by O'Toole (2011) with the following modifications. The overnight culture was grown in LB medium. The microtiter plates were incubated for 18 h at 30°C. Each strain was analyzed by triplicate and we used the strains P. aeruginosa PAO1 and E. coli MC4100 as positive and negative controls, respectively. After incubation, plates were gently washed with water, and subsequently 125 μl of a 0.1% solution of crystal violet in water was added to each well. The plates were incubated at room temperature for 15 min and then washed with water. 125 μl of 30% acetic acid in water was added to each well to solubilize the crystal violet and it was incubated for 15 min. The volume was transferred to a new flat-bottomed microtiter dish. Absorbance was read with a Synergy HT plate reader (BioTek, Winooski, VT) at 550 nm.
To determine the statistical significance of differences in biofilm formation between the different treatments and the initial condition, we organized the data as contingency tables including the number of positive and negative strains for biofilm formation for the initial time and each treatment, and applied Barnard's exact tests. To correct multiple testing, a Bonferroni correction was applied.
To characterize the changes in diversity of cultivated Proteobacteria related to pseudomonads associated with nutrient enrichment in a water system, we performed an in situ mesocosm experiment with three experimental manipulations, adding phosphate (P), phosphate and nitrogen (NP), and phosphate and excess nitrogen (NNP), as well as an un-enriched control. The experiment is part of a bigger project that is described in detail elsewhere (Lee et al., in press). Cultures were obtained from surface water and sediment in two sampling events, prior to the experimental manipulations, and after 21 days, and a total of 960 isolates were analyzed for interaction phenotypes such as biofilm formation, and antibiotic resistance. A subsample of the isolates was characterized by 16S rRNA sequence.
Nutrient conditions during the experiment were characterized by low concentrations of P but relatively high concentrations of N that increased in NP and NNP treatments. Despite nutrient enrichment, N:P ratios remained quite stoichiometrically imbalanced. Except for the NNP treatment, total dissolved phosphorus and soluble reactive phosphorus were significantly different respect to control. Total phosphorus, total dissolved phosphorus and the N:P ratio were significantly different in all enriched treatments. Total dissolved nitrogen was significantly different only in the nitrogen amended treatments (Table 2). Details of dynamics and fate of N and P in the mesocosms are given elsewhere (Lee et al., in press).
Phylogenetic Diversity and Responses to Treatments
The phylogenetic relationships of partial 16S rRNA gene sequences (700 bases) of the 152 strains were determined by queries against the Ribosomal Database Project (Cole et al., 2009) using SeqMatch and Classifier (Wang et al., 2007). The phylogenetic results showed that the strain collection was dominated by Proteobacteria (145 isolates/95.4%), with members from the gamma- (137/90.1%) and alpha- (8/5.3%) subdivisions (Table 1), consistent with the use of the PIA medium for isolation. We also found several Actinobacteria (7/4.6%). As described by King et al. (1954), this medium includes Irgasan®, a broad spectrum antimicrobial substance that is not active against Pseudomonas spp. The medium also enhances the formation of pigments by Pseudomonas by adding magnesium chloride and potassium sulfate. Pseudomonas species were found in all un-enriched water cultures as well as in almost all enriched samples, with the exception of the NP treatment. Pseudomonas was the most abundant genus for the water isolates and overall for all samples (35.5% of the isolates). Halomonas was the most abundant genus among the sediment strains, representing 30.9% of the isolates.
Using the aligned sequences, a maximum likelihood tree was constructed (Figure 1). For any given treatment, groups of closely related, almost identical strains dominated the samples, even though they were obtained from different replicate mesocosms. Previous to the experimental manipulations, only gamma-proteobacteria were retrieved from the water samples, which were dominated by Pseudomonas. Other genera were present, including Aeromonas, Shewanella, Citrobacter, Rheinheimera, and Halomonas. With the addition of phosphorus (P treatment), many closely related Pseudomonas were retrieved. In the NP treatment, all the obtained isolates belonged to the genus Aeromonas. For the high N:P treatment (NNP), a group of Pseudomonas became abundant and two Aeromonas isolates were obtained.
Figure 1. Phylogenetic relationships among 16S ribosomal RNA genes. Sequences from the 152 isolates (highlighted in colors) and 21 reference strains. The coded name of the strains corresponds to the treatments and the environment from which they were isolated. Water samples: #[A-E]#, mesocosm T0, XC#, control treatment, XP#, phosphorus enrichment, XNP#, N:P 16:1 treatment, XNNP#, and N:P 75:1 treatment. Sediment samples: JPXs_#, mesocosm T0, XCs_#, control treatment, XPs_#, phosphorus treatment, XNPs_#, N:P 16:1 treatment, XNNPs_#, and N:P 75:1 treatment, where # is a number and X a letter from A to E, each one represents a replicate.
Phylogenetic differences were observed between the abundant strains obtained from the different habitats (water vs. sediment) but also between the strains corresponding to the different treatments. At the beginning of the experiment, the sediment samples were dominated by another gamma-proteobacterium, Stenotrophomonas, as well as an alpha-proteobacterium, Rhizobium. Some Actinobacteria were also isolated from these initial samples, belonging to the genus Brachybacterium. With P enrichment, Rhizobium was also found, as well as an Actinobacterium, from the genus, Kocuria. However, a larger change in the sampled community was found in all treatments including nitrogen along with P (NP and NNP), where Halomonas (gamma-proteobacteria) was largely dominant.
We analyzed the prevalence of antibiotic resistances among the 960 isolates. The antibiotics analyzed were Carbenicillin, Kanamycin, Tetracycline, Streptomycin and Gentamicin (Supplementary Table 1; Supplementary Figures 1, 2). Carbenicillin was the most common form of resistance among the isolates. The Carbenicillin resistance was found in 55.2% of the water isolates at time zero and in 68.3% of the total water isolates after the experiment, and in 28.6 and 90.5% of the total sediment isolates, at time zero and after the experiment, respectively. Tetracycline was the least common form of resistance, being represented in 10.4 and 0.8% of the water isolates and 7.1 and 1.0% of the sediment isolates before and after the experiment, respectively. Kanamycin resistance was not observed in sediment samples after the experiment. Overall, there was a statistically significant decrease in the number of resistant strains after the experiment (Bernard's exact test; p < 0.05; Figure 2; Table 3), in all three fertilization treatments and in the control, with the marked exception of Carbenicillin, for which resistance increased significantly in NNP and P treatments in water and in all treatments in sediment.
Figure 2. Prevalence of antibiotic resistance among 960 isolates. The isolates analyzed for each treatment were as follows. T0—163, Control—87, NNP—186, NP—176, P—201, T0S—42, CS—11, NNPS—38, NPS—31, and PS—25.
To further characterize the ecological traits of the studied strains, biofilm formation was analyzed for 923 isolates. In the water samples, the enriched isolates exhibited an increase in the proportion of the ability to form biofilm (Figure 3; Supplementary Figure 1). In a structured environment, interactions can be enhanced by biofilm formation. As expected, sediment strains generally had a greater tendency to produce biofilm that those isolated from water. However, in sediment, differences between before and after the treatments were not significant (Barnard's exact test; p > 0.05; Table 4), except for the NNP treatment, for which no biofilm-producing strains were detected (Figure 3; Supplementary Figure 2). In the case of water, a significant increase in the percentage of biofilm-producing strains was observed for all treatments involving nutrient additions, namely P, NP, and NNP (Barnard's exact test; p < 0.05; Table 4).
Figure 3. Relative frequency of biofilm forming strains among 923 isolates. The isolates analyzed for each treatment were as follows. T0—157, Control—84, NNP—174, NP—167, P—196, T0S—42, CS—11, NNPS—38, NPS—29, and PS—25.
In this study, the effect of experimental nutrient enrichment in a shallow, nutrient-deficient pond was analyzed to assess changes in the composition of the cultivable microbial community as a function of the N:P ratio of enrichment, as well as the modifications in the interaction network among cultivable isolates. The selective medium PIA was found to be highly selective to gamma-proteobacteria, a class previously reported as abundant at Churince as well as other sites at CCB (Souza et al., 2006; Escalante et al., 2009; Bonilla-Rosso et al., 2012) and it has not been analyzed in detail. Unfortunately, a great amount of the isolates could not be identified, as it seems that Cuatro Ciénegas bacteria produce compounds that inhibit the PCR reaction, compounds that remain despite the use of several DNA extraction methods. To analyze the community response to the nutrient enrichment without the cultivation bias, 16S libraries from water and sediment samples, from May to June 2011, were also sequenced. The analysis of the community showed, in agreement with this study, a shift in the community composition with the disappearance of several bacterial genera after the increased nutrient availability (Elser et al., 2014, in preparation).
Our hypothesis that Lagunita is P-limited was supported, as it was found that all added P that remained in the water column was immobilized into seston (Table 2). This is reflected by the fact that TDP and SRP were not significantly different respect to the control, but TP and phosphorus in seston were. This difference was also observed in the N:P ratio. With respect to total dissolved nitrogen, it was found that both nitrogen-amended treatments were significantly different but the phosphorus-amended treatment was not. This supports that our observations were not due to an isolation bias.
Based on the 16S rRNA sequences from the cultivated strains, we found that two different groups of gamma-proteobacteria responded to increases of nutrient availability in water and sediment environments, Pseudomonas and Halomonas, respectively. The strains characterized differed among the different treatments in the physiological characteristics analyzed, biofilm formation and antibiotic resistance. In contrast, genera such as Shewanella, Citrobacter, and Rheinheimera in water, as well as Stenotrophomonas and Brachybacterium in sediment were not isolated from the fertilized mesocosms, although they were present at time zero. Moreover, some rare genera, such as Kocuria, which was not documented at time zero, was found after the enrichment. However, we note that this rare genus of Actinobacteria has been obtained in past work at CCB (Cerritos et al., 2011), so it is not foreign to this environment. Several genera isolated in this study, had not been previously reported to grow in PIA medium (McCaig et al., 2001; Rajkowski and Rice, 2001; Falcone-Dias et al., 2012; Weiser et al., 2014).
All these changes in bacterial groups, which were isolated during different times in the experiment, suggest that an important community change at the bacterial level took place after the enrichment, however it should be taken carefully given the sampling method. It can be speculated that the original community in Lagunita was dominated by K-strategists, which rely on long-term survival on limited resources (Pianka, 1970; Fuchs et al., 2000; Singer et al., 2011). After the enrichment, the microbial community was dominated by faster-growing gamma-proteobacteria, which can be considered r-strategists, rapidly exploiting nutrient patches and then dying or becoming dormant after substrate exhaustion. Indeed, Pseudomonas has been previously characterized in general as a r-strategist, as pseudomonads rapidly colonize and grow on nutrient-rich environments (Juteau et al., 1999; Margesin et al., 2003), due to its metabolic versatility (Clarke, 1982; Hallsworth et al., 2003; Domínguez-Cuevas et al., 2006). This shift is important as it has been reported that K-strategists are expected to allocate more energy and interact in a broader way with their environment, for example developing strategies to cope with their environment, than to grow (Fontaine et al., 2003).
Although the experiment showed a reduction in diversity of the isolated strains and interactions as a whole, details of these responses of strains to the individual treatments were largely idiosyncratic and seemed independent of their phylum. For example, the group that includes strain AP29, affiliated with Pseudomonas, responded to the P treatment, while another group of Pseudomonas, which includes strain ENNP10, responded to the NNP treatment. However, another Pseudomonas group was abundant at the beginning of the experiment and remained abundant in the control, as well as in the P treatment (group including strains AP9, AC10, and 1C11) (Figure 1). The sudden availability of P in P-deficient environment could favor the dominance of the Pseudomonas genus, which has been reported to solubilize phosphate (Park et al., 2009; Parani and Saha, 2012). In sediment samples it was observed in both N+P treatments that Halomonas exhibited the greatest response. This observation is consistent with the fact that the Halomonas genus has been identified as capable of denitrification and may have taken advantage of the added NO3 in the NP and NNP treatments (Mormille et al., 1999; Guo et al., 2013).
The production of chemical compounds as bacteriocins and/or antibiotics is a common mechanism of antagonism among microorganisms (Riley and Gordon, 1999; Lenski and Riley, 2002; Riley and Wertz, 2002; Kirkup and Riley, 2004; Hibbing et al., 2010; Kohanski et al., 2010; Majeed et al., 2011, 2013; Pérez-Gutiérrez et al., 2013; Aguirre-von-Wobeser et al., 2014). In this study, we observed a general decrease in the antibiotic resistance in both the water and sediment environment, for all antibiotics except Carbenicillin, after the experiment (Table 3). This decrease likely reflects modifications in microbial survival strategies under different conditions, including the control were wind movement of both water and sediment, was restricted by the mesocosm tubes. Given the antagonistic network previously documented among Pseudomonas bacterial isolates from CCB (Aguirre-von-Wobeser et al., 2014), this antagonism could be due to the competition for resources. Following “microbial market logic” (Werner et al., 2014), without the acute nutrient limitation, the cost of producing antibiotics to repel competitors for a scarce resource is no longer beneficial in the case of increased P and the ideal N:P ratio, while in the NNP treatment, the community is so perturbed by the further limitation of P in relation to N that antagonism or cooperation through biofilm formation is no longer an economic option. On the other hand, in a rock—paper—scissors (RPS) model behavior scenario, that can be applied only for structured environment (Kirkup and Riley, 2004; Nahum et al., 2011) such as sediment or biofilm, the strains that produce toxins (C) kill sensitive strains (S), which outcompete resistant strains (R), which in turn outcompete C (Czarán et al., 2002; Kirkup and Riley, 2004). In this RPS game, the resistant and producer strains spend resources to keep the resistance, which in an enriched environment may be no longer needed, thus those strains are outcompeted by the sensitive strains (Kerr et al., 2002). Interestingly, this shift from a collection of isolates with a high prevalence of resistance to several antibiotic tested, to a collection dominated by sensitive strains was observed in both habitats: water and sediment (Figure 2), suggesting that the fact of predicting the neighborhood, as required for the rock-paper-scissor model, is not a requisite for antagonism, at least in the analyzed system.
Market logic suggests that local environments determine trade connections (Werner et al., 2014), and biofilm formation is a way to ensure a local environment both in water and sediment. The mechanisms of bacterial biofilm formation are processes by which single cells coordinate and implement the formation of complex surface-attached communities (Davey and O'Toole, 2000). Bacteria that reside within the biofilm are to some extent isolated from environmental stresses, such as desiccation or nutrient limitation (Danhorn et al., 2004). Complex interactions, in particular mutualistic behaviors, are expected in a biofilm, since the secretion of the matrix that forms the biofilm is a form of public good that will increase the survival of the coexisting partners and will avoid the presence of cheaters (Werner et al., 2014).
Biofilm-forming strains were found in our experiment in both habitats at time zero and in most of the enrichment treatments, suggesting that this cooperation strategy is more ingrained in the whole CCB microbial community than antagonistic interactions promoted by antibiotics, where rare species cannot afford to pay its costs. Biofilm formation was not present only in the sediments of the NNP treatment where P was further limited in relation to N. This treatment was dominated by strain of Halomonas. The NP treatment was also dominated by closely related strains of Halomonas that were biofilm-forming (Supplementary Figure 2). This differential response suggests that a wide range of strategies to cope with environmental limitations is present in CCB microbes, even within a single genus.
As the Black Queen Hypothesis (BQH) suggests, certain biological functions are not only expensive, but are also broadly distributed in the community, since they are public goods. Hence, the majority of the community can afford to lose those functions if they are at least a proportion of helpers that produce such public good (Morris et al., 2012). As a consequence, for a microbial market, we also need “price” differences and supply, as well as “demand” variation. In this study, we observed a shift in the composition of the bacterial community, from a diverse community prior the enrichment to a community dominated by few genera such as Pseudomonas in water or Halomonas in sediment after the addition of N+P. As previously stated, the change was not only related to the community composition, but also to its physiological characteristics as antibiotic resistances and the ability of biofilm formation. These changes could be explained with this microbial market theory as the supplies changed with the experiment the balance between cost and benefit.
According to the market theory, the benefit of trade depends not only on the interacting partners but also on the available supply of commodities from other sources and it is expected that biotic or abiotic conditions influence the demand for a particular service (Werner et al., 2014). Based on an analysis of Illumina tags of 16S rRNA (Elser et al., 2014, in preparation), our study site exhibits a large bacterial diversity dominated by alpha-proteobacteria and bacteroidetes. In the original community, we expected to find a strong competition among the members of this community. Under the original condition of extreme low P and N availabilities, it is predicted that specialization will be favored. All this considering, it is not surprising that when this particularly diverse and fragile network of interactions was perturbed by the mesocosm conditions and the nutrient input, the community structure and its biological market equilibrium changed, reducing its overall diversity, not only in the few cultured genera that we could follow, but in the overall community (Elser et al., 2014), suggesting that the resilience of this extremely oligotrophic oasis depends precisely on the permanence of such unbalanced stoichiometry.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This paper constitutes a partial fulfillment of the Graduate Program in Biological Sciences of the National Autonomous University of México (UNAM). GYP-S acknowledges the scholarship and financial support provided by the National Council of Science and Technology (CONACyT; Grant no. 412697/269818), and UNAM. This study was conducted with financial support from PAPIIT project IN203712-3, FCS-Alianza WWF and CONACYT. We thank Abigail González-Valdéz, Victoria Grosso-Becerra, Mirna Vázquez-Rosas-Landa, Erika Aguirre-Planter and Laura Espinosa-Asuar for technical assistance. JJE acknowledges support from the U. S. National Science Foundation (DEB-0950175) and National Aeronautics and Space Administration Astrobiology (NAI5-0018). All samples were collected under the permit SGPA/DGVS/00614/13 FAUT – 0230 granted to VS by the “Secretaría de Medio Ambiente y Recursos Naturales (SEMARNAT).”
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2015.00246/abstract
1. ^DifcoTM and BBLTM Manual, 2nd Edition. Pseudomonas Isolation Agar.
Aguirre-von-Wobeser, E., Soberón-Chávez, G., Eguiarte, L. E., Ponce-Soto, G. Y., Vázquez-Rosas-Landa, M., and Souza, V. (2014). Two-role model of an interaction network of free-living γ-proteobacteria from an oligotrophic environment. Environ. Microbiol. 16, 1366–1377. doi: 10.1111/1462-2920.12305
Bennie, J., Hill, M., Baxter, R., and Huntley, B. (2006). Influene of slope aspect on long-term vegetation change in British chalk grasslands. J. Ecol. 94, 355–368. doi: 10.1111/j.1365-2745.2006.01104.x
Bevivino, A., Paganin, P., Bacci, G., Florio, A., Pellicer, M. S., Papaleo, M. C., et al. (2014). Soil Bacterial community response to differences in agricultural management along with seasonal changes in a mediterranean region. PLoS ONE 9:e105515. doi: 10.1371/journal.pone.0105515
Bonilla-Rosso, G., Peimbert, M., Alcaraz, L. D., Hernández, I., Eguiarte, L. E., Olmedo-Alvarez, G., et al. (2012). Comparative metagenomics of two microbial mats at Cuatro Cienegas Basin II: community structure and composition in oligotrophic environments. Astrobiology 12, 659–673. doi: 10.1089/ast.2011.0724
Brauer, V. S., Stomp, M., and Huisman, J. (2012). The nutrient-load hypothesis: patterns of resource limitation and community structure driven by competition for nutrients and light. Am. Nat. 179, 721–740. doi: 10.1086/665650
Breitbart, M., Hoare, A., and Nitti, A. (2009). Metagenomic and stable isotopic analyses of modern freshwater microbialites in Cuatro Cienegas, Mexico. Environ. Microbiol. 11, 16–34. doi: 10.1111/j.1462-2920.2008.01725.x
Cerritos, R., Eguiarte, L. E., Avitia, M., Siefert, J., Travisano, M., Rodríguez-Verdugo, A., et al. (2011). Diversity of culturable thermo-resistant aquatic bacteria along and environmental gradient in Cuatro Cienegas, Coahuila, México. Antonie Van Leeuwenhoek 99, 303–318. doi: 10.1007/s10482-010-9490-9
Christofoli, S., Monty, A., and Mahy, G. (2010). Historical landscape structure affects plant species richness in wet heathlands with complex landscape dynamics. Landsc. Urban Plann. 98, 92–98. doi: 10.1016/j.landurbplan.2010.07.014
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The Ribosomal Database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37 D141–D145. doi: 10.1093/nar/gkn879
Cordero, O. X., Wildschutte, H., Kirkup, B., Proehl, S., Ngo, L., Hussain, F., et al. (2012). Ecological populations of bacteria act as socially cohesive units of antibiotic production and resistance. Science 337, 1228–1231. doi: 10.1126/science.1219385
Danhorn, T., Hentzer, M., Givskov, M., Parsek, M. R., and Fuqua, C. (2004). Phosphorus limitation enhances biofilm formation of the plant pathogen Agrobacterium tumefaciens through the PhoR-PhoB regulatory system. J. Bacteriol. 186, 4492–4501. doi: 10.1128/JB.186.14.4492-4501.2004
Domínguez-Cuevas, P., González-Pastor, J. E., Marqués, S., Ramos, J. L., and de Lorenzo, V. (2006). Transcriptional tradeoff between metabolic and stress-response programs in Pseudomonas putida KT2440 cells exposed to toluene. J. Biol. Chem. 281, 11981–11991. doi: 10.1074/jbc.M509848200
Elser, J. J., Elser, T. J., Carpenter, S. R., and Brock, W. A. (2014). Regime shift in fertilizer commodities indicates more turbulence ahead for food security. PLoS ONE 9:e93998. doi: 10.1371/journal.pone.0093998
Elser, J. J., Schampel, J. H., Garcia-Pichel, F., Wade, B. D., Souza, V., Eguiarte, L., et al. (2005). Effects of phosphorus enrichment and grazing snails on modern stromatolitic communities. Freshw. Biol. 50, 1808–1825. doi: 10.1111/j.1365-2427.2005.01451.x
Elser, J. J., Watts, J., Schampel, J. H., and Farmer, J. (2006). Early food webs on a trophic knife-edge? Experimental data from modern microbialite-based ecosystem. Ecol. Lett. 9, 295–303. doi: 10.1111/j.1461-0248.2005.00873.x
Escalante, A., Caballero-Mellado, J., Martínez-Aguilar, L., Rodríguez-Verdugo, A., González-González, A., Toribio-Jiménez, J., et al. (2009). Pseudomonas cuatrocienegasensis sp. nov., isolated from an evaporating lagoon in the Cuatro Cienegas valley in Coahuila, Mexico. Int. J. Syst. Evol. Microbiol. 59, 1416–1420. doi: 10.1099/ijs.0.006189-0
Falcone-Dias, M. F., Vaz-Moreira, I., and Manaia, C. M. (2012). Bottled mineral water as a potential source of antibiotic resistan bacteria. Water Res. 46, 3612–3622. doi: 10.1016/j.watres.2012.04.007
Fuchs, B. M., Zubkov, M. V., Sahm, K., Burkill, P. H., and Amann, R. (2000). Changes in community composition during dilution cultures of marine bacterioplankton as assessed by flow cytometric and molecular biological techniques. Environ. Microbiol. 2, 190–201. doi: 10.1046/j.1462-2920.2000.00092.x
Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0). Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010
Jansson, M., Bergstöm, A-K., Lymer, D., Vrede, K., and Karlsson, J. (2006). Bacterioplankton growth and nutrient use efficiencies under variable organic carbon and inorganic phosphorus ratios. Microb. Ecol. 52, 358–364. doi: 10.1007/s00248-006-9013-4
Juteau, P., Larocque, R., Rho, D., and LeDuy, A. (1999). Analysis of the relative abundance of different types of bacteria capable of toluene degradation in a compost biofilter. Appl. Microbiol. Biotechnol. 52, 863–868. doi: 10.1007/s002530051604
Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi: 10.1093/bioinformatics/btm404
Lee, Z. M.-P., Steger, L., Corman, J. R., Neveu, M., Poret-Peterson, A., Souza, V., et al. (in press). Stoichiometrically imbalanced shallow pond in the Cuatro Ciénegas Basin, Mexico, is highly sensitive to nutrient enrichment. PLoS ONE.
Makino, W., and Cotner, J. B. (2005). Elemental stoichiometry of a heterotrophic bacterial community in a freshwater lake: implications for growth- and resource-dependent variations. Aquat. Microbiol. Ecol. 34, 33–41. doi: 10.3354/ame034033
Margesin, R., Labbé, D., Schinner, F., Greer, C. W., and Whyte, L. G. (2003). Characterizaton of hydrocarbon-degrading microbial populations in contaminated and pristine alpine soils. Appl. Environ. Microbiol. 69, 3085–3092. doi: 10.1128/AEM.69.6.3085-3092.2003
Marini, L., Scotton, M., Klimek, S., Isselsteinc, J., and Pecile, A. (2007). Effects of local factors on plant species richness and composition of Alpine meadows. Agric. Ecosyst. Environ. 119, 281–288. doi: 10.1016/j.agee.2006.07.015
McCaig, A. E., Grayston, S. J., Prosser, J. I., and Glover, A. (2001). Impact of cultivation on characterisation of species composition of soil bacterial communities. FEMS Microbiol. Ecol. 35, 37–48. doi: 10.1111/j.1574-6941.2001.tb00786.x
Mormille, M. R., Romine, M. F., Garcia, M. T., Ventosa, A., Bailey, T. J., and Peyton, B. M. (1999). Halomonas campisalis sp. nov. a denitrifying moderately haloalkaliphilic bacterium. Syst. Appl. Microbiol. 22, 551–558. doi: 10.1016/S0723-2020(99)80008-3
Nelson, C. E., and Carlson, C. A. (2011). Differential response of high-elevation planktonic bacterial community structure and metabolism to experimental nutrient enrichment. PLoS ONE 6:e18320. doi: 10.1371/journal.pone.0018320
Park, K. H., Lee, C. Y., and Son, H.-J. (2009). Mechanisms of insoluble phosphate solubilization by Pseudomonas fluorescens RAF15 isolated from ginseng rhizosphere and its plant growth-promoting activities. Lett. Appl. Microbiol. 49, 222–228. doi: 10.1111/j.1472-765X.2009.02642.x
Peimbert, M., Alcaraz, L. D., Bonilla-Rosso, G., Olmedo-Alvarez, G., García-Oliva, F., Segovia, L., et al. (2012). Comparative metagenomics of two microbial mats at Cuatro Cienegas Basin I: ancient lessons on how to cope with an environment under severe nutrient stress. Astrobiology 12, 648–658. doi: 10.1089/ast.2011.0694
Pérez-Gutiérrez, R. A., López-Ramírez, V., Islas, Á., Alcaraz, L. D., Hernández-González, I., Olivera, B. C., et al. (2013). Antagonism influences assembly of a Bacillus guild in a local community and is depicted as a food-chain network. ISME J. 7, 487–497. doi: 10.1038/ismej.2012.119
Reischl, U., Linde, H. J., Metz, M., Leppmeier, B., and Lehn, N. (2000). Rapid identification of methicillin-resistant Staphylococcus aureus and simultaneous species confirmation using real-time fluorescence PCR. J. Clin. Microbiol. 38, 2429–2433. doi: 10.1128/JCM.40.7.2392-2397.2002
Rodríguez-Verdugo, A., Souza, V., Eguiarte, L. E., and Escalante, A. E. (2012). Diversity across seasons of culturable Pseudomonas from a desiccation lagoon in Cuatro Cienegas, México. Int. J. Microbiol. 2012:201389. doi: 10.1155/2012/201389
Schäfer, H., Bernard, L., Courties, C., Lebaron, P., Servais, P., Pukall, R., et al. (2001). Microbial community dynamics in Mediterranean nutrient-enriched seawater mesocosms: changes in the genetic diversity of bacterial populations. FEMS Microbiol. Ecol. 34, 243–253. doi: 10.1111/j.1574-6941.2001.tb00775.x
Singer, E., Webb, E. A., Nelson, W. C., Heidelberg, J. F., Ivanova, N., Pati, A., et al. (2011). Genomic potential of Marinobacter aquaeolei, a biogeochemical “Opportunitroph.” Appl. Environ. Microbiol. 77, 2763–2771. doi: 10.1128/AEM.01866-10
Souza, V., Espinosa-Asuar, L., Escalante, A. E., Eguiarte, L. E., Farmer, J., Forney, L., et al. (2006). An endangered oasis of aquatic microbial biodiversity in the Chihuahuan desert. Proc. Natl. Acad. Sci. U.S.A. 103, 6565–6570. doi: 10.1073/pnas.0601434103
Souza, V., Siefert, J. L., Escalante, A. E., Elser, J. J., and Eguiarte, L. E. (2012). The Cuatro Ciénegas Basin in Coahuila, Mexico: an astrobiological Precambrian Park. Astrobiology 12, 641–647. doi: 10.1089/ast.2011.0675
Tilman, D., Mattson, M., and Langer, S. (1981). Competition and nutrient kinetics along a temperature gradient: an experimental test of a mechanistic approach to niche theory. Limnol. Oceanogr. 26, 1020–1033. doi: 10.4319/lo.19184.108.40.2060
Validov, S., Mavrodi, O., De La Fuente, L., Boronin, A., Weller, D., Thomashow, L., et al. (2005). Antagonistic activity among 2,4-diacetylphloroglucinolproducing fluorescent Pseudomonas spp. FEMS Microbiol. Lett. 242, 249–256. doi: 10.1016/j.femsle.2004.11.013
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Weiser, R., Donoghue, D., Weightman, A., and Mahenthiralingam, E. (2014). Evaluation of five selective media for the detection of Pseudomonas aeruginosa using a strain panel from clinical, environmental and industrial sources. J. Microbiol. Methods. 99, 8–14. doi: 10.1016/j.mimet.2014.01.010
Keywords: mesocosm, nutrient enrichment, Cuatro Cienegas Basin, proteobacteria interactions, community structure
Citation: Ponce-Soto GY, Aguirre-von-Wobeser E, Eguiarte LE, Elser JJ, Lee ZM-P and Souza V (2015) Enrichment experiment changes microbial interactions in an ultra-oligotrophic environment. Front. Microbiol. 6:246. doi: 10.3389/fmicb.2015.00246
Received: 15 October 2014; Accepted: 13 March 2015;
Published: 01 April 2015.
Edited by:Michael Travisano, University of Minnesota, USA
Reviewed by:Scott Clingenpeel, Department of Energy Joint Genome Institute, USA
Jillian Petersen, Max Planck Institute for Marine Microbiology, Germany
Copyright © 2015 Ponce-Soto, Aguirre-von-Wobeser, Eguiarte, Elser, Lee and Souza. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Valeria Souza, Laboratorio de Ecología Molecular y Experimental, Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Circuito Exterior S/N, AP 70-275, Coyoacán DF 04510, México firstname.lastname@example.org