Phenotypic plasticity,but not adaptive tracking,underlies seasonal variation in post‐cold hardening freeze tolerance of Drosophila melanogaster

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short‐lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post‐cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm‐caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post‐cold hardening freeze tolerance. Comparing indoor and field‐caught flies and their laboratory‐reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field‐ and laboratory‐based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post‐cold hardening freeze tolerance results from environmental influences and not genetic changes.     

Contrasting acclimation responses to elevated CO2 and warming between an evergreen and a deciduous boreal conifer

Rising atmospheric carbon dioxide (CO₂) concentrations may warm northern latitudes up to 8°C by the end of the century. Boreal forests play a large role in the global carbon cycle, and the responses of northern trees to climate change will thus impact the trajectory of future CO₂ increases. We grew two North American boreal tree species at a range of future climate conditions to assess how growth and carbon fluxes were altered by high CO₂ and warming. Black spruce (Picea mariana, an evergreen conifer) and tamarack (Larix laricina, a deciduous conifer) were grown under ambient (407 ppm) or elevated CO₂ (750 ppm) and either ambient temperatures, a 4°C warming, or an 8°C warming. In both species, the thermal optimum of net photosynthesis (T_optA) increased and maximum photosynthetic rates declined in warm‐grown seedlings, but the strength of these changes varied between species. Photosynthetic capacity (maximum rates of Rubisco carboxylation, V_cmax, and of electron transport, J_max) was reduced in warm‐grown seedlings, correlating with reductions in leaf N and chlorophyll concentrations. Warming increased the activation energy for V_cmax and J_max (E_aV and E_aJ, respectively) and the thermal optimum for J_max. In both species, the T_optA was positively correlated with both E_aV and E_aJ, but negatively correlated with the ratio of J_max/V_cmax. Respiration acclimated to elevated temperatures, but there were no treatment effects on the Q₁₀ of respiration (the increase in respiration for a 10°C increase in leaf temperature). A warming of 4°C increased biomass in tamarack, while warming reduced biomass in spruce. We show that climate change is likely to negatively affect photosynthesis and growth in black spruce more than in tamarack, and that parameters used to model photosynthesis in dynamic global vegetation models (E_aV and E_aJ) show no response to elevated CO₂.     

Physiological and regulatory underpinnings of geographic variation in reptilian cold tolerance across a latitudinal cline

Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CT_min), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co‐expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity‐limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Genetic mapping shows intraspecific variation and transgressive segregation for caterpillar‐induced aphid resistance in maize

Plants in nature have inducible defences that sometimes lead to targeted resistance against particular herbivores, but susceptibility to others. The metabolic diversity and genetic resources available for maize (Zea mays) make this a suitable system for a mechanistic study of within‐species variation in such plant‐mediated interactions between herbivores. Beet armyworms (Spodoptera exigua) and corn leaf aphids (Rhopalosiphum maidis) are two naturally occurring maize herbivores with different feeding habits. Whereas chewing herbivore‐induced methylation of 2,4‐dihydroxy‐7‐methoxy‐1,4‐benzoxazin‐3‐one glucoside (DIMBOA‐Glc) to form 2‐hydroxy‐4,7‐dimethoxy‐1,4‐benzoxazin‐3‐one glucoside (HDMBOA‐Glc) promotes caterpillar resistance, lower DIMBOA‐Glc levels favour aphid reproduction. Thus, caterpillar‐induced DIMBOA‐Glc methyltransferase activity in maize is predicted to promote aphid growth. To test this hypothesis, the impact of S. exigua feeding on R. maidis progeny production was assessed using seventeen genetically diverse maize inbred lines. Whereas aphid progeny production was increased by prior caterpillar feeding on lines B73, Ki11, Ki3 and Tx303, it decreased on lines Ky21, CML103, Mo18W and W22. Genetic mapping of this trait in a population of B73 × Ky21 recombinant inbred lines identified significant quantitative trait loci on maize chromosomes 1, 7 and 10. There is a transgressive segregation for aphid resistance, with the Ky21 alleles on chromosomes 1 and 7 and the B73 allele on chromosome 10 increasing aphid progeny production. The chromosome 1 QTL coincides with a cluster of three maize genes encoding benzoxazinoid O‐methyltransferases that convert DIMBOA‐Glc to HDMBOA‐Glc. Gene expression studies and benzoxazinoid measurements indicate that S. exigua ‐induced responses in this pathway differentially affect R. maidis resistance in B73 and Ky21.     

