Patterns and drivers of intraspecific variation in avian life history along elevational gradients: a meta‐analysis

Elevational gradients provide powerful natural systems for testing hypotheses regarding the role of environmental variation in the evolution of life‐history strategies. Case studies have revealed shifts towards slower life histories in organisms living at high elevations yet no synthetic analyses exist of elevational variation in life‐history traits for major vertebrate clades. We examined (i) how life‐history traits change with elevation in paired populations of bird species worldwide, and (ii) which biotic and abiotic factors drive elevational shifts in life history. Using three analytical methods, we found that fecundity declined at higher elevations due to smaller clutches and fewer reproductive attempts per year. By contrast, elevational differences in traits associated with parental investment or survival varied among studies. High‐elevation populations had shorter and later breeding seasons, but longer developmental periods implying that temporal constraints contribute to reduced fecundity. Analyses of clutch size data, the trait for which we had the largest number of population comparisons, indicated no evidence that phylogenetic history constrained species‐level plasticity in trait variation associated with elevational gradients. The magnitude of elevational shifts in life‐history traits were largely unrelated to geographic (altitude, latitude), intrinsic (body mass, migratory status), or habitat covariates. Meta‐population structure, methodological issues associated with estimating survival, or processes shaping range boundaries could potentially explain the nature of elevational shifts in life‐history traits evident in this data set. We identify a new risk factor for montane populations in changing climates: low fecundity will result in lower reproductive potential to recover from perturbations, especially as fewer than half of the species experienced higher survival at higher elevations.     

A unified framework for diversity gradients: the adaptive trait continuum

Aim Adaptive trait continua are axes of covariation observed in multivariate trait data for a given taxonomic group. These continua quantify and summarize life‐history variation at the inter‐specific level in multi‐specific assemblages. Here we examine whether trait continua can provide a useful framework to link life‐history variation with demographic and evolutionary processes in species richness gradients. Taking an altitudinal species richness gradient for Mediterranean butterflies as a study case, we examined a suite of traits (larval diet breadth, adult phenology, dispersal capacity and wing length) and species‐specific habitat measures (temperature and aridity breadth). We tested whether traits and species‐specific habitat measures tend to co‐vary, whether they are phylogenetically conserved, and whether they are able to explain species distributions and spatial genetic variation in a large number of butterfly assemblages. Location Catalonia, Spain. Methods We formulated predictions associated with species richness gradients and adaptive trait continua. We applied principal components analyses (PCAs), structural equation modelling and phylogenetic generalized least squares models. Results We found that traits and species‐specific habitat measures covaried along a main PCA axis, ranging from multivoltine trophic generalists with high dispersal capacity to univoltine (i.e. one generation per year), trophic specialist species with low dispersal capacity. This trait continuum was closely associated with the observed distributions along the altitudinal gradient and predicted inter‐specific differences in patterns of spatial genetic variability (F_ST and genetic distances), population responses to the impacts of global change and local turnover dynamics. Main conclusions The adaptive trait continuum of Mediterranean butterflies provides an integrative and mechanistic framework to: (1) analyse geographical gradients in species richness, (2) explain inter‐specific differences in population abundances, spatial distributions and demographic trends, (3) explain inter‐specific differences in patterns of genetic variation (F_ST and genetic distances), and (4) study specialist–generalist life‐history transitions frequently involved in butterfly diversification processes.     

Evolving possibilities: postembryonic axial elongation in salamanders with biphasic (Eurcyea cirrigera,Eurycea longicauda,Eurycea quadridigitata) and paedomorphic life cycles (Eurycea nana and Ambystoma mexicanum)

