Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Seasonal and altitudinal variation in dorsal skin reflectance and thermic rates in a high-altitude montane lizard

Temperature is one of the most important factors in the life histories of ectotherms, as body temperature has an undeniable effect on growth, activity, and reproduction. Lizards have a wide variety of strategies to acquire and maintain body temperature in an optimal range. The "Thermal Melanism Hypothesis" proposes that individuals with lower skin reflectance can heat up faster as a result of absorbing more solar radiation compared to lighter conspecifics. Therefore, having a darker coloration might be advantageous in cold habitats. Dorsal skin reflectance has been found to change rapidly with body temperature in several lizard species, and it can also vary over longer, seasonal time scales. These variations may be important in thermoregulation, especially in lizards that inhabit areas with a large temperature variation during the year. Here, we study how dorsal reflectance fluctuates with body temperature and varies among seasons. We compared dorsal skin reflectance at three body temperature treatments, and measured thermal rates (i.e., heat and cool rate, thermic lapse, and net heat gain) by elevation (2500-4100 m) and seasons (spring, summer, and autumn) in the mesquite lizard, Sceloporus grammicus. Our results show that lizards were darker at high elevations and during the months with the lowest environmental temperatures. The rate of obtaining and retaining heat also varied during the year and was highest during the reproductive season. Our results indicate that the variation of dorsal skin reflectance and thermal rates follows a complex pattern in lizard populations and is affected by both elevation and season.

Fitness difference between two synonymous mutations of Phytophthora infestans ATP6 gene

BACKGROUND

Sequence variation produced by mutation provides the ultimate source of natural selection for species adaptation. Unlike nonsynonymous mutation, synonymous mutations are generally considered to be selectively neutral but accumulating evidence suggests they also contribute to species adaptation by regulating the flow of genetic information and the development of functional traits. In this study, we analysed sequence characteristics of ATP6, a housekeeping gene from 139 Phytophthora infestans isolates, and compared the fitness components including metabolic rate, temperature sensitivity, aggressiveness, and fungicide tolerance among synonymous mutations.

RESULTS

We found that the housekeeping gene exhibited low genetic variation and was represented by two major synonymous mutants at similar frequency (0.496 and 0.468, respectively). The two synonymous mutants were generated by a single nucleotide substitution but differed significantly in fitness as well as temperature-mediated spatial distribution and expression. The synonymous mutant ending in AT was more common in cold regions and was more expressed at lower experimental temperature than the synonymous mutant ending in GC and vice versa.

CONCLUSION

Our results are consistent with the argument that synonymous mutations can modulate the adaptive evolution of species including pathogens and have important implications for sustainable disease management, especially under climate change.

Fish diversity of Colombian Andes-Amazon streams at the end of conflict is a reference for conservation before increased land use

Reference conditions are difficult to find in the Anthropocene but essential for effective biodiversity conservation. Aquatic ecosystems in the Andes-Amazon transition zone of Colombia are now at high risk due to expanded human activities after peace agreements in 2016 ended armed conflict because lands formerly controlled by FARC and other armed groups are now prone to agricultural and urban expansion. Particularly, expanding human land use may reduce fish diversity across the altitudinal gradient, especially in the premontane streams (i.e., <500 m a.s.l.) because lands are more amenable to human use than at greater altitudes. We evaluated fish α-diversity (measured as species richness, total abundance, and effective species number) and β-diversity (spatial and temporal) in 12 sites over 8 years bracketing the end of armed conflict. All α-diversity and β-diversity analyses were evaluated relative to categorical altitude (< or >500 m) and continuous altitude. Strong differences in fish community structure among sites occurred as a function of altitude. Fish communities exhibit altitudinal biodiversity gradients that are consistent in space and time, and that need to be accounted for conservation and management considerations. Our results provide a reference to identify short- and long-term changes due to impending human land use at a critical moment for the conservation of tropical fish diversity. Similar studies in other areas of the upper Amazon Basin are needed to evaluate effects of subsequent human activities on diversity patterns and our study area to compare to reference conditions reported here.

Thermal limits along tropical elevational gradients: Poison frog tadpoles show plasticity but maintain divergence across elevation

Temperature is arguably one of the most critical environmental factors impacting organisms at molecular, organismal, and ecological levels. Temperature variation across elevation may cause divergent selection in physiological critical thermal limits (CTMAX and CTMIN). Generally, high elevation populations are predicted to withstand lower environmental temperatures than low elevation populations. Organisms can also exhibit phenotypic plasticity when temperature varies, although theory and empirical evidence suggest that tropical ectotherms have relatively limited ability to acclimate. To study the effect of temperature variation along elevational transects on thermal limits, we measured CTMAX and CTMIN of 934 tadpoles of a poison frog species, Epipedobates anthonyi, along two elevational gradients (200-1700 m asl) in southwestern Ecuador to investigate their thermal tolerance across elevation. We also tested if tadpoles could plastically shift their critical thermal limits in response to exposure to different temperatures representing the range of temperatures they experience in nature (20 °C, 24 °C, and 28 °C). Overall, we found that CTMAX did not change across elevation. In contrast, CTMIN was lower at higher elevations, suggesting that elevational variation in temperature influences this thermal trait. Moreover, all populations shifted their CTMAX and CTMIN according to treatment temperatures, demonstrating an acclimation response. Overall, trends in CTMIN among high, mid, and low elevation populations were maintained despite plastic responses to treatment temperature. These results demonstrate that, for tadpoles of E. anthonyi across tropical elevational gradients, temperature acts as a selective force for CTMIN, even when populations show acclimation abilities in both, CTMAX and CTMIN. Our findings advance our understanding on how environmental variation affects organisms' evolutionary trajectories and their abilities to persist in a changing climate in a tropical biodiversity hotspot.