Circannual variation in blood parasitism in a sub‐Saharan migrant passerine bird,the garden warbler

Knowing the natural dynamics of pathogens in migratory birds is important, for example, to understand the factors that influence the transport of pathogens to and their transmission in new geographical areas, whereas the transmission of other pathogens might be restricted to a specific area. We studied haemosporidian blood parasites of the genera Plasmodium, Haemoproteus and Leucocytozoon in a migratory bird, the garden warbler Sylvia borin. Birds were sampled in spring, summer and early autumn at breeding grounds in Sweden, on migration at Capri, Italy and on arrival and departure from wintering staging areas in West Africa: mapping recoveries of garden warblers ringed in Fennoscandia and Capri showed that these sites are most probably on the migratory flyway of garden warblers breeding at Kvismaren. Overall, haemosporidian prevalence was 39%, involving 24 different parasite lineages. Prevalence varied significantly over the migratory cycle, with relatively high prevalence of blood parasites in the population on breeding grounds and at the onset of autumn migration, followed by marked declines in prevalence during migration both on spring and autumn passage. Importantly, we found that when examining circannual variation in the different lineages, significantly different prevalence profiles emerged both between and within genera. Our results suggest that differences in prevalence profiles are the result of either different parasite transmission strategies or coevolution between the host and the various parasite lineages. When separating parasites into common vs. rare lineages, we found that two peaks in the prevalence of rare parasites occur; on arrival at Swedish breeding grounds, and after the wintering period in Africa. Our results stress the importance of appropriate taxonomic resolution when examining host‐parasite interactions, as variation in prevalence both between and within parasite genera can show markedly different patterns.     

Ecological and evolutionary factors of intraspecific variation in inducible defenses: Insights gained from Daphnia experiments

Phenotypic variation among individuals and species is a fundamental principle of natural selection. In this review, we focus on numerous experiments involving the model species Daphnia (Crustacea) and categorize the factors, especially secondary ones, affecting intraspecific variations in inducible defense. Primary factors, such as predator type and density, determine the degree to which inducible defense expresses and increases or decreases. Secondary factors, on the other hand, act together with primary factors to inducible defense or without primary factors on inducible defense. The secondary factors increase intraspecies variation in inducible defense, and thus, the level of adaptation of organisms varies within species. Future research will explore the potential for new secondary factors, as well as the relative importance between factors needs to be clarified.     

Geographic variation in fitness‐related traits of the bladderwrack Fucus vesiculosus along the Baltic Sea‐North Sea salinity gradient

In the course of the ongoing global intensification and diversification of human pressures, the study of variation patterns of biological traits along environmental gradients can provide relevant information on the performance of species under shifting conditions. The pronounced salinity gradient, co‐occurrence of multiple stressors, and accelerated rates of change make the Baltic Sea and its transition to North Sea a suitable region for this type of study. Focusing on the bladderwrack Fucus vesiculosus, one of the main foundation species on hard‐bottoms of the Baltic Sea, we analyzed the phenotypic variation among populations occurring along 2,000 km of coasts subjected to salinities from 4 to >30 and a variety of other stressors. Morphological and biochemical traits, including palatability for grazers, were recorded at 20 stations along the Baltic Sea and four stations in the North Sea. We evaluated in a common modeling framework the relative contribution of multiple environmental drivers to the observed trait patterns. Salinity was the main and, in some cases, the only environmental driver of the geographic trait variation in F. vesiculosus. The decrease in salinity from North Sea to Baltic Sea stations was accompanied by a decline in thallus size, photosynthetic pigments, and energy storage compounds, and affected the interaction of the alga with herbivores and epibiota. For some traits, drivers that vary locally such as wave exposure, light availability or nutrient enrichment were also important. The strong genetic population structure in this macroalgae might play a role in the generation and maintenance of phenotypic patterns across geographic scales. In light of our results, the desalination process projected for the Baltic Sea could have detrimental impacts on F. vesiculosus in areas close to its tolerance limit, affecting ecosystem functions such as habitat formation, primary production, and food supply.     

Climate‐dependent variation in cold tolerance of weedy rice and rice mediated by OsICE1 promoter methylation