Vaglia, JL., White, K, and Case, A. 2012. Evolving possibilities: postembryonic axial elongation in salamanders with biphasic (Eurcyea cirrigera, Eurycea longicauda, Eurycea quadridigitata) and paedomorphic life cycles (Eurycea nana and Ambystoma mexicanum). —Acta Zoologica (Stockholm) 93 : 2–13. Typically, the number of vertebrae an organism will have postembryonically is determined during embryogenesis via the development of paired somites. Our research investigates the phenomenon of postembryonic vertebral addition in salamander tails. We describe body and tail growth and patterns of postsacral vertebral addition and elongation in context with caudal morphology for four plethodontids (Eurycea) and one ambystomatid. Eurycea nana and Ambystoma mexicanum have paedomorphic life cycles; Eurcyea cirrigera, Eurycea longicauda and Eurycea quadridigitata are biphasic. Specimens were collected, borrowed and/or purchased, and cleared and stained for bone and cartilage. Data collected include snout‐vent length (SVL), tail length (TL), vertebral counts and centrum lengths. Eurycea species with biphasic life cycles had TLs that surpassed SVL following metamorphosis. Tails in paedomorphic species elongated but rarely exceeded body length. Larger TLs were associated with more vertebrae and longer vertebrae in all species. We observed that rates of postsacral vertebral addition varied little amongst species. Regional variation along the tail becomes prominent following metamorphosis in biphasic developers. In all species, vertebrae in the posterior one‐half of the tail taper towards the tip. We suggest that a developmental link might exist between the ability to continually add vertebrae and regeneration in salamanders.     

EVOLUTION OF PAEDOMORPHOSIS IN PLETHODONTID SALAMANDERS: ECOLOGICAL CORRELATES AND RE‐EVOLUTION OF METAMORPHOSIS

Life‐history modes can profoundly impact the biology of a species, and a classic example is the dichotomy between metamorphic (biphasic) and paedomorphic (permanently aquatic) life‐history strategies in salamanders. However, despite centuries of research on this system, several basic questions about the evolution of paedomorphosis in salamanders have not been addressed. Here, we use a nearly comprehensive, time‐calibrated phylogeny of spelerpine plethodontids to reconstruct the evolution of paedomorphosis and to test if paedomorphosis is (1) reversible; (2) associated with living in caves; (3) associated with relatively dry climatic conditions on the surface; and (4) correlated with limited range size and geographic dispersal. We find that paedomorphosis arose multiple times in spelerpines. We also find evidence for re‐evolution of metamorphosis after several million years of paedomorphosis in a lineage of Eurycea from the Edwards Plateau region of Texas. We also show for the first time using phylogenetic comparative methods that paedomorphosis is highly correlated with cave‐dwelling, arid surface environments, and small geographic range sizes, providing insights into both the causes and consequences of this major life history transition.     

Contrasting latitudinal patterns of life‐history divergence in two genera of new world thrushes (Turdinae)

Several long‐standing hypotheses have been proposed to explain latitudinal patterns of life‐history strategies. Here, we test predictions of four such hypotheses (seasonality, food limitation, nest predation and adult survival probability) by examining life‐history traits and age‐specific mortality rates of several species of thrushes (Turdinae) based on field studies at temperate and tropical sites and data gathered from the literature. Thrushes in the genus Catharus showed the typical pattern of slower life‐history strategies in the tropics while co‐occuring Turdus thrushes differed much less across latitudes. Seasonality is a broadly accepted hypothesis for latitudinal patterns, but the lack of concordance in latitudinal patterns between co‐existing genera that experience the same seasonal patterns suggests seasonality cannot fully explain latitudinal trait variation in thrushes. Nest‐predation also could not explain patterns based on our field data and literature data for these two genera. Total feeding rates were similar, and per‐nestling feeding rates were higher at tropical latitudes in both genera, suggesting food limitation does not explain trait differences in thrushes. Latitudinal patterns of life histories in these two genera were closely associated with adult survival probability. Thus, our data suggest that environmental influences on adult survival probability may play a particularly strong role in shaping latitudinal patterns of life‐history traits.     

Natural selection drives clinal life history patterns in the perennial sunflower species, Helianthus maximiliani

