Exposing organisms to a particular stressor may enhance tolerance to a subsequent stress, when protective mechanisms against the two stressors are shared. Such cross-tolerance is a common adaptive response in dynamic multivariate environments and often indicates potential co-evolution of stress traits. Many aquatic insects in inland saline waters from Mediterranean-climate regions are sequentially challenged with salinity and desiccation stress. Thus, cross-tolerance to these physiologically similar stressors could have been positively selected in insects of these regions. We used adults of the saline water beetles Enochrus jesusarribasi (Hydrophilidae) and Nebrioporus baeticus (Dytiscidae) to test cross-tolerance responses to desiccation and salinity. In independent laboratory experiments, we evaluated the effects of (i) salinity stress on the subsequent resistance to desiccation and (ii) desiccation stress (rapid and slow dehydration) on the subsequent tolerance to salinity. Survival, water loss and haemolymph osmolality were measured. Exposure to stressful salinity improved water control under subsequent desiccation stress in both species, with a clear cross-tolerance (enhanced performance) in N. baeticus. In contrast, general negative effects on performance were found under the inverse stress sequence. The rapid and slow dehydration produced different water loss and haemolymph osmolality dynamics that were reflected in different survival patterns. Our finding of cross-tolerance to salinity and desiccation in ecologically similar species from distant lineages, together with parallel responses between salinity and thermal stress previously found in several aquatic taxa, highlights the central role of adaption to salinity and co-occurring stressors in arid inland waters, having important implications for the species' persistence under climate change.

The persistence of animal populations in dynamic and multivariate environments greatly depends on their ability to deal with the interactive effects of different stressors occurring simultaneously or sequentially over short time scales (Gunderson et al., 2016). Exposure to a particular stressor might enhance tolerance to a subsequent stress if the physiological protective mechanisms against the two stressors are shared (cross-tolerance, e.g. Elnitsky et al., 2009; Holmstrup et al., 2002) or, conversely, can cause the organisms to be more susceptible to the second stress (cross-susceptibility) (Sinclair et al., 2013; Todgham and Stillman, 2013).

Studies on multiple stressors have received increasing attention for their potential to reveal interesting information, which would be difficult to predict based on single stressor approaches (DeBiasse and Kelly, 2016; Gunderson et al., 2016), as well as to increase our understanding of responses to global change in natural multivariate environments (Hewitt et al., 2016). However, these approaches are still scarce in the literature and are mostly focused on the combined effects of temperature with other factors (e.g. Huth and Place, 2016; Pansch et al., 2012; Todgham et al., 2005).

Inland saline waters in arid and semi-arid regions represent a template for the evolution of mechanisms to deal with multiple sources of stress, especially high and fluctuating temperatures and water salinity levels coupled with seasonal wet and dry periods (Gasith and Resh, 1999; Millán et al., 2011). In temporary saline waterbodies, droughts are often preceded by an increase in salinity as the water level drops (Hershkovitz and Gasith, 2013); therefore, many aquatic organisms living in these environments are sequentially and repeatedly challenged with salinity and desiccation stress. Some aquatic insects face droughts in situ in microrefuges (e.g. Stubbington et al., 2016), whereas others with flying adults disperse to more favourable wet habitats (Bilton, 2014; Strachan et al., 2015). During dispersal, aquatic insects may experience important water losses associated with exposure to the desiccating aerial medium and flight activity (Dudley, 2000). The persistence of these species greatly depends on the ability of adults to deal with such osmotic and dehydration stress, especially for those with no desiccation-resistant stages (eggs or larvae), like most true water beetles (sensuJäch and Balke, 2008) (Millán et al., 2014).

Both salinity and desiccation lead to dehydration and osmotic stress, which is a critical problem at the cellular level (Bradley, 2009; Cohen, 2012; Evans, 2008). Therefore, salinity and desiccation stress in insects trigger common physiological mechanisms, mainly aimed at increasing water content (e.g. drinking from the medium), avoiding its loss (e.g. control of cuticle permeability) and maintaining ionic homeostasis (e.g. activity of Malpighian tubules and specialized parts of the hindgut) (Bradley, 2009; Dow and Davies, 2006; Gibbs and Rajpurohit, 2010; Larsen et al., 2014). However, the efficiency of these mechanisms in terms of water and ionic balance under sequential exposure to both stressors, as frequently occurs in nature, has not been studied in saline aquatic insects. If the stressors occur close enough in time, so that exposure to the second stressor takes place while the physiological response to the first stressor is still being mounted, interactive effects that result in cross-tolerance or cross-susceptibility are likely to occur (Gunderson et al., 2016).

List of abbreviations
     
  • CD

    control desiccation

  •  
  • M0

    initial mass

  •  
  • Md

    fresh mass after desiccation pre-treatment

  •  
  • Ms

    fresh mass after salinity pre-treatment

  •  
  • OS

    optimum salinity

  •  
  • RD

    rapid desiccation

  •  
  • SD

    slow desiccation

  •  
  • SLS

    sublethal salinity

  •  
  • WC0

    initial water content

  •  
  • WCd

    water content after desiccation pre-treatment

  •  
  • WCs

    water content after salinity pre-treatment

  •  
  • WLR

    water loss rate (mean, maximum, final)

Cross-tolerance is not necessarily adaptive per se (i.e. it might be the consequence of general responses to stressors that are not experienced together in nature), but in many cases it appears to be under selection in response to synchronous or sequential stressors (e.g. Chen and Stillman, 2012; Kumlu et al., 2010; Sánchez-Fernández et al., 2010; Todgham et al., 2005). Therefore, identifying cross-tolerance responses could offer significant information on the evolutionary history of interactions among stressors (Bubliy et al., 2012; Sinclair et al., 2013). In a recent study on a lineage of water beetles (Enochrus species of the subgenus Lumetus), Arribas et al. (2014) found that transitions from fresh to saline waters occurred in periods of global aridification and showed a positive correlation between the salinity and aridity of species' habitats. From these results, the authors hypothesized a correlated evolution of salinity and desiccation tolerance in this group, potentially as an exaptation process due to linked physiological mechanisms to deal with the two stressors. The positive association found between desiccation resistance and salinity tolerance across species of this genus (Pallarés et al., 2016) also points in that direction. In light of this potential evolutionary link between salinity and desiccation tolerance and their common physiological basis, cross-tolerance might have been selected in species living in saline inland waters where salinity and desiccation stress co-occur.

Here, we tested cross-tolerance to salinity and desiccation in two Iberian water beetle species with a clear habitat preference for saline waters: Enochrus jesusarribasi and Nebrioporus baeticus. They belong to two representative coleopteran families in inland waters (Hydrophilidae and Dytiscidae) that have colonized and diversified across aquatic habitats independently (Hunt et al., 2007). We examined survival and two key measures of osmoregulation and water conservation capacity (water loss and haemolymph osmolality) under controlled-sequential exposure to salinity and desiccation in the laboratory. We predicted that the activation of mechanisms for ionic and water control during exposure to stressful but non-lethal levels of either salinity or desiccation would help specimens to deal with further osmotic-dehydration stress, improving performance under a subsequent exposure to the other stressor.

