Experimental simulation of environmental warming selects against pigmented morphs of land snails

Abstract In terrestrial snails, thermal selection acts on shell coloration. However, the biological relevance of small differences in the intensity of shell pigmentation and the associated thermodynamic, physiological, and evolutionary consequences for snail diversity within the course of environmental warming are still insufficiently understood. To relate temperature‐driven internal heating, protein and membrane integrity impairment, escape behavior, place of residence selection, water loss, and mortality, we used experimentally warmed open‐top chambers and field observations with a total of >11,000 naturally or experimentally colored individuals of the highly polymorphic species Theba pisana (O.F. MÜller, 1774). We show that solar radiation in their natural Mediterranean habitat in Southern France poses intensifying thermal stress on increasingly pigmented snails that cannot be compensated for by behavioral responses. Individuals of all morphs acted neither jointly nor actively competed in climbing behavior, but acted similarly regardless of neighbor pigmentation intensity. Consequently, dark morphs progressively suffered from high internal temperatures, oxidative stress, and a breakdown of the chaperone system. Concomitant with increasing water loss, mortality increased with more intense pigmentation under simulated global warming conditions. In parallel with an increase in mean ambient temperature of 1.34°C over the past 30 years, the mortality rate of pigmented individuals in the field is, currently, about 50% higher than that of white morphs. A further increase of 1.12°C, as experimentally simulated in our study, would elevate this rate by another 26%. For 34 T. pisana populations from locations that are up to 2.7°C warmer than our experimental site, we show that both the frequency of pigmented morphs and overall pigmentation intensity decrease with an increase in average summer temperatures. We therefore predict a continuing strong decline in the frequency of pigmented morphs and a decrease in overall pigmentation intensity with ongoing global change in areas with strong solar radiation.