The mechanisms by which weedy rice (Oryza sativa f. spontanea) has adapted to endure low‐temperature stress in northern latitudes remain unresolved. In this study, we assessed cold tolerance of 100 rice varieties and 100 co‐occurring weedy rice populations, which were sampled across a broad range of climates in China. A parallel pattern of latitude‐dependent variation in cold tolerance was detected in cultivated rice and weedy rice. At the molecular level, differential cold tolerance was strongly correlated with relative expression levels of CBF cold response pathway genes and with methylation levels in the promoter region of OsICE1, a regulator of this pathway. Among all methylated cytosine sites of the OsICE1 promoter, levels of CHG and CHH methylation were found to be significantly correlated with cold tolerance among accessions. Furthermore, within many of the collection locales, weedy rice shared identical or near‐identical OsICE1 methylation patterns with co‐occurring cultivated rice. These findings provide new insights on the possible roles that methylation variation in the OsICE1 promoter may play in cold tolerance, and they suggest that weedy rice can rapidly acquire cold tolerance via methylation patterns that are shared with co‐occurring rice cultivars.     

Geographic variation in thermal sensitivity of early life traits in a widespread reptile

Taxa with large geographic distributions generally encompass diverse macroclimatic conditions, potentially requiring local adaptation and/or phenotypic plasticity to match their phenotypes to differing environments. These eco‐evolutionary processes are of particular interest in organisms with traits that are directly affected by temperature, such as embryonic development in oviparous ectotherms. Here we examine the spatial distribution of fitness‐related early life phenotypes across the range of a widespread vertebrate, the painted turtle (Chrysemys picta). We quantified embryonic and hatchling traits from seven locations (in Idaho, Minnesota, Oregon, Illinois, Nebraska, Kansas, and New Mexico) after incubating eggs under constant conditions across a series of environmentally relevant temperatures. Thermal reaction norms for incubation duration and hatchling mass varied among locations under this common‐garden experiment, indicating genetic differentiation or pre‐ovulatory maternal effects. However, latitude, a commonly used proxy for geographic variation, was not a strong predictor of these geographic differences. Our findings suggest that this macroclimatic proxy may be an unreliable surrogate for microclimatic conditions experienced locally in nests. Instead, complex interactions between abiotic and biotic factors likely drive among‐population phenotypic variation in this system. Understanding spatial variation in key life‐history traits provides an important perspective on adaptation to contemporary and future climatic conditions.     

Climate‐driven reduction of genetic variation in plant phenology alters soil communities and nutrient pools

We examined the hypothesis that climate‐driven evolution of plant traits will influence associated soil microbiomes and ecosystem function across the landscape. Using a foundation tree species, Populus angustifolia, observational and common garden approaches, and a base population genetic collection that spans 17 river systems in the western United States, from AZ to MT, we show that (a) as mean annual temperature (MAT) increases, genetic and phenotypic variation for bud break phenology decline; (b) soil microbiomes, soil nitrogen (N), and soil carbon (C) vary in response to MAT and conditioning by trees; and (c) with losses of genetic variation due to warming, population‐level regulation of community and ecosystem functions strengthen. These results demonstrate a relationship between the potential evolutionary response of populations and subsequent shifts in ecosystem function along a large temperature gradient.     

INDUCER OF CBF EXPRESSION 1 integrates cold signals into FLOWERING LOCUS C‐mediated flowering pathways in Arabidopsis

Plants constantly monitor changes in photoperiod and temperature throughout the year to synchronize flowering with optimal environmental conditions. In the temperate zones, both photoperiod and temperature fluctuate in a somewhat predictable manner through the seasons, although a transient shift to low temperature is also encountered during changing seasons, such as early spring. Although low temperatures are known to delay flowering by inducing the floral repressor FLOWERING LOCUS C (FLC), it is not fully understood how temperature signals are coordinated with photoperiodic signals in the timing of seasonal flowering. Here, we show that the cold signaling activator INDUCER OF CBF EXPRESSION 1 (ICE1), FLC and the floral promoter SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) constitute an elaborate signaling network that integrates cold signals into flowering pathways. The cold‐activated ICE1 directly induces the gene encoding FLC, which represses SOC1 expression, resulting in delayed flowering. In contrast, under floral promotive conditions, SOC1 inhibits the binding of ICE1 to the promoters of the FLC gene, inducing flowering with a reduction of freezing tolerance. These observations indicate that the ICE1‐FLC‐SOC1 signaling network contributes to the fine‐tuning of flowering during changing seasons.     

Phenology and species determine growing‐season albedo increase at the altitudinal limit of shrub growth in the sub‐Arctic

