bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
Behavioural adaptations in egg laying ancestors facilitate evolutionary
transitions to live birth
Amanda K. Pettersen1*, Nathalie Feiner1, Daniel W.A. Noble2, Geoffrey M. While3, Charlie
K. Cornwallis1# & Tobias Uller1#
1. Department of Biology, Lund University, Lund, 22 362, Sweden
2. Division of Ecology and Evolution, Research School of Biology, The Australian
National University, Canberra, 2600, Australia
3. School of Natural Sciences, University of Tasmania, Hobart, 7005, Australia
#Joint senior author
*Corresponding author: Amanda K. Pettersen. Email: [email protected]
Keywords: Viviparity, Oviparity, Thermal plasticity, Behavioural plasticity, Reproductive
mode
1 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
1 Abstract
2 Live birth is a key innovation that has evolved from egg laying over 100 times in reptiles.
3 One significant feature in this transition is the thermal conditions experienced by developing
4 embryos. Adult lizards and snakes often have preferred body temperatures that can be lethal
5 to developing embryos and should prevent egg retention: how has viviparity repeatedly
6 evolved in the face of this pervasive mismatch? Here we resolve this paradox by conducting
7 phylogenetic analyses using data on thermal preference from 224 species. Thermal
8 mismatches between mothers and offspring are widespread but resolved by gravid females
9 behaviourally down-regulating their body temperature towards the thermal optimum of
10 embryos. Importantly, this thermoregulatory behaviour evolved in ancestral egg-laying
11 species before the evolutionary emergence of live birth. Maternal thermoregulatory behaviour
12 therefore bypasses constraints imposed by a slowly evolving thermal physiology and is likely
13 to have been a key requirement for repeated transitions to live birth.
2 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
14 Introduction
15 The evolution of live bearing is an important life-history transition in vertebrates (1–3). The
16 ecological conditions that favour the transition from egg laying (oviparity) to live birth
17 (viviparity) are relatively well understood, especially in reptiles, with particularly strong
18 support for the adaptive value of viviparity in cool climates (3, 4). By retaining embryos
19 throughout development, mothers can buffer offspring from suboptimal nest temperatures,
20 thereby ensuring higher hatching success and increased offspring viability (4). Live bearing
21 also allows females to maintain offspring at higher temperatures than those of potential nest
22 sites, effectively increasing developmental rates to facilitate early-season hatching (5, 6).
23 This has allowed reptile species to diversify and persist in cool climates globally (7).
24 Despite its adaptive advantage, it is unclear how viviparity initially evolved. A major
25 challenge to the evolution of viviparity is that adult lizards and snakes typically have
26 preferred body temperatures that exceed the upper lethal limit of embryos (8). For example,
27 the average nest temperature for the Iberian emerald lizard, Lacerta schreiberi, is 24°C,
28 rarely exceeding 30°C, while the preferred body temperature of females is 33°C (9). Since
29 embryos are well adapted to the temperatures they typically experience in the nest, they
30 generally have limited capacity to sustain development at temperatures outside a narrow
31 thermal window (10, 11). Prolonged exposure to temperatures beyond this range can result in
32 offspring abnormalities and even death (12). A mismatch between thermal optima of embryos
33 and females should therefore prevent the maternal retention of embryos throughout
34 development, and therefore evolutionary transitions to live birth (8).
35 Despite this apparent constraint, live birth has evolved over 100 times in squamate
36 reptiles (13, 14). How can the repeated evolution of viviparity be reconciled with the
37 widespread thermal mismatch between embryos and adults? There are two potential scenarios
38 by which viviparity can evolve in the presence of maternal-embryo thermal mismatches.
3 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
39 First, viviparity may only evolve in species where there is no mismatch, that is, adult body
40 temperature (Pbt) and embryo thermal optima (Topt) are well aligned. Second, females may
41 behaviourally adjust their body temperature when pregnant to close the gap between adult
42 and embryo thermal optima, even when substantial mismatches exist. Such plasticity may
43 come at a cost to female performance, but temporarily shifting body temperature while gravid
44 to match the thermal optima of embryos could eliminate thermal barriers to the evolution of
45 viviparity.
46 Elucidating how the mismatch between female and offspring thermal optima have
47 been resolved in viviparous lineages is challenging. Retracing the evolutionary steps that
48 have led to viviparity requires the reconstruction of thermal preferences of adults and
49 embryos at the very origins of these lineages. Here, we collate data on the thermal
50 preferences of female adults (163 species) and embryos (50 species) of squamate reptiles, and
51 use phylogenetic analyses to reconstruct the evolutionary history of live birth (Table S1).
52 First, we estimated the ancestral states of preferred body temperature of females and embryos
53 to quantify the thermal conflict between mothers and embryos at the origins of viviparity.
54 Second, we tested the prediction that transitions to viviparity have been restricted to
55 oviparous lineages with a low degree of maternal-embryo mismatch. This was assessed by
56 calculating if viviparity has originated more frequently from oviparous ancestors where
57 female Pbt overlaps with the thermal optima of embryos. Finally, we examined if females
58 adjust their body temperature when gravid to reduce the mismatch between adult and embryo
59 thermal optima, and if so, whether this behaviour evolved prior to the emergence of
60 viviparity, or afterwards in response to such mismatch.
4 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
61 Results and Discussion
62 Wide-spread mismatch in thermal optima of oviparous mothers and their embryos
63 Across oviparous lizards and snakes, there were many species with a significant disparity in
64 preferred body temperature (Pbt) and the temperature that maximizes hatching success (Topt)
65 (Multi-response BPMM (MR-BPMM): posterior mode (PM) of difference = 5.60, 95%
66 Credible Interval (CI): -0.60, 10.18, pMCMC = 0.06. Figure 1 & 2; Table S3). The female Pbt
67 of oviparous species is on average 4 °C higher than Topt (raw data mean ± SD Pbt 102 species:
68 31.40 ± 4.31 ºC; Topt 50 species: 27.25 ± 2.42 ºC). Incubation experiments have shown that
69 such a difference can significantly reduce hatching success, illustrating that the retention of
70 offspring throughout development may jeopardise embryo survival (15, 16); reviewed in
71 (17). Pbt and Topt of oviparous species were both estimated to have high phylogenetic inertia
2 2 72 (Figure 1; MR-BPMM: Pbt phylo H : 0.99, CI: 0.91, 1.0. Topt phylo H : 0.93, CI: 0.72, 0.98.
73 Table S3) and female Pbt and offspring Topt only show weak evidence of coevolving (Figure
74 2; phylogenetic correlation (MR-BPMM): PM = 0.57, CI: -0.09, 0.97, pMCMC = 0.08. Table
75 S4), suggesting that this mismatch may be difficult to resolve.
76
77 Transitions to live birth occur despite mismatches between females and embryos
78 We found evidence that preferred body temperature was lower in viviparous, compared with
79 oviparous, lineages (MR-BPMM: PM = -0.20, CI: -0.43, -0.02, pMCMC = 0.02; Table S6).
80 This is consistent with the hypothesis that viviparity is an adaptation to cool climates (2, 3,
81 18). However, reconstructions of ancestral preferred body temperatures of adults indicated
82 that in the majority of lineages, lower preferred body temperatures evolved after, not before,
83 transitions to live bearing. This was evident from the lack of difference in Pbt between
84 ancestors of egg-laying and live-bearing lineages (MR-BPMM, PM = 0.04, CI: -3.96, 4.27,
85 pMCMC = 0.42; Table S7). Additionally, live-bearing lineages with egg-laying ancestors
5 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
86 showed higher Pbt than those with live-bearing ancestors (MR-BPMM, PM = 6.78, CI: -0.64,
87 8.74, pMCMC = 0.06; Table S7).
88 Ancestral reconstructions of Pbt and Topt of oviparous species showed that transitions
89 from egg laying to live bearing were no more likely to occur when adult and embryo thermal
90 optima were similar (% of lineages with thermal mismatches: ancestors of live bearing 84%,
91 ancestors of egg laying species 85%. Table S5). Only 5 of the estimated 31 evolutionary
92 origins of live birth had adult and embryo thermal optima that overlapped, which was similar
93 to the proportion of lineages where live bearing did not occur (Table S5). In the remaining 26
94 egg-laying ancestors of live-bearing species, the thermal optimum of embryos was
95 significantly lower than the preferred body temperature of adults (there were no cases of Topt
96 being significantly higher than Pbt; Table S5). In summary, the degree of similarity between
97 adult Pbt and embryo Topt does not predict transitions towards viviparity, and the lower Pbt
98 observed for viviparous taxa appears to have evolved after transitions to live birth, not before.
99
100 Female behavioural plasticity resolves the constraints imposed by thermal physiology
101 Next, we tested the hypothesis that the thermal mismatch between females and embryos
102 during the evolution of viviparity is resolved by females adjusting their body temperature
103 when gravid to meet the thermal requirements of their embryos (19). We collected data on
104 the preferred body temperature of gravid (Pbt-g) and non-gravid (Pbt-ng) adult females from 52
105 species of lizards and snakes were collected from published field and experimental studies
106 (live bearing n = 32, egg laying n = 20. Table S1). Using these data, an effect size of female
107 behavioural plasticity, Hedges’ g was calculated (positive values = higher Pbt when gravid,
108 negative values = lower Pbt when gravid).
109 Across those 52 species, there was a negative phylogenetic correlation between Pbt
110 and the degree to which females reduced their body temperature when gravid (MR-BPMM:
6 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
111 PM = -0.29, CI: -0.62, -0.02, pMCMC = 0.02. Figure 4; Table S8). Females with high Pbt
112 have therefore evolved to behaviourally reduce their body temperature when gravid to a
113 greater extent than females with low Pbt. This relationship was not restricted to viviparous
114 lineages, which is expected if female behavioural plasticity is an adaptation that evolves in
115 response to selection for embryo retention. Instead, comparable phylogenetic correlations
116 between Hedges’ g and Pbt were evident in both oviparous and viviparous taxa, showing that
117 behavioural adjustment of temperature when gravid occurs to a similar extent in both egg
118 laying and live bearing species (Egg laying: PM = -0.61, CI: -0.92, -0.08, pMCMC = 0.02.
119 Live bearing: PM = -0.88, CI: -0.96, -0.30, pMCMC = 0.005. Table S9).
120 To further investigate the role of female plasticity in the evolution of live bearing, we
121 used phylogenetic models to reconstruct values of Hedges’ g for the ancestors of egg-laying
122 and live-bearing species. Consistent with the results from our correlational analyses, we
123 found that, in the ancestors of both egg laying and live bearing species, thermal mismatches
124 between adults and embryos were resolved by females adjusting their body temperature when
125 gravid (Figure 4A; MR-BPMM: Egg-laying Hedges’ g PM = -0.81, CI = -1.43, 0.05,
126 pMCMC = 0.03. Live bearing Hedges’ g PM = -0.75, CI = -1.88, -0.09, pMCMC = 0.03.
127 Table S10). As a result, Hedges’ g did not differ between the ancestors of oviparous and
128 viviparous species (Figure 4B; Table S9-S11), suggesting that behavioural plasticity was
129 likely present prior to the emergence of live birth.
130 Our results illustrate that an ancient and evolutionarily conserved ability of egg laying
131 females to adjust their body temperature during pregnancy provides a general resolution to
132 the conflict over adult and embryo thermal optima that paves the way for the evolution of live
133 bearing. The down-regulation of female body temperature in egg-laying species while gravid
134 may appear surprising considering that most of these species lay their eggs within the early
135 stages of development (commonly around the time of limb bud formation (20)). However,
7 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
136 early developmental stages, involving gastrulation, neurulation and organogenesis, are
137 potentially even more sensitive to temperature stress than later stages, which are
138 predominantly associated with growth (8, 10). The temperature sensitivity of early-stage
139 embryos therefore suggests that resolutions to mother-offspring thermal conflicts are required
140 in both egg-laying and live-bearing species. If true, the key innovation of live birth may owe
141 its evolutionary origin to mechanisms of temperature regulation put in place long before live
142 birth emerged. We suggest that this behavioural flexibility of female lizards and snakes
143 enabled the gradual evolution of live birth by helping traverse potential fitness valleys
144 associated with embryo retention. In many live-bearing lineages, temperature mismatches
145 have subsequently been eliminated through the evolution of a lower preferred body
146 temperature in females. However, behavioural flexibility is frequently maintained in
147 viviparous lineages that have colonized cool climates, possibly as a means for mothers to
148 maintain embryos at temperatures significantly warmer than the external environment (2, 5,
149 6, 21, 22).
150 In conclusion, female thermoregulatory behaviour is a key adaptation that overcomes
151 the constraints imposed by thermal mismatches between mothers and embryos, and has
152 facilitated both the evolution of live birth, and expansion of reptiles into a variety of
153 environments.
8 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
154 Methods
155 Literature search and data collection
156 To investigate the relationship between maternal behavioural plasticity and embryo thermal
157 sensitivity, we collected existing data on reproductive mode (oviparous versus viviparous),
158 the female preferred body temperature (Pbt) (“Pbt dataset”), the temperature at which hatching
159 success was maximised (Topt) (“Topt dataset”), and female body temperature when gravid (Pbt-
160 g) and not gravid (Pbt-ng) (“Hedges’ g dataset”). Complete data are presented in Table S1. Data
161 were collected for each of the variables from published literature. We conducted a literature
162 search using ISI Web of Science (v.5.30) with search terms specific to each dataset.
163 The Pbt dataset provided preferred body temperatures (Pbt) of adult females measured at
164 one point in time. We used the search terms ‘body temperature*’, along with one of the
165 following: ‘squamat*’, ‘lizard*’, ‘snake*’ which yielded a total of 1075 papers. We only
166 used data from studies where Pbt for females was stated explicitly, unless pooled male/female
167 data stated no significant effect of sex. Data on females that were described to be gravid or
168 data collected during the reproductive season were excluded. This provided a final dataset of
169 female preferred body temperature for 163 species that was independent of the gravid and
170 non-gravid measures used to calculate Hedges’ g (live bearing: n = 61 and egg laying: n =
171 102).
172 The Topt dataset quantified the temperature that maximises embryonic survival (measured
173 as hatching success; Topt) in 50 egg laying species. Hatching success data were taken from the
174 Reptile Developmental Database (23) and supplemented with additional studies using the
175 same method described for the Pbt dataset (search terms: ‘temperature* AND incubat* AND
176 hatch* OR surv*’, along with one of the following: ‘squamat*’, ‘lizard*’, ‘snake*’), yielding
177 a total of 671 papers. We only included studies where three or more constant temperature
178 treatments were used under controlled laboratory conditions, resulting in 661 papers being
9 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.07.430163; this version posted February 8, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
179 rejected due to irrelevance or overlap with the Reptile Development Database. Given that Topt
180 had a strong phylogenetic signature (Phylogenetic Heritability (H2) = 91.93%, 95% CI: 71.65
181 – 98.16%; Table S3) we fitted a Bayesian Phylogenetic Mixed Effects Model (BPMM) to
182 hatching success data for all species jointly. Including phylogenetic information allowed for
183 the Topt of each species to be estimated with greater accuracy and precision given that the
184 range and number of temperatures across species varied (range: 10-40 ºC, mean number of
185 temperatures per species ± SD: 7.62 ± 4.29). Compared to non-phylogenetic models, Topt
186 estimates produced from phylogenetic models (BPMM) showed smaller sampling error and
187 avoided convergence problems in estimating model parameters. Importantly, this approach
188 was not used to estimate Topt data for species without data, only to better predict Topt values
189 for species for which there were data. The Topt BPMM model was run for 1,100,000 iterations
190 with a burn-in of 10,000 iterations and thinning rate of 500, leaving us with 2,000 samples
191 from the posterior distribution. Autocorrelation was low (lag values < 0.1) and trace plots
192 showed chains mixed well for all parameters. Our model included temperature as a fixed
193 effect (estimating both a linear and quadratic slopes) and random slopes of temperature
194 (linear and quadratic slopes) fitted at the phylogenetic level. From our BPMM we estimated
195 Topt, and its corresponding sampling variance, using the posterior distribution of fixed effects
196 and best linear unbiased predictors (BLUPs) for the random slopes (linear and quadratic) for
197 each species as follows: