BACKGROUND
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  <cit> .
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  <cit> .
the molecular characterization of small  eukaryotic assemblages has highlighted an unexpected phylogenetic and functional diversity , 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  <cit> .
however, combined effects of physical factors, such as water temperature and uvb radiation  are still poorly investigated.
it is recognized that either temperature or uvbr increases can modify microbial dynamics and structure at various levels  .
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  <cit> , 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   <cit> .
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  and uvbr , at two different nutrients levels.
here, we investigate short-term responses of both pigmented and non-pigmented small eukaryotes  to these simulated climatic conditions by using morphological and molecular methods .
the increases in temperature and uvbr tested in this study correspond to the mean temperature increase expected in the mediterranean region by 2080– <dig>  and the high-uvbr scenario for the european region during spring in future years  <cit> .
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.
methods
experimental design
the four-day experiment  was conducted on the mediterranean platform for marine ecosystems experimental research  located in sète  on the shore of the thau lagoon .
the experimental platform was composed of submerged enclosures  which allowed the isolation of up to  <dig>  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 .
four enclosures, filled with lagoon surface-water at random, were used as incubators for the 2 l experimental bags  in which microbial communities were isolated.
the factorial experimental design constituted eight different treatments : c: control, c + nut: control with nutrient addition, uv: uvbr increase , uv + nut: uvbr increase  and nutrient addition, t: temperature increase , t + nut: temperature increase  and nutrient addition, tuv: temperature  and uvbr  increases, tuv + nut: temperature  and uvbr increases  and nutrient addition .
in order to fill the  <dig> whirl pack bags, 100 l subsurface lagoon water was pumped and pre-filtered through 6-μm-pore-size polycarbonate membranes  in order to isolate the smallest planktonic fraction.
this water sample  was equally distributed into  <dig> sterile whirl pack® polyethylene bags.
<dig> 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  and phosphate in order to maintain a substrate c:n:p molar ratio close to that of marine bacteria  <cit>  as described in bouvy et al.
<cit> .
the bags with and without nutrient addition exhibited concentrations of  <dig>  μm and  <dig>  μm of po <dig>  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   <cit> ).
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  and experimental volume  were chosen to be consistent with the generation time of microorganisms, as validated in other experimental studies  <cit> .
in the present study, the composition of unicellular eukaryotes was studied at t <dig> and t96h.
the data provided by bouvy et al.
<cit>  regarding the evolution of abundances of the main biological communities  at  <dig> sampling times  under the same experimental conditions as ours, informed this choice.
measurement of abiotic parameters
temperature was continuously measured using thermistor probes .
incident uvbr  was constantly monitored by a uvb radiometer .
during the experiment, temperature varied between  <dig> °c and  <dig> °c , while incident uvb radiations , which were measured around local zenith time, varied between  <dig> and  <dig> μwcm- <dig> .
at t <dig> and t96 h, samples were taken for abiotic analysis.
a volume of 80 ml of water was filtered on pre-combusted glass fiber filters  and stored at −20°c until nitrate and phosphate concentrations were measured, following standard nutrient analysis methods  <cit> .
bacterial and viral counting by flow cytometry
at t <dig> 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  equipped with an air-cooled laser providing 15 mw at 488 nm.
for photosynthetic-cells  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  was added, generally for  <dig> to 8 minutes in order to obtain > <dig>  events.
for the analysis of bacteria and viruses, 1 ml fixed  sub-samples were incubated with sybr green i  at a final concentration of 1/ <dig>  for 15 min at room temperature in the dark.
the cytometry flow counts were performed as described in brussard et al.
<cit> .
small eukaryotes microscopy observation
for enumeration of non-pigmented and pigmented eukaryotes, water samples  taken at t <dig> and t96h were fixed with glutaraldehyde  and stored at 4°c for 24 h.  <dig> to 25 ml of each preserved water sample was stained with dapi  for 15 min, filtered onto a black nuclepore filter , stored at −20°c, and counted under an epifluorescence microscope with uv excitation .
under uv light , 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  <dig>  μm pore size polycarbonate membranes  under very low vacuum  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.
<cit>  and extracts were stored at −20°c until analysis.
capillary electrophoresis – single strand conformation polymorphism
nucleic acids from each sample were used as templates for pcr amplification of the 18s rrna gene with primers uni1392r  labelled at the 5’-end with phosphoramidite  <cit>  and euk1209f   <cit> .
each  <dig> μl reaction mixture contained 50 μm of each primer, 1x pfu reaction buffer, 20 mm dntps,  <dig>  u of pfu dna polymerase  and  <dig>  μg of template dna.
pcr amplification was performed with a rob cycler  under the following conditions: an initial denaturation step of 94°c for 2 min, followed by  <dig> touchdown cycles of denaturation at 94°c for 1 min, annealing at 65°c  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  and diluted in sterile te  in order to obtain around 10 ng ml– <dig> of pcr product.
one μl of the dilution was mixed with  <dig>  μl of formamide  and  <dig>  μl of the internal size standard gene-scan- <dig> rox , denatured at 94°c for 5 minutes, and immediately cooled on ice for 10 minutes before electrokinetic injection  into a capillary tube  filled with  <dig> % of gene scan polymer in a abi prism  <dig> genetic analyser .
electrophoresis was carried out and data were collected as described in sauret et al.
<cit> .
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  and ek-1520r  were used for pcr amplification  <cit> .
the pcr mixture  contained about 10 ng of environmental dna, 200 μm of each deoxynucleoside triphosphate, 2 mm mgcl <dig>   <dig> pmol of each primer,  <dig>  u of taq dna polymerase  and the pcr buffer supplied with the enzyme.
reactions were carried out in an automated thermocycler  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  from the triplicate samples of each experimental condition were pooled, precipitated with ethanol–sodium acetate and re-suspended in  <dig> μl of sterile water.
clone libraries were constructed for the t <dig> control and for each of the eight treatments at t96 h using a topo ta cloning kit  with pcr vector  <dig>  according to the manufacturer’s instructions.
phylogenetic analysis
dotur was used to determine operational taxonomic units  from 18s sequences data  <cit>  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   <cit> .
secondly, the otus were aligned with complete sequences in an arb database using the latter’s automatic alignment tool   <cit> .
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 .
the resulting tree was pruned to retain the closest relatives, sequences representative of eukaryotic evolution and our clones .
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 .
the relative distribution of otus in the library was used to calculate coverage values   <cit>  and the non-parametric richness estimator chao <dig>  <cit>  and ace  <cit>  which are the most appropriate indices for microbial clone libraries  <cit> .
statistical analysis
univariate analysis
we tested the homogeneity of the main biological parameters in experimental bags at the initial point  of the experiment using an anova test.
to test the effects of temperature, uv and nutrients on the abundance of all biological groups , we used a three-way anova test .
equality of the variances and normality of the residuals were tested by bartlett and shapiro-wilk tests.
the software sigmastattm  <dig>  was used for all analyses.
multivariate analysis
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- <dig> rox using the safum software  <cit> .
similarity matrices based on bray-curtis distances, dendrograms  and ordination by non-metric multidimensional scaling  were then obtained by using the primer  <dig> software .
one-way analysis of similarity  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  of each sample.
in order to compare eukaryotic communities from the  <dig> genetic libraries unifrac  metrics were used to perform a principal coordinate analysis .
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  and physic-chemical and biological parameters, we used direct multivariate canonical correspondence analysis   <cit> .
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 .
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  at t96 h. eigen values for site scores, biplot and diversity data were plotted to illustrate the associations between these data  <cit> .
RESULTS
initial conditions
biological and chemical parameters
at t <dig>  conditions were considered as homogeneous in all experimental bags.
the statistical analysis showed no significant difference between experimental bags in terms of biological parameters .
means values  are presented for the two sets of experimental microcosms  at t <dig>  for nitrogen and phosphorus compounds, bacteria, viruses, picocyanobacteria, non-pigmented and pigmented small eukaryotes.
* data obtained by flow-cytometry.
abundances and structure of the small eukaryotic community
the microscope counts showed that the eukaryotic community was largely dominated by pigmented cells .
their mean abundance was  <dig>  x <dig> cells ml- <dig> and  <dig> of the  <dig> otus identified at t <dig> from sequencing results were affiliated to pigmented groups .
mamiellophyceae was the dominant group  and they were represented by  <dig> otus affiliated to micromonas pusilla and ostreococcus tauri .
the microscope observations allowed detection of other viridiplantae at low densities.
in particular, some pyramimonadales  were observed but were not recorded among sequences at t <dig>  the mean relative abundance of cryptophyceae  was  <dig> %, while very low relative abundances of bacillariophyceae  and prymnesiophyceae  were found by microscopy  and sequencing.
finally, dinophyceae  accounted for only 3% of total pigmented eukaryotes abundance, and was represented by  <dig> otu .
the mean abundance of non-pigmented eukaryotes was  <dig> cells ml- <dig> at t <dig>  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  <dig> % of total clones .
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,  <dig> different otus could be attributed to non-pigmented groups: cercozoa , stramenopiles affiliated to hyphochytrids , syndiniales affiliated to amoebophrya , uncultured alveolates , and choanoflagellida  .
changes in the small eukaryotes structure under uvbr, temperature and/or nutrient increase
changes in abundances
at t96h, pigmented eukaryotes had abundances varying between  <dig>  x <dig> cells ml- <dig>  and  <dig>  x <dig> cells ml- <dig>  according to the treatment, and they still dominated small eukaryotes regardless of the treatment .
all treatments with increased temperature were characterised by a significant increase in the density of pigmented eukaryotes .
p values obtained for the effects of temperature , uvbr , nutrient addition  and the interactions between the three factors are presented.
+ and – signs indicate the direction of the effect .
bold font corresponds to significant values, where p <  <dig> , while normal font corresponds to a lower significance .
ns is the code for a non-significant effect.
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 .
inversely, for cryptophyceae, a general negative impact of nutrient addition  counteracted the positive impact of temperature increase .
the relative abundance of mamiellophyceae  decreased from t <dig> to t96h in all treatments, and they represented only between  <dig>  and  <dig> % of pigmented eukaryotes at the end of the experiment .
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  tended to increase in abundance in all conditions.
the highest values were found in tuv + nut treatments , however, the 3-way anova did not reveal any significant impact of the manipulated factors .
changes in small eukaryotes structure
a multidimensional scaling  plot generated from bray-curtis similarity indices based on the 18s rdna ce-sscp profiles, showed that all samples from temperature increase simulation  grouped in a separate cluster .
the small eukaryotic community structures of all other treatments  had closer similarity to initial conditions.
overall, ce-sscp profiles generated from all experimental bags showed good reproducibility within triplicate of each treatment , except for one replicate of the uvbr condition which had an atypical profile.
mds ordination plot stress value was low  which indicated good ordination without misleading interpretation  <cit> .
the same trends were found with the upgma  analysis .
changes in small eukaryotes phylogenetic composition
a total of  <dig> otus were identified  .
during the incubation, the richness detected by molecular analyses showed a general decrease in  <dig>  treatments .
tuv + nut was the only treatment characterised by a clear increase in the richness , whereas the greatest decrease was recorded in the c + nut treatment .
even though no general trend was observed in the responses of small eukaryotes in terms of overall richness, the beta-diversity  studied from unifrac metrics revealed a clear association between all treatments with increased temperature .
this highlights the significant structuring impact of increased temperature, while on axis  <dig>  nutrient addition appeared as the second-most important factor in shaping the eukaryotic composition .
these observations were confirmed by analyzing the correlations between coordinates on the pca axis and environmental parameters: coordinates on axis  <dig> were indeed significantly correlated to temperature values  while coordinates on axis  <dig> were significantly correlated to inorganic nutrients concentrations .
the p-values matrix that compares each sample to each other sample showed significant differences in the phylogenetic composition of eukaryotes between t, t + nut, tuv on the one hand and c + nut on the other .
thus, ce-sscp profiles and unifrac analysis led to the same general pattern of changes in the small eukaryote structure.
the otus affiliated to non-pigmented taxa generally dominated the clone libraries .
among them, ciliates and uncultured alveolates were generally well represented .
however, the increase of non-pigmented group proportions within most of the libraries  was mainly linked to the emergence of taxa affiliated to parasitic groups: hyphochytrids and genus pirsonia , and amoebophrya .
the proportion of these sequences clearly increased during the incubation in all types of treatment.
parasitic taxa related to amoebophrya particularly emerged in treatments with the highest temperatures , while hyphochytrids were strongly associated with all other treatments  .
the cca plot illustrates the significant link between the increase in temperature and the presence of numerous sequences affiliated to amoebophrya, while sequences affiliated to hyphochytrides have an opposite position in the plot .
the potential hosts of amoebophrya are primarily found within the class dinophyceae, and it is noticeable that we observed a large number of pigmented dinophyceae cells infected by parasites  at t96 h in all types of treatment .
pigmented dinophyceae were indeed favored by the temperature increase but were also strongly positively affected by nutrient addition and uvbr increase .
pigmented dinophyceae and amoebophrya were represented by  <dig> different otus each.
even though the presence/absence of these otus varied according to the treatments, no association between the abundance of host and parasite otus was observed.
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 .
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 t <dig>  but were not detected at t96 h .
in contrast, pyramimonadales sequences  appeared at t96 h in  <dig> out of the  <dig> types of treatment.
overall, the analysis of the otus dynamics  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 .
the cca showed that  <dig> % of the total variation in the eukaryotic structure was explained by temperature, whereas, uvbr and nutrients explained 11% and  <dig> %, respectively.
discussion
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  <cit> , 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  <cit> .
our goal was to reveal the potential rapid responses of small eukaryote assemblage  during the productive spring season when plankton may be particularly vulnerable to elevated temperature and uvbr  <cit> .
molecular analyses revealed the presence of various phylogenetic groups within the “black box” of small eukaryotes, especially non-pigmented eukaryotes .
some limitations in the pcr-based methods are recognized, for instance, the over-representation of alveolata  in 18s rrna gene clone libraries due their high ssu rrna gene copy number  <cit> .
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  and the comparison of phylogenetic composition obtained for all treatments .
the impact of temperature nutrients and uvbr explained  <dig> %,  <dig> % and  <dig> % of the variance of the small eukaryotes structure respectively.
while bouvy et al.
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  under the influence of temperature, uvbr and nutrient treatments.
among the  <dig> 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  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  effects.
multi-factorial interactions were also apparent for cryptophyceae which experienced antagonistic effects of nutrient addition  and temperature .
in addition to the manipulated factors , 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  in all treatments might be linked to  manipulation effects since these fragile cells might have been affected by filtration steps,  limitation by inorganic nutrients under the rather low orthophosphate concentrations at t96h ,  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  <cit>  and were shown to exert a strong control of bacterioplankton during the study period  <cit> .
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  was observed through sequencing results on taxa affiliated to parasitic and saprotroph groups .
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  <cit> .
the ‘fungi-like’ hyphochytrids possess many morphological and ecological similarities to chytrids  <cit> , and their role as saprotrophs and/or parasites is unclear  <cit> , whereas the amoebophrya are well recognized as a widely distributed parasitic order within the dinophyceae  <cit> .
amoebophrya and hyphochytrids emerged in clone libraries at t96 h and were presumably present among the rare species at t <dig>  the taxa found to be phylogenetically close to amoebophrya particularly emerged in treatments with increased temperature , along with their hosts .
this observation supports guillou et al.’s  <cit>  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 .
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  <cit> .
the responses of saprotrophs to these drivers may result from direct and/or indirect effects as demonstrated in soils  <cit> ; further research is probably needed on the saprotrophs in aquatic systems since changes in their assemblages may influence organic matter decomposition and nutrient cycling.
CONCLUSIONS
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 .
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.
authors’ contributions
all authors have made substantial intellectual contributions to the study.
they read and approved the final manuscript.
tb was the principal investigator of this study.
tb, id, mb, sj, jpt, yb, fv, bm, el, ef participated in the experimental design.
bm, el, tb supervised the operational realisation of the experiment.
id, hm, cb, ef, el realised chemical  and biological analyses , sj performed the flow cytometric analysis.
jfg performed and interpreted the ce-sscp analysis.
cl, id, dd performed the molecular analyses and the post sequencing analysis, ak contributed with cl id and dd to the statistical analysis.
writing was mainly prepared by id, cl, dd and mb, helped by ak, jfg, sj, fv, bm, yb, jpt, tb.
supplementary material
additional file 1
figure s <dig>  maximum parsimony tree showing phylogenetic relationships of the partial 18s rrna gene sequences.
the tree was constructed with the  <dig> sequences generated in this study and sequences from genbank.
only one representative sequence per otu per library is presented in this phylogenetic tree.
the labels show the origin of each sequence .
values in brackets correspond to the otu numbers as presented in figure  <dig> and additional file 2: table s <dig>
click here for file
additional file 2
table s <dig>  composition of the nine 18s rrna genes clone libraries in terms of otus at t <dig> and t96h, the affiliation to phylogenetic groups is specified for each otu.
* the number associated to each otu corresponds to numbers used in figure  <dig> and in the phylogenetic tree .
table s <dig>  unifrac metrics: the grey area  corresponds to the distance matrix obtained from the comparison of each pair of samples.
bold text denotes values in the upper quartile .
the white area  corresponds to the p-values obtained by comparing each sample to each other sample.
all p-values have been corrected for multiple comparisons by multiplying the calculated p-value by the number of comparisons made .
bold text denotes significant p values.
click here for file
acknowledgements
we gratefully acknowledge jean nouguier and yvan vergne for their technical help during the experiment.
this study was supported by the french program pnec  and the anr aquaphage .
this work was also supported by the ‘groupement de recherches   <dig> réseaux trophiques pélagiques.
the experimental platform for mediterranean ecosystem research  was funded by umr  <dig> ecolag, cnrs-inee, institut fédératif de recherche  <dig> a. sabatier, gdr  <dig> réseaux trophiques aquatiques, région languedoc roussillon.
we thank joseph kirkman for improving the text.

 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. 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. the molecular characterization of small eukaryotic assemblages has highlighted an unexpected phylogenetic and functional diversity, 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. however, combined effects of physical factors, such as water temperature and uvb radiation are still poorly investigated. it is recognized that either temperature or uvbr increases can modify microbial dynamics and structure at various levels. 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, 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. 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 and uvbr, at two different nutrients levels. here, we investigate short-term responses of both pigmented and non-pigmented small eukaryotes to these simulated climatic conditions by using morphological and molecular methods. the increases in temperature and uvbr tested in this study correspond to the mean temperature increase expected in the mediterranean region by 2080– 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. methods experimental design the four-day experiment was conducted on the mediterranean platform for marine ecosystems experimental research located in sète on the shore of the thau lagoon. the experimental platform was composed of submerged enclosures which allowed the isolation of up to  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. four enclosures, filled with lagoon surface-water at random, were used as incubators for the 2 l experimental bags in which microbial communities were isolated. the factorial experimental design constituted eight different treatments: c: control, c + nut: control with nutrient addition, uv: uvbr increase, uv + nut: uvbr increase and nutrient addition, t: temperature increase, t + nut: temperature increase and nutrient addition, tuv: temperature and uvbr increases, tuv + nut: temperature and uvbr increases and nutrient addition. in order to fill the whirl pack bags, 100 l subsurface lagoon water was pumped and pre-filtered through 6-μm-pore-size polycarbonate membranes in order to isolate the smallest planktonic fraction. this water sample was equally distributed into sterile whirl pack® polyethylene bags. 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 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  μm and  μm of po 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 ). 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 and experimental volume were chosen to be consistent with the generation time of microorganisms, as validated in other experimental studies. in the present study, the composition of unicellular eukaryotes was studied at t and t96h. the data provided by bouvy et al. regarding the evolution of abundances of the main biological communities at sampling times under the same experimental conditions as ours, informed this choice. measurement of abiotic parameters temperature was continuously measured using thermistor probes. incident uvbr was constantly monitored by a uvb radiometer. during the experiment, temperature varied between °c and °c, while incident uvb radiations, which were measured around local zenith time, varied between and μwcm-. at t and t96 h, samples were taken for abiotic analysis. a volume of 80 ml of water was filtered on pre-combusted glass fiber filters and stored at −20°c until nitrate and phosphate concentrations were measured, following standard nutrient analysis methods. bacterial and viral counting by flow cytometry at t 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 equipped with an air-cooled laser providing 15 mw at 488 nm. for photosynthetic-cells 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 was added, generally for to 8 minutes in order to obtain > events. for the analysis of bacteria and viruses, 1 ml fixed sub-samples were incubated with sybr green i at a final concentration of 1/ 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 taken at t and t96h were fixed with glutaraldehyde and stored at 4°c for 24 h. to 25 ml of each preserved water sample was stained with dapi for 15 min, filtered onto a black nuclepore filter, stored at −20°c, and counted under an epifluorescence microscope with uv excitation. under uv light, 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  μm pore size polycarbonate membranes under very low vacuum 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 nucleic acids from each sample were used as templates for pcr amplification of the 18s rrna gene with primers uni1392r labelled at the 5’-end with phosphoramidite and euk1209f. each μl reaction mixture contained 50 μm of each primer, 1x pfu reaction buffer, 20 mm dntps, u of pfu dna polymerase and  μg of template dna. pcr amplification was performed with a rob cycler under the following conditions: an initial denaturation step of 94°c for 2 min, followed by touchdown cycles of denaturation at 94°c for 1 min, annealing at 65°c 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 and diluted in sterile te in order to obtain around 10 ng ml– of pcr product. one μl of the dilution was mixed with μl of formamide and μl of the internal size standard gene-scan- rox, denatured at 94°c for 5 minutes, and immediately cooled on ice for 10 minutes before electrokinetic injection into a capillary tube filled with % of gene scan polymer in a abi prism genetic analyser. 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 and ek-1520r were used for pcr amplification. the pcr mixture contained about 10 ng of environmental dna, 200 μm of each deoxynucleoside triphosphate, 2 mm mgcl pmol of each primer, u of taq dna polymerase and the pcr buffer supplied with the enzyme. reactions were carried out in an automated thermocycler 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 from the triplicate samples of each experimental condition were pooled, precipitated with ethanol–sodium acetate and re-suspended in μl of sterile water. clone libraries were constructed for the t control and for each of the eight treatments at t96 h using a topo ta cloning kit with pcr vector according to the manufacturer’s instructions. phylogenetic analysis dotur was used to determine operational taxonomic units 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. secondly, the otus were aligned with complete sequences in an arb database using the latter’s automatic alignment tool. 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. the resulting tree was pruned to retain the closest relatives, sequences representative of eukaryotic evolution and our clones. 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. the relative distribution of otus in the library was used to calculate coverage values and the non-parametric richness estimator chao and ace which are the most appropriate indices for microbial clone libraries. statistical analysis univariate analysis we tested the homogeneity of the main biological parameters in experimental bags at the initial point of the experiment using an anova test. to test the effects of temperature, uv and nutrients on the abundance of all biological groups, we used a three-way anova test. equality of the variances and normality of the residuals were tested by bartlett and shapiro-wilk tests. the software sigmastattm was used for all analyses. multivariate analysis 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- rox using the safum software. similarity matrices based on bray-curtis distances, dendrograms and ordination by non-metric multidimensional scaling were then obtained by using the primer software. one-way analysis of similarity 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 of each sample. in order to compare eukaryotic communities from the genetic libraries unifrac metrics were used to perform a principal coordinate analysis. 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 and physic-chemical and biological parameters, we used direct multivariate canonical correspondence analysis. 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. 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 at t96 h. eigen values for site scores, biplot and diversity data were plotted to illustrate the associations between these data.  initial conditions biological and chemical parameters at t conditions were considered as homogeneous in all experimental bags. the statistical analysis showed no significant difference between experimental bags in terms of biological parameters. means values are presented for the two sets of experimental microcosms at t for nitrogen and phosphorus compounds, bacteria, viruses, picocyanobacteria, non-pigmented and pigmented small eukaryotes. * data obtained by flow-cytometry. abundances and structure of the small eukaryotic community the microscope counts showed that the eukaryotic community was largely dominated by pigmented cells. their mean abundance was x cells ml- and of the otus identified at t from sequencing results were affiliated to pigmented groups. mamiellophyceae was the dominant group and they were represented by otus affiliated to micromonas pusilla and ostreococcus tauri. the microscope observations allowed detection of other viridiplantae at low densities. in particular, some pyramimonadales were observed but were not recorded among sequences at t the mean relative abundance of cryptophyceae was %, while very low relative abundances of bacillariophyceae and prymnesiophyceae were found by microscopy and sequencing. finally, dinophyceae accounted for only 3% of total pigmented eukaryotes abundance, and was represented by otu. the mean abundance of non-pigmented eukaryotes was cells ml- at t 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 % of total clones. 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, different otus could be attributed to non-pigmented groups: cercozoa, stramenopiles affiliated to hyphochytrids, syndiniales affiliated to amoebophrya, uncultured alveolates, and choanoflagellida. changes in the small eukaryotes structure under uvbr, temperature and/or nutrient increase changes in abundances at t96h, pigmented eukaryotes had abundances varying between x cells ml- and x cells ml- according to the treatment, and they still dominated small eukaryotes regardless of the treatment. all treatments with increased temperature were characterised by a significant increase in the density of pigmented eukaryotes. p values obtained for the effects of temperature, uvbr, nutrient addition and the interactions between the three factors are presented. + and – signs indicate the direction of the effect. bold font corresponds to significant values, where p < , while normal font corresponds to a lower significance. ns is the code for a non-significant effect. 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. inversely, for cryptophyceae, a general negative impact of nutrient addition counteracted the positive impact of temperature increase. the relative abundance of mamiellophyceae decreased from t to t96h in all treatments, and they represented only between and % of pigmented eukaryotes at the end of the experiment. 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 tended to increase in abundance in all conditions. the highest values were found in tuv + nut treatments, however, the 3-way anova did not reveal any significant impact of the manipulated factors. changes in small eukaryotes structure a multidimensional scaling plot generated from bray-curtis similarity indices based on the 18s rdna ce-sscp profiles, showed that all samples from temperature increase simulation grouped in a separate cluster. the small eukaryotic community structures of all other treatments had closer similarity to initial conditions. overall, ce-sscp profiles generated from all experimental bags showed good reproducibility within triplicate of each treatment, except for one replicate of the uvbr condition which had an atypical profile. mds ordination plot stress value was low which indicated good ordination without misleading interpretation. the same trends were found with the upgma analysis. changes in small eukaryotes phylogenetic composition a total of otus were identified. during the incubation, the richness detected by molecular analyses showed a general decrease in treatments. tuv + nut was the only treatment characterised by a clear increase in the richness, whereas the greatest decrease was recorded in the c + nut treatment. even though no general trend was observed in the responses of small eukaryotes in terms of overall richness, the beta-diversity studied from unifrac metrics revealed a clear association between all treatments with increased temperature. this highlights the significant structuring impact of increased temperature, while on axis nutrient addition appeared as the second-most important factor in shaping the eukaryotic composition. these observations were confirmed by analyzing the correlations between coordinates on the pca axis and environmental parameters: coordinates on axis were indeed significantly correlated to temperature values while coordinates on axis were significantly correlated to inorganic nutrients concentrations. the p-values matrix that compares each sample to each other sample showed significant differences in the phylogenetic composition of eukaryotes between t, t + nut, tuv on the one hand and c + nut on the other. thus, ce-sscp profiles and unifrac analysis led to the same general pattern of changes in the small eukaryote structure. the otus affiliated to non-pigmented taxa generally dominated the clone libraries. among them, ciliates and uncultured alveolates were generally well represented. however, the increase of non-pigmented group proportions within most of the libraries was mainly linked to the emergence of taxa affiliated to parasitic groups: hyphochytrids and genus pirsonia, and amoebophrya. the proportion of these sequences clearly increased during the incubation in all types of treatment. parasitic taxa related to amoebophrya particularly emerged in treatments with the highest temperatures, while hyphochytrids were strongly associated with all other treatments. the cca plot illustrates the significant link between the increase in temperature and the presence of numerous sequences affiliated to amoebophrya, while sequences affiliated to hyphochytrides have an opposite position in the plot. the potential hosts of amoebophrya are primarily found within the class dinophyceae, and it is noticeable that we observed a large number of pigmented dinophyceae cells infected by parasites at t96 h in all types of treatment. pigmented dinophyceae were indeed favored by the temperature increase but were also strongly positively affected by nutrient addition and uvbr increase. pigmented dinophyceae and amoebophrya were represented by different otus each. even though the presence/absence of these otus varied according to the treatments, no association between the abundance of host and parasite otus was observed. 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. 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 t but were not detected at t96 h. in contrast, pyramimonadales sequences appeared at t96 h in out of the types of treatment. overall, the analysis of the otus dynamics 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. the cca showed that % of the total variation in the eukaryotic structure was explained by temperature, whereas, uvbr and nutrients explained 11% and %, respectively. discussion 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. our goal was to reveal the potential rapid responses of small eukaryote assemblage 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. some limitations in the pcr-based methods are recognized, for instance, the over-representation of alveolata in 18s rrna gene clone libraries due their high ssu rrna gene copy number. 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 and the comparison of phylogenetic composition obtained for all treatments. the impact of temperature nutrients and uvbr explained %, % and % of the variance of the small eukaryotes structure respectively. while bouvy et al. 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 under the influence of temperature, uvbr and nutrient treatments. among the 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 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 effects. multi-factorial interactions were also apparent for cryptophyceae which experienced antagonistic effects of nutrient addition and temperature. in addition to the manipulated factors, 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 in all treatments might be linked to manipulation effects since these fragile cells might have been affected by filtration steps, limitation by inorganic nutrients under the rather low orthophosphate concentrations at t96h, 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 was observed through sequencing results on taxa affiliated to parasitic and saprotroph groups. 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. the ‘fungi-like’ hyphochytrids possess many morphological and ecological similarities to chytrids, and their role as saprotrophs and/or parasites is unclear, 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 t the taxa found to be phylogenetically close to amoebophrya particularly emerged in treatments with increased temperature, along with their hosts. 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. 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. 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. authors’ contributions all authors have made substantial intellectual contributions to the study. they read and approved the final manuscript. tb was the principal investigator of this study. tb, id, mb, sj, jpt, yb, fv, bm, el, ef participated in the experimental design. bm, el, tb supervised the operational realisation of the experiment. id, hm, cb, ef, el realised chemical and biological analyses, sj performed the flow cytometric analysis. jfg performed and interpreted the ce-sscp analysis. cl, id, dd performed the molecular analyses and the post sequencing analysis, ak contributed with cl id and dd to the statistical analysis. writing was mainly prepared by id, cl, dd and mb, helped by ak, jfg, sj, fv, bm, yb, jpt, tb. supplementary material additional file 1 figure s maximum parsimony tree showing phylogenetic relationships of the partial 18s rrna gene sequences. the tree was constructed with the sequences generated in this study and sequences from genbank. only one representative sequence per otu per library is presented in this phylogenetic tree. the labels show the origin of each sequence. values in brackets correspond to the otu numbers as presented in figure and additional file 2: table s click here for file additional file 2 table s composition of the nine 18s rrna genes clone libraries in terms of otus at t and t96h, the affiliation to phylogenetic groups is specified for each otu. * the number associated to each otu corresponds to numbers used in figure and in the phylogenetic tree. table s unifrac metrics: the grey area corresponds to the distance matrix obtained from the comparison of each pair of samples. bold text denotes values in the upper quartile. the white area corresponds to the p-values obtained by comparing each sample to each other sample. all p-values have been corrected for multiple comparisons by multiplying the calculated p-value by the number of comparisons made. bold text denotes significant p values. click here for file acknowledgements we gratefully acknowledge jean nouguier and yvan vergne for their technical help during the experiment. this study was supported by the french program pnec and the anr aquaphage. this work was also supported by the ‘groupement de recherches réseaux trophiques pélagiques. the experimental platform for mediterranean ecosystem research was funded by umr ecolag, cnrs-inee, institut fédératif de recherche a. sabatier, gdr réseaux trophiques aquatiques, région languedoc roussillon. we thank joseph kirkman for improving the text.
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. for most of the phylogenetic groups, some rearrangements occurred at the otus level even when their relative proportion did not change significantly. we coupled molecular tools with microscope-based identification and counting to experimentally investigate the short-term responses of small eukaryotes to a warming treatment and uvb radiation increases at two different nutrient levels. temperature explained almost 20% of the total variance of the small eukaryote community structure. interestingly, the increase in temperature resulted in higher pigmented eukaryotes abundances and in community structure changes clearly illustrated by molecular analyses. small eukaryotesmolecular diversitytemperatureuvb radiationmicrocosms experimentmediterranean lagoon 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 and nutrient availability is still poorly investigated.
 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 and nutrient availability is still poorly investigated.  we coupled molecular tools with microscope-based identification and counting to experimentally investigate the short-term responses of small eukaryotes to a warming treatment and uvb radiation increases 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 did not change significantly. temperature explained almost 20% of the total variance of the small eukaryote community structure. 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. small eukaryotesmolecular diversitytemperatureuvb radiationmicrocosms experimentmediterranean lagoon 