Unraveling Adaptive Evolutionary Divergence at Microgeographic Scales

AbstractStriking examples of local adaptation at fine geographic scales are increasingly being documented in natural populations. However, the relative contributions made by natural selection, phenotype-dependent dispersal (when individuals disperse with respect to a habitat preference), and mate preference in generating and maintaining microgeographic adaptation and divergence are not well studied. Here, we develop quantitative genetics models and individual-based simulations (IBSs) to uncover the evolutionary forces that possibly drive microgeographic divergence. We also perform Bayesian estimation of the parameters in our IBS using empirical data on habitat-specific variation in bill morphology in the island scrub-jay (Aphelocoma insularis) to apply our models to a natural system. We find that natural selection and phenotype-dependent dispersal can generate the patterns of divergence we observe in the island scrub-jay. However, mate preference for a mate with similar bill morphology, even though observed in the species, does not play a significant role in driving divergence. Our modeling approach provides insights into phenotypic evolution occurring over small spatial scales relative to dispersal ranges, suggesting that adaptive divergence at microgeographic scales may be common across a wider range of taxa than previously thought. Our quantitative genetic models help to inform future theoretical and empirical work to determine how selection, habitat preference, and mate preference contribute to local adaptation and microgeographic divergence.

Genome-wide DNA methylation patterns in bumble bee (Bombus vosnesenskii) populations from spatial-environmental range extremes

Unraveling molecular mechanisms of adaptation to complex environments is crucial to understanding tolerance of abiotic pressures and responses to climatic change. Epigenetic variation is increasingly recognized as a mechanism that can facilitate rapid responses to changing environmental cues. To investigate variation in genetic and epigenetic diversity at spatial and thermal extremes, we use whole genome and methylome sequencing to generate a high-resolution map of DNA methylation in the bumble bee Bombus vosnesenskii. We sample two populations representing spatial and environmental range extremes (a warm southern low-elevation site and a cold northern high-elevation site) previously shown to exhibit differences in thermal tolerance and determine positions in the genome that are consistently and variably methylated across samples. Bisulfite sequencing reveals methylation characteristics similar to other arthropods, with low global CpG methylation but high methylation concentrated in gene bodies and in genome regions with low nucleotide diversity. Differentially methylated sites (n = 2066) were largely hypomethylated in the northern high-elevation population but not related to local sequence differentiation. The concentration of methylated and differentially methylated sites in exons and putative promoter regions suggests a possible role in gene regulation, and this high-resolution analysis of intraspecific epigenetic variation in wild Bombus suggests that the function of methylation in niche adaptation would be worth further investigation.

Influences of Maternal Weight and Geographic Factors on Offspring Traits of the Edible Dormouse in the NE of the Iberian Peninsula

The main goal of this study was to analyze the reproductive patterns of edible dormouse (Glis glis) populations in the northeast of the Iberian Peninsula using an 18-year period of data obtained from nest boxes collected between 2004 and 2021. The average litter size in Catalonia (Spain) was 5.5 ± 1.60 (range: 2-9, n = 131), with litter sizes between 5 and 7 pups as the more frequent. The overall mean weight in pink, grey and open eyes pups was 4.8 g/pup, 11.7 g/pup and 23.6 g/pup, respectively. No differences in offspring weights between sexes were found in any of the three age groups. Maternal body weight was positively associated with mean pup weight, whereas no correlation between the weight of the mothers and litter size was found. The trade-off between offspring number and size was not detected at birth. Regarding litter size variation across the geographic gradient (and their climatic gradient associated) from the southernmost populations of the Iberian Peninsula located in Catalonia to the Pyrenees region in Andorra, no evidence to suggest that geographic variables affect litter size was found, discarding (1) an investment in larger litters to compensate shorter seasons related to higher altitudes or northern latitudes, and (2) variation in litter size related to weather changes (e.g., temperature and precipitation) along latitudinal and/or altitudinal gradients.

Patterns and drivers of soil surface-dwelling Oribatida diversity along an altitudinal gradient on the Changbai Mountain, China

Distribution patterns of biodiversity and environmental interactions are dominant themes in ecology. In montane ecosystems, biodiversity is closely associated with altitudinal gradients. However, studies of biodiversity in montane ecosystems are focused on plants and vertebrates, with relatively less on invertebrates. Here, the present study used a Vortis arthropod suction sampler to explore the biodiversity patterns of soil surface-dwelling Oribatida and their drivers along an altitudinal gradient (600, 800, 1600, 2000, and 2300 m) from typical temperate forests, evergreen coniferous forests, subalpine birch forests to alpine tundra on the north slope of Changbai Mountain, Northeast China. Trichoribates berlesei, Platynothrus peltifer, and Oribatula tibialis were the dominant soil surface-dwelling species on Changbai Mountain. Generally, alpha diversity and beta diversity of soil surface-dwelling Oribatida decreased with the rising altitude, with a peaking density value at 2000 m. The result of beta diversity showed that the structures of community were more influenced by the species turnover component than the nestedness component. Nonmetric multidimensional scaling (NMDS) ordination showed that the community structure of soil surface-dwelling Oribatida varied significantly along the altitudinal gradient. The variance partitioning showed that the elevation and climatic conditions determined the soil surface-dwelling Oribatida community. Spatial filtering represented by geographic and elevation distances was particularly associated with soil surface-dwelling Oribatida community variation between altitudes on Changbai Mountain. However, the variation of the Oribatida community between adjacent altitudes was only associated with geographic distance. Our study provides supportive evidence for the biodiversity analyzing of soil surface-dwelling Oribatida in montane ecosystems along an altitudinal gradient.

