# An interspecific assessment of Bergmann’s rule in 22 mammalian families

- Jostein Gohli
^{1}Email authorView ORCID ID profile and - Kjetil L. Voje
^{2}

**16**:222

https://doi.org/10.1186/s12862-016-0778-x

© The Author(s). 2016

**Received: **30 May 2016

**Accepted: **29 September 2016

**Published: **19 October 2016

## Abstract

### Background

Bergmann’s rule proposes that animals in cold habitats will be larger than those in warm habitats. This prediction has been tested thoroughly at the intraspecific level, but few studies have investigated the hypothesis with interspecific data using phylogenetic comparative approaches. Many clades of mammals have representatives in numerous distinct biomes, making this order highly suitable for a large-scale interspecific assessment of Bergmann’s rule. Here, we evaluate Bergmann’s rule within 22 mammalian families—with a dataset that include ~35 % of all described species—using a phylogenetic comparative approach. The method is based on an Ornstein-Uhlenbeck model of evolution that allows for joint estimation of adaptation and constraints (phylogenetic inertia) in the evolution of a trait. We use this comparative method to investigate whether body mass evolves towards phenotypic optima that are functions of median latitude, maximum latitude or temperature. We also assess the closely related Allen’s rule in five families, by testing if relative forelimb length evolves as a function of temperature or latitude.

### Results

Among 22 mammalian families, there was weak support for Bergmann’s rule in one family: A decrease in temperature predicted increased body mass in Canidae (canids). We also found latitude and temperature to significantly predict body mass in Geomyidae (pocket gophers); however, the association went in the opposite direction of Bergmann’s predictions. Allen’s rule was supported in one of the five examined families (Pteropodidae; megabats), but only when forelimb length evolves towards an optimum that is a function of maximum latitude, not median latitude or temperature.

### Conclusions

Based on this exhaustive assessment of Bergmann’s rule, we conclude that factors other than latitude and temperature are the major drivers of body mass evolution at the family level in mammals.

## Keywords

## Background

Biologists have always sought ecological and evolutionary generalizations that structure the rich and complex diversity of life. Bergmann’s rule is such an eco-evolutionary generalization stating that animals will be larger in cold climates and smaller in warm climates [1]. Given latitude’s association with temperature, Bergmann’s hypothesis predicts a positive association between body mass and latitude. One might think that such a simple idea—increased body size with latitude—would be easily confirmed or discarded. However, Bergmann’s hypothesis continues to receive a great deal of attention more than 150 years after it was originally proposed.

Bergmann proposed his rule as applying to closely related species within the same genus (see Blackburn et al. [2] for a review of definitions of Bergmann’s rule), and it seems that Bergmann’s rule can operate at both the intraspecific and interspecific level [3, 4]. However, it can be argued that the expected correlation between body size and temperature—given that Bergmann’s hypothesis is correct—may differ in strength among species or populations. For example, there might be few biological differences among populations of the same species at different latitudes. Thus, if temperature clines have an effect on body mass evolution, they may be readily detectable among populations, since most other variables are the same. Conversely, temperature may describe a comparatively smaller amount of interspecific variation among species, since (many) species are different in their general biology, and will thus differ in important constraints or selective pressures on body mass. Larger amounts of residual variation in interspecific studies of body mass and temperature might therefore be expected, but species specific median or maximum latitudes may still explain a substantial amount of variation in body mass.

If we consider only endotherms—which some deem appropriate [4]—Bergmann’s rule has been thoroughly evaluated on the intraspecific level: Studies have shown that temperature or latitude is associated with body mass within species of both mammals e.g., [5–7] and birds e.g., [6, 8]. These observed positive correlations between latitude and body size within species suggests that a similar pattern could also be found among species. However, we are only aware of a handful of studies that evaluate Bergmann’s rule interspecifically in endotherms, and their collective results are inconclusive: While some studies find strong evidence [9–11], others find no evidence for Bergmann’s rule across species [12, 13].

One reason for the conflicting evidence for Bergmann’s rule at the interspecific level may be due to the different statistical models used to investigate the hypothesis. Blackburn and Hawkins [10] and Rodríguez et al. [9] investigated Bergmann’s rule using standard linear regression models which assumes that all species constitute independent data points and that shared evolutionary history among species has no effect on the data. Kamilar et al. [13] and Blackburn and Gaston [11] used Felsenstein’s [14] independent contrasts approach in their evaluation of Bergmanns’s rule in Malagasy Strepsirhines, a method that controls for phylogenetic effects by assuming that any statistical influence of the phylogeny is the result of ancestry [15]. Another comparative study by Diniz-Filho and colleges [16] use what they call Phylogenetic Eigenvector Regression (PVR) to estimate how much variation in body mass that can be due to ecological factors. None of these tests of Bergmann’s rule assume an underlying model of adaptive evolution that can disentangle how much of the estimated phylogenetic effect that stems from shared ancestry and how much that is due to adaptive evolution of body mass towards niches that are placed non-randomly on the phylogeny [17].

