Global patterns of body size evolution are driven by precipitation in legless amphibians

Body size shapes ecological interactions across and within species, ultimately influencing the evolution of large‐scale biodiversity patterns. Therefore, macroecological studies of body size provide a link between spatial variation in selection regimes and the evolution of animal assemblages through space. Multiple hypotheses have been formulated to explain the evolution of spatial gradients of animal body size, predominantly driven by thermal (Bergmann's rule), humidity (‘water conservation hypothesis’) and resource constraints (‘resource rule’, ‘seasonality rule’) on physiological homeostasis. However, while integrative tests of all four hypotheses combined are needed, the focus of such empirical efforts needs to move beyond the traditional endotherm–ectotherm dichotomy, to instead interrogate the role that variation in lifestyles within major lineages (e.g. classes) play in creating neglected scenarios of selection via analyses of largely overlooked environment–body size interactions. Here, we test all four rules above using a global database spanning 99% of modern species of an entire Order of legless, predominantly underground‐dwelling amphibians (Gymnophiona, or caecilians). We found a consistent effect of increasing precipitation (and resource abundance) on body size reductions (supporting the water conservation hypothesis), while Bergmann's, the seasonality and resource rules are rejected. We argue that subterranean lifestyles minimize the effects of aboveground selection agents, making humidity a dominant selection pressure – aridity promotes larger body sizes that reduce risk of evaporative dehydration, while smaller sizes occur in wetter environments where dehydration constraints are relaxed. We discuss the links between these principles with the physiological constraints that may have influenced the tropically‐restricted global radiation of caecilians.     

The geography of body size – challenges of the interspecific approach

Recent compilations of large-scale data bases on the geographical distributions and body sizes of animals, coupled with developments in spatial statistics, have led to renewed interest in the geographical distribution of animal body sizes and the interspecific version of Bergmann's rule. Standard practice seems to be an examination of mean body sizes within higher taxa on gridded maps, with little regard to species richness or phylogeny. However, because the frequency distribution of body sizes is typically highly skewed, average size within grid cells may differ significantly between species-rich and species-poor cells even when the median and modal sizes remain constant. Species richness influences body size patterns because species are not added to communities at random in relation to their size: areas of low diversity are characterized by a higher range of body sizes than is expected by chance. Finally, a consideration of phylogenetic structure within taxa is necessary to elucidate whether patterns in the geography of size result from turnover between or within intermediate taxonomic levels. We suggest that the highest and lowest quantiles of body size distribution be mapped in order to expose possible physiological or ecological limitations on body size.     

Latitudinal variation in the shape of the species body size distribution: an analysis using freshwater fishes

Many taxonomic and ecological assemblages of species exhibit a right-skewed body size-frequency distribution when characterized at a regional scale. Although this distribution has been frequently described, factors influencing geographic variation in the distribution are not well understood, nor are mechanisms responsible for distribution shape. In this study, variation in the species body size-frequency distributions of 344 regional communities of North American freshwater fishes is examined in relation to latitude, species richness, and taxonomic composition. Although the distribution of all species of North American fishes is right-skewed, a negative correlation exists between latitude and regional community size distribution skewness, with size distributions becoming left-skewed at high latitudes. This relationship is not an artifact of the confounding relationship between latitude and species richness in North American fishes. The negative correlation between latitude and regional community size distribution skewness is partially due to the geographic distribution of families of fishes and apparently enhanced by a nonrandom geographic distribution of species within families. These results are discussed in the context of previous explanations of factors responsible for the generation of species size-frequency distributions related to the fractal nature of the environment, energetics, and evolutionary patterns of body size in North American fishes.     

Climate history,human impacts and global body size of Carnivora (Mammalia: Eutheria) at multiple evolutionary scales

Aim One of the longest recognized patterns in macroecology, Bergmann’s rule, describes the tendency for homeothermic animals to have larger body sizes in cooler climates than their phylogenetic relatives in warmer climates. Here we provide an integrative process‐based explanation for Bergmann’s rule at the global scale for the mammal order Carnivora. Location Global. Methods Our database comprises the body sizes of 209 species of extant terrestrial Carnivora, which were analysed using phylogenetic autocorrelation and phylogenetic eigenvector regression. The interspecific variation in body size was partitioned into phylogenetic (P) and specific (S) components, and mean P‐ and S‐components across species were correlated with environmental variables and human occupation both globally and for regions glaciated or not during the last Ice Age. Results Three‐quarters of the variation in body size can be explained by phylogenetic relationships among species, and the geographical pattern of mean values of the P‐component is the opposite of the pattern predicted by Bergmann’s rule. Partial regression revealed that at least 43% of global variation in the mean phylogenetic component is explained by current environmental factors. In contrast, the mean S‐component of body size shows large positive deviations from ancestors across the Holarctic, and negative deviations in southern South America, the Sahara Desert, and tropical Asia. There is a moderately strong relationship between the human footprint and body size in glaciated regions, explaining 19% of the variance of the mean P‐component. The relationship with the human footprint and the P‐component is much weaker in the rest of the world, and there is no relationship between human footprint and S‐component in any region. Main conclusions Bergmannian clines are stronger at higher latitudes in the Northern Hemisphere because of the continuous alternation of glacial–interglacial cycles throughout the late Pliocene and Pleistocene, which generated increased species turnover, differential colonization and more intense adaptive processes soon after glaciated areas became exposed. Our analyses provide a unified explanation for an adaptive Bergmann’s rule within species and for an interspecific trend towards larger body sizes in assemblages resulting from historical changes in climate and contemporary human impacts.     

Natural selection on body size is mediated by multiple interacting factors: a comparison of beetle populations varying naturally and experimentally in body size

Body size varies considerably among species and among populations within species, exhibiting many repeatable patterns. However, which sources of selection generate geographic patterns, and which components of fitness mediate evolution of body size, are not well understood. For many animals, resource quality and intraspecific competition may mediate selection on body size producing large-scale geographic patterns. In two sequential experiments, we examine how variation in larval competition and resource quality (seed size) affects the fitness consequences of variation in body size in a scramble-competing seed-feeding beetle, Stator limbatus. Specifically, we compared fitness components among three natural populations of S. limbatus that vary in body size, and then among three lineages of beetles derived from a single base population artificially selected to vary in size, all reared on three sizes of seeds at variable larval density. The effects of larval competition and seed size on larval survival and development time were similar for larger versus smaller beetles. However, larger-bodied beetles suffered a greater reduction in adult body mass with decreasing seed size and increasing larval density; the relative advantage of being large decreased with decreasing seed size and increasing larval density. There were highly significant interactions between the effects of seed size and larval density on body size, and a significant three-way interaction (population-by-density-by-seed size), indicating that environmental effects on the fitness consequences of being large are nonadditive. Our study demonstrates how multiple ecological variables (resource availability and resource competition) interact to affect organismal fitness components, and that such interactions can mediate natural selection on body size. Studying individual factors influencing selection on body size may lead to misleading results given the potential for nonlinear interactions among selective agents.     

Interspecific patterns of species richness, geographic range size, and body size among New World venomous snakes

Many higher taxa exhibit latitudinal gradients in species richness, geographic range size, and body size. However, these variables are often interdependent, such that examinations of univariate or bivariate patterns alone may be misleading. Therefore, I examined latitudinal gradients in, and relationships between, species richness, geographic range size, and body size among 144 species of New World venomous snakes [families Elapidae (coral snakes) and Viperidae (pitvipers)]. Both lineages are monophyletic, collectively span 99° of latitude, and are extremely variable in body size and geographic range sizes. Coral snakes exhibit highest species richness near the equator, while pitviper species richness peaks in Central America. Species – range size distributions were strongly right-skewed for both families. There was little support for Bergmann's rule or Rapoport's rule for snakes of either family, as neither body size nor range size increased significantly with latitude. However, range area and median range latitude were positively correlated above 15° N, indicating a possible "Rapoport effect" at high northern latitudes. Geographic range size was positively associated with body size. Available continental area strongly influenced range size. Comparative (phylogenetically-based) analyses revealed that shared history is a poor predictor of range size variation within clades. Among vipers, trends in geographic range sizes may have been structured more by historical biogeography than by macroecological biotic factors.     

Global patterns of geographic range size in birds

Large-scale patterns of spatial variation in species geographic range size are central to many fundamental questions in macroecology and conservation biology. However, the global nature of these patterns has remained contentious, since previous studies have been geographically restricted and/or based on small taxonomic groups. Here, using a database on the breeding distributions of birds, we report the first (to our knowledge) global maps of variation in species range sizes for an entire taxonomic class. We show that range area does not follow a simple latitudinal pattern. Instead, the smallest range areas are attained on islands, in mountainous areas, and largely in the southern hemisphere. In contrast, bird species richness peaks around the equator, and towards higher latitudes. Despite these profoundly different latitudinal patterns, spatially explicit models reveal a weak tendency for areas with high species richness to house species with significantly smaller median range area. Taken together, these results show that for birds many spatial patterns in range size described in geographically restricted analyses do not reflect global rules. It remains to be discovered whether global patterns in geographic range size are best interpreted in terms of geographical variation in species assemblage packing, or in the rates of speciation, extinction, and dispersal that ultimately underlie biodiversity.     

Environmental factors, regional body size distributions and spatial variation in body size of local avian assemblages

Aim To determine how well variation in median body size of avian assemblages is predicted by (1) the environmental models usually employed in analyses of Bergmann's rule and (2) random sampling from the regional body size frequency distribution. If body size frequency distributions of local assemblages represent a random sample of a regional frequency distribution, then geographical variation in body sizes of assemblages might be a consequence of the determinants of spatial variation in species richness rather than direct influences on body size per se. Location Southern Africa. Methods Median body masses (as a measure of body size) of avian assemblages were calculated for quarter‐degree grid cells across South Africa and Lesotho. The relationship between median body mass and four environmental variables (minimum and maximum monthly temperatures, precipitation and seasonality in the normalized difference vegetation index, as a measure of seasonality in productivity) was examined using general linear models first without taking spatial autocorrelation into account, and then accounting for it by fitting an exponential spatial covariance structure. Model fit was assessed using the Akaike information criterion and Akaike weights. At each species richness value, random assemblages were sampled by either drawing species randomly from the regional body mass frequency distribution, or drawing species from the regional body mass frequency distribution with a probability proportional to their geographical distribution in the area. The ability of randomizations to predict actual body masses was examined using two‐tailed Fisher exact tests. Results Seasonality in productivity was the only environmental variable that remained a significant predictor of body mass variation in spatially explicit models, though the positive relationship was weak. When species richness was included in the models it remained the only significant predictor of size variation. Randomizations predicted median body mass poorly at low species richness, but well at high richness. Main conclusions Environmental models that have previously been proposed explain little of the variation in body mass across avian assemblages in South Africa. However, much of the variation in the median mass of assemblages could be predicted by randomly drawing species from the regional body mass frequency distribution, particularly using randomizations in which all species were drawn from the regional body mass frequency distribution with equal probability and at high species richness values. This outcome emphasizes the need to consider null expectations in investigations of the geographical variation in body size together with the probable environmental mechanisms underlying spatial variation in average size. Moreover, it suggests that in the South African avifauna, spatial variation in the body sizes of assemblages may be determined indirectly by the factors that influence geographical variation in species richness.     

Body size–climate relationships of European spiders

Aim Geographic body size patterns of mammals and birds can be partly understood under the framework of Bergmann’s rule. Climatic influences on body size of invertebrates, however, appear highly variable and lack a comparable, generally applicable theoretical framework. We derived predictions for body size–climate relationships for spiders from the literature and tested them using three datasets of variable spatial extent and grain. Location Europe. Methods To distinguish climate from space, we compared clines in body size within three datasets with different degrees of co‐variation between latitude and climate. These datasets were: (1) regional spider faunas from 40 European countries and large islands; (2) local spider assemblages from standardized samples in 32 habitats across Europe; and (3) local spider assemblages from Central European habitats. In the latter dataset climatic conditions were determined more by habitat type than by geographic position, and therefore this dataset provided a non‐spatial gradient of various microclimates. Spider body size was studied in relation to latitude, temperature and water availability. Results In all three datasets the mean body size of spider assemblages increased from cool/moist to warm/dry environments. This increase could be accounted for by turnover from small‐bodied to large‐bodied spider families. Body size–climate relationships within families were inconsistent. Main conclusions Starvation resistance and accelerated maturation can be ruled out as explanations for the body size clines recorded, because they predict the inverse of the observed relationship between spider body size and temperature. The relationship between body size and climate was partly independent of geographic position. Thus, the restriction of large‐bodied spiders to their glacial refugia owing to dispersal limitations can be excluded. Our results are consistent with mechanisms invoking metabolic rate, desiccation resistance and community interactions to predict a decrease in body size from warm and dry to cool and moist conditions.     

Interspecific pairwise relationships among body size, clutch size and latitude: deconstructing a macroecological triangle in birds

Aim Ecogeographical ‘rules’, large‐scale patterns in ecological variables across geographical space, can provide important insights into the mechanisms of evolution and ecological assembly. However, interactions between rules could obscure both the observation of large‐scale patterns and their interpretation. Here, we examine a system of three variables interrelated by ecogeographical rules – the latitudinal increase in body size within closely related homeotherms (Bergmann’s rule), the negative allometry of clutch size (Calder’s rule) and the latitudinal increase in clutch size (Lack’s rule) – in a global dataset of birds. Location Global. Methods We used linear regressions and meta‐analysis techniques to quantify the three rules across clades and through the taxonomic hierarchy. Path analysis was used to quantify interactions between rules at multiple taxonomic levels, as a function of both phylogenetic inheritance of traits and indirect feedbacks between the three rules. Independent contrasts analyses were performed on four clades with available phylogenies, and the taxonomic partitioning of variation in each trait was quantified. Results Standardizing across all clades, Lack’s and Bergmann’s rules were supported at all taxonomic levels, with Calder’s rule being supported at the order level. Lack’s rule was consistently stronger and more often detected than the other two rules. Path analysis showed that the indirect effects often outweighed the direct effects of Calder’s rule at the genus level and Bergmann’s rule at the order level. Strong interactions between Calder’s and Bergmann’s rules led to a trade‐off between the rules depending on taxonomic resolution. Main conclusions We found strong interactions between Bergmann’s, Lack’s and Calder’s rules in birds, and these interactions varied in strength and direction over the taxonomic hierarchy and among avian clades. Ecogeographical rules may be masked by feedbacks from other, correlated variables, even when the underlying selective mechanism is operating. The apparently conflicting pairwise relationships among clutch size, body size and latitude illustrate the difficulty of interpreting individual pairwise correlations without recognition of interdependence with other variables.     

Resources and global avian assemblage structure in forests

Explaining spatial variation in a number of bird species, particularly from temperate to tropical regions, has been a longstanding challenge. We test at a global scale whether species‐rich forest assemblages are associated with division of a larger resource pool, a finer division of that pool, or some combination of the two. Species richness increases with increasing assemblage abundance, biomass and energy use. As assemblage abundance, biomass and energy use increase with increasing energy availability, and as per species numbers of individuals, biomass and energy use do not decrease with increasing energy availability, we provide direct evidence that the avian species–energy relationship in forests is associated foremost with an increase in the size of the resource pool and not with a finer level of its subdivision.     

Intra-annual patterns of saproxylic beetle assemblages inhabiting Mediterranean oak forests

The temporal factor has been widely neglected in studies related to species distribution. However, knowing the species phenology throughout the year could be a key element in ecological studies, more so considering the current climate change framework. The intra-annual temporal patterns of the saproxylic beetle assemblage inhabiting the oak forests in the western Iberian Peninsula were evaluated in terms of species richness and composition. Moreover, we analysed the temporal patterns of species richness and body size, with respect to the abiotic factors of temperature and relative humidity. The marked seasonality of the Mediterranean weather governed, to a great extent, the temporal patterns of the assemblage, as both species richness and body size showed a positive correlation with temperature and a negative correlation with relative humidity. The high temperatures reached in August limited the activity of certain species, causing a marked decrease in species richness. The results revealed that species body size decreased in the colder months of the year and vice versa, which could fit with the inverse of Bergmann’s rule; nevertheless, studies are needed to address the contribution of species physiology and autoecology to the observed temporal patterns. Regarding species composition, the partition of beta diversity in its two components showed a constant process of species turnover throughout the year. We conclude that, due to global warming, the study of temporal patterns of assemblages is becoming all the more important, since it could significantly affect the spatio-temporal distribution of species and the interactions between them.     

Genome size variation and species diversity in salamanders

Salamanders (Urodela) have among the largest vertebrate genomes, ranging in size from 10 to 120 pg. Although changes in genome size often occur randomly and in the absence of selection pressure, nonrandom patterns of genome size variation are evident among specific vertebrate lineages. Several reports suggest a relationship between species richness and genome size, but the exact nature of that relationship remains unclear both within and across different taxonomic groups. Here, we report (a) a negative relationship between haploid genome size (C‐value) and species richness at the family taxonomic level in salamander clades; (b) a correlation of C‐value and species richness with clade crown age but not with diversification rates; (c) strong associations between C‐value and both geographic area and climatic‐niche rate. Finally, we report a relationship between C‐value diversity and species diversity at both the family‐ and genus‐level clades in urodeles.     

Spatial turnover in the global avifauna

Despite its wide implications for many ecological issues, the global pattern of spatial turnover in the occurrence of species has been little studied, unlike the global pattern of species richness. Here, using a database on the breeding distributions of birds, we present the first global maps of variation in spatial turnover for an entire taxonomic class, a pattern that has to date remained largely a matter of conjecture, based on theoretical expectations and extrapolation of inconsistent patterns from different biogeographic realms. We use these maps to test four predictions from niche theory as to the form that this variation should take, namely that turnover should increase with species richness, towards lower latitudes, and with the steepness of environmental gradients and that variation in turnover is determined principally by rare (restricted) species. Contrary to prediction, we show that turnover is high both in areas of extremely low and high species richness, does not increase strongly towards the tropics, and is related both to average environmental conditions and spatial variation in those conditions. These results are closely associated with a further important and novel finding, namely that global patterns of spatial turnover are driven principally by widespread species rather than the restricted ones. This complements recent demonstrations that spatial patterns of species richness are also driven principally by widespread species, and thus provides an important contribution towards a unified model of how terrestrial biodiversity varies both within and between the Earth's major land masses.     

Regional variation in the historical components of global avian species richness

Aim Using a global data base of the distribution of extant bird species, we examine the evidence for spatial variation in the evolutionary origins of contemporary avian diversity. In particular, we assess the possible role of the timing of mountain uplift in promoting diversification in different regions. Location Global. Methods We mapped the distribution of avian richness at four taxonomic levels on an equal‐area 1° grid. We examined the relationships between richness at successive taxonomic levels (e.g. species richness vs. genus richness). We mapped the residuals from linear regressions of these relationships to identify areas that are exceptional in the number of lower taxa relative to the number of higher taxa. We use generalized least squares models to test the influence of elevation range and temperature on lower‐taxon richness relative to higher‐taxon richness. Results Peaks of species richness in the Neotropics were congruent with patterns of generic richness, whilst peaks in Australia and the Himalayas were congruent with patterns of both genus and family richness. Hotspots in the Afrotropics did not reflect higher‐taxon patterns. Regional differences in the relationship between richness at successive taxonomic levels revealed variation in patterns of taxon co‐occurrence. Species and genus co‐occurrence was positively associated with elevational range across much of the world. Taxon occurrence in the Neotropics was associated with a positive interaction between elevational range and temperature. Conclusions These results demonstrate that contemporary patterns of richness show different associations with higher‐taxon richness in different regions, which implies that the timing of historical effects on these contemporary patterns varies across regions. We suggest that this is due to dispersal limitation and phylogenetic constraints on physiological tolerance limits promoting diversification. We speculate that diversification rates respond to long‐term changes in the Earth's topography, and that the role of tropical mountain ranges is implicated as a correlate of contemporary diversity, and a source of diversification across avian evolutionary history.     

A phylogenetic perspective on the evolutionary processes of floristic assemblages within a biodiversity hotspot in eastern Asia

Abstract How to maximize the conservation of biodiversity is critical for conservation planning, particularly given rapid habitat loss and global climatic change. The importance of preserving phylogenetic diversity has gained recognition due to its ability to identify some influences of evolutionary history on contemporary patterns of species assemblages that traditional taxonomic richness measures cannot identify. In this study, we evaluate the relationship between taxonomic richness and phylogenetic diversity of angiosperms at genus and species levels and explore the spatial pattern of the residuals of this relationship. We then incorporate data on historical biogeography to understand the process that shaped contemporary floristic assemblages in a global biodiversity hotspot, Yunnan Province, located in southwestern China. We identified a strong correlation between phylogenetic diversity residuals and the biogeographic affinity of the lineages in the extant Yunnan angiosperm flora. Phylogenetic diversity is well correlated with taxonomic richness at both genus and species levels between floras in Yunnan, where two diversity centers of phylogenetic diversity were identified (the northwestern center and the southern center). The northwestern center, with lower phylogenetic diversity than expected based on taxonomic richness, is rich in temperate‐affinity lineages and signifies an area of rapid speciation. The southern center, with higher phylogenetic diversity than predicted by taxonomic richness, contains a higher proportion of lineages with tropical affinity and seems to have experienced high immigration rates. Our results highlight that maximizing phylogenetic diversity with historical interpretation can provide valuable insights into the floristic assemblage of a region and better‐informed decisions can be made to ensure different stages of a region's evolutionary history are preserved.     

An Intercontinental Analysis of Climate-Driven Body Size Clines in Reptiles: No Support for Patterns, No Signals of Processes

Climatic gradients impose clinal selection on animal ecological and physiological performance, often promoting geographic body size clines. Bergmann’s rule predicts that body size increases with decreasing environmental temperatures given the need to retain body-heat through adjustments of body-mass-to-surface-area ratio. This prediction generally holds for endotherms, but remains controversial for ectotherms. An alternative interpretation, the ‘resource rule’, suggests that food abundance, primary productivity and precipitation (which, unlike temperature, do not necessarily correlate with geography), drive body size clines. We investigate geographic variation in body size within 65 species of lizards and snakes (squamates) based on an intercontinental dataset (6,500+ specimens belonging to 56 Israeli species, and multiple populations of nine Liolaemus species from Argentina and Chile). Bergmann’s rule is only rarely supported by our data (in four species, 6 %), whereas six species (9 %) follow its converse (hence, it is unsupported in 94 % of cases). Similarly, size increases with resource abundance in only 12 species (18 %). Therefore, although neither of the rules is supported, factors suggested by the resource rule are better predictors of body size than temperature. Surprisingly, we show that some measures of the extent of a species’ climatic envelope do not affect the likelihood of it showing a size-climate relationship. We conclude that negative size-temperature associations are an exception rather than a generality among squamates.     

Abundance, diversity and body size: patterns from a grassland arthropod community

1. The empirical relationships among body size, species richness and number of individuals may give insight into the factors controlling species diversity and the relative abundances of species. To determine these relationships, we sampled the arthropods of grasslands and savannahs at Cedar Creek, MN using sweep nets (90 525 individuals of 1225 species) and pitfall traps (12 721 individuals of 92 species). Specimens were identified, enumerated and measured to determine body size.
2. Both overall and within abundant taxonomic orders, species richness and numbers of individuals peaked at body sizes intermediate for each group. Evolution could create unimodal diversity patterns by random diversification around an ancestral body size or from size-dependent fitness differences. Local processes such as competition or predation could also create unimodal diversity distributions.
3. The average body size of a species depended significantly on its taxonomic order, but on contemporary trophic role only within the context of taxonomic order.
4. Species richness ( S _i) within size classes was related to the number of individuals ( I _i) as S _i =  I i^0·5. This relationship held across a 100 000-fold range of body sizes. Within size classes, abundance distributions of size classes were all similar power functions. A general rule of resource division, together with similar minimum population sizes, is sufficient to generate the relationship between species richness and number of individuals.
5. Smaller bodied species had slightly shallower abundance distributions and may, in general, persist at lower densities than larger species.
6. Our results suggest there may be fewer undescribed small arthropod species than previously thought and that most undescribed species will be smaller than arthropods.     

Shallow water ray-finned marine fishes follow Bergmann’s rule

Understanding the interspecific variation in body size across macroclimatic gradients has been of paramount importance to naturalists and biogeographers. Bergmann’s rule, which describes a trend of increasing body size polewards, is arguably the best-known ecogeographical rule in terrestrial environments but remains largely unexplored in the marine realm. In this study we tested Bergmann's rule in marine ray-finned fishes (Pisces, Actinopterygii), analyzing the relationship between body size and latitude in 5662 species. To examine possible underlying mechanisms, we adopted a cross-species approach to evaluate the association of body size with four predictors: Sea Surface Temperature, Net Primary Productivity, Salinity, and Human impact. We analyzed the relationships between body size and environmental and anthropogenic variables building mixed linear models, which considered the taxonomic structure in the data. We conducted complementary analyses dividing the data into five latitudinal bands. Actinopterygii showed a clear Bergmannian pattern, with the largest species observed in temperate regions, being the first global analysis on ray-finned fishes showing a pattern consistent with Bergmann’s rule. Sea Surface Temperature and Net Primary Productivity were the best predictors, in accordance with the time to sexual maturity and resource availability hypotheses. Our analyses based on latitudinal bands showed a differential response of body size to the environment, with temperature, salinity and human impact more strongly associated with size variation at cold environments. These results agree with previous studies on Bergmann’s rule for terrestrial ectothermic, freshwater and marine fishes. Our findings suggest that temperature rise in the ocean and growing human impact may have effects on the distribution of body size, thus altering ecosystem functioning. Fundamental differences often assumed to exist between marine and terrestrial systems are not so evidently reflected in the emergence of large-scale body size gradients.