Short-term responses of unicellular planktonic eukaryotes to increases in temperature and UVB radiation
- Isabelle Domaizon1Email author,
- Cécile Lepère2,
- Didier Debroas2,
- Marc Bouvy3,
- Jean Francois Ghiglione4,
- Stephan Jacquet1,
- Yvan Bettarel3,
- Corinne Bouvier3,
- Jean Pascal Torréton3,
- Francesca Vidussi3,
- Behzad Mostajir3,
- Amy Kirkham1,
- Emilie LeFloc’h3,
- Eric Fouilland3,
- Helene Montanié5 and
- Thierry Bouvier3
© Domaizon et al.; licensee BioMed Central Ltd. 2012
Received: 22 March 2012
Accepted: 27 August 2012
Published: 11 September 2012
Small size eukaryotes play a fundamental role in the functioning of coastal ecosystems, however, the way in which these micro-organisms respond to combined effects of water temperature, UVB radiations (UVBR) and nutrient availability is still poorly investigated.
We coupled molecular tools (18S rRNA gene sequencing and fingerprinting) with microscope-based identification and counting to experimentally investigate the short-term responses of small eukaryotes (<6 μm; from a coastal Mediterranean lagoon) to a warming treatment (+3°C) and UVB radiation increases (+20%) at two different nutrient levels. Interestingly, the increase in temperature resulted in higher pigmented eukaryotes abundances and in community structure changes clearly illustrated by molecular analyses. For most of the phylogenetic groups, some rearrangements occurred at the OTUs level even when their relative proportion (microscope counting) did not change significantly. Temperature explained almost 20% of the total variance of the small eukaryote community structure (while UVB explained only 8.4%). However, complex cumulative effects were detected. Some antagonistic or non additive effects were detected between temperature and nutrients, especially for Dinophyceae and Cryptophyceae.
This multifactorial experiment highlights the potential impacts, over short time scales, of changing environmental factors on the structure of various functional groups like small primary producers, parasites and saprotrophs which, in response, can modify energy flow in the planktonic food webs.
KeywordsSmall eukaryotes Molecular diversity Temperature UVB radiation Microcosms experiment Mediterranean lagoon
Small-sized plankton plays critical roles in aquatic systems, mostly as major contributors to production and biomass, and as key players driving carbon and nutrient cycles [1, 2]. The study of the gene coding for 18S rRNA has brought opportunities to investigate the eukaryotic composition in the smallest size fraction in various aquatic systems, independently of morphological identification and cultivation [3–7]. The molecular characterization of small (pico and/or nano) eukaryotic assemblages has highlighted an unexpected phylogenetic and functional diversity (e.g. [8–11]), and many important questions are now emerging about the in situ dynamics of diverse eukaryotic groups, and the regulatory factors that drive changes in their structure.
A few studies have investigated the effects of structuring factors on the molecular diversity of small eukaryotes, and shown that trophic status, predation by met zooplankton, and/or viral lytic activity are involved in the regulation of the eukaryotic microbial assemblage [5, 12–15]. However, combined effects of physical factors, such as water temperature and UVB radiation (UVBR: 280–320 nm) are still poorly investigated. It is recognized that either temperature or UVBR increases can modify microbial dynamics and structure at various levels (species, population, trophic network) (e.g. [16–20]). Nevertheless, previous investigations have generally focused on only one specific stressor and little is known about the combined effects of climatic and anthropogenic stressors on diversity and food web structure. Since these stressors are expected to exert complex interactive effects [21–23], multi-factorial studies are required to improve the understanding of the mechanistic basis underlying ecological responses of planktonic food webs to these regulatory factors. A series of enclosure experiments using natural microbial communities from the Mediterranean Thau lagoon were recently performed to assess the response of microbial communities to top-down and bottom-up control under various simulated climatic conditions (temperature and UVBR) . This study showed a much larger effect of temperature than UVBR on bacterial dynamics. In addition to this study, in order to describe the composition of small eukaryotes and potentially to observe changes in their structure, we used a similar microcosm experiment to tease apart the effects of single and combined increase of temperature (+3°C) and UVBR (+20%), at two different nutrients levels. Here, we investigate short-term responses of both pigmented and non-pigmented small eukaryotes (size fraction <6 μm) to these simulated climatic conditions by using morphological and molecular methods (18S rRNA gene sequencing and a fingerprint technique: Capillary Electrophoresis Single Strand Conformation Polymorphism CE-SSCP).
The increases in temperature and UVBR tested in this study correspond to the mean temperature increase expected in the Mediterranean region by 2080–2099 (IPCC 2007) and the high-UVBR scenario for the European region during spring in future years . This approach enables us to describe the short term responses of eukaryotic community assemblages when exposed to these drivers during the productive spring season. The changes induced by these regulatory factors could be detected at different taxonomic levels thanks to the coupling of morphological and molecular approaches.
The four-day experiment (20–23 April 2006) was conducted on the Mediterranean platform for Marine Ecosystems Experimental Research (MEDIMEER) located in Sète (France) on the shore of the Thau lagoon (43°24’49”N, 3°41’19”E). The experimental platform was composed of submerged enclosures (1.2 m diameter and 2 m depth) which allowed the isolation of up to 2,000 L and the simulation of UVBR and temperature increases in order to study the responses of pelagic communities to these manipulated factors simultaneously. The regulations of UVBR and temperature are performed with high frequency monitoring following the in situ temperature and natural incident UVBR (see details in supplementary data; full description in Nouguier et al. ).
In order to fill the 24 Whirl Pack bags, 100 L subsurface lagoon water was pumped and pre-filtered through 6-μm-pore-size polycarbonate membranes (47 mm in diameter) in order to isolate the smallest planktonic fraction. This water sample (<6 μm) was equally distributed into 24 sterile Whirl Pack® polyethylene bags. 12 of these experimental bags received nutrients addition at time zero, while the others were kept without nutrient addition. The set bags which represented the enriched nutrient conditions were obtained by addition of a mixture of leucine (C and N) and phosphate in order to maintain a substrate C:N:P molar ratio close to that of marine bacteria  as described in Bouvy et al.. The bags with and without nutrient addition exhibited concentrations of 0.20 μM and 0.07 μM of PO4, respectively. The two levels of P concentration mimicked natural fluctuations in coastal lagoon waters. These concentrations were chosen to be relevant to phosphorus concentrations recently measured in Thau lagoon (a general decrease over the past 30 years has led to low values of soluble reactive phosphorus: i.e. from 3 μM to undetectable values (<0.03 μM in winter) ).
Since nutrients usually refer to inorganic nutrients, it should be noted that in this study, “nutrients” actually refer to “nutrients and organic source of C and N”.
The incubation time (4 days) and experimental volume (2 L) were chosen to be consistent with the generation time of microorganisms, as validated in other experimental studies [28–31]. In the present study, the composition of unicellular eukaryotes was studied at T0 and T96h. The data provided by Bouvy et al. regarding the evolution of abundances of the main biological communities (i.e. bacteria, viruses, heterotrophic flagellates) at 3 sampling times (T0, T48h, T96h) under the same experimental conditions as ours, informed this choice.
Measurement of abiotic parameters
Environmental conditions (temperature, salinity, chlorophyll a concentration, natural UVBR intensities) during the four days experiment
Environmental conditions during the 4 days of study
Spring (18–24 April)
In situ Temperature
15.7°C to 17.2°C
In situ Salinity
In situ Chl a
Approx. 1 μg/L
In situ maximum UVBR incidentsN (local zenith time)
150 to 185 μW/cm2
Bacterial and viral counting by flow cytometry
At T0 and T96h, 5 ml of water was collected from each of the polyethylene bags for flow cytometry counts. Picocyanobacteria, heterotrophic bacteria and viruses were counted using a FACSCalibur flow cytometry (Becton Dickinson) equipped with an air-cooled laser providing 15 mW at 488 nm. For photosynthetic-cells (i.e. picocyanobacteria) neither fixative nor fluorochrome were used. Samples were stored at <4°C until analysis, which was performed within 2 h of sampling in field laboratories. Analysis was therefore performed on fresh samples, to which a suspension of 1-μm beads (Molecular probes) was added, generally for 4 to 8 minutes in order to obtain >20,000 events. For the analysis of bacteria and viruses, 1 mL fixed (glutaraldehyde 0.5% final concentration) sub-samples were incubated with SYBR Green I (Molecular Probes, Eugene, OR, USA) at a final concentration of 1/10,000 for 15 min at room temperature in the dark. The cytometry flow counts were performed as described in Brussard et al. .
Small eukaryotes microscopy observation
For enumeration of non-pigmented and pigmented eukaryotes, water samples (100 mL) taken at T0 and T96h were fixed with glutaraldehyde (1% final concentration) and stored at 4°C for 24 h. 20 to 25 ml of each preserved water sample was stained with DAPI (final concentration, 15 μg mL−1) for 15 min, filtered onto a black Nuclepore filter (0.8 μm-pore-size), stored at −20°C, and counted under an epifluorescence microscope with UV excitation (modified from Boenigk et al.). Under UV light (350/461 nm), the eukaryotic cell nucleus appears as a separate organelle, while prokaryotic organisms appear as cells uniformly stained without visible nuclei. The blue and green light excitations were used to reveal pigmented cells.
Molecular analysis of small eukaryotes
Sampling and preservation
Water samples from each treatment were taken at the beginning and at the end of the experiment. The microbial biomass was collected on 0.2 μm pore size polycarbonate membranes (Millipore) under very low vacuum (<20 mbar) to prevent cell damage. Filters were then stored at −80°C until nucleic acid extraction.
Nucleic acid extraction
Nucleic acid extraction was performed as described by Lefranc et al.  and extracts were stored at −20°C until analysis.
Capillary electrophoresis – single strand conformation polymorphism (CE-SSCP)
Nucleic acids from each sample were used as templates for PCR amplification of the 18S rRNA gene with primers Uni1392r (5’-ACG-GGC-GGT-GTG-TRC-3’) labelled at the 5’-end with phosphoramidite  and Euk1209f (5’-CAG-GTC-TGT-GAT-GCC-CGC-3’) . Each 25 μL reaction mixture contained 50 μM of each primer, 1X Pfu reaction buffer, 20 mM dNTPs, 1.0 U of Pfu DNA polymerase (Promega) and 0.1 μg of template DNA. PCR amplification was performed with a Rob cycler (Stratagene) under the following conditions: an initial denaturation step of 94°C for 2 min, followed by 10 touchdown cycles of denaturation at 94°C for 1 min, annealing at 65°C (with the temperature decreasing 1°C each cycle) for 1 min, and extension at 72°C for 1 min, followed by 15 cycles of 94°C for 1 min, 55°C for 1 min and 72°C for 1 min, and a final elongation step at 72°C for 10 min. The TET-labelled PCR products were quantified by visualization in ethidium bromide-stained agarose gels (2%) and diluted in sterile TE (10 mM Tris, 1 mM EDTA) in order to obtain around 10 ng mL–1 of PCR product. One μL of the dilution was mixed with 18.9 μL of formamide (Applera Corp. Norwalk, Connecticut) and 0.1 μL of the internal size standard Gene-Scan-400 Rox (Applied Biosystems), denatured at 94°C for 5 minutes, and immediately cooled on ice for 10 minutes before electrokinetic injection (5 s, 12 kV) into a capillary tube (47 cm x 50 μm) filled with 5.6% of Gene Scan polymer in a ABI Prism 310 Genetic analyser (Applied Biosystems). Electrophoresis was carried out and data were collected as described in Sauret et al. .
Eukaryotic rRNA genetic libraries
Environmental DNA extracts were also used to construct the 18S rRNA gene clone libraries. The eukaryote-specific primers Ek-1 F (5’-CTG-GTT-GAT-CCT-GCC-AG-3’) and Ek-1520R (5-CYG-CAG-GTT-CAC-CTA-C-3’) were used for PCR amplification . The PCR mixture (50 μL) contained about 10 ng of environmental DNA, 200 μM of each deoxynucleoside triphosphate, 2 mM MgCl2, 10 pmol of each primer, 1.5 U of Taq DNA polymerase (Eurobio) and the PCR buffer supplied with the enzyme. Reactions were carried out in an automated thermocycler (MJ Research PTC 200-cycler) with the following cycle: initial denaturation at 95°C for 5 min, 30 cycles of denaturation at 95°C for 1 min, annealing at 57°C for 1 min, and extension at 72°C for 1 min 30 s, and a final extension at 72°C for 10 min. PCR products (at least four 50 μL samples) from the triplicate samples of each experimental condition were pooled, precipitated with ethanol–sodium acetate and re-suspended in 50 μL of sterile water. Clone libraries were constructed for the T0 control and for each of the eight treatments at T96 h using a TOPO TA cloning kit (Invitrogen, Carlsbad, CA) with PCR vector 2.1 according to the manufacturer’s instructions.
DOTUR was used to determine operational taxonomic units (OTUs) from 18S sequences data  with a cut-off of 97% sequence similarity. To determine the phylogenetic affiliation, each sequence was first compared with sequences available in public databases using BLAST (National Center for Biotechnology Information and the Ribosomal Database Project) . Secondly, the OTUs were aligned with complete sequences in an ARB database using the latter’s automatic alignment tool (http://www.arb-home.de) . The resulting alignments were checked and corrected manually. Sequences were inserted into an optimised tree according to the maximum parsimony criteria without allowing any changes to the existing tree topology (ARB software). The resulting tree was pruned to retain the closest relatives, sequences representative of eukaryotic evolution and our clones (Additional file 1: Figure S1). The sequences were screened for potential chimeric structures by using Chimera check from Ribosomal Database project II and by performing fractional treeing of the 5' and 3' ends of the sequenced DNA fragments. The sequences reported in this paper have been deposited into Genbank (accession numbers: HQ393974 to HQ394162).
The relative distribution of OTUs in the library was used to calculate coverage values (Good’s coverage)  and the non-parametric richness estimator Chao1  and ACE  which are the most appropriate indices for microbial clone libraries .
We tested the homogeneity of the main biological parameters in experimental bags at the initial point (T0) of the experiment using an ANOVA test.
To test the effects of temperature, UV and nutrients on the abundance of all biological groups (bacteria, picocyanobacteria, viruses, heterotrophic flagellates and pigmented eukaryote abundances at T96 h), we used a three-way ANOVA test (with Bonferroni adjustment). Equality of the variances and normality of the residuals were tested by Bartlett and Shapiro-Wilk tests. The software SigmastatTM 3.1 was used for all analyses.
Indirect multivariate analysis was used to compare CE-SSCP fingerprinting. Total fingerprinting area was normalized between the different CE-SSCP profiles using the internal size standard Gene-Scan-400 Rox using the SAFUM software . Similarity matrices based on Bray-Curtis distances, dendrograms (complete linkage clustering) and ordination by non-metric multidimensional scaling (MDS) were then obtained by using the PRIMER 5 software (PRIMER-E, Ltd., UK). One-way analysis of similarity (ANOSIM, Primer-E) was performed on the same distance matrix to test the null hypothesis that there was no difference between eukaryotic communities from replicate samples of each condition.
Statistics applied to phylogenetic information
From the sequencing results, the beta-diversity was studied from the Unifrac distance (fraction of the total branch length in the phylogeny that is unique to each environment) of each sample. In order to compare eukaryotic communities from the 9 genetic libraries Unifrac (http://bmf2.colorado.edu/unifrac/index.psp; ) metrics were used to perform a principal coordinate analysis (PCA). The P-values matrix that compares each sample to each other sample was also performed from UNIFRAC metrics.
To investigate the relationships between changes in the eukaryote community structure (number of clones affiliated to each OTUs within main phylogenetic groups) and physic-chemical and biological parameters, we used direct multivariate canonical correspondence analysis (CCA) . In addition to temperature values, UVB radiation, and nutrient concentrations, we considered the abundances of bacteria, picocyanobacteria, viruses, pigmented eukaryotes and heterotrophic flagellates as explanatory variables. CCA was calculated for the T96 h dataset using the Vegan package within the R software (http://cran.rproject.org/). A minimal set of explanatory variables associated with variation in eukaryote community structure was identified, allowing us to exclude the most redundant explanatory variables. Forward selection was performed to identify environmental variables that could explain a significant portion of the variation in small eukaryote structure (P < 0.05) at T96 h. Eigen values for site scores, biplot and diversity data were plotted to illustrate the associations between these data .
Biological and chemical parameters
Initial conditions for chemical and biological parameters
Chemical and biological parameters in experimental bags at T0
No nutrient addition
Bacteria 106 cell mL -1*
Virus 108 cell mL-1*
Picocyanobacteria 103 cell mL-1*
Non-pigmented Euk. 102 cell mL -1
Pigmented Euk. 103 cell mL -1
Abundances and structure of the small eukaryotic community
The mean abundance of non-pigmented eukaryotes was 776 cells mL-1 at T0, accounting for about 15% of total eukaryotes. In comparison to microscope counting, the proportion of typical non-pigmented eukaryotes was over-estimated in the clone library, accounting for 43.2% of total clones (such over-representation of non-pigmented groups in 18S rRNA gene clone libraries has been discussed previously e.g.[50–52]). The diversity of these non-pigmented groups cannot be discriminated by classical microscopy due to a lack of distinct morphological features and/or their small size. However, from cloning-sequencing results, 11 different OTUs could be attributed to non-pigmented groups: Cercozoa (2 OTUs), Stramenopiles affiliated to Hyphochytrids (1 OTU), Syndiniales affiliated to Amoebophrya (2 OTUs), uncultured alveolates (4 OTUs), and Choanoflagellida (2 OTUs) (Figure 2 Additional file 2: Table S1).
Changes in the small eukaryotes structure under UVBR, temperature and/or nutrient increase
Changes in abundances
Results of the three-way ANOVA performed from T96h abundance values
Anova results (P)
Temp x UV
Temp x Nut
Temp x UV
Temp x UV x Nut
Pigmented eukaryotes (total) cells mL -1
Non-pigmented eukaryotes cells mL -1
Bacteria cell mL -1
Virus particles mL -1
Picocyanobacteria cells mL -1
Some major changes were observed in the relative proportions of the main taxonomic groups. The abundance of pigmented Dinophyceae increased in all treatments, with the highest increases where nutrients were added. Indeed, the 3-way ANOVA showed a significant effect of nutrients (p = 0.028, Table 3). Inversely, for Cryptophyceae, a general negative impact of nutrient addition (p < 0.001) counteracted the positive impact of temperature increase (Table 3, Figure 2). The relative abundance of Mamiellophyceae (Micromonas and Ostreococcus) decreased from T0 to T96h in all treatments, and they represented only between 0.1 and 14.8% of pigmented eukaryotes at the end of the experiment (depending on the treatment). Pyramimonadales seemed to take advantage of the general reduction of Mamiellophyceae densities and developed strongly, especially in treatments with increased UVBR. The 3-way ANOVA showed a positive impact of UVBR on Pyramimonadales abundance.
Non-pigmented eukaryotes (mainly free flagellated forms) tended to increase in abundance in all conditions. The highest values were found in TUV + Nut treatments (mean abundance: 2.5 x103 cells mL-1), however, the 3-way ANOVA did not reveal any significant impact of the manipulated factors (Table 3).
Changes in small eukaryotes structure (CE SSCP)
Changes in small eukaryotes phylogenetic composition (sequencing)
The OTUs affiliated to non-pigmented taxa generally dominated the clone libraries (from 67.6% in C + Nut to 85.3% in UV + Nut; Figure 4 and Additional file 2: Table S1). Among them, Ciliates and uncultured Alveolates were generally well represented (accounting from 14 to 32% of total OTUs, and from 13 to 37% of clones, according to the treatments). However, the increase of non-pigmented group proportions within most of the libraries (compared to T0) was mainly linked to the emergence of taxa affiliated to parasitic groups: Hyphochytrids and genus Pirsonia (Heterokonta), and Amoebophrya (Alveolata).
The CCA supported the conclusion obtained from the UNIFRAC analysis, clearly showing that all treatments with increased temperature grouped together. Furthermore, the highest abundances of bacteria, picocyanobacteria, and pigmented groups such as Cryptophyceae and Bacillariophyceae were tightly associated with treatments receiving an increased temperature (Figure 5).
The CCA plot also illustrates the strong negative impact of experimental conditions on Mamiellophyceae in general. Mamiellophyceae represented 28% of sequences in the clone library at T0, but were not detected at T96 h (except 1 OTU detected in the C treatments). In contrast, Pyramimonadales sequences (2 OTUs) appeared at T96 h in 6 out of the 8 types of treatment.
Overall, the analysis of the OTUs dynamics (either generally or for specific phylogenetic groups) showed that, even when the abundance of a given group did not change significantly from one treatment to another, some rearrangements could occur at the OTUs level (Additional file 2: Table S1). The CCA showed that 18.8% of the total variation in the eukaryotic structure was explained by temperature, whereas, UVBR and nutrients explained 11% and 8.4%, respectively.
The Thau lagoon, characterised by a high abundance of small eukaryotes and by recent in situ changes in phytoplankton structure due to water temperature increase , is an interesting ecosystem to investigate the responses of small eukaryotes to climatic and anthropogenic regulatory factors. Our experimentation does not intend to predict the impact of long-term global change on the structure of small planktonic eukaryotes. Indeed, only a combination of approaches including laboratory studies on model microbes, microcosm and mesocosm experiments, and in situ comparative studies would help to provide realistic predictions of the effects of environmental changes [23, 54]. Our goal was to reveal the potential rapid responses of small eukaryote assemblage (using molecular and morphological methods) during the productive spring season when plankton may be particularly vulnerable to elevated temperature and UVBR .
Molecular analyses revealed the presence of various phylogenetic groups within the “black box” of small eukaryotes, especially non-pigmented eukaryotes (poorly discriminated by microscopy). Some limitations in the PCR-based methods are recognized, for instance, the over-representation of Alveolata (particularly Dinoflagellates and Ciliates) in 18S rRNA gene clone libraries due their high SSU rRNA gene copy number [50–52]. However, the molecular methods used here enabled the description of the diversity within dominant eukaryotic populations, and allowed examination of the effects of regulatory factors by considering both the dynamics of OTUs (using the sequencing and fingerprinting datasets) and the comparison of phylogenetic composition obtained for all treatments (using the sequencing data). The impact of temperature nutrients and UVBR explained 18.8%, 11.0% and 8.4% of the variance of the small eukaryotes structure respectively. While Bouvy et al. (2011) could not detect any significant responses of pico- or nano-eukaryotic plankton in the same experimental conditions, we demonstrated here, at a different taxonomic resolution, that small eukaryotes community structure was actually affected by this multi-factorial pressure.
The simultaneous use of molecular and morphological methods was therefore essential to provide evidence of rapid shifts that occur at various taxonomic levels (abundance of large groups or community composition at OTU level) under the influence of temperature, UVBR and nutrient treatments.
Among the 3 regulatory factors tested, both sequencing and CE-SSCP demonstrated that increased temperature had the greatest influence on the small eukaryote community structure and composition. The single effect of temperature (without any significant interaction with UVBR and nutrients) on total pigmented eukaryote abundance was observed by microscopy. Considering the different phylogenetic groups within pigmented eukaryotes, complex interaction effects were also suggested. For instance, our results showed that under multi-factorial environmental changes, the general impact on the molecular diversity and abundance of pigmented Dinophyceae resulted from complex interactive (non-additive) effects. Multi-factorial interactions were also apparent for Cryptophyceae which experienced antagonistic effects of nutrient addition (significantly negative impact) and temperature (positive impact on relative abundance).
In addition to the manipulated factors (temperature, UVBR and nutrients), some biotic interactions such as predation, viral lysis and competition, are involved in the responses observed in this experiment. For example, the general reduction of Mamiellophyceae (Micromonas and Ostreococcus) in all treatments might be linked to (i) manipulation effects since these fragile cells might have been affected by filtration steps, (ii) limitation by inorganic nutrients under the rather low orthophosphate concentrations at T96h (from 0.05 to 0.08 μM of PO4), (iii) the grazing impact of heterotrophic flagellates: these microorganisms are known to play a significant role in the regulation of Ostreococcus populations in the Thau lagoon  and were shown to exert a strong control of bacterioplankton during the study period . We could not detect a link between the dynamics of Micromonas/Ostreococcus and viruses. Since biological descriptors can explain some of the variance of small eukaryote phylogenetic structure, the observations made here regarding the effect of temperature, UVBR, and nutrients have to be considered in view of possible biological effects. Predation by zooplankton and competition with larger phytoplanktonic species were not considered in our size fractionated approach and should be taken into account, especially if long-term extrapolation of in situ responses of small eukaryotes is considered.
Our data provide further illustration of the need to consider the taxonomic and functional diversity of heterotrophic flagellates. The lack of discrimination between heterotrophic bacterivores and parasitic/saprotrophic zoospores within the non-pigmented flagellates can lead to misinterpretation of the functioning and responses of planktonic food webs. Indeed, while microscope observations did not allow us to detect changes in the abundance and structure of non-pigmented eukaryotes, a structuring impact of manipulated factors (especially temperature) was observed through sequencing results on taxa affiliated to parasitic and saprotroph groups (particularly Syndiniales and Hyphochytrids). The existence of eukaryotic parasites among small-size plankton was recently re-discovered by molecular environmental surveys, and the ecological significance of these groups has been highlighted by several authors [57, 58]. The ‘Fungi-like’ Hyphochytrids possess many morphological and ecological similarities to chytrids [58, 59], and their role as saprotrophs and/or parasites is unclear [60, 61], whereas the Amoebophrya are well recognized as a widely distributed parasitic order within the Dinophyceae . Amoebophrya and Hyphochytrids emerged in clone libraries at T96 h and were presumably present among the rare species at T0. The taxa found to be phylogenetically close to Amoebophrya particularly emerged in treatments with increased temperature (Figure 5), along with their hosts (pigmented Dinoflagellates). This observation supports Guillou et al.’s  suggestion that warming could promote rapid infection cycles of Amoebophrya. However, broad extrapolation would need to take into account various aspects of the host-parasite relationships, such as the mechanisms underlying the parasitic specificity. In contrast to the Amoebophrya, hyphochytrids were associated with all treatments except those with increased temperature (Figure 5). From our results, we hypothesized that not only parasite communities, but also saprotroph communities would be shaped by temperature and UVBR conditions, as already described in other ecosystems . The responses of saprotrophs to these drivers may result from direct and/or indirect effects as demonstrated in soils ; further research is probably needed on the saprotrophs in aquatic systems since changes in their assemblages may influence organic matter decomposition and nutrient cycling.
Even though caution should be exercised when applying the results of small-scale experiments to larger-scale systems, these results can be treated as an insight into ecological interactions that may occur in larger natural systems with more complex planktonic assemblages. Our results indicate that these ecosystem drivers, which are associated with climate change, and their interactions may cause changes in small eukaryotic community abundance and structure involving various functional groups including the small primary producers, parasites and saprotrophs. Notably, temperature tends to have a much greater effect on the community composition of small eukaryotes compared to UVBR (at least at the level tested in our experiment). Due to their strong link with other communities within the food web, the small eukaryotes variability may have potential consequences in food webs structure and energy flow. Currently, our knowledge of the potential for plankton in general and small eukaryotes in particular to adapt genetically and phenotypically to multifactorial physico-chemical climate drivers is poor. To improve our understanding, additional experimental investigations in other types of ecosystems and over longer periods of warming and UVBR exposure are required before generalization may be confidently applied. Future investigations should be based on the coupling of methods such as microscopy, flow cytometry, molecular analyses targeting several gene markers or fluorescence in situ hybridization in order to analyse the responses of the microbial community structure to multiple stressors at various taxonomic levels.
We gratefully acknowledge Jean Nouguier and Yvan Vergne for their technical help during the experiment. This study was supported by the French program PNEC (10301705 to TB) and the ANR AQUAPHAGE (ANR07 BIODIV 015–02 to TB). This work was also supported by the ‘Groupement De Recherches (GDR) 2476 Réseaux Trophiques Pélagiques. The experimental platform for Mediterranean Ecosystem Research (MEDIMEER) was funded by UMR 5119 ECOLAG, CNRS-INEE, Institut Fédératif de Recherche 129 A. Sabatier, GDR 2476 Réseaux Trophiques Aquatiques, Région Languedoc Roussillon.
We thank Joseph Kirkman for improving the text.
- Caron DA, Peele ER, Lim EL, Dennett MR: Picoplankton and nanoplankton and their trophic coupling in the surface waters of the Sargasso Sea south of Bermuda. Limnol Oceanogr. 1999, 44: 259-272. 10.4319/lo.1999.44.2.0259.View ArticleGoogle Scholar
- Li WKW: Primary production of prochlorophytes, cyanobacteria and eucaryotic ultraphytoplankton: measurements from flow cytometric sorting. Limnol Oceanogr. 1994, 39: 169-175. 10.4319/lo.1994.39.1.0169.View ArticleGoogle Scholar
- Moreira D, Lόpez-García P: The molecular ecology of microbial eukaryotes unveils a hidden world. Trends Microbiol. 2002, 10: 31-38. 10.1016/S0966-842X(01)02257-0.PubMedView ArticleGoogle Scholar
- Šlapeta J, Moreira D, Lopez-Garcia P: Global dispersal and ancient cryptic species in the smallest marine eukaryotes. Mol Biol Evol. 2006, 23 (1): 23-29.PubMedView ArticleGoogle Scholar
- Lepère C, Boucher D, Jardillier L, Domaizon I, Debroas D: Succession and regulation factors of small eukaryote community composition in a lacustrine ecosystem (Lake Pavin). Appl Environ Microb. 2006, 72: 2971-2981. 10.1128/AEM.72.4.2971-2981.2006.View ArticleGoogle Scholar
- Jardillier L, Zubkov MV, Pearman J, Scanlan DJ: Significant CO2 fixation by small prymnesiophytes in the subtropical and tropical northeast Atlantic Ocean. ISME J. 2010, 4: 1180-1192. 10.1038/ismej.2010.36.PubMedView ArticleGoogle Scholar
- Cuvelier ML, Allen AE, Monier A, McCrow JP, Messie M: Targeted metagenomics and ecology of globally important uncultured eukaryotic phytoplankton. PNAS. 2010, 107 (33): 14679-14684. 10.1073/pnas.1001665107.PubMedPubMed CentralView ArticleGoogle Scholar
- Lepère C, Domaizon I, Debroas D: Unexpected importance of potential parasites in the composition of the freshwater small-eukaryote community. Appl Environ Microb. 2008, 74: 2940-2949. 10.1128/AEM.01156-07.View ArticleGoogle Scholar
- Amaral-Zettler LA, McCliment EA, Ducklow HW, Huse SM: A method for studying protistan diversity using massively parallel sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA genes. PLoS One. 2009, 4 (7): e6372-10.1371/journal.pone.0006372.PubMedPubMed CentralView ArticleGoogle Scholar
- Massana R, Unrein F, Rodriguez-Martinez R, Forn I, Lefort T, Pinhassi J, Not F: Grazing rates and functional diversity of uncultured heterotrophic flagellates. ISME J. 2009, 3 (5): 588-596. 10.1038/ismej.2008.130.PubMedView ArticleGoogle Scholar
- Shi XL, Marie D, Jardillier L, Scanlan DJ, Vaulot D: Groups without cultured representatives dominate eukaryotic picophytoplankton in the oligotrophic South East Pacific Ocean. PLoS One. 2009, 4 (10): e7657-10.1371/journal.pone.0007657.PubMedPubMed CentralView ArticleGoogle Scholar
- Evans C, Archer SD, Jacquet S, Wilson WH: Direct estimates of the contribution of viral lysis and microzooplankton grazing to the decline of a Micromonas spp population. Aquat Microb Ecol. 2003, 30: 207-219.View ArticleGoogle Scholar
- Baudoux AC, Veldhuis MJW, Noordeloos AAM, Vann Noort G, Brussaard CPD: Estimates of virus-vs grazing induced mortality of picophytoplankton in the North Sea during summer. Aquat Microb Eco. 2008, 52: 69-82.View ArticleGoogle Scholar
- Chen M, Chen F, Zhao B, Wu QL, Kong FX: Genetic diversity of eukaryotic microorganisms in Lake Taihu, a large shallow subtropical lake in China. Microb Ecol. 2008, 56 (3): 572-583. 10.1007/s00248-008-9377-8.PubMedView ArticleGoogle Scholar
- Masquelier S, Foulon E, Jouenne F, Ferréol M, Brussard CPD, Vaulot D: Distribution of eukaryotic in the English Channel and North Sea in summer. J Sea Res. 2011, 66: 111-122. 10.1016/j.seares.2011.05.004.View ArticleGoogle Scholar
- Petchey OL, McPhearson PT, Casey TM, Morin PJ: Environmental warming alters food-web structure and ecosystem function. Nature. 1999, 402: 69-72. 10.1038/47023.View ArticleGoogle Scholar
- Mostajir B, Sime-Ngando T, Demers S, Belzile C: Ecological implications of changes in cell size and photosynthetic capacity of marine Prymnesiophyceae induced by ultraviolet-B radiation. Mar Ecol Prog Ser. 1999, 187: 89-100.View ArticleGoogle Scholar
- Sommaruga R, Hofer JS, Alonso-Saez L, Gasol JM: Differential Sunlight Sensitivity of Picophytoplankton from Surface Mediterranean Coastal Waters. Appl Environ Microb. 2005, 71 (4): 2154-2157. 10.1128/AEM.71.4.2154-2157.2005.View ArticleGoogle Scholar
- Ferreyra GA, Mostajir B, Schloss IR, Chatila K, Ferrario ME, Sargian P, Roy S, Prod'homme J, Demers S: Ultraviolet-B radiation effects on the structure and function of lower trophic levels of the marine planktonic food web. Photochem Photobiol. 2006, 82 (4): 887-897. 10.1562/2006-02-23-RA-810.PubMedView ArticleGoogle Scholar
- Conan P, Joux F, Torréton JP, Pujo-Pay M, Rochelle-Newall E, Mari X: Impact of solar ultraviolet radiation on bacterio- and phytoplankton activity in a large coral reef lagoon (SW New Caledonia). Aquat Microb Ecol. 2008, 52: 83-98.View ArticleGoogle Scholar
- Christensen MR, Graham MD, Vinebrooke RD, Findlay DL, Paterson MJ, Turner MA: Multiple anthropogenic stressors cause ecological surprises in boreal lakes. Global Change Biol. 2006, 12 (12): 2316-2322. 10.1111/j.1365-2486.2006.01257.x.View ArticleGoogle Scholar
- Vidussi F, Mostajir B, Fouilland E, Le Floc’h E: Effects of experimental warming and increased ultraviolet B radiation on the Mediterranean plankton food web. Limnol Oceanogr. 2011, 56 (1): 206-218. 10.4319/lo.2011.56.1.0206.View ArticleGoogle Scholar
- Doyle SA, Saros JE, Williamson CE: Interactive effects of temperature and nutrient limitation on the response of alpine phytoplankton growth to ultraviolet radiation. Limnol Oceanogr. 2005, 50 (5): 1362-1367. 10.4319/lo.2005.50.5.1362.View ArticleGoogle Scholar
- Bouvy M, Bettarel Y, Bouvier C, Domaizon I, Jacquet S, LeFloc’h E, Montanié H, Mostajir B, Sime-Ngando T, Torréton JP, Vidussi F, Bouvier T: Trophic interactions between viruses, bacteria, and nanoflagellates under various nutrient conditions and simulated climate change. Environ Microbiol. 2011, 13 (7): 1842-1857. 10.1111/j.1462-2920.2011.02498.x.PubMedView ArticleGoogle Scholar
- Nouguier J, Mostajir B, Le Floc'h E, Vidussi F: An automatically operated system for simulating global change temperature and ultraviolet B radiation increases: application to the study of aquatic ecosystem responses in mesocosm experiments. Limnol Oceanog Methods. 2007, 5: 269-279.View ArticleGoogle Scholar
- Goldman JC, Caron DA, Dennet MR: Regulation of gross efficiency and ammonium regeneration in bacteria by C:N ratio. Limnol Oceanogr. 1987, 32: 1239-1252. 10.4319/lo.19188.8.131.529.View ArticleGoogle Scholar
- Collos Y, Gagne C, Laabir M, Vaquer A, Cecchi P, Souchu P: Nitrogenous nutrition of Alexandrium catenella (Dinophyceae) in cultures and in Thau lagoon, southern France. J Phycol. 2004, 40 (1): 96-103. 10.1046/j.1529-8817.2004.03034.x.View ArticleGoogle Scholar
- Jardillier L, Boucher D, Personnic S, Jacquet S, Thenot A, Sargos D, Amblard C, Debroas D: Relative importance of nutrients and mortality factors on prokaryotic community composition in two lakes of different trophic status: Microcosm experiments. FEMS Microbiol Ecol. 2005, 53 (3): 429-443. 10.1016/j.femsec.2005.01.011.PubMedView ArticleGoogle Scholar
- Brussard CPD, Marie D, Bratbak G: Flow cytometric detection of viruses. J. Virological methods. 2000, 85: 175-182. 10.1016/S0166-0934(99)00167-6.View ArticleGoogle Scholar
- Pradeep RAS, Sime-Ngando T: Functional responses of prokaryotes and viruses to grazer effects and nutrient additions in freshwater microcosms. ISME J. 2008, 2: 498-509. 10.1038/ismej.2008.15.View ArticleGoogle Scholar
- Bonilla-Findji O, Herndl GJ, Gattuso JP, Weinbauer MG: Viral and flagellate control of prokaryotic production and community structure in Offshore Mediterranean Waters. Appl Environ Microb. 2009, 75 (14): 4801-4812. 10.1128/AEM.01376-08.View ArticleGoogle Scholar
- Tréguer P, LeCorre P: Manuel d'analyse des sels nutritifs dans l'eau de mer. Utilisation de l'AutoAnalyser II Technicon. 1975, Univ. Bretagne Occidentale, Laboratoire de Chimie marine, Brest, France, 2Google Scholar
- Boenigk J, Stadler P, Wiedlroither A, Hahn MW: Strain-specific differences in the grazing sensitivities of closely related ultramicrobacteria affiliated with the Polynucleobacter Cluster. Appl Environ Microb. 2004, 70 (10): 5787-5793. 10.1128/AEM.70.10.5787-5793.2004.View ArticleGoogle Scholar
- Lefranc M, Thénot A, Lepère C, Debroas D: Genetic diversity of small eukaryotes in lakes differing by their trophic status. Appl Environ Microb. 2005, 71: 5935-5942. 10.1128/AEM.71.10.5935-5942.2005.View ArticleGoogle Scholar
- Lane DJ, Pace B, Olsen GJ, Stahl DA, Sogin ML, Pace NR: Rapid determination of 16S ribosomal RNA sequences for phylogenetic analyses. PNAS. 1985, 82: 6955-6959. 10.1073/pnas.82.20.6955.PubMedPubMed CentralView ArticleGoogle Scholar
- Giovannoni SJ, DeLong EF, Olsen GJ, Pace NR: Phylogenetic group-specific oligodeoxynucleotide probes for identification of single microbial cell. J Bacteriol. 1988, 170: 720-726.PubMedPubMed CentralGoogle Scholar
- Sauret C, Christaki U, Moutsaki P, Hatzianestis I, Gogou A, Ghiglione JF: Influence of pollution history on the response of coastal bacterial and nanoeukaryote communities to crude oil and biostimulation assays. Mar Environ Res. 2012, 79: 70-78.PubMedView ArticleGoogle Scholar
- Lopez-Garcia P, Philippe H, Gail F, Moreira D: Autochthonous eukaryotic diversity in hydrothermal sediment and experimental microcolonizers at the Mid-Atlantic Ridge. PNAS. 2003, 100 (2): 697-702. 10.1073/pnas.0235779100.PubMedPubMed CentralView ArticleGoogle Scholar
- Schloss PD, Handelsman J: Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl Environ Microb. 2005, 71 (3): 1501-1506. 10.1128/AEM.71.3.1501-1506.2005.View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Mille W, Lipman DJ: Gapped BLAST and PSIBLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Ludwig W, Strunk O, Westram R, Richter L: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32: 1363-1371. 10.1093/nar/gkh293.PubMedPubMed CentralView ArticleGoogle Scholar
- Hughes JB, Hellmann JJ, Ricketts TH, Bohannan BJ: Counting the uncountable: Statistical approaches to estimating microbial diversity. Appl Environ Microbiol. 2001, 67 (10): 4399-4406. 10.1128/AEM.67.10.4399-4406.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Chao A: Nonparametric estimation of the number of classes in a population. Scandinavian J Stat. 1984, 11: 265-270.Google Scholar
- Chao A, Lee SM: Estimating the number of classes via sample coverage. J Am Stat Assoc. 1992, 87: 210-217. 10.1080/01621459.1992.10475194.View ArticleGoogle Scholar
- Kemp PF, Aller JY: Estimating prokaryotic diversity: When are 16S rDNA libraries large enough?. Limnol Oceanogr: Methods. 2004, 2: 114-125.View ArticleGoogle Scholar
- Zemb O, Haegeman B, Delgenes JP, Lebaron P, Godon JJ: Safum: statistical analysis of SSCP fingerprints using PCA projections, dendrograms and diversity estimators. Mol Ecol Notes. 2007, 7: 767-770. 10.1111/j.1471-8286.2007.01882.x.View ArticleGoogle Scholar
- Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microb. 2005, 71 (12): 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.View ArticleGoogle Scholar
- ter Braak CJF: Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis. Ecology. 1986, 67: 1167-1179. 10.2307/1938672.View ArticleGoogle Scholar
- Legendre P, Legendre L: Numerical ecology. 1998, 2nd English edition Elsevier Science BV, AmsterdamGoogle Scholar
- Not F, del Campo J, Balagué V, de Vargas C, Massana R: New Insights into the Diversity of Marine Picoeukaryotes. PLoS ONE. 2009, 4: e7143-10.1371/journal.pone.0007143.PubMedPubMed CentralView ArticleGoogle Scholar
- Shi XL, Lepère C, Scanlan DJ, Vaulot D: Plastid 16S rRNA Gene Diversity among Eukaryotic Picophytoplankton Sorted by Flow Cytometry from the South Pacific Ocean. PLoS ONE. 2011, 6 (4): e18979-10.1371/journal.pone.0018979.PubMedPubMed CentralView ArticleGoogle Scholar
- Lepère C, Masquelier S, Mangot JF, Debroas D, Domaizon I: Vertical distribution of small eukaryote diversity in lakes: a quantitative approach. The ISME Journal. 2010, 4: 1509-1519. 10.1038/ismej.2010.83.PubMedView ArticleGoogle Scholar
- Clarke KR, Warwick R: Change in Marine Communities: An Approach to Statistical Analysis and Interpretation. 2001, 2nd edition: PRIMER-E, Plymouth, UKGoogle Scholar
- Joint I, Donay SC, Karl DM: Will ocean acidification affect marine microbes?. The ISME J. 2011, 5: 1-7. 10.1038/ismej.2010.79.PubMedView ArticleGoogle Scholar
- Joint I, Jordan MB: Effect of short-term exposure to UVA and UVB on potential phytoplanlton production in UK coastal waters. J Plankton Res. 2008, 3052: 199-210.Google Scholar
- Bec B, Husseini-Ratrema J, Collos Y, Souchu P, Vaquer A: Phytoplankton seasonal dynamics in a Mediterranean coastal lagoon: emphasis on the picoeukaryote community. J Plankton Res. 2005, 27 (9): 881-894. 10.1093/plankt/fbi061.View ArticleGoogle Scholar
- Guillou L, Alves-de Souza C, Siano R, Gonzalez H: The ecological significance of small eukaryotic parasites in marine ecosystems. Microbiol Today. 2010, 92-95.http://www.sgm.ac.uk/pubs/micro_today/about.cfm,Google Scholar
- Lefèvre E, Roussel B, Amblard C, Simé-Ngando T: The molecular diversity of freshwater picoeukaryotes reveals high occurrence of putative parasitoids in the plankton. PLoS ONE. 2008, 3: 2324-2333. 10.1371/journal.pone.0002324.View ArticleGoogle Scholar
- Hausner G, Inglis G, Yanke LJ, Kawchuk LM, McAllister TA: Analysis of restriction fragment length polymorphisms in the ribosomal DNA of a selection of anaerobic chytrids. Can J Bot. 2000, 78 (7): 917-927.Google Scholar
- Alster A, Zohary T: Interactions between the bloom-forming dinoflagellate Peridinium gatunense and the chytrid fungus Phlyctochytrium sp. Hydrobiologia. 2007, 578 (1): 131-139. 10.1007/s10750-006-0439-y.View ArticleGoogle Scholar
- Ibelings B, Arnout De Bruin W, Kagami M, Rijkeboer M, Brehm M, Van D, Ibelings B, Arnout De Bruin W, Kagami M, Rijkeboer M, Brehm M, Van Donk E: Host parasite interactions between freshwater phytoplankton and chytrid fungi (chytridiomycota). J Phycol. 2004, 40: 437-453. 10.1111/j.1529-8817.2004.03117.x.View ArticleGoogle Scholar
- Guillou L, Viprey M, Chambouvet A, Welsh RM, Kirkham AR, Massana R, Scanlan DJ, Worden AZ: Widespread occurrence and genetic diversity of marine parasitoids belonging to Syndiniales (Alveolata). Environ Microiol. 2008, 10 (12): 3349-3365. 10.1111/j.1462-2920.2008.01731.x.View ArticleGoogle Scholar
- Reuder J, Dameris M, Koepke P: Future UVradiation in Central Europe modeled from ozone scenarios. J Photoch Photobio B. 2001, 61: 94-105. 10.1016/S1011-1344(01)00143-9.View ArticleGoogle Scholar
- Duguay KJ, Kliromonos JN: Direct and indirect effects of enhanced UV-B radiation on the decomposing and competitive abilities of saprobic fungi. Applied Soil Ecol. 2000, 14 (2): 157-164. 10.1016/S0929-1393(00)00049-4.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.