- Open Access
The role of thermal niche selection in maintenance of a colour polymorphism in redback salamanders (Plethodon cinereus)
Frontiers in Zoologyvolume 3, Article number: 10 (2006)
In eastern North America two common colour morphs exist in most populations of redback salamanders (Plethodon cinereus). Previous studies have indicated that the different morphs may be adapted to different thermal niches and the morphological variation has been linked to standard metabolic rate at 15°C in one population of P. cinereus. It has therefore been hypothesized that a correlated response to selection on metabolic rate across thermal niches maintains the colour polymorphism in P. cinereus. This study tests that hypothesis.
We found that the two colour morphs do sometimes differ in their maintenance metabolic rate (MMR) profiles, but that the pattern is not consistent across populations or seasons. We also found that when MMR profiles differ between morphs those differences do not indicate that distinct niches exist. Field censuses showed that the two colour morphs are sometimes found at different substrate temperatures and that this difference is also dependent on census location and season.
While these morphs sometimes differ in their maintenance energy expenditures, the differences in MMR profile in this study are not consistent with maintenance of the polymorphism via a simple correlated response to selection across multiple niches. When present, differences in MMR profile do not indicate the existence of multiple thermal niches that consistently mirror colour polymorphism. We suggest that while a relationship between colour morph and thermal niche selection appears to exist it is neither simple nor consistent.
Maintenance of variation in populations is a fundamental topic in evolutionary ecology . Mechanisms that maintain variation can be selective, non-selective or a mixture of both. Previous studies have identified factors that maintain colour polymorphisms in many taxa (e.g. [2–6]). However, the factors explaining the coexistence of two common morphs of the Redback Salamander,Plethodon cinereus (Plethodontidae) remain poorly understood despite considerable study [7–11].
Natural selection can maintain variation in numerous ways [3, 12]. For instance variation can be maintained across multiple niches where different variants are favoured in different microhabitats or at different times [13–15]. Predation, parasitism, competition or sexual selection could result in the maintenance of variation through frequency or density-dependent selection [16–19]. Non-selective processes could also maintain variation. Frequencies of selectively neutral polymorphisms will fluctuate across populations via genetic drift and will generally lead to fixation in some populations . If neutral variation is pleiotropically or chromosomally linked to traits under selection, then the neutral variation could also be maintained via selection on those traits [18, 21].
P. cinereus are fully terrestrial, lungless, woodland salamanders, which exist as three relatively distinct colour morphs, of which two are commonly found. The striped (redback) morph has a reddish dorsal stripe that extends from the base of the head to the tail, darkly pigmented sides, and a mottled black and white venter. The black (leadback or unstriped) morph is identical to the striped morph except that it has a darkly pigmented dorsum. This variation is heritable, appears to be polygenic in nature [9, 22], and is not known to be environmentally plastic. While the variation is polygenic, it is not quantitative since the unstriped morph is qualitatively distinguished by changes in the cellular architecture of chromatophores including a complete loss of erythrophores .
The two common colour morphs exist in most P. cinereus populations in eastern North America . The relative rarity of monomorphic populations throughout most of the range suggests that genetic drift alone is unlikely to maintain this widespread polymorphism. Previous studies have found that the frequency of the two morphs sometimes varies with climate either geographically , or temporally [8, 10], indicating that the different morphs may be adapted to different thermal niches. Colour phenotype has also been linked to standard metabolic rate (SMR) at 15°C in one population of P. cinereus . Because SMR is an estimate of the maintenance costs that represent a large fraction of an ectotherm's energy budget [25, 26] such linkage of colour and SMR provides a mechanism by which the polymorphism might be maintained. This mechanism assumes that differences in metabolic rate translate to differences in fitness since the energy available for survival and reproduction depends on the energy remaining after maintenance costs are paid [25, 27].
The purpose of this study was to test a current hypothesis that the colour polymorphism in P. cinereus could be maintained via a correlated response to thermal niche selection on metabolic rate . Although a link between metabolic rate and colour morph is known to exist in one population , demonstrating the existence of distinct thermal niches requires that each morph specialize on a particular temperature. Only then would thermal niche specialization maintain the polymorphism. This study tested for differences in optimal thermal niche between the morphs by determining their maintenance metabolic rates (MMR) across a range of typical field temperatures (MMR profile) in two different populations. It further tested for differences between the two colour morphs in the temperatures of the microhabitats they inhabit. If thermal niche specialization maintains the polymorphism, the morphs should display distinct thermal optima. Furthermore, in the field, the morphs should segregate into thermal niches corresponding to their particular thermal optima.
Maintenance Metabolic Rate
Metabolic rate data, log (ln) transformed to linearize the relationship between oxygen consumption and mass [28, 29], were examined using repeated-measures MANCOVA. This allowed us to test for significant influence of location, season, and phenotype on MMR and on the response profile of MMR across a thermal range. Log transformed mass was included as a covariate to appropriately adjust for variation in the size of animals .
In the repeated measures analysis 'within individual' variation in MMR indicates variation in maintenance costs across temperatures (i.e. variation in MMR profile). Phenotype was a significant source of variation in MMR profile (temperature X phenotype effect in Table 1) indicting that phenotypes vary in their response to temperature (see Figure 1). Location was also a significant source of variation in the MMR profile primarily through its significant interaction with phenotype (temperature X location X phenotype effect in Table 1).
Between individuals comparisons in the repeated measures analysis tested for differences in mean MMR across all temperatures. Season was the only significant source of variation in MMR between individuals (Table 2). This indicates that the MMR profile in the autumn (Figure 2) was shifted up relative to the summer. Differences between phenotypes or locations were not significant (Table 2).
In independent contrasts for each population at each temperature, phenotype was a significant source of variation only in salamanders from VPA at 10°C (p < 0.01). The significant phenotypic influence on MMR profile (Table 1) is attributable in part to this effect (Figure 1). The mean whole body MMR (O2 consumption ml/h) for each morph within each location and season is listed in Table 3. Estimates of MMR (Table 3) are similar to reported estimates that have been standardized to ml/h O2 consumption  for temperatures in the 10°C, 15°C and 20°C ranges .
Field distribution across thermal environments
Seasonal change was operationally defined by the point of temperature stabilization after a rapid and consistent decline in surface and substrate temperature across both sites for several weeks beginning 5 September and ending 19 September. Location, season and the interaction between location and season were significant sources of variation (3 way ANOVA) in the substrate temperature where the salamanders were found (Table 4). Phenotype was a significant source of variation in substrate temperature through its three-way interaction with season and location (Table 4). This significant interaction appears to be attributable in part to the marginally significant difference between the morphs from VPA in the autumn (p = 0.079; t-test). No differences between the morphs were detected in the other three comparisons (p > 0.35 in all cases; t-test). Figure 3 summarises the mean substrate temperature at which each morph within each location and season was found.
If morphs differ in their MMR profiles, they may also differ in their net allocatable energy available for growth and reproduction, and hence in their fitness[21, 25, 31]. This in turn could constitute a mechanism for the maintenance of a colour polymorphism. We found that colour phenotype significantly contributes to variation in MMR profile (Table 1). This pattern of colour contributing to variation in metabolic rate and, in particular, the pattern of the black morph from VPA having a significantly lower metabolic rate than the striped morph at 10°C (Figure 1b) agrees with a study by Moreno , who found that the black morph had a lower SMR than the striped morph at 15°C. For thermal niche specialization to directly maintain the polymorphism however, both morphs must have a niche where they are favoured, and our study did not find evidence of this. Within the range of temperatures that we studied we did not find a consistent pattern across populations in either metabolic response or in thermal niche choice.
The data collected in this study provided the ability to differentiate morphs with even subtle differences in metabolic rate and thermal profile. While differences between morph could be detected we did not find evidence that these differences support a hypothesis that morphs are adapted for different thermal optima. We only found significant differences in MMR between the morphs at 10°C for animals from VPA and in that population there was no evidence of separate thermal niches. For separate niches to exist each morph would have to have a thermal niche in which it is favoured. This pattern of response was evident in neither the O'Neil nor VPA population. Furthermore, given that the influence of phenotype on MMR profile varied significantly between locations, thermal niche specialization is unlikely to be a uniform force directly maintaining the polymorphism in these populations. The existence of significant variation between morphs in MMR profile across temperature (Table 1), does however indicate a linkage of some sort between morph type and thermal optima. The existence of significant interactions between morph phenotype and location in the determination of thermal profile also indicated that while potentially important, any relationship between phenotype and thermal specialization is not consistent across populations.
This study demonstrates that black and striped morphs of P. cinereus are sometimes found in different thermal microhabitats. The two morphs seem to segregate into different thermal habitat in some locations (marginally significant for VPA in Figure 3). However, the difference is once again dependent upon the census location and season (Table 4). Differences between morphs in their thermal microhabitat distribution are consistent with data suggesting that morph frequencies can vary with air temperature  and with temperature proxies such as climate  and season . Furthermore, since we find that patterns of thermal segregation are not consistent, we confirm the contradictory patterns of niche preference in these previous studies. Lotter and Scott  found striped morphs in lower frequencies in regions defined as "warm", and Moreno  found that striped morph frequency decreased with increasing temperatures. In contrast, Test  found that striped morph frequency increased with increasing temperatures through the summer. We conclude from these previous studies that our results accurately reflect a relationship between thermal habitat segregation and colour morph that is inconsistent from population to population.
This study, as well as Moreno's study  found the black morph to have a lower metabolic rate under certain conditions. Either chromosomal or pleiotropic linkage between colour phenotype and metabolism could be responsible for this relationship. If pleiotropy is responsible for the relationship between these traits then there should be a consistent relationship across all populations. However, if chromosomal linkage is responsible, linkage could be inconsistent from population to population since linkage pairing could vary from population to population. Because we find differences in the relationship between MMR and morph type between locations, it appears that chromosomal linkage is more likely to be responsible for the association of colour and metabolism than pleiotropy. Correlated response to selection on colour morphs should, in this case, vary from location to location depending on the particular linkage within a given population or region.
Alternative explanations for the maintenance of the polymorphism given evidence (here and ) that the black morph has at times a lower maintenance expenditure remain unexplored. Possible explanations for maintenance of such colour polymorphisms are extensive and have been thoroughly reviewed by Roulin  and Hoffman and Blouin. Roulin has pointed out that proximate adaptive explanations for the existence of a colour polymorphism may fall into three classes: historical, direct, or indirect effects . A historical association between colour and other fitness related traits, due to divergent selection during allopatry, might be carried into the present. The apparent mosaic of various associations between colour morph and metabolic rate indicate that no single historical event is responsible for the polymorphism in P. cinereus.
There may be some direct selection on the colour polymorphism in P. cinereus. Differences in solar heating are unlikely to be of any importance in this species as is common in many species  since the species is not active during daylight. It is possible however that colour is important in defense or mating. Such direct selection on colour may support the polymorphism or it may balance indirect selective effects due to correlated response to selection on thermal niche. The potential for selection on the all red 'eurythristic' morph as a Batesian mimic of red spotted newts has been demonstrated [33, 34]. However the eurythristic form is too geographically restricted to account for the more common striped vs. black polymorphism. Striped morphs could be more cryptically coloured than the black morphs or, since predators often prey on the most common prey form, frequency-dependent selection can be maintaining the polymorphism even without mimicry . No studies, however, have yet documented direct selection supporting the common striped/black polymorphism.
This and previous studies [7, 10, 11] have found evidence for thermal niche differentiation of the two different colour morphs. This provides the opportunity for a correlated response of to selection on temperature specific metabolic rate. While this study fails to find evidence for unique thermal optima for the two morphs over typical summer field temperatures, seasonal dynamics may provide a more complex selective regime than examined here. The differences between populations in the way in which the colour morphs respond to temperature indicate that while correlated response to selection may be an important part of the dynamics of the polymorphism it is not a consistent one. We cannot eliminate the possiblity of plastic response of metabolic rate to environmental factors upon which morphs might be segregating. This type of influence has been documented in Tawny Owls  although in this case we have no evidence supporting such a plastic response. We speculate insead that some sort of indirect response to selection on thermal niche is probably important in maintenance of the polymorphism, but that historical linkage associations between the morphs and the genes under selection may be locally variable.
The two colour morphs of P. cinereus sometimes differ in their MMR profile and in the temperature microhabitat in which they are found. These differences are dependent on location and season. However these differences in MMR profile suggest that only one morph could enjoy a fitness advantage via this mechanism in these populations under the conditions studied. This is not consistent with the hypothesis that a simple correlated response to selection across multiple niches maintains the polymorphism. Instead, a correlated response to selection in conjunction with another selective regime(s) may maintain the polymorphism. Alternatively incomplete linkage between colour morph and the genes under selection may create a geographic mosaic of different correlated selective responses.
Two study sites, 6.7 kilometres apart were used. One site is located in O'Neil Woods Metro Park (O'Neil), Bath, Ohio (N41° 10.1120, W081° 35.448) and is 1.9 km2. The second site, the Valley Picnic Area on Riverview Road (VPA), is located in the Cuyahoga Valley National Park, Peninsula, Ohio (N41° 13.4410, W081° 33.7020) and is 6.3 km2. Both sites are hilly beech-maple forests, with well-drained soil, limited underbrush, and many cover objects.
Collection and maintenance of organisms
Animals (n = 101) from the two sites were collected for metabolic rate measurements. Salamanders were held for 5–6 days and housed individually in plastic containers with moistened paper towels. Salamanders were held at 10°C for the first 24 hours of captivity and then at the next day's test temperature for each subsequent night. Salamanders were kept on a light:dark cycle of 14:10 h for the summer trials and 12:12 h for the autumn trials. Each was fed six fruit flies (Drosophila melanogaster) every evening beginning two days before testing so that all animals were in a digestive state similar to free-ranging salamanders.
Metabolic rate measures
Temperature-dependent metabolic rate was measured at 10.0°C, 15.0°C, and 20.0°C. Metabolic measurements were conducted while salamanders were resting inside respirometry chambers under full-spectrum lighting, all housed within a digitally controlled incubator. Salamanders did not have access to radiant or other external sources of heat, insuring that their body temperatures were identical to the controlled temperature of the cabinet during all trials. The order of testing for the three temperatures was randomized for each individual. Measurement of MMR, the O2 consumption during the scotophase (07:00 to 20:00 hours local time) of fed animals , allowed inclusion of relevant sources of variation such as nutritional status and variation in metabolic rate not included in SMR [37, 38].
Automated closed system respirometry (Sable TR-3; Sable Systems International, Henderson, NV) was used to simultaneously record oxygen (Ametek S-3A, Pittsburgh, PA) and carbon dioxide (LiCor 6251, Lincoln, NE) levels. Five animals measured at each temperature over five evenly spaced automated cycles of 6 minutes each during any given day. The average O2 consumption rate (ml/h) of each animal was calculated using equations derived from Withers . During the trials, the animals sat on moistened sponges to prevent dehydration in 5 cc cylindrical chambers.
Field distributions across thermal environments
To determine correlations between temperatures and colour morph location, censuses were taken within each site. Sites were searched between June 12 and October 23, 2004. Five-hour censuses were conducted biweekly at each site by turning over potential cover objects. Data collected during each census included date, location, time, substrate temperature, temperature 5 cm above the substrate, colour phenotype, snout-vent length, and recapture status. Both substrate and air temperatures were measured simultaneously using type K thermocouples (Fluke 54 II Thermometer; 1999 Fluke Corporation, Everett, WA). Individuals (n = 993) were marked via toe clipping to prevent resampling of individuals.
Metabolic rate data were natural log (ln) transformed to linearize the relationship between oxygen consumption and mass [28, 29]. Repeated-measures MANCOVA was used to test for significant influence of location, season, and phenotype on MMR. Log transformed mass was included as a covariate to appropriately adjust for variation in the size of animals . Sequential model reduction resulted in a final model that included ln mass, location, season, phenotype, and phenotype X location as sources of variation in the MMR. Inspection of model residuals gave no indication of violation of MANCOVA assumptions.
Endler JA: Natural Selection in the Wild. Monographs in Population Biology. Edited by: May RM. 1986, Princeton, Princeton University Press, 21:
Nielsen MG, Watt WB: Interference competition and sexual selection promote polymorphism in Colias (Lepidoptera, Pieridae). Functional Ecology. 2000, 14: 718-730. 10.1046/j.1365-2435.2000.00472.x.
Roulin A: The evolution, maintenance and adaptive function of genetic colour polymorphism in birds. Biological Reviews. 2004, 79: 1-34. 10.1017/S1464793104006487.
Hoekstra HE, Krenz JG, Nachman MW: Local adaptation in the rock pocket mouse (Chaetodipus intermedius): natural selection and phylogenetic history of populations. Heredity. 2005, 94: 217-228. 10.1038/sj.hdy.6800600.
Whitney KD, Stanton ML: Insect seed predators as novel agents of selection on fruit color. Ecology. 2004, 85: 2153-2160.
Wlasiuk G, Garza JC, Lessa EP: Genetic and geographic differentiation in the rio negro tuco-tuco (Ctenomys rionegrensis): Inferring the roles of migration and drift from mulitple genetic markers. Evolution. 2003, 57: 913-926. 10.1554/0014-3820(2003)057[0913:GAGDIT]2.0.CO;2.
Lotter F, Scott NJJ: Correlation between climate and distribution of the color morphs of the salamander Plethodon cinereus. Copeia. 1977, 1977: 681-690. 10.2307/1443166.
Gibbs JP: Genetic structure of reback salamander Plethodon Cinereus populations in continuous and fragmented forests. Biological Conservation. 1998, 86: 77-81. 10.1016/S0006-3207(97)00173-0.
Highton R: Geographic variation in genetic dominance of the color morphs of the red-backed salamander, Plethodon cinereus. Genetics. 1975, 80: 363-374.
Test FH: Seasonal differences in populations of the red-backed salamander in southeastern Michigan. Papers of the Michigan Academy of Science, Arts, and Letters. 1955, 40: 137-153.
Moreno G: Behavioral and physiological differentiation between the color morphs of the salamander, Plethodon cinereus. Journal of Herpetology. 1989, 23: 335-341. 10.2307/1564043.
Hoffman EA, Blouin MS: A review of colour and pattern polymorphisms in anurans. Biological Journal of the Linnean Society. 2000, 70: 633-665. 10.1006/bijl.1999.0421.
Levene H: Genetic equilibrium when more than one ecological niche is available. American Naturalist. 1953, 87: 331-333. 10.1086/281792.
Levins R: Evolution in changing environments: Some theoretical explorations. 1968, Princeton, NJ, Princeton University Press
Gillespie J: Polymorphism in patchy environments. American Naturalist. 1974, 108: 145-151. 10.1086/282894.
Ayala FJ, Campbell CA: Frequency-dependent selection. Annual Review of Ecology and Systematics. 1974, 5: 115-138. 10.1146/annurev.es.05.110174.000555.
Horth L: Predation and the persistence of melanic male mosquitofish (Gambusia holbrooki). Journal of Evolutionary Biology. 2004, 17: 672-679. 10.1111/j.1420-9101.2004.00710.x.
Mayr E, Stresemann E: Polymorphism in the chat genus Oenanthe (Aves). Evolution. 1950, 4: 291-300. 10.2307/2405596.
Clarke BC: The evolution of genetic diversity. Proceedings of the Royal Society of London, Series B, Biological Sciences. 1979, 205: 453-474.
Wright S: Evolution and the Genetics of Populations, Volume 4: Variability Within and Among Natural Populations. 1978, Chicago, University of Chicago Press, 4:
Mather K: Polymorphism as an outcome of disruptive selection. Evolution. 1955, 9: 52-61. 10.2307/2405357.
Highton R: The inheritance of the color phases of Plethodon cinereus. Copeia. 1959, 1959: 33-37. 10.2307/1440097.
Bagnara JT, Taylor JD: Differences in pigment-containing organelles between color forms of the red-backed salamander, Plethodon cinereus. Cell and Tissue Research. 1970, 106: 412-417.
Petranka JW: Salamanders of the United States and Canada. 1998, Washington, D.C., Smithsonian Institution Press
Congdon JD, Dunham AE, Tinkle DW: Energy budgets and life histories of reptiles. Biology of the Reptila. Edited by: Pough FH. 1982, New York, Academic Press, 13:
Spotila JR, Standora EA: Energy budgets of ectothermic vertebrates. American Zoologist. 1985, 25: 973-986.
Dunham AE, Grant BW, Overall KL: Interfaces between biophysical and physiological ecology and the population ecology of terrestrial vertebrate ectotherms. Physiological Zoology. 1989, 62: 335-355.
Beaupre SJ, Duvall D: Variation in oxygen consumption of the western diamondback rattlesnake (Crotalus atrox): Implications for sexual size dimorphism. Journal of Comparative Physiology B Biological Sciences. 1998, 168: 497-506. 10.1007/s003600050170.
Currens CR, Niewiarowski PH, Whiteman HH: Effects of temperature and time of day on the resting metabolic rates of paedomorphic and metamorphic mole salamanders, Ambystoma talpoideum. Copeia. 2002, 2002: 489-495. 10.1643/0045-8511(2002)002[0489:EOTATO]2.0.CO;2.
Packard GC, Boardman TJ: The use of percentages and size specific indices to normalize physiological data for variation in body size: wasted time, wasted effort?. Comparative Biochemistry and Physiology. 1999, 122A: 37–44.-
Gatten REJ, Miller K, Full RJ: Energetics at rest and during locomotion. Environmental Physiology of the Amphibians. Edited by: Burggren WW. 1992, Chicago, University of Chicago Press
Feder ME: Acclimation to constant and variable temperatures in plethodontid salamanders: I. Rates of oxygen consumption. Comparative Biochemistry and Physiology. 1985, 81: 673-682. 10.1016/0300-9629(85)91046-1.
Brodie EDJ, Brodie EDIII: Differential avoidance of mimetic salamanders by freee-ranging birds. Science. 1980, 208: 181-182.
Tilley SG, Lundrigan BL, Brower LP: Erythrism and mimicry in the salamander Plethodon cinereus. Herpetologica. 1982, 38: 409-417.
Roulin A, Bize P, Tzaud N, Bianchi M, Ravussin PA, Christe P: Oxygen consumption in offspring tawny owls Strix aluco is associated with colour morph of foster mother. Journal of Ornithology. 2005, 146: 390-394. 10.1007/s10336-005-0096-3.
Anderson RA, Karasov WH: Energetics of the lizard Cnemidophorus tigris and life history consequences of food-acquisition mode. Ecological Monographs. 1988, 58: 79-110. 10.2307/1942462.
Niewiarowski PJ, Waldschmidt SR: Variation in metabolic rates of a lizard: Use of SMR in ecological contexts. Functional Ecology. 1992, 6: 15-22. 10.2307/2389766.
Beaupre SJ: An ecological study of oxygen consumption in the mottled rock rattlesnake, Crotalus lepidus lepidus, and the black-tailed rattlesnake, Crotalus molossus molossus, from two populations. Physiological Zoology. 1993, 66: 437-454.
Withers PC: Measurement of VO2, VCO2, and evaporative water loss with a flow-thorugh mask. Journal of Applied Physiology. 1977, 42: 120-123.
We thank C. Dostal for fieldwork and S.C. Weeks for comments. The Cuyahoga Valley National Park and the Summit County Metro Parks provided for study permits. This research was supported by grants to EEP from the Department of Biology at the University of Akron.
EEP designed the study, executed the experiments, collaborated on analysis of the data and produced the first draft of the manuscript. FBGM gave advice in the design of the study, took a large role in the statistical analyses and took primary responsibility for the final manuscript. PHN provided the expertise and equipment necessary for metabolic rate measures and assisted in drafting the manuscript.