Coping with Temperature at the Warm Edge – Patterns of Thermal Adaptation in the Microbial Eukaryote Paramecium caudatum

Background Ectothermic organisms are thought to be severely affected by global warming since their physiological performance is directly dependent on temperature. Latitudinal and temporal variations in mean temperatures force ectotherms to adapt to these complex environmental conditions. Studies investigating current patterns of thermal adaptation among populations of different latitudes allow a prediction of the potential impact of prospective increases in environmental temperatures on their fitness. Methodology/Principal Findings In this study, temperature reaction norms were ascertained among 18 genetically defined, natural clones of the microbial eukaryote Paramecium caudatum. These different clones have been isolated from 12 freshwater habitats along a latitudinal transect in Europe and from 3 tropical habitats (Indonesia). The sensitivity to increasing temperatures was estimated through the analysis of clone specific thermal tolerances and by relating those to current and predicted temperature data of their natural habitats. All investigated European clones seem to be thermal generalists with a broad thermal tolerance and similar optimum temperatures. The weak or missing co-variation of thermal tolerance with latitude does not imply local adaptation to thermal gradients; it rather suggests adaptive phenotypic plasticity among the whole European subpopulation. The tested Indonesian clones appear to be locally adapted to the less variable, tropical temperature regime and show higher tolerance limits, but lower tolerance breadths. Conclusions/Significance Due to the lack of local temperature adaptation within the European subpopulation, P. caudatum genotypes at the most southern edge of their geographic range seem to suffer from the predicted increase in magnitude and frequency of summer heat waves caused by climate change.


Introduction
Temperature is one of the most important environmental factors determining a variety of ecosystem elements, e.g. species ecophysiology, abundance and distribution, as well as species diversity and population dynamics [1][2][3][4]. Due to the current climate change, scientists started to re-evaluate the impact of elevated temperatures on the ecology of species. Here, ectothermic organisms are of special interest as their physiological performance is highly dependent on environmental temperature.
To make predictions of organisms' and population responses to global warming, studies on genetic and phenotypic diversity over a species' geographic range are important. Such investigations can unveil patterns of evolutionary temperature adaptation to the current thermal heterogeneity on Earth by determining which ectotherms have a high acclimatisation capacity and which only occur at specific temperatures. Adaptive phenotypic plasticity, for instance, may cause a higher tolerance to changing thermal conditions [5,6], while local temperature adaptation might be detrimental.
Several studies could show a co-variation of latitude and thermal tolerance (e.g. [7,8]) suggesting that organisms are adapted to the mean temperatures of their environment, but others failed (e.g. [9,10]). Climate change is supposed to affect both climate averages and variability [11] and it has been shown that the thermal tolerance of many organisms is proportional to the magnitude of variation they are exposed to [12]. Organisms are also expected to be adapted to the thermal heterogeneity of their particular environment. This thermal heterogeneity increases with latitude. Therefore, organisms from variable climates, such as the temperate zone, should evolve a broad thermal tolerance resulting in thermal generalists. In contrast, tropical ectotherms, experiencing less variation in temperature, should be selected for narrow thermal niches resulting in thermal specialists [13,14]. Consequently, analysing thermal niches of different populations along a latitudinal transect is necessary to understand the process of adaptation to novel thermal environments.
It has been shown that populations and individuals at the edge of the species range may suffer the most from increasing temperatures, because they often live close to the limit of their species' physiological thermal tolerance [15]. Therefore, it is not only important to investigate the intraspecific variation in species' thermal tolerance, but also to consider populations from the margins of their current distribution range. Especially if one would expect a thermal generalist pattern for ubiquitous species, populations at the 'warm edge' (such as the tropics or subtropics) might be most at risk due to global warming (cf. [16,17]).
Because of the anticipated increasing risk of more intense, more frequent and longer-lasting heat waves during summer [18], species heat resistances are of particular significance [15,19,20,21]. Here, thermal safety margins as well as the maximum warming tolerance are suitable characters to qualitatively elucidate the impact of climate change effects across latitude on different populations. These indicators are based on an organism's thermal tolerance and its relation to the local temperature regime [17]. Studies investigating species' current thermal adaptation patterns with respect to present-day and future environmental temperatures therefore allow predictions of species' and population responses to elevated temperatures.
Beside these patterns of evolutionary temperature adaptation obviously related to climate change, many other patterns are important in thermal adaptation with respect to species evolution and ecology. For example, the warmer is better hypothesis [22,23,24], which predicts a positive correlation between an organism's optimal temperature and its maximum performance; or the Jack-of-all-temperatures is a master of none hypothesis [25], which assumes an evolutionary trade-off between the performance breadth and the maximal performance of an organism, are controversially discussed. These patterns are relevant in a climate change context, too, but only few investigators have experimentally tested these basic ideas of evolutionary temperature adaptation [26][27][28][29].
For the investigation of such elementary hypotheses, the determination of thermal performance curves ( Figure 1) provides a suitable framework to evaluate an organism's thermal tolerance [30]. Thermal performance curves (TPCs) allow estimations on how basic physiological functions are influenced by environmental temperature [17]. Furthermore, TPCs permit the calculation of ecophysiological key characteristics like the lower and upper critical thermal limits (CT min and CT max ) as well as the optimum temperature (T opt ) and the maximum performance ( [31]; cf. Figure 1]). Such key parameters are useful indicators for the thermal tolerance or thermal niche as well as for a potential environmental adaptation of different genotypes. As before mentioned, these ecophysiological characteristics can show a covariation with latitude in metazoan species (e.g. [12,32]), although other studies unveiled that the upper thermal limits of ectotherms vary little with latitude (e.g. [33,34]). However, this has never been critically evaluated for microbial eukaryotes, which are not only important for aquatic ecosystems [35,36], but also constitute well suited organisms for experimental evolution [37,38].
While some recent studies have investigated the response of protozoan species to increasing temperatures (e.g. [39][40][41][42]), little is known about thermal adaptation patterns of globally distributed eukaryotic microbes and how temperature might affect the genetic diversity of natural populations. Furthermore, investigations on the intraspecific variation in species' thermal tolerance by considering populations from the margins of their current distribution range are rare as well.
In the present study, the microbial eukaryote Paramecium caudatum was used to investigate the intraspecific variation in temperature reaction norms of different genotypes. These were isolated from natural habitats along a latitudinal transect in Europe, while three genotypes from tropical habitats (Indonesia) served as a genetic and phenotypic outgroup. This globally distributed ciliate species inhabits the mud-water interface of littoral freshwater environments, which are considerably affected by atmospheric temperature changes [43,44]. Consequently, P. caudatum has to cope with large temporal and spatial variations in temperature. It therefore constitutes a suitable ectotherm to test hypotheses in thermal adaptation as well as the consequences of climate change on such ubiquitous protists. Here, we performed temperature dependent growth experiments to (i) test for a hypothesized local temperature adaptation of different P. caudatum genotypes; (ii) investigate thermal constraints resulting from evolutionary temperature adaptation; and (iii) understand the sensitivity of P. caudatum to predicted future temperatures.