Elevation and blood traits in the mesquite lizard: Are patterns repeatable between mountains?

Ecogeographical patterns describe predictable variation in phenotypic traits between ecological communities. For example, high-altitude animals are expected to show elevated hematological values as an adaptation to the lower oxygen pressure. Mountains act like ecological islands and therefore are considered natural laboratories. However, the majority of ecophysiological studies on blood traits lack replication that would allow us to infer if the pattern reported is a local event or whether it is a widespread pattern resulting from larger-scale ecological processes. In lizards, in fact, the increase of hematological values at high altitudes has received mixed support. Here, for the first time, we compare blood traits in lizards along elevational gradients with replication. We tested the repeatability of blood traits in mesquite lizards between different elevations in three different mountains from the Trans-Mexican Volcanic Belt. We measured hematocrit, hemoglobin concentration, mean corpuscular hemoglobin concentration, and erythrocyte size in blood samples of low, medium, and high-elevation lizards. We obtained similar elevational patterns between mountains, but the blood traits differed among mountains. Middle-altitude populations had greater oxygen-carrying capacity than lizards from low and high altitudes. The differences found between mountain systems could be the result of phenotypic plasticity or genetic differentiation as a consequence of abiotic factors not considered.

Temperature and livestock grazing trigger transcriptome responses in bumblebees along an elevational gradient

Climate and land-use changes cause increasing stress to pollinators but the molecular pathways underlying stress responses are poorly understood. Here, we analyzed the transcriptomic response of Bombus lucorum workers to temperature and livestock grazing. Bumblebees sampled along an elevational gradient, and from differently managed grassland sites (livestock grazing vs unmanaged) in the German Alps did not differ in the expression of genes known for thermal stress responses. Instead, metabolic energy production pathways were upregulated in bumblebees sampled in mid- or high elevations or during cool temperatures. Extensive grazing pressure led to an upregulation of genetic pathways involved in immunoregulation and DNA-repair. We conclude that widespread bumblebees are tolerant toward temperature fluctuations in temperate mountain environments. Moderate temperature increases may even release bumblebees from metabolic stress. However, transcriptome responses to even moderate management regimes highlight the completely underestimated complexity of human influence on natural pollinators.

•

Upregulation of energy metabolism pathways in Bombus lucorum with increasing elevation

•

Genes known for thermal stress responses did not change with increased elevation

•

Bombus lucorum are tolerant toward relatively broad temperature fluctuations

•

Grazing lead to an upregulation in genetic information processes in B. lucorum

Relationship between brain size and digestive tract length support the expensive-tissue hypothesis in Feirana quadranus

The brain is among the most energetically costly organs in the vertebrate body, while the size of the brain varies within species. The expensive-tissue hypothesis (ETH) predicts that increasing the size of another costly organ, such as the gut, should compensate for the cost of a small brain. Here, the ETH was tested by analyzing the relationship between brain size variation and digestive tract length in a Swelled-vented frog (Feirana quadranus). A total of 125 individuals across 10 populations ranging from 586 to 1,702 m a.s.l. from the Qinling-Daba Mountains were sampled. With the increase in altitude, the brain size decreases and the digestive tract length increases. Different brain regions do not change their relative size in a consistent manner. The sizes of telencephalon and cerebellum decrease with the increase in altitude, while the olfactory nerve increases its size at high altitudes. However, the olfactory bulb and optic tectum have no significant relationship with altitude. After controlling for snout-vent length (SVL), a significant negative correlation could be found between brain size and digestive tract length in F. quadranus. Therefore, the intraspecific variation of brain size follows the general patterns of ETH in this species. The results suggest that annual mean temperature and annual precipitation are environmental factors influencing the adaptive evolution of brain size and digestive tract length. This study also suggests that food composition, activity times, and habitat complexity are the potential reasons driving the adaptive evolution of brain size and digestive tract length.

Challenging ecogeographical rules: Phenotypic variation in the Mountain Treeshrew (Tupaia montana) along tropical elevational gradients

