How do climate-linked sex ratios and dispersal limit range boundaries?
© Boyle et al.; licensee BioMed Central Ltd. 2014
Received: 6 April 2014
Accepted: 11 June 2014
Published: 10 July 2014
Skip to main content
© Boyle et al.; licensee BioMed Central Ltd. 2014
Received: 6 April 2014
Accepted: 11 June 2014
Published: 10 July 2014
Geographic ranges of ectotherms such as reptiles may be determined strongly by abiotic factors owing to causal links between ambient temperature, juvenile survival and individual sex (male or female). Unfortunately, we know little of how these factors interact with dispersal among populations across a species range. We used a simulation model to examine the effects of dispersal, temperature-dependent juvenile survival and sex determining mechanism (temperature-dependent sex determination (TSD) and genotypic sex determination (GSD)) and their interactions, on range limits in populations extending across a continuous range of air temperatures. In particular, we examined the relative importance of these parameters for population persistence to recommend targets for future empirical research.
Dispersal influenced the range limits of species with TSD to a greater extent than in GSD species. Whereas male dispersal led to expanded species ranges across warm (female-producing) climates, female dispersal led to expanded ranges across cool (male-producing) climates. Two-sex dispersal eliminated the influence of biased sex ratios on ranges.
The results highlight the importance of the demographic parameter of sex ratio in determining population persistence and species range limits.
For many organisms, particularly ectothermic animals, individual traits and demographic parameters are linked to temperature, such that ranges could be largely determined by abiotic factors [1, 2]. Abiotic factors are often incorporated into the most simple modelling approaches (for example, climate envelope modelling) used to predict the future distributions of organisms . However, there is increasing evidence [3–8] that fast evolutionary changes, dispersal and population dynamics are among the factors that are equally as important as abiotic factors in determining species ranges [3, 9]. For example, it has been demonstrated that population dynamics and dispersal can alter a species responses to shifting climates . In particular there can be a considerable lag between climate becoming suitable in any one place and colonisation by a species through dispersal . Across a species’ range, dispersal can play a key role by expanding ranges outside of a thermal limit (i.e. creating sinks,), [Harts A, Schwanz L, Kokko H: Demography can favour female-advantageous alleles. P R SOC B. submitted], and preventing local adaptation . These issues have gained new imperative for understanding range change in the context of climatic warming.
Colonisation and extinction rates may change spatially due to environmental gradients in habitat characteristics (for example, temperature, rainfall, or humidity) . The role of temperature in demography and range limits may be especially important for species with temperature-dependent sex determination (TSD; [12, 13]). In these species, individual sex (male or female) is permanently determined during embryonic development by incubation temperature . The proportion of new hatchlings in a population that are male (the Cohort Sex Ratio, or CSR) often shows wide geographical variation associated with climatic or microclimatic variation and their effects on nest temperatures [15, 16]. Numerous studies have predicted that climate warming will skew the cohort sex ratios of reptiles towards females (or males in the tuatara) with negative consequences for population persistence (e.g. [12, 15, 17–22]). Yet we know very little about how sex ratios in populations of reptiles with TSD limit range boundaries. There has been at least one attempt to explain the geographical distribution of TSD species based on sex ratios, suggesting that imbalanced sex ratios (and marginal habitats) limit population growth at range boundaries . It was proposed that balanced sex ratios are located at range centres (or favourable habitats), and skewed sex ratios at more marginal habitats .
The importance of biased sex ratios in limiting the ranges of species with TSD will depend on a few interacting factors. First, climatically-linked juvenile survival may also play a role in determining range boundaries, as eggs only survive at certain temperatures [16, 23] [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted]. Second, the extent to which female-biased sex ratios limit population persistence will depend on how ‘limiting’ males are (i.e. the dependence of female fecundity on male abundance, [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted] . If male density is not strongly limiting to female fecundity because of polygyny, sperm storage and sex-specific breeding intervals [25–27], female-biased populations will have enhanced population growth and persistence [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted] .
Finally, dispersal among populations may bring the rare sex into biased populations, rescuing them from the risk of demographic collapse and potentially leading to expanded ranges. Dispersal by male hatchlings is thought to have an important role in facilitating population persistence in increasingly feminised populations . However, most studies of the effects of climate warming on sex ratios are of isolated populations, and no consideration has been given to the role of dispersal, or the interaction between climatically-linked sex ratios and sex-specific dispersal in driving population dynamics. For example, in many TSD species dispersal may be inefficient for range expansion under climate change . This suggests that the sensitivity of range limitation to dispersal distance should be examined.
Dispersal in reptiles is thought to be primarily male-biased, although data are sparse and alternative dispersal tendencies (female-biased dispersal or dispersal by both sexes) have not been well investigated [30, 31]. For example, male-biased dispersal has been demonstrated in many species of lizards [32, 33], snakes [31, 34, 35], marine turtles [36–39], and in some freshwater turtles [40–42]. Female dispersal has been demonstrated only in an alpine lizard (Niveoscincus microlepidotus) .
There are no documented studies of two-sex dispersal in reptiles. A theoretical analysis of two-sex dispersal tendency found it enhanced population persistence in marginal habitats as gene flow through dispersal occurred in both directions . This analysis included benefits associated with choosing habitats to maximize environmental quality and a reduction in inbreeding. However, because there was no population variation in sex ratio, it did not explore the importance of dispersal in recruiting a rare sex into a population.
In this paper, we explore the role of dispersal in determining the extent of population persistence and range limits in species of reptiles with the common TSD pattern 1A (females produced at higher temperatures) and compare it to results for species with genotypic sex determination (GSD). We do this across populations that vary spatially, but not temporally, in air temperature. Importantly, we include in this analysis 1) the potential for male limitation on female fecundity when males are rare, and 2) the additional effect of temperature on juvenile survival.
A matrix (100 rows by 100 columns) of 10,000 populations was distributed across a continuous air temperature gradient, with each column assigned a temperature from 18°C to 33°C. Each population in the matrix was initiated with 100 adult males and 100 adult females (for both TSD and GSD species). We projected these populations in a simulation to determine range limits under different scenarios of dispersal and temperature-dependence of sex ratios. Within each population, simulated population operations approximated logistic growth, including density-dependent juvenile survival [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted].
In each time step, offspring were first produced. The total number of offspring produced was the product of the number of adult females in the population and per-female fecundity, which was subject to male limitation (see Male Limitation). The number of offspring in each population that were male was sampled using a random binomial distribution (and, hence sex ratio was stochastically determined), given the total number of offspring and the population-specific sex ratio probability of producing a male, p. This sex ratio probability depended on population air temperature for TSD scenarios, but was a constant 0.5 for GSD scenarios (see Cohort sex ratio). The number of offspring that recruited as adults into their natal population or a non-natal population (survived to and bred at age 1) depended on temperature-dependent juvenile survival (see Juvenile survival) followed by dispersal (see Dispersal function). The number of adults surviving to the next time step was sampled using a random binomial distribution given the initial number of adults and a constant probability of survival (s = 0.95 for both males and females; average life expectancy of 20 years). There was no adult dispersal.
The population cycle was iterated 1,000 times. Throughout the simulation, populations were counted as extinct if they had either zero adult males or zero adult females, and were excluded from reproduction in that round. Each simulation was replicated ten times and the averaged numbers of extant populations were plotted in ten temperature intervals, each of 1.5°C in width, in the range of temperatures in the population matrix. Simulations were performed with MATLAB R2012b.
The number of juvenile males and females surviving the embryonic stage was sampled using a random binomial distribution based on the number of offspring produced of each sex and the juvenile survival rate. Juvenile survival rate was density dependent, a e-cN, where N is the total number of adult males and females in the population, c is the density-dependent constant (set to 0.001), and the baseline survival of juveniles (a) depended on temperature according to a normal distribution (the temperature-dependent embryonic survival (TS) curve, [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted].
The TS curve had a maximum baseline embryonic survival of a max =0.015 (at 25°C) and a minimum baseline survival value of zero. The range of temperatures that produced non-zero baseline juvenile survival probabilities was 18 to 33°C. Survival rates are based on estimates from published values . The temperature ranges for juvenile survival were similar to those reported for turtles with GSD  and TSD [15, 46–48].
Several dispersal scenarios were evaluated owing to insufficient empirical information available to accurately parameterise dispersal in reptile species and variation among species in their dispersal tendencies. We explored three levels of dispersal: none, ‘small’ and ‘large’. However, these terms are relative as there is no readily available information to estimate what small or large dispersal would be in these species. For the ‘small’ amount of dispersal, juveniles had a higher probability of not dispersing and a lower probability of travelling large distances compared to juveniles with a ‘large’ amount of dispersal (Figure 1b). ‘Small’ and ‘large’ amounts of dispersal had three separate sex-based tendencies, i.e. male only, female only and two-sex dispersal. Under two-sex dispersal, males and females had equal probabilities of dispersing. Due to the stochastic natures of the juvenile sex ratio in each population and the dispersal function (below), the sex ratio of dispersers among cells was therefore also stochastic.
Pr(dispersal distance) is the probability that an individual moves a certain distance, Dij, the distance moved between cells or populations. Dij includes a distance of zero, thereby including the probability of not dispersing. Parameter A defines the distribution of dispersal differences, 1/A is the average dispersal distance, and Z is a shape parameter for the dispersal curve . Pr(dispersal distance) represents the relative probability that an individual moves a certain distance between cells compared to other possible distances. This was converted to an actual or ‘real’ probability of dispersing each distance by dividing each relative probability by the sum of the values for all dispersal distances in the matrix (range of distances is 0–198 cells, see below).
Distances were calculated between cells (populations) as a von Neumann neighborhood, which counts the 4 cells immediately to each side of the focal cell as a distance of one, and the diagonal 4 cells as a distance of two . Firstly, dispersal distance for each juvenile in each population was chosen based on the probabilities specified by the dispersal kernel. If a distance of zero was chosen the individual did not disperse and recruited to its natal population. If the distance was greater than zero, the target cell for dispersal was chosen randomly from all cells of the specified distance. Thus, dispersal from the cell in which the individual was born occurred in any direction on the matrix with equal probability. The edge cells of the matrix were defined as the ‘boundaries’ beyond which an individual could move no further. If a dispersal distance was chosen to which no cells corresponded (i.e. off the matrix), a new distance was chosen. Note that maximum possible dispersal distance varied among cells. Only the four corner cells of the matrix had a non-zero probability of dispersing a distance of the maximum 198 cells. Cells along the edge had higher maximum possible dispersal distances, but were constrained to disperse in fewer directions compared to cells in the center, which were more limited in distance but could disperse in more directions. We do not think this variation introduces much bias into the model as the parameter values we chose for the dispersal function specify exceptionally low probabilities of dispersing beyond 20–50 cells. An additional, small probability (0.1) of dispersal-related mortality in juveniles was incorporated into the dispersal algorithm. Dispersal related mortality did not increase with distance moved.
When temperature influences the primary sex ratio, the results share some similarities with GSD, but also have several important differences. For both TSD and GSD species, there are large numbers of surviving populations around the central temperature intervals (from 24 to 28.5°C, Figures 3 and 2). In TSD species this 4.5°C temperature range represents a transition from a slightly male-biased sex ratio (around 60% male at 24°C) to a very female-biased sex ratio (around 80% female at 28.5°C). Without dispersal, TSD species have smaller ranges compared to GSD species. The range edge at warm temperatures is reduced modestly owing to decreased female fecundity where males are limiting. The range edge at cool temperatures is strongly reduced owing to low numbers of females, reducing overall offspring production.Under male dispersal, when males move into warmer areas they encounter increasingly larger numbers of females, and populations are no longer male limited. The importance of male limitation on female fecundity in influencing range limits is evident by comparing the results with and without male dispersal (Figure 2b, c, d). In the absence of dispersal (Figure 2b) warm populations are limited by the lack of males. In contrast, small and large amounts of male dispersal result in considerably greater population persistence in warmer areas (Figure 2c, d). Populations in the temperature intervals from 28.5 to 32.5°C (Figure 2c, d) benefit most from dispersal by males, as these very female-biased populations are able to persist in warmer areas where they do not in the absence of dispersal. Male dispersal cannot expand the range at the cold edge, where only males are produced, because dispersing males mostly encounter other males.Female dispersal produces the opposite result to male dispersal, expanding the range in the cool climates. This is because females dispersing to colder areas encounter a higher proportion of males than in populations of reptiles with GSD (Figure 3f). In contrast, if females disperse to warmer areas they encounter only other females, leading to consequent declines in reproduction and population persistence.Two sex dispersal (Figure 2g, h) increases ranges, similar to GSD species (Figure 3g, h). The shapes of the distributions of TSD and GSD species with two sex dispersal are very similar.
There is increasing evidence that factors such as population dynamics and dispersal strongly interact with abiotic factors such as climate in determining species ranges [3–8]. Species range limits can be strongly influenced by climatic effects on demographic parameters such as juvenile or adult survival, growth and reproduction [1, 2]. Dispersal can play a key role in rescuing failing populations or creating “sinks” at range edges. Along range margins if the rate of colonisation via dispersal exceeds local extinctions, the range will expand. By contrast if local extinction events become more frequent than colonisations owing to extreme climates then the species range will shrink . Here, we demonstrate the importance of the sex ratio in limiting species ranges, and explore its interaction with sex-specific dispersal.
Consistent with our previous work and Kallimanis , we found that, in the absence of dispersal, species with TSD are restricted in the climates in which they can persist compared to those with GSD owing to biased sex ratios in both warm and cool climates. In contrast to Kallimanis , we found that the restriction was greater in cool (male-producing) climates than in warm (female-producing) climates. The difference is owing to assumptions regarding male limitation [Boyle M, Hone J, Schwanz L, Georges A: Under what conditions do climate-driven sex ratios enhance versus diminish population persistence? Ecology and Evolution. submitted] – we assume that female fecundity is only reduced when males are less than 10% of the adult population.
In addition, we had several novel findings: 1) when both sexes disperse, the existence of biased sex ratios at extreme climates has almost no effect on limiting ranges; 2) male-only dispersal led to ranges covering areas where sex ratios were slightly male-biased to where they were heavily female-biased; 3) female-only dispersal led to ranges covering areas where sex ratios were slightly female-biased to where they were heavily male-biased; 4) dispersal scenario was more influential in driving range boundaries in TSD species compared to those with GSD. Our findings provide a clear distinction between populations of reptiles with GSD and TSD in the effect of dispersal at range boundaries, and by inference their responses to climate change. Understanding the sex-specific tendencies of dispersal will be imperative for predicting the possibility of range change for species with TSD.
Female-biased populations have been described as more likely to experience growth than populations with even sex ratios . Furthermore, dispersal by male hatchlings is thought to have an important role in facilitating population persistence in increasingly feminised TSD populations under climatic warming . We have demonstrated that male dispersal increases population persistence in female-biased populations located at warmer areas towards range margins (Figure 2c and d).
We show that population persistence at the colder (male-producing) edge of the range occurs only through female dispersal. Yet female dispersal has only been demonstrated in one GSD reptile, the alpine skink . Skinks are considered to be inefficient at dispersal, and dispersal may occur over very short distances of a few metres . In contrast, male marine turtles are considered to be very effective dispersers and may travel across oceans [36, 37, 39, 52].
Two-sex dispersal has been poorly investigated across species . Dispersal of both sexes in our model resulted in a potential for large ranges for both TSD and GSD species compared to no dispersal (Figure 3g, 3h, 2g, 2h). This expanded range led to populations persisting even in locations where juvenile survival was almost zero. In these scenarios, populations located at range edges are likely sinks, where the rate of production is below replacement level, and without sufficient immigration may become extinct . Hence, under a climate change scenario, they may have a poor ability to produce their own migrants capable of expanding the range outwards.
Given the strong impact of the sex-bias and distance of dispersal on the results, these will be key behavioural traits to determine empirically when testing our theoretical predictions or making predictions about climate change. In addition, the distance of dispersal in a species must be considered with respect the spatial scale of its populations (e.g. distance is measured in number of populations that an individual can pass when dispersing). Our model makes the unique and testable prediction that, when comparing species of reptiles with pattern 1A of TSD, those species that have ranges where the sex ratios tend to be female biased will also to have male-mediated dispersal, while those that persist where the sex ratios tend to be male biased will have female-mediated dispersal.
While TSD species may persist in the short-term with biased sex ratios at range margins, an important caveat for our results is that populations at range margins (in extreme climates) may evolve in response to frequency-dependent selection on sex to produce more even sex ratios . That is, local adaptation in the pivotal temperature (temperature at which a 50:50 sex ratio is produced) or maternal nesting behaviours could occur and reduce geographic variation in cohort sex ratios [54, 55]. Existing geographical variation in nesting behaviour (for example, timing, nest depth and shade) may ameliorate the effects of local climate on nest temperatures, thus reduce spatial variation in survival and sex ratio . However, recent work on the co-evolution of the pivotal temperature and dispersal suggests that local adaptation is often limited . Owing to the demographic dominance of females, robust populations that produce a lot of dispersing offspring (i.e. populations with even or female-biased sex ratios) send their locally-adapted genes into other environments where the genes are maladapted and lead to biased sex ratios .
A critical assumption of our model is that each sex has an increasing benefit to population growth as the sex ratio becomes biased towards the opposite sex. However, as females become rarer this assumption may not hold anymore. For example, a low density of females has been demonstrated to result in male aggression and population collapse .
For our theoretical model to have explanatory power for the distributions of species with TSD, sex ratios must vary spatially. There is some evidence that sex ratios may be constant across the ranges of the painted turtle and the water dragon [54, 55]. In contrast, sex ratios vary geographically in some sea turtles and crocodiles, although the variation is idiosyncratic and not linked to a continuous spatial component such as latitude [15, 16]. A recent paper proposes that TSD reptiles do not exhibit balanced sex-ratios at the centre of their geographic distributions, and biased sex-ratios at range margins . Future theoretical work should incorporate evolution of TSD traits as well as spatial variation in sex ratios that is more idiosyncratic across the range.
More broadly, testing our model and understanding the influence of climate of the ranges of species with TSD requires additional empirical data on the key demographic parameters across the ranges of TSD and GSD species. While data on dispersal are hard to collect, inferences about dispersal can be made from genetic analyses. Genetic analyses of sea turtles suggest that gene-flow is male-mediated and hence male-biased dispersal is the likely pattern [37–39, 52]. Crude estimates of dispersal could be made by combining fixation index (FST) (a measure of population differentiation owing to genetic structure) and demographic parameters. There is also some behavioural data which suggest that male-biased dispersal is the likely pattern in freshwater turtles [40–42], and this could be further investigated using techniques such as radio telemetry.
Our theoretical model revealed that dispersal was more influential in determining the range limits of species with TSD than GSD. Male dispersal led to expanded species ranges across warm climates, where more females were produced, and female dispersal led to expanded ranges across cool climates, where more males were produced. Two-sex dispersal eliminated the influence of biased sex ratios on ranges.
This study was supported by an Australian Postgraduate Research Award to Maria Boyle, and supplementary funding by the Institute for Applied Ecology at the University of Canberra. Lisa Schwanz was funded by a US National Science Foundation International Research Fellowship and a University of Canberra Postdoctoral Fellowship. We thank the attendees at the Institute journal club for their vigorous discussion of ideas leading up to this paper, and J Lynch for providing comment on an earlier draft. Denis O’Meally established and maintained the high performance computer on which the modelling was done.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.