Sampling sites and Organisms
Paramecium caudatum cells were isolated from freshwater samples of 12 different natural habitats along a north-south transect in Europe as well as from three tropical habitats in Indonesia, Sulawesi (see Figure 2 and Table 1 for specifications). No specific permits were required for the described field studies. In Europe and Indonesia, work with Paramecium does not require specific permission and samples were not taken from water bodies where private property was indicated or from nature reserves where sampling is prohibited. The field studies did not involve endangered or protected species.
The food bacteria Enterobacter aerogenes were obtained from the American Type Culture Collection (ATCC 35028) and the kanamycin-resistant strain Pseudomonas fluorescens SBW25 EeZY-6KX [45] was acquired from the University of Oxford. Relationship between environmental temperature and a physiological rate of an ectotherm expressed as a thermal performance curve (grey line). The optimum temperature (T opt ) specifies the temperature at maximum performance. The ecophysiological key characteristics critical thermal minimum (CT min ) and maximum (CT max ) delimit an organism's thermal tolerance. doi:10.1371/journal.pone.0030598.g001 Paramecium Stock MAINTENANCE The investigated P. caudatum stock cultures were maintained in a 0.25% CEROPHYL infusion, prepared according to the methods of Sonneborn [46] with minor modifications [31]. Isolated cells were separated in 1 ml of filtrated habitat water in 24-well tissue culture plates (TPPH AG) to establish clonal cultures. Afterwards, cells were washed and maintained at 22uC in a CEROPHYL infusion inoculated with Enterobacter aerogenes to establish mass cultures. Later, cultures were kept at 10uC, lowering the growth and ageing of P. caudatum. Previous to the start of the experiments, monoxenic P. caudatum cultures were established at 22uC in a CEROPHYL infusion with Pseudomonas fluorescens serving as the only food bacteria (for details see [31]).

Temperature Dependent Growth Experiments
Cells from exponentially growing, monoxenic P. caudatum cultures were transferred to tissue culture flat tubes and acclimatised to experimental temperatures between 7uC and 35.5uC in steps of 61.5 K d 21 . Cultures were kept in exponential growth phase (500-1000 cells ml 21 ) during the acclimation period by doubling the culture volume with Pseudomonas fluorescens inoculated CEROPHYL infusion (CMP; pH 7.0) as appropriate (1-5 ml per day). Due to the different acclimation phases from 22uC up to 35.5uC or down to 7uC, respectively, temperaturedependent experiments were conducted time-delayed. All experiments were performed in microprocessor-controlled, cooled incubators obtained from BINDER GmbH (Type KB 53).
Before the experimental start, acclimatised P. caudatum precultures were adjusted to ,250 cells ml 21 with CMP. Two millilitres of these starting cultures were added to each microcosm containing 18 ml CMP and resulting in an initial abundance of ,25 cells ml 21 . Growth experiments were run in triplicate in 60ml tissue culture flasks with filter lids (TPPH AG) over two to eight days depending on the experimental growth temperature. The bacterial start density was regulated to a saturating prey level of about 2N10 8 cells ml 21 . If the P. caudatum pre-culture densities were below 250 cells ml 21 because of growth-limiting temperatures (e.g. 7uC or $34uC), initial cell abundance was adjusted to the highest possible cell number ($10 cells ml 21 ).
Paramecium cell abundance was estimated by sampling 1 ml every nine to 41 hours depending on the experimental growth temperature. This resulted in five to eight samples per replicate. For precise counting, cells were fixed by the addition of Bouin's solution [47] to a final concentration of 1%. Cell numbers were  enumerated microscopically by threefold counting 100 ml to 300 ml subsamples using a dark field stereoscopic microscope (Olympus GmbH). The population growth rate (m, d 21 ) for each replicate and at each experimental temperature was calculated over the period of exponential increase using the slope of the linear regression of log e -transformed cell densities versus time (t).