In plants, ecologically important life history traits often display clinal patterns of population divergence. Such patterns can provide strong evidence for spatially varying selection across environmental gradients but also may result from nonselective processes, such as genetic drift, population bottlenecks and spatially restricted gene flow. Comparison of population differentiation in quantitative traits (measured as Q(ST) ) with neutral molecular markers (measured as F(ST) ) provides a useful tool for understanding the relative importance of adaptive and nonadaptive processes in the formation and maintenance of clinal variation. Here, we demonstrate the existence of geographic variation in key life history traits in the diploid perennial sunflower species Helianthus maximiliani across a broad latitudinal transect in North America. Strong population differentiation was found for days to flowering, growth rate and multiple size-related traits. Differentiation in these traits greatly exceeds neutral predictions, as determined both by partial Mantel tests and by comparisons of global Q(ST) values with theoretical F(ST) distributions. These findings indicate that clinal variation in these life history traits likely results from local adaptation driven by spatially heterogeneous environments.     

An invasive species imposes selection on life‐history traits of a native frog

As well as their direct ecological impacts on native taxa, invasive species can impose selection on phenotypic attributes (morphology, physiology, behaviour, etc.) of the native fauna. In anurans, body size at metamorphosis is a critical life‐history trait: for most challenges faced by post‐metamorphic anurans, larger size at metamorphosis probably enhances survival. However, our studies on Australian frogs (Limnodynastes convexiusculus) show that this pattern can be reversed by the arrival of an invasive species. When metamorph frogs first encounter invasive cane toads (Bufo marinus), they try to eat the toxic invader and, if they are able to do so, are likely to die from poisoning. Because frogs are gape‐limited predators, small metamorphs cannot ingest a toad and thus survive long enough to disperse away from the natal pond (and thus from potentially deadly toads). These data show that larger size at metamorphosis can reduce rather than increase anuran survival rates, because larger metamorphs are more easily able to ingest (and thus be poisoned by) metamorph cane toads. Our results suggest that patterns of selection on life‐history traits of native taxa (such as size and age at metamorphosis, seasonal timing of breeding and duration of pondside aggregation prior to dispersal) can be modified by the arrival of an invasive species. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 329–336.     

CLINES IN CUTICULAR HYDROCARBONS IN TWO DROSOPHILA SPECIES WITH INDEPENDENT POPULATION HISTORIES

We took a comparative approach utilizing clines to investigate the extent to which natural selection may have shaped population divergence in cuticular hydrocarbons (CHCs) that are also under sexual selection in Drosophila. We detected the presence of CHC clines along a latitudinal gradient on the east coast of Australia in two fly species with independent phylogenetic and population histories, suggesting adaptation to shared abiotic factors. For both species, significant associations were detected between clinal variation in CHCs and temperature variation along the gradient, suggesting temperature maxima as a candidate abiotic factor shaping CHC variation among populations. However, rainfall and humidity correlated with CHC variation to differing extents in the two species, suggesting that response to these abiotic factors may vary in a species‐specific manner. Our results suggest that natural selection, in addition to sexual selection, plays a significant role in structuring among‐population variation in sexually selected traits in Drosophila.     

Picturing local adaptation: Spectral and structural traits from drone remote sensing reveal clinal responses to climate transfer in common-garden trials of interior spruce (Picea engelmannii × glauca)