Bergmann’s and Allen’s rules were defined to describe macroecological patterns across latitudinal gradients. Bergmann observed a positive association between body size and latitude for endothermic species while Allen described shorter appendages as latitude increases. Almost two centuries later, there is still ongoing discussion about these patterns. Temperature, the common variable in these two rules, varies predictably across both latitude and elevation. Although these rules have been assessed extensively in mammals across latitude, particularly in regions with strong seasonality, studies on tropical montane mammals are scarce. We here test for these patterns and assess the variation of several other locomotory, diet-associated, body condition, and thermoregulatory traits across elevation in the Mountain Treeshrew (Tupaia montana) on tropical mountains in Borneo. Based on morphological measurements from both the field and scientific collections, we found a complex pattern: Bergmann’s rule was not supported in our tropical mountain system, since skull length, body size, and weight decreased from the lowest elevations (<1000 m) to middle elevations (2000–2500 m), and then increased from middle elevations to highest elevations. Allen’s rule was supported for relative tail length, which decreased with elevation, but not for ear and hindfoot length, with the former remaining constant and the latter increasing with elevation. This evidence together with changes in presumed diet-related traits (rostrum length, zygomatic breadth and upper tooth row length) along elevation suggest that selective pressures other than temperature, are playing a more important role shaping the morphological variation across the distribution of the Mountain Treeshrew. Diet, food acquisition, predation pressure, and/or intra- and inter-specific competition, are some of the potential factors driving the phenotypic variation of this study system. The lack of variation in body condition might suggest local adaptation of this species across its elevational range, perhaps due to generalist foraging strategies. Finally, a highly significant temporal effect was detected in several traits but not in others, representing the first phenotypic variation temporal trends described on treeshrews.

Anna's hummingbird (Calypte anna) physiological response to novel thermal and hypoxic conditions at high elevations

Many species have not tracked their thermal niches upslope as predicted by climate change, potentially because higher elevations are associated with abiotic challenges beyond temperature. To better predict whether organisms can continue to move upslope with rising temperatures, we need to understand their physiological performance when subjected to novel high-elevation conditions. Here, we captured Anna's hummingbirds - a species expanding their elevational distribution in concordance with rising temperatures - from across their current elevational distribution and tested their physiological response to novel abiotic conditions. First, at a central aviary within their current elevational range, we measured hovering metabolic rate to assess their response to oxygen conditions and torpor use to assess their response to thermal conditions. Second, we transported the hummingbirds to a location 1200 m above their current elevational range limit to test for an acute response to novel oxygen and thermal conditions. Hummingbirds exhibited lower hovering metabolic rates above their current elevational range limit, suggesting lower oxygen availability may reduce performance after an acute exposure. Alternatively, hummingbirds showed a facultative response to thermal conditions by using torpor more frequently and for longer. Finally, post-experimental dissection found that hummingbirds originating from higher elevations within their range had larger hearts, a potential plastic response to hypoxic environments. Overall, our results suggest lower oxygen availability and low air pressure may be difficult challenges to overcome for hummingbirds shifting upslope as a consequence of rising temperatures, especially if there is little to no long-term acclimatization. Future studies should investigate how chronic exposure and acclimatization to novel conditions, as opposed to acute experiments, may result in alternative outcomes that help organisms better respond to abiotic challenges associated with climate-induced range shifts.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Evidence for genetic isolation and local adaptation in the field cricket Gryllus campestris

Understanding how species can thrive in a range of environments is a central challenge for evolutionary ecology. There is strong evidence for local adaptation along large-scale ecological clines in insects. However, potential adaptation among neighbouring populations differing in their environment has been studied much less. We used RAD sequencing to quantify genetic divergence and clustering of ten populations of the field cricket Gryllus campestris in the Cantabrian Mountains of northern Spain, and an outgroup on the inland plain. Our populations were chosen to represent replicate high and low altitude habitats. We identified genetic clusters that include both high and low altitude populations indicating that the two habitat types do not hold ancestrally distinct lineages. Using common-garden rearing experiments to remove environmental effects, we found evidence for differences between high and low altitude populations in physiological and life-history traits. As predicted by the local adaptation hypothesis, crickets with parents from cooler (high altitude) populations recovered from periods of extreme cooling more rapidly than those with parents from warmer (low altitude) populations. Growth rates also differed between offspring from high and low altitude populations. However, contrary to our prediction that crickets from high altitudes would grow faster, the most striking difference was that at high temperatures, growth was fastest in individuals from low altitudes. Our findings reveal that populations a few tens of kilometres apart have independently evolved adaptations to their environment. This suggests that local adaptation in a range of traits may be commonplace even in mobile invertebrates at scales of a small fraction of species' distributions.

Brain size variation along altitudinal gradients in the Asiatic Toad (Bufo gargarizans)

Size changes in brain and brain regions along altitudinal gradients provide insight into the trade-off between energetic expenditure and cognitive capacity. We investigated the brain size variations of the Asiatic Toad (Bufo gargarizans) across altitudes from 700 m to 3,200 m. A total of 325 individuals from 11 sites and two transects were sampled. To reduce confounding factors, all sampling sites within each transect were within a maximum distance of 85 km and an altitudinal difference close to 2,000 m. Brains were dissected, and five regions were both measured directly and with 3D CT scan. There is a significant negative correlation between the relative whole-brain volume (to snout-vent length) and altitude. Furthermore, the relative volumes (to whole-brain volume) of optic tectum and cerebellum also decrease along the altitudinal gradients, while the telencephalon increases its relative volume along the gradients. Therefore, our results are mostly consistent with the expensive brain hypothesis and the functional constraint hypothesis. We suggest that most current hypotheses are not mutually exclusive and data supporting one hypothesis are often partially consistent with others. More studies on mechanisms are needed to explain the brain size evolution in natural populations.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Elevational Changes in Mormon Cricket Life Histories: Minimum Temperature for Nymphal Growth Declines With Elevation

