- Open Access
Environment, but not genetic divergence, influences geographic variation in colour morph frequencies in a lizard
BMC Evolutionary Biologyvolume 15, Article number: 156 (2015)
Identifying the causes of intraspecific phenotypic variation is essential for understanding evolutionary processes that maintain diversity and promote speciation. In polymorphic species, the relative frequencies of discrete morphs often vary geographically; yet the drivers of spatial variation in morph frequencies are seldom known. Here, we test the relative importance of gene flow and natural selection to identify the causes of geographic variation in colour morph frequencies in the Australian tawny dragon lizard, Ctenophorus decresii.
Populations of C. decresii are polymorphic for male throat coloration and all populations surveyed shared the same four morphs but differed in the relative frequencies of morphs. Despite genetic structure among populations, there was no relationship between genetic similarity or geographic proximity and similarity in morph frequencies. However, we detected remarkably strong associations between morph frequencies and two environmental variables (mean annual aridity index and vegetation cover), which together explained approximately 45 % of the total variance in morph frequencies.
Spatial variation in selection appears to play an important role in shaping morph frequency patterns in C. decresii. Selection associated with differences in local environmental conditions, combined with relatively low levels of gene flow, is expected to favour population divergence in morph composition, but may be counteracted by negative frequency-dependent selection favouring rare morphs.
Understanding the processes generating intraspecific phenotypic variation, and how they relate to population divergence and speciation, is a fundamental goal for evolutionary biologists. Polymorphic species, in which multiple, discrete phenotypic variants coexist within a population , exhibit extreme, easily measurable, intraspecific phenotypic variation that often varies geographically in the composition and frequency of morphs. Such geographic variation in polymorphism may act as a precursor to speciation if the processes generating geographic variation in morph frequencies also promote phenotypic and genetic differentiation leading to reproductive isolation among populations [2–4]. Indeed, colour polymorphism has been shown to accelerate speciation in birds .
One of the most important processes generating morph frequency variation is spatial variation in selection and local adaptation . This has been shown to influence population morph frequencies in multiple species (e.g. [2, 6]), including associations between morph frequencies and broad scale environmental gradients (e.g. [7, 8]). In many cases, however, spatial variation in polymorphism cannot be explained by selection alone [9, 10], and it is therefore likely that genetic drift and gene flow also play a role. Despite a rich history of research on polymorphic species, the relative importance of selection, genetic drift and gene flow for explaining geographic variation in polymorphism remains poorly understood (reviewed in ).
Gene flow can be an important determinant of geographic variation in morph frequencies (e.g. ) and is expected to constrain local adaptation by introducing individuals from populations experiencing different selective pressures [12, 13]. Therefore, high levels of gene flow may favour the maintenance of polymorphism and homogenise morph frequencies across populations . Conversely, when gene flow is low, geographic variation in polymorphism may reflect local selective pressures . Notably, when local selection is particularly strong, population divergence in morph frequencies can be maintained despite ongoing gene flow [16, 17]. However, local adaptation driven by strong directional selection may be opposed by negative frequency-dependent selection (NFDS), where rare morphs are favoured, generating cyclical changes in morph frequencies over time [18–20]. For example, female preference for novel mates and predator avoidance of rare, unfamiliar colour morphs, likely maintain extreme colour polymorphism in male guppies . Therefore, frequency-dependent selection, spatial variation in selection driven by environmental factors, genetic drift and gene flow can interact in complex ways to influence morph frequencies (e.g. ). Consequently, in order to determine the processes generating geographic variation in polymorphism, and their relative contribution, it is necessary to assess the relationship between morph frequencies and multiple factors, including local environmental variables, broad environmental gradients and gene flow.
Although many polymorphic species exhibit population differences in morph composition and frequencies, few studies have examined geographic variation in polymorphism with specific focus on its causes and evolutionary consequences (reviewed in ). Here, we evaluate processes that may influence morph frequencies in the Australian tawny dragon lizard, Ctenophorus decresii, in which males exhibit striking throat colour variation, both within and among populations . The species comprises two genetically and phenotypically distinct lineages, characterised by unique throat coloration, and corresponding to a polymorphic ‘northern’ lineage and a monomorphic ‘southern’ lineage [22, 23]. There are four discrete throat colour morphs within the ‘northern’ lineage: orange, yellow, orange surrounded by yellow (orange-yellow) and grey (which lacks both orange and yellow; ). Data from field mark-recapture studies (M. Yewers and D. Stuart-Fox, unpublished data), and long-term captive breeding experiments [25–27], show that throat coloration is highly heritable, fixed at sexual maturity and does not vary with body size, such that all morphs are observed in the same age/size class . Previous work on C. decresii has also detected considerable morph frequency variation among populations ; however, this geographic variation is yet to be explicitly documented.
We quantify and compare colour morph frequencies in eight populations of C. decresii and test the potential roles of gene flow and selection in shaping morph frequency variation. This system is well suited to studying how the interaction between gene flow and selection influences morph frequencies. Firstly, C. decresii is a habitat specialist, restricted to rocky outcrops throughout South Australia, which may form discrete ‘island’ populations with varying degrees of connectivity between populations. Secondly, throat coloration has an important signalling function , as males display their throats during both territorial and courtship displays , and throat coloration is the main feature distinguishing C. decresii from its closest relatives . Additionally, throat coloration in the ‘northern’ and ‘southern’ lineage is locally adapted to different signalling environments. Lineage specific throat colours are more conspicuous to conspecific lizards when viewed against native rock and lichen backgrounds . Finally, colour polymorphism is associated with alternative behavioural strategies. The orange morph exhibits high levels of aggression towards territory intruders in the field, while the grey morph is the least aggressive and the least bold when exposed to a simulated predator attack. The remaining morphs (yellow and orange-yellow males) appear to be intermediates and show aggression that is conditional on the colour morph of the intruder, reacting more aggressively to orange or yellow males than to grey (M. Yewers and D. Stuart-Fox, unpublished data). Therefore, environmentally-driven selection, acting directly on colour or on correlated traits, is expected to influence throat coloration in C. decresii. By quantifying the relationship between morph frequencies and genetic and environmental factors, our study provides insights into processes generating and maintaining geographic variation in polymorphism.
Sites and data collection
Field work was conducted in South Australia during spring and summer between 2010 and 2012. We sampled eight populations within the polymorphic ‘northern’ lineage of Ctenophorus decresii including; Aroona (31°16′39.82" S, 138°36′19.61" E), Wilpena (31°32′24.59" S, 138°36′4.46" E), Yourambulla Caves (31°57′12.99" S, 138°22′24.79" E), Warren Gorge (32°11′15.34" S, 138°00′32.90" E), Devil’s Peak (32°24′57.19" S, 137°59′27.69" E), Mt Remarkable (32°50′39.89" S, 138°03′7.52" E) and Telowie Gorge (33°01′21.43" S, 138°06′24.35" E) in the Flinders Ranges, and Bimbowrie Station (31°58′40.61" S, 140°12′37.65" E) in the Olary Ranges. Within each population, we captured approximately 50 lizards (at least 20 adult males; see Fig. 1 for sample sizes) and took a tissue sample (10 mm tail clip) for subsequent molecular analysis.
Males were photographed using a Canon PowerShot SX1-IS digital camera (saved in RAW format; ) and each photograph included a Micnova grey card. Photographs were then linearised for radiance and calibrated (relative to the grey card) to account for differences in illumination between images . A previous study by Teasdale et al.  used discriminant function analysis of colour and pattern variables and spectral reflectance data to demonstrate that C. decresii males can be reliably classified into discrete throat colour morphs based on the presence or absence of orange and yellow. Using the standardised photographs and these criteria, we visually allocated males into one of four morph types: orange, orange-yellow, yellow or grey. To ensure that morph frequencies did not fluctuate temporally during our study period, two populations (Telowie Gorge and Warren Gorge) were sampled across consecutive years. Morph frequencies were compared among populations and across years using χ 2 tests in R (chisq.test, chisq.post.hoc; ).
Genetic structure and connectivity
Genomic DNA was extracted from tissue with proteinase-K and a GenCatch (TM) Blood and Tissue Genomic Mini-Prep Kit (Epoch Life Sciences, Sugar Land, TX, USA). We amplified eight microsatellite loci previously developed for C. decresii (Ctde03, Ctde05, Ctde08, Ctde12, Ctde21 and Ctde45; ), or for the congeneric Ctenophorus pictus (CP10, CP11; [32, 33]), using published polymerase chain reaction (PCR) protocols . PCR products were sent to Macrogen (Korea) for genotyping and fragment sizes were called using Peak Scanner ver. 1.0 (Applied Biosystems, Foster City, CA, USA). For each locus within each population, we calculated observed and expected heterozygosities and Hardy-Weinberg equilibrium in Arlequin ver. 3.5 (Additional file 1: Table S1 ). All loci were checked for the presence of null alleles using Microchecker ver. 2.2.3  and tested for linkage disequilibrium using GenePop ver. 4.2 (Additional file 1: Table S1 ). In total we genotyped 30 individuals from each of the eight populations sampled.
Pairwise FST between populations were calculated in MSA ver. 4.05 . We derived P values based on 104 permutations and applied Bonferonni correction for multiple tests. Isolation by distance was assessed using a simple Mantel test performed in the zt software , with pairwise FST as the independent variable and Euclidian distance between sites as the dependent variable. To investigate genetic structure among populations, we performed a Bayesian clustering analysis in STRUCTURE ver. 2.3.3 . The program uses Markov Chain Monte Carlo (MCMC) sampling to probabilistically assign individuals to populations and determine the most likely number of populations (K). We performed 20 independent runs, with random starting seeds, for values of K ranging from 1 to 10 with a burn-in length of 1 × 105 iterations and MCMC length of 1 × 106 iterations. We employed the admixture model with correlated allele frequencies and sampling location as a prior. The most likely value of K was determined by calculating the rate of change in likelihood between successive K values (ΔK; ) in Structure Harvester . Given that several hierarchical levels of population structure may exist, and the ΔK method detects only the uppermost level , we also performed separate STRUCTURE analyses on each of the multi-population clusters inferred in the previous step. This was repeated until the inferred number of clusters (K) was one (e.g. ). We performed the STRUCTURE analyses with the full microsatellite dataset, and excluding locus CP10 which showed departures from Hardy Weinberg equilibrium across most populations (Additional file 1: Table S1); however, results were qualitatively the same.
To investigate potential environmental correlates of geographic variation in colour morph frequencies, we acquired 0.01 degree (approximately 1 km) resolution climatic and topographic layers from the Atlas of Living Australia (ALA; http://www.ala.org.au; ). Given the modest number of populations in our dataset, it was important to minimise the number of variables considered. Therefore, we selected two layers which are likely to be important to C. decresii given that the species is a semi-arid, rocky habitat specialist: mean annual aridity index (the monthly ratio of precipitation to potential evaporation and an indicator of dryness) and topographic relief (the range in elevation). Topographic relief provides a measure of the ruggedness of a site, and given that C. decresii inhabits rocky outcrops and gorges, this is likely to be a better predictor of habitat quality than elevation alone. Within each 1 km cell, it is calculated as the average change in elevation between adjacent sub-cells of a grid with spacing of approximately 250 m (9 s in latitude and longitude). Using the ALA website, we extracted values for these two variables for each lizard, and calculated the average value within our eight study sites (Additional file 2: Figure S1). Given the resolution of the climatic and topographic layers used, the values obtained were the same for most lizards within a site (lizards within 1 km2 of each other); however, the area covered within sites ranged between 0.2–2 km2. Furthermore, for each individual captured, the proportion cover of vegetation < 1 m in height and rocks > 50 cm in diameter were visually estimated within a circular quadrat of 3 m radius centred on the point where the lizard was first sighted. These microhabitat variables were assumed to influence territory quality as vegetation structure can affect thermoregulation, predation risk and prey availability (e.g. [43, 44]), and large rocks are used as lookouts and refuges by territorial lizards, including the congeneric ornate dragon, Ctenophorus ornatus . With the exception of geographic distance and aridity index (Pearson correlation coefficient = 0.636; PROC CORR, SAS ver. 9.3), there was minimal correlation between the environmental variables examined (Additional file 3: Table S2).
Determinants of morph frequency variation
We assessed whether geographic variation in morph frequencies in C. decresii can be explained by genetic and/or environmental factors using multiple matrix regression as implemented in FSTAT ver. 126.96.36.199 . Although using measures of distance can reduce statistical power, and must be interpreted in relation to the raw variables , matrix regressions on distance measures are a powerful method to infer the relative contribution of individual covariates to population differentiation . We included morph frequencies as the response variable and all possible combinations of pairwise FST, geographic distance (to account for spatial autocorrelation), aridity index (aridity), topographic relief (topography), rock cover (rock) and vegetation cover (vegetation) as predictor variables in the models; however, we did not consider models that contained both geographic distance and aridity due to the significant relationship detected between these variables . P values for slope parameters were calculated using 104 permutations.
With the exception of pairwise FST, distance matrices for each of the independent variables were calculated as the Euclidean distance between population pairs . For the dependent variable (morph frequencies) Euclidean distance (between populations A and B) was calculated as:
dist(A, B) = √((A O – B O)2 + (A OY – B OY)2 + (A Y – B Y)2 + (A G – B G)2),
which incorporates population differences in the frequencies of each of the four colour morphs (O: orange, OY: orange-yellow, Y: yellow, and G: grey). The best fitting models for the data were determined using the corrected Akaike’s information criterion (AICc; ). We considered models with AICC differences (ΔAICC) within five units of the best possible model to be strongly supported. Based on this subset of competing models, we measured the relative importance of each variable by summing the normalised AICC weights across all models in which a given variable was present (Table 1; [49, 51]). To further aid interpretation of the relationship between morph frequencies and the predictor variables identified in the best model, we performed separate multiple linear regression analyses in R  using the raw variables (rather than distances). The frequency of each of the morphs was the dependent variable and the independent variables were the best predictors identified from matrix regressions. Morph frequency data were arcsine transformed and environmental data were log transformed prior to analysis to meet model assumptions.
Morph frequency variation
Within Telowie Gorge and Warren Gorge there was no significant difference in morph frequencies across years (Telowie Gorge: χ 2 = 6.98, df = 3, P = 0.073; Warren Gorge: χ 2 = 1.62, df = 3, P = 0.655), which is not surprising given that longevity is approximately 5 years in the wild (Yewers and Stuart-Fox, unpublished data). Furthermore, the relative rank of morphs was consistent; for example in both years yellow was the most common morph at Telowie Gorge, followed by grey, orange-yellow and then orange. Conversely, at Warren Gorge orange-yellow was the most common morph, followed by yellow, orange and then grey (Fig. 1; Additional file 4: Table S3). Each of the eight populations sampled contained all four colour morphs; however, the frequencies of morphs varied significantly among populations (χ 2 = 243.16, df = 21, P < 0.0001; Fig. 1), ranging from 10–53.3 % for orange, 9.5–57.7 % for orange-yellow, 13.3–61.9 % for yellow, and 3.8–36.4 % for grey (Fig. 1). Pairwise comparisons revealed similar morph frequencies in Mt Remarkable and Telowie Gorge (P = 0.111) and in Wilpena, Yourambulla Caves, Warren Gorge and Devil’s Peak (P values > 0.05), with the exception of Warren Gorge and Devil’s Peak within this group (P = 0.027; Additional file 5: Table S4). Morph frequencies at Mt Remarkable/Telowie Gorge, Aroona and Bimbowrie Station were significantly different from all other populations (all P values < 0.0001; Additional file 5: Table S4).
Genetic structure and connectivity
There was genetic subdivision between populations of C. decresii. Although pairwise FST values were low to moderate (ranging from 0.013–0.068), they were statistically significant for all population pairs (with the exception of Mt Remarkable and Telowie Gorge; Additional file 6: Table S5). Genetic similarity did not necessarily reflect geographic proximity as population structure did not conform to a model of isolation by distance (r partial = 0.274, P = 0.171; Additional file 2: Figure S1). The highest ΔK value from the initial STRUCTURE analysis was for K = 2 (ΔK = 22.01, mean likelihood of K = −9901), with Yourambulla Caves forming a separate cluster from all other populations (Fig. 2a). The subsequent analysis excluding Yourambulla Caves inferred three clusters (K = 3: ΔK = 4.33, mean likelihood of K = −8484; Fig. 2b): Aroona/Warren Gorge/Devil’s Peak, Mt Remarkable/Telowie Gorge and Wilpena/Bimbowrie Station. With the exception of Mt Remarkable/Telowie Gorge, these clusters were further separated into individual populations in the final STRUCTURE analyses (Fig. 2c, d).
Factors affecting morph frequencies
Multiple matrix regression analyses revealed that genetic distance (pairwise FST) was not a significant predictor of divergence in colour morph frequencies in C. decresii (r partial = 0.026, P = 0.470; Additional file 7: Figure S2). Of the other variables included in our model, mean annual aridity index (aridity) and proportion cover of vegetation <1 m in height (vegetation) were the most important predictor variables (RIV = 0.90 and 0.89 respectively) for the observed geographic variation in morph frequencies (aridity: r partial = 0.568, P = 0.004, vegetation: r partial = 0.211, P = 0.180; Table 1). The model with the highest AICC weight combined these two variables (R2 = 45.46), and, with the exception of one model (geographic distance; model 7, Table 1), all models with ΔAICC within five units of the best AICc model included aridity or vegetation or both (Table 1). Despite the moderate correlation between aridity and geographic distance, aridity was a better predictor of morph frequency divergence with RIV = 0.90 compared with RIV = 0.10 for geographic distance (Table 1).
Linear regressions of the relationship between the frequencies of each morph and the two most informative predictors from the multiple matrix regression models (aridity and vegetation) revealed a consistent pattern of association underpinning the distance-based relationships. The predictor variables explained between 57 and 84 % of the variance in morph frequencies (orange: R2 = 0.842, F2, 5 = 13.31, P = 0.01; orange-yellow: R2 = 0.676, F2, 5 = 5.22, P = 0.06; yellow: R2 = 0.756, F2, 5 = 7.74, P = 0.029; grey: R2 = 0.574, F2, 5 = 3.36, P = 0.119). Specifically, an increase in aridity was associated with an increase in the proportion of orange individuals (standardised coefficient = −0.892, P = 0.004) and a decrease in the proportion of yellow individuals (standardised coefficient = 0.858, P = 0.012; Fig. 3; Additional file 8: Table S6). Bimbowrie Station had a particularly high frequency of orange individuals (Grubbs’ outlier test: G = 2.32, P = 0.002). The relationship between orange morph frequencies and aridity remained significant when Bimbowrie Station was removed from the analysis (standardised coefficient = −0.876, P = 0.024). Additionally, there was a trend of higher frequencies of the orange-yellow morph (standardised coefficient = 0.830, P = 0.023), and lower frequencies of the grey morph (standardised coefficient = −0.744, P = 0.053; Fig. 3; Additional file 8: Table S6), in habitats with a greater proportion cover of vegetation <1 m in height (Fig. 3; Additional file 8: Table S6); however, these relationships were not significant after false discovery rate correction for multiple tests .
We quantified colour morph frequencies in eight populations of C. decresii and examined their relationship with geographic distance, genetic divergence, broad scale environmental gradients and microhabitat features, to elucidate the causes of geographic variation in polymorphism. Our analyses revealed genetic structure among populations, but no relationship between geographic proximity or genetic similarity and similarity in morph frequencies. However, we found a strong association between divergence in morph frequencies and divergence in environmental variables (aridity index and vegetation cover <1 m in height), with the proportion of variance explained being remarkably high for the relatively simple models and modest number of populations considered here (accounting for approximately 45 % of the variance in morph frequencies). This indicates a potentially important role for environmentally driven selection in shaping morph frequencies in C. decresii, as is also observed in multiple other species, including land snails, Cepaea nemoralis , damselflies, Ischnura senegalans , rock pocket mice, Chaetodipus intermedius  and owls . Conversely, environmentally driven selection has been shown to be an unlikely explanation for geographic variation in morph frequencies in several polymorphic lizard species (Podarcis gaigae , Anolis sagrei , Urosaurus ornatus ).
Increasing aridity was associated with a higher proportion of orange males, and fewer yellow males. Broad environmental gradients such as aridity could influence morph frequencies directly by affecting the relative fitness of morphs (e.g. [6, 56]) and/or indirectly by altering intra- and inter-specific interactions (e.g. [57, 58]). Given that discrete colour morphs are often genetically correlated with other traits such as behaviour, physiology and life-history, selection may act on these correlated traits which make certain morphs better adapted to different environments. In C. decresii, colour morphs appear to differ in behavioural strategies. In particular, a study performed at Yourambulla Caves found that orange males react more aggressively to model intruders (M. Yewers and D. Stuart-Fox, unpublished data). Therefore, the orange morph may adopt a dominant behavioural strategy, as is also observed in other polymorphic lizard species (Uta stansburiana, ; Podarcis melisellensis, ; Sceloporus consobrinus, ), and being more aggressive may provide a competitive advantage in more arid habitats where resources are scarce. Additionally, morphs may also differ in other physiological and/or life history traits that are yet to be quantified. For example, colour morphs of the common wall lizard, Podarcis muralis, differ in their immune response to parasites  and sensitivity to stress . Further work is needed to determine the full suite of correlated traits in C. decresii and the complex interaction between morphs and their environment.
We also detected a trend of fewer grey males, and more orange-yellow males, in populations with a greater proportion cover of vegetation. The morphs do not differ in microhabitat preferences, including vegetation cover ; however, vegetation cover may affect signal efficacy by altering light conditions  or obscuring colour signals. For example, in the lizard Anolis cristatellus, dewlap colour is adapted to habitat type, and populations from more mesic habitats have brighter dewlaps . In C. decresii, the grey morph is the least conspicuous, whereas the orange and orange-yellow morphs are the most conspicuous against natural rock backgrounds to the visual system of a conspecific lizard [24, 29]. Furthermore, orange and yellow coloration should maximize conspicuousness in small patches of sunlight where the light environment is rich in long wavelengths . Consequently, orange-yellow throat coloration may be a more effective colour signal than grey throat coloration in highly vegetated habitats, as it is more conspicuous to females and rival males. Vegetation cover could also alter the abundance and species composition of predators, which may interact with selection for signal efficacy to influence the relative fitness of colour morphs [65, 66]. Given the strong relationship between morph frequencies and both aridity and vegetation cover identified here, it would be illuminating to document associated geographic variation in predator species and the relative conspicuousness of morphs in different habitats.
In contrast to the strong relationship between morph frequencies and environmental factors, there was no relationship with genetic divergence or geographic distance. Differences in morph frequencies between geographically adjacent or closely related populations may arise when neighbouring populations differ strongly in their environment, and local conditions select against immigrants (a process termed ‘Isolation By Adaptation’, and a critical component of ecological speciation [15, 67, 68]). Notably, neighbouring populations of C. decresii differed markedly in aridity (e.g. Aroona and Wilpena), vegetation cover (e.g. Mt Remarkable and Devil’s Peak), or both (e.g. Warren Gorge and Devil’s Peak). Furthermore, local directional selection favouring common morphs may be countered by NFDS favouring rare morphs [18–20]. NFDS could also potentially reduce or magnify apparent gene flow by favouring or selecting against immigrants. Specifically, NFDS will favour immigrants belonging to locally rare morphs and select against locally common morphs, which are expected to constitute the largest proportion of migrants between populations with similar morph ranks. For example, Wilpena and Yourambulla Caves have almost identical morph frequencies and are geographically proximate with no obvious barriers to gene flow; however, these populations are genetically distinct. The potential effect of NFDS on morph-specific migration will ultimately depend on the magnitude of differences in morph frequencies between populations and the strength of NFDS in relation to local directional selection. Lastly, morph frequencies in neighbouring populations may be affected by phenotypic or behavioural differences between morphs (e.g. ). For instance, morphs may differ in dispersal, such as in the colour polymorphic side-blotched lizard, Uta stansburiana . Therefore, future research on the C. decresii system should focus on the role of NFDS in maintaining colour polymorphism and whether morph-specific behavioural strategies are associated with dispersal patterns.
Assessing how gene flow and selection interact to shape morph frequencies within a species has broad implications for understanding whether or not polymorphism is likely to promote speciation (reviewed in ). Spatial variation in ecological selection coupled with low or absent gene flow between populations (or selection against migrants) is expected to facilitate population divergence through the eventual fixation of different morphs in different populations . In C. decresii, despite conditions potentially favouring divergence associated with local environmental conditions, populations did not differ in morph composition (i.e. the number and type of morphs). Rare morphs may be maintained by migration between populations; however, gene flow between populations varied, with some populations (e.g. Yourambulla Caves) appearing to have low or absent gene flow. Hence, an alternative explanation is that morphs are under NFDS, a powerful evolutionary force maintaining polymorphism (by favouring rare morphs, including in lizards e.g. Uta stansburiana, ; Anolis sagrei, ), thus countering the influence that local environmental conditions may have on morph frequencies. Assessing morph frequency changes over time would provide insight into whether NFDS also maintains polymorphism in C. decresii.
Morph frequency patterns in C. decresii did not reflect genetic divergence or geographic distance. By contrast, we detected strong relationships between colour morph frequencies and a broad environmental gradient (mean annual aridity index) and microhabitat feature (vegetation cover). Furthermore, while morph frequencies differed significantly across populations, all populations surveyed contained the same four morph types. Selection associated with geographic variation in environmental factors is expected to favour local adaptation and population divergence in morph composition, but may be counteracted by NFDS. Ultimately, it is the interplay between NFDS and divergent natural selection that is likely to be the most important determinant of the potential for speciation in geographically variable polymorphic species.
The data set supporting the results of this article is available in the LabArchives repository, http://dx.doi.org/10.6070/H42805MZ.
Akaike’s information criterion
Atlas of living Australia
Negative frequency-dependent selection
Relative importance of variables
Ford EB. Polymorphism. Heredity. 1945;20:73–88.
Corl A, Davis AR, Kuchta SR, Sinervo B. Selective loss of polymorphic mating types is associated with rapid phenotypic evolution during morphic speciation. Proc Natl Acad Sci U S A. 2010;107(9):4254–9. doi:10.1073/pnas.0909480107.
Hugall AF, Stuart-Fox D. Accelerated speciation in colour polymorphic birds. Nature. 2012;485:631–4. doi:10.1038/nature11050.
McLean CA, Stuart-Fox D. Geographic variation in animal colour polymorphisms and its role in speciation. Biol Rev. 2014;89(4):860–73. doi:10.1111/brv.12083.
Endler JA. Geographic variation, speciation, and clines. Princeton: Princeton University Press; 1977.
Hoekstra HE, Drumm KE, Nachman MW. Ecological genetics of adaptive color polymorphism in pocket mice: geographic variation in selected and neutral genes. Evolution. 2004;58(6):1329–41.
Oxford GS. Genetic drift within a protected polymorphism: Enigmatic variation in color-morph frequencies in the candy-stripe spider, Enoplognatha ovata. Evolution. 2005;59(10):2170–84.
Roulin A, Burri R, Antoniazza S. Owl melanin-based plumage redness is more frequent near than away from the equator: implications on the effect of climate change on biodiversity. Biol J Linn Soc Lond. 2011;102:573–82.
Calsbeek B, Bonvini L, Cox RM. Geographic variation, frequency-dependent selection, and the maintenance of a female-limited polymorphism. Evolution. 2009;64:116–25.
Van Gossum H, Beirinckx K, Forbes MR, Sherratt TN. Do current hypotheses explain continental and seasonal variation in female morph frequencies of the damselfly, Nehalennia irene? Biol J Linn Soc Lond. 2007;90:501–8.
Runemark A, Hansson B, Pafilis P, Valakos ED, Svensson EI. Island biology and morphological divergence of the Skyros wall lizard Podarcis gaigeae: a combined role for local selection and genetic drift on color morph frequency divergence? Evol Biol. 2010;10:269.
Garcia-Ramos G, Kirkpatrick M. Genetic models of adaptation and gene flow in peripheral populations. Evolution. 1997;51:21–8.
Storfer A, Cross J, Rush V, Caruso J. Adaptive coloration and gene flow as a constraint to local adaptation in the streamside salamander, Ambystoma barbouri. Evolution. 1999;53:889–98.
Merilaita S. Habitat heterogeneity, predation and gene flow: colour polymorphism in the isopod, Idotea baltica. Evol Ecol. 2001;15:103–16.
Nosil P. Ecological Speciation. New York: Oxford University Press, Inc.; 2012.
de Leon LF, Bermingham E, Podos J, Hendry AP. Divergence with gene flow as facilitated by ecological differences: within-island variation in Darwin’s finches. Philos Trans R Soc Lond B Biol Sci. 2010;365:1041–52.
Saccheri IJ, Rousset F, Watts PC, Brakefield PM, Cook LM. Selection and gene flow on a diminishing cline of melanic peppered moths. Proc Natl Acad Sci U S A. 2008;105:16212–7.
Svensson EI, Abbott J, Hardling R. Female polymorphism, frequency dependence, and rapid evolutionary dynamics in natural populations. Am Nat. 2005;165(5):568–76.
Takahashi Y, Yoshimura J, Morita S, Watanabe M. Negative frequency-dependent selection in female color polymorphism of a damselfly. Evolution. 2010;64:3620–8.
Sinervo B, Lively CM. The rock-paper-scissors game and the evolution of alternative male strategies. Nature. 1996;380(6571):240–3.
Olendorf R, Rodd FH, Punzalan D, Houde AE, Hurt C, Reznick DN, et al. Frequency-dependent survival in natural guppy populations. Nature. 2006;441(7093):633–6.
McLean CA, Moussalli A, Sass S, Stuart-Fox D. Taxonomic assessment of the Ctenophorus decresii (Reptilia: Agamidae) complex reveals a new species of dragon lizard from western New South Wales. Rec Aus Mus. 2013;65:51–63.
McLean CA, Stuart-Fox D, Moussalli A. Phylogeographic structure, demographic history and morph composition in a colour polymorphic lizard. J Evol Biol. 2014;27(10):2123–37. doi:10.1111/jeb.12464.
Teasdale LC, Stevens M, Stuart-Fox D. Discrete colour polymorphism in the tawny dragon lizard (Ctenophorus decresii) and differences in signal conspicuousness among morphs. J Evol Biol. 2013;26:1035–46.
Rankin K. Genetic and environmental influences on male throat colouration in the colour polymorphic tawny dragon lizard, Ctenophorus decresii. Melbourne, Australia: The University of Melbourne; 2014. master’s thesis.
Stuart-Fox D. Evolution of colour variation and species richness in agamid lizards. Brisbane: Ph.D. thesis. Brisbane, Australia: The University of Queensland; 2003.
Osborne L. Male contest behaviour and information content of signals used by the Australian Tawny Dragon: Ph.D. thesis. Canberra, Australia: Australian National University; 2004.
Gibbons JRH. The hind leg pushup display of the Amphibolurus decresii species complex (Lacertilia: Agamidae). Copeia. 1979;1979(1):29–40.
McLean CA, Moussalli A, Stuart-Fox D. Local adaptation and divergence in colour signal conspicuousness between monomorphic and polymorphic lineages in a lizard. J Evol Biol. 2014;27(12):2654–64. doi:10.1111/jeb.12521.
Stevens M, Parraga CA, Cuthill IC, Partridge JC, Troscianko TS. Using digital photography to study animal coloration. Biol J Linn Soc Lond. 2007;90:211–37.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/. 2010.
Schwartz TS, Warner DA, Beheregaray LB, Olsson M. Microsatellite loci for Australian agamid lizards. Mol Ecol Notes. 2007;7:528–31.
McLean CA, Stuart-Fox D, Moussalli A. Environment, but not genetic divergence, influences geographic variation in colour morph frequencies in a lizard. LabArchives 2015, http://dx.doi.org/10.6070/H42805MZ.
Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Resour. 2004;4:535–8.
Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.
Dieringer D, Schlötterer C. Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003;3(1):167–9.
Bonnet R, Ven De Peer Y. Zt: a software tool for simple and partial Mantel tests. J Stat Softw. 2002;7:1–12.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.
Earl DA, VonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno mothod. Conserv Genet Resour. 2012;4:359–60.
Coulon A, Fitzpatrick JW, Bowman R, Stith BM, Makarewich CA, Stenzler LM, et al. Congruent population structure inferred from dispersal behaviour and intensive genetic surveys of the threatened Florida scrub-jay (Aphelocoma cœrulescens). Mol Ecol. 2008;17:1685–701.
Martín J, López P. Shifts in microhabitat use by the lizard Psammodromus algirus: responses to seasonal changes in vegetation structure. Copeia. 1998;1998:780–6.
Adolph SC. Influence of behavioral thermoregulation on microhabitat use by two Sceloporus lizards. Ecology. 1990;71:315–27.
Lebas N. Microsatellite determination of male reproductive success in a natural population of the territorial ornate dragon lizard, Ctenophorus ornatus. Mol Ecol. 2001;10:193–203.
Goudet J. FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Available from http://www2.unil.ch/popgen/softwares/fstat.htm. 2001.
Legendre P, Fortin M-J. Comparison of the Mantel test and alternative approaches for detecting complex multivariate relationships in the spatial analysis of genetic data. Mol Ecol Resour. 2010;10:831–44.
Wang IJ. Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution. 2013;67:3403–11.
Symonds MRE, Moussalli A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav Ecol Sociobiol. 2001;65:13–21.
Akaike H. Information theory as an extension of the maximum likelihood principle. In: Petrov BN, Csaki F, editors. Proceedings of the Second International Symposium in Information Theory. Budapest: Akademiai Kiado; 1973.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.
Verhoeven KJF, Simonsen KL, McIntyre LM. Implementing false discovery rate control: increasing your power. Oikos. 2005;108:643–7.
Ozgo M. Rapid evolution in unstable habitats: a success story of the polymorphic land snail Cepaea nemoralis (Gastropoda: Pulmonata). Biol J Linn Soc Lond. 2011;102:251–62.
Takahashi Y, Morita S, Yoshimura J, Watanabe M. A geographic cline induced by negative frequency-dependent selection. BMC Evol Biol. 2011;11:256. doi:10.1186/1471-2148-11-256.
Hews DK, Thompson CW, Moore IT, Moore MC. Population frequencies of alternative male phenotypes in tree lizards: geographic variation and common-garden rearing studies. Behav Ecol Sociobiol. 1997;41(6):371–80.
Brisson JA, De Toni DC, Duncan I, Templeton AR. Abdominal pigmentation variation in Drosophila polymorpha: Geographic variation in the trait, and underlying phylogeography. Evolution. 2005;59:1046–59.
Gosden TP, Svensson EI. Density-dependent male mating harassment, female resistance, and male mimicry. Am Nat. 2009;173(6):709–21.
Gosden TP, Stoks R, Svensson E. Range limits, large-scale biogeographic variation, and localized evolutionary dynamics in a polymorphic damselfly. Biol J Linn Soc Lond. 2011;102:775–85.
Huyghe K, Vanhooydonck B, Herrel A, Tadic Z, Van Damme R. Morphology, performance, behavior and ecology of three color morphs in males of the lizard Podarcis melisellensis. Integr Comp Biol. 2007;47(2):211–20.
Rand M. Courtship and aggressive behavior in male lizards exhibiting two different sexual colorations. Am Zool. 1988;28:153A.
Sacchi R, Rubolini D, Gentilli A, Pupin F, Razzetti E, Scali S, et al. Morph-specific immunity in male Podarcis muralis. Amphibia-Reptilia. 2007;28:408–12.
Galeotti P, P-R D, Sacchi R, Gentilli A, Pupin F, Rubolini D, et al. Sex-, morph- and size-specific susceptibility to stress measured by haematological variable in captive common wall lizard Podarcis muralis. Comp Biochem Physiol. 2010;157:354–63.
Endler JA. The color of light in forests and its implications. Ecol Monogr. 1993;63:1–27.
Leal M, Fleishman LJ. Differences in visual signal design and detectability between allopatric population of Anolis lizards. Am Nat. 2004;163:26–39.
Ahnesjo J, Forsman A. Differential habitat selection by pygmy grasshopper color morphs; interactive effects of temperature and predator avoidance. Evol Ecol. 2006;20:335–257.
Macedonia JM, Brandt Y, Clark DL. Sexual dichromatism and differential conspicuousness in two populations of the common collared lizard (Crotaphytus collaris) from Utah and New Mexico, USA. Biol J Linn Soc Lond. 2002;77:67–85.
Hendry AP. Selection against migrants contributes to the rapid evolution of ecologically dependent reproductive isolation. Evol Ecol Res. 2004;6:1219–36.
Shafer ABA, Wolf JBW. Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology. Ecol Lett. 2013;16:940–50. doi:10.1111/ele.12120.
Socha R, Zemek R. Wing morph-related differences in the walking pattern and dispersal in a flightless bug, Pyrrhocoris apterus (Heteroptera). Oikos. 2003;100:35–42.
Sinervo B, Clobert J. Morphs, dispersal behavior, genetic similarity, and the evolution of cooperation. Science. 2003;300(5627):1949–51.
We are grateful to Maik Fiedel, Adam Elliott, Bryant Turffs, Tom McLean and Danial Nordinson for assistance in the field, and thank Peter Watkins and Leah and Richard Khoe for their hospitality. Fieldwork was conducted with approval from the South Australia Department of Environment, Water and Natural Resources (permit: E25861-1) and ethics approval was obtained from the South Australian Wildlife Ethics Committee (18/2010) and The University of Melbourne Animal Ethics Committee (1011760.1). Research funding was provided by the Australian Research Council to DS-F (DP1092908), and Nature Foundation South Australia and the Holsworth Wildlife Research Endowment to CM.
The authors declare that they have no competing interests.
CM designed the study, conducted the fieldwork, performed the laboratory work, analysed and interpreted the data and wrote the manuscript. DS-F conceived of the study, contributed to data interpretation and edited the manuscript. AM helped design the study, contributed to data analysis and interpretation and edited the manuscript. All authors gave final approval for publication.
Characteristics of eight microsatellite loci, screened across thirty individuals from eight populations of Ctenophorus decresii. Significant linkage disequilibrium was detected between Ctde08 and Ctde45 across all populations. Statistically significant values are bold and italicised, * indicates possible null alleles. (PDF 197 kb)
Mean annual aridity index and topographic relief across the range of northern C. decresii with values for Aroona (A), Bimbowrie Station (BS), Devil’s Peak (DP), Mt Remarkable (MR), Telowie Gorge (TG), Wilpena (W), Warren Gorge (WG) and Yourambulla Caves (YC). (PDF 311 kb)
Correlations among geographic distance (distance), mean annual aridity index (aridity), topographic relief (topography), proportion cover of vegetation <1 m in height (vegetation), and proportion cover of rocks >50 cm in diameter. Pearson correlation coefficients are below the diagonal and corresponding P values are above the diagonal. Statistically significant relationships after false discovery rate correction for multiple tests are bold and italicised . (PDF 85 kb)
Relative frequencies of male colour morphs at Telowie Gorge and Warren Gorge in 2010 and 2011. (PDF 84 kb)
Test for significant differences in male morph frequencies among all pairs of populations in a chi-square test (A: Aroona, W: Wilpena, YC: Yourambulla Caves, WG: Warren Gorge, DP: Devil’s Peak, BS: Bimbowrie Station, MR: Mt Remarkable, TG: Telowie Gorge). Statistically significant values after false discovery rate correction for multiple tests are bold and italicised . (PDF 146 kb)
Pairwise corrected FST between Aroona (A), Wilpena (W), Yourambulla Caves (YC), Warren Gorge (WG), Devil’s Peak (DP), Bimbowrie Station (BS), Mt Remarkable (MR) and Telowie Gorge (TG) calculated from eight microsatellite loci (below the diagonal) and P values after sequential Bonferonni correction for multiple tests (above the diagonal). (PDF 139 kb)
Figures showing weak relationships between pairwise FST and geographic distance (a), and pairwise FST and divergence in colour morph composition (b). (PDF 152 kb)
Results of four separate multiple regression analyses with the frequency of each colour morph (orange, orange-yellow, yellow and grey) as the dependent variables and mean annual aridity index (aridity) and proportion cover of vegetation <1 m in height (vegetation) as the predictor variables. Statistically significant relationships after false discovery rate correction for multiple tests are bold and italicised . (PDF 88 kb)