| INTRODUC TI ON
To better understand and predict the impact of climate warming on organisms requires multitrait approaches (Debecker & Stoks, 2019).
Recognizing the potential impacts of global warming, the physiological mechanisms that limit heat tolerance in ectotherms regained interest around the turn of the millennium (Pörtner, 2002). More recently, the entirety of phenotypic traits including behavior, physiology, and life cycle parameters that correlate in concert with environmental variation, has been identified and studied as syndromes (Boyle et al., 2015;Killen et al., 2013), helping to gain further understanding into the mechanisms behind the selection pressures generated by climate change (Brook et al., 2008;Stevens et al., 2013).
Among the best-studied examples for phenotypic variation in response to climate and habitat are shell pigmentation polymorphisms of helicoid land snails, which have been investigated in detail for more than 50 years (reviewed in Ożgo, 2014;Schilthuizen, 2012;Schweizer et al., 2019). Numerous studies on these animals have addressed single morphological features, physiological traits, or biochemical parameters that are modified by thermal stress. However, the concurrent syndromic responses to thermal stress have rarely been addressed, but are of crucial importance for a mechanistic understanding of stress responses across biological levels of increasing complexity (Schweizer et al., 2019).
In addition, two resurvey studies in England and the Netherlands revealed an increase in individuals with high-albedo shells at the cost of darker ones after 43 years in the garden snail, Cepaea nemoralis (Cameron et al., 2013;, which was interpreted to be adaptive to environmental temperature increase. Stine (1989) reported the same effect for C. nemoralis after several decades for Virginia, USA, Johnson (2011) for Theba pisana in Western Australia. Another large study that specifically addressed the effects of climate change-induced temperature increases in Europe revealed an increase in yellow individuals in sand dunes, the most exposed habitat type, but did not find general support for selective advantages of high shell albedo across habitats in C. nemoralis (Silvertown et al., 2011). Although land surface temperatures in Europe have increased by an average 1.3°C during the 20th century until 2009 (European Environment Agency, 2018), data obtained in this study made evident that garden snails did not show any increase in the frequency of the lightest (yellow) morphs over time but rather an unexpected decrease in the frequency of unbanded shells (independent of the basic coloration, which can be yellow, pink, or brown) | 1113 KÖHLER Et aL. and an increase in the "mid-banded" morph. Although precise shell pigmentation intensity had not been quantified, Silvertown et al. (2011) and Cameron and Cook (2012) concluded that banding was probably not the most important trait that selection had acted upon. They suggested that factors other than climate change-induced temperature increases, such as changing predation pressure or habitat changes affecting microclimatic conditions, might have asserted stronger selection pressure. Indeed, many other studies have suggested alternative explanations to thermal selection for variation in the banding frequency of land snails. For example, the selection on visual crypsis through predation by birds or rodents has often been hypothesized as a key component limiting morph frequencies (Allen, 2004;Bond, 2007;Cain & Currey, 1963;Cain & Sheppard, 1954;Clarke, 1969;Cook, 2005;Endler, 1978;Heller & Gadot, 1984;Moreno-Rueda, 2009;Punzalan et al., 2005;Rosin et al., 2011). Furthermore, while a recent study on C. nemoralis from urban "heat island" areas in the Netherlands found a higher frequency of the palest morph (yellow), these yellow snails were also more likely to carry bands on the umbilical side of the shell. As a result, the authors proposed a possible role of such pigmentation in a yet unknown way (Kerstes et al., 2019), adding further complexity to the theory of selection due to thermal stress.
The manifold results and conclusions of these studies suggest that the role of thermal selection in the evolutionary ecology of land snails in response to global change, and its relevance for differently pigmented morphs, remains only partly understood. This is especially true in snails with extremely variable shell coloration such as the Mediterranean species T. pisana. In this species, shell color and banding patterns are primarily determined by alleles of at least three loci (Cain, 1984;Cowie, 1984), but the interplay of genetic variation and phenotypic plasticity results in a situation in which almost every individual displays a unique pigmentation pattern and intensity (Cain, 1984;Cowie, 1984;Johnson, 2012;Köhler et al., 2013). Furthermore, physical integrity and survival over the hot season are of particular evolutionary relevance in T. pisana because this species reproduces at the end of its lifetime. In this context, combining field observations with common garden experiments, the analysis of trait syndromes can shed light on the causal effects of thermal selection.
The Mediterranean basin is globally among the regions which are most vulnerable to climate change, with annual mean surface air temperature increasing from 0.19 to 0.25°C per decade between 1960 and 2005 (Mariotti et al., 2015) and 1.1°C from 1985 to 2015 (Cramer et al., 2018). Depending on the climate scenario and the season, a rise in temperature versus pre-industrial time from 2 to 6°C by 2100 is expected in the Mediterranean (Giorgi & Lionelli, 2008;Jacob et al., 2014). Using T. pisana, we here experimentally simulated warming in the expected range under global climate change and artificially altered shell pigmentation in a Mediterranean field plot.
In a multitrait approach, we combined our experiments with field measurements of thermodynamics, behavioral studies on snail competition, biochemical markers for protein and lipid destruction, and measurements of transpiration and survival. The goal of this study thus was to identify a mechanistically plausible interplay between thermodynamics, biochemistry, physiology, and anticipated population changes in the course of global warming at a given site in one of the prospectively most threatened regions on earth. Connecting field measurements with passive open-top chambers (OTCs) to manipulate air and soil surface temperatures Marion et al., 1997;Welshofer et al., 2018;Yan et al., 2015), we aimed to test the following hypotheses: • [H1] Even faint differences in shell pigmentation intensity result in different warming of morphs.
• [H2] Simulation of global warming causes increasing protein degradation and oxidative stress in increasingly pigmented snails.
• [H3] As T. pisana often forms clusters at resting sites, morphs of different pigmentation intensity interact (i.e., actively compete or act jointly) in their behavior when escaping the hot soil surface in a morph-specific way.
• [H4] Water loss due to environmental warming differs among different morphs.
• [H5] Survival rates of morphs decrease with increasing pigmentation intensity under simulated global warming conditions. The use of OTCs enabled us to change the expression of fluctuations in environmental temperature rather than just to increase a constant mean temperature. This approach thus allowed close to nature acclimation independent of the basal thermotolerance (Gerken et al., 2015) and avoided shifts in the animals' performance optima that may occur under constant conditions (Dowd et al., 2015). Using our OTC-derived data, we also aimed to predict future morph frequencies and average shell pigmentation intensity in our experimental field site with ongoing global warming. We tested our predictions by field observations on morph frequencies in other Mediterranean sites that already are subject to such thermal conditions that are expected for the experimental site.  (Figure 1). The site is in the vicinity of Avignon climate station from which we obtained climate data (global solar irradiation, air temperature in the shade at 1 m above the ground, and precipitation) for the time of our experiments (9 August-6 September 2017; temporal resolution 15 min). Furthermore, we obtained long-term temperature data (monthly maximum, monthly mean) for this site over a 30-year period (January 1988-December 2017, resolution 1 month).

| Field site and population abundance
Individuals of T. pisana from the field site were classified according to the pigmentation of their shells in two categories, "naturally white" (W) and "naturally pigmented" (P), the latter comprising individuals with one or more (even faint or dotted) bands. To quantify the frequency of both morphs at the field site, 24 samples of randomly collected snails with a total number of 3,339 individuals were taken from the vegetation and the snails allocated to the categories W and P. The scarcity of vegetation at this site prevented potential camouflage of snails and thus bias in sampling, as all individuals were easily visible. We refrained from sampling empty shells from the ground to circumvent uncertainty in respect to postmortem shell bleaching. To quantify the mortality of the different morphs in the field under natural conditions, we randomly collected 2,160 individuals that were attached to the vegetation by dried mucus surrounding the aperture, allocated them to the categories W and P, and checked whether individuals were dead or alive.

| Manipulation of snail pigmentation
The purpose of shell color manipulation was to provide a series of morphs with gradually increasing pigmentation intensity from white to black. As most individuals at the Montfavet site were naturally white, we exclusively applied black ink on them and refrained from painting pigmented ones white. All experiments were thus conducted with four morphs of T. pisana. Apart from the two natural morph classes ("naturally white," W; and "naturally pigmented," P), we experimentally colored white individuals either by adding a rather thick band to the apical side of the shell ("artificially striped," S) or by painting both sides of the shell entirely black ("artificially blackened," B) using a carbon-based tattoo ink on water basis without solvents and free from heavy metals, phenol, and polycyclic aromatic hydrocarbons (Deep Colours! Black, DC-TP Europe LLC). Using scanning electron microscopy (SEM), we found that the artificial black color only added a very thin layer of pigment to the exterior of the shell. The fine structure of the shell's surface was smoothened in the painted areas. This smoothening, however, rather resembles the natural surface of pigmented shell parts (bands) in T. pisana, as bands exhibit a smoother surface than pale shell parts in SEM pictures (Schweizer et al., 2019).

