Studying the geographic variation of phenotypic traits can provide key information about the potential adaptive function of alternative phenotypes. Gloger's rule posits that animals should be dark- vs. light-colored in warm and humid vs. cold and dry habitats, respectively. The rule is based on the assumption that melanin pigments and/or dark coloration confer selective advantages in warm and humid regions. This rule may not apply, however, if genes for color are acting on other traits conferring fitness benefits in specific climes. Covariation between coloration and climate will therefore depend on the relative importance of coloration or melanin pigments and the genetically correlated physiological and behavioral processes that enable an animal to deal with climatic factors. The Barn Owl (Tyto alba) displays three melanin-based plumage traits, and we tested whether geographic variation in these traits at the scale of the North American continent supported Gloger's rule. An analysis of variation of pheomelanin-based reddish coloration and of the number and size of black feather spots in 1,369 museum skin specimens showed that geographic variation was correlated with ambient temperature and precipitation. Owls were darker red in color and displayed larger but fewer black feather spots in colder regions. Owls also exhibited more and larger black spots in regions where the climate was dry in winter. We propose that the associations between pigmentation and ambient temperature are of opposite sign for reddish coloration and spot size vs. the number of spots because selection exerted by climate (or a correlated variable) is plumage trait–specific or because plumage traits are genetically correlated with different adaptations.
INTRODUCTION
The study of geographic variation in phenotypic traits is key to exploring the adaptive function of alternative phenotypes. At a large geographic scale, such as across a continent, ecological and climatic factors vary greatly among regions, promoting the evolution of specific phenotypes adapted to these specific local conditions (Huey et al. 2000). Local adaptation requires populations to differ genetically (Blanquart et al. 2013), potentially leading to speciation even in the presence of gene flow (Barton 2010). Although studying geographic variation in phenotypic traits has a long history (Wright 1943, Mayr 1956), this approach still remains essential to identifying the potential adaptive function of alternative phenotypes in changing environments (Newton 2003).
As early as 1833, Constantin Gloger noticed that humidity and temperature accounted for geographic variation in bird coloration (Gloger 1833, Caro 2005, Kamilar and Bradley 2011). This led to the development of “Gloger's rule,” which states that animals should be darkly colored in warm and humid environments and paler in cold and dry habitats. This rule applies to a large number of animals, including mammals and insects (Mayr 1963, Stoner et al. 2003, Kamilar and Bradley 2011), and in 50 out of 54 (96%) studies of North American birds, populations were more pigmented in more humid areas (Zink and Remsen 1986). A number of selective processes can explain this pattern (Burtt and Ichida 2004). First, predator–prey interactions may be at the heart of Gloger's rule if dark and pale individuals are cryptic in different habitats (Miller and Miller 1951, Burtt 1981, Zink and Remsen 1986). Second, animals may be darker in wetter habitats because melanin pigments play a role in thermoregulation, or because heavily pigmented feathers repel water better, allowing the body to dry more rapidly (Heppner 1970, Walsberg et al. 1978, Burtt 1981, Clusella-Trullas et al. 2008). Third, melanin pigments may protect the integument against ectoparasites (Ruiz-De-Castañeda et al. 2012) and bacteria (Burtt and Ichida 2004, Goldstein et al. 2004) that thrive in warm and humid habitats (Guernier et al. 2004). Finally, melanin pigments improve hardness of the feathers and hence protect them against abrasion (Bonser and Witter 1993), which might be more pronounced in warm and humid areas.
Recently, it has been proposed that genes that participate in the production of melanin pigments may pleiotropically regulate other phenotypes (Ducrest et al. 2008). These pleiotropic genes may therefore account for the covariation between melanin-based coloration, behaviour, and physiology (Roulin and Ducrest 2011), regardless of whether melanin-based coloration is involved in camouflage or whether melanin pigments participate in thermoregulation or in the protection of the integument against parasites and abrasion. Therefore, not only selection exerted on coloration but also selection exerted on genetically correlated physiological or behavioral traits may account for geographic variation. As a consequence, in contrast to Gloger's rule, we may find that individuals are heavily pigmented in colder and dryer habitats if genes encoding for the production of melanin pleiotropically regulate physiological and/or behavioral processes adapted to such habitats. This mechanism may explain why dark, melanic Pied Flycatchers (Ficedula hypoleuca; Sirkiä et al. 2010) and Tawny Owls (Strix aluco; Karell et al. 2011) are more sensitive to the cold than their pale conspecifics.
Here, we studied geographic variation in three melanin-based plumage traits in North American Barn Owls (Tyto alba pratincola). On museum skin specimens, we counted and measured the size of black eumelanic feather spots, as well as measuring pheomelanin-based coloration, which varies from white to dark red. In European Barn Owls (T. a. alba, T. a. guttata), these three plumage traits are genetically correlated with behavioral and physiological processes (Roulin and Ducrest 2011). For instance, smaller-spotted and redder pheomelanic owls have a higher appetite (Dreiss et al. 2010), and larger-spotted owls are more resistant to parasitism and to various sources of stress that induce a rise in glucocorticoids (Roulin et al. 2000, Almasi et al. 2012). Because the genetics of melanogenesis are very well conserved across vertebrates (Schiöth et al. 2005), plumage traits may be associated with similar phenotypes in North American as in European Barn owls. In this paper, we first describe geographic variation in the three plumage traits and then examine whether rainfall, as a proxy for humidity, and temperature account for geographic variation as predicted by Gloger's rule. According to this rule, we would expect Barn Owls to be darker red and to display more and larger black spots in warm and humid compared with cold and dry regions. However, if plumage traits are genetically correlated with physiological and behavioral traits that are useful for coping with variation in climate or in any correlated environmental variable, plumage traits may evolve through selection exerted by climate on these genetically correlated traits rather than on the plumage traits themselves. Depending on how melanin-based plumage traits are genetically correlated with these physiological and behavioral adaptations, geographic variation in melanin-based traits may or may not follow Gloger's rule.
METHODS
We obtained data on 1,369 barn owl skins (636 males and 733 females) collected between the years 1844 and 2013 (mean year is 1948) in North America (Appendix A) by 59 natural history museums (Appendix B). Plumage traits were assessed on the breast and belly, and mean values over these two body parts were used in the statistical analyses. These methods were reliable, as shown by measuring the same individuals twice (repeatability was 0.90; Roulin 1999, 2004a). A. Roulin compared the degree of reddish coloration against eight chips ranging from −8 for white to −1 for dark red (Roulin et al. 1998); we have previously shown that coloration does not change in museum specimens (Roulin 2003) and that coloration scored by eye is strongly correlated (Pearson's correlation: r = 0.78, n = 1,107, P < 0.001) with reflectance in the range 400–700 nm as measured with an Ocean Optics S2000 spectrometer (Ocean Optics, Dunedin, Florida, USA) and a PX-2 xenon lamp (Ocean Optics; Dreiss and Roulin 2010). Coloration of skin specimens was not measured with a spectrophotometer because our study started in 1996 when we did not have access to a spectrophotometer. Black spots were counted within a 60 × 40 mm frame and their diameter was measured to the nearest 0.1 mm using calipers. We scored reddish coloration in all 1,369 individuals, and we counted and measured the size of spots in 1,327 and 1,368 individuals, respectively.
The geographic coordinates of the sites where each skin specimen was collected allowed us to extract long-term monthly sums of precipitation and monthly daily mean 2-m ambient temperature from the Worldclim dataset for the period 1950–2000 ( http://www.worldclim.org/, Hijmans et al. 2005) at a 30-arc-second spatial resolution. For each individual we had one value per day, for each year we calculated a mean monthly value, and for each month we computed a mean value over the entire period 1950–2000. Therefore, for each skin specimen and each month, we had a mean precipitation level value and a mean ambient temperature value.
We log10-transformed monthly precipitation values to generate a normal distribution for the dataset. To reduce the number of weather variables and to obtain integrative and independent indices of monthly precipitation and ambient temperature, we performed a principal components analysis (PCA; Table 1). We extracted four components (eigenvalues >1.0) and used them to examine how melanin-based plumage traits were associated with temperature and precipitation.
TABLE 1.
Principal components analysis of monthly sums of precipitation (log-transformed) and of monthly means of daily mean ambient temperatures in North America. Only the components with eigenvalues larger than 1 are reported.
Statistical Procedures
All statistical analyses were performed with the software JMP version 9.0.0 (SAS Institute, Cary, North Carolina, USA). In order to compare the extent to which each plumage trait was associated with geographic coordinates and the principal components of weather, we standardized each plumage trait with the formula ([value − mean]/standard deviation). To test for phenotypic correlations between plumage traits, we performed ANCOVAs. A significant interaction between the factor sex and plumage trait (covariate) on the phenotypically correlated plumage trait (dependent variable) indicates that the magnitude of phenotypic correlation between 2 plumage traits differs between males and females. To compare sexual dimorphism in the 3 plumage traits, we standardized plumage trait values across the 2 sexes. We then ran linear mixed models with standardized plumage traits as dependent variables and individual identity as a random variable, since each individual appears 3 times in the analysis. Sex and trait identity (reddish coloration, number of spots, spot diameter) were 2 factors, and a significant interaction between them indicated that sexual dimorphism differed between traits. To investigate geographic variation in the degree of melanin-based plumage traits, we ran an ANCOVA for each trait including sex as a factor and latitude and longitude plus their interaction as covariates; we added latitude2 and longitude2 to test for nonlinear relationships. To examine the relationships between temperature, rain, and plumage traits, we introduced the 4 principal components of monthly precipitation and ambient temperature as independent variables in ANCOVAs, including sex as a factor and each plumage trait as a dependent variable in separate analyses. Tests are 2-tailed, and P-values lower than 0.05 are considered significant.
RESULTS
Phenotypic Correlation between Melanin-based Plumage Traits
The phenotypic correlation between spot diameter and number of spots was significantly stronger in males than in females (Pearson's r = 0.61 and 0.27, respectively; ANCOVA with spot diameter as dependent variable, sex: F1,1323 = 253.13, P < 0.001; number of spots: F1,1323 = 308.79, P < 0.001; interaction: F1,1323 = 35.50, P < 0.001), as was the phenotypic correlation between pheomelanin-based color and spot diameter (r = 0.41 and 0.27, respectively; interaction: F1,1364 = 6.37, P = 0.012). In contrast, the phenotypic correlation between coloration and number of spots was not significantly different between males and females (r = 0.49 and 0.36, respectively; interaction: F1,1323 = 1.54, P = 0.22).
Sexual Dimorphism in Melanin-based Plumage Traits
The most strongly sexually dimorphic trait was reddish coloration, followed by spot diameter and number of spots (Table 2). Compared with males, females were on average darker red (females: −4.23 ± 0.05; males: −6.36 ± 0.06) and displayed more black spots (32.0 ± 0.4 vs. 27.2 ± 0.4), which were also larger (2.41 ± 0.02 vs. 1.93 ± 0.02 mm). In a linear mixed model with individual identity as random variable (explaining 41% of the variation), there was a significant interaction between sex and trait (sex: F1,1370 = 481.70, P < 0.001; trait: F1,2706 = 0.67, P = 0.51; interaction: F1,2706 = 93.43, P < 0.001), indicating that sexual dimorphism differed between plumage traits. By performing similar analyses with pairs of plumage traits (number of spots × reddish coloration, number of spots × spot diameter, spot diameter × reddish coloration), we found that the degree of sexual dimorphism differed significantly among all traits (P-values for all three interactions between sex and trait were < 0.001; Figure 1).
TABLE 2.
Results of ANCOVAs testing for geographic variation in three melanin-based plumage traits in North American Barn Owls. Analyses were performed on standardized trait values (i.e. [value − mean]/standard deviation) to facilitate comparisons of estimates. Significant results are highlighted in bold font.
Geographic Variation in Melanin-based Plumage Traits
In the northern parts of the North American continent, Barn Owls were darker red and displayed more and larger black spots than in the southern parts of the continent (Table 2, Figures 2 and 3). Owls were darker red and displayed fewer spots on the West Coast than on the East Coast (Table 2, Figures 2 and 3). These relationships were not linear (Table 2). With respect to latitude, owls were paler red and displayed fewer and smaller black spots in the center of the continent than in the southern and northern parts (Figure 3). With respect to longitude, owls exhibited more and larger black spots in the center of the continent than on the West and East coasts (Figure 3).
Climate
Reddish coloration was significantly negatively related to the first principal component of weather (Figure 4A), and there was a tendency for a negative association with the fourth principal component (Tables 1 and 3). This indicates that owls were darker red in colder regions, particularly in summer.
TABLE 3.
Covariation between indices of climate and 3 melanin-based plumage traits in North American Barn Owls. Indices of climate are given by the first 4 principal components of a principal components analysis including monthly values of ambient temperature and amount of precipitation (Table 1). Analyses were performed on standardized plumage trait values across the 2 sexes (i.e. [value − mean]/standard deviation) to facilitate comparisons of estimates, because none of the 2-way interactions between sex and the principal components of weather were significant (i.e. the relationships between the principal components of weather and plumage traits did not differ between males and females). Significant results are highlighted in bold font.
The number of black spots was significantly positively associated with the first principal component of weather and was negatively associated with the third component (Table 3, Figure 4). This means that owls displayed more black spots in regions where it was warmer and where it rained less in winter but more in summer.
Spot size was negatively associated with the first and third principal components of weather (Table 3, Figure 4). This indicates that owls displayed larger black spots in regions where it was colder, particularly in winter, and where it rained less in winter but more in summer.
For all 3 plumage traits, interactions between sex and the 4 principal components of weather were not significant (P-values > 0.06), and hence are not reported in the final models (Table 3). Furthermore, when we modeled all 3 plumage traits together rather than separately, we obtained qualitatively similar results (not shown). This indicates that the relationship between climate and a given plumage trait was not confounded by the other phenotypically correlated plumage traits.
DISCUSSION
In North American Barn Owls, the three melanin-based plumage traits were associated with ambient temperature and precipitation. Owls were darker red and displayed larger but fewer black spots in colder regions. Furthermore, owls tended to display fewer and larger black spots in regions where the climate was rather cold (Figure 4A) and dry in winter but wet in summer. Two sources of noise may have prevented us from detecting stronger relationships, or may have hidden existing relationships. First, because skin specimens were collected between 1844 and 2013, whereas climatic values were available for the period 1950–2000, for each specimen we assigned the mean climatic value for the period 1950–2000 instead of the actual value for the year when an individual was collected. Second, in northern parts of the North American continent, Barn Owls appear to be (partially) migratory (Stewart 1952, Bolen 1978, Bunn et al. 1982). Unfortunately, we cannot evaluate the impact of these two drawbacks on our analyses, and hence further work is required on geographic variation in North American Barn Owls.
Sexual Dimorphism and Phenotypic Correlations
As in Europe, sexual dimorphism of Barn Owls in North America was more pronounced in pheomelanin-based coloration than in the size of eumelanic spots, which was itself more pronounced than sexual dimorphism in the number of eumelanic spots (A. Roulin personal observation). A major issue is whether sexual dimorphism is related to the degree of sex-specific selection, i.e. whether selection on females to be more melanistic than males is much more pronounced with respect to pheomelanin-based coloration than spot size and spot number. Because none of the interactions between sex and the principal components of weather were significant (Table 3), it seems that climate (or a correlated variable) has similar selective impact on plumage for males as for females. As in Europe, phenotypic correlations between plumage traits were stronger in males than in females on the northern American continent (Roulin 2004a). These results suggest that the quantitative and molecular genetics of melanin-based plumage traits are globally similar in European and North American Barn owls.
Climate and Melanin-based Plumage Traits
A number of studies have demonstrated a link between melanin-based plumage traits and survival along a gradient of temperature. In two species of owl, Eurasian Scops-Owl (Otus scops) and Tawny Owl, the dark reddish morph increased in frequency during the last century in Europe, apparently as a consequence of climate warming (Galeotti et al. 2009, Karell et al. 2011). These observations suggest that less melanistic individuals survive better in cold regions and more melanistic individuals in warm regions. This interpretation is consistent with a comparative study of owls at the worldwide scale. In the northern hemisphere, owls of the genus Strix, Glaucidium, Megascops, and Otus are indeed paler poleward and darker red near the equator, as predicted by Gloger's rule (Roulin et al. 2011). Color-specific performance under different ambient temperatures has also been found in Pied Flycatchers: in cold ambient temperatures for rearing, nestlings had higher survival when their father was pale rather than dark, whereas in warm temperatures the opposite tendency was detected (Sirkiä et al. 2010). Thus, fitness differences between dark and pale flycatchers appear to be stronger in cold than warm ambient temperatures. It remains unclear whether these results are due to selection being directly exerted on plumage traits or because of genetically correlated traits. For instance, genes of the melanocortin system not only regulate melanogenesis, but also other phenotypes such as energy homeostasis and thermoregulation (Ducrest et al. 2008).
Barn Owls were darker red (and displayed larger black spots) in the northern parts of the American continent (Figure 3), in contrast to the pattern observed among owl species (Roulin et al. 2011). This may suggest that selective agents exerted on reddish coloration in the Barn Owl are not the same as those exerted on other owl species. If cold ambient temperatures may indeed select darker reddish over lighter reddish Barn Owls (Figure 4A), predator–prey interactions may also play a role. In both Switzerland and Israel, dark reddish owls are mainly found in fields, and white owls close to forests, something that could explain why dark and pale reddish Barn Owls have a different diet in these two countries (Roulin 2004b, Charter et al. 2012, Dreiss et al. 2012). Because Barn Owls in the USA consume more shrews and fewer rodents in regions where it rains more (Clark and Bunck 1991), it is possible that the association between reddish coloration and climate may be driven by predator–prey relationships. If the predator–prey relationship is the main selective agent in the Barn Owl, whereas climate exerts the most selective pressure in other owls, this could explain why cline variation in the degree of pheomelanin-based coloration in the Barn Owl does not follow the general pattern observed in other owls.
Interestingly, we detected the opposite relationship between the first principal component of weather and the number of black spots vs. reddish coloration and spot size (Figure 4A). In colder regions, Barn Owls displayed fewer spots but larger spots and a darker reddish plumage, which is particularly remarkable given that the number of spots was positively correlated with reddish coloration and spot size. These phenotypic correlations imply that selection toward displaying fewer spots in colder regions must be particularly strong, in order to counteract indirect selection due to dark reddish and large-spotted owls being favored in cold climates. A potential scenario to explain this finding is that genes that generate variation in the number of spots pleiotropically regulate processes associated with thermoregulation other than plumage traits, and hence selection on these correlated thermoregulatory traits constrains selection on the number of spots due to Gloger's rule. Indeed, unpublished results for European Barn Owls have shown that body temperature and thermoregulation are associated with the number of black spots, but not with reddish coloration or spot size. Whatever the reason for spot size and spot number not being similarly correlated with climate, our results suggest that the adaptive function of these two plumage traits is not the same.
Conclusion and Perspectives
In North American Barn Owls, geographic variation in melanin-based plumage traits is consistent with Gloger's rule with respect to the number of black spots, but not with regard to pheomelanin-based coloration and spot size. Patterns of geographic variation may not always follow Gloger's rule because selective factors other than climate may act on melanin-based plumage traits, or because genes for plumage traits may also regulate physiological or behavioral processes related to thermoregulation other than plumage characteristics. Experiments are now required to measure how differently colored Barn Owls cope with various temperatures. For instance, Mosher and Henny (1976) measured oxygen consumption in differently colored Eastern Screech-Owls (Megascops asio) at different ambient temperatures. A similar design could be used with hatchlings swapped between nests, in order to measure the performance of individuals raised by foster parents in relation to plumage traits assessed in the nestlings themselves and in their biological and foster parents. In this way, it would be possible to disentangle genetic from environmental effects on the relationship between plumage traits and the ability to cope with cold and warm temperatures. A relationship with the coloration measured in the foster parents would indicate that nestlings raised by differently colored parents experienced different rearing conditions that altered their performance under different temperatures. A relationship with the coloration measured in the biological parents would indicate that plumage traits are genetically correlated with physiological or behavioral process related to thermoregulation. Our study should therefore be considered as just a first attempt to identify the adaptive value of melanin-based plumage traits in North American Barn Owls.
ACKNOWLEDGMENTS
We warmly thank all curators who opened their museums to us to measure Barn Owl skins for their hospitality. We hope that this study helps to demonstrate how important museum collections are to tackling scientific questions. We are grateful to Reto Burri, Sylvain Antoniazza, and two anonymous reviewers for useful comments on a previous version of the manuscript.
Funding statement. This study was financed by the Swiss National Science Foundation, the American Natural History Museum in New York, the Swiss Academy of Sciences, the Basler Stiftung für biologische Forschung, and the “Fondation du 450ème anniversaire” and “Fondation Agassiz” of the University of Lausanne. No funders had any input into the manuscript content, nor required approval of the manuscript prior to publication.