TY - JOUR AU1 - Dupoué,, Andréaz AU2 - Sorlin,, Mahaut AU3 - Richard,, Murielle AU4 - Le Galliard, Jean, François AU5 - Lourdais,, Olivier AU6 - Clobert,, Jean AU7 - Aubret,, Fabien AB - Abstract Parent-offspring conflicts are widespread given that resources are often limited. Recent evidence has shown that availability of water can trigger such conflict during pregnancy in viviparous squamate species (lizards and snakes) and thus questions the role of water in the evolution of reproductive modes. Here, we examined the impact of water restriction during gravidity in the oviparous form of the bimodal common lizard (Zootoca vivipara), using a protocol previously used on the viviparous form. Females were captured in early gravidity from six populations along a 600 m altitudinal gradient to investigate whether environmental conditions (altitude, water access and temperature) exacerbate responses to water restriction. Females were significantly dehydrated after water restriction, irrespective of their reproductive status (gravid vs. non-reproductive), relative reproductive effort (relative clutch mass), and treatment timing (embryonic development stage). Female dehydration, together with reproductive performance, varied with altitude, probably due to long term acclimation or local adaptation. This moderate water-based intergenerational conflict in gravid females contrasts sharply with previous findings for the viviparous form, with implications to the evolutionary reversion from viviparity to oviparity. It is likely that oviparity constitutes a water-saving reproductive mode which might help mitigate intensive temperature-driven population extinctions at low altitudes. altitude, dehydration, ectotherm, mother-offspring conflicts, parity mode, reproduction INTRODUCTION Life history trade-offs, by shaping the allocation of resources within and between generations, are central to the evolution of reproductive strategies (Stearns, 1992; Harshman & Zera, 2007; Kölliker et al., 2015). Intergenerational trade-offs refer to parental allocation strategies to the offspring and may lead to parent-offspring conflicts (POCs) when resources become scarce (Trivers, 1974). POCs have been used to model the evolution of parental care, optimal parental investment, life history strategies (Godfray, 1995; Haig, 2010; Kölliker et al., 2015), and even the evolution of reproductive modes (Zeh & Zeh, 2000; Crespi & Semeniuk, 2004; Pollux et al., 2014). In animals, parity modes range from oviparity (egg-laying) to viviparity (live-bearing) and are unimodal in many animal families and orders (e.g. birds or eutherian mammals). Yet, squamate reptiles (i.e. lizards and snakes) display exceptional diversity in reproductive strategies by covering the whole spectrum of reproductive modes from some species laying eggs in early developmental stages to others bearing embryos with complex placental structures (Blackburn, 2006; Van Dyke et al., 2014). The different selective forces involved in the evolution of viviparity in squamates remain under debate, as shown by recent phylogenetic and comparative studies (Pyron & Burbrink, 2014; Blackburn, 2015; Shine, 2015). An accepted view is that viviparity evolved from oviparity more than 110 times independently (Blackburn, 2006) given the general selective advantage associated with extended maternal care (Shine, 2014). In ectotherms, it is often advantageous for the mother to control stable and optimal thermal conditions of development through behavioural means (Li et al., 2009; Lorioux et al., 2013a, b; Foucart et al., 2018). However, prolonged egg retention may generate an arms race between maternal investment and embryonic resource acquisition [the ‘Viviparity conflict hypothesis’ (Crespi & Semeniuk, 2004)]. Such elevated costs of pregnancy may preclude the transition from oviparity to viviparity and explain intermediate stages of retention as observed in many oviparous squamates (Andrews, 2004; Blackburn, 2015). Additionally, pregnancy costs may explain reverse transitions from viviparity to oviparity as recently hypothesized (Recknagel et al., 2018; Gao et al., 2019; Horreo et al., 2020). That said, how water shapes the evolution of reproductive strategies remains an open question, and our objective was to determine whether a mother-offspring conflict for water would differ between reproductive modes. The vast majority of reptile species are lecithotrophic, meaning that mothers invest nutrients into the yolk prior to ovulation (Fig. 1). Embryos then rely on this store of nutrients for their growth and development (Blackburn & Stewart, 2011). In lecithotrophic species, embryos cannot manipulate energy allocation after ovulation, so there is a release from the mother-offspring conflict for energy. However, mothers must also control thermal conditions during pregnancy and supply a substantial amount of water to their developing embryos (Packard, 1991). Maternal water supply is essential for embryos to convert vitellus into embryonic tissues and, in addition, water demand of embryos increases with exponential somatic growth (Packard, 1991; Shine & Thompson, 2006; Lourdais et al., 2015). Yet, in most squamate species from temperate regions, pregnancy occurs during summer to maximize opportunities for thermoregulation at a time when water is potentially scarce. Water restriction during pregnancy can therefore trigger intergenerational conflicts between mother and offspring in viviparous lizards and snakes (Dupoué et al., 2015a, 2018a). Remarkably, water-based mother-offspring conflicts have been associated with higher offspring mortality in a viviparous lizard (Dupoué et al., 2018a), suggesting that water-limiting environments may challenge the benefits of prolonged egg retention. Figure 1. Open in new tabDownload slide Representation of water restriction protocols in the oviparous form (this study) or the viviparous form (Dupoué et al., 2018a) of the bimodal lizard Z. vivipara, over a standard reproductive cycle (Foucart et al., 2014). In both reproductive modes, females were exposed in early June to the same experimental 14-day period of water restriction. The associated range of treatment timing (oviparous: blue arrow; viviparous: orange arrow, scaled in the figure) covered mid-late vitellogenesis to early gravidity in oviparous females (32 to 4 days before oviposition), while water restriction occurred relatively later in the viviparous form (35 to 3 days before parturition). We hypothesized water constraints to be relatively low in oviparous females given that most of embryo water demand occurs during somatic growth (Lourdais et al., 2015), once oviparous females have laid their eggs. Figure 1. Open in new tabDownload slide Representation of water restriction protocols in the oviparous form (this study) or the viviparous form (Dupoué et al., 2018a) of the bimodal lizard Z. vivipara, over a standard reproductive cycle (Foucart et al., 2014). In both reproductive modes, females were exposed in early June to the same experimental 14-day period of water restriction. The associated range of treatment timing (oviparous: blue arrow; viviparous: orange arrow, scaled in the figure) covered mid-late vitellogenesis to early gravidity in oviparous females (32 to 4 days before oviposition), while water restriction occurred relatively later in the viviparous form (35 to 3 days before parturition). We hypothesized water constraints to be relatively low in oviparous females given that most of embryo water demand occurs during somatic growth (Lourdais et al., 2015), once oviparous females have laid their eggs. Species with bimodal reproduction (i.e. distinct reproductive modes between populations) offer the ideal opportunity to clarify the evolution of reproductive strategies. For instance, the bimodal European common lizard (Z. vivipara) is arguably a model species to examine costs and benefits of parity modes (Foucart et al., 2014; Recknagel & Elmer, 2019), and to understand the factors leading to the transition to viviparity (Surget-Groba et al., 2006; Rodríguez-Díaz & Braña, 2012). This lizard is oviparous in the southern margin of distribution range (Pyrenean Mountains, Northern Spain and locally in the Alps and Balkans) but viviparous in the rest of its Eurasian distribution (Heulin et al., 2000). The two reproductive forms share similar affinities for relatively wet habitats, suggesting high water dependence for reproduction (Lorenzon et al., 1999; Marquis et al., 2008; Le Galliard et al., 2012). Two independent oviparous clades have been identified (Recknagel et al., 2018): while oviparity is an ancestral trait in the Eastern oviparous populations (Z. vivipara carniolica) it likely results from reversal in Western oviparous ones (Z. vivipara louislantzi). Females lay slightly calcified eggs at a relatively advanced stage in embryonic development [stage 30–35 sensuDufaure & Hubert (1961)], and embryo stage at oviposition can vary with altitude (Heulin et al., 1997; Rodríguez-Díaz & Braña, 2012) and reproductive effort (Foucart et al., 2017). In the present study, we tested if water-based conflict between mother and offspring occurred in the Western oviparous form of the common lizard (Z. vivipara louislantzi), as previously found in the viviparous form and using the same protocol (Dupoué et al., 2018a). We compared the physiological responses (dehydration rate) of gravid and non-reproductive oviparous females exposed to a 14-day period of water restriction. Treatment exposure occurred in early June, relatively soon in female reproductive cycle compared to the viviparous females (Fig. 1). We characterized the two assumptions of water-based conflict by determining: i) if gravid females paid an extra water cost compared to non-reproductive ones after water restriction, and ii) whether female dehydration correlated to relative clutch mass (RCM, the residuals extracted from an independent linear mixed model of clutch mass against female body size) (Supporting Information, Fig. S1) (Dupoué et al., 2015a, 2018a). We also examined the consequences of water restriction on reproductive performance including oviposition date, RCM, incubation time, hatching success and hatching traits (offspring size and mass). Water restriction may (Dauphin-Villemant & Xavier, 1986; Brusch et al., 2018) or may not (Dupoué et al., 2015a, 2018a) impact reproductive traits such as clutch mass at laying (oviparous species) or litter mass at parturition (viviparous species). Furthermore, a competition for water between developing siblings may exist (Bonnet et al., 2017). These intergenerational and intrauterine conflicts may explain why hatchlings of oviparous species are generally larger than neonates of viviparous ones [see Table 1 in Bonnet et al. (2017)]. Indeed after oviposition, the physiological constraints on female water balance are released (Fig. 1).We therefore hypothesized mother-offspring conflict for water to be modest in the oviparous common lizard due to relatively low hydric investment toward the eggs. We predicted similar responses to water restriction between gravid and non-reproductive females, and that dehydration of gravid females should also be uncorrelated with RCM or treatment timing due to lower water demand for reproduction. Females were captured in the Pyrenean Mountains from six natural populations, and we further examined the effect of population altitude, water availability (permanent access in peat-bogs vs. periodic in underwood and dry meadows) and average thermal conditions. We hypothesized that females would be locally adapted to the environmental conditions associated with altitude, as seen in a previous comparative study of thermal preferences in the same geographic area (Trochet et al., 2018). In response to water restriction, we expected gravid females from low altitudes (i.e. relatively hot and dry habitats), to exhibit stronger resistance to dehydration and greater reproductive performance compared with those from high altitudes. MATERIAL AND METHODS Study species, populations and husbandry The European common lizard (Z. vivipara) is a small (adult snout-vent length (SVL) ~50–75 mm), widespread lacertid typically found in cold humid peat bogs and heathland habitats from Western Europe to Scandinavia and Eastern Russia (Heulin et al., 2000). Between the 15th and 31st of May 2018, we caught a total of 134 adult females (105 gravid and 29 non-reproductive) from six populations distributed throughout the Pyrenees mountain range (Supporting Information, Table S1). The reproductive strategy of the Western oviparous form substantially changes with increasing SVL across altitude (Supporting Information, Fig. S2), and earlier and multiple clutches as well as earlier age at maturation observed in warmer habitats at lowland elevations (< 300 m) compared to single reproductive event at higher altitudes (Heulin et al., 1997). In our study, the altitudinal range (990–1580 m) was strong enough to examine altitudinal variation of reproductive performance (see below), while narrow enough to avoid strong differences in life history strategies. Females were captured by hand and then transferred to the laboratory where they were housed in individual terraria (18 x 12 x 12 cm) with sterilized soil, a shelter, and basking heat until parturition. Each individual was provided a 20–35 °C thermal gradient for 6 h per day (09:00-12:00 and 14:00-17:00) using a 25 W incandescent light bulb placed over one end of each terrarium. They had ad libitum access to water in a petri dish and we further provided water three times per day at 09:00, 13:00 and 17:00. We further provided them with three mealworms (Tenebrio molitor) every 2 days. Females were kept in these standard conditions until oviposition except during the water restriction experiment (see Experimental Design). We characterized environmental conditions for each population, with altitude (Supporting Information, Table S1), presence of permanent vs. periodic water sources and air temperature, because these measures were likely to shape local adaptations in the regulation of the water balance as documented previously in the viviparous form (Dupoué et al., 2017a). We recorded air temperature using three data loggers (iButtons, Maxim Integrated Products, Sunnyvale, CA, USA, ± 0.5 °C) per population placed at locations where we found most lizards within vegetation at ground level and completely shaded to avoid the effect of radiation. Air temperature was recorded every hour, and we standardized the sampling period from the June 30th to July 25th 2018 to enable population comparisons (Dupoué et al., 2017a). During this sampling period, we extracted the average daily minimum and maximum temperatures (Tmin and Tmax, respectively) to assess the thermal conditions of each population (Supporting Information, Table S1). Experimental design A few days following capture, we randomly assigned females within each population to two experimental treatments following the exact same protocol we previously used on viviparous females from the same species (Dupoué et al., 2018a). In the water-restricted treatment, we removed the water bowl and reduced the misting frequency to once per day occurring in the morning. In the control treatment, lizards had permanent access to the water bowl and were misted three times per day. The control treatment mimics conditions in which lizards find permanent access to water (e.g. peat bog, marsh). Instead, the water-restricted treatment reflects summer conditions in dry habitats where in the absence of precipitation, morning dew is the only source of drinking water. Water restriction lasted for 14 days and occurred on a range of treatment timing from 32 to 4 days before oviposition, a time when females were between mid-late vitellogenesis to early gravidity (Fig. 1). After the period of water restriction, all females returned to the control water conditions, having permanent access to water in a water bowl and being misted three times per day. We released non-reproductive females within 2 weeks following experiments after controlling their body mass (BM) trajectories and palpation to confirm their non-reproductive status. Gravid females were released within 3 days post-laying at their capture location. After oviposition, clutches (range = 1 to 9 eggs) were weighed (± 1 mg) and placed in individual plastic cups on water-saturated vermiculite to maintain hydric conditions and incubated at Tset = 25 °C to optimise incubation time without risk of overheating for embryo development (Rodríguez-Díaz et al., 2010; Foucart et al., 2018). We used three incubators (Novital Covatutto Eco, Italy) and we randomly distributed the clutches in the three incubators (N = 35 per incubator) that remained within a 1 °C range of observed temperatures (mean ± SD, incubator 1: Tobs = 24.38 ± 0.26 °C; incubator 2: Tobs = 24.05 ± 0.39 °C; incubator 3: Tobs = 24.63 ± 0.30 °C). Juveniles were released within 3 days post-hatching in their respective populations. Female water balance Females were weighed (BM, ± 1 mg) every 3 days throughout the water restriction period to assess effects of water restriction on short-term changes since BM is an indicator of hydration state (Lillywhite et al., 2012). However, for simplicity, we only examined BM changes (ΔBM = BMfinal - BMinitial) from the initiation to the end of the water restriction period. We also measured absolute changes in plasma osmolality (ΔOsmo = Osmofinal - Osmoinitial), which is a rigorous measure of whole-body hydration in vertebrates (Peterson, 2002). For these assessments, females were bled at the onset and at the end of the water restriction period using a standard protocol (Meylan et al., 2003). Blood samples (c. 40 µL whole blood) were collected from the post-orbital sinus and centrifuged for 5 min at 11 000 rpm. Plasma was then separated from the blood cells and kept frozen at -30 °C in airtight tubes until analyses were performed. Plasma osmolality (± 1 mOsm.kg-1) was determined using a vapour pressure osmometer (Model 5500, Wescor, Logan, UT, USA) according to the protocol previously established in this species (Dupoué et al., 2017a). Before analyses, plasma was diluted (1:1) in a physiological serum (304 mOsm.kg-1) so that plasma osmolality could be determined from 10 µL duplicates (intra-individual coefficient of variation: 1.17%). Reproductive performance We examined the effects of water restriction on different reproductive traits. We checked females daily to compare oviposition date. For each female we subtracted the oviposition date from the last day in hydric treatment to calculate the treatment timing. As an index of reproductive effort, we estimated size-adjusted relative clutch mass (RCM) at oviposition (i.e. residuals from the linear relationship between clutch mass and body size, CM—SVL: t103 = 7.2, P < 0.001) (Supporting Information, Fig. S1) to determine how much females deviate from the reproductive effort predicted by their body size (Bonnet et al., 2003). We used SVL as the denominator in this linear regression instead of post-laying BM because: i) SVL is fixed and less biased than post-laying BM (Bonnet et al., 2003), ii) SVL-adjustments might better explain reproductive changes than BM-adjustments (Dupoué & Lourdais, 2014; Foucart et al., 2014), and iii) because SVL was a better predictor of female clutch mass (r2 = 0.34) than post-laying BM (r2 = 0.21). We calculated the incubation time as the day difference between hatching and oviposition dates. We determined hatching success as the number of live juveniles within the clutch against stillborn and undeveloped eggs. All alive offspring were then counted, weighed (BM, ± 1 mg), measured (SVL, ± 1 mm), and sexed by counting ventral scales on the medioventral lines (Lecomte et al., 1992). This method assumes a discriminant relationship between the number of ventral scales and phenotypic sex, based on sexual dimorphism in body size at birth, which has been recently shown to be independent from water restriction (Dupoué et al., 2019). Z. vivipara from oviparous populations differ in body size, shape and life history from viviparous populations (A. Dupoué, pers. obs.). Thus, we used the adult females sampled for this study and males with obvious secondary sexual characters (N = 79) captured only to count ventral scales (released the same day) to fit a discriminant function on ventral scales and sex (mean ± SE, left side, females: 29.61 ± 0.09, males: 26.16 ± 0.13, χ 1,220 = 31.8, P < 0.001; right side, females: 29.52 ± 0.10, males: 26.39 ± 0.14; χ 1,221 = 189.8, P < 0.001) with high determination success (96.5%) (Supporting Information, Table S2). Scalation may vary across altitude (Thorpe & Baez, 1993), but here scale number differences between sexes was independent of geographic locality (interaction term between population and ventral scales, left side: χ 3,168 = 5.3, P = 0.153, right side: χ 3,171 = 0.9, P = 0.835). Statistical analyses All analyses were performed using R software (version 3.2.0, R Core Team 2016, https://www.r-project.org/). First, we checked whether initial BM (BMini) and osmolality (Osmoini) differed between gravid and non-reproductive females according to treatment affiliation, using linear mixed models [package lme4 (Bates et al., 2015)]. Models included fixed effects of hydric treatment (control vs. water restriction), reproductive state (gravid vs. non-reproductive) and their first-order interaction. We set the population origin as a random effect to control for non-independence of females within populations. ΔBM and ΔOsmo were then analysed using similar model construction, including the fixed effects of hydric treatment, reproductive state, their first-order interaction, and the random effects of population. We added the effects of initial physiological value (at the onset of experiment) as linear covariate. In gravid females, we analysed the relationships between dehydration rate (ΔBM or ΔOsmo) and RCM using similar design except that we replaced the reproductive status by RCM as described in first models. Preliminary analyses showed that female SVL had no effect on dehydration rate, which was confirmed in further analyses (Supporting Information, Tables S3-S4). We examined whether female dehydration rate could be further related to treatment timing and environmental conditions. For each response variable (ΔBM or ΔOsmo), we used the Akaike information criterion corrected for small sample size [AICc, package AICcmodavg (Mazerolle, 2019)], to compare a set of models with different environmental measures (Supporting Information, Table S3: all females; Supporting Information, Table S4: gravid females only). In all models, population was set as a random factor to account for the non-independence of females within a population. Our set of models included: i) a null model with only the intercept and random factor, ii) simple models (initial values of BM, SVL, osmolality, reproductive status or hydric treatment alone), iii) additive models with initial value, SVL, reproductive status, hydric treatment and additive effects of each environmental conditions (altitude, water access, and Tmin and Tmax treated separately; see all models in Supporting Information, Table S3), and iv) interactive models with initial value, SVL, reproductive status, hydric treatment, first- and second-order interactive effects between hydric treatment, reproductive status and each environmental condition treated separately (altitude, water access, Tmin and Tmax) (see Supporting Information, Table S3). In gravid females specifically, we used a similar procedure to test the effects of initial physiological state, SVL, treatment timing, hydric treatment, environmental conditions and their interaction on ΔBM and ΔOsmo (see all models in Supporting Information, Table S4). In some models, altitude was treated as explanatory covariates in addition with SVL or treatment timing, a procedure that may induce multicollinearity given the positive correlation between those (respectively: t6.1 = 3.1, P = 0.022, t106.0 = 6.9, P < 0.001). Despite this relationship, the Pearson’s correlation coefficients between SVL and altitude (r = 0.32) or treatment timing and altitude (r = 0.55) were below the threshold when multicollinearity bias may occur (r = 0.70) (Dormann et al., 2013). We used similar approaches to determine whether reproductive performance was correlated with treatment timing and whether it differed between hydric treatments, alone, in addition to, or in interaction with environmental conditions (altitude, water access, and Tmin and Tmax treated separately; see all models in Supporting information, Table S5). In all models, population was set as a random factor. Oviposition date (Julian date), RCM and incubation time were analysed with linear mixed models, and hatching success was analysed with mixed-effect logistic regressions including a logit link and binomial error term (number of viable vs. number of failed offspring). Offspring SVL and BM were analysed with similar linear mixed models, with offspring sex as an additive or interactive factor (see all models in Supporting information, Table S5) and random factors included both the mother identity and population to account for non-independence between siblings and within populations. RESULTS Physiological responses to water restriction At the onset of experiment, BMini was higher in gravid than non-gravid females (t135.6 = -2.6, P = 0.010), irrespective of treatment affiliation (control: 3.200 ± 0.087, water restriction 3.192 ± 0.085 g, t132.9 = -0.5, P = 0.632). Osmoini did not differ either between reproductive statuses (t135.0 = -0.5, P = 0.635) or between treatment affiliation (control: 319.1 ± 2.2, water restriction 319.4 ± 2.7 mOsm.kg-1, t129.4 = 0.2, P = 0.845). At the end of exposure to hydric treatment, changes in body mass (ΔBM) were positively correlated with BMini (t134.0 = 4.3, P < 0.001) and lower in water-restricted females compared to controls (t134.0 = -5.5, P < 0.001) (Fig. 2A), irrespective of their SVL and reproductive state (Supporting information, Table S3). In gravid females, ΔBM was positively correlated with RCM (t105.0 = 3.1, P = 0.003) (Fig. 2B), and negatively correlated to treatment timing (t105.0 = -2.1, P = 0.041) (Fig. 2C), irrespective of hydric treatment (Supporting Information, Table S4). At the same time, the changes in plasma osmolality (ΔOsmo) were negatively correlated with Osmoini (t133.2 = -8.3, P < 0.001) and significantly increased following water restriction compared to control conditions (t129.2 = 6.6, P < 0.001) (Fig. 2D) irrespective of their SVL and reproductive state (Supporting Information, Table S3). In gravid females, ΔOsmo was neither correlated with RCM (t104.9 = -0.8, P = 0.410) (Fig. 2E), nor with treatment timing (t105.0 = -0.1, P = 0.957) (Fig. 2F), in both hydric treatments (Supporting Information, Table S4). Figure 2. Open in new tabDownload slide Examination of physiological responses to water restriction in gravid and non-reproductive female oviparous common lizards (N = 134). In response to a 2-week period of water restriction (filled circles), (A) females experienced greater loss in body mass (ΔBM) compared to control conditions (open circles), independent of their reproductive state. In addition, (B) ΔBM were positively correlated with RCM, and (C) negatively correlated with treatment timing (number of days between the end of hydric treatment and oviposition). Although water-restricted females were more dehydrated (higher plasma osmolality, ΔOsmo) compared to controls, this was independent of (D) reproductive state, (E) RCM, and (F) treatment timing. Significant differences between control and water-restricted females are symbolized: (***) = P < 0.001. Figure 2. Open in new tabDownload slide Examination of physiological responses to water restriction in gravid and non-reproductive female oviparous common lizards (N = 134). In response to a 2-week period of water restriction (filled circles), (A) females experienced greater loss in body mass (ΔBM) compared to control conditions (open circles), independent of their reproductive state. In addition, (B) ΔBM were positively correlated with RCM, and (C) negatively correlated with treatment timing (number of days between the end of hydric treatment and oviposition). Although water-restricted females were more dehydrated (higher plasma osmolality, ΔOsmo) compared to controls, this was independent of (D) reproductive state, (E) RCM, and (F) treatment timing. Significant differences between control and water-restricted females are symbolized: (***) = P < 0.001. Females from lower altitude populations exhibited higher sensitivity to water restriction since they experienced a greater loss of BM (t105.0 = 2.9, P = 0.005) (Fig. 3A) and higher increase in plasma osmolality (t105.0 = -3.2, P = 0.002) (Fig. 3B) than those from highlands, whereas no significant altitudinal variation of hydration state occurred in control females (both, P > 0.800) (Fig. 3A, B), irrespective of reproductive status (Supporting information, Table S3). Figure 3. Open in new tabDownload slide Relationships between physiological deviation of water balance in gravid females (N = 105) and altitude of the population. After a 2-week period of water restriction (filled circles, solid lines), females from the lower altitudes experienced greater dehydration, whereas no-relation between changes in hydration state and altitude occurred in control females (open circles). This is illustrated by (A) higher loss of BM and (B) increased plasma osmolality in water-restricted females from lower altitudes. Data are represented by mean ± SE of body mass changes (Δ BM) and osmolality changes (Δ Osmo) within population. Figure 3. Open in new tabDownload slide Relationships between physiological deviation of water balance in gravid females (N = 105) and altitude of the population. After a 2-week period of water restriction (filled circles, solid lines), females from the lower altitudes experienced greater dehydration, whereas no-relation between changes in hydration state and altitude occurred in control females (open circles). This is illustrated by (A) higher loss of BM and (B) increased plasma osmolality in water-restricted females from lower altitudes. Data are represented by mean ± SE of body mass changes (Δ BM) and osmolality changes (Δ Osmo) within population. Effects of water restriction on reproductive output Water restriction had no effect on reproductive performance, including oviposition date (mean ± SE, control mothers: 29th June ± 1 day, water-restricted mothers: 30th June ± 1 day, t100.3 = 0.8, P = 0.415), clutch mass (control mothers: 1.366 ± 0.059 g, water-restricted mothers: 1.370 ± 0.065 g, t102.7 = 0.0, P = 0.967), RCM (control mothers: 0.005 ± 0.044, water-restricted mothers: -0.005 ± 0.055, t102.6 = -0.1, P = 0.891), incubation time (control mothers: 19.0 ± 0.2 days, water-restricted mothers: 19.1 ± 0.2 days, t73.0 = 0.7, P = 0.479), hatching success (control mothers: 52.3 ± 5.4%, water-restricted mothers: 59.3 ± 5.8%, z = 1.7, P = 0.097), or on offspring SVL (control mothers: 20.07 ± 0.11 mm, water-restricted mothers: 20.19 ± 0.10 mm, t72.5 = 0.5, P = 0.632) or BM (control mothers: 225.1 ± 2.7 mg, water-restricted mothers: 229.2 ± 2.7 mg, t69.6 = 0.5, P = 0.615). Female SVL did not influence oviposition date, RCM, hatching success, and offspring SVL and BM (Supporting information, Table S5); however, it was negatively correlated with incubation time (t78.0 = -2.1, P = 0.037). Effects of environmental conditions on reproductive output Regarding the effects of environmental conditions, we found that altitude and treatment timing significantly explained variation in reproductive performance (Supporting information, Table S5). Gravid females from lower altitude populations laid eggs earlier than those from higher altitude (t106.0 = 8.2, P < 0.001) (Fig. 4A), while the incubation time decreased with altitude (t78.0 = -4.3, P < 0.001) (Fig. 4B). RCM was negatively correlated with treatment timing (t106.0 = -5.7, P < 0.001), so that females laying eggs sooner also had greater reproductive effort than those with late oviposition (Fig. 4C). Hatching success was negatively correlated with treatment timing (z = -3.0, P = 0.002) (Fig. 4D) and interactively impacted by mother treatment and altitude (z = -4.1, P < 0.001), since it decreased with altitude in clutches from water-restricted mothers (z = -3.2, P = 0.002) (Fig. 4E), while slightly increasing in those from control mothers (z = 2.1, P = 0.033) (Fig. 4E). Both morphometric measures of offspring (SVL and BM) were positively correlated with altitude (offspring SVL: t7.4 = 4.0, P = 0.005; offspring BM: t7.1 = 3.6, P = 0.008), irrespective of mother treatment (Supporting information, Table S5). In addition, female offspring had longer SVL than males (t289.3 = 6.0, P < 0.001) (Fig. 4F), and sex differences in BM depended upon mother treatment (interaction term: t269.5 = -3.1, P = 0.002). That is, daughters of control mothers had lower BM than their brothers (t268.6 = -2.9, P = 0.004) (Fig. 4G), whereas there was no difference between offspring of water-restricted mothers (t269.3 = 1.5, P = 0.137) (Fig. 4G). Figure 4. Open in new tabDownload slide Variation of reproductive performance in water-restricted (filled circles, solid lines) and control (open circles, dashed lines) females (N = 105) and their offspring (N = 343). (A) oviposition date was positively correlated with altitude, while altitude negatively correlated with (B) incubation time. Treatment timing (number of days between the end of hydric treatment and oviposition) negatively correlated with (C) relative clutch mass (RCM) and (D) hatching success. (E) hatching success was also negatively correlated with altitude in water-restricted females (solid line) but not in control ones (dashed line). In offspring that survived, (F) female offspring SVL was longer than males, and (G) sex differences in BM between daughter and sons depended on mother treatment. Data are represented by mean ± SE of reproductive outputs within population and significant relationships are symbolised by trend lines (predictions of final models). Significant effects of interaction terms between offspring sex and mother treatments are symbolized: {n.s.} = non-significant and {**} = P < 0.01. Figure 4. Open in new tabDownload slide Variation of reproductive performance in water-restricted (filled circles, solid lines) and control (open circles, dashed lines) females (N = 105) and their offspring (N = 343). (A) oviposition date was positively correlated with altitude, while altitude negatively correlated with (B) incubation time. Treatment timing (number of days between the end of hydric treatment and oviposition) negatively correlated with (C) relative clutch mass (RCM) and (D) hatching success. (E) hatching success was also negatively correlated with altitude in water-restricted females (solid line) but not in control ones (dashed line). In offspring that survived, (F) female offspring SVL was longer than males, and (G) sex differences in BM between daughter and sons depended on mother treatment. Data are represented by mean ± SE of reproductive outputs within population and significant relationships are symbolised by trend lines (predictions of final models). Significant effects of interaction terms between offspring sex and mother treatments are symbolized: {n.s.} = non-significant and {**} = P < 0.01. DISCUSSION Following a 2-week water restriction period, both gravid and non-reproductive females responded the same way, suggesting no additional hydric cost of reproduction [sensuLourdais et al. (2017)]. This was further confirmed by an absence of a relationship between dehydration rate and both relative reproductive effort and treatment timing. These results probably reflect the low water requirements of embryos at early developmental stages, and contrast with recent findings in an oviparous snake, where water restriction applied during the entire gravidity period (3 weeks) did affect females’ hydration state and resulted in lower egg mass (Brusch et al., 2018). Interestingly, we found that independently of reproductive state, female physiological responses to water restriction were shaped by altitude. We had expected habitat water access or temperature to better reflect phenotypic variation in water balance regulation (Guillon et al., 2014; Cox & Cox, 2015; Dupoué et al., 2017a). We assessed microclimatic conditions within a restricted spatiotemporal window (i.e. sensor in shade under vegetation over 3 summer weeks), following previous methodology (Dupoué et al., 2017a). Due to logistical constraints, we were unable to measure other ecologically relevant components of thermal (e.g. operative temperature or thermal heterogeneity) or hydric (e.g. deficit in water vapour pressure) environments as potential determinants of local constraints for a heliothermic ectotherm. Alternatively, altitude is a fixed environmental measure that integrates many factors from macro- to microclimatic conditions, vegetation type, and snow cover, and may therefore represent a better descriptor of annual conditions than punctual measures of temperatures. Here, females from lower altitudes faced higher dehydration rates probably because they were exposed to water restriction when embryos were more developed and with higher water needs (Lourdais et al., 2015). In addition, females at high altitudes may be locally adapted and less permeable to water loss due to body shape [e.g. lower surface-to-volume ratio (Dupoué et al., 2015b)], or more prompt to initiate water-saving strategies (e.g. lower thermal preferences), to limit evaporative water loss and dehydration risk (Köhler et al., 2011; Rozen‐Rechels et al., 2019). In support of this last hypothesis, we recently documented thermal preference negatively correlated with altitude in females from the exact same populations (Trochet et al., 2018). Altitude also correlated with oviposition date (sooner at lower altitudes) and incubation duration (faster at higher altitudes). This is consistent with expectations and previous evidence on phenology variation along elevation gradients (Heulin et al., 1997; Rodríguez-Díaz & Braña, 2012; Rutschmann et al., 2016). RCM and hatching success increased with treatment timing irrespective of hydric treatment, thus suggesting that reproduction follows a “sooner is better” pattern as found in another lizard species (Le Henanff et al., 2013). Once corrected for treatment timing, more unexpected was the interactive impact of water restriction and altitude on hatching success. In control conditions, hatching success increased by 6% in high altitude populations whereas it increased by a strong 40% in lower altitude populations when females faced higher dehydration. This result is surprising, because dehydration during gravidity or incubation is usually associated with lower reproductive success (Packard, 1991). Due to the orographic effects of altitude, females from low altitudes regularly face periods of more severe water restriction (summer droughts) compared to those from high altitudes with more frequent rain episodes. This contrast might have resulted in some degree of local adaptation with enhanced reproductive fitness for low altitude females in the drier conditions (water restriction treatment) and in high altitude females in the wetter conditions (control treatment). Additionally, greater hatching success in water-restricted females from low altitudes might have indirectly resulted from locally adapted thermoregulation. As stated before, females from low altitude populations have higher preferred body temperatures (Trochet et al., 2018). In this species, small differences in thermoregulation during gravidity may strongly impact hatching success (Foucart et al., 2018). Surviving offspring from higher altitudes had larger body size and body mass than those from lower altitudes. Overall, daughters were longer than sons as classically documented in this species (Lecomte et al., 1992; Le Galliard et al., 2006; Dupoué et al., 2019), and sex differences in mass depended upon mother treatment. That is, daughters from control mothers had lower body mass than sons, whereas offspring from water-restricted mothers had similar mass. Altogether, our results therefore suggest multi-level trade-offs since females from low altitudes laid eggs sooner and despite prolonged egg incubation time in standard conditions, they had greater hatching success but produced smaller offspring. These trade-offs were relatively independent of punctual dehydration but they were shaped by long-term acclimation or local adaptation to altitude. In the Western viviparous lineages of the common lizard, population extinction risk is increasing at low altitudes where lizards are exposed to higher temperatures (Massot et al., 2008; Sinervo et al., 2010) and their associated costs (Dupoué et al., 2017b, 2018b). Viviparous species are generally predicted to be more vulnerable to climate change than oviparous ones (Sinervo et al., 2010, 2018). In support of this hypothesis, the oviparous common lizard remains present (in low abundance) at sea level and only in those humid landscapes (peat bogs, forest marshes) that support its presence (Berroneau, 2014). Although causal factors remain elusive, we showed here that water restriction early in gravidity had relatively low physiological impact on the oviparous form and this contrasts with previous findings in the viviparous form (Dupoué et al., 2018a). In fact, our results suggest that the oviparous form at low altitudes might even gain in fitness (hatching success) when exposed to water restriction. Future work is now critically needed to experimentally manipulate both temperature and water to measure the adaptive significance and long-term response in this bimodal species. To conclude, the evolution of reproductive strategies constitutes the core part of life-history theory (Stearns, 1992). Given their outstanding variability and flexibility in reproductive modes, squamates remain key models to investigate the causes and consequences of transitioning to viviparity (Blackburn, 2006; Laird et al., 2019). Here, we repeated an experimental procedure in both forms of a reproductively bimodal species. Although some fitness consequences may also appear over time, the lack of changes in females physiology depending on reproductive effort suggest that they did not pay immediate water costs contrary to their viviparous relative (Dupoué et al., 2018a). This implies that water demand associated with late pregnancy stages might represent a barrier to prolonged egg retention when evolving in water limiting environments. Additionally, water constraints may also favour a reverse transition from viviparity to oviparity, if the hydric costs of reproduction become too high (Lourdais et al., 2017; Dupoué et al., 2019a). In support of this innovative hypothesis, phylogenomic analyses unravelled that the Western lineages populations sampled here likely originated from a viviparous ancestor and evolved back to an oviparous reproductive mode (Recknagel et al., 2018; Horreo et al., 2020). Hence, egg-laying might have evolved during the last Pleistocene and Holocene as an adaptation to high habitat aridity. Our study therefore calls for future work to test this hypothesis and include water constraints when studying the environmental drivers of reproductive strategies. SUPPORTING INFORMATION Additional Supporting Information may be found in the online version of this article at the publisher's web-site: FigureS1. Positive relationship between reproductive effort (clutch mass) and female body size (SVL). Figure S2. Positive relationship between female body size (SVL) and altitude (t6.4 = 3.1, P = 0.022). Table S1. Localizations, altitude, climatic conditions and water access (permanent in peatbog type habitat, periodic in dry meadow) in the six natural populations of the oviparous form of common lizard (Z. vivipara) from the Pyrenees Mountain range. Table S2. Sex determination table using scale counting along the medioventral (left and right side) as previously described (Lecomte et al., 1992). Table S3. AICc based model selection in all females (gravid and non-reproductive) comparing a null model (intercept only) to models testing relationships between female dehydration indexes (changes in body mass ΔBM and changes in plasma osmolality ΔOsmo) and reproductive status, hydric treatment and environmental conditions. Table S4. AICc based model selection in gravid females only comparing a null model (intercept only) to models testing relationships between female dehydration indexes (changes in body mass ΔBM and changes in plasma osmolality ΔOsmo) and environmental conditions. Table S5. AICc based model selection comparing a null model (intercept only) to models testing relationships between female reproductive performance and environmental conditions. ACKNOWLEDGMENTS We thank Pierre-Louis Duarte and Loïse Bessaguet for their help in the field and Pauline Blaimont for English revising and helping during experiments. We are also grateful to the ‘Office Nationale des Forêts’, the ‘Parc National des Pyrénées’ and the region Occitanie for allowing us to sample lizards. We thank Prof. Daniel G. Blackburn and three anonymous reviewers for greatly improving the early drafts of this manuscript. This study received the financial and technical support from the Centre National de la Recherche Scientifique (CNRS), and was funded by the Interreg POCTEFA under the project ‘Ectopyr’ (EFA031/15 to FA), by the Agence Nationale de la Recherche under the ‘Aquatherm’ project (ANR-17-CE02-0013 to JFLG), and by the Region Nouvelle Aquitaine under the ‘Aquastress’ project (2018-1R20214 to OL). AD was financially supported by the Australian Research Council (DP180100058 to A/Prof Anne Peters) at the time of writing. AD conceived the study, captured individuals, led the experiments, performed the analyses, and drafted the manuscript. MS and MR helped in capture and performed the experiments. All authors contributed critically to the conception, the redaction and gave final approval for publication. This study was performed in accordance with laws relative to capture, transport and experiments on Z. vivipara (DREAL Occitanie permit #2017-s-02) and further approved by an independent ethical committee (Apafis #14651-2018041215394197 v5). The authors declare no conflicts of interest. REFERENCES Andrews RM . 2004 . Patterns of embryonic development. In: Deeming DC, ed. Reptilian incubation: environment, evolution and behaviour . Nottingham : Nottingham University Press , 75 – 102 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Bates D , Maechler M, Bolker B, Walker S. 2015 . Fitting linear mixed-effects models using {lme4} . Journal of Statistical Software 67 : 1 – 48 . Google Scholar Crossref Search ADS WorldCat Berroneau M . 2014 . Lézard vivipare Zootoca vivipara (Lichtenstein, 1823). In: Nature C, ed. Atlas des amphibiens et reptiles d’Aquitaine. Le Haillan, France: Association Cistude Nature, 106 – 109 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Blackburn DG . 2006 . Squamate reptiles as model organisms for the evolution of viviparity . Herpetological Monographs 20 : 131 – 146 . Google Scholar Crossref Search ADS WorldCat Blackburn DG . 2015 . Evolution of viviparity in squamate reptiles: reversibility reconsidered . Journal of Experimental Zoology Part B: Molecular and Developmental Evolution 324 : 473 – 486 . Google Scholar Crossref Search ADS WorldCat Blackburn DG , Stewart JR. 2011 . Viviparity and placentation in snakes. In: Aldrich RD, Sever DM, eds. Reproductive biology and phylogeny of snakes . Enfield, New Hampshire : Science Publishers , 119 – 181 . Google Scholar Crossref Search ADS Google Scholar Google Preview WorldCat COPAC Bonnet X , Naulleau G, Shine R. 2017 . The evolutionary economics of embryonic-sac fluids in squamate reptiles . The American Naturalist 189 : 333 – 344 . Google Scholar Crossref Search ADS PubMed WorldCat Bonnet X , Shine R, Lourdais O, Naulleau G. 2003 . Measures of reproductive allometry are sensitive to sampling bias . Functional Ecology 17 : 39 – 49 . Google Scholar Crossref Search ADS WorldCat Brusch GA IV , Lourdais O, Kaminsky B, DeNardo DF. 2018 . Muscles provide an internal water reserve for reproduction . Proceedings of the Royal Society B: Biological Sciences 285 : 20180752 . Google Scholar Crossref Search ADS WorldCat Cox CL , Cox RM. 2015 . Evolutionary shifts in habitat aridity predict evaporative water loss across squamate reptiles . Evolution 69 : 2507 – 2516 . Google Scholar Crossref Search ADS PubMed WorldCat Crespi B , Semeniuk C. 2004 . Parent-offspring conflict in the evolution of vertebrate reproductive mode . The American Naturalist 163 : 635 – 653 . Google Scholar Crossref Search ADS PubMed WorldCat Dauphin-Villemant C , Xavier F. 1986 . Adrenal activity in the females Lacerta vivipara Jacquin: possible involvement in the success of gestation. In: Assemacher I, Boissin J, eds. Endocrine regulation as adaptive mechanism to environment . Paris : CNRS , 241 – 250 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Dormann CF , Elith J, Bacher S, Carré GCG, García Márquez JR, Gruber B, Lafourcade B, Leitao PJ, Münkemüller T, McClean CJ, Osborne PE, Reneking B, Schröder B, Skidmore AK, Zurell D, Lautenbach S. 2013 . Collinearity: a review of methods to deal with it and a simulation study evaluating their performance . Ecography 36 : 027 – 046 . Google Scholar Crossref Search ADS WorldCat Dufaure JP , Hubert J. 1961 . Table de développement du lézard vivipare: Lacerta (Zootoca) vivipara . Archives d’Anatomie Microscopique et de Morphologie Experimentale 50 : 309 – 328 . OpenURL Placeholder Text WorldCat Dupoué A , Blaimont P, Rozen‐Rechels D, Richard M, Meylan S, Clobert J, Miles DB, Martin R, Decencière B, Agostini S, Le Galliard J‐F. 2020 . Water availability and temperature induce changes in oxidative status during pregnancy in a viviparous lizard . Functional Ecology . 34 : 475 – 485 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Lourdais O, Meylan S, Brischoux F, Angelier F, Rozen‐Rechels D, Marcangeli Y, Decencière B, Agostini S, Le Galliard J‐F. 2019 . Some like it dry: Water restriction overrides heterogametic sex determination in two reptiles . Ecology and Evolution 9 : 6524 – 6533 . Google Scholar PubMed OpenURL Placeholder Text WorldCat Dupoué A , Le Galliard J-F, Josserand R, DeNardo DF, Decencière B, Agostini S, Haussy C, Meylan S. 2018a . Water restriction causes an intergenerational trade-off and delayed mother-offspring conflict in a viviparous lizard . Functional Ecology 32 : 676 – 686 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Rutschmann A, Le Galliard J-F, Clobert J, Blaimont P, Sinervo B, Miles DB, Haussy C, Meylan S. 2018b . Reduction of baseline corticosterone secretion correlates with climate warming and drying across wild lizard populations . Journal of Animal Ecology 87 : 1331 – 1341 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Rutschmann A, Le Galliard J-F, Miles DB, Clobert J, DeNardo DF, Brusch GA IV, Meylan S. 2017a . Water availability and environmental temperature correlate with geographic variation in water balance in common lizards . Oecologia 185 : 561 – 571 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Rutschmann A, Le Galliard J-F, Clobert J, Angelier F, Marciau C, Ruault S, Miles D, Meylan S. 2017b . Shorter telomeres precede population extinction in wild lizards . Scientific Reports 7 : 16976 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Brischoux F, Angelier F, DeNardo DF, Wright CD, Lourdais O. 2015a . Intergenerational trade-off for water may induce a mother-offspring conflict in favour of embryos in a viviparous snake . Functional Ecology 29 : 414 – 422 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Stahlschmidt ZR, Michaud B, Lourdais O. 2015b . Physiological state influences evaporative water loss and microclimate preference in the snake Vipera aspis . Physiology & Behavior 144 : 82 – 89 . Google Scholar Crossref Search ADS WorldCat Dupoué A , Lourdais O. 2014 . Relative reproductive effort drives metabolic changes and maternal emaciation during pregnancy in a viviparous snake . Journal of Zoology 293 : 49 – 56 . Google Scholar Crossref Search ADS WorldCat Foucart T , Heulin B, Lourdais O. 2017 . Clutch size influences embryonic stages at oviposition in a lizard with prolonged egg retention . Amphibia Reptilia 38 : 557 – 561 . Google Scholar Crossref Search ADS WorldCat Foucart T , Heulin B, Lourdais O. 2018 . Small changes, big benefits: testing the significance of maternal thermoregulation in a lizard with extended egg retention . Biological Journal of the Linnean Society 125 : 280 – 291 . Google Scholar Crossref Search ADS WorldCat Foucart T , Lourdais O, DeNardo DF, Heulin B. 2014 . Influence of reproductive mode on metabolic costs of reproduction: insight from the bimodal lizard Zootoca vivipara . Journal of Experimental Biology 217 : 4049 – 4056 . Google Scholar Crossref Search ADS PubMed WorldCat Le Galliard J-F , Massot M, Baron J-P, Clobert C. 2012 . Ecological effects of climate change on European reptiles. In: Brodie J, Post E, Doak D, eds. Wildlife conservation in a changing climate . Chicago : University of Chicago Press , 179 – 203 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Le Galliard J-F , Massot M, Landys MM, Meylan S, Clobert J. 2006 . Ontogenic sources of variation in sexual size dimorphism in a viviparous lizard . Journal of Evolutionary Biology 19 : 690 – 704 . Google Scholar Crossref Search ADS PubMed WorldCat Gao W , Sun YB, Zhou WW, Xiong Z-J, Chen L, Li H, Fu T-T, Xu K, Xu W, Ma L, Chen Y-J, Xiang X-Y, Zhou L, Zeng T, Zhang S, Jin J-Q, Chen H-M, Zhang G, Hillis DM, Ji X, Zhang Y-P, Che J. 2019 . Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity . Proceedings of the National Academy of Sciences of the United States of America 116 : 3646 – 3655 . Google Scholar Crossref Search ADS PubMed WorldCat Godfray HCJ . 1995 . Evolutionary theory of parent-offspring conflict . Nature 376 : 133 – 138 . Google Scholar Crossref Search ADS PubMed WorldCat Guillon M , Guiller G, DeNardo DF, Lourdais O. 2014 . Microclimate preferences correlate with contrasted evaporative water loss in parapatric vipers at their contact zone . Canadian Journal of Zoology 92 : 81 – 86 . Google Scholar Crossref Search ADS WorldCat Haig D . 2010 . Transfers and transitions: parent-offspring conflict, genomic imprinting, and the evolution of human life history . Proceedings of the National Academy of Sciences of the United States of America 107 : 1731 – 1735 . Google Scholar Crossref Search ADS PubMed WorldCat Harshman LG , Zera AJ. 2007 . The cost of reproduction: the devil in the details . Trends in Ecology & Evolution 22 : 80 – 86 . Google Scholar Crossref Search ADS PubMed WorldCat Le Henanff M , Meylan S, Lourdais O. 2013 . The sooner the better: reproductive phenology drives ontogenetic trajectories in a temperate squamate (Podarcis muralis) . Biological Journal of the Linnean Society 108 : 384 – 395 . Google Scholar Crossref Search ADS WorldCat Heulin B , Guillaume CP, Vogrin N, Surget-Groba Y, Tadic Z. 2000 . Further evidence of the existence of oviparous populations of Lacerta (Zootoca) vivipara in the NW of the Balkan Peninsula . Comptes Rendus de l’Académie des Sciences. Séries III, Sciences de la Vie 323 : 461 – 468 . OpenURL Placeholder Text WorldCat Heulin B , Osenegg-Leconte K, Michel D. 1997 . Demography of a bimodal reproductive species of lizard (Lacerta vivipara): survival and density characteristics of oviparous populations . Herpetologica 53 : 432 – 444 . OpenURL Placeholder Text WorldCat Horreo JL , Suarez T, Fitze PS. 2020 . Reversals in complex traits uncovered as reticulation events: Lessons from the evolution of parity-mode, chromosome morphology, and maternal resource transfer . Journal of Experimental Zoology Part B: Molecular and Developmental Evolution 334 : 5 – 13 . Google Scholar Crossref Search ADS WorldCat Köhler A , Sadowska J, Olszewska J, . 2011 . Staying warm or moist? Operative temperature and thermal preferences of common frogs (Rana temporaria) and effects on locomotion . Herpetological Journal 21 : 17 – 26 . OpenURL Placeholder Text WorldCat Kölliker M , Boos S, Wong JWY, Röllin L, Stucki D, Raveh S, Wu M, Meunier J. 2015 . Parent-offspring conflict and the genetic trade-offs shaping parental investment . Nature Communications 6 : 6850 . Google Scholar Crossref Search ADS PubMed WorldCat Laird MK , Thompson MB, Whittington CM. 2019 . Facultative oviparity in a viviparous skink (Saiphos equalis) . Biology Letters 15 : 20180827 . Google Scholar Crossref Search ADS PubMed WorldCat Lecomte J , Clobert J, Massot M. 1992 . Sex identification in juveniles of Lacerta vivipara . Amphibia-Reptilia 13 : 21 – 25 . Google Scholar Crossref Search ADS WorldCat Li H , Qu Y-F, Hu R-B, Ji X. 2009 . Evolution of viviparity in cold-climate lizards: testing the maternal manipulation hypothesis . Evolutionary Ecology 23 : 777 – 790 . Google Scholar Crossref Search ADS WorldCat Lillywhite HB , Brischoux F, Sheehy CM, Pfaller JB. 2012 . Dehydration and drinking responses in a pelagic sea snake . Integrative and Comparative Biology 52 : 227 – 234 . Google Scholar Crossref Search ADS PubMed WorldCat Lorenzon P , Clobert J, Oppliger A, John-Alder H. 1999 . Effect of water constraint on growth rate, activity and body temperature of yearling common lizard (Lacerta vivipara) . Oecologia 118 : 423 – 430 . Google Scholar Crossref Search ADS PubMed WorldCat Lorioux S , Lisse H, Lourdais O. 2013a . Dedicated mothers: predation risk and physical burden do not alter thermoregulatory behaviour of pregnant vipers . Animal Behaviour 86 : 401 – 408 . Google Scholar Crossref Search ADS WorldCat Lorioux S , Vaugoyeau M, DeNardo DF, Clobert J, Guillon M, Lourdais O. 2013b . Stage dependence of phenotypical and phenological maternal effects: insight into squamate reptile reproductive strategies . The American Naturalist 182 : 223 – 233 . Google Scholar Crossref Search ADS WorldCat Lourdais O , Dupoué A, Guillon M, Guiller G, Michaud B, DeNardo DF. 2017 . Hydric ‘costs’ of reproduction: pregnancy increases evaporative water loss in the snake Vipera aspis . Physiological and Biochemical Zoology 90 : 663 – 672 . Google Scholar Crossref Search ADS PubMed WorldCat Lourdais O , Lorioux S, Dupoué A, Wright C, DeNardo DF. 2015 . Embryonic water uptake during pregnancy is stage- and fecundity-dependent in the snake Vipera aspis . Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 189 : 102 – 106 . Google Scholar Crossref Search ADS WorldCat Marquis O , Massot M, Le Galliard J-F. 2008 . Intergenerational effects of climate generate cohort variation in lizard reproductive performance . Ecology 89 : 2575 – 2583 . Google Scholar Crossref Search ADS PubMed WorldCat Massot M , Clobert J, Ferrière R. 2008 . Climate warming, dispersal inhibition and extinction risk . Global Change Biology 14 : 461 – 469 . Google Scholar Crossref Search ADS WorldCat Mazerolle MJ . 2019 . AICcmodavg: Model selection and multimodel inference based on (Q)AIC(c) . R package version 2.2-2. https://cran.r-project.org/package=AICcmodavg. Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Meylan S , Dufty AM, Clobert J. 2003 . The effect of transdermal corticosterone application on plasma corticosterone levels in pregnant Lacerta vivipara . Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 134 : 497 – 503 . Google Scholar Crossref Search ADS WorldCat Packard GC . 1991 . The physiological and ecological importance of water to embryos of oviparous reptiles. In: Deeming DC, Ferguson MWJ, eds. Egg incubation: its effect on embryonic development in birds and reptiles . Cambridge : Cambridge University Press , 213 – 228 . Google Scholar Crossref Search ADS Google Scholar Google Preview WorldCat COPAC Peterson CC . 2002 . Temporal, population, and sexual variation in hematocrit of free-living desert tortoises: correlational tests of causal hypotheses . Canadian Journal of Zoology 80 : 461 – 470 . Google Scholar Crossref Search ADS WorldCat Pollux BJA , Meredith RW, Springer MS, Garland T, Reznick DN. 2014 . The evolution of the placenta drives a shift in sexual selection in livebearing fish . Nature 513 : 233 – 236 . Google Scholar Crossref Search ADS PubMed WorldCat Pyron RA , Burbrink FT. 2014 . Early origin of viviparity and multiple reversions to oviparity in squamate reptiles . Ecology Letters 17 : 13 – 21 . Google Scholar Crossref Search ADS PubMed WorldCat R Development Core Team (2016). R: A language and environment for statistical computing . Vienna, Austria : R Foundation for Statistical Computing . Recknagel H , Elmer KR. 2019 . Differential reproductive investment in co-occurring oviparous and viviparous common lizards (Zootoca vivipara) and implications for life-history trade-offs with viviparity . Oecologia 190 : 85 – 98 . Google Scholar Crossref Search ADS PubMed WorldCat Recknagel H , Kamenos NA, Elmer KR. 2018 . Common lizards break Dollo’s law of irreversibility: Genome-wide phylogenomics support a single origin of viviparity and re-evolution of oviparity . Molecular Phylogenetics and Evolution 127 : 579 – 588 . Google Scholar Crossref Search ADS PubMed WorldCat Rodríguez-Díaz T , Braña F. 2012 . Altitudinal variation in egg retention and rates of embryonic development in oviparous Zootoca vivipara fits predictions from the cold-climate model on the evolution of viviparity . Journal of Evolutionary Biology 25 : 1877 – 1887 . Google Scholar Crossref Search ADS PubMed WorldCat Rodríguez-Díaz T , González F, Ji X, Braña F. 2010 . Effects of incubation temperature on hatchling phenotypes in an oviparous lizard with prolonged egg retention: are the two main hypotheses on the evolution of viviparity compatible? Zoology 113 : 33 – 38 . Google Scholar Crossref Search ADS PubMed WorldCat Rozen‐Rechels D , Dupoué A, Meylan S, Qitout K, Decencière B, Agostini S, Le Galliard J-F. 2020 . Acclimation to water restriction implies different paces for behavioral and physiological responses in a lizard species . Physiological and Biochemical Zoology 93 : 160 – 174 . Google Scholar Crossref Search ADS PubMed WorldCat Rutschmann A , Miles DB, Le Galliard J-F, Richard M, Moulherat S, Sinervo B, Clobert J. 2016 . Climate and habitat interact to shape the thermal reaction norms of breeding phenology across lizard populations . Journal of Animal Ecology 85 : 457 – 466 . Google Scholar Crossref Search ADS PubMed WorldCat Shine R . 2014 . Evolution of an evolutionary hypothesis: a history of changing ideas about the adaptive significance of viviparity in reptiles . Journal of Herpetology 48 : 147 – 161 . Google Scholar Crossref Search ADS WorldCat Shine R . 2015 . The evolution of oviparity in squamate reptiles: An adaptationist perspective . Journal of Experimental Zoology Part B: Molecular and Developmental Evolution 324 : 487 – 492 . Google Scholar Crossref Search ADS WorldCat Shine R , Thompson MB. 2006 . Did embryonic responses to incubation conditions drive the evolution of reproductive modes in squamate reptiles? Herpetological Monographs 20 : 159 – 171 . Google Scholar Crossref Search ADS WorldCat Sinervo B , Méndez-de-la-Cruz F, Miles DB, Heulin B, Bastiaans E, Cruz MVS, Lara-Resendiz R, Martínez-Méndez N, Calderón-Espinosa ML, Meza-Lázaro RN, Gadsden H, Avila LJ, Morando M, De la Riva IJ, Sepulveda PV, Rocha CFD, Ibargüengoytía N, Puntriano CA, Massot M, Lepetz V, Oksanen TA, Chapple DG, Bauer AM, Branch WR, Clobert J, Sites JW Jr. 2010 . Erosion of lizard diversity by climate change and altered thermal niches . Science 328 : 894 – 899 . Google Scholar Crossref Search ADS PubMed WorldCat Sinervo B , Miles DB, Wu Y, Méndez-de-la-Cruz F, Kirchhof S, Qi Y. 2018 . Climate change, thermal niches, extinction risk and maternal-effect rescue of toad-headed lizards, Phrynocephalus, in thermal extremes of the Arabian Peninsula to the Qinghai-Tibetan Plateau . Integrative Zoology 13 : 450 – 470 . Google Scholar Crossref Search ADS PubMed WorldCat Stearns SC . 1992 . The evolution of life histories . Oxford : Oxford University Press . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Surget-Groba Y , Heulin B, Guillaume CP, Puky M, Semenov D, Orlova V, Kupriyanova L, Ghira I, Smajda B. 2006 . Multiple origins of viviparity, or reversal from viviparity to oviparity? The European common lizard (Zootoca vivipara, Lacertidae) and the evolution of parity . Biological Journal of the Linnean Society 87 : 1 – 11 . Google Scholar Crossref Search ADS WorldCat Thorpe RS , Baez M. 1993 . Geographic variation in scalation of the lizard Gallotia stehlini within the island of Gran Canaria . Biological Journal of the Linnean Society 48 : 75 – 87 . Google Scholar Crossref Search ADS WorldCat Trivers RL . 1974 . Parent-offspring conflict . American Zoologist 14 : 249 – 264 . Google Scholar Crossref Search ADS WorldCat Trochet A , Dupoué A, Souchet J, Bertrand R, Deluen M, Murarasu S, Calvez O, Martinez-Silvestre A, Verdaguer-Foz I, Darnet E, Le Chevalier H, Mossoll-Torres M, Guillaume O, Aubret F. 2018 . Variation of preferred body temperatures along an altitudinal gradient: a multi-species study . Journal of Thermal Biology 77 : 38 – 44 . Google Scholar Crossref Search ADS PubMed WorldCat Van Dyke JU , Brandley MC, Thompson MB. 2014 . The evolution of viviparity: molecular and genomic data from squamate reptiles advance understanding of live birth in amniotes . Reproduction 147 : R15 – R26 . Google Scholar Crossref Search ADS PubMed WorldCat Zeh DW , Zeh JA. 2000 . Reproductive mode and speciation: the viviparity-driven conflict hypothesis . BioEssays 22 : 938 – 946 . Google Scholar Crossref Search ADS PubMed WorldCat Author notes " Co-senior authors © 2020 The Linnean Society of London, Biological Journal of the Linnean Society This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model) TI - Mother-offspring conflict for water and its mitigation in the oviparous form of the reproductively bimodal lizard, Zootoca vivipara JF - Biological Journal of the Linnean Society DO - 10.1093/biolinnean/blaa012 DA - 2020-03-28 UR - https://www.deepdyve.com/lp/oxford-university-press/mother-offspring-conflict-for-water-and-its-mitigation-in-the-Y6jRpRgx6a SP - 888 VL - 129 IS - 4 DP - DeepDyve ER -