Intraspecific Differentiation
To identify and distinguish the individual, clonal cultures of natural P. caudatum genotypes from different geographic regions, the mitochondrial cytochrome c oxidase subunit I (COI) gene was sequenced following the protocol of Barth et al. [48]. Five cells from each stock culture were washed four times in sterile Eau de VolvicHand then incubated overnight with 100 ml of 10% ChelexHsolution and 10 ml Proteinase K (10 mg ml 21 ) at 56uC. Afterwards, the mixture was boiled for 20 min and frozen at 220uC; the supernatant was used for subsequent PCR reactions. Each PCR reaction mix contained 10 ml of ChelexH extracted genomic DNA, 10 pmol of each primer, 1 U Taq-polymerase (SIGMA, Taufkirchen, Germany), 16PCR buffer with 2 mM MgCl 2 and 200 mM dNTPs in a total volume of 50 ml. PCR conditions were as follows: 5 min initial denaturation (95uC); 35 cycles of 1 min at 95uC, 1 min at 50uC and 45 s at 72uC; and a final extension step of 5 min (72uC). Using the primers CoxL11058 and CoxH10176 (see [48]), an 880-bp fragment of the mitochondrial COI gene was amplified. After purification with the Rapid PCR Purification System (Marligen Bioscience, Ijamsville, USA), PCR products were directly sequenced. Sequencing reactions were performed in both directions and analysed on an ABI 3100 Genetic Analyzer (Applied Biosystems).

Thermal Performance Curves
The calculation of thermal performance curves (TPCs) was performed to describe the temperature dependent growth rate data of the individual P. caudatum clones and to determine clone specific key ecophysiological characteristics. TPCs have a common general shape with a gradual increase from a lower critical temperature (CT min ) to a thermal optimum (T opt ) where the investigated biological function reaches its maximum. With a further increase in temperature above T opt the TPCs show a rapid decline towards a critical temperature maximum (CT max ; Figure 1). It was shown that the nonlinear Lactin-2 optimum function [49] can adequately describe the temperature -growth rate relationship of Paramecium caudatum resulting in typically skewed TPCs with a right-shift towards warmer temperatures [31].
The TPC estimation was done by fitting nonlinear mixed-effects models [50] simultaneously to the whole data set. The mixedeffects models were compared with AIC based model selection at three hierarchical levels; the whole data set with common fixed effects for all 18 clones (null models nm0a and nm0b, cf. Table 2), with separate fixed effects for the two regions, Europe and Indonesia (model nm2, cf. Table 2) and with separate fixed effects for the four regions northern, central, southern Europe and Indonesia (model nm4, cf. Table 2). In all cases, all four original parameters of the Lactin-2 function (r, T max , D and l) were used as fixed effects. In model nm0a, all four parameters were also used as random effects while for models nm0b, nm2, and nm4 only T max , D and l were used because of the high correlation between the parameters r and l resulting in a low model convergence. The decision which of the two parameters had to be omitted for nm2 and nm4 was made by comparing the respective AIC values (not shown). Model nm0b is shown for comparison only (cf. Table 2).
Then, the ecophysiological characteristics CT min and CT max were derived numerically as the intersection points of the resulting thermal performance curve with the temperature axis (m = 0). The maximum growth rate (m max, cal ) was calculated analytically as the growth rate (m) at T opt using the Lactin-2 function (Eq.1), while T opt was calculated using its first derivative (Eq.2) as follows: where, the parameter r is a constant influencing m max and the slope of the low-temperature branch, T max is the maximum temperature, and D defines the temperature range of the thermal inhibition above T opt . Parameter l is an intercept parameter that forces the curve to intersect the abscissa at low temperatures and allows the estimation of CT min .
Finally, standard errors for both, the original Lactin-2 function parameters (see Table S1) and the derived ecophysiological key parameters (Table 3) were estimated by nonparametric residual bootstrapping [51] with 1000 bootstrap replicates. For all further analyses based on these key ecophysiological characteristics, estimated data derived from the nm0a mixed-effects model fitting and bootstrapping procedure were used if not otherwise stated.