Common-garden trials of forest trees provide phenotype data used to assess growth and local adaptation; this information is foundational to tree breeding programs, genecology, and gene conservation. As jurisdictions consider assisted migration strategies to match populations to suitable climates, in situ progeny and provenance trials provide experimental evidence of adaptive responses to climate change. We used drone technology, multispectral imaging, and digital aerial photogrammetry to quantify spectral traits related to stress, photosynthesis, and carotenoids, and structural traits describing crown height, size, and complexity at six climatically disparate common-garden trials of interior spruce (Picea engelmannii × glauca) in western Canada. Through principal component analysis, we identified key components of climate related to temperature, moisture, and elevational gradients. Phenotypic clines in remotely sensed traits were analyzed as trait correlations with provenance climate transfer distances along principal components (PCs). We used traits showing clinal variation to model best linear unbiased predictions for tree height (R² = .98–.99, root mean square error [RMSE] = 0.06–0.10 m) and diameter at breast height (DBH, R² = .71–.97, RMSE = 2.57–3.80 mm) and generated multivariate climate transfer functions with the model predictions. Significant (p < .05) clines were present for spectral traits at all sites along all PCs. Spectral traits showed stronger clinal variation than structural traits along temperature and elevational gradients and along moisture gradients at wet, coastal sites, but not at dry, interior sites. Spectral traits may capture patterns of local adaptation to temperature and montane growing seasons which are distinct from moisture-limited patterns in stem growth. This work demonstrates that multispectral indices improve the assessment of local adaptation and that spectral and structural traits from drone remote sensing produce reliable proxies for ground-measured height and DBH. This phenotyping framework contributes to the analysis of common-garden trials towards a mechanistic understanding of local adaptation to climate.     

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min^?1) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life‐stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.     

Land‐use legacies in a central Appalachian forest: differential response of trees and herbs to historic agricultural practices

Question: Are contemporary herb and tree patterns explained by historic land use practices? If so, are observed vegetation patterns associated with life‐history characteristics, soil properties, or other environmental variables? Location: Southeastern Ohio, USA. Methods: Using archival records, currently forested sites were identified with distinct land use histories: cultivated, pasture (but not plowed), and reference sites which appear to have never been cleared. Trees were recorded by size and species on twenty 20 m × 20 m plots; percent cover was estimated for each herb species in nested 10 m × 10 m plots. Environmental characteristics were noted, and soil samples analysed for nutrient availability and organic matter. Nonmetric multidimensional scaling ordination was performed separately on both tree and herb datasets to graphically characterize community composition among plots. Life‐history traits were investigated to explain observed compositional differences. Results: Vegetation patterns were explained by current environmental gradients, especially by land‐use history. Cultivated and pasture sites had similar tree composition, distinct from reference sites. Herb composition of pasture and reference sites was similar and distinct from cultivated sites, suggesting the ‘tenacity’ of some forest herbs on formerly cleared sites. Tilling removes rhizomatous species, and disfavors species with unassisted dispersal. These life‐history traits were underrepresented on cultivated sites, although ant‐dispersed species were not. Conclusions: Historic land‐use practices accounted for as much variation in species composition as environmental gradients. Furthermore, trees and herbs responded differently to past land‐use practices. Life‐history traits of individual species interact with the nature of disturbance to influence community composition.     

HABITAT DURATION,LENGTH OF LARVAL PERIOD,AND THE EVOLUTION OF A COMPLEX LIFE CYCLE OF A SALAMANDER,AMBYSTOMA TEXANUM

In many organisms, genotypic selection may be a less effective means of adapting to unpredictable environments than is selection for phenotypic plasticity. To determine whether genotypic selection is important in the evolution of complex life cycles of amphibians that breed in seasonally ephemeral habitats, we examined whether mortality risk from habitat drying in natural populations of small-mouthed salamanders (Ambystoma texanum) corresponded to length of larval period when larvae from the same populations were grown in a common laboratory environment. Comparisons were made at two levels of organization within the species: 1) among geographic races that are under strongly divergent selection regimes associated with the use of pond and stream habitats and 2) among populations within races that use the same types of breeding habitats. Morphological evidence indicates that stream-breeding A. texanum evolved from pond-breeding populations that recently colonized streams. Larvae in streams incur heavy mortality from stream drying, so the upper bound on length of larval period is currently set by the seasonal duration of breeding sites. We hypothesized that selection would reduce length of larval period of pond-breeders that colonize streams if their larval periods are inherently longer than those of stream-breeders. The results of laboratory experiments support this hypothesis. When grown individually in a common environment, larvae from stream populations had significantly shorter larval periods than larvae from pond populations. Within races, however, length of larval period did not correlate significantly with seasonal duration of breeding sites. When males of both races were crossed to a single pond female, offspring of stream males had significantly shorter larval periods than offspring of pond males. Collectively, these data suggest that differences in complex life cycles among pond and stream-breeders are due to genotypic selection related to mortality from habitat drying. Stream larvae in the common-environment experiment were significantly smaller at metamorphosis than pond larvae. Yet, the evolution of metamorphic size cannot be explained readily by direct selection: there are no intuitively obvious advantages of being relatively small at metamorphosis in streams. A positive phenotypic correlation was observed between size at metamorphosis and length of larval period in most laboratory populations. A positive additive genetic correlation between these traits was demonstrated recently in another amphibian. Thus, we suspect that metamorphic size of stream-breeders evolved indirectly as a consequence of selection to shorten length of larval period.     