Our results on the comparative physiology of two ecologically similar species from different lineages of water beetles show the importance of the adaptation to both stressors in inland saline waters and provide new insights into the processes of colonization and diversification in these systems. This information could also be highly relevant to understanding how aquatic insects respond to the ongoing aridification in Mediterranean inland aquatic ecosystems, where more extreme and prolonged droughts and increased salinity levels are predicted (Bonada and Resh, 2013; Filipe et al., 2013; Sala et al., 2000).

Target species, specimen collection and maintenance

Adult specimens of the water beetle species Enochrus jesusarribasi Arribas and Millán 2013 (suborder Polyphaga, family Hydrophilidae) and Nebrioporus baeticus (Schaum 1864) (suborder Adephaga, family Dytiscidae) were used as models for this cross-tolerance study. Both species have been proposed as ‘vulnerable’ because of their strong habitat specificity and endemic character, restricted to inland saline streams in southeastern Spain (Arribas et al., 2013; Millán et al., 2014; Sánchez-Fernández et al., 2008).

These species have been shown to be effective euryhaline osmoregulators in laboratory assays (Céspedes et al., 2013; Pallarés et al., 2015) and are mainly found in meso-hypersaline waters in nature constituting high-abundance populations (Velasco et al., 2006). Dispersal by flying in the adult stage is the main strategy for coping with seasonal droughts in these species, whose larvae and eggs are desiccation sensitive.

Adult beetles were collected in two intermittent saline streams located in Murcia (SE Spain): Rambla Salada (E. jesusarribasi) and Río Chícamo (N. baeticus). Specimens were held in the laboratory at 20°C with substratum and water taken from the collection sites, at conductivities of 65 and 12 mS cm−1 for E. jesusarribasi and N. baeticus, respectively, until they were used for the experiments. During this time, beetles were fed with macrophytes (E. jesusarribasi) or chironomid larvae (N. baeticus).

Cross-tolerance experiments

We conducted two independent experiments to assess the effects of (i) exposure to stressful salinity on the subsequent resistance to desiccation and (ii) exposure to desiccation stress on the subsequent tolerance to salinity. The stress treatments represented sublethal conditions, which were adapted to the specific tolerance ranges of each species (see specific conditions in Fig. 1) according to pilot trials and previous studies (Arribas et al., 2012b; Céspedes et al., 2013; Pallarés et al., 2015, 2016; Sánchez-Fernández et al., 2010). Pilot trials were also conducted to determine the number of replicates (i.e. specimens) needed to ensure adequate power to detect the effect of the stressors.

Fig. 1.

Experimental conditions and variables obtained for each experimental step. (A) Effect of salinity on desiccation resistance. (B) Effect of desiccation on salinity tolerance. Experimental conditions [relative humidity (RH) osmolality and exposure time] that differ between the species are indicated with species symbols. OS, optimum salinity; SLS, sublethal salinity; CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

Fig. 1.

Experimental conditions and variables obtained for each experimental step. (A) Effect of salinity on desiccation resistance. (B) Effect of desiccation on salinity tolerance. Experimental conditions [relative humidity (RH) osmolality and exposure time] that differ between the species are indicated with species symbols. OS, optimum salinity; SLS, sublethal salinity; CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

The experiments were conducted at a constant temperature (20°C) and light:dark cycle (12 h:12 h) in an environmental chamber with humidity control (CLIMACELL-404, MMM Medcenter Einrichtungen GmbH, Germany). At the end of each experiment, dry mass (Mdry) of all the tested specimens was measured (after drying at 50°C for 48 h) with an electronic high-precision balance (±0.00001 g) and beetles were sexed by examination of genitalia.

Effect of salinity on desiccation resistance

Groups of 30–40 specimens of each species were randomly assigned to the following pre-treatments, for 1 week (Fig. 1): (i) optimum salinity (OS), i.e. the most frequent salinity levels of each species' habitat or (ii) a higher sublethal salinity (SLS). Solutions were prepared by dissolving the appropriate amount of marine salt (Ocean Fish, Prodac, Citadella Pd, Italy) in distilled water. Food was provided for the first 5 days and removed 48 h prior to desiccation exposure. After the salinity pre-treatments, surviving specimens of each salinity that showed no sign of critical stress (i.e. were able to normally move) were used to obtain fresh mass (Ms; mg) and their survival and water loss under desiccation were investigated. For this, specimens were gently dried on blotting paper, held for 10 min at room temperature until the cuticle surface was totally dry and then individually placed in open glass vials and subjected to 40±5% relative humidity (RH) in an environmental chamber. Survival was checked and specimens were re-weighed every 2 h. Using such fresh mass measures and Mdry of each specimen (see above), we estimated their water content after salinity pre-treatment (WCs; mg), the rate of total water loss (i.e. cuticular, respiratory and excretory) under desiccation (WLR; mg h−1) and the final water content remaining at death (WCf; mg).

Kaplan–Meier survivorship curves (Altman, 1992) were used to compare survival to desiccation between salinity pre-treatments, specifying right censored data for those individuals that were alive at the end of the experiment. Differences in mean survival between optimum and stressful salinity pre-treatments were tested using the log-rank (Mantel–Cox) test (e.g. Folguera et al., 2011; Kefford et al., 2012). Generalized linear models (GLMs) were used to test for differences between the two salinity pre-treatments in WCs and WLR under the subsequent desiccation. Because the change in WLR was not linear over desiccation exposure, this variable was estimated for each individual as (i) the mean of the rates measured every 2 h (WLRmean), (ii) the maximum rate (WLR­max) and (iii) the final rate (WLRf; i.e. the rate measured at the interval prior to dead). To correct the analyses for the individual variation in mass and water status and to account for potential sex-specific differences in the response variables (see Chown and Nicolson, 2004; Le Lagadec et al., 1998), initial Ms and sex were included as covariates, plus WCs in the analysis of WLRs. In this case, as Ms and WCs are highly correlated, their effects were evaluated separately to avoid statistical problems of collinearity between predictors.

Effect of desiccation on salinity tolerance

Groups of 50 specimens of each species were randomly assigned to the following pre-treatments (specific exposure times for each treatment and species are shown in Fig. 1): (i) non-desiccation control (CD; at RH>90%), (ii) rapid desiccation (RD; at 10% RH) or (iii) slow desiccation (SD; at 40% RH). For the SD pre-treatment, specimens were placed in individual open glass tubes in the environmental chamber set at 40±5% RH. For the RD pre-treatment, 2 g of silica gel were added to individual glass tubes and separated from the specimen using a piece of foam. The tubes were covered with Parafilm to maintain a low RH (10±5%) (e.g. Bazinet et al., 2010; Lyons et al., 2014). In the control, the open glass tubes were introduced into a 7 l plastic aquarium with deionized water in the base (approximately 2 cm) enclosed with plastic film, reaching RH levels close to saturation (i.e. >90%).