Climate Data
Investigating an organism's local temperature adaptation or its extinction risk due to climate change requires specific knowledge about the thermal conditions within its natural habitats. Here, we used site-specific temperature data to compare the clonal specific ecophysiological characteristics T opt and CT max with the current The null models nm0a and nm0b were fitted with common fixed effects for all regions, model nm2 with separate fixed effects for the tropical and the European region and nm4 with separate fixed effects for northern, central, southern Europe and the tropical region. For fixed effects the complete set of parameters of the Lactin-2 model (r, T max , D and l; cf. Eq.1) was used in all cases. In model nm0a, all four parameters were also used as random effects, while for models nm0b, nm2, and nm4 only T max , D and l were used. doi:10.1371/journal.pone.0030598.t002 Table 3. Ecophysiological characteristics of individual Paramecium caudatum clones and the two regions, Europe and Indonesia. climate conditions of the specific habitats. Climate data were obtained from nearby meteorological stations (Table S2) or derived from the WorldClim database [52] using the program DIVA-GIS. In case of the meteorological station data, daily mean and maximum air temperature data of the years 2000-2011 (if available) were used to calculate the mean surface air temperature (T hab, mean ) as well as the mean maximal surface air temperature (T hab, max ), both for the warmest three months of the specific habitat. The WorldClim database is a set of interpolated global climate layers considering monthly precipitation as well as mean, minimum and maximum temperatures of the years ,1950-2000.
The database also provides 19 derived bioclimatic variables. The 2.5 arc-minutes resolution database was used to obtain the habitat mean temperature of the warmest quarter (bioclimatic variable 10 o T hab, mean ) and to calculate the habitat mean maximum temperature of the warmest three months (T hab, max ). Furthermore, we used climate change data (2.5 arc-minutes) provided by DIVA-GIS to calculate future conditions for the respective habitats of the investigated P. caudatum clones. These data were derived from high-resolution simulations of global warming [53] using the CCM3 model and assuming a CO 2 doubling until 2100 (see Table S2).
The use of such temperature data has proven controversial and it is well known that local and microhabitat temperature extremes and fluctuations can differ significantly from the regional average [54,55]. However, it could be demonstrated that the summer lake surface water temperature of shallow lakes clearly correlates with the local air temperature [56]. Therefore, the use of local air temperature data seems to be a valid approach to estimate an aquatic ectotherm's performance temperature such as for P. caudatum that inhabits the littoral zone of freshwater environments.

Statistical Analyses
Correlation analyses between each key ecophysiological characteristic (CT min , T opt , CT max ) and the latitude of the respective natural habitats were performed to compare these thermal adaptation indicators with the geographical origin of the different clonal P. caudatum cultures. We used a subset of data containing all investigated European clones and a second dataset including also the Indonesian paramecia. Additionally, the latitudes of the European habitats were corrected for the altitude by assuming that 100 m elevation translates into a ,100 km latitudinal increment within the temperate zone [57]. This correction was done to circumvent altitude effects on the thermal tolerance -latitude dependency. A correction for the Indonesian clones was disclaimed due to nonlinear and extraordinary steep elevational temperature gradients in the tropics [58]. The correlation analyses were performed with absolute values for the altitude-corrected latitude, where Spearman's correlation coefficient and the respective p-values were estimated by Pearson correlation on ranks.
Thermal safety margins (TSM = T opt 2T hab, mean ) as well as the warming tolerances of maximum temperatures (MWT = CT max 2 T hab, max ) were calculated to correlate the clonal specific thermal adaptation indicators and the local climate conditions with the altitude-corrected latitudinal gradient (according to [17]).
A Spearman rank correlation analysis between all calculated maximum growth rates (m max, cal ) and T opt values as well as all m max, cal data and thermal tolerance breadths (TTB = CT max 2 CT min ) was used to test for the so-called warmer is better hypothesis [23] and the Jack-of-all-temperatures is a master of none hypothesis [25], respectively.
Additionally, the coefficients of variation (CV) of the mean growth rate of all investigated clonal P. caudatum cultures were analysed to estimate the intraspecific variation among all tested European clones as well as to assess intra-populational divergence. The calculation for each experimental temperature was as follows: where s m T is the standard deviation of all growth rates (m i ) at the investigated temperatures (T i ) and its arithmetic mean (m T ).
To assess multivariate correlations between the genetic and the ecophysiological as well as geographic distances of the investigated P. caudatum clones, Mantel tests [59] were performed for the whole dataset as well as the European subset. Please note that this does not test for a causal link between the genetic variation of the COI gene and differing thermal tolerances as this gene was chosen to estimate the genetic differentiation within P. caudatum. It is a so called barcoding gene and facilitates comparisons of genetic variation within and among species [48,60].
All statistical analyses were performed using the R system for statistical computing [61] with the add-on package nlme [62] for mixed-effects modelling and package vegan [63,64] for Mantel tests.