Developmental acclimation affects clinal variation in stress resistance traits in Drosophila buzzatii

Patterns of clinal genetic variation in Drosophila are often characterized after rearing at constant temperatures. However, clinal patterns might change after acclimation if populations differ in their plastic response to fluctuating environments. We studied longevity, starvation and heat knock‐down resistance after development at either constant or fluctuating temperatures in nine Drosophila buzzatii populations collected along an altitudinal gradient in Tenerife, Spain. Flies that developed at fluctuating temperatures had higher stress resistance despite experiencing a slightly lower average temperature than those at constant temperatures. Genetic variation along the gradient was found in both stress‐resistance traits. Because Q_ST values greatly exceeded F_ST values, genetic drift could not explain this diversification. In general, differences among populations were larger after rearing at fluctuating temperatures, especially in heat knock‐down, for which clinal patterns disappeared when flies were reared at constant temperatures. This result emphasizes the importance of determining whether populations originating from different environments differ in their plastic responses to stress.     

INTRASPECIFIC LIFE HISTORY VARIATION IN A POND SNAIL: THE ROLES OF POPULATION DIVERGENCE AND PHENOTYPIC PLASTICITY

Two approaches were used to determine the degree of divergence in life histories among populations of the pond snail, Lymnaea elodes. Juvenile snails were reciprocally transferred between ponds differing in permanence and productivity, and the resulting variation in life history traits was recorded. In a second experiment, parents and their offspring from both a vernal and a permanent pond population were reared in the same pond. Proximal factors had by far the greatest effects on life history traits in the transfer experiment, with snails reared in a more productive pond showing earlier reproduction at a larger size, higher fecundity, and longer life cycle length. Snails from the more uncertain pond in terms of drying date did reproduce at an earlier age and smaller size and grew less in each pond. However, these population differences, for the most part, disappeared when snails were reared for two generations in the same environment. Much of the intraspecific variation in life histories seen in this species must therefore be considered the result of phenotypic plasticity. I argue that the plasticity in life histories itself may be adaptive to this inhabitant of unpredictable, vernal ponds.     

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an F_ST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.     

Breeding site characteristics regulating life history traits of the brown tree frog, Litoria ewingii

Past research has determined the habitat requirements of amphibian species predominantly from presence/absence studies. This study tested the hypothesis that relationships between breeding site habitat components, life history traits and fitness may provide a higher resolution of biological data relating to the habitat requirements of amphibian species. We tested this novel approach by using Litoria ewingii as our model species. We correlated larval and metamorph life history traits with habitat variables at 28 small to medium sized ponds within a commercially logged forest in southern Tasmania, Australia. To avoid larval mortality due to pond desiccation, L. ewingii laid eggs and metamorphosed earlier in smaller ponds. Snout vent length at metamorphosis increased with elevation and metamorphosis was earlier in less shaded ponds. Breeding ponds that maximised the fitness of L. ewingii were higher elevation ponds with reduced shading, steeper bank slopes and reduced pond isolation. The findings of the study equip land managers with a greater ecological understanding of ecosystem function in relation to specific species. The methodological approach has broad application to conservation biology where an awareness of the specific habitat requirements of amphibians is critical to successful ecosystem management.     