As the mean temperature and the duration of the growing season decline with elevation, growth of immature insects should initiate at a lower temperature, but it should also be faster to complete development prior to season's end. Although flightless, Mormon crickets migrate in large aggregations across broad spatial and elevational distances that might limit adaptations to local environments. In addition, selection to be active at cooler temperatures might limit selection to maximize growth rate. I measured growth rate in controlled environments for nymphs from three populations that vary in altitude (87-2,688 m) but are similar in latitude (43.2-45.7°N). Growth rate increased significantly with mean rearing temperature between 22 and 30°C. The intercept of the regression of growth rate on temperature increased with elevation, whereas the slope did not change significantly. For any given rearing temperature, growth rate increased with elevation, which suggests that selection to initiate growth at cooler temperatures did not compromise growth rate. Body mass did not differ between the two lower elevations, whereas the highest elevation population had smaller hatchlings and adults. Critical thermal minimum (base temperature) declined with elevation (0.7°C per 1,000 m), and the degree days were 509 across all elevations. For pest management, a base temperature from midelevation of 15.3°C (60°F) and growing degree days of 509 (equivalent to 916 Fahrenheit-based degree days) are reasonable estimates for applications from sea level to 2,700 m.

To be small and dark is advantageous for gaining heat in mezquite lizards, Sceloporus grammicus (Squamata: Phrynosomatidae)

Body temperature is important in determining individual performance in ectotherms such as lizards. Environmental temperature decreases with increasing altitude, but nevertheless many lizards inhabit high-altitude environments. The ‘thermal melanism hypothesis’ proposes that a dark dorsal coloration enables darker individuals to heat up faster because they absorb more solar radiation and thus being darker may be advantageous in cold habitats. The aim of the present study is to evaluate how heating rate, cooling rate and net heat gain vary with body size and dorsal skin coloration in Sceloporus grammicus lizards along an altitudinal gradient. We measured these traits multiple times in the same individuals with a radiation heat source and spectrophotometry under laboratory conditions. Our results showed that S. grammicus lizards are smaller and darker at high elevations than at low elevations. In addition, the smallest and darkest lizards showed the greatest heating rate and net heat gain. Therefore, in S. grammicus, we suggest that small body size and dark dorsal coloration provide thermoregulatory benefits in high-altitude environments. Hence, this study supports the thermal melanism hypothesis in a lizard species under varied thermal environments.

Anthropogenic and natural barriers affect genetic connectivity in an Alpine butterfly

Dispersal is a key biological process serving several functions including connectivity among populations. Habitat fragmentation caused by natural or anthropogenic structures may hamper dispersal, thereby disrupting genetic connectivity. Investigating factors affecting dispersal and gene flow is important in the current era of anthropogenic global change, as dispersal comprises a vital part of a species' resilience to environmental change. Using finescale landscape genomics, we investigated gene flow and genetic structure of the Sooty Copper butterfly (Lycaena tityrus) in the Alpine Ötz valley system in Austria. We found surprisingly high levels of gene flow in L. tityrus across the region. Nevertheless, ravines, forests, and roads had effects on genetic structure, while rivers did not. The latter is surprising as roads and rivers have a similar width and run largely in parallel in our study area, pointing towards a higher impact of anthropogenic compared with natural linear structures. Additionally, we detected eleven loci potentially under thermal selection, including ones related to membranes, metabolism, and immune function. This study demonstrates the usefulness of molecular approaches in obtaining estimates of dispersal and population processes in the wild. Our results suggest that, despite high gene flow in the Alpine valley system investigated, L. tityrus nevertheless seems to be vulnerable to anthropogenically-driven habitat fragmentation. With anthropogenic rather than natural linear structures affecting gene flow, this may have important consequences for the persistence of species such as the butterfly studied here in altered landscapes.

Adaptation to elevation but limited local adaptation in an amphibian

We performed a reciprocal transplant experiment to estimate "parallel" adaptation to elevation and "unique" adaptation to local sites at the same elevation, using the frog Rana temporaria in the Swiss Alps. It is important to distinguish these two processes because they have different implications for population structure and ecological specialization. Larvae were reared from hatching to metamorphosis within enclosures installed in their pond of origin, in three foreign ponds at the same elevation, and in four ponds at different elevation (1500-2000 m higher or lower). There were two source populations from each elevation, and adults were held in a common environment for 1 year before they were crossed to produce offspring for the experiment. Fitness was a measure that integrated larval survival, development rate, and body size. Parallel adaptation to elevation was indicated by an advantage at the home elevation (11.5% fitness difference at low elevation and 47% at high elevation). This effect was stronger than that observed in most other studies, according to a survey of previous transplant experiments across elevation (N = 8 animal species and 71 plants). Unique local adaptation within elevational zones was only 0.3-0.7 times as strong as parallel adaptation, probably because gene flow is comparatively high among nearby wetlands at the same elevation. The home-elevation advantage may reduce gene flow across the elevational gradient and enable the evolution of habitat races specialized on elevation.

Myco- and photobiont associations in crustose lichens in the McMurdo Dry Valleys (Antarctica) reveal high differentiation along an elevational gradient