Food was removed from the maintenance aquaria 48 h prior to desiccation and withheld for the duration of the experiments. All specimens were initially weighed (M0) and after desiccation, survival and ability to move were checked, discarding individuals that showed movement difficulties. A subgroup of 20–30 specimens were re-weighed (fresh mass after desiccation, Md) and used for estimation of WC0 and WCd (before and after the desiccation pre-treatment, respectively) as well as WLR. These same specimens were subsequently used for the survival assay under stressful high salinity. For this purpose, specimens were placed in individual plastic containers with 40 ml of the specific SLS solution for each species (Fig. 1). Survival was checked every hour for the first 12 h and subsequently at 12 h intervals, for 1 week.

Survival during the salinity exposure was compared among the three desiccation pre-treatments following the procedure explained for the effect of salinity (see above). In this case, interval censored data were specified in Kaplan–Meier curves for those deaths registered between 12 h intervals. WLR and WCd were compared among pre-treatments using GLMs and Bonferroni post hoc pairwise comparisons in order to quantify the magnitude of desiccation stress and determine how it affected survival under the subsequent salinity exposure. M0, sex and WC0 were included as covariates.

In parallel, another subgroup of 18–24 specimens per species exposed to the same three desiccation pre-treatments was used to obtain haemolymph samples immediately after desiccation (time 0) and after 8 and 24 h of salinity exposure (N=6–8 specimens per species per time) (Fig. 1). Haemolymph extraction was conducted following the procedures described in Pallarés et al. (2015) and osmolality was obtained using a calibrated nanolitre osmometer (Otago Osmometers, Dunedin, New Zealand) (see details of the same measurement procedure in Williams et al., 2004). Some haemolymph samples were discarded because the volume was insufficient or they showed a dark colour indicating potential oxidation, which resulted in a small sample size for some groups (n<5). Therefore, non-parametric tests (Kruskall–Wallis and Dunn's post hoc multiple comparison test) were used to compare the osmotic concentration of haemolymph after desiccation among treatments and its temporal variation during salinity exposure.

The number of replicates in statistical analyses was equal to the number of tested specimens in individual vials (desiccation exposure) or containers (salinity exposure). Gaussian distribution and identity link function were assumed in all GLMs. Models were validated by graphical inspection of residuals versus fitted values to verify homogeneity and QQ plots of the residuals for normality (Zuur et al., 2009). All analyses were implemented in R v.3.2.2 using the packages stats, phia and survival.

All the measurements from our experiments as well as the variables estimated for analyses (as both percentages and absolute units in the case of WC and WLR) are supplied in Tables S1–S3.

Effect of salinity stress on desiccation resistance

Most of the individuals exposed to the sublethal salinities survived (around 70% in both species) and showed no signs of critical stress (i.e. were able to move normally).

Specimens' performance under desiccation did not differ between salinity pre-treatments in E. jesusarribasi (log rank test: χ2=2.2, P=0.135), although survivorship curves showed a tendency of higher survival in individuals from the SLS pre-treatment than those from the OS pre-treatment (Fig. 2). The mean body water content of specimens after salinity pre-treatments (WCs) was significantly higher in the sublethal (5.66±0.14 mg) than in the OS pre-treatment group (4.92±0.14 mg; Table 1). WLR tended to decrease during exposure to desiccation (Fig. S1). Mean WLR did not differ between salinity pre-treatments (0.113±0.004 and 0.117±0.006 mg h−1 in the SLS and OS groups, respectively). WLRmax and WLRfinal did not differ either (Table 2; Table S4). The mean water content at death (WCf) was 2.58±0.07 mg (i.e. approximately 52% of WCs).

Fig. 2.

Kaplan–Meir survivorship curves for exposure to desiccation following salinity pre-treatment. Each data point represents survival probability (mean±s.e.m.). Numbers in parentheses indicate the number of specimens in each group. OS, optimum salinity; SLS, sublethal salinity.

Fig. 2.

Kaplan–Meir survivorship curves for exposure to desiccation following salinity pre-treatment. Each data point represents survival probability (mean±s.e.m.). Numbers in parentheses indicate the number of specimens in each group. OS, optimum salinity; SLS, sublethal salinity.

Table 1.

GLM results on the differences in water content between salinity pre-treatments, and influenceof initial body mass (Ms) and sex

GLM results on the differences in water content between salinity pre-treatments, and influence of initial body mass (Ms) and sex
GLM results on the differences in water content between salinity pre-treatments, and influence of initial body mass (Ms) and sex
Table 2.

GLM results on the differences in mean, maximum and final water loss rate under desiccation between salinity pre-treatments, and influence of initial water content and sex

GLM results on the differences in mean, maximum and final water loss rate under desiccation between salinity pre-treatments, and influence of initial water content and sex
GLM results on the differences in mean, maximum and final water loss rate under desiccation between salinity pre-treatments, and influence of initial water content and sex

A clear cross-tolerance response was observed in N. baeticus; individuals exposed to SLS showed higher survival than those exposed to OS (χ2=6.5, P=0.011; Fig. 2). Similar to E. jesusarribasi, WCs after the SLS exposure (5.17±0.10 mg) was higher than that following OS exposure (4.83±0.08 mg; Table 1). The change in WLR with exposure time was not linear, reaching a maximum at 4 h in both pre-treatments (Fig. S1). In agreement with survival patterns, WLRmean and WLRf during desiccation were slightly but significantly lower in individuals from the SLS pre-treatment group (WLRmean: 0.226±0.012 mg h−1) when compared with those from the OS group (WLRmean: 0.262±0.069 mg h−1). However, WLRmax did not differ between treatments (Table 2; Table S4). The mean water content at death (WCf) was 2.53±0.06 mg (i.e. approximately 59% of WCs).

Effect of desiccation stress on salinity tolerance

Desiccation pre-treatments did not cause significant mortality in either species; only a few specimens of N. baeticus (<10% of the exposed individuals) died during the SD pre-treatment.

In E. jesusarribasi, survival under stressful salinity conditions after SD exposure showed a rapid decline when compared with that of individuals exposed to RD or not desiccated (CD) (log rank test CD versus SD: χ2=14.4, P<0.001; Fig. 3). Although differences in mean survival between the RD and CD pre-treatments were not significant (log rank test: χ2=1.2, P=0.277), Kaplan–Meier survivorship curves showed a better performance in individuals previously subjected to RD during the first 12 h of salinity exposure. However, after 72 h this pre-treatment showed a higher mortality than the control (Fig. 3). The highest WLRs were recorded in the RD group (in accordance with the nature of this treatment), but specimens lost a significantly higher amount of water (lower WCd) in the longer SD pre-treatment group (Table 3, Fig. 4; Table S5). Haemolymph osmolality differed among treatments (χ2=10.1, P=0.006) and exposure time (χ2=7.5, P=0.023). Specimens' osmotic concentration at time 0 (i.e. after the pre-treatment) in the RD pre-treatment group was higher than that in both the CD and the SD group (P<0.05 in Dunn's post hoc comparison), but remained stable during salinity exposure. In contrast, in specimens from the SD pre-treatment and the CD group, haemolymph osmolality significantly increased during salinity exposure (P<0.05) (Fig. 5).

