- Research article
- Open Access
Increasing melanism along a latitudinal gradient in a widespread amphibian: local adaptation, ontogenic or environmental plasticity?
© Alho et al; licensee BioMed Central Ltd. 2010
Received: 29 June 2010
Accepted: 21 October 2010
Published: 21 October 2010
The thermal benefits of melanism in ectothermic animals are widely recognized, but relatively little is known about population differentiation in the degree of melanism along thermal gradients, and the relative contributions of genetic vs. environmental components into the level of melanism expressed. We investigated variation in the degree of melanism in the common frog (Rana temporaria; an active heliotherm thermoregulator) by comparing the degree of melanism (i) among twelve populations spanning over 1500 km long latitudinal gradient across the Scandinavian Peninsula and (ii) between two populations from latitudinal extremes subjected to larval temperature treatments in a common garden experiment.
We found that the degree of melanism increased steeply in the wild as a function of latitude. Comparison of the degree of population differentiation in melanism (P ST ) and neutral marker loci (F ST ) revealed that the P ST > F ST , indicating that the differences cannot be explained by random genetic drift alone. However, the latitudinal trend observed in the wild was not present in the common garden data, suggesting that the cline in nature is not attributable to direct genetic differences.
As straightforward local adaptation can be ruled out, the observed trend is likely to result from environment-driven phenotypic plasticity or ontogenetic plasticity coupled with population differences in age structure. In general, our results provide an example how phenotypic plasticity or even plain ontogeny can drive latitudinal clines and result in patterns perfectly matching the genetic differences expected under adaptive hypotheses.
Melanins and carotenoids are two classes of pigment compounds that are responsible for much of the variation in animal coloration . Melanins have been shown or suggested to be involved in a wide range of vital adaptive functions in animals, including signaling [2, 3], crypsis [4, 5], thermoregulation [6–10], protection from ultraviolet radiation [11, 12] and immune function . Whether melanism has a significant role in ectotherm thermoregulation has been subject to a lot of research over many decades , and the thermoregulatory explanation has been typically favored over other hypotheses. However, while the hypothesis about the thermally adaptive value of melanism has a lot of intuitive appeal because of its simplicity, the relationship between phenotype and genotype is complex [15–17], and other adaptive and non-adaptive causes should also be considered. Further, even if melanism in some circumstances is thermally adaptive, the question of whether it is a purely genetic or environmentally driven plastic trait is important, especially in light of the recent rapid anthropogenic changes in the environment . Namely, changes in degree of melanism in a population could be driven by direct environmental induction (i.e. plastic response) as well as by selection acting on heritable genetic variation (i.e. adaptation).
Intraspecific variation in melanism, when present, can express itself as distinct color morphs including completely dark, or melanistic, individuals [6, 18, 19], or as a continuum of animals with different amounts of dark pigmentation [20, 21]. The possible genetic basis of this variation has received considerable attention in literature. Several researchers have found significant family effects in common garden experiments, raising the possibility of genetic variation in the degree of melanization [14, 22]. There are also numerous molecular genetic studies that have investigated the association between melanism and sequence variation e.g. in the melanocortin-1 receptor gene [23–27] and its antagonist, the agouti signaling protein (agouti) [28, 29]. These molecular analyses have often, although not always, found a correspondence between certain genotypes and phenotypes. However, few studies have attempted to establish an adaptive basis for variation in melanization by linking the genetic variation to pigmentation clines observed in nature (but see e.g. [25, 30–33]).
Since melanism is genetically based in numerous taxa [5, 23, 24, 26] and has been shown to confer thermally adaptive advantage [6, 8–10], one would expect to observe intraspecific pigmentation clines along latitudinal or altitudinal gradients in the wild, as these gradients typically correspond to ambient temperature. However, such clines have been rarely reported in other taxa than insects  and the genetic basis, and thus the adaptive nature, of melanism clines remain little explored. To show that melanism can be thermally adaptive, one needs to simultaneously demonstrate both the existence and the genetic basis of a melanism cline along an environmental gradient. A failure to find genetic basis for an observed trend would indicate that the cline is caused either by environmentally induced plasticity (which can itself be either adaptive or non-adaptive), or by ontogenetic plasticity coupled with geographical differences in age or longevity. For instance, given that the degree of melanism in at least some anuran amphibians increases with age , and average age of individuals in a given population tends to increase with increasing latitude and altitude [34, 35], changes in population age structure could also be driving clines in the degree of dark pigmentation.
Here, we studied the existence and possible genetic basis of a latitudinal cline in the degree of melanization in the common frog (Rana temporaria). Amphibians are a particularly well-suited group of organisms to study the genetic basis of pigmentation due to the large color variation they exhibit within and between species [22, 36]. The common frog makes an especially interesting model since (i) it has one of the widest distribution ranges among amphibians including both high altitudes and latitudes , (ii) it is an active heliotherm thermoregulator , (iii) exhibits variation in the level of melanization [10, 27], and (iv) the degree of melanization in the species has been indicated to correlate positively with heating rate . Furthermore, many traits in this species have previously been shown to have undergone adaptive divergence along a latitudinal gradient [38–41].
The aims of our study were twofold. First, we investigated the presence of a latitudinal trend in the degree of melanization in wild-caught adult common frogs. Second, we used a common garden experiment to test whether the observed cline in melanism was genetically based, or reflected either plastic response to environmental heterogeneity or ontogenetic plasticity with geographical differences in age structure. To these ends, we compared the level of dorsal melanism in twelve populations along over 1500 km latitudinal gradient across the Scandinavian Peninsula, and conducted a common garden study applying larval thermal treatments using two latitudinally extreme populations. We hypothesized that a latitudinally ordered thermal adaptation might have occurred in Rana temporaria favouring elevated heating rates towards higher latitudes. Based on this hypothesis we predicted that the level of melanism would increase towards north and that the pattern would at least partially be genetically based and independent of rearing environment.
Data from the wild
The locations of the study populations and sample sizes (N) in the analyses
55° 30' N, 13° 25' E
55° 42' N, 13° 26' E
56° 19' N, 14° 42' E
59° 28' N, 13° 31' E
59° 51' N, 17° 14' E
60° 07' N, 16° 56' E
61° 16' N, 17° 11' E
62° 37' N, 17° 59' E
63° 49' N, 20° 14' E
65° 54' N, 16° 18' E
67° 51' N, 21° 02' E
69° 03' N, 20° 47' E
Common garden data
55° 42' N, 13° 26' E
68° 24' N, 18° 40' E
Common garden data
A common garden experiment was conducted for two populations along the latitudinal gradient across the Scandinavian Peninsula (Figure 2, Table 1). In the case of the southern population adult frogs were collected from spawning sites and kept in 4 °C until used in artificial fertilizations. Crosses were performed following North Carolina type I design . Briefly, the experiment utilized twenty females and ten males from the southern population, and each male was crossed with two females. This resulted in twenty full-sib and ten paternal half-sib families. Artificial fertilization was performed following  with some modifications. Males were first injected with hormone, the cloaca was rinsed with Ringer's solution and the sperm solution was separated into two containers. The eggs were squeezed out from females' cloaca to the sperm solution with the help of blunt forceps. In the case of the northern population we collected eggs from the wild. These naturally laid eggs were subsequently handled in similar manner to those from artificial crosses after fertilization.
The eggs and hatchlings from each full-sib family were raised in two 3 L vials (ca. 200 eggs in each) until Gosner stage 25 . Temperature in the laboratory room was 18 °C and the light rhythm 18L:6D. After reaching stage 25, ten tadpoles from each cross were raised individually in 1 L vials at two temperature treatments, 18 °C (high) and 14 °C (low), in two separate climate-controlled rooms with 18L:6D light rhythm. There were thus five tadpoles per cross per treatment. During this time the tadpoles were fed slightly boiled spinach ad libitum. To assure homogenous water quality, the animals were reared throughout the aquatic stage in reconstituted soft water (RSW)  renewed every third day. After metamorphosis, the tadpoles were moved to 1 L boxes with the bottom covered with moss and with a dripping system maintaining moisture in a climate room maintained at 16 °C. At this stage, the juvenile frogs were fed with unlimited amount of Drosophila flies and small crickets. In late November, the juveniles were moved to another climate chamber where the conditions were gradually changed to resemble wintering conditions (0L:24D, 4 °C). They were allowed to overwinter until late February, when the conditions were gradually changed to correspond to the summer conditions again. The juveniles were then maintained under similar conditions as above. Nine months after metamorphosis, all juveniles were photographed for the estimation of the degree of melanism (note that this date varies between individuals due to the difference in phenology between the populations and in development rate between the larval temperature treatments). Snout-vent length of the juveniles was measured from the digital photographs.
Mortality limited the available data in the southernmost population mostly to full-sib families. Since also the data on the northernmost population came from full-sib families, we restricted our analysis below to family effects, and did not attempt to estimate sire and dam effects.
Estimation of level of melanism
Estimation of the relative area covered by the permanent melanistic spots was done similarly with the wild-caught and common garden samples. The individuals' (thawed carcasses in the wild-caught and live specimens in the common garden sample) dorsal sides were gently pressed against a transparent plexiglass panel to obtain a nearly flat surface for photographing. Digital photographs were taken under similar circumstances (different sets for wild-caught and common garden samples). A ruler or millimeter paper was placed in every photograph for scaling. Color calibration was not employed, but melanistic and non-melanistic areas were clearly distinguishable in the images, as shown by the relatively high repeatability estimate for the degree of melanism (see Results).
The digital photographs were processed with Image Pro Plus 4.5 software (Media Cybernetics Inc., Bethesda, USA). We first defined the area (the dorsal torso without head and appendages; Figure 1b) and then the color of the spots of interest. Then, with the aid of the scale on each photograph, the software automatically calculated the area of each spot and the remaining area. The degree of melanism was defined as the proportion of dorsal surface area covered by pigmentation spots larger than 1 mm2. 24 images were measured twice in order to estimate the repeatability for length and the degree of melanism.
Since the degree of melanism was a proportion, it was arcsine transformed for the statistical analyses, i.e. we calculated the arcsine of the square root of the proportion. Back-transformation, when needed, was done by taking the square of the sine of the value.
where MS A was the among individuals mean squares, MS W the within individuals means squares, and n the number of measurements per individual, i.e. two. 95% confidence intervals for the repeatability estimates were obtained by non-parametric bootstrap, resampling the data 5000 times.
where was the variance between populations, the variance within populations, and h2 the heritability. P ST was estimated for two scenarios, with heritability values h 2 = 1 and h 2 = 0.5. The variance components were estimated using the lmer function. We fitted a linear mixed model with sex and snout-vent length as fixed effects and population as a random effect. 95% confidence intervals for the P ST estimates were obtained by non-parametric bootstrap, resampling data 5000 times.
An estimate of F ST - describing the degree of neutral genetic divergence [51, 52] - calculated based on six of our twelve study populations was available from a previous study . We calculated P ST separately both for the six populations and for all twelve populations.
where was the family variance and residual variance. Family and residual variances were assumed to be equal between the two populations. 95% confidence interval for the heritability estimate was obtained by non-parametric bootstrap, resampling family data 5000 times.
The repeatability for the juveniles reared in common garden was 1.00 (F23,24 = 427.26; 95% CI: 0.98-1.00) for snout-vent length and 0.81 (F23,24 = 9.35; 95% CI: 0.25-1.00) for the degree of melanism. The snout-vent length was 16.9-27.1 mm in the common garden juveniles and 53.3-90.6 mm in the adults from the wild.
Fixed effect estimates for the degree of melanism in wild-caught adult common frogs
Latitude × Sex
Snout-vent length × Sex
Fixed effect estimates for degree of melanism in juvenile common frogs reared in common garden
Population × Temperature
The thermoregulation hypothesis, stating that dark individuals have a fitness advantage in low temperatures as compared to light individuals due to their ability to heat up faster , is perhaps the most prominent adaptive explanation for melanism in ectotherms. However, few studies have tested whether the pigmentation clines observed in nature and explained as thermal adaptations are genetically based as assumed by the adaptive hypothesis (but see e.g. [30, 32, 33]). Here we detected an increasing latitudinal trend across the Scandinavian Peninsula in the degree of melanism in common frogs, but found no evidence of corresponding genetic differentiation in a common garden experiment with different larval temperature treatments. In fact, the observed genetic differences between populations were in the opposite direction of the latitudinal trend observed in the wild. As the evidence thus points against direct local adaptation as an explanation for the patterns we found in the wild, it seems plausible that environmental plasticity and/or ontogenetic plasticity together with geographic differences in age structure are responsible for the observed latitudinal trend. Interestingly, the trend persisted even after correcting for differences in snout-vent length - a rough proxy for age  - suggesting that the cline in nature could be more than a simple reflection of differences in age structure. This interpretation was supported by a significant environmental effect in the common garden experiment, with the juveniles originating from the northern population exhibiting higher degree of melanism in the colder larval rearing environment (Table 3). While previous studies on amphibians have found that both larval and adult coloration is affected by temperature and developmental stage [54–56], this is, to our knowledge, the first to suggest that the effects of larval environment on coloration can be carried over to the terrestrial stage. In light of these results, it is not necessarily surprising that an earlier study  did not find an association between MC1R sequence variation and variation in the degree of melanism in our study species.
Comparison of quantitative and neutral marker differentiation, as measured by Q ST and F ST , respectively, is often used to infer the relative roles of natural selection and random drift in among-population divergence [51, 52]. The argument is that when Q ST exceeds F ST there is evidence of divergent selection, when the reverse is true there is support for stabilizing selection, and when Q ST and F ST are not significantly different, random genetic drift as the only cause of phenotypic divergence cannot be excluded . When Q ST estimates have not been available, P ST values calculated from purely phenotypic data have been used as surrogates [50, 57]. In our study, the degree of phenotypic divergence, measured as P ST , was significantly higher than the F ST estimate published in  for our populations, when the more realistic P ST value estimated assuming h 2 = 0.5 was used. Hence, looking naïvely at the phenotypic trend and divergence, without a further genetic analysis or study of population demography, the results suggest that the populations might have diverged under natural selection consistent with thermal adaptation hypothesis of melanism. However, when looking at the genetic population effects, not only do we find that there is no direct genetic basis for the latitudinal trend, but also that the observed population effects are to the opposite direction of the cline. Our results are thus in line with the criticism of the use of P ST in studies of wild populations .
The problem of disentangling genetic and environmental effects in studies of geographic or temporal phenotypic variation is a general one [16, 59, 60]. The relevance of this issue is likely to increase with ongoing rapid anthropogenic changes in the environment and the increasing integration of evolutionary biology to nature conservation in predicting the long-term effects of these changes . Our study cautions once again against interpreting phenotypic trends - even if making perfect evolutionary sense - as evidence of local adaptation without explicit genetic analyses. In other words, environmental plasticity or even plain ontogeny can produce patterns that conform precisely to those expected under adaptive hypotheses (see also e.g. [59, 61, 62]).
Environmental plasticity itself can be either adaptive or non-adaptive . Although somewhat speculative with only two populations, it is interesting to note that while larval rearing temperature had a significant effect on the degree of melanism in the common garden juveniles in the northern population (Table 3), there was also a significant interaction between population and temperature, with the southern population predicted to reverse the pattern observed in the north. This difference suggests that if melanism serves a thermoregulatory function in the north, natural selection might have acted on environmental plasticity, causing the individuals in the north to respond to cold temperature with increased melanistic coloration. However, the limited number of populations and the lack of accurate age data from the wild prevent any firm conclusions. The observed pattern in the wild might also reflect plain age differences or adaptive environmental plasticity serving other functions, such as protection from ultraviolet radiation , crypsis [4, 5], disease resistance , or multiple simultaneous functions . Finally, the pattern might arise from a correlation with an adaptive trait, or from entirely non-adaptive environmental plasticity . Without measuring age and selection, the adaptive significance of the degree of melanism in the common frog is, however, speculative and the realm of future studies.
Our study comes with a number of caveats. The sample size in the common garden data was small in terms of populations, families and offspring. This was reflected e.g. in the fact that the confidence interval for the upper limit for heritability was very wide and prohibited any firm conclusions on heritability. However, the small sample size is taken into account in the confidence intervals and highest posterior density intervals, and hence the results are robust e.g. as to finding that in the common garden there was no population difference similar to that observed in the wild. The low number of offspring in the common garden analysis resulted partly from low survival in the southern population, but although mortality was thus high, there is no reason to suspect that it would have influenced the results. It is also true that the age structure of the wild and common garden data did not overlap, a fact reflected in the difference between the ranges of snout-vent lengths (see Results). Finally, while environmental conditions can be controlled in the laboratory, this also often unavoidably distorts them from natural conditions. For example, while we had two larval temperature treatments, neither of them arguably fully corresponded to the normal temperature conditions common frogs experience in the wild - while temperatures in the experiment were constant, in the wild they fluctuate daily and seasonally and vary spatially. In addition, the temperature treatments were limited to larval stage, ending after metamorphosis. Although these issues should be kept in mind when interpreting our results, the conclusions we have drawn should be justified in the light of the available evidence and information.
Latitudinal and altitudinal clines in the degree of melanism provide an opportunity to study both the evolution of coloration and evolutionary processes in general. A genetic basis is necessary for any adaptive explanation for such a cline to be applicable. While we observed a positive correlation between latitude and the degree of melanistic pigmentation in the common frog, there was no straightforward genetic foundation for it. This emphasizes the potential role of environment in the degree of melanism, while leaving population demographic explanations and the possible genetic basis of environmental plasticity to be addressed in future studies.
We thank John Loehr and three anonymous reviewers for helpful comments, Maarit Haataja, Ane Laugen, Katja Räsänen and Mattias Sterner for helping with the fieldwork, Jonas Andersson for assisting in the common garden experiment, and Heidi Kristanjankroon for measuring the degree of melanism in the frogs reared in the common garden. Our research was supported by the Swedish Research Council (KIJ, AL, JM), Academy of Finland (GH, AL, JM), Finnish Ministry of Education (JM) and LUOVA Graduate School (JSA).
- Fox DL: Animal biochromes and structural colours: physical, chemical, distributional & physiological features of coloured bodies in the animal world. 1976, Berkeley, California: University of California Press, 2Google Scholar
- Jawor JM, Breitwisch R: Melanin ornaments, honesty, and sexual selection. Auk. 2003, 120: 249-265. 10.1642/0004-8038(2003)120[0249:MOHASS]2.0.CO;2.View ArticleGoogle Scholar
- Griffith SC, Parker TH, Olson VA: Melanin-versus carotenoid-based sexual signals: is the difference really so black and red?. Anim Behav. 2006, 71: 749-763. 10.1016/j.anbehav.2005.07.016.View ArticleGoogle Scholar
- Hoekstra HE, Nachman MW: Different genes underlie adaptive melanism in different populations of rock pocket mice. Mol Ecol. 2003, 12: 1185-1194. 10.1046/j.1365-294X.2003.01788.x.View ArticleGoogle Scholar
- Hoekstra HE, Hirschmann RJ, Bundey RA, Insel PA, Crossland JP: A single amino acid mutation contributes to adaptive beach mouse color pattern. Science. 2006, 313: 101-104. 10.1126/science.1126121.View ArticleGoogle Scholar
- Gibson AR, Falls JB: Thermal biology of the common garter snake Thamnophis sirtalis (L.) II. The effects of melanism. Oecologia. 1979, 43: 99-109. 10.1007/BF00346675.View ArticleGoogle Scholar
- de Jong PW, Gussekloo SWS, Brakefield PM: Differences in thermal balance, body temperature and activity between non-melanic and melanic two-spot ladybird beetles (Adalia bipunctata) under controlled conditions. J Exp Biol. 1996, 199: 2655-2666.Google Scholar
- Forsman A: Heating rates and body temperature variation in melanistic and zigzag Vipera berus: does colour make a difference?. Ann Zool Fenn. 1995, 32: 365-374.Google Scholar
- Forsman A: Thermal capacity of different colour morphs in the pygmy grasshopper Tetrix subulata. Ann Zool Fenn. 1997, 34: 145-149.Google Scholar
- Vences M, Galán P, Vieites DR, Puente M, Oetter K, Wanke S: Field body temperatures and heating rates in a montane frog population: the importance of black dorsal pattern for thermoregulation. Ann Zool Fenn. 2002, 39: 209-220.Google Scholar
- Kollias N, Sayre RM, Zeise L, Chedekel MR: Photoprotection by melanin. J Photochem Photobiol B. 1991, 135-160. 10.1016/1011-1344(91)80147-A.Google Scholar
- Hofer R, Mokri C: Photo-protection in tadpoles of the common frog, Rana temporaria. J Photochem Photobiol B. 2001, 59: 48-53. 10.1016/S1011-1344(00)00136-6.View ArticleGoogle Scholar
- Wilson K, Cotter SC, Reeson AF, Pell JK: Melanism and disease resistance in insects. Ecol Lett. 2001, 4: 637-649. 10.1046/j.1461-0248.2001.00279.x.View ArticleGoogle Scholar
- Clusella Trullas S, van Wyk JH, Spotila JR: Thermal melanism in ectotherms. J Therm Biol. 2007, 32: 235-245. 10.1016/j.jtherbio.2007.01.013.View ArticleGoogle Scholar
- Merilä J, Kruuk LEB, Sheldon BC: Cryptic evolution in a wild bird population. Nature. 2001, 412: 76-79. 10.1038/35083580.View ArticleGoogle Scholar
- Gienapp P, Teplitsky C, Alho JS, Mills JA, Merilä J: Climate change and evolution: disentangling environmental and genetic responses. Mol Ecol. 2008, 17: 167-178. 10.1111/j.1365-294X.2007.03413.x.View ArticleGoogle Scholar
- Ozgul A, Tuljapurkar S, Benton TG, Pemberton JM, Clutton-Brock TH, Coulson T: The dynamics of phenotypic change and the shrinking sheep of St. Kilda. Science. 2009, 325: 464-467. 10.1126/science.1173668.View ArticleGoogle Scholar
- King RB: Color pattern polymorphism in the Lake Erie water snake, Nerodia sipedon insularium. Evolution. 1987, 41: 241-255. 10.2307/2409135.View ArticleGoogle Scholar
- Andrén C, Nilson G: Reproductive success and risk of predation in normal and melanistic colour morphs of the adder, Vipera berus. Biol J Linn Soc Lond. 1981, 15: 235-246. 10.1111/j.1095-8312.1981.tb00761.x.View ArticleGoogle Scholar
- Riobó A, Rey J, Puente M, Miramontes C, Vences M: Ontogenetic increase of black dorsal pattern in Rana temporaria. British Herpetological Society Bulletin. 2000, 70: 1-6.Google Scholar
- Cheviron ZA, Hackett SJ, Brumfield RT: Sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R) is not associated with plumage variation in the blue-crowned manakin (Lepidopthrix coronata). Proc R Soc Lond B. 2006, 273: 1613-1618. 10.1098/rspb.2006.3499.View ArticleGoogle Scholar
- Hoffman EA, Blouin MS: A review of colour and pattern polymorphisms in anurans. Biol J Linn Soc Lond. 2000, 70: 633-665. 10.1111/j.1095-8312.2000.tb00221.x.View ArticleGoogle Scholar
- Ritland K, Newton C, Marshall HD: Inheritance and population structure of the white-phased "Kermode" bear. Curr Biol. 2001, 11: 1468-1472. 10.1016/S0960-9822(01)00448-1.View ArticleGoogle Scholar
- Theron E, Hawkins K, Bermingham E, Ricklefs RE, Mundy NI: The molecular basis of an avian plumage polymorphism in the wild: A melanocortin-1 receptor point mutation is perfectly associated with the melanic plumage morph of the bananaquit, Coereba flaveola. Curr Biol. 2001, 11: 550-557. 10.1016/S0960-9822(01)00158-0.View ArticleGoogle Scholar
- Mundy NI, Badcock NS, Hart T, Scribner K, Janssen K, Nadeau NJ: Conserved genetic basis of a quantitative plumage trait involved in mate choice. Science. 2004, 303: 1870-1873. 10.1126/science.1093834.View ArticleGoogle Scholar
- Rosenblum EB, Hoekstra HE, Nachman MW: Adaptive reptile colour variation and the evolution of the MC1R gene. Evolution. 2004, 58: 1794-1808.Google Scholar
- Herczeg G, Matsuba C, Merilä J: Sequence variation in the melanocortin-1 receptor gene (MC1R) does not explain variation in the degree of melanism in a widespread amphibian. Ann Zool Fenn. 2010, 47: 37-45.View ArticleGoogle Scholar
- Eizirik E, Yuhki N, Johnson WE, Menotti-Raymond M, Hannah SS, O' Brien SJ: Molecular genetics and evolution of melanism in the cat family. Curr Biol. 2003, 13: 448-453. 10.1016/S0960-9822(03)00128-3.View ArticleGoogle Scholar
- Kingsley EP, Manceau M, Wiley CD, Hoekstra HE: Melanism in Peromyscus is caused by independent mutations in Agouti. PLoS One. 2009, 4: e6435-10.1371/journal.pone.0006435.PubMed CentralView ArticleGoogle Scholar
- Dearn JM: Latitudinal cline in a colour pattern polymorphism in the Australian grasshopper Phaulacridium vittatum. Heredity. 1981, 47: 111-119. 10.1038/hdy.1981.64.View ArticleGoogle Scholar
- Loehr J, Worley K, Moe J, Carey J, Coltman DW: MC1R variants correlate with thinhorn sheep colour cline but not individual colour. Can J Zool. 2008, 86: 147-150. 10.1139/Z07-130.View ArticleGoogle Scholar
- Rajpurohit S, Parkash R, Ramniwas S, Singh S: Variations in body melanisation, ovariole number and fecundity in highland and lowland populations of Drosophila melanogaster from the Indian subcontinent. Insect Sci. 2008, 15: 553-561. 10.1111/j.1744-7917.2008.00245.x.View ArticleGoogle Scholar
- Parkash R, Sharma V, Kalra B: Correlated changes in thermotolerance traits and body color phenotypes in montane populations of Drosophila melanogaster: analysis of within- and between-population variations. J Zool. 2010, 280: 49-59. 10.1111/j.1469-7998.2009.00641.x.View ArticleGoogle Scholar
- Miaud C, Guyétant R, Elmberg J: Variations in life-history traits in the common frog Rana temporaria (Amphibia: Anura): a literature review and new data from the French Alps. J Zool. 1999, 249: 61-73. 10.1111/j.1469-7998.1999.tb01060.x.View ArticleGoogle Scholar
- Laugen AT, Laurila A, Jönsson KI, Söderman F, Merilä J: Do common frogs (Rana temporaria) follow Bergmann's rule?. Evol Ecol Res. 2005, 7: 717-731.Google Scholar
- Richards CM, Nace GW: Dark pigment variants in anurans: classification, new descriptions, color changes and inheritance. Copeia. 1983, 979-990. 10.2307/1445100.Google Scholar
- Gasc J-P, Cabela A, Crnobrnja-Isailovic J, Dolmen D, Grossenbacher K, Haffner P, Lescure J, Martens H, Martínez Rica JP, Maurin H, Oliveira ME, Sofianidou TS, Veith M, Zuiderwijk A: Atlas of amphibians and reptiles in Europe. 1997, Paris: Societas Europea Herpetologica & Museum National d'Histoire Naturelle (IEGB/SPN)Google Scholar
- Laugen AT, Laurila A, Räsänen K, Merilä J: Latitudinal countergradient variation in the common frog (Rana temporaria) development rates - evidence for local adaptation. J Evol Biol. 2003, 16: 996-1005. 10.1046/j.1420-9101.2003.00560.x.View ArticleGoogle Scholar
- Palo JU, O'Hara RB, Laugen AT, Laurila A, Primmer CR, Merilä J: Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data. Mol Ecol. 2003, 12: 1963-1978. 10.1046/j.1365-294X.2003.01865.x.View ArticleGoogle Scholar
- Hettyey A, Laurila A, Herczeg G, Jönsson KI, Kovács T, Merilä J: Does testis weight decline towards the Subarctic? A case study on the common frog, Rana temporaria. Naturwissenschaft. 2005, 92: 188-192. 10.1007/s00114-005-0607-3.View ArticleGoogle Scholar
- Jönsson KI, Herczeg G, O'Hara RB, Söderman F, ter Schure AFH, Larsson P, Merilä J: Sexual patterns of prebreeding energy reserves in the common frog Rana temporaria along a latitudinal gradient. Ecography. 2009, 32: 831-839. 10.1111/j.1600-0587.2009.05352.x.View ArticleGoogle Scholar
- Fog K, Schmedes A, Rosenørn de Lasson D: Nordens padder og krybdyr. 1997, Copenhagen: Gads ForlagGoogle Scholar
- Laugen AT, Laurila A, Merilä J: Maternal and genetic contributions to geographical variation in Rana temporaria larval life-history traits. Biol J Linn Soc Lond. 2002, 76: 61-70. 10.1111/j.1095-8312.2002.tb01714.x.View ArticleGoogle Scholar
- Lynch M, Walch B: Genetics and Analysis of Quantitative Traits. 1998, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Berger L, Rybacki M, Hotz H: Artificial fertilization of water frogs. Amphibia-Reptilia. 1994, 15: 108-113.Google Scholar
- Gosner KL: A simplified table for staging anuran embryos and larvae with notes on identification. Copeia. 1960, 183-190.Google Scholar
- American Public Health Association (APHA): Standard methods for the examination of water and wastewater. 1985, Washington, 16Google Scholar
- Lessells CM, Boag PT: Unrepeatable repeatabilities - a common mistake. Auk. 1987, 104: 116-121.View ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2009, Vienna, Austria: R Foundation for Statistical Computing, [http://www.r-project.org]Google Scholar
- Leinonen T, Cano JM, Makinen H, Merilä J: Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks. J Evol Biol. 2006, 19: 1803-1812. 10.1111/j.1420-9101.2006.01182.x.View ArticleGoogle Scholar
- Merilä J, Crnokrak P: Comparison of genetic differentiation at marker loci and quantitative traits. J Evol Biol. 2001, 14: 892-903. 10.1046/j.1420-9101.2001.00348.x.View ArticleGoogle Scholar
- Leinonen T, O'Hara RB, Cano JM, Merilä J: Comparative studies of quantitative trait and neutral marker divergence: a meta-analysis. J Evol Biol. 2008, 21: 1-17.Google Scholar
- Guarino FM, Di Già I, Sindaco R: Age structure in a declining population of Rana temporaria from Northern Italy. Acta Zool Academ Sci Hung. 2008, 54: 99-112.Google Scholar
- Duellman WE, Trueb L: Biology of amphibians. 1986, Baltimore, Maryland: Johns Hopkins University PressGoogle Scholar
- Garcia TS, Straus R, Sih A: Temperature and ontogenetic effects on color change in the larval salamander species Ambystoma barbouri and Ambystoma texanum. Can J Zool. 2003, 81: 710-715. 10.1139/z03-036.View ArticleGoogle Scholar
- Stegen JM, Ginger CM, Sun L: The control of color change in the Pacific tree frog, Hyla regilla. Can J Zool. 2004, 82: 889-896. 10.1139/z04-068.View ArticleGoogle Scholar
- Merilä J: Quantitative trait and allozyme divergence in the Greenfinch (Carduelis chloris, Aves: Fringillidae). Biol J Linn Soc Lond. 1997, 61: 243-266.View ArticleGoogle Scholar
- Pujol B, Wilson AJ, Ross RIC, Pannell JR: Are QST-FST comparisons for natural populations meaningful?. Mol Ecol. 2008, 17: 4782-4785. 10.1111/j.1365-294X.2008.03958.x.View ArticleGoogle Scholar
- Teplitsky C, Mills JA, Alho JS, Yarrall JW, Merilä J: Bergmann's rule and climate change revisited: Disentangling environmental and genetic responses in a wild bird population. Proc Natl Acad Sci USA. 2008, 105: 13492-13496. 10.1073/pnas.0800999105.PubMed CentralView ArticleGoogle Scholar
- Merilä J: It's the genotype, stupid!. J Anim Breed Genet. 2009, 126: 1-2. 10.1111/j.1439-0388.2008.00795.x.View ArticleGoogle Scholar
- Przybylo R, Sheldon BC, Merilä J: Climatic effects on breeding and morphology: evidence of phenotypic plasticity. J Anim Ecol. 2000, 69: 395-403. 10.1046/j.1365-2656.2000.00401.x.View ArticleGoogle Scholar
- Réale D, McAdam AG, Boutin S, Berteaux D: Genetic and plastic responses of a northern mammal to climate change. Proc R Soc Lond B. 2003, 270: 591-596. 10.1098/rspb.2002.2224.View ArticleGoogle Scholar
- Ghalambor CK, McKay JK, Carroll SP, Reznick DN: Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007, 21: 394-407. 10.1111/j.1365-2435.2007.01283.x.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.