Climatically extreme regions such as the polar deserts of the McMurdo Dry Valleys (78° S) in Continental Antarctica are key areas for a better understanding of changes in ecosystems. Therefore, it is particularly important to analyze and communicate current patterns of biodiversity in these sensitive areas, where precipitation mostly occurs in form of snow and liquid water is rare. Humidity provided by dew, clouds, and fog are the main water sources, especially for rock-dwelling crustose lichens as one of the most common vegetation-forming organisms. We investigated the diversity and interaction specificity of myco-/photobiont associations of 232 crustose lichen specimens, collected along an elevational gradient (171–959 m a.s.l.) within the McMurdo Dry Valleys. The mycobiont species and photobiont OTUs were identified by using three markers each (nrITS, mtSSU, RPB1, and nrITS, psbJ-L, COX2). Elevation, positively associated with water availability, turned out to be the key factor explaining most of the distribution patterns of the mycobionts. Pairwise comparisons showed Lecidea cancriformis and Rhizoplaca macleanii to be significantly more common at higher elevations and Carbonea vorticosa and Lecidea polypycnidophora at lower elevations. Lichen photobionts were dominated by the globally distributed Trebouxia OTU, Tr_A02 which occurred at all habitats. Network specialization resulting from myco-/photobiont bipartite network structure varied with elevation and associated abiotic factors. Along an elevational gradient, the spatial distribution, diversity, and genetic variability of the lichen symbionts appear to be mainly influenced by improved water relations at higher altitudes.

Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee

Global temperature changes have emphasized the need to understand how species adapt to thermal stress across their ranges. Genetic mechanisms may contribute to variation in thermal tolerance, providing evidence for how organisms adapt to local environments. We determine physiological thermal limits and characterize genome-wide transcriptional changes at these limits in bumble bees using laboratory-reared Bombus vosnesenskii workers. We analyze bees reared from latitudinal (35.7-45.7°N) and altitudinal (7-2154 m) extremes of the species' range to correlate thermal tolerance and gene expression among populations from different climates. We find that critical thermal minima (CTMIN) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CTMAX) was invariant among populations. Concordant patterns are apparent in gene expression data, with regional differentiation following cold exposure, and expression shifts invariant among populations under high temperatures. Furthermore, we identify several modules of co-expressed genes that tightly correlate with critical thermal limits and temperature at the region of origin. Our results reveal that local adaptation in thermal limits and gene expression may facilitate cold tolerance across a species range, whereas high temperature responses are likely constrained, both of which may have implications for climate change responses of bumble bees.

Machine-learning-based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine-learning method (i.e., random forest) besides traditional distance-based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non-neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non-neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.

Distinct Patterns of Selective Sweep and Polygenic Adaptation in Evolve and Resequence Studies

In molecular population genetics, adaptation is typically thought to occur via selective sweeps, where targets of selection have independent effects on the phenotype and rise to fixation, whereas in quantitative genetics, many loci contribute to the phenotype and subtle frequency changes occur at many loci during polygenic adaptation. The sweep model makes specific predictions about frequency changes of beneficial alleles and many test statistics have been developed to detect such selection signatures. Despite polygenic adaptation is probably the prevalent mode of adaptation, because of the traditional focus on the phenotype, we are lacking a solid understanding of the similarities and differences of selection signatures under the two models. Recent theoretical and empirical studies have shown that both selective sweep and polygenic adaptation models could result in a sweep-like genomic signature; therefore, additional criteria are needed to distinguish the two models. With replicated populations and time series data, experimental evolution studies have the potential to identify the underlying model of adaptation. Using the framework of experimental evolution, we performed computer simulations to study the pattern of selected alleles for two models: 1) adaptation of a trait via independent beneficial mutations that are conditioned for fixation, that is, selective sweep model and 2) trait optimum model (polygenic adaptation), that is adaptation of a quantitative trait under stabilizing selection after a sudden shift in trait optimum. We identify several distinct patterns of selective sweep and trait optimum models in populations of different sizes. These features could provide the foundation for development of quantitative approaches to differentiate the two models.

Geographic variation in adult and embryonic desiccation tolerance in a terrestrial-breeding frog

Intraspecific variation in the ability of individuals to tolerate environmental perturbations is often neglected when considering the impacts of climate change. Yet this information is potentially crucial for mitigating deleterious effects of climate change on threatened species. Here we assessed patterns of intraspecific variation in desiccation tolerance in the frog Pseudophryne guentheri, a terrestrial-breeding species experiencing a drying climate. Adult frogs were collected from six populations across a rainfall gradient and their dehydration and rehydration rates were assessed. We also compared desiccation tolerance of embryos and hatchlings originating from within-population parental crosses from four of the populations. Embryos were reared on soil at three soil-water potentials and their desiccation tolerance was assessed across a range of traits. We found significant and strong patterns of intraspecific variation in almost all traits, both in adults and first-generation offspring. Adult frogs exhibited clinal variation in their water balance responses, with populations from drier sites both dehydrating and rehydrating more slowly compared to frogs from more mesic sites. Similarly, desiccation tolerance of first-generation offspring was significantly greater in populations from xeric sites. Our findings suggest that populations within this species will respond differently to the regional reduction in rainfall predicted by climate change models.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Paper wasps are darker at high elevation

High mountains are harsh environments in which colder temperatures and higher levels of UV-B radiation are common. These abiotic conditions strongly affect animals' biology, often constraining their survival and reproduction. As a result, adaptations to live in such habitats are expected to evolve. Body color is thought to be adaptive to the environment that animals experience. Tegument melanization improves heat gain and provides photoprotection. Therefore, at high elevation, ectotherms are expected to be darker (well-melanized). We test this prediction in the paper wasp Agelaia pallipes (Hymenoptera: Vespidae), a species distributed across an elevational gradient in the Colombian Andes. We used Malaise traps and sampled a total of 146 wasps along nine elevations, ranging from 2,600-3,380 m above sea level. Standard digital photography was used to measure the body luminance and colour patterning in different body parts of dry-preserved specimens. There was striking variation in body luminance (darker and lighter), color patterning (patched, smoothed, homogeneous) and surface texture (shiny and matte), but the kind and degree of variation depended on the body part examined. Wasps from higher elevations had darker thoraces, confirming our prediction. Besides, at high elevation, the frequency of wasps with a matte rather than a shiny face strongly increased. Overall, our findings support the thermal melanism hypothesis and suggest that intraspecific color variation might be an adaptation to the environment of paper wasps.

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life-history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome-wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short-winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short-winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high-elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.