Intraspecific Variation in Thermal Performance
Using the Lactin-2 model [49] to describe thermal performance curves (TPCs) of individual Paramecium caudatum clones resulted in typical left-skewed TPCs (Figure 3). These clonal specific TPCs allowed the calculation of ecophysiological key characteristics, which were qualitatively distinguishable between the different P. caudatum clones. For example, the Swedish clone SWV-2A possessed the lowest heat tolerance (CT max = 32.1360.05uC), while the clone from Greece (GRL-1) showed the highest (CT max = 35.1260.06uC) among all investigated European clones. Conversely, the Spanish clone (ESH-2) possessed a comparatively low heat tolerance (CT max = 32.9660.06uC) compared to another Swedish clone (CT max, SWL-1A = 33.8460.07uC). Not only the CT max values showed high differences among the tested European clones, the calculated CT min values were also considerably different (DCT min, EU = 8.4763.13uC, cf. Table 3). On the other hand, all European clones showed their highest growth rates (m max ) at the same experimental temperature of 28uC, while the calculated optimum temperatures (T opt ) of the fitted TPCs ranged from 27.9460.23uC (ESH-2) to 29.4960.15uC (GRL-1).
The comparison between relative differences of intrinsic growth rate data of all European clones showed considerably larger differences with increasing distance from T opt . More precisely, large relative variations (CV) were obtained at low temperatures [CV(m 7uC ) = 21.14%] and especially at temperatures above T opt [CV(m 32.5uC ) = 91.41%], compared to the low relative variation at T opt [CV(m 28uC ) = 7.50%]. While clones from the same habitat (GMA-1A, GMA-1B, GMA-2 and GMA-3) belonging to the same as well as to different COI haplotypes (see Table 1) showed seemingly similar reaction norms, we could detect some variation [CV(m 7uC ) = 14.41%, CV(m 28uC ) = 7.59%, CV(m 32.5uC ) = 17.91%]. Performing one-way ANOVAs of clone specific growth rates at each experimental temperature suggests significant differences for the lowest and the highest temperatures tested (7uC, 32.5uC, 34uC; df = 3, p,0.001, with Bonferroni correction). This result indicates the existence of an intra-populational variation even though we found no significant differences for the intermediate temperatures.
While we could obviously detect only slight differences among the performance curves of different European P. caudatum clones (cf. Figure 3A-C), the Indonesian clones showed considerably different reaction norms compared to the European clones. For these paramecia we could identify significantly higher CT min , T opt , CT max as well as m max values compared to the key characteristics of all European P. caudatum clones (Mann-Whitney U-test; U = 0, p = 0.002, n 1 = 15, n 2 = 3). In addition, an AIC based comparison of nonlinear mixed-effects models with different levels of spatial aggregation (Table 2) identified the model with separate fixed effects for the two regions, Europe and Indonesia (model nm2) as the second best model. While the null model nm0a with common fixed effects for all 18 individual clones was the most parsimonious and significantly superior model, the model with two regions was significantly better than the four regions model (nm4, cf. Table 2). This indicates a clear separation between the European and Indonesian clones, but not within the European paramecia.

Correlation between Ecophysiology and Latitude
As illustrated in Figure 4, we obtained significant negative correlations for CT min (r s = 20.795, p,0.001) and CT max (r s = 20.596, p,0.01) with latitude using the complete dataset. When the correlation analyses were restricted to the European subset, we could only detect a significant correlation for CT min (r s = 20.647; p,0.01; see Figure 4). These results were also supported by Mantel tests for the correlation between geographic and ecophysiological (CT min , T opt , CT max ) distances. Here, analyses on the complete dataset revealed highly significant correlations for all tests (distance matrices of CT min , T opt and CT max vs. geographic distance matrix), while Mantel tests on the European subset showed non-significant relationships at all (cf .  Table S3).
Our results further showed that thermal safety margins increased with latitude as well as altitude ( Figure 5A-C). All low-latitudinal European paramecia from low altitude possessed considerably smaller thermal safety margins than all other clones tested. Interestingly, thermal safety margins derived from nearby weather stations were on average 1.3660.65 K lower than data derived from the WorldClim database (cf. Figure 5A,B). This is potentially due to differently estimated time scales (years 2000-2011 versus ,1950-2000) indicated by the fact that the first decade of the 21 st century apparently was the warmest since climate records began and by holding at least two summers most likely been the warmest in Europe since year 1500 [65,66]. Taking this into account, the lowest observed thermal safety margin was 2.7960.45 K for clone POE-1 from Portugal and the highest for clone NOE-1 from Norway (16.5160.16 K).
As shown in Figure 4C, we could not detect a significant decrease in CT max with increasing latitude by analysing the European dataset. Therefore, the decrease in maximum air temperature with increasing latitude was considerably higher than the decrease in CT max of the respective P. caudatum clones. This resulted in a steep increase in maximum warming tolerance with increasing latitude for the European P. caudatum clones ( Figure 5D-F). Here, the low-latitude clones from Spain and Portugal showed the lowest tolerance window for extreme temperatures ( Figure 5D). Analysing maximum warming tolerances derived from modelled climate change scenarios further revealed that these two lowlatitudinal European clones (ESH-2, POE-1) would show negative maximum warming tolerances ( Figure 5F). The three investigated Indonesian clones, which served as a tropical outgroup, showed higher critical maximum temperatures than all tested European clones ( Figure 4C). Hence, these clones showed positive maximum  warming tolerances even for the predicted future maximum temperatures in their current habitats ( Figure 5D-F).