Fig. 3.

Kaplan–Meir survivorship curves for exposure to salinity following desiccation pre-treatment. Each data point represents survival probability (mean±s.e.m.). Numbers in parentheses indicate the number of specimens in each group. CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

Fig. 3.

Kaplan–Meir survivorship curves for exposure to salinity following desiccation pre-treatment. Each data point represents survival probability (mean±s.e.m.). Numbers in parentheses indicate the number of specimens in each group. CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

Table 3.

GLM results on variation in WLR and final water content between desiccation pre-treatments and influence of initial water content and sex

GLM results on variation in WLR and final water content between desiccation pre-treatments and influence of initial water content and sex
GLM results on variation in WLR and final water content between desiccation pre-treatments and influence of initial water content and sex
Fig. 4.

WLR and WCd following each desiccation pre-treatment. WLR, water loss rate; WCd, water content after desiccation pre-treatment. Box length represents the interquartile range (IQR) of the data and whiskers are 1.5 times the IQR. Data outside this range are represented as points. Letters indicate significant differences between treatments (Bonferroni post hoc comparisons, P<0.05) and numbers in parentheses indicate the number of specimens in each group.

Fig. 4.

WLR and WCd following each desiccation pre-treatment. WLR, water loss rate; WCd, water content after desiccation pre-treatment. Box length represents the interquartile range (IQR) of the data and whiskers are 1.5 times the IQR. Data outside this range are represented as points. Letters indicate significant differences between treatments (Bonferroni post hoc comparisons, P<0.05) and numbers in parentheses indicate the number of specimens in each group.

Fig. 5.

Haemolymph osmolality measured immediately after desiccation (time 0) and after 8 and 24 h of exposure to sublethal salinity. Salinity was 2470 and 1580 mOsmol kg−1 for E. jesusarribasi and N. baeticus, respectively. Bars represent means±s.e.m. and numbers in parentheses indicate the number of samples in each group. CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

Fig. 5.

Haemolymph osmolality measured immediately after desiccation (time 0) and after 8 and 24 h of exposure to sublethal salinity. Salinity was 2470 and 1580 mOsmol kg−1 for E. jesusarribasi and N. baeticus, respectively. Bars represent means±s.e.m. and numbers in parentheses indicate the number of samples in each group. CD, control desiccation; SD, slow desiccation; RD, rapid desiccation.