With the recent accumulation of high quality phylogenies, the stage is set for a large scale interspecific assessment of Bergmann’s rule in mammals. Our hypothesis is that of Bergmann: that body mass will increase as an adaptation to increased absolute latitude, and increase in response to decreased temperature. We test the hypothesis for 22 mammalian families using a phylogenetic comparative approach that allows for the joint estimation of phylogenetic inertia and adaptation in body mass based on an Ornstein-Uhlenbeck process [18–20]. More specifically, we investigate whether body mass evolves as a consequence of changes in an optimum that is modeled as a function of various predictor variable (average temperature and median and max absolute latitude).

For five mammalian families, data also allowed us to test Allen’s rule [21], an ecogeographic pattern related to Bergmann’s rule, which states that the relative size of body extremities of endotherms (e.g., limbs, tails, ears etc.) should be smaller in colder environments in order to reduce thermoregulatory costs. The level of support for Allen’s rule in the literature is comparable to Bergmann’s: Allen’s rule has been found to hold for some intraspecific temperature clines in mammals e.g., [22, 23], poikiloterms [24] and birds [25, 26], while we are only aware of two interspecific analyses, on birds [27, 28], which both found support for the rule.

## Methods

Results from regressing ln (body mass) on maximum-, median latitude, or temperature

Family | Predictor |
| Phylogentic half-life | Stationary variance |
| Optimal regression slope | AICc | AICc- |
---|---|---|---|---|---|---|---|---|

− | 31.50 | 30.00 | − | − | 186.25 | − | ||

Bovidae | median latitude | 81 | 32.50 | 29.00 | 0.030 | 0.983 | 186.11 | −0.14 |

max latitude | 32.50 | 29.00 | 0.040 | 1.037 | 185.29 | −0.96 | ||

temperature | 31.00 | 29.00 | 0.019 | −1.180 | 186.99 | 0.73 | ||

Median latitude - temperature regression: | ||||||||

− | 2.41 | 1.40 | − | − | 52.12 | − | ||

Canidae | median latitude | 25 | 11.61 | 5.30 | 0.076 | 0.390 | 53.06 | 0.94 |

max latitude | 4.41 | 2.10 | 0.085 | 0.103 | 52.72 | 0.60 | ||

temperature | 10.81 | 3.70 | 0.233 | −0.935 | 48.27 | −3.85 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 200.10 | − | − | 15.20 | − | ||

Cebidae | median latitude | 35 | ∞ | 180.10 | 0.040 | 23.736 | 16.49 | 1.29 |

max latitude | ∞ | 180.10 | 0.035 | 20.376 | 16.66 | 1.46 | ||

temperature | ∞ | 180.10 | 0.037 | −56.352 | 16.58 | 1.38 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 195.10 | − | − | 55.34 | − | ||

Cercopithecidae | median latitude | 71 | ∞ | 104.10 | 0.023 | 8.854 | 55.94 | 0.61 |

max latitude | ∞ | 152.10 | 0.076 | 25.809 | 51.91 | −3.42 | ||

temperature | ∞ | 136.10 | 0.039 | −36.621 | 54.73 | −0.60 | ||

Median latitude - temperature regression: | ||||||||

− | 0.41 | 0.80 | − | − | 62.38 | − | ||

Cervidae | median latitude | 29 | 0.46 | 0.85 | 0.009 | −0.005 | 64.90 | 2.51 |

max latitude | 0.31 | 0.70 | 0.031 | 0.007 | 64.41 | 2.02 | ||

temperature | 0.31 | 0.70 | 0.036 | −0.020 | 64.26 | 1.88 | ||

Median latitude - temperature regression: | ||||||||

− | 2.30 | 1.45 | − | − | 616.80 | − | ||

Cricetidae | median latitude | 415 | 2.70 | 1.70 | 0.000 | −0.005 | 618.81 | 2.01 |

max latitude | 3.00 | 1.85 | 0.004 | −0.024 | 617.22 | 0.42 | ||

temperature | 2.70 | 1.70 | 0.000 | −0.002 | 618.89 | 2.08 | ||

Median latitude - temperature regression: | ||||||||

− | 1.15 | 1.80 | − | − | 136.22 | − | ||

Dasyuridae | median latitude | 51 | 1.38 | 2.04 | 0.010 | 0.033 | 138.06 | 1.84 |

max latitude | 1.74 | 2.40 | 0.032 | 0.063 | 137.02 | 0.80 | ||

temperature | 1.29 | 1.86 | 0.041 | −0.142 | 136.46 | 0.24 | ||

Median latitude - temperature regression: | ||||||||

− | 1.67 | 0.59 | − | − | 69.93 | − | ||

Echimyidae | median latitude | 49 | 1.22 | 1.22 | 0.008 | −0.029 | 71.98 | 2.05 |

max latitude | 3.68 | 1.15 | 0.010 | 0.076 | 71.91 | 1.98 | ||

temperature | 1.46 | 0.52 | 0.014 | 0.109 | 71.63 | 1.70 | ||

Median latitude - temperature regression: | ||||||||

− | 43.40 | 12.10 | − | − | 55.29 | − | ||

Emballonuridae | median latitude | 34 | 43.80 | 12.10 | 0.004 | −0.587 | 57.72 | 2.44 |

max latitude | 44.60 | 12.10 | 0.019 | 0.964 | 57.19 | 1.91 | ||