Thermal Constraints
Analyses of all optimum temperatures (T opt ) and calculated maximum growth rates (m max, cal ) revealed a significant positive correlation of T opt and m max, cal (n = 18, r s = 0.756, p,0.001; Figure 6A). This pattern was not significant among the European subset.
A negative correlation between the maximum performance (m max, cal ) and thermal tolerance breadth was detected using the complete dataset and performing Spearman's rank correlation (n = 18, r s = 20,546, p,0.05; Figure 6B). This indicates that higher maximum growth rates at higher optimum temperatures resulted in a narrower thermal tolerance (Indonesian clones), while a broader thermal tolerance was connected with lower maximum performance (European clones).

Genetic and Ecophysiological Distances
Performing Mantel tests [59] on the complete dataset revealed highly significant correlations between ecophysiological (CT min , T opt , CT max ) and genetic distances. Excluding the Indonesian clones from the analyses resulted in non-significant relationships (see Table 4). This indicates a missing correlation of the ecophysiology and the genetic distance within the investigated European P. caudatum clones using the mitochondrial COI gene as a phylogenetic marker and temperature dependent population growth rates as a fitness component. However, please note that the Indonesian clones possess both high geographical distances (10,715 km-13,261 km) as well as high genetic distances (0.073-0.087 substitutions per site) compared to the European P. caudatum clones, which exhibit comparatively low genetic distances among each other (0-0.02 substitutions per site). Performing Mantel tests for the genetic and geographic distances resulted in significant correlations for both datasets (Table S4).

Latitude-dependent Ecophysiology
Generally, high-latitude Paramecium caudatum populations encounter lower temperatures than low-latitude populations, which should select for higher growth rates at low or high temperatures, respectively. Northern European P. caudatum clones should therefore exhibit lower critical minimum temperatures (CT min ) to which the southern European clones were barely exposed to, while southern European clones should possess higher critical maximum temperatures (CT max ). Our data partly support this general assumption of a co-variation of the critical thermal limits (CT min , CT max ) and the optimum temperature (T opt ) with the latitude when analysing the European dataset. Performing correlation analyses on this dataset revealed a significant dependence of CT min with altitude-corrected latitudes, although we partly obtained high standard errors based on the used fitting procedure (cf. Figure 4A, Table 3). In addition, we could not detect such a correlation for T opt and CT max using the European dataset ( Figure 4B,C). Our data therefore indicate a potential thermal adaptation of the northern European P. caudatum clones to lower winter temperatures only. These results are in accordance with a number of investigations on terrestrial ectotherms, which revealed that the lower critical temperatures significantly decline with increasing latitude while the upper thermal limits do not or are less variable [12,33,34,67,68]. This might be due to different costs for cold and warm tolerance, but up to now there is no clear evidence whether cold or warm adaptation came at a higher cost or if high heat tolerance is the ancestral state with low retaining costs (cf. [33,69,70,71]).
When using the complete dataset that includes experimental data of all tested European as well as the three Indonesian clones, significant co-variations of CT min as well as CT max with latitude were obtained ( Figure 4A,C). Additionally, the results of the model selection approach (Table 2) identified the model nm2 with two regions (Europe and Indonesia) as significantly superior compared to the four regions model nm4 (with regions for northern, central, southern Europe and Indonesia). These results imply a stronger impact of latitude on the thermal performance of the eukaryotic microbe P. caudatum when comparing data on a large intercontinental geographic scale, in contrast to the analyses on the intra-continental European scale. Further, significant positive correlations of all ecophysiological distances with genetic distances could be identified only when using the complete dataset, but not among the European subset (Table 4). Hence, the detected high genetic distances and significantly different ecophysiological characters of the Indonesian compared to the European clones as well as the results of the model selection approach suggest a large-scale biogeographic diversification within Paramecium caudatum on the phenotypic as well as the genetic level.