| Escape behavior and snail temperature
For the studies on escape behavior, we alternately placed five indi- and (c) reach 20 cm height to the wall's top (rank 1 quickest, rank 10 slowest). Rank sums were calculated for each morph and each behavioral criterion (a-c). Subsequently, for each morph pair, a Δ rank sum was calculated (the rank sum of the darker morph minus the rank sum of the lighter morph), which was related to the difference in pigmentation intensity calculated for the corresponding pair of snails (for pigmentation measurements see below). We recorded randomly selected runs with a high-resolution thermocamera (Micro-Epsilon Thermo Imager 400, equipped with a macrolens) and measured shell and soft body temperature of the snails on the ground (n = 31 for each morph) and at heights of 6.5 cm (n = 18-25), 10 cm (n = 17-21), and 20 cm (n = 14-17) above ground level for every individual possible.

| Morph-specific body and shell temperatures
Individuals from all four morphs were kept in the laboratory at 22°C in a box, which was transported outside right before measuring body and shell temperatures. Two individuals from different morphs were randomly placed within 30 s on two needle thermometers that had been calibrated against one another at a height of 50 cm above ground in the full sunlight (90-100 klux). After piercing the body at the aperture, next to the parietal wall with the needle, the internal temperatures of the two individuals were measured at 30-to 60-s intervals in the sun. At the same time, the shell temperature (apical side) was measured using high-resolution thermography. We calculated ΔT (T of the darker morph minus T of the lighter morph) for both the internal and shell temperatures, and for each pair of snails.

F I G U R E 1 Theba pisana in its natural habitat at Montfavet, Southern France
In total, 44 individuals were measured. All snails survived the piercing of their skin.

| Open-top chamber experiments on climbing ability, water loss, and mortality
Passive open-top chambers were constructed for the purpose of locally raising the temperature in the snails' habitat, simulating the conditions of environmental warming (Marion et al., 1997;Welshofer et al., 2018). showing that this experiment simulated environmental warming of 1.12 ± 0.94°C above ambient temperatures. Control snails were kept in the laboratory at a constant temperature of 25°C and a natural day/night regime.
Climbing behavior and mortality were recorded in each OTC on days 1, 2, 3, 5, 7, 9, 12, 14, 16, 19, 21, 23, and 26. For climbing, we counted the number of morphs that had climbed ≥10 cm on a skewer, as well as the number and identity of the morph that had climbed highest, but at least 10 cm ("top" position). After counting, all individuals were placed on the ground again to avoid inactivity (and thus double counting) at an elevated position for more than one day. Mortality (i.e., the number of dead snails relative to the starting number) was recorded on the same days for each OTC and each morph. All snails were weighed at day 0 (prior to the introduction into the OTCs) and at day 8. As it was impossible to mark the snails individually, we concentrated on the average fresh weight of surviving snails for every morph and OTC.
Climbing, mortality, and weight measurements were statistically analyzed using generalized liner mixed models in R version 3.4.2 (R Development Core Team, 2014) using package "lme4" (Bates et al., 2015). Statistical significance was calculated by model simplification, with multicomparison tests and estimates performed using least-square means (emmeans package, Lenth, 2018). Intercepts of OTC identity were used as random variables in all models.
Climbing had three response variables: (a) "number of a morph higher than 10 cm" per OTC-modeled as a either a Poisson distribution (log-link function) with the number of individuals per morph that climbed higher than 10 cm, or modeled as (b) a binomial distri- We also used two response variables for mortality: "Death rate"-modeled as a binomial distribution (logit link function), with the number of living individuals out of the total starting number (20)-and "time to D50"-modeled as a gamma distribution (log-link function) as the time (days) it took for 50% (10 individuals per morph) to die (note that if 10 deaths were reached between observational recordings, fractional days were calculated on a linear scale).
For both climbing and mortality, models were constructed to investigate OTCs that only contained the individual pairwise comparisons of morph type (i.e., B vs. W or B vs. S) "single-competition" OTCs. Models were also constructed to investigate overall comparisons of morph type acting as either the "competitive response" of an individual morph (the value of the response variable being dependent on the morph of its neighbor) or "competitive effect" of an individual morph (the value of the neighboring morph in response to any impact a morph may have). Tests were made for each morph to test for differences in the effect or response with all three alternative morph combinations. General patterns were also calculated for competitive response and effect using like-for-like comparisons, for example, in a comparison of the competitive responses of B versus W morphs, we also included B versus S and B versus P comparisons to W versus S and W versus P, in order to remove bias from W versus B OTCs and B versus W OTCs.
To test for differences in morph response and effect on climbing, solely early observational time points at 2, 3, and 5 days were used with a categorical explanatory variable "time" included in the model, to exclude impact of death at later time points from the model.
Death rate also included time as a continuous variable from points 1 to 9, as later time points included too many zeros for a reliable analysis. Here, the slope and therefore interaction of morph identity x time were the important factor of analysis. Finally, weight at time 0 and time point 8 were analyzed separately to determine differences in weight dependent on morph type. Weights were log-transformed