Arctic warming is resulting in reduced snow cover and increased shrub growth, both of which have been associated with altered land surface–atmospheric feedback processes involving sensible heat flux, ground heat flux and biogeochemical cycling. Using field measurements, we show that two common Arctic shrub species (Betula glandulosa and Salix pulchra), which are largely responsible for shrub encroachment in tundra, differed markedly in albedo and that albedo of both species increased as growing season progressed when measured at their altitudinal limit. A moveable apparatus was used to repeatedly measure albedo at six precise spots during the summer of 2012, and resampled in 2013. Contrary to the generally accepted view of shrub‐covered areas having low albedo in tundra, full‐canopy prostrate B. glandulosa had almost the highest albedo of all surfaces measured during the peak of the growing season. The higher midsummer albedo is also evident in localized MODIS albedo aggregated from 2000 to 2013, which displays a similar increase in growing‐season albedo. Using our field measurements, we show the ensemble summer increase in tundra albedo counteracts the generalized effect of earlier spring snow melt on surface energy balance by approximately 40%. This summer increase in albedo, when viewed in absolute values, is as large as the difference between the forest and tundra transition. These results indicate that near future (<50 years) changes in growing‐season albedo related to Arctic vegetation change are unlikely to be particularly large and might constitute a negative feedback to climate warming in certain circumstances. Future efforts to calculate energy budgets and a sensible heating feedback in the Arctic will require more detailed information about the relative abundance of different ground cover types, particularly shrub species and their respective growth forms and phenology.     

Epigenetic variation predicts regional and local intraspecific functional diversity in a perennial herb

The ecological significance of epigenetic variation has been generally inferred from studies on model plants under artificial conditions, but the importance of epigenetic differences between individuals as a source of intraspecific diversity in natural plant populations remains essentially unknown. This study investigates the relationship between epigenetic variation and functional plant diversity by conducting epigenetic (methylation‐sensitive amplified fragment length polymorphisms, MSAP) and genetic (amplified fragment length polymorphisms, AFLP) marker–trait association analyses for 20 whole‐plant, leaf and regenerative functional traits in a large sample of wild‐growing plants of the perennial herb Helleborus foetidus from ten sampling sites in south‐eastern Spain. Plants differed widely in functional characteristics, and exhibited greater epigenetic than genetic diversity, as shown by per cent polymorphism of MSAP fragments (92%) or markers (69%) greatly exceeding that for AFLP ones (41%). After controlling for genetic structuring and possible cryptic relatedness, every functional trait considered exhibited a significant association with at least one AFLP or MSAP marker. A total of 27 MSAP (13.0% of total) and 12 AFLP (4.4%) markers were involved in significant associations, which explained on average 8.2% and 8.0% of trait variance, respectively. Individual MSAP markers were more likely to be associated with functional traits than AFLP markers. Between‐site differences in multivariate functional diversity were directly related to variation in multilocus epigenetic diversity after multilocus genetic diversity was statistically accounted for. Results suggest that epigenetic variation can be an important source of intraspecific functional diversity in H. foetidus, possibly endowing this species with the capacity to exploit a broad range of ecological conditions despite its modest genetic diversity.     

Differential remodeling of the lipidome during cold acclimation in natural accessions of Arabidopsis thaliana

Freezing injury is a major factor limiting the geographical distribution of plant species and the growth and yield of crop plants. Plants from temperate climates are able to increase their freezing tolerance during exposure to low but non‐freezing temperatures in a process termed cold acclimation. Damage to cellular membranes is the major cause of freezing injury in plants, and membrane lipid composition is strongly modified during cold acclimation. Forward and reverse genetic approaches have been used to probe the role of specific lipid‐modifying enzymes in the freezing tolerance of plants. In the present paper we describe an alternative ecological genomics approach that relies on the natural genetic variation within a species. Arabidopsis thaliana has a wide geographical range throughout the Northern Hemisphere with significant natural variation in freezing tolerance that was used for a comparative analysis of the lipidomes of 15 Arabidopsis accessions using ultra‐performance liquid chromatography coupled to Fourier‐transform mass spectrometry, allowing the detection of 180 lipid species. After 14 days of cold acclimation at 4°C the plants from most accessions had accumulated massive amounts of storage lipids, with most of the changes in long‐chain unsaturated triacylglycerides, while the total amount of membrane lipids was only slightly changed. Nevertheless, major changes in the relative amounts of different membrane lipids were also evident. The relative abundance of several lipid species was highly correlated with the freezing tolerance of the accessions, allowing the identification of possible marker lipids for plant freezing tolerance.     