Crypsis Decreases with Elevation in a Lizard

Predation usually selects for visual crypsis, the colour matching between an animal and its background. Geographic co-variation between animal and background colourations is well known, but how crypsis varies along elevational gradients remains unknown. We predict that dorsal colouration in the lizard Psammodromus algirus should covary with the colour of bare soil—where this lizard is mainly found—along a 2200 m elevational gradient in Sierra Nevada (SE Spain). Moreover, we predict that crypsis should decrease with elevation for two reasons: (1) Predation pressure typically decreases with elevation, and (2) at high elevation, dorsal colouration is under conflicting selection for both crypsis and thermoregulation. By means of standardised photographies of the substratum and colourimetric measurements of lizard dorsal skin, we tested the colour matching between lizard dorsum and background. We found that, along the gradient, lizard dorsal colouration covaried with the colouration of bare soil, but not with other background elements where the lizard is rarely detected. Moreover, supporting our prediction, the degree of crypsis against bare soil decreased with elevation. Hence, our findings suggest local adaptation for crypsis in this lizard along an elevational gradient, but this local adaptation would be hindered at high elevations.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

The role of dispersal, selection, and timing of sampling on the false discovery rate of loci under selection during geographic range expansion

Identifying adaptive loci is important to understand the evolutionary potential of species undergoing range expansion. However, in expanding populations, spatial demographic processes such as allele surfing can create spatial patterns of neutral genetic variation that appear similar to those generated through adaptive processes. As a result, the false discovery rate of adaptive loci may be inflated in landscape genomic analyses. Here, we take a simulation modelling approach to investigate how range expansion affects our ability to correctly distinguish between neutral and adaptive genetic variation, using the mountain pine beetle outbreak system as a motivating example. We simulated the demographic and population genetic dynamics of populations undergoing range expansion using an individual-based genetic model CDMetaPOP. We investigated how the false discovery rate of adaptive loci is affected by (i) dispersal capacity, (ii) timing of sampling, and (iii) the strength of selection on an adaptive reference locus. We found that a combination of weak dispersal, weak selection, and early sampling presents the greatest risk of misidentifying loci under selection. Expanding populations present unique challenges to the reliable identification of adaptive loci. We demonstrate that there is a need for further methodological development to account for directional demographic processes in landscape genomics.

Divergent selection along elevational gradients promotes genetic and phenotypic disparities among small mammal populations

Species distributed along mountain slopes, facing contrasting habitats in short geographic scale, are of particular interest to test how ecologically based divergent selection promotes phenotypic and genetic disparities as well as to assess isolation-by-environment mechanisms. Here, we conduct the first broad comparative study of phenotypic variation along elevational gradients, integrating a large array of ecological predictors and disentangling population genetic driver processes. The skull form of nine ecologically distinct species distributed over a large altitudinal range (100-4200 m) was compared to assess whether phenotypic divergence is a common phenomenon in small mammals and whether it shows parallel patterns. We also investigated the relative contribution of biotic (competition and predation) and abiotic parameters on phenotypic divergence via mixed models. Finally, we assessed the population genetic structure of a rodent species (Niviventer confucianus) via analysis of molecular variance and FST along three mountain slopes and tested the isolation-by-environment hypothesis using Mantel test and redundancy analysis. We found a consistent phenotypic divergence and marked genetic structure along elevational gradients; however, the species showed mixed patterns of size and skull shape trends across mountain zones. Individuals living at lower altitudes differed greatly in both phenotype and genotype from those living at high elevations, while middle-elevation individuals showed more intermediate forms. The ecological parameters associated with phenotypic divergence along elevation gradients are partly related to species' ecological and evolutionary constraints. Fossorial and solitary animals are mainly affected by climatic factors, while terrestrial and more gregarious species are influenced by biotic and abiotic parameters. A novel finding of our study is that predator richness emerged as an important factor associated with the intraspecific diversification of the mammalian skull along elevational gradients, a previously overlooked parameter. Population genetic structure was mainly driven by environmental heterogeneity along mountain slopes, with no or a week spatial effect, fitting the isolation-by-environment scenario. Our study highlights the strong and multifaceted effects of heterogeneous steep habitats and ecologically based divergent selective forces in small mammal populations.

Prey availability, prey selection, and trophic niche width in the lizard Psammodromus algirus along an elevational gradient