Phenotypic Plasticity AND Thermal Adaptation
In general, the thermal performance of all investigated European Paramecium caudatum clones is indicative for a high phenotypic plasticity of this freshwater ciliate. Exemplified by the northern European clones, they showed a higher physiological optimum (around 29uC) compared to the temperatures they experienced in their natural habitats ( Figure 5A). The optimum temperatures and the shapes of the thermal performance curves (TPCs) were reasonably similar for all genetically distinct clones from Europe (cf. Figure 3, 4B). Further, all European clones showed a general broad thermal tolerance and non-covarying CT max as well as T opt values with latitude ( Figure 4B,C). These results disapprove the hypothesised latitudinal clines for thermal adaptation indicators such as heat tolerance and optimum temperature in European P. caudatum genotypes, but are in Figure 5. Latitudinal trends in thermal safety margin and maximum warming tolerance. Thermal safety margins (A-C) and maximum warming tolerances (D-F) of the investigated Paramecium caudatum clones were calculated using habitat temperatures from global climate layers (A+D), near-by meteorological station data (B+E) and from climate change projections (C+F). Symbols represent the mean 6 standard error of the mean. Latitudes of the European clones were corrected for altitude assuming a 100 km increase in latitude for a 100 m increase in altitude. Arrows indicate P. caudatum clones from high altitude (cf. agreement with the climatic variability or the seasonal variability hypotheses [13,72,73]. These hypotheses claim that greater environmental variability at higher latitudes, for example due to seasonal changes, select for a more 'generalist' climatic tolerance and favours phenotypic plasticity within populations. Such a scenario would be also supported by the low genetic differentiation among the investigated European P. caudatum clones. This fact could be interpreted as an increased gene flow facilitated by the high phenotypic plasticity and the resulting low dispersal costs [13]. In that case, the high gene flow and high dispersal rates would also limit a potential local adaptation to specific habitat temperatures. As our data revealed a comparatively broad thermal tolerance for all European clones with exceptional high T opt and CT max values, we would argue for an adaptive phenotypic plasticity among the investigated European populations resulting in thermal generalists. On the other hand, arguments for a potential thermal adaptation of the European P. caudatum to local microhabitat conditions are obvious due to the co-varying CT min values with latitude ( Figure 4A), but also because of large intraspecific growth rate variations at temperatures above T opt (Figure 3). These differences at low and high temperatures indicate the existence of various ecotypes and a potential microadaptation to the local microclimate within the European P. caudatum clones. Temperature adaptation to microhabitat conditions has also been shown by several studies for ectothermic metazoans [74][75][76][77].
Along the three tested Indonesian clones, it was remarkable that in comparison to the European clones the averaged lower critical temperature of these tropical paramecia was shifted by more than +6 K, whereas T opt and CT max were only shifted by approx. +1.5 or +2.1 K, respectively. This means a reduction of their thermal tolerance breadth compared to the European P. caudatum clones. Given that in general the Indonesian clones were hardly ever stressed by temperatures below 10uC in their natural habitats, these tropical populations have not needed to adapt to lower temperatures. They have either lost the low-temperature tolerance or have never had this ability, depending on whether the tropical or the European clones represent the most ancestral phenotype. On the other hand, the Indonesian populations from low altitudes experience higher frequencies of hotter daily maximum temperatures. For example, the maximum air temperature for the natural habitat of clone INP-3 could reach temperatures of up to 43uC. However, the mean maximum air temperature, which corresponds arguably better to the maximum water temperature because of the buffering capacity of water [78], is around 33uC. Here, we could show that the CT max values of all investigated Indonesian P. caudatum clones (cf. Table 3) were somewhat beyond this temperature regardless of the elevation of their natural habitats. This result is indicative for a common thermal adaptation of the investigated Indonesian P. caudatum clones to the tropical temperature regime. Figure 6. Thermal constraints. A) Relationship between calculated optimum temperature (T opt ) and maximum growth rate (m max, cal ) for all investigated P. caudatum clones supporting warmer is better. Symbols represent the mean 6 standard error of the mean derived from the nonlinear mixed-effects model nm0a with residual bootstrapping. Significance was tested with Spearman's rank correlation (n = 18, r s = 0.775, p,0.001). B) Trade-off between calculated maximum growth rate (m max, cal ) and thermal tolerance breadth (TTB = CT max 2CT min ). Symbols represent the mean 6 standard error of the mean derived from the nonlinear mixed-effects model nm0a with residual bootstrapping. Significance was tested by using Spearman's rank correlation (n = 18, r s = 20,554, p,0.05). doi:10.1371/journal.pone.0030598.g006 These findings support the general assumption that acclimatization and a high phenotypic plasticity is more likely in temperate species or populations at higher latitudes because the overall temperature variation increases with latitude. Thus, the evolution of broad thermal tolerances is needed that temperate organisms can cope with a large seasonal variation while organisms from less variable tropical environments should have evolved narrow thermal tolerances and reduced acclimation responses ( [13,32,33]; but see [79]), as shown by our study. To generalise this pattern for P. caudatum, more clones from other tropical habitats need to be investigated.

Thermal Constraints
In this study, a significant correlation between optimum temperature and maximum population growth rates could be shown ( Figure 6A), which supports the so-called warmer is better hypothesis [22,23,24]. This thermodynamic-constraint hypothesis argues for a dependence of maximum performance on optimum temperature because of the thermodynamic properties of biochemical and physiological systems [24,80]. Our results document that P. caudatum clones with higher optimum temperatures (T opt ) have generally higher maximum growth rates (m max ) giving evidence for warmer is better concerning the growth performance as an important component of overall fitness. Nevertheless, we could not observe a significant correlation of T opt and m max for the European subset, which indicates that the scale at which investigations are performed is of importance. Further studies on geographically well described organisms are necessary to understand if this finding is of general relevance for a variety of ectotherm species, for microbial eukaryotes only, or just for Paramecium caudatum.
The present study also provides an indication of an evolutionary trade-off between the performance breadth and maximum performance, which is known as the Jack-of-all-temperatures is a master of none hypothesis [25]. This result indicates that the selection for a broad thermal tolerance could result in a lower peak performance [81] while the selection for greater performance at a higher temperature would cause a correlated decrease in performance at lower temperatures [20]. Several studies provided mixed or no support of such specialist-generalist tradeoffs that constrain TPCs (e.g. [82][83][84][85]). Here, we could demonstrate such a negative correlation between maximum growth rates (m max ) and thermal tolerance breadths ( Figure 6B). Especially the three tested Indonesian paramecia possessed higher m max values but narrower thermal tolerances than the European clones, mainly caused by a considerably larger shift of CT min to higher temperatures compared to CT max . This result supports the above formulated suggestion that the European P. caudatum clones could be indicated as thermal generalists with a maximised performance breadth while the Indonesian clones seem to be thermal specialists with a maximised peak performance at higher temperatures.