Effects of conditionally expressed phenotypes and environment on amphibian dispersal in nature

Individuals vary greatly in the distance they disperse, and in doing so, strongly affect ecological and evolutionary processes. Dispersal, when viewed as a component of phenotype, can be affected independently or jointly by environment. However, among taxa with complex life cycles that occupy different habitats over ontogeny, the effects of environment on dispersal and the interaction between environment and phenotype remains poorly understood. Here, we conducted a field experiment to measure how dispersal distance was affected by phenotype, environment experienced before and after metamorphosis, and their interaction. We manipulated the environment encountered by a pond‐breeding salamander Ambystoma annulatum during the aquatic larval stage and again as dispersing terrestrial juveniles. After assaying juvenile phenotype (exploration behavior, body condition, and morphology), we then measured the initial distance dispersed by juveniles. The distance moved by dispersing salamanders was affected by attributes of both larval and juvenile habitat, with salamanders that encountered low quality habitat in either life stage moving the farthest. However, we did not find support for an interactive effect of phenotype and environment affecting the distance moved by dispersers. Interestingly, exploration behavior explained the distance moved by philopatric animals but not dispersing ones. Our findings indicate that the environment experienced before metamorphosis can affect juvenile dispersal behavior, and demonstrates the need to consider dispersal in species with complex life cycles to understand the coupling between local and regional population dynamics.     

Amphibians live longer at higher altitudes but not at higher latitudes

Why and how organisms differ in life‐history strategies across their range is a long‐standing topic of interest to evolutionary ecologists. Although many studies have addressed this issue for several life‐history traits, such as body size and clutch size, very few have been made for some others traits, including longevity. In the present study, we performed a comparative study aiming to develop general patterns of geographical variation in longevity of urodele and anuran amphibians using published information on demographic age derived from skeletochronology. We conducted within‐species meta‐analyses using datasets of two (ten urodele and 12 anuran species) and multiple (two urodele and nine anuran species) spatially‐separated populations and found that maturation, mean, and maximum age all increased with altitude but not with latitude in each sex of both amphibian groups. This geographical pattern held true across 33 urodele and 86 anuran species at common body sizes, independent of phylogeny. It is likely that metabolic rate, reproductive investment, and mortality risk, which are the key factors that affect longevity as suggested by ageing theory, vary systemically along altitudinal gradients but not along latitudinal gradients. The evolutionary causes behind these puzzling patterns deserve further investigation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 623–632.     

Patterns of differential introgression in a salamander hybrid zone: inferences from genetic data and ecological niche modelling

Hybrid zones have yielded considerable insight into many evolutionary processes, including speciation and the maintenance of species boundaries. Presented here are analyses from a hybrid zone that occurs among three salamanders –Plethodon jordani, Plethodon metcalfi and Plethodon teyahalee– from the southern Appalachian Mountains. Using a novel statistical approach for analysis of non‐clinal, multispecies hybrid zones, we examined spatial patterns of variation at four markers: single‐nucleotide polymorphisms (SNPs) located in the mtDNA ND2 gene and the nuclear DNA ILF3 gene, and the morphological markers of red cheek pigmentation and white flecks. Concordance of the ILF3 marker and both morphological markers across four transects is observed. In three of the four transects, however, the pattern of mtDNA is discordant from all other markers, with a higher representation of P. metcalfi mtDNA in the northern and lower elevation localities than is expected given the ILF3 marker and morphology. To explore whether climate plays a role in the position of the hybrid zone, we created ecological niche models for P. jordani and P. metcalfi. Modelling results suggest that hybrid zone position is not determined by steep gradients in climatic suitability for either species. Instead, the hybrid zone lies in a climatically homogenous region that is broadly suitable for both P. jordani and P. metcalfi. We discuss various selective (natural selection associated with climate) and behavioural processes (sex‐biased dispersal, asymmetric reproductive isolation) that might explain the discordance in the extent to which mtDNA and nuclear DNA and colour‐pattern traits have moved across this hybrid zone.