Local environment,not local adaptation,drives leaf‐out phenology in common gardens along an elevational gradient in Acadia National Park,Maine

Premise of the Study

Climate‐driven changes in phenology are substantially affecting ecological relationships and ecosystem processes. The role of variation among species has received particular attention; for example, variation among species’ phenological responses to climate can disrupt trophic interactions and can influence plant performance. Variation within species in phenological responses to climate, however, has received much less attention, despite its potential role in ecological interactions and local adaptation to climate change.

Methods

We constructed three common gardens across an elevation gradient on Cadillac Mountain in Acadia National Park, Maine, to test population‐level responses in leaf‐out phenology in a reciprocal transplant experiment. The experiment included three native species: low bush blueberry (Vaccinium angustifolium), sheep's laurel (Kalmia angustifolia), and three‐toothed cinquefoil (Sibbaldiopsis tridentata).

Key Results

Evidence for local adaptation of phenological response to temperature varied among the species, but was weak for all three. Rather, variation in phenological response to temperature appeared to be driven by local microclimate at each garden site and year‐to‐year variation in temperature.

Conclusions

Population‐level adaptations in leaf‐out phenology appear to be relatively unimportant for these species in Acadia National Park, perhaps a reflection of strong genetic mixing across elevations, or weak differences in selection on phenological response to spring temperatures at different elevations. These results concur with other observational data in Acadia and highlight the utility of experimental approaches to understand the importance of annual and local site variation in affecting phenology both among and within plant species.     

Behavioral constraints on local adaptation and counter‐gradient variation: Implications for climate change

Resource allocation to growth, reproduction, and body maintenance varies within species along latitudinal gradients. Two hypotheses explaining this variation are local adaptation and counter‐gradient variation. The local adaptation hypothesis proposes that populations are adapted to local environmental conditions and are therefore less adapted to environmental conditions at other locations. The counter‐gradient variation hypothesis proposes that one population out performs others across an environmental gradient because its source location has greater selective pressure than other locations. Our study had two goals. First, we tested the local adaptation and counter‐gradient variation hypotheses by measuring effects of environmental temperature on phenotypic expression of reproductive traits in the burying beetle, Nicrophorus orbicollis Say, from three populations along a latitudinal gradient in a common garden experimental design. Second, we compared patterns of variation to evaluate whether traits covary or whether local adaptation of traits precludes adaptive responses by others. Across a latitudinal range, N. orbicollis exhibits variation in initiating reproduction and brood sizes. Consistent with local adaptation: (a) beetles were less likely to initiate breeding at extreme temperatures, especially when that temperature represents their source range; (b) once beetles initiate reproduction, source populations produce relatively larger broods at temperatures consistent with their local environment. Consistent with counter‐gradient variation, lower latitude populations were more successful at producing offspring at lower temperatures. We found no evidence for adaptive variation in other adult or offspring performance traits. This suite of traits does not appear to coevolve along the latitudinal gradient. Rather, response to selection to breed within a narrow temperature range may preclude selection on other traits. Our study highlights that N. orbicollis uses temperature as an environmental cue to determine whether to initiate reproduction, providing insight into how behavior is modified to avoid costly reproductive attempts. Furthermore, our results suggest a temperature constraint that shapes reproductive behavior.     

Colony‐level behavioural variation correlates with differences in expression of the foraging gene in red imported fire ants

Among social insects, colony‐level variation is likely to be widespread and has significant ecological consequences. Very few studies, however, have documented how genetic factors relate to behaviour at the colony level. Differences in expression of the foraging gene have been associated with differences in foraging and activity of a wide variety of organisms. We quantified expression of the red imported fire ant foraging gene (sifor) in workers from 21 colonies collected across the natural range of Texas fire ant populations, but maintained under standardized, environmentally controlled conditions. Colonies varied significantly in their behaviour. The most active colonies had up to 10 times more active foragers than the least active colony and more than 16 times as many workers outside the nest. Expression differences among colonies correlated with this colony‐level behavioural variation. Colonies with higher sifor expression in foragers had, on average, significantly higher foraging activity, exploratory activity and recruitment to nectar than colonies with lower expression. Expression of sifor was also strongly correlated with worker task (foraging vs. working in the interior of the nest). These results provide insight into the genetic and physiological processes underlying collective differences in social behaviour. Quantifying variation in expression of the foraging gene may provide an important tool for understanding and predicting the ecological consequences of colony‐level behavioural variation.