temperature | 43.80 | 12.10 | 0.010 | 4.156 | 57.53 | 2.24 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 1488.10 | − | − | 91.82 | − | ||

Felidae | median latitude | 29 | ∞ | 1440.10 | 0.030 | −43.326 | 93.71 | 1.88 |

max latitude | ∞ | 1320.10 | 0.055 | 53.462 | 92.99 | 1.17 | ||

temperature | ∞ | 1520.10 | 0.002 | 17.582 | 94.47 | 2.64 | ||

Median latitude - temperature regression: | ||||||||

− | 0.39 | 0.46 | − | − | 64.03 | − | ||

Geomyidae | median latitude | 32 | 0.01 | 0.23 | 0.366 | −0.035 | 54.29 | −9.74 |

max latitude | 0.01 | 0.23 | 0.380 | −0.030 | 53.55 | −10.47 | ||

temperature | 0.01 | 0.29 | 0.237 | 0.055 | 60.21 | −3.82 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 390.10 | − | − | 70.63 | − | ||

Heteromyidae | median latitude | 51 | ∞ | 375.10 | 0.009 | −20.960 | 72.54 | 1.91 |

max latitude | ∞ | 360.10 | 0.050 | −38.684 | 70.39 | −0.24 | ||

temperature | ∞ | 345.10 | 0.045 | 71.262 | 70.65 | 0.02 | ||

Median latitude - temperature regression: | ||||||||

− | 1.36 | 0.55 | − | − | 57.47 | − | ||

Leporidae | median latitude | 45 | 1.45 | 0.58 | 0.010 | −0.014 | 59.42 | 1.95 |

max latitude | 1.81 | 0.66 | 0.037 | −0.029 | 58.19 | 0.71 | ||

temperature | 1.21 | 0.50 | 0.011 | 0.024 | 59.41 | 1.94 | ||

Median latitude - temperature regression: | ||||||||

− | 0.55 | 0.60 | − | − | 98.63 | − | ||

Macropodidae | median latitude | 46 | 1.21 | 0.90 | 0.068 | −0.068 | 98.59 | −0.05 |

max latitude | 0.49 | 0.58 | 0.003 | 0.007 | 100.94 | 2.31 | ||

temperature | 0.97 | 0.74 | 0.071 | 0.098 | 98.01 | −0.62 | ||

Median latitude - temperature regression: | ||||||||

− | 0.01 | 0.50 | − | − | 94.02 | − | ||

Molossidae | median latitude | 41 | 0.01 | 0.50 | 0.016 | −0.009 | 95.86 | 1.83 |

max latitude | 0.01 | 0.50 | 0.017 | −0.009 | 95.80 | 1.78 | ||

temperature | 0.01 | 0.50 | 0.017 | 0.033 | 95.80 | 1.78 | ||

Median latitude - temperature regression: | ||||||||

− | 0.75 | 1.60 | − | − | 786.97 | − | ||

Muridae | median latitude | 324 | 0.80 | 1.70 | 0.000 | 0.002 | 789.02 | 2.05 |

max latitude | 0.80 | 1.70 | 0.000 | 0.005 | 788.90 | 1.93 | ||

temperature | 0.80 | 1.70 | 0.001 | 0.016 | 788.71 | 1.75 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 990.70 | − | − | 112.15 | − | ||

Mustelidae | median latitude | 42 | ∞ | 975.70 | 0.002 | −3.630 | 114.53 | 2.38 |

max latitude | ∞ | 975.70 | 0.001 | 2.613 | 114.55 | 2.41 | ||

temperature | ∞ | 975.70 | 0.003 | 9.162 | 114.48 | 2.33 | ||

Median latitude - temperature regression: | ||||||||

− | ∞ | 945.70 | − | − | 97.62 | − | ||

Nesomyidae | median latitude | 32 | ∞ | 900.70 | 0.004 | −22.286 | 100.15 | 2.52 |

max latitude | ∞ | 850.70 | 0.048 | −51.197 | 98.82 | 1.19 | ||

temperature | ∞ | 925.70 | 0.007 | −42.047 | 100.04 | 2.42 | ||

Median latitude - temperature regression: | ||||||||

− | 1.71 | 1.20 | − | − | 184.97 | − | ||

Phyllostomidae | median latitude | 101 | 1.66 | 1.15 | 0.000 | 0.006 | 187.13 | 2.15 |

max latitude | 1.96 | 1.30 | 0.013 | −0.050 | 185.82 | 0.85 | ||

temperature | 1.66 | 1.15 | 0.009 | −0.157 | 186.23 | 1.25 | ||

Median latitude - temperature regression: | ||||||||

− | 0.15 | 1.18 | − | − | 134.12 | − | ||

Pteropodidae | median latitude | 52 | 0.13 | 1.02 | 0.069 | 0.025 | 136.36 | 2.24 |

max latitude | 0.16 | 0.96 | 0.143 | 0.034 | 132.27 | −1.85 | ||

temperature | 0.13 | 1.02 | 0.065 | −0.126 | 136.71 | 2.59 | ||

Median latitude - temperature regression: | ||||||||

− | 1.41 | 2.50 | − | − | 418.47 | − | ||

Sciuridae | median latitude | 196 | 1.21 | 2.20 | 0.005 | 0.018 | 419.71 | 1.24 |