In N. baeticus, the specimens exposed to the RD and SD pre-treatments showed lower survival to salinity than those from the CD group (log rank test CD versus SD: χ2=4.1, P=0.043; CD versus RD: χ2=6.5, P=0.011). Such a decline in performance occurred in both treatment groups after 12 h of exposure to salinity (Fig. 3). The RD pre-treatment produced a significantly higher WLR than the SD and CD pre-treatments, but the WCd was similar between the two desiccation pre-treatment groups (Table 3, Fig. 4; Table S5). Haemolymph osmolality followed a similar temporal variation pattern in all treatments, remaining stable during salinity exposure (χ2=4.5, P=0.106) but differed in magnitude among treatments (χ2=13.4, P=0.001). Specimens previously exposed to SD showed a higher osmotic concentration than those from the CD and RD pre-treatments (P<0.001 in Dunn's post hoc comparisons) (Fig. 5).

We found similar interactive effects of salinity and desiccation stressors in two species from main representative coleopteran families inhabiting saline inland waters. Exposure to stressful salinity had beneficial effects on the regulation of water balance under a subsequent desiccation stress. In contrast, a negative synergistic effect on performance was found when the order of exposure to the stressors was inverted. These results are clear evidence of the mechanistic links between tolerance to these co-occurring stressors in water beetles, which could have played a key role in the colonization of these systems and may have important implications in the context of climate change.

Physiological mechanisms linking tolerance to salinity and desiccation

In our first experiment, specimens of the two studied species showed higher body water content after exposure to the stressful salinities than those held at their respective optimum salinities, and also reduced water loss under the subsequent desiccation exposure in the case of N. baeticus. Such adjustments contributed to extended survival time in N. baeticus and had a smaller, but still positive, effect on the performance of E. jesusarribasi.

Although knowledge on the specific mechanisms of osmoregulation in aquatic beetles is still too poor to provide a mechanistic explanation for the observed cross-tolerance responses, different non-mutually exclusive processes might underlie this pattern. The specimens could have increased their drinking rates during exposure to hyperosmotic conditions in order to compensate for water loss by osmosis, a common behaviour in other saline aquatic insects such as mosquito larvae (e.g. Bradley and Phillips, 1975; Patrick and Bradley, 2000). This could account for the increase in water content and subsequent higher desiccation resistance observed here, but would also have the obvious collateral effect of ingestion of a substantial amount of salts from the medium. However, the species studied here have been shown to be able to osmoregulate over a wide hyperosmotic range including the stressful salinities tested here (Pallarés et al., 2015). In Coleoptera, excretion of salts and water reabsorption are mainly achieved through the activity of Malpighian tubules and specialized parts of the hindgut, such as the rectal pads (Crowson, 1981; Elliott and King, 1985; Machin and O'Donnell, 1991; Ramsay, 1964). The pre-activation of these osmoregulatory organs and tissues to maintain water and osmotic balance during salinity exposure probably contributed to minimize water loss during the subsequent desiccation exposure. The control of cuticle permeability is also one of the main mechanisms used to prevent water loss in terrestrial insects and has been shown to be a phenotypically plastic trait (e.g. Stinziano et al., 2015; Terblanche et al., 2010), although its role in aquatic insects has been less well explored (e.g. Alarie et al., 1998; Jacob and Hanssen, 1986). In the case of the species studied here, which use plastron or air bubbles for aquatic respiration, the cuticle surface exposed to water is reduced during water immersion. In consequence, modulation of cuticle permeability could make a smaller contribution to the control of water loss in the aquatic than in the aerial medium and therefore would have a relatively minor influence in the cross-tolerance between salinity and desiccation. Further studies on the relative contributions of these different mechanisms to the maintenance of water and ionic balance would improve understanding of the physiological basis of the cross-tolerance pattern found here.

In the second experiment, contrary to our expectations, exposure to either slow desiccation at a moderate relative humidity or rapid extreme desiccation reduced performance under a subsequent salinity stress in the two species. However, it cannot be discarded that less severe conditions than those tested here could elicit cross-tolerance to salinity, especially in E. jesusarribasi, which showed a short-term survival improvement after rapid desiccation exposure.

Homeostasis of the extracellular fluid is highly plastic in insects; different types of dehydration (fast versus slow) can induce different homeostatic processes (Beyenbach, 2016) and molecular responses (e.g. Lopez-Martinez et al., 2009). This seems to be the case in E. jesusarribasi, which showed clearly different responses after the two desiccation pre-treatments tested here in relation to the WLR and total water loss. After a slow but intense desiccation, this species maintained its haemolymph osmolality close to the values measured under control conditions and around the typical osmotic concentration of osmoregulatory insects, i.e. 300–400 mOsmol kg−1 (Bradley, 2009). Such an ability to display strict osmotic regulation under extreme conditions of dehydration has been observed in desert beetles (e.g. Naidu and Hattingh, 1988; Naidu, 2001) and Drosophila (Albers and Bradley, 2004). Nevertheless, osmolality increased rapidly when beetles were transferred to the hyperosmotic medium (second stress), probably as a consequence of a large intake of saline water to compensate for the large quantity of water previously lost and a disruption of osmoregulatory mechanisms. In contrast, when E. jesusarribasi was subjected to rapid dehydration, specimens showed a high haemolymph osmolality, suggesting that the osmotic concentration could have been sacrificed in this case in order to preserve extracellular and intracellular fluid volume under desiccation (Beyenbach, 2016). Although this water conservation strategy is apparently less energetically costly than active osmoregulation (Evans, 2008; Peña-Villalobos et al., 2016), maintaining high haemolymph concentrations was detrimental under the subsequent salinity stress in the long term for this species. These differences in osmoregulatory responses between the different desiccation conditions were not so evident in N. baeticus, probably because this species experienced similar water loss under the two desiccation pre-treatments, resulting in a similar decline of performance under the subsequent salinity stress.

Taken together, the results of our two experiments show that the effects on fitness of the combination of salinity and desiccation differ drastically depending on the order of the stress sequence (i.e. which stress precedes the other) as well as on their intensity and duration. Plastic osmoregulatory and water balance responses have both costs and benefits, and these are determined by the time scale and magnitude of variation in environmental conditions (e.g. Kleynhans et al., 2014; Todgham et al., 2005).

Ecological and evolutionary implications of linked salinity and desiccation tolerance

To our knowledge, only one case of cross-tolerance between salinity and desiccation has previously been reported in insects, specifically in larvae of the Antarctic midge Belgica antarctica (Elnitsky et al., 2009). However, there is evidence of cross-tolerance between salinity and thermal tolerance in diverse saline-tolerant taxa in inland waters, as for the branchiopod Daphnia pulex (Chen and Stillman, 2012), the water boatmen Tricocorixa verticalis and Sigara lateralis (Coccia et al., 2013), as well as several water beetle species, including those studied here (Arribas et al., 2012b; Botella-Cruz et al., 2016; Sánchez-Fernández et al., 2010). The finding of these patterns in ecologically similar species from different lineages highlights the ecological relevance of the cross-tolerance phenomenon (Hochachka and Somero, 2002; Kültz, 2005) and the central role of adaption to salinity and co-occurring stressors in arid inland waters.

During drought events in saline water, drying is often preceded by an increase in salinity levels. Acclimatization to such increasing salinities might allow insect populations showing cross-tolerance to enhance their resistance to the high temperatures and dehydration stress that they face during dispersal to wet refuges. This would also imply that adults from generations that emerge and develop in different seasons (i.e. spring–summer versus autumn–winter) could show different desiccation and thermal resistance in relation to the seasonal salinity levels of their habitat (e.g. Kalra and Parkash, 2016).

Considering the predicted intensification of droughts and water salinization across the Mediterranean region (Bonada and Resh, 2013; Filipe et al., 2013; Sala et al., 2000), cross-tolerance to salinity and desiccation could provide a significant physiological advantage for saline species to deal with such changes over related freshwater species in the same climatic area. However, as the combined effects of these stressors greatly depend on the intensity and relative timing of each stressor and also considering that most of the fauna in saline inland waters already inhabit conditions that are close to their physiological limits (Arribas et al., 2012a), persistent droughts may strongly limit the potential for salinity acclimation of these endemic species (Arribas et al., 2012b; Sánchez-Fernández et al., 2010), compromising their persistence in their current localities.

Salinity and desiccation play essential roles in the distribution and diversification of aquatic lineages. Tolerance to these stressors could have co-evolved in water beetle lineages as an exaptation process (Arribas et al., 2014; Pallarés et al., 2016). However, the exaptation hypothesis assumes that salinity and desiccation tolerance are mechanistically linked, something which had not previously been demonstrated. The cross-tolerance found here provides a solid, experimentally based trace of a potential parallel evolution of these traits in water beetles, offering a new frame to interpret diversification in inland waters. Under this scenario, global aridification events, which are broadly recognized as drivers of diversification in aquatic taxa (e.g. Pinceel et al., 2013; Dorn et al., 2014), could have been particularly relevant in the case of saline environments as one of the main forces for the colonization and further diversification in these systems (Arribas et al., 2014). Further studies on the temporal sequence of evolution of desiccation resistance and salinity tolerance in aquatic lineages could provide important insights for understanding the role of these mechanisms in driving evolution in inland waters.

We thank Almudena Gutiérrez and technicians from the SAI (Servicio de Apoyo a la Investigación de la Universidad de Murcia) for their technical support.

Author contributions

All authors conceived the idea and revised the manuscript; J.V., P.A., M.B.-C. and S.P. designed the experiments; J.V., A.M., M.B.-C. and S.P. conducted specimen collection in field; S.P. and M.B.-C. performed the experiments and statistical analyses; S.P. wrote the manuscript and the final version was revised and approved by all authors.

Funding

This research was funded by the I+D+i project CGL2013-48950-C2-2-P [Ministerio de Economía y Competitividad cofinanced with FEDER (Fonds Européen de Développement Économique et Régional) funds], a PhD grant from Universidad de Murcia to S.P. and two post-doctoral grants from the Royal Society UK (Newton International Program) and from Ministerio de Economía y Competitividad (Juan de la Cierva Formación Program) to P.A.

Alarie
,
Y.
,
Joly
,
H.
and
Dennie
,
D.
(
1998
).
Cuticular hydrocarbon analysis of the aquatic beetle Agabus anthracinus Mannerheim (Coleoptera: Dytiscidae)
.
Can. Entomol.
130
,
615
-
629
.
Albers
,
M. A.
and
Bradley
,
T. J.
(
2004
).
Osmotic regulation in adult Drosophila melanogaster during dehydration and rehydration
.
J. Exp. Biol.
207
,
2313
-
2321
.
Altman
,
D. G.
(
1992
).
Analysis of survival times
. In
Practical Statistics for Medical Research
(ed.
D. G.
Altman
), pp.
365
-
393
.
London
:
Chapman and Hall
.
Arribas
,
P.
,
Abellán
,
P.
,
Velasco
,
J.
,
Bilton
,
D. T.
,
Millán
,
A.
and
Sánchez-Fernández
,
D.
(
2012a
).
Evaluating drivers of vulnerability to climate change: a guide for insect conservation strategies
.
Glob. Change Biol.
18
,
2135
-
2146
.
Arribas
,
P.
,
Velasco
,
J.
,
Abellán
,
P.
,
Sánchez-Fernández
,
D.
,
Andújar
,
C.
,
Calosi
,
P.
,
Millán
,
A.
,
Ribera
,
I.
and
Bilton
,
D. T.
(
2012b
).
Dispersal ability rather than ecological tolerance drives differences in range size between lentic and lotic water beetles (Coleoptera: Hydrophilidae)
.
J. Biogeogr.
39
,
984
-
994
.
Arribas
,
P.
,
Andújar
,
C.
,
Sánchez-Fernández
,
D.
,
Abellán
,
P.
and
Millán
,
A.
(
2013
).
Integrative taxonomy and conservation of cryptic beetles in the Mediterranean region (Hydrophilidae)
.
Zool. Scr.
42
,
182
-
200
.
Arribas
,
P.
,
Andújar
,
C.
,
Abellán
,
P.
,
Velasco
,
J.
,
Millán
,
A.
and
Ribera
,
I.
(
2014
).
Tempo and mode of the multiple origins of salinity tolerance in a water beetle lineage
.
Mol. Ecol.
23
,
360
-
373
.
Bazinet
,
A. L.
,
Marshall
,
K. E.
,
MacMillan
,
H. A.
,
Williams
,
C. M.
and
Sinclair
,
B. J.
(
2010
).
Rapid changes in desiccation resistance in Drosophila melanogaster are facilitated by changes in cuticular permeability
.
J. Insect Physiol.
56
,
2006
-
2012
.
Beyenbach
,
K. W.
(
2016
).
The plasticity of extracellular fluid homeostasis in insects
.
J. Exp. Biol.
219
,
2596
-
2607
.
Bilton
,
D. T.
(
2014
).
Dispersal in Dytiscidae
. In
Ecology, Systematics, and Natural History of Predaceous Diving Beetles (Coleoptera: Dytiscidae)
(ed.
D. A.
Yee
), pp.
387
-
407
.
New York
:
Springer
.
Bonada
,
N.
and
Resh
,
V. H.
(
2013
).
Mediterranean-climate streams and rivers: geographically separated but ecologically comparable freshwater systems
.
Hydrobiologia
719
,
1
-
29
.
Botella-Cruz
,
M.
,
Carbonell
,
J. A.
,
Pallarés
,
S.
,
Millán
,
A.
and
Velasco
,
J.
(
2016
).
Plasticity of thermal limits in the aquatic saline beetle Enochrus politus (Küster 1849). under changing environmental conditions
.
Limnetica
35
,
131
-
142
.
Bradley
,
T.
(
2009
).
Animal Osmoregulation
.
New York
:
Oxford University Press
.
Bradley
,
T. J.
and
Phillips
,
J. E.
(
1975
).
The secretion of hyperosmotic fluid by the rectum of a saline-water mosquito larva, Aedes taeniorhynchus
.
J. Exp. Biol.
63
,
331
-
342
.
Bubliy
,
O. A.
,
Kristensen
,
T. N.
,
Kellermann
,
V.
and
Loeschcke
,
V.
(
2012
).
Plastic responses to four environmental stresses and cross-resistance in a laboratory population of Drosophila melanogaster
.
Funct. Ecol.
26
,
245
-
253
.
Céspedes
,
V.
,
Pallarés
,
S.
,
Arribas
,
P.
,
Millán
,
A.
and
Velasco
,
J.
(
2013
).
Water beetle tolerance to salinity and anionic composition and its relationship to habitat occupancy
.
J. Insect Physiol.
59
,
1076
-
1084
.
Chen
,
X.
and
Stillman
,
J. H.
(
2012
).
Multigenerational analysis of temperature and salinity variability affects on metabolic rate, generation time, and acute thermal and salinity tolerance in Daphnia pulex
.
J. Therm. Biol.
37
,
185
-
194
.
Chown
,
S. L.
and
Nicolson
,
S.
(
2004
).
Insect Physiological Ecology: Mechanisms and patterns
.
Oxford
:
Oxford University Press
.
Coccia
,
C.
,
Calosi
,
P.
,
Boyero
,
L.
,
Green
,
A. J.
and
Bilton
,
D. T.
(
2013
).
Does ecophysiology determine invasion success? A comparison between the invasive boatman Trichocorixa verticalis verticalis and the native Sigara lateralis (Hemiptera, Corixidae) in South-West Spain
.
PLoS ONE
8
,
e63105
.
Cohen
,
E.
(
2012
).
Roles of aquaporins in osmoregulation, desiccation and cold hardiness in insects
.
Entomol. Ornithol. Herpetol.
S1
,
1
-
17
.
Crowson
,
R. A.
(
1981
).
The biology of the Coleoptera
.
London
:
Academic Press
.
DeBiasse
,
M. B.
and
Kelly
,
M. W.
(
2016
).
Plastic and evolved responses to global change: what can we learn from comparative transcriptomics?
J. Hered.
107
,
71
-
81
.
Dorn
,
A.
,
Musilová
,
Z.
,
Platzer
,
M.
,
Reichwald
,
K.
and
Cellerino
,
A.
(
2014
).
The strange case of East African annual fishes: aridification correlates with diversification for a savannah aquatic group?
BMC Evol. Biol.
14
,
210
.
Dow
,
J. A. T.
and
Davies
,
S. A.
(
2006
).
The Malpighian tubule: rapid insights from post–genomic biology
.
J. Insect Physiol.
52
,
365
-
378
.
Dudley
,
R.
(
2000
).
The Biomechanics of Insect Flight
.
Princenton
:
Princenton University Press
.
Elliott
,
P.
and
King
,
P. E.
(
1985
).
A comparison of rectal pad structure in some terrestrial and intertidal Carabid beetles (Insecta)
.
J. Nat. Hist.
19
,
1209
-
1216
.
Elnitsky
,
M. A.
,
Benoit
,
J. B.
,
Lopez-Martinez
,
G.
,
Denlinger
,
D. L.
and
Lee
,
R. E.
Jr
. (
2009
).
Osmoregulation and salinity tolerance in the Antarctic midge, Belgica antarctica: seawater exposure confers enhanced tolerance to freezing and dehydration
.
J. Exp. Biol.
212
,
2864
-
2871
.
Evans
,
D. H.
(
2008
).
Osmotic and Ionic Regulation: Cells and Animals
.
Boca Ratón
:
CRC Press
.
Filipe
,
A. F.
,
Lawrence
,
J. E.
and
Bonada
,
N.
(
2013
).
Vulnerability of stream biota to climate change in mediterranean climate regions: a synthesis of ecological responses and conservation challenges
.
Hydrobiologia
719
,
331
-
351
.
Folguera
,
G.
,
Bastías
,
D. A.
,
Caers
,
J.
,
Rojas
,
J. M.
,
Piulachs
,
M.-D.
,
Bellés
,
X.
and
Bozinovic
,
F.
(
2011
).
An experimental test of the role of environmental temperature variability on ectotherm molecular, physiological and life–history traits: implications for global warming
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
159
,
242
-
246
.
Gasith
,
A.
and
Resh
,
V. H.
(
1999
).
Streams in Mediterranean climate regions: abiotic influences and biotic responses to predictable seasonal events
.
Annu. Rev. Ecol. Syst.
30
,
51
-
81
.
Gibbs
,
A. G.
and
Rajpurohit
,
S.
(
2010
).
Cuticular lipids and water balance
. In
Insect Hydrocarbons: Biology, Biochemistry, and Chemical Ecology
(ed.
G. J.
Blomquist
and
A. G.
Bagnères
), pp.
100
-
120
.
Cambridge
:
Cambridge University Press
.
Gunderson
,
A. R.
,
Armstrong
,
E. J.
and
Stillman
,
J. H.
(
2016
).
Multiple stressors in a changing world: the need for an improved perspective on physiological responses to the dynamic marine environment
.
Annu. Rev. Mar. Sci.
8
,
357
-
378
.
Hershkovitz
,
Y.
and
Gasith
,
A.
(
2013
).
Resistance, resilience, and community dynamics in mediterranean-climate streams
.
Hydrobiologia
719
,
59
-
75
.
Hewitt
,
J. E.
,
Ellis
,
J. I.
and
Thrush
,
S. F.
(
2016
).
Multiple stressors, nonlinear effects and the implications of climate change impacts on marine coastal ecosystems
.
Glob. Change Biol.
22
,
2665
-
2675
.
Hochachka
,
P. W.
and
Somero
,
G. N.
(
2002
).
Biochemical Adaptation
.
New York
:
Oxford University Press
.
Holmstrup
,
M.
,
Hedlund
,
K.
and
Boriss
,
H.
(
2002
).
Drought acclimation and lipid composition in Folsomia candida: implications for cold shock, heat shock and acute desiccation stress
.
J. Insect Physiol.
48
,
961
-
970
.
Hunt
,
T.
,
Bergsten
,
J.
,
Levkanicova
,
Z.
,
Papadopoulou
,
A.
,
John
,
O. S.
,
Wild
,
R.
,
Hammond
,
P. M.
,
Ahrens
,
D.
,
Balke
,
M.
,
Caterino
,
M. S.
, et al. 
(
2007
).
A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation
.
Science
318
,
1913
.
Huth
,
T. J.
and
Place
,
S. P.
(
2016
).
Transcriptome wide analyses reveal a sustained cellular stress response in the gill tissue of Trematomus bernacchii after acclimation to multiple stressors
.
BMC Genomics
17
,
453
.
Jäch
,
M. A.
and
Balke
,
M.
(
2008
).
Global diversity of water beetles (Coleoptera) in freshwater
.
Hydrobiologia
595
,
419
-
442
.
Jacob
,
J.
and
Hanssen
,
H.-P.
(
1986
).
Distribution and variability of cuticular hydrocarbons within the Coleoptera
.
Biochem. Syst. Ecol.
14
,
207
-
210
.
Kalra
,
B.
and
Parkash
,
R.
(
2016
).
Effects of saturation deficit on desiccation resistance and water balance in seasonal populations of the tropical drosophilid Zaprionus indianus
.
J. Exp. Biol.
219
,
3237
-
3245
.
Kefford
,
B. J.
,
Hickey
,
G. L.
,
Gasith
,
A.
,
Ben-David
,
E.
,
Dunlop
,
J. E.
,
Palmer
,
C. G.
,
Allan
,
K.
,
Choy
,
S. C.
and
Piscart
,
C.
(
2012
).
Global scale variation in the salinity sensitivity of riverine macroinvertebrates: eastern Australia, France, Israel and South Africa
.
PLoS ONE
7
,
e35224
.
Kleynhans
,
E.
,
Clusella-Trullas
,
S.
and
Terblanche
,
J. S.
(
2014
).
Impacts of environmental variability on desiccation rate, plastic responses and population dynamics of Glossina pallidipes
.
J. Evol. Biol.
27
,
337
-
348
.
Kültz
,
D.
(
2005
).
Molecular and evolutionary basis of the cellular stress response
.
Annu. Rev. Physiol.
67
,
225
-
257
.
Kumlu
,
M.
,
Kumlu
,
M.
and
Turkmen
,
S.
(
2010
).
Combined effects of temperature and salinity on critical thermal minima of pacific white shrimp Litopenaeus vannamei (Crustacea: Penaeidae)
.
J. Therm. Biol.
35
,
302
-
304
.
Larsen
,
E. H.
,
Deaton
,
L. E.
,
Onken
,
H.
,
O'Donnell
,
M.
,
Grosell
,
M.
,
Dantzler
,
W. H.
and
Weihrauch
,
D.
(
2014
).
Osmoregulation and Excretion
.
Compr. Physiol.
4
,
405
-
573
.
Le Lagadec
,
M. D.
,
Chown
,
S. L.
and
Scholtz
,
C. H.
(
1998
).
Desiccation resistance and water balance in southern African keratin beetles (Coleoptera, Trogidae): the influence of body size and habitat
.
J. Comp. Physiol. B Biochem. Syst. Environ. Physiol.
168
,
112
-
122
.
Lopez-Martinez
,
G.
,
Benoit
,
J. B.
,
Rinehart
,
J. P.
,
Elnitsky
,
M. A.
,
Lee
,
R. E.
, Jr
and
Denlinger
,
D. L.
(
2009
).
Dehydration, rehydration, and overhydration alter patterns of gene expression in the Antarctic midge, Belgica antarctica
.
J. Comp. Physiol. B.
179
,
481
-
491
.
Lyons
,
C. L.
,
Coetzee
,
M.
,
Terblanche
,
J. S.
and
Chown
,
S. L.
(
2014
).
Desiccation tolerance as a function of age, sex, humidity and temperature in adults of the African malaria vectors Anopheles arabiensis and Anopheles funestus
.
J. Exp. Biol.
217
,
3823
-
3833
.
Machin
,
J.
and
O'Donnell
,
M. J.
(
1991
).
Rectal complex ion activities and electrochemical gradients in larvae of the desert beetle, Onymacris: comparisons with Tenebrio
.
J. Insect Physiol.
37
,
829
-
838
.
Millán
,
A.
,
Velasco
,
J.
,
Gutiérrez-Cánovas
,
C.
,
Arribas
,
P.
,
Picazo
,
F.
,
Sánchez-Fernández
,
D.
and
Abellán
,
P.
(
2011
).
Mediterranean saline streams in southeast Spain: what do we know?
J. Arid Environ.
75
,
1352
-
1359
.
Millán
,
A.
,
Sánchez–Fernández
,
D.
,
Abellán
,
P.
,
Picazo
,
F.
,
Carbonell
,
J. A.
,
Lobo
,
J. M.
and
Ribera
,
I.
(
2014
).
Atlas de los coleópteros acuáticos de España peninsular
.
Madrid
:
Ministerio de Agricultura, Alimentación y Medio Ambiente
.
Naidu
,
S. G.
(
2001
).
Water balance and osmoregulation in Stenocara gracilipes, a wax–blooming tenebrionid beetle from the Namib Desert
.
J. Insect Physiol.
47
,
1429
-
1440
.
Naidu
,
S. G.
and
Hattingh
,
J.
(
1988
).
Water balance and osmoregulation in Physadesmia globosa, a diurnal tenebrionid beetle from the Namib desert
.
J. Insect Physiol.
34
,
911
-
917
.
Pallarés
,
S.
,
Arribas
,
P.
,
Bilton
,
D. T.
,
Millán
,
A.
and
Velasco
,
J.
(
2015
).
The comparative osmoregulatory ability of two water beetle genera whose species span the fresh–hypersaline gradient in inland waters (Coleoptera: Dytiscidae, Hydrophilidae)
.
PLoS ONE
10
,
e0124299
.
Pallarés
,
S.
,
Velasco
,
J.
,
Millán
,
A.
,
Bilton
,
D. T.
and
Arribas
,
P.
(
2016
).
Aquatic insects dealing with dehydration: do desiccation resistance traits differ in species with contrasting habitat preferences?
PeerJ
4
,
e2382
.
Pansch
,
C.
,
Nasrolahi
,
A.
,
Appelhans
,
Y. S.
and
Wahl
,
M.
(
2012
).
Impacts of ocean warming and acidification on the larval development of the barnacle Amphibalanus improvisus
.
J. Exp. Mar. Biol. Ecol.
420-421
,
48
-
55
.
Patrick
,
M. L.
and
Bradley
,
T. J.
(
2000
).
The physiology of salinity tolerance in larvae of two species of Culex mosquitoes: the role of compatible solutes
.
J. Exp. Biol.
203
,
821
-
830
.
Peña-Villalobos
,
I.
,
Narváez
,
C.
and
Sabat
,
P.
(
2016
).
Metabolic cost of osmoregulation in a hypertonic environment in the invasive African clawed frog Xenopus laevis
.
Biol. Open
5
,
955
-
961
.
Pinceel
,
T.
,
Brendonck
,
L.
,
Larmuseau
,
M. H. D.
,
Vanhove
,
M. P. M.
,
Timms
,
B. V.
and
Vanschoenwinkel
,
B.
(
2013
).
Environmental change as a driver of diversification in temporary aquatic habitats: does the genetic structure of extant fairy shrimp populations reflect historic aridification?
Freshwater Biol.
58
,
1556
-
1572
.
Ramsay
,
J. A.
(
1964
).
The rectal complex of the mealworm Tenebrio molitor, L. (Coleoptera, Tenebrionidae)
.
Philos. Trans. R. Soc. Lond. B Biol. Sci.
248
,
279
-
314
.
Sala
,
O. E.
,
Chapin
,
F. S.
, III
,
Armesto
,
J. J.
,
Berlow
,
E.
,
Bloomfield
,
J.
,
Dirzo
,
R.
,
Huber-Sanwald
,
E.
,
Huenneke
,
L. F.
,
Jackson
,
R. B.
,
Kinzig
,
A.
, et al. 
(
2000
).
Global biodiversity scenarios for the year 2100
.
Science
287
,
1770
-
1774
.
Sánchez-Fernández
,
D.
,
Bilton
,
D. T.
,
Abellán
,
P.
,
Ribera
,
I.
,
Velasco
,
J.
and
Millán
,
A.
(
2008
).
Are the endemic water beetles of the Iberian Peninsula and the Balearic Islands effectively protected?
Biol. Conserv.
141
,
1612
-
1627
.
Sánchez-Fernández
,
D.
,
Calosi
,
P.
,
Atfield
,
A.
,
Arribas
,
P.
,
Velasco
,
J.
,
Spicer
,
J. I.
,
Millán
,
A.
and
Bilton
,
D. T.
(
2010
).
Reduced salinities compromise the thermal tolerance of hypersaline specialist diving beetles
.
Physiol. Entomol.
35
,
265
-
273
.
Sinclair
,
B. J.
,
Ferguson
,
L. V.
,
Salehipour-Shirazi
,
G.
and
MacMillan
,
H. A.
(
2013
).
Cross–tolerance and cross–talk in the cold: relating low temperatures to desiccation and immune stress in insects
.
Integr. Comp. Biol.
53
,
545
-
556
.
Stinziano
,
J. R.
,
Sové
,
R. J.
,
Rundle
,
H. D.
and
Sinclair
,
B. J.
(
2015
).
Rapid desiccation hardening changes the cuticular hydrocarbon profile of Drosophila melanogaster
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
180
,
38
-
42
.
Strachan
,
S. R.
,
Chester
,
E. T.
and
Robson
,
B. J.
(
2015
).
Freshwater invertebrate life history strategies for surviving desiccation
.
Springer Sci. Rev.
3
,
57
-
75
.
Stubbington
,
R.
,
Gunn
,
J.
,
Little
,
S.
,
Worrall
,
T. P.
and
Wood
,
P. J.
(
2016
).
Macroinvertebrate seedbank composition in relation to antecedent duration of drying and multiple wet–dry cycles in a temporary stream
.
Freshwater Biol.
61
,
1293
-
1307
.
Terblanche
,
J. S.
,
Clusella-Trullas
,
S.
and
Chown
,
S. L.
(
2010
).
Phenotypic plasticity of gas exchange pattern and water loss in Scarabaeus spretus (Coleoptera: Scarabaeidae): deconstructing the basis for metabolic rate variation
.
J. Exp. Biol.
213
,
2940
-
2949
.
Todgham
,
A. E.
and
Stillman
,
J. H.
(
2013
).
Physiological responses to shifts in multiple environmental stressors: relevance in a changing world
.
Integr. Comp. Biol.
53
,
539
-
544
.
Todgham
,
A. E.
,
Schulte
,
P. M.
and
Iwama
,
G. K.
(
2005
).
Cross–tolerance in the tidepool sculpin: the role of heat shock proteins
.
Physiol. Biochem. Zool.
78
,
133
-
144
.
Velasco
,
J.
,
Millán
,
A.
,
Hernández
,
J.
,
Gutiérrez
,
C.
,
Abellán
,
P.
,
Sánchez
,
D.
and
Ruiz
,
M.
(
2006
).
Response of biotic communities to salinity changes in a Mediterranean hypersaline stream
.
Saline Systems
2
,
1
-
15
.
Williams
,
C. M.
,
Poulin
,
R.
and
Sinclair
,
B. J.
(
2004
).
Increased haemolymph osmolality suggests a new route for behavioural manipulation of Talorchestia quoyana (Amphipoda: Talitridae). by its mermithid parasite
.
Funct. Ecol.
18
,
685
-
691
.
Zuur
,
A. F.
,
Ieno
,
E. N.
,
Walker
,
N. J.
,
Saveliev
,
A. A.
and
Smith
,
G. M.
(
2009
).
Mixed Effects Models and Extensions in Ecology with R
.
New York
:
Spring Science and Business Media
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information