| OTC experiments on snail temperature, oxidative stress, and proteotoxicity
Sixteen open-top transparent polyethylene chambers (OTCs) with a ground area of 30 × 20 cm and a height of 20 cm were arranged in a 4 × 4 block design on the field ground mentioned before (INRA Avignon territory, 3 blocks) and in the laboratory (control at 25°C, 1 block). The OTCs were equipped as described above. In each of the four blocks, the four OTCs housed a single morph (W, P, S, or B) with 20 individuals. Animals were exposed to ambient conditions Supernatants were stored at -80°C until further analysis. The assay was conducted in 96-well plates. In each well (except for the blank),

TA B L E 1 (Continued)
in which 200 is the total sample volume, V is the added sample supernatant volume (15 μl), and 2 is the dilution factor with methanol (1:5).

| Proteotoxicity
As a proxy for proteotoxicity, the level of the stress protein 70 family (Hsp70) relative to the control was used. Frozen snails were homogenized individually on ice in extraction buffer (80 mM potassium acetate, 5 mM magnesium acetate, 20 mM HEPES, and 2% protease inhibitor at pH 7.5) according to their body mass (2 μl buffer/mg snail) and centrifuged for 10 min at 20,000 g and 4°C.
To determine the total protein content of each sample, the protein-dye binding assay of Bradford (1976)

| Morph frequencies and shell pigmentation intensity in 35 Mediterranean sites
To compare the morph communities in large populations of To determine morph frequencies, the "naturally white" (W) and the "naturally pigmented" (P) morphs were separated and counted for every population. Pigmentation analysis, which provides more detailed information on the shell's pigmentation intensity, was conducted as follows. Animals were killed in 70% ethanol, the shell

| RE SULTS
In summary, climate data revealed that the magnitude of global warming was simulated realistically in the OTCs. All biological data obtained for the different morphs showed that gradually more intense pigmentation poses increasingly strong physiological disadvantages under the tested conditions in the Mediterranean climate, particularly under conditions of simulated warming.

| Climate data, exposure conditions, and snail characteristics
Climate data during exposure in the OTCs are displayed in Figure 2.

| Warming of snails
Data on the warming of snails obtained in all three experiments (escape behavior, OTCs, and exposure experiments at 50 cm above the ground) all showed that darker shells heat up more than paler ones.
In the escape behavior experiments, darker shells became significantly warmer than paler ones, prior to moving and on their way upwards (Figure 3b, left panel). Even in contact with the hot ground, at least parts of the shell of all (naturally and artificially) pigmented morphs became significantly warmer than those of the white ones.
When moving upwards, shells of darker snails reached higher temperatures at all heights. At 20 cm above ground level (occasionally also at 6.5 cm), the soft body of artificially painted morphs (S, B) was also warmer than that of natural ones (W, P).
In addition, thermography of snails that escaped from the soil

| Behavior
The escape behavior experiments with field-collected, non-OTCexposed animals showed that the ability of snails to escape from a hot soil surface did not differ consistently among the four morphs.
The rank sum analyzes for the parameters "start to move," "start to F I G U R E 3 Theba pisana morphs and corresponding internal/shell temperatures in different situations in the field or the global warming experiment. a: Representative individuals of the morphs "naturally white" (W), "naturally pigmented" (P), "artificially striped" (S), and "artificially blackened" (B). b: Significant differences in temperature of the four morphs exposed to different situations, either at the bottom or top of the shell, at the shell surface or at/inside of the soft body. Red wedges depict differences at p ≤ .05 with the wedges' bigger part symbolizing the, respectively, higher temperature. Left panel: situation at different heights in open-top chambers, on the ground or at different heights when moving up the chambers' walls. Right panel: Situation in an inactive state on vertically arranged wooden sticks inside the open-top chambers or at 50 cm above ground on a needle thermometer in the field. Curved wedges are based on significances obtained by pairwise comparisons, and straight wedges are based on ANOVAs. When lying on hot ground thermal conduction heated up the shell (bottom) independent of its coloration. When snails flee, the heat by moving upward the chambers' walls darker shells tend to become warmer than white ones, which also correspond to the respective body surface temperatures when the snails reach positions more distant to the soil. This phenomenon was corroborated by data obtained for individuals, which deliberately climbed wooden sticks inside the opentop chambers or which were measured under field situations at a constant 50 cm distance from soil. c: Two individuals of morphs W and P of the same size and origin in direct vicinity to one another at identical distance from the soil (top), and thermographic images when exposed to solar radiation around noon (12.30 hr at 27°C air temperature, middle) or in the evening (22.15 hr at 22°C, bottom). The shell surface of the pigmented morph displayed higher temperatures at all parts, compared with the white one, both in the middle of the day and in the evening. As well, in both morphs, the side directed to the ground was warmer than the top side indicating the importance of indirect thermal radiation reflected by the ground. Increasing temperatures are depicted in the order blue < purple < orange < dark yellow < bright yellow We did not find significant differences of overall competitive effect of morph coloration on the snails for either ">10 cm" or "top position occupancy." Thus, morph color did not alter the climbing behavior of neighbors.