Mountains imply enormous environmental variation, with alpine habitats entailing harsh environments, especially for ectotherms such as lizards. This environmental variability also may imply variation in prey availability. However, little is known about how lizard trophic ecology varies with elevation. In this study, we analyze diet, prey availability, prey selection, and trophic niche width in the lacertid lizard Psammodromus algirus along a 2,200-m elevational gradient in the Sierra Nevada (SE Spain). The analysis of fecal samples has shown that Orthoptera, Formicidae, Hemiptera, Coleoptera, and Araneae are the main prey, although, according to their abundance in pitfall traps, Formicidae and Coleoptera are rejected by the lizard whereas Orthoptera, Hemiptera, and Araneae are preferred. Prey abundance and diversity increase with elevation and diet subtly varies along with the elevational gradient. The consumption of Coleoptera increases with elevation probably as a consequence of the lizard foraging more in open areas while basking. The electivity for Araneae increases with elevation. Araneae are rejected in the lowlands—where they are relatively abundant—whereas, at high elevation, this lizard positively selects them, despite they being less abundant. The lizard trophic niche width expands with elevation due to concomitant greater prey diversity and hence this lizard feeds on more prey types in highlands. Although no sex difference in diet has been found, the trophic niche is broader in females than males. As a whole, alpine lizards show a trophic niche similar to that found at lower elevations, suggesting that P. algirus is well adapted to the harsh environment found in alpine areas.

Morphological homeostasis in the fossil record

Morphological homeostasis limits the extent to which genetic and/or environmental variation is translated into phenotypic variation, providing generation-to-generation fitness advantage under a stabilizing selection regime. Depending on its lability, morphological homeostasis might also have a longer-term impact on evolution by restricting the variation-and thus the response to directional selection-of a trait. The fossil record offers an inviting opportunity to investigate whether and how morphological homeostasis constrained trait evolution in lineages or clades on long timescales (thousands to millions of years) that are not accessible to neontological studies. Fossils can also reveal insight into the nature of primitive developmental systems that might not be predictable from the study of modern organisms. The ability to study morphological homeostasis in fossils is strongly limited by taphonomic processes that can destroy, blur, or distort the original biological signal: genetic data are unavailable; phenotypic data can be modified by tectonic or compaction-related deformation; time-averaging limits temporal resolution; and environmental variation is hard to study and impossible to control. As a result of these processes, neither allelic sensitivity (and thus genetic canalization) nor macroenvironmental sensitivity (and thus environmental canalization) can be unambiguously assessed in the fossil record. However, homeorhesis-robustness against microenvironmental variation (developmental noise)-can be assessed in ancient developmental systems by measuring the level of fluctuating asymmetry (FA) in a nominally symmetric trait. This requires the analysis of multiple, minimally time-averaged samples of exquisite preservational quality. Studies of FA in fossils stand to make valuable contributions to our understanding of the deep-time significance of homeorhesis. Few empirical studies have been conducted to date, and future paleontological research focusing on how homeorhesis relates to evolutionary rate (including stasis), species survivorship, and purported macroevolutionary trends in evolvability would reap high reward.

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta-analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.

Opposed elevational variation in prevalence and intensity of endoparasites and their vectors in a lizard

Studying the causes of parasite geographic distribution is relevant to understand ecological and evolutionary processes that affect host populations as well as for species conservation. Temperature is one of the most important environmental variables affecting parasite distribution, as raising temperatures positively affect development, reproduction, and rate of transmission of both endo- and ectoparasites. In this context, it is generally accepted that, in mountains, parasite abundance decreases with elevation. However, empirical evidence on this topic is limited. In the present study, we analyzed the elevational variation of hemoparasites and ectoparasites of a lizard, Psammodromus algirus, along a 2,200-m elevational gradient in Sierra Nevada (SE Spain). As predicted, ectoparasite (mites, ticks, mosquitoes, and sandflies) abundance decreased with elevation. However, hemoparasite prevalence and intensity in the lizard augmented with altitude, showing a pattern contrary to their vectors (mites). We suggest that tolerance to hemoparasites may increase with elevation as a consequence of lizards at high altitudes taking advantage of increased body condition and food availability, and reduced oxidative stress. Moreover, lizards could have been selected for higher resistance against hemoparasites at lowlands (where higher rates of replication are expected), thus reducing hemoparasite prevalence and load. Our findings imply that, in a scenario of climate warming, populations of lizards at high elevation may face increased abundance of ectoparasites, accompanied with strong negative effects.

Temperature-driven selection on metabolic traits increases the strength of an algal-grazer interaction in naturally warmed streams

Trophic interactions are important determinants of the structure and functioning of ecosystems. Because the metabolism and consumption rates of ectotherms increase sharply with temperature, there are major concerns that global warming will increase the strength of trophic interactions, destabilizing food webs, and altering ecosystem structure and function. We used geothermally warmed streams that span an 11°C temperature gradient to investigate the interplay between temperature-driven selection on traits related to metabolism and resource acquisition, and the interaction strength between the keystone gastropod grazer, Radix balthica, and a common algal resource. Populations from a warm stream (~28°C) had higher maximal metabolic rates and optimal temperatures than their counterparts from a cold stream (~17°C). We found that metabolic rates of the population originating from the warmer stream were higher across all measurement temperatures. A reciprocal transplant experiment demonstrated that the interaction strengths between the grazer and its algal resource were highest for both populations when transplanted into the warm stream. In line with the thermal dependence of respiration, interaction strengths involving grazers from the warm stream were always higher than those with grazers from the cold stream. These results imply that increases in metabolism and resource consumption mediated by the direct, thermodynamic effects of higher temperatures on physiological rates are not mitigated by metabolic compensation in the long term, and suggest that warming could increase the strength of algal-grazer interactions with likely knock-on effects for the biodiversity and productivity of aquatic ecosystems.