Sensitivity to Predicted Increasing Temperatures
The impact of the ongoing climate change on an ectotherm's fitness depends on numerous factors, including the community responses or the resource availability and the temperature-specific resource demand [86]. However, due to the expected temperature rise of up to 5.8uC by the year 2100 and the anticipated increase in diurnal variability of summer temperatures in the Northern Hemisphere [18], temperature is one of the most important factors which can drive shifts within the structure of natural populations [1,2,3]. Therefore, climate change effects are often predicted to give rise to species extinctions [82,87]. The consequence of this temperature effect will considerably depend on the genotype specific T opt and CT max values relative to the mean and extreme habitat temperatures [88]. Thermal safety margins as well as maximum warming tolerances are suitable characters to reveal the potential extinction risk of organisms due to increasing temperatures [17].
In the present study, positive thermal safety margins reflect the temperature range at which the respective P. caudatum clone may benefit from a future warming of the habitat due to increasing growth rates. Negative values are rather a measure of potential risk. Our analysis revealed considerably smaller thermal safety margins for the low-latitudinal and -altitudinal European and tropical paramecia compared to all other P. caudatum clones tested ( Figure 5A). However, all investigated clones currently live in environments that are on average cooler than their physiological thermal optimum. All tested clones seem to potentially benefit from future increasing temperatures ( Figure 5C) at least initially.
Nevertheless, a second important key characteristic especially in consideration of the expected increasing intensity, frequency and duration of summer heat waves is the maximum warming tolerance. This character illustrates the average increase in maximum temperatures which Paramecium can tolerate before harmful growth conditions will be reached. Unexpectedly, analyses of the maximum warming tolerance showed that the European paramecia from low latitude and altitude currently experience near-or even above-lethal temperatures during summer ( Figure 5E). Furthermore, the predicted future maximum temperatures of the low-latitudinal habitats (Portugal and Spain) are higher than their critical maximum temperatures ( Figure 5F). Climate change models further predict the highest warming rates for low-latitude European habitats based on our dataset. Consequently, if these genotypes cannot adapt to the expected higher temperatures in their current natural habitats, they will potentially suffer from global warming. This seems not to be the case for the Indonesian and all high-latitudinal and -altitudinal European clones, which possess higher maximum warming tolerances ( Figure 5D-F).
Our results, therefore, only partly support the hypothesis that tropical ectotherms are most at risk due to novel as well as disappearing climates in the tropics and subtropics in consequence of climate change [17,89,90]. In terms of the microbial eukaryote Paramecium caudatum, only the most southern European populations seem to be adversely affected by global warming. While some of the high-latitude European populations may actually benefit from increased temperatures by an enhanced population growth, the investigated Indonesian clones seem also not to suffer from the expected temperature rise. This is due to their adaptation to higher temperatures as well as the fact that temperature increase in the tropics is expected to be less intensive compared to temperate habitats [18]. However, tropical P. caudatum populations may become also affected by global warming, because the predicted changing environmental conditions such as temperature and precipitation seem to be comparatively heterogeneous across latitude [90,91] and we have only a small dataset of tropical paramecia. Further, how climate change affects species abundances, distribution or diversity depends on the multigenerational response of their survival and reproduction within ecosystems [92]. The magnitude of temperature effects on species hinge on several factors such as food-web interactions, altered competition as well as the species specific acclimation or adaptation capacity [88,93,94]. Hence, experimental selection and micro-evolutionary studies as well as competition experiments of artificial populations pose an interesting research outline for future studies.

Conclusions
All investigated European clones showed a broad thermal tolerance with high upper thermal limits as well as similar optimum temperatures, which is indicative for thermal generalists. These genetically closely related P. caudatum clones may have been selected to perform well in thermally variable environments utilising a high acclimation capacity. In contrast, the investigated Indonesian genotypes showed significantly higher optimum temperatures, maximum performances and critical thermal limits. This suggests local temperature adaptation of these tested clones to the less thermally variable tropical temperature regime. Further studies on additional 'continental subpopulations' of P. caudatum are needed to generalise the suggested thermal generalist and thermal specialist patterns in this ubiquitous microbial eukaryote.
Relating ecophysiological key characteristics such as optimum and maximum temperatures to current and predicted temperatures of all investigated P. caudatum clones suggest that only lowlatitude European clones would be sensitive to global warming. During summer, they currently perform closer to their thermal limits than high-latitude European and tropical clones. Increasing temperature extremes, therefore, may severely affect their performance and fitness. Future investigations on experimental evolution would be useful to examine whether the different genotypes tested in this study can adapt to the predicted increase in temperatures and extreme events, or not.