| Stress physiology and biochemistry
After a 3-day exposure in the OTCs, the Hsp70 level of the snails was significantly linearly correlated with the shell pigmentation intensity of the morphs and Hsp70 content (Figure 4a). The more intense the shell pigmentation, the higher was the destruction of the capacity of the Hsp70 system (linear regression, ANOVA). This correlation held true for datasets either including moribund individuals with zero Hsp70 levels (p < .0001, F = 40.1127) or excluding them (p = .0003, F = 13.2879). The Hsp70 level of laboratory controls and field-collected individuals did not differ much from OTC-exposed morphs W, P, and S, but B morphs showed a comparatively (though insignificantly) lower level of stress protein.
Concomitant with the breakdown of the Hsp70 chaperoning system, black morphs showed significantly increased lipid peroxidation after 3 days in the OTCs (Figure 4b). This effect was evident with respect to the other OTC-exposed morphs (B > [W = P = S], p between .0001 and .0039) and vs. blacks of the laboratory control (p < .0001, F = 5.6172, ANOVA and Tukey-Kramer). Also for this parameter, the laboratory controls and field-collected individuals did not differ much from OTC-exposed morphs W, P, and S, but B morphs showed a significantly elevated lipid peroxidation level (p < .0001, F = 5.6172, ANOVA and Tukey-Kramer).

F I G U R E 4
Levels of the protein-stabilizing Hsp70 stress protein family (a: optical volume relative to standard) and of lipid peroxidation (b: cumol hydroperoxide (CHP) equivalents per mg wet wt., y-axis truncated) after 3 days, and weight loss (c: % of initial wet weight) after 7 days of exposure to simulated global warming conditions in open-top chambers in the two natural and the two artificial morphs of Theba pisana (means ± SD, respectively). Average shell pigmentation intensity (±SD) of the morphs was quantified on a grayscale from 0 [black] to 255 [white]. Data obtained for laboratory controls (without exposure to direct sun at a constant 26°C) or field specimens in their natural environment at Montfavet are displayed at the right margin of the graphs (each bar represents a sample, and shaded area indicates the range between min and max). Linear regression analysis and 95% confidence intervals. Data showed a significant decrease in the chaperoning system (ANOVA, p = .0158) and a significant increase in weight (presumably water) loss (ANOVA, p = .0149) with increasing shell pigmentation. *The lipid peroxidation level in artificially blackened individuals was significantly increased (0.0001 ≤ p ≤ .0039 vs. the other morphs in the open-top chambers and p < .0001 vs. the laboratory control) In the OTC experiment, morphs differed significantly with respect to their weight loss at day 8 (Figure 4c). ANOVA/Tukey-Kramer analysis showed S and B morphs to lose more weight than W and P. Black morphs lost significantly more weight than S morphs ([W = P] < S < B, p between <.0001 and .0113, F = 46.7770, Figure 4). In the mixed models, all pairwise comparisons involving black-shelled snails were highly significant, that is, B individuals on day 8 weighed far less than snails of any other color (p < .001 for all comparisons). Pairwise comparisons between other shell colors were nonsignificant (P vs. S: p = .344; W vs. P: p = .977; W vs. S: p = .535). Generally, all OTCexposed snails had lost significantly more weight than the laboratory controls (p < .0001, F = 46.7770, ANOVA and Tukey-Kramer).

| Mortality
In the OTCs, the four morphs differed a lot in rates of mortality ( Figure 5). Already at days 3 and 5, the survival rate of white morphs Mortality in naturally pigmented P individuals was also higher than in W when exposed to the OTCs for a week (f OTC : on the average 1.87 times the W mortality with a very slight incline, Figure 5). At the field site Montfavet, the mortality of individuals, indicated by the percentage of dead snails detached from plants (not on the soil surface), was higher in the P morph (5.93% dead snails on plants, n = 540) than in W (4.01%, n = 1,620), leading to a factor f field of 1.48 times the mortality in P versus W under current climate conditions in 2017.
The time to 50% mortality was significantly different in all mul- quently followed by S (which was similar to B but significantly different), P (very different to S by 2 days), and W (2.5 days more than P).
Mixed models also revealed significant differences in the lethal competitive "response" of morph over the first nine days in the OTCs ("mortality rate," binomial GLMMs): The slopes of the models showed that B snails died quickest, followed by S (which was significantly different to B: z = −5.58; p < .0001), then W and P were the slowest, both significantly slower than S (W vs. S z = −4.422, p < .0001; P vs. S z = −4.276, p < .0001). In the mixed models, the death rate response of the P morph was a little quicker (like in the regression analyzes) but not significantly different from W (p = .269).
In terms of competitive effects, snail color hardly played a role in the mortality of neighboring snails. Most pairwise comparisons were nonsignificant. There were no significant differences for time to 50% mortality, and only one significant result for mortality rate (B had more effect than P, p = .0315). Therefore, the impact on a neighbors' mortality was generally independent of morph pigmentation.
In respect to the relevance of size for survival, we did not find early mortality of particularly small or large individuals in any morph (Welch's test, Bonferroni-corrected).

| Morph frequencies and pigmentation intensity in Mediterranean populations
In  (Table 1). Both parameters showed a significant correlation with the mean temperature of the hottest quarter of the year (Figure 6b,c). The frequency of the white morph increased with increasing ΔT and followed a saturation curve approaching almost 100% whites (p = .0015, W = 7.9873, ANOVA). In addition, T. pisana from sites warmer than Montfavet were found to be on average paler with increasing ΔT (p < .0001, W = 35.1335, ANOVA).
Based on the observed mortality rates in the two natural morphs, W and P, in the OTCs and the field site at Montfavet, we extrapolated the morph frequencies and the mean pigmentation intensity at this site assuming an additional temperature increase of 1.12°C.
According to the calculated values for f OTC and f field, the average mortality of the P morph (vs. W) will further increase by the factor of 1.264 (f warming ). Therefore, an increase in the average temperature by 1.12°C will result in a shift in the frequency of naturally pigmented individuals from today's 20.58 ± 2.74% to 16.28 ± 2.17% and for the white morph from 79.42 ± 2.74% to 83.72 ± 2.17%. Thus, with 1.12°C warming we predict a loss of 20.9% of pigmented individuals at this site and a slight increase in the average pigmentation value by 0.48% toward "white" (arrows in Figure 6b,c).

| D ISCUSS I ON
The role of shell pigmentation in the physiology and evolution of land snails is among the best-studied and most discussed topics in animal evolutionary ecology. Whereas numerous publications have addressed aspects of the impact of thermal stress, they mainly focused on hν absorption and heat transfer or on the distribution of morphs in particular habitats (reviewed by Schweizer et al., 2019).
In our current study, we highlight three innovative aspects that have not been previously addressed: • The joint analysis of multiple parameters, collated in a syndrome across different levels of biological organization (thermodynamics, Temperature data for (b) and (c) were obtained from WorldClim.org. Nonlinear or linear regression curves and 95% confidence intervals (red) display a strong trend toward a dominance of paler/white individuals at sites with increasing summer temperatures. The green circle displays the shift of the population at Montfavet as calculated from the mortality results obtained in this study for the simulated warming of 1.12°C. The rather exceptional current position of the Montfavet data point (which will shift to the average trend with ongoing global change) is in accordance with the exceptionally high degree of thermal pressure change this population has been facing during the last 30 years: +1.34°C, in comparison with an average of +1.1°C in the entire Mediterranean biochemistry, physiology, behavior, survival rates, and the population frequencies of shell banding morphs), allowing us to trace thermal effects from the molecular to the population level.
• An experimental approach to mechanistically show the effects of environmental warming on land snails, which has not been conducted before in such complexity.
• The extrapolation of our results to predict the consequences of environmental warming for future morph frequencies in a field population in the Mediterranean.
Our results, obtained in experiments tailored to suit our specific research questions, consistently suggest that T. pisana individuals with pigmented shells, in both the current Mediterranean climate and under simulated environmental warming, become warmer than white ones, thus experiencing additional thermal stress. This warming effect was most pronounced in the artificially painted individuals, but even the shell and body of the naturally pigmented morph were significantly warmer than those of nonpigmented individuals. Previous studies investigating heat stress on Theba, Cepaea, and Littoraria (Cook & Freeman, 1986;Hazel & Johnson, 1990;Heath, 1975;Jones, 1973;Knigge et al., 2017) generally compared fully banded with unbanded morphs. Our results not only confirm the findings from those studies, but also elaborate on them to include strong experimental manipulations of color that define the mechanism of the warming. Another recent study has found the same trend for faintly pigmented T. pisana individuals versus white ones, even under moderate illumination in the laboratory (Tull et al., 2020). At least under Mediterranean climate conditions, at present and in the future, there is no doubt that dark individuals will heat up more intensely than pale ones.
Evidence is thus mounting that differently pigmented morphs experience different shell and body temperature in a given situation: Darker morphs warm up more quickly and more intensely. It has long been known that snails respond to elevated body temperature by selecting preferred substrate temperatures (Chapperon & Seuront, 2011;Riddle, 1990), or by climbing tall objects to escape from substrates with extreme temperatures (Cowie, 1985;Di Lellis et al., 2012;McQuaid et al., 1979). These responses may to a certain degree buffer the physiological impacts of heat (Ng et al., 2017).
Climbing objects often results in aggregates of multiple individuals. Depending on the climate, the aggregate centers are either warmer (in cold climate, Chapperon & Seuront, 2012) or cooler (in hot climate, McQuaid et al., 1979) than the aggregate surfaces. An alteration in temperature due to climbing or within aggregates implies that there may be intraspecific competition to find the "best place" for temperature regulation. However, Chapperon and Seuront (2012) reported that all snails of a cluster benefit from the aggregate.
The question thus arises whether either fighting for an optimal resting place or joint action (i.e., "co-operation") between individuals has been evolutionary favored, or whether any trade-offs in other traits exist. Particularly in the Mediterranean region, differently colored snails may need to possess different abilities to escape the soil surface. However, we did not find any morph-specific behavior. While darker morphs heated up more quickly, they did not flee the hot surface any quicker than pale morphs, nor did they climb more quickly, either during observed "racing" experiments, or as determined from the number of top positions occupied on the skewers in the OTCs.
Yet, by analyzing the earlier time points-which takes into account mortality-we found minor evidence in the opposite direction to the "fleeing heat" hypothesis. In general, while very few results were significant, the evidence suggested that paler individuals generally reached higher positions on the skewers than darker individuals, but with little effect on overall mortality. There was even less evidence to suggest that morph coloration had any consistent competitive effect (in terms of neither climbing nor mortality). Based on our results, we thus conclude that competition among morphs is at most indirect and may just be a consequence of limited space, supporting no particular morph.
Among the main biochemical effects of elevated body temperature is the impaired integrity of intracellular proteins and cell membranes. As protein integrity is mandatory for a cell's survival, a biochemical chaperoning machinery that assists in protein (re) folding has been established early in the evolution of life. A central component of this machinery is the Hsp70 protein family, whose induction by heat and drought has been well characterized in land snails Di Lellis et al., 2014;Dieterich et al., 2013Dieterich et al., , 2015Gaitán-Espitia et al., 2013;Köhler et al., 2009;Kotsakiozi et al., 2015;Mizrahi et al., 2010Mizrahi et al., , 2011Mizrahi et al., , 2012aMizrahi et al., , 2012bMizrahi et al., , 2015Mizrahi et al., , 2016Scheil et al., 2011;. As long as an induction of Hsp70 may compensate for the proteotoxic impact of body temperature elevation, organisms can limit pathologies and ensure survival. A breakdown of the Hsp70-based refolding machinery reflected by declining levels of this chaperone, however, is indicative for oxygen limitations in the body fluids at high temperatures, resulting in oxidative stress (Pörtner, 2002) and subsequent destructive processes in tissues that, on the long run, lead to elevated morbidity (Eckwert et al., 1997). In our OTC experiment, we found that simulated environmental warming gradually decreased the Hsp70 level in morphs with increasing pigmentation intensity, symbolizing a collapsing chaperoning system, particularly in the blackened snails but also, tentatively, in the other morphs. Concomitant with this temperature effect on proteins, the lipid peroxidation level increased with increasing pigmentation intensity of the morphs as a consequence of elevated respiration and thus oxygen consumption in snails with elevated body temperature (Pörtner, 2002). The disruption of redox homeostasis, however, generally is a key phenotype of pathological conditions in many organisms and cell types, as excessive lipid peroxidation alters the physical properties of cellular membranes and can result in covalent modification of proteins and nucleic acids (Gaschler & Stockwell, 2017). The association between high temperatures and lipid peroxidation in land snails has experimentally been shown before in Xeropicta derbentina individuals that have been exposed to 38-48°C for 8 hr (Dieterich et al., 2015;, but it is noteworthy that, in our study, even a rather moderate increase in elevated temperatures induced lipid peroxidation in a morph-specific manner. In addition to the biochemical stress posed by elevated body temperature, all animals exposed to simulated environmental warming lost more weight than the controls. Weight loss in this context is almost exclusively due to the loss of water, and the impact of catabolic processes is neglectable within a rather short time frame (Reuner et al., 2008). Particularly, the artificially colored morphs B and S suffered from water loss >30% (wet wt.) within a few days, confirming the results of Moreno-Rueda (2008) who has painted Sphincterochila candidissima black. Highly effective water retention, however, seems to be essential for land snails living in arid and semi-arid environments, which makes for instance the desert snail Sphincterochila boissieri to lose only 0.45 mg water per day, allowing survival of up to 4 years in the drought (Schmidt-Nielsen et al., 1971). Kotsakiozi et al. (2015) showed that only one of six Codringtonia species lose water during summer estivation, and Cantareus apertus, which estivated for 6 months had only 14% less water than active individuals (Reuner et al., 2008). Plausibly, the higher internal temperature in the intensely (and artificially) colored morphs of our study caused a higher need for evaporative cooling, which considerably decreased the water content of the snails' soft bodies.
The syndrome of thermal stress, which comprises higher internal temperature, weakening of the protein chaperoning system, enhanced peroxidation of membranes, and increased loss of water in the pigmented morphs, feasibly resulted in differences in the survival rate among morphs. This was particularly obvious under conditions of environmental warming in the OTCs, but was already visible in the field. Even though, in our experiments, inactivity of climbed snails was disturbed once a day for a few minutes-an additional stress factor which is potentially relevant also in the natural habitat-the relative difference in water loss and mortality among the morphs was striking: The OTC experiments revealed that, in a situation that simulated the global warming effect, the brighter the shell of an individual the higher is the chance for its survival. This was independent of its size or neighboring morph, provided that there was enough space to climb up and escape the hot soil. Thus, the capacity of a morph to occupy top positions on a vertical object did not determine survival but was rather a consequence of it. The experimental design using artificial painting helped to reveal the general association between darker pigmentation and higher mortality, but this effect was not restricted to the artificially colored snails, but also visible in the naturally pigmented morph P, which had a 1.87 times higher mortality rate (and reached 50% mortality 2.5 days earlier) than the unpigmented morph W under simulated environmental warming conditions.
Even in Montfavet's climate of 2017, which was a snapshot in the ongoing process of global warming, the naturally pigmented morphs displayed a 1.48 times higher mortality than their white conspecifics. As T. pisana reproduces at the end of its lifetime in autumn and early dying individuals do not reproduce at all, such a difference in survival is of crucial relevance for the population structure. Hence, in their entirety, our results imply that the association between shell pigmentation intensity and thermal effects reflects causality for thermal selection pressure on pigmentation. To establish evidence for cause-consequence relationships between these parameters and across the different levels of biological organization (Triebskorn et al., 2003), we applied Hill's criteria of causation (Hill, 1965). The following criteria strengthened our interpretation.
• Strength of association: Pigmentation and a high number of thermal effects were significantly correlated at p-values as low as <.0001.
• Coherence: Thermal effects all pointed toward the same direction: The darker the shell, the higher the internal and surface temperatures, the physiological challenge on protein and membrane integrity, the water supply, and the survival. For some of these parameters, also different methods of data analysis resulted in the same findings. Furthermore, results from OTC experiments and field observations were coherent.
The consequences of this selection pressure are already reflected in the large field population at Montfavet that shows a higher mortality rate in naturally pigmented individuals than in white ones, even though the pigmentation intensity in the P morph was not very intense, compared with other Mediterranean populations.
Evolutionary aspects of thermal selection on shell coloration in land pulmonates have been addressed in the past, however, without investigating the thermodynamic, biochemical, and physiological fundamentals. There is consensus that snail populations in shaded habitats, less dominated by solar radiation, evolve toward darker shells than those in open habitats, and this effect has been reported for a number of species (Cain & Sheppard, 1954;Heller, 1981;Ożgo & Komorowska, 2009;Schilthuizen, 2013), despite remaining uncertainties about the (local) role of shell banding in, for example, thermoregulation (Cameron & Cook, 2012;Kerstes et al., 2019;Silvertown et al., 2011). The debate on the underlying mechanism for this phenomenon-thermal selection versus visual predation pressure by birds and rodents-is still ongoing (reviewed in Schweizer et al., 2019), but known geographic clines of particular phenotypes of land snails (Cowie, 1990;Jones et al., 1977) signify the relevance of climate parameters in this respect.
What will be the consequences for T. pisana and their different morphs at Montfavet and, eventually, the entire Mediterranean? It is highly unlikely that this species will go extinct in Southern France provided that environmental warming will be limited to about 2.5°C. This geographical area is not located at the edge of the environmental limits of T. pisana, and it has recently been shown (for 10,000 Arabidopsis thaliana genome-wide SNPs) that the strongest climate-driven selection is currently experienced by populations living at such edges (Exposito-Alonso et al., 2019). Furthermore, in 136 case studies of possible climate change-related extinctions, none of the proximate causes of global change have shown a relationship between local extinction and limited tolerance to high temperatures (Cahill et al., 2013), possibly because life history, behavior, and physiology can be altered by microevolution in response to environmental warming (Debecker & Stoks, 2019). The only proposed extinction of a snail, Rhachistia aldabrae, as a result of climate change had to be revised because this species was later rediscovered (Gerlach, 2007).
The variation in T. pisana's shell pigmentation will probably persist in the future, either as a result of year-to-year change in the intensity of physiological pressure (Johnson, 2011(Johnson, , 2012 or as a result of buffering mechanisms within the limits of phenotypic plasticity . However, we expect a decrease in the frequency of pigmented individuals in the Montfavet population with increasing temperatures. This prediction is based upon our results on the higher physiological pressure on the pigmented morph and the higher frequency of white individuals, concomitant with an average decrease in pigmentation intensity in 34 other T. pisana populations in the Mediterranean, which is related to increasing temperatures in the hottest quarter of the year. In our model, we solely used the selection against darker morphs to predict the future average pigmentation intensity of this population, even though phenotypic plasticity will most likely modify morphological variation also in a future population at this site. In this context, the broadness of heritable ranges of phenotypic plasticity in different genotypes will determine the number of genotypes that will go extinct in a warming environment and thus possibly limits the ability of the population to restore the full spectrum of phenotypic variation . Quantitative information on the broadness of morphological reaction norms of given phenotypes, unfortunately, is completely lacking for T. pisana, which prevents inclusion of this aspect in a mathematical model. However, the close fit of the selection-based prediction for the Montfavet population to the pigmentation data for already warmer sites indicates that the selection of phenotypes very probably outweighs phenotypic plasticity in a warming Mediterranean. Even in the colder climate of Wales, two thirds of the variation in the pigmentation pattern of T. pisana, documented by Cowie (1992), could be explained by the genetic basis and only about one third by environmental factors .
We also propose that the selection against pigmented individuals at Montfavet has already been ongoing in the last years, because the pigmentation intensity in this population is comparably low and current differences in the morphs' mortality rates are evident.
In addition, the Avignon region has experienced exceptionally high warming during the last three centuries compared with the average in the Mediterranean, which already is a hot spot of global warming.
Consequently, further temperature increases will probably not lead to the extinction of T. pisana in this area but will predictably reduce intraspecific diversity.

CO N FLI C T O F I NTE R E S T
None declared.

AUTH O R CO NTR I B UTI O N S
Heinz-R. Köhler: Conceptualization (equal); Data curation (equal);