max latitude | 1.21 | 2.20 | 0.000 | 0.001 | 420.58 | 2.11 | ||

temperature | 1.61 | 2.80 | 0.000 | 0.002 | 420.58 | 2.11 | ||

Median latitude - temperature regression: | ||||||||

− | 0.16 | 0.45 | − | − | 155.53 | − | ||

Vespertilionidae | median latitude | 91 | 0.16 | 0.45 | 0.041 | 0.002 | 158.77 | 3.24 |

max latitude | 0.16 | 0.45 | 0.042 | −0.003 | 158.66 | 3.13 | ||

temperature | 0.16 | 0.45 | 0.042 | −0.006 | 158.72 | 3.20 | ||

Median latitude-temperature regression: |

We investigate whether log transformed body mass (Bergmann’s rule) and relative forelimb length (forelimb length controlled for body mass; Allen’s rule) have evolved towards optima that are influenced by latitude or temperature within different families of mammals. We did this using a phylogenetic comparative approach implemented in the R program SLOUCH, designed to study adaptive evolution of a trait along a phylogenetic tree [17–20, 34]. The output of the model can be summarized by a regression, which includes information on whether the analyzed trait is evolving towards the estimated optima, how fast (or slow) the trait approaches the optimum, and how much of the trait variation is explained by adaptation towards the optimum. The model of evolution in SLOUCH is based on an Ornstein–Uhlenbeck model and assumes that the trait (body mass and relative forelimb length in our case) has a tendency to evolve towards a ‘primary’ optimum *Θ*, defined as the average optimal state that species will reach in a given environment when ancestral constraints have disappeared [18]. The primary optimum is modeled as a linear function of the predictor variable, which evolves as if by a Brownian-motion process. Lag in adaptation towards primary optima is quantified by a half-life parameter, *t*
_{
1/2
} = ln (2)/*α*, which can be interpreted as the average time it takes a species to evolve half the distance from the ancestral phenotype towards the predicted optimal phenotype. A half-life of zero means there is no evolutionary lag, while a half-life above zero indicates that adaptation is not immediate.

The SLOUCH model returns an “optimal regression”, which is the best fit of the estimated primary optimum on the response variable (e.g. ln body mass). In tests of Bergmann’s rule, this optimal regression describes the expected relationship between ln (body mass) and the predictor in the model (e.g. latitude) if there were no constraints on the evolution of body mass towards the optimal state (instantaneous adaptation). If Bergmann’s rule applies, the optimal regression coefficient would be positive in models of log body mass and the two latitude variables, while it will be negative for the optimal regression of log body mass and temperature. The optimal regression is contrasted with an “evolutionary regression”, which is the best fit of the predictor variable on the response variable. The evolutionary regression represents the observed relationship between the variables and is shallower than the optimal regression whenever there is a lag in adaptation (i.e. when the half-life is not zero). The model of evolution implemented in SLOUCH also includes a stochastic component with standard deviation *σ,* which can be interpreted as evolutionary changes in body mass due to unmeasured selective forces and genetic drift. This component of the model is reported as v_{y} = *σ*
^{2}/2*α*, which signify the expected residual variance when adaptation and stochastic changes have come to an equilibrium. Generalized least squares is used for estimation of the regression parameters (i.e., the influence of the predictor on the primary optimum) and maximum likelihood for estimation of *α* and *σ*
^{2} in an iterative procedure. For a full description of the model implemented in SLOUCH, see Hansen et al. [20].

Models that include a predictor variable are referred to as adaptation models since we in these models test whether body mass (in tests of Bergmann’s rule) or relative forelimb length (in tests of Allen’s rule) evolves towards optima influenced by latitude or temperature. The adaptation models are contrasted with an intercept-only model without a predictor variable. The half-life parameter in such intercept-only models is a measure of the phylogenetic effect in the response variable, which is an estimate of how well the phylogeny explains the distribution of body mass or relative forelimb length in the investigated family. A half-life of zero in such a model means the response variable is not phylogenetically structured, while a half-life > 0 indicates that there exists an influence of phylogeny on the data. A half-life value larger than 30 times the length of the phylogeny is reported as infinity, as the OU model reduces to a Brownian motion when the half-life is very large. A phylogenetic effect can be due to slowness of adaptation, adaptation towards phylogenetically structured optima, or a combination of both. The adaptation models can determine what proportion of the phylogenetic effect within each family that can be accounted for by adaptation towards optima influenced by latitude and temperature, respectively. Adaptation models are compared to the intercept-only models using the small sample-size corrected version of Akaike information criterion (AICc); adaption models that have AICc scores that are two or more units lower than their respective intercept-only models are considered substantially better [35]. To judge whether better models (according to the AICc score) support Bergmann’s rule, we interpret the slope of the optimal regression together with the amount of variation in body mass that the optimal regression explains. All statistical analyses were done in R v3.1.3 [36].

## Results

### Tests of Bergmann’s rule

The phylogenetic effect in body mass varied among the examined families, but were generally large (Table 1). 18 out of 22 families had very strong phylogenetic effects in body mass (half-life > 0.5). We found median latitude and temperature to be strongly negatively correlated in all 22 families, except for Nesomyidae (Table 1).

*r*

^{2}, and AICc-

*θ*AICc values are given in Fig. 2. Many models had steep slopes due to very large half-life values, which means body mass shows no tendency to evolve as a response to changes in the predictor in these models and that the slope estimate is not meaningful. Also, the majority of these steep slopes have a 95 % CI covering zero due to large standard errors. The negative association between temperature and body mass in Canidae (canines) is the strongest evidence in favor of Bergmann’s rule of all the models we tested (

*r*

^{2}= 0.23; AICc-

*θ*AICc = −3.9; Fig. 3; Fig. 2; Table 1). The half-life estimate of this model was more than ten times the length of the phylogeny, which means there is no tendency for body mass to evolve towards the optimum and that the model residuals changes similar to a Brownian motion, possibly with a trend [18, 37]. Body mass in Cercopithecidae (Old World monkeys) also showed a positive relationship with maximum latitude, but the coefficient of determination was small (

*r*

^{2}= 0.08), which means that maximum latitude has marginal predictive power on body size variation within this family. The strongest relationship between body mass and a predictor variable was found for Geomyidae (pocket gophers); the model using maximum latitude as a predictor explained 39 % of the variation in body mass (AICc-

*θ*AICc = −10.5; Fig. 3; Table 1), and the model using median latitude as predictor performed qualitatively similar. Importantly, Geomyidae exhibited a negative association between body mass and latitude. Half-life estimates were low for both models (<0.06), indicating rapid adaptation in body mass towards the optimum. The optimal regression slope of Geomyidae’s maximum latitude model was − 0.030, which predicts a 3 % decrease in body mass for every increase of one latitudinal degree.

### Tests of Allen’s rule

*r*

^{2}= 0.16; AICc-

*θ*AICc = −3.0; Fig. 4; Table 2). Half-life estimates indicate that adaptation of forelimb length as a function of maximum latitude was relatively rapid in Pteropodidae (half-life = 0.15). For each increase of one latitudinal degree, the model predicted a decrease of ~2.5 % in forelimb length. In all the five datasets (families) where we tested Allen’s rule, median latitude and temperature was strongly and negatively correlated (Table 2).

Results from regressing ln (relative forelimb length) on maximum-, median latitude, and temperature

Family | Predictor |
| Phylogentic half-life | Stationary variance |
| Optimal regression slope | AICc | AICc- |
---|---|---|---|---|---|---|---|---|

− | ∞ | 168.10 | − | − | 38.64 | − | ||

Emballonuridae | median latitude | 34 | ∞ | 160.10 | 0.004 | 9.693 | 41.09 | 2.45 |

max latitude | ∞ | 160.10 | 0.025 | −18.349 | 40.36 | 1.71 | ||

temperature | ∞ | 160.10 | 0.015 | −88.544 | 40.70 | 2.05 | ||

Median latitude - temperature regression: | ||||||||

− | 0.01 | 0.25 | − | − | 66.48 | − | ||

Molossidae | median latitude | 40 | 0.01 | 0.25 | 0.020 | 0.008 | 68.13 | 1.65 |

max latitude | 0.01 | 0.25 | 0.082 | 0.014 | 65.52 | −0.96 | ||

temperature | 0.01 | 0.25 | 0.048 | −0.039 | 66.95 | 0.47 | ||

Median latitude - temperature regression: | ||||||||

− | 0.73 | 0.43 | − | − | 141.60 | − | ||

Phyllostomidae | median latitude | 97 | 0.66 | 0.40 | 0.015 | 0.021 | 142.31 | 0.71 |

max latitude | 0.86 | 0.45 | 0.031 | 0.033 | 140.81 | −0.79 | ||

temperature | 0.66 | 0.40 | 0.001 | −0.017 | 143.72 | 2.12 | ||

Median latitude - temperature regression: | ||||||||

− | 0.13 | 0.50 | − | − | 93.80 | − | ||

Pteropodidae | median latitude | 52 | 0.13 | 0.45 | 0.068 | −0.012 | 95.95 | 2.15 |

max latitude | 0.15 | 0.40 | 0.163 | −0.025 | 90.80 | −3.00 | ||

temperature | 0.13 | 0.45 | 0.074 | 0.100 | 95.72 | 1.92 | ||

Median latitude - temperature regression: | ||||||||

− | 0.10 | 0.22 | − | − | 105.53 | − | ||

Vespertilionidae | median latitude | 88 | 0.06 | 0.18 | 0.100 | −0.004 | 105.80 | 0.27 |

max latitude | 0.08 | 0.20 | 0.070 | 0.001 | 106.63 | 1.11 | ||

temperature | 0.06 | 0.18 | 10.11 | 0.010 | 105.69 | 0.17 | ||

Median latitude - temperature regression: |

## Discussion

Here, we have performed tests of Bergmann’s rule in 22 mammal families. To do this, we used a phylogenetic comparative method to test if log body mass evolved as a response to changes in median latitude, maximum latitude or temperature. The dataset includes ~35 % of all described species and covers eight orders, which makes this a comprehensive interspecific test of Bergmann’s rule in mammals. The general pattern was one of very little support for Bergmann’s rule, i.e., of latitude or temperature being important predictors of body mass evolution. Most of the models in SLOUCH that included a predictor variable were not better than a model without a predictor according to their AICc score and most did not have a slope estimate that differed significantly from zero (Fig. 2). Also, most models explained very little variation in body mass (Fig. 2; Table 1). We also performed tests of Allen’s rule in 5 families, were we found a similar pattern of little statistical and biological significance in the tested predictor variables (Table 2).

Only one of the 22 families exhibited a somewhat clear negative relationship between temperature and body mass, thus supporting Bergmann’s rule. In this family, Canidae, the temperature model explained about 23 % of the variation in body mass, however, this adaption model was only 4 AICc units better than the model without any predictors. Furthermore, the estimated rate of adaptation in body mass was extremely slow, which indicates that changes in temperature are not necessarily followed by any actual change in body size. Furthermore, given that we ran 66 distinct regression analyses in our tests of Bergmann’s rule, this result may be the product of a type I error. Thus, we caution that further analyses (including all 34 species) of this relationship in Canidae should be conducted before concluding that species within this family follows the prediction from Bergmann’s rule.

We did find a strong association between body mass and latitude in one of the 22 families, namely Geomyidae. Interestingly, the relationship was in the opposite direction of Bergmann’s prediction. The negative association between body mass and latitude in Geomyidae is congruent with Medina et al. [12]: In their study of a rodent genus (*Ctenomys*), which is not included in our data set, body size decreased with increasing latitude. This congruence may be explained by the fact that both *Ctenomys* and Geomyidae are subterranean rodents; the burrowing lifestyle will probably involve some constraints on body mass not found in other rodents or mammals in general. Also, the importance of overland temperature is likely less for animals that spend a lot of time underground [38, 39]. A likely causal driver for the correlation between latitude and body mass in Geomyidae is soil humidity, which seem to affects burrowing; larger species are found in areas with dry, sandy, and brittle soil [12]. This fits the observed pattern since humidity correlates negatively with overland temperature in the latitudinal range of both Geomyidae and *Ctenomys*. However, resource availability is also a possible driver of body mass evolution in Geomyidae and *Ctenomys* [12], which could be a confounding variable in our analyses of body mass in Geomyidae.

While there was no reliable support for Bergmann’s rule among the 22 families (with the potential exception of Canidae), we did detect weak support for Allen’s rule in one family of bats—namely Pteropodidae (megabats; Fig. 4; Table 2). Bats are unique, being the only truly flying mammals, and while they may use crowding and rolling up into a spherical shape in order to limit heat loss while roosting [40], the way in which they are subjected to the environment when flying is fundamentally different to other mammals. Many mammals live underground or under vegetation, and even those that do not, benefit from the insulation of the ground, terrain or vegetation [41]. The wings of bats are excellent tools for dissipating excess heat during flight in warm habitats [42], but are inherently poor at retaining heat. This is in contrast to birds, whose wings are covered in highly insulating feathers. There is some support for Allen’s rule across bird species, but only for featherless limbs [27, 28]. In fact, migratory birds that breed further north/south are expected to have longer wings due to longer migration distance [43]. With their extreme exposure to the elements, and lack of insulation on their wings, and the fact that they have relatively much longer forelimbs than any other mammal group, it is not surprising that increased latitude (i.e., lower temperatures) would inflict a strong selective pressure on wing length in Pteropodidae.

Interestingly, we found no support for Allen’s rule in a second bat family included in our study (Emballonuridae; sac-winged bats), which suggests that the result found for Pteropodidae should be interpreted with care. With large species such as the flying fox (*Pteropus vampyrus*) most commonly representing Pteropodidae, one might intuitively think that Pteropodidae’s exceedingly large wingspan (1.7 m in *P. vampyrus*) is what sets it apart from other groups of bats. And while it is true that Pteropodidae has larger average forelimb length than Emballonuridae (96.6 mm, and 54.9 mm respectively), the average forelimb length to body mass ratio is actually larger in Emballonuridae (5.13, compared to 1.02 in Pteropodidae). The two families also have similar latitudinal distributions, so this variable offers no explanation for the lack of support for Allen’s rule in Emballonuridae. One major difference between the two is that Emballonuridae species are mainly insectivorous and hawk flying insects in flight. It may be that this foraging behavior inflicts strong constraints on wing morphology—that counteracts a selective pressure towards decreasing wing span with increased latitude—which is absent in the frugivorous Pteropodidae.

Bergmann’s original formulation of his rule was about how body size variation in a group of closely related species was related to temperature (see Blackburn et al. [2]). Our investigation of whether Bergmann’s rule applies at the family level may therefore be argued to be outside the taxonomic scope of how the rule was originally formulated. However, analyzing families separately allow us to investigate whether Bergmann’s rule holds true across mammals in general. Our results indicate that temperature and latitude are not universal drivers of body mass variation among mammals.

Our main result, that temperature and latitude do not represent important factors affecting the adaptive landscape of body mass or limb length evolution within families of mammals, comes with some caveats. Performing comparative analyses of interspecific data on species means is associated with several levels of uncertainty. The members of species may show substantial variation in phenotypes and ecology, which implies that the analyzed species median values may not be very representative for all populations of a given species. This could partly have been accounted for in our analyses if variance measures for all variables were available, but this was not the case. For latitude, this issue is largely remedied by the inclusion of both median and maximum latitude in separate tests. However, the issue remains for body mass and temperature.

There is also the issue of whether our data variables are appropriate and exhaustive. For instance, we did not account for altitude in this study. It is likely that altitude and latitude have similar effects on temperature, which would entail that species living only at high elevations experience a climate similar to that at higher latitudes and lower altitudes. The inclusion of temperature as a predictor of body mass acts to control for this issue. Another potential issue is the use of mean temperature, which isn’t necessarily the most relevant metric when investigating temperature’s effect on body mass evolution; perhaps the most extreme temperature experienced by species is the most essential factor affecting body mass evolution. However, there might be a correlation between the extreme and the mean temperatures within the geographic range of most species, which means that some of the relevant variation may be captured by the mean values. Also, extreme temperature values may to some extent be represented by the maximum latitude variable.

## Conclusion

In summation, we found no reliable support for Bergmann’s rule among the 22 examined mammalian families. Further, we found weak support for Allen’s rule in only one very atypical mammalian family, the megabats. We conclude that neither Bergmann’s rule nor Allen’s rule are important interspecific phenomena in mammals at the family level.

## Declarations

### Acknowledgments

We thank Becky Cramer, Inger Maren Rivrud, José Alexandre Felizola Diniz-Filho, and an anonymous reviewer for helpful comments on the manuscript.

### Availability of data and materials

A dataset supporting the conclusions of this article-the PanTHERIA database-is available at http://esapubs.org/archive/ecol/E090/184/. Further data supporting this article are included in the additional file.

### Authors’ contributions

The study was conceived by KLV and designed by JG and KLV. JG carried out the analyses and drafted the first version of the manuscript. JG and KLV finalized the manuscript together. Both authors read and approved the final manuscript.

### Competing interests

The authors declare that they have no competing interests.

### Consent for publication

Not applicable.

### Ethics approval and consent to participate

Not applicable.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Bergmann C. Uber die Verhaltnisse der Warmeokonomie der Thiere zu ihrer Grosse. Gottinger Studien. 1847;3:595–708.Google Scholar
- Blackburn TM, Gaston KJ, Loder N. Geographic gradients in body size: a clarification of Bergmann’s rule. Divers Distrib. 1999;5(4):165–74.View ArticleGoogle Scholar
- Meiri S. Bergmann’s Rule-what’s in a name? Global Ecol Biogeogr. 2011;20(1):203–7.View ArticleGoogle Scholar
- Pincheira-Donoso D. The balance between predictions and evidence and the search for universal macroecological patterns: taking Bergmann’s rule back to its endothermic origin. Theor Biosci. 2010;129(4):247–53.View ArticleGoogle Scholar
- Ashton KG, Tracy MC, Ad Q. Is Bergmann’s rule valid for mammals? Am Nat. 2000;156(4):390–415.View ArticleGoogle Scholar
- Meiri S, Dayan T. On the validity of Bergmann’s rule. J Biogeogr. 2003;30(3):331–51.View ArticleGoogle Scholar
- Rodríguez MÁ, Olalla-Tárraga MÁ, Hawkins BA. Bergmann’s rule and the geography of mammal body size in the Western Hemisphere. Global Ecol Biogeogr. 2008;17(2):274–83.View ArticleGoogle Scholar
- Ramirez L, Diniz-Filho JAF, Hawkins BA. Partitioning phylogenetic and adaptive components of the geographical body-size pattern of New World birds. Global Ecol Biogeogr. 2008;17(1):100–10.Google Scholar
- Rodríguez MÁ, López-Sañudo IL, Hawkins BA. The geographic distribution of mammal body size in Europe. Global Ecol Biogeogr. 2006;15(2):173–81.View ArticleGoogle Scholar
- Blackburn TM, Hawkins BA. Bergmann’s rule and the mammal fauna of northern North America. Ecography. 2004;27(6):715–24.View ArticleGoogle Scholar
- Blackburn TM, Gaston KJ. Spatial patterns in the body sizes of bird species in the New World. Oikos. 1996;77(3):436–46.View ArticleGoogle Scholar
- Medina AI, Martí DA, Bidau CJ. Subterranean rodents of the genus Ctenomys (Caviomorpha, Ctenomyidae) follow the converse to Bergmann’s rule. J Biogeogr. 2007;34(8):1439–54.View ArticleGoogle Scholar
- Kamilar JM, Muldoon KM, Lehman SM, Herrera JP. Testing Bergmann’s rule and the resource seasonality hypothesis in Malagasy primates using GIS-based climate data. Am J Phys Anthropol. 2012;147(3):401–8.View ArticlePubMedGoogle Scholar
- Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.View ArticleGoogle Scholar
- Hansen TF, Orzack SH. Assessing current adaptation and phylogenetic inertia as explanations of trait evolution:The need for controlled comparisons. Evolution. 2005;59(10):2063–72.PubMedGoogle Scholar
- Diniz-Filho JAF, Bini LM, Rodríguez MÁ, Rangel TFLVB, Hawkins BA. Seeing the forest for the trees: partitioning ecological and phylogenetic components of Bergmann’s rule in European Carnivora. Ecography. 2007;30(4):598–608.View ArticleGoogle Scholar
- Labra A, Pienaar J, Hansen T. Evolution of thermal physiology in liolaemus lizards: Adaptation, phylogenetic inertia, and niche tracking. Am Nat. 2009;174(2):204–20.View ArticlePubMedGoogle Scholar
- Hansen TF. Stabilizing selection and the comparative analysis of adaptation. Evolution. 1997;51(5):1341–51.View ArticleGoogle Scholar
- Butler MA, King AA. Phylogenetic comparative analysis: A modeling approach for adaptive evolution. Am Nat. 2004;164(6):683–95.View ArticleGoogle Scholar
- Hansen TF, Pienaar J, Orzack SH. A comparative method for studying adaptation to a randomly evolving environment. Evolution. 2008;62(8):1965–77.PubMedGoogle Scholar
- Allen G. Physiological aesthetics. New York: D. Appleton; 1877.Google Scholar
- Griffing JP. Body measurements of black-tailed jackrabbits of southeastern New Mexico with implications of Allen’s rule. J Mammal. 1974;55(3):674–8.View ArticlePubMedGoogle Scholar
- Lindsay SL. Geographic Size and Non-Size Variation in Rocky Mountain
*Tamiasciurus hudsonicus*: Significance in Relation to Allen’s Rule and Vicariant Biogeography. J Mammal. 1987;68(1):39–48.View ArticleGoogle Scholar - Ray C. The application of Bergmann’s and Allen’s rules to the poikilotherms. J Morphol. 1960;106(1):85–108.View ArticlePubMedGoogle Scholar
- Johnston RF, Selander RK. Evolution in the house sparrow. II. Adaptive differentiation in North American populations. Evolution. 1971;25:1–28.View ArticleGoogle Scholar
- Danner RM, Greenberg R. A critical season approach to Allen’s rule: bill size declines with winter temperature in a cold temperate environment. J Biogeogr. 2015;42(1):114–20.View ArticleGoogle Scholar
- Nudds RL, Oswald SA. An interspecific test of allen’s rule: evolutionary implications for endothermic species. Evolution. 2007;61(12):2839–48.View ArticlePubMedGoogle Scholar
- Cartar RV, Guy Morrison RI. Metabolic correlates of leg length in breeding arctic shorebirds: the cost of getting high. J Biogeogr. 2005;32(3):377–82.View ArticleGoogle Scholar
- Fabre P-H, Hautier L, Dimitrov D, P Douzery E. A glimpse on the pattern of rodent diversification: a phylogenetic approach. BMC Evol Biol. 2012;12(1):88.View ArticlePubMedPubMed CentralGoogle Scholar
- Bibi F. A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics. BMC Evol Biol. 2013;13(1):166.View ArticlePubMedPubMed CentralGoogle Scholar
- Agnarsson I, Zambrana-Torrelio CM, Flores-Saldana NP, May-Collado LJ. A time-calibrated species-level phylogeny of bats (Chiroptera, Mammalia). PLoS Cur. 2011;3.Google Scholar
- Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Grenyer R, et al. The delayed rise of present-day mammals. Nature. 2007;446(7135):507–12.View ArticlePubMedGoogle Scholar
- Jones KE, Bielby J, Cardillo M, Fritz SA, O’Dell J, Orme CDL, et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals: Ecological Archives E090-184. Ecology. 2009;90(9):2648.View ArticleGoogle Scholar
- Bartoszek K, Pienaar J, Mostad P, Andersson S, Hansen TF. A phylogenetic comparative method for studying multivariate adaptation. J Theor Biol. 2012;314:204–15.View ArticlePubMedGoogle Scholar
- Burnham KP, Anderson DR. Multimodel inference understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33(2):261–304.View ArticleGoogle Scholar
- R Development Core Team. R: a language and environment for statistical computing. In
*.*Vienna, Austria, http://www.r-project.org; 2014. - Grabowski M, Voje KL, Hansen TF. Evolutionary modeling and correcting for observation error support a 3/5 brain-body allometry for primates. J Hum Evol. 2016;94:106–16.View ArticlePubMedGoogle Scholar
- Barlow JC. Land mammals from Uruguay: ecology and zoogeography. 1965.Google Scholar
- Weir BJ. The tuco-tuco and plains viscacha. In: Symposia of the Zoological Society of London. 1974. p. 113–30.Google Scholar
- Hannon B, Ruth M. Bat Thermo-Regulation. In: Hannon B, Ruth M, editors. Modeling Dynamic Biological Systems. Cham: Springer International Publishing; 2014. p. 147–50.Google Scholar
- Mysterud A, Østbye E. Cover as a habitat element for temperate ungulates: effects on habitat selection and demography. Wildlife Soc B. 1999;27(2):385–94.Google Scholar
- Kluger MJ, Heath JE. Vasomotion in the bat wing: A thermoregulatory response to internal heating. Comp Biochem Physiol. 1970;32(2):219–26.View ArticlePubMedGoogle Scholar
- Marchetti K, Price T, Richman A. Correlates of wing morphology with foraging behaviour and migration distance in the genus Phylloscopus. J Avian Biol. 1995;26(3):177–81.View ArticleGoogle Scholar