Life history patterns in female moose (Alces alces): the relationship between age,body size,fecundity and environmental conditions

I examined the relationship between age, body size and fecundity in 833 female moose (Alces alces) from 14 populations in Sweden sampled during 1989–1992. Data on population density, food availability and climatic conditions were also collected for each population. Age and body mass were both significantly positively related to fecundity, measured as ovulation rate, among female moose. The relationship between the probability of ovulation and body mass was dependent on age with (1) a higher body mass needed in younger females for attaining a given fecundity, and (2) body mass having a stronger effect on fecundity in yearling (1.5 year) than in older (2.5 year) females. Thus, a 40 kg increase in yearling body mass resulted in a 42% increase in the probability of ovulation as compared to a 6% increase in older females. The lower reproductive effort per unit body mass, and the relatively stronger association between fecundity and body mass in young female moose compared to older ones, is likely to primarily represent a mechanism that trades off early maturation against further growth, indicating a higher cost of reproduction in young animals. In addition to age and body mass, population identity explained a significant amount of the individual variation in fecundity, showing that the relationship between body mass and fecundity was variable among populations. This variation was in turn related to the environment, in terms of climatic conditions forcing female moose living in relatively harsh/more seasonal climatic conditions to attain a 22% higher body mass to achive the same probability of multiple ovulation (twinning) as females living in climatically milder/less seasonal environments. The results suggests that the lower fecundity per unit body mass in female moose living in climatically harsh/more seasonal environments may be an adaptive response to lower rates of juvenile survival, compared to females experiencing relatively milder/less seasonal climatic conditions.     

Moose body mass variation revisited: disentangling effects of environmental conditions and genetics

Large-scale geographical variation in phenotypic traits within species is often correlated to local environmental conditions and population density. Such phenotypic variation has recently been shown to also be influenced by genetic structuring of populations. In ungulates, large-scale geographical variation in phenotypic traits, such as body mass, has been related to environmental conditions and population density, but little is known about the genetic influences. Research on the genetic structure of moose suggests two distinct genetic lineages in Norway, structured along a north-south gradient. This corresponds with many environmental gradients, thus genetic structuring provides an additional factor affecting geographical phenotypic variation in Norwegian moose. We investigated if genetic structure explained geographical variation in body mass in Norwegian moose while accounting for environmental conditions, age and sex, and if it captured some of the variance in body mass that previously was attributed to environmental factors. Genetic structuring of moose was the most important variable in explaining the geographic variation in body mass within age and sex classes. Several environmental variables also had strong explanatory power, related to habitat diversity, environmental seasonality and winter harshness. The results suggest that environmental conditions, landscape characteristics, and genetic structure should be evaluated together when explaining large-scale patterns in phenotypic characters or life history traits. However, to better understand the role of genetic and environmental effects on phenotypic traits in moose, an extended individual-based study of variation in fitness-related characters is needed, preferably in an area of convergence between different genetic lineages.     

The roles of pre‐ and post‐hatching growth rates in generating a latitudinal cline of body size in the eastern fence lizard (Sceloporus undulatus)

Countergradient variation in norms of reaction can dampen the direct effects of environmental influences on phenotypic traits, allowing phenotypic similarity among populations despite exposure to different environmental conditions. Such norms of reaction may occur at any phase of the life‐history (e.g. growth rates during both embryonic and postembryonic stages may influence geographical variation in adult body size). We collected gravid female lizards (Sceloporus undulatus) from northern (Indiana), central (Mississippi), and southern (Florida) populations, spanning almost the full latitudinal range of the species. Adult females from the southern population were smaller. Intrinsic growth rates of hatchlings were higher for the central population than for the other two populations. This pattern does not parallel the countergradient variation previously found in embryonic developmental rates among these populations. Earlier hatching enhanced survival rates of juveniles to a similar degree among populations, although juvenile survival rates in the field generally increase with latitude in this species. Our data reveal geographical variation in the ways in which intrinsic developmental/growth rates and survival shift during ontogeny, and suggest that latitudinal patterns in adult body size (such as Bergmann's rule) can result from both faster growth, and longer periods of growth. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 202–209.     

Multiple selection pressures generate adherence to Bergmann's rule in a Neotropical migratory songbird

Aim Bergmann's rule, the tendency for body size to be positively correlated with latitude, is widely accepted but the mechanisms behind the patterns are still debated. Bergmann's originally conceived mechanism was based on heat conservation; other proposed mechanisms invoke phylogeny, migration distance and resource seasonality. With the goal of examining these mechanisms, we quantified morphological variation across the breeding range of a Neotropical migratory songbird, the cerulean warbler (Dendroica cerulea). Location Deciduous forests of eastern North America. Methods We sampled nine cerulean warbler populations, spanning the species’ breeding range. We captured 156 males using targeted playback and model presentation, and included 127 adult males in our analyses of morphological variation. We used an information‐theoretical approach to identify climatic variables associated with geographical variation in body size. Results Cerulean warbler body size adheres to Bergmann's rule: individuals in northern populations are larger than those in southern populations. Variation in body size is best explained by variation in dry and wet‐bulb temperature and actual evapotranspiration. Main conclusions Adherence to Bergmann's rule by the cerulean warbler appears to be linked to thermodynamics (heat conservation in the north, evaporative cooling in the south) and resource seasonality. Multiple selection pressures can interact to generate a single axis of morphological geographical variation, and even subtle fluctuations in climatic variables can exert significant selection pressures. We suggest that the influence of selection pressures on migrants might be enhanced by migratory connectivity, providing further support for the important role played by this phenomenon in the ecology, evolution and population dynamics of migratory songbirds.     

Environmental phenology and geographical gradients in moose body mass

Intraspecific body mass in ungulates has often been shown to increase with latitude. The biological basis for such latitudinal gradients is, however, poorly known. Here we examined whether satellite-derived indices of environmental phenology, based on the normalised difference vegetation index (NDVI), as well as variables derived from meteorological stations, altitude, and population density, can explain latitudinal gradients and regional variation in body mass of Norwegian moose. The best model gave a considerably better fit than latitude alone, and included all explanatory environmental variables. Accordingly, heavy moose were found in areas with short and intense summers that were followed by long, cold winters, at low altitude relative to the tree-limit, and with low population density relative to the available plant biomass. This relationship was stronger for yearlings than for calves, except for the effect of population density. This indicates that differences in the characteristics of the vegetation quality and environmental phenology, as well as winter harshness and population density, are important factors that shape both the latitudinal and other geographical gradients in moose body mass.     

Behavioural Responses to Thermal Conditions Affect Seasonal Mass Change in a Heat-Sensitive Northern Ungulate

Background

Empirical tests that link temperature-mediated changes in behaviour (activity and resource selection) to individual fitness or condition are currently lacking for endotherms yet may be critical to understanding the effect of climate change on population dynamics. Moose (Alces alces) are thought to suffer from heat stress in all seasons so provide a good biological model to test whether exposure to non-optimal ambient temperatures influence seasonal changes in body mass. Seasonal mass change is an important fitness correlate of large herbivores and affects reproductive success of female moose.

Methodology/Principal Findings

Using GPS-collared adult female moose from two populations in southern Norway we quantified individual differences in seasonal activity budget and resource selection patterns as a function of seasonal temperatures thought to induce heat stress in moose. Individual body mass was recorded in early and late winter, and autumn to calculate seasonal mass changes (n = 52 over winter, n = 47 over summer). We found large individual differences in temperature-dependent resource selection patterns as well as within and between season variability in thermoregulatory strategies. As expected, individuals using an optimal strategy, selecting young successional forest (foraging habitat) at low ambient temperatures and mature coniferous forest (thermal shelter) during thermally stressful conditions, lost less mass in winter and gained more mass in summer.

Conclusions/Significance

This study provides evidence that behavioural responses to temperature have important consequences for seasonal mass change in moose living in the south of their distribution in Norway, and may be a contributing factor to recently observed declines in moose demographic performance. Although the mechanisms that underlie the observed temperature mediated habitat-fitness relationship remain to be tested, physiological state and individual variation in thermal tolerance are likely contributory factors. Climate-related effects on animal behaviour, and subsequently fitness, are expected to intensify as global warming continues.     

Body size variation and sexual size dimorphism across climatic gradients in the widespread treefrog Scinax fuscovarius (Anura,Hylidae)

Variation in body size represents one of the crucial raw materials for evolution. However, at present, it is still being debated what is the main factor affecting body size or if the final body size is the consequence of several factors acting synergistically. To evaluate this, widespread species seem to be suitable models because the different populations occur along a geographical gradient and under contrasted climatic and environmental conditions. Here we describe the spatial pattern of variation in body size and sexual size dimorphism in the snouted treefrog Scinax fuscovarius (Anura, Hylidae) along a 10° range in latitude, 25° longitude, and 2000 m in altitude from Argentina, Brazil and Paraguay using an information‐theoretic approach to evaluate the support of the data for eight a priori hypotheses proposed in the literature to account for geographical body size, and three hypotheses for sexual size dimorphism variation. Body size of S. fuscovarius varied most dramatically with longitude and less so with latitude; frogs were largest in the northwestern populations. Body size was positively related with precipitation seasonality, and negatively with annual precipitation. Furthermore, the degree of sexual size dimorphism was greatest in the western populations with less annual precipitation, as the increase in body size was stronger for females. Our results on body size variation are consistent with two ecogeographical hypotheses, the starvation resistance and the water availability hypotheses, while our results on sexual size dimorphism in S. fuscovarius supports the differential‐plasticity hypothesis but the inverse to Rensch's rule and the parental investment hypothesis. Due to the weak association between environmental variables and body size and sexual size dimorphism variation, we stress that there are other factors, mainly those related to the life history, driving the geographical variation of S. fuscovarius.     

Variation in juvenile growth rates among and within latitudinal populations of the medaka

In ectotherms, lower temperatures at high latitudes would theoretically reduce annual growth rates of individuals. If slower growth and resulting smaller body size reduce fitness, individuals at high latitudes may evolve compensatory growth. This study compares individual growth rates among and within 12 latitudinal populations of the medaka (Oryzias latipes). Growth rates during juvenile stage were measured in a common, temperature-controlled (28°C) environment. The results revealed that juvenile growth rates differed significantly among the populations. Growth rates were, moreover, significantly correlated with latitudes of source populations, such that higher-latitude individuals grew faster. Significant variation in growth rates among full-sib families within populations was also demonstrated. The results strongly suggest that higher-latitude O. latipes have acquired a greater capacity for growth as an adaptation to shorter growing seasons (which would reduce annual growth rates), thus refuting probability processes, i.e., genetic drift, founder, or bottleneck effects, as a cause of the among-population variation.     

Patterns and causes of demographic variation in a harvested moose population: evidence for the effects of climate and density-dependent drivers

1.?Better understanding of the mechanisms affecting demographic variation in ungulate populations is needed to support sustainable management of harvested populations. While studies of moose Alces alces L. populations have previously explored temporal variation in demographic processes, managers responsible for populations that span large heterogeneous landscapes would benefit from an understanding of how demography varies across biogeographical gradients in climate and other population drivers. Evidence of thresholds in population response to manageable and un-manageable drivers could aid resource managers in identifying limits to the magnitude of sustainable change. 2.?Generalized additive models (GAMs) were used to evaluate the relative importance of population density, habitat abundance, summer and winter climatic conditions, primary production, and harvest intensity in explaining spatial variation in moose vital rates in Ontario, Canada. Tree regression was used to test for thresholds in the magnitudes of environmental predictor variables that significantly affected population vital rates. 3.?Moose population growth rate was negatively related to moose density and positively related to the abundance of mixed deciduous habitat abundant in forage. Calf recruitment was negatively related to a later start of the growing season and calf harvest. The ratio of bulls to cows was related to male harvest and hunter access, and thresholds were evident in predictor variables for all vital rate models. 4.?Findings indicate that the contributions of density-dependent and independent factors can vary depending on the scale of population process. The importance of density dependence and habitat supply to low-density ungulate populations was evident, and management strategies for ungulates may be improved by explicitly linking forest management and harvest. Findings emphasize the importance of considering summer climatic influences to ungulate populations, as recruitment in moose was more sensitive to the timing of vegetation green-up than winter severity. The efficacy of management decisions for harvested ungulates may require regional shifts in targets where populations span bioclimatic gradients. The use of GAMs in combination with recursive partitioning was demonstrated to be an informative analytical framework that captured nonlinear relationships common in natural processes and thresholds that are relevant to population management in diverse systems.     

Structure of allozyme variation in Nordic Silene nutans (Caryophyllaceae): population size, geographical position and immigration history

We investigated allozyme variation in 34 populations of the perennial herb Silene nutans from Sweden and northern Finland, areas that were ice-covered during the last (Weichselian) glaciation. The present geographical structure of genetic variation in S. nutans in Sweden and northern Finland appears to have been mainly shaped by ancient historical processes. Patterns of variation in allele frequencies suggest two major postglacial immigration routes into Sweden, with populations entering the area from both the south and the east and forming a contact zone with admixed populations in central Sweden. While estimates of within-population genetic diversity and allelic richness are significantly correlated with present population size and geographical position (latitude), population size is not correlated with latitude. Low genetic diversity in the northern populations is more likely to have resulted from ancient stochastic events during the process of immigration than from recent population fragmentation. F _IS values are high and increase with latitude. Evidence of recent bottlenecks was detected in several southern Swedish populations: these can be interpreted in terms of population fragmentation as a result of anthropogenic disturbance. Soil pH is uncorrelated with population size and position.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 81 , 357–371.     

Body size clines in the European badger and the abundant centre hypothesis

Aim To test the abundant centre hypothesis by analysing the physical and climatic factors that influence body size variation in the European badger (Meles meles). Location Data were compiled from 35 locations across Europe. Methods We used body mass, body length and condylo‐basal length (CBL) as surrogates of size. We also compiled data on latitude, several climatic variables, habitat type and site position relative to the range edge. We collapsed all continuous climatic variables into independent vectors using principal components analysis (PCA), and used a general linear model to explain the morphometric variation in badger populations across the species’ range. Results Body mass and body length were nonlinearly and significantly related to latitude. In contrast, CBL was linearly related to latitude. Body mass changed nonlinearly along the temperature (PC1) gradient, with the highest values observed at mid‐range. Furthermore, body mass, body length and CBL differed significantly among habitats, with badgers showing larger size in temperate habitats and core areas relative to peripheral zones. Main conclusions Our analysis supports the nonlinear pattern predicted by the abundant centre hypothesis only for body mass and body length. These results imply that individuals are largest and heaviest at the centre of the climatic range of badger distribution. Variation of CBL with latitude follows a linear trend, consistent with Bergmann’s rule. Our results provide mixed support for the abundant centre hypothesis, and suggest food availability/quality to be the main mechanism underlying body size clines in this species.     

Population characteristics predict responses in moose body mass to temporal variation in the environment

1. A general problem in population ecology is to predict under which conditions stochastic variation in the environment has the stronger effect on ecological processes. By analysing temporal variation in a fitness-related trait, body mass, in 21 Norwegian moose Alces alces (L.) populations, we examined whether the influence of temporal variation in different environmental variables were related to different parameters that were assumed to reflect important characteristics of the fundamental niche space of the moose. 2. Body mass during autumn was positively related to early access to fresh vegetation in spring, and to variables reflecting slow phenological development (low June temperature, a long spring with a slow plant progression during spring). In contrast, variables related to food quantity and winter conditions had only a minor influence on temporal variation in body mass. 3. The magnitude of the effects of environmental variation on body mass was larger in populations with small mean body mass or living at higher densities than in populations with large-sized individuals or living at lower densities. 4. These results indicate that the strongest influence of environmental stochasticity on moose body mass occurs towards the borders of the fundamental niche space, and suggests that populations living under good environmental conditions are partly buffered against fluctuations in environmental conditions.     

Ecogeographic analysis of morphological and life-history variation in the Italian treefrog

We describe the spatial pattern of variation in body size, body shape, and adult population age structure of the Italian treefrog (Hyla intermedia). By means of a non-parametric test of matrix association (the Partial Mantel Test), we investigate the role of climate and geography in explaining the patterns of phenotypic variation. Body size is the most important source of morphometric variation, both at intra- and inter-populational levels. Body shape, despite its secondary role in explaining variation, still maintains high discriminative power among populations. Age structure differs significantly among populations, since body size correlates positively to age, these differences partly explain geographic variation of body size. The patterns of variation in both age and body size are congruent with climate variation, suggesting causal relationship. However, we find no evidence for adaptive explanations and suggest that climate-induced demographically differential survivorship might be the proximate cause responsible for the associations observed. Unlike body size, body shape varies congruently with the pattern of between-population geographic distances, which we show to be highly correlated with their genetic distances, thus suggesting isolation by distance as the proximate cause of the association.     

Changes in cuticular lipids, water loss and desiccation resistance in a tropical drosophilid: analysis of variation between and within populations

We investigated within as well as between population variability in desiccation resistance, cuticular lipid mass per fly and cuticular water loss in nine geographical populations of a tropical drosophilid, Zaprionus indianus. Interestingly, the amount of cuticular lipids and desiccation resistance in this non-melanic species are significantly higher as compared with melanic Drosophila melanogaster. On the basis of isofemale line analysis, within population trait variability in cuticular lipid mass per fly is positively correlated with desiccation resistance and negatively correlated with cuticular water loss but show lack of correlation with body size. We observed geographical variation in the amount of cuticular lipid mass per fly in Zaprionus indianus but no such divergence was found in D.melanogaster. In both the species, geographical variations in desiccation resistance are negatively correlated with cuticular water loss but the underlying mechanisms for changes in cuticular permeability are quite different. Thus, we may suggest that body melanisation and cuticular lipids may represent alternative strategies for coping with dehydration stress in melanic versus non-melanic drosophilids. For both the species, desiccation resistance and cuticular water loss are correlated with regular increase in aridity in the northern subtropical localities as compared with southern peninsular humid tropical localities. The role of climatic selection is evident from multiple regression analysis with seasonal changes in temperature and humidity (Tcv and RHcv) of the sites of origin of populations of Zaprionus indianus along latitude.     

Latitudinal variation of morphological characteristics in the swallowtail Sericinus montelus Gray, 1798 (Lepidoptera: Papilionidae)

Patterns of latitudinal variation in the phenotype or genotype of an organism may provide evidence for natural selection. In this study, we investigated seven populations of swallowtail Sericinus montelus Gray, 1798 (Lepidoptera: Papilionidae), a non‐migratory species, to explore the latitudinal variation of morphological characteristics in adults. The results showed that body size and the development of dark pigmentation on wings in this species responded strongly to latitude. The body size of both male and female adult of S. montelus was negatively correlated with latitude. These findings provided solid evidence to support the converse Bergmann's rule. We considered that the observed variation in morphological characteristics was most likely mediated by the seasonal length and thermoperiod to adapt to different latitudinal environment (e.g. shortened developmental time of immature stages for smaller body size at higher latitude). Moreover, the tendency towards progressively darker colour patterns (only in adult males) at increasingly low latitudes was consistent with Gloger's rule. We suggested that the observed colour variation was most likely associated with thermoregulation. Slight variation in the morphology of the W‐shaped stripe on the forewing of adult females was also found, and we presumed that the functions of sexual preferences, mimicry and thermoregulation might be involved.     

Genetic variation in relation to demography of peripheral pool frog populations (Rana lessonae)

Summary Genetic variation in an isolated northern metapopulation of the pool frog (Rana lessonae) in Sweden was compared to that of Central European populations using enzyme electrophoresis and literature data. Of the 31 loci scored, two (EST-2 andIDH-2) were polymorphic while no variation occurred in seven of the eight loci which are polymorphic in Central European populations.The heterozygosity level of the Swedish pool frogs is very low compared to that of other anuran populations, but their mean proportion of fertilized eggs within egg masses (97.5%) was not lower than in more heterozygous species, and their body size-specific fecundity did not differ from that of Polish conspecifics. The low genetic variability of the Swedish pool frogs is discussed in relation to features of the local populations such as size (N), calculated effective size (N _e ) reproductive success and probable history. It is concluded that long-term strong fluctuations in population size caused by reporductive failure in cold years have contributed more to the low genetic variability than could a single founder event due to a recent introduction by man.     

Evidence of divergent growth rates among populations of the lizard Anolis carolinensis based on experimental manipulations of egg size

Geographic variation in body size is of special interest because it affects nearly all aspects of an organism’s life. I examined whether differences in body size among four populations of the green anole lizard, Anolis carolinensis, were attributable to maternal investment in egg size and/or growth rates of embryos and juveniles. Larger body size and larger egg size relative to female size in the northern part of the range have been documented in this species, and suggested to be adaptive responses to more extreme winters. The current study confirmed the trends in adult size and egg size in the north, but rejected the trend of larger egg size relative to body size in the south. To control for differences in maternal investment in egg size among populations, I performed yolk removals on eggs from two northern populations to produce comparably sized eggs relative to one southern population. This manipulation was designed to minimize the confounding effect of maternal investment in yolk, the primary energy reserves for eggs, so that intrinsic differences in embryonic growth due to metabolism could be investigated. I found that differences in juvenile and, potentially, embryonic growth rates existed among populations of A. carolinensis, both due to and independent of differences in egg size. Juveniles from the northernmost population were bigger not only due to larger egg size, but also due to faster juvenile growth and possibly differences in developmental stage of oviposition or conversion of egg mass to hatchling mass. Larger body size may hold a number of advantages in northern populations of this species, including starvation resistance through winters and better competitive access to food resources and warmer microhabitats.     

Role of adaptive and non-adaptive mechanisms forming complex patterns of genome size variation in six cytotypes of polyploid Allium oleraceum (Amaryllidaceae) on a continental scale

Background and Aims

Although the large variation in genome size among different species is widely acknowledged, the occurrence and extent of variation below the species level are still controversial and have not yet been satisfactorily analysed. The aim of this study was to assess genome size variation in six ploidy levels (2n = 3x–8x) of the polyploid Allium oleraceum over a large geographical gradient and to search for potential interpretations of the size variation.

Methods

The genome sizes of 407 individuals of A. oleraceum collected from 114 populations across Europe were determined by flow cytometry using propidium iodide staining. The genome size variation was correlated with spatial, climatic and habitat variables.

Key Results

The mean holoploid genome size (2C DNA) was 42·49, 52·14, 63·34, 71·94, 85·51 and 92·12 pg at the tri-, tetra-, penta-, hexa-, hepta- and octoploid levels, respectively. Genome size varied from a minimum of 2·3 % in the octoploids to a maximum of 18·3 % in the tetraploids. Spatial structuring of genome size was observed within the tetra- and pentaploids, where 2C DNA significantly increased with both latitude and longitude, and correlated with several climatic variables, suggesting a gradient of continentality. Genome size in hexaploids showed low variation, weak correlation with climatic variables and no spatial structuring. Downsizing in monoploid genome size was observed between all cytotypes except for heptaploids. Splitting populations into western and eastern European groups resulted in strong differences in monoploid genome size between groups in tetra- and pentaploids but not in hexaploids. The monoploid genome sizes of the cytotypes were similar in the western group but diverged in the eastern group.

Conclusions

Complex patterns of holoploid and monoploid genome size variation found both within and between A. oleraceum cytotypes are most likely the result of several interacting factors, including different evolutionary origins of cytotypes via hybridization of parental combinations with different genome sizes in the south-western and south-eastern part of Europe, introgression between cytotypes, and antropic dispersal. The role of broad-scale and fine-scale environmental variables in shaping genome size is probably of minor importance in A. oleraceum.     

Temperature and clinal variation in larval growth efficiency in Drosophila melanogaster

Geographic clines in ectotherm species including Drosophila melanogaster have been found throughout the world, with genetically larger body size and shorter development time occurring at high latitudes. Temperature is thought to play a major role in the evolution of this clinal variation. Laboratory thermal selection has effects similar to those seen in geographical clines. Evolution at low temperatures results in more rapid development to larger adult flies. This study investigated the effects of geographical origin and experimental temperature on larval growth efficiency in D. melanogaster. Larvae from populations that had evolved at high latitudes were found to use limited food more efficiently, so that the overall adult body size achieved was larger. Larvae reared at a lower experimental temperature (18 °C) used food more efficiently than those reared at a higher temperature (25 °C). The increases in growth efficiency found in populations from high latitudes could explain their increased body size and more rapid development.     

Variation in body size and metamorphic traits of Iberian spadefoot toads over a short geographic distance

Determinants of geographic variation in body size are often poorly understood, especially in organisms with complex life cycles. We examined patterns of adult body size and metamorphic traits variation in Iberian spadefoot toad ( Pelobates cultripes ) populations, which exhibit an extreme reduction in adult body size, 71.6% reduction in body mass, within just about 30 km at south-western Spain. We hypothesized that size at and time to metamorphosis would be predictive of the spatial pattern observed in adult body size. Larvae from eight populations were raised in a common garden experiment at two different larval densities that allow to differentiate whether population divergence was genetically based or was simply a reflection of environmental variation and, in addition, whether this population divergence was modulated by differing crowding larval environments. Larger adult size populations had higher larval growth rates, attaining larger sizes at metamorphosis, and exhibited higher survival than smaller-sized populations at both densities, although accentuated at a low larval density. These population differences appeared to be consistent once embryo size variation was controlled for, suggesting that this phenotypic divergence is not due to maternal effects. Our results suggest considerable genetic differentiation in metamorphic traits that parallels and may be a causal determinant of geographic variation in adult body size.