Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

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.     

Lower thermal tolerance in nocturnal than in diurnal ants: a challenge for nocturnal ectotherms facing global warming

1. The thermal adaptation hypothesis proposes that because thermoregulation involves a high metabolic cost, thermal limits of organisms must be locally adapted to temperatures experienced in their environments. There is evidence that tolerance to high temperatures decreases in insects inhabiting colder habitats and microclimates. However, it is not clear if thermal limits of ectotherms with contrasting temporal regimes, such as diurnal and nocturnal insects, are also adapted to temperatures associated with their circadian activities. 2. This study explores differences in heat tolerance among diurnal and nocturnal ant species in four ecosystems in Mexico: tropical montane, tropical rainforest, subtropical dry forests, and high‐elevation semi‐desert. 3. The critical thermal maximum (CT_max), i.e. the temperature at which ants lost motor control, was estimated for diurnal and nocturnal species. CT_max for 19 diurnal and 12 nocturnal ant species distributed among 45 populations was also estimated. 4. Semi‐desert and subtropical dry forest ants displayed higher tolerances to high temperatures than did ants in tropical rainforest. The lowest tolerance to high temperatures was recorded in tropical montane forest ants. In general, among all habitats, the CT_max of nocturnal ants was lower than that of diurnal ants. 5. An increase in nocturnal temperatures, combined with lower tolerance to high temperatures, may represent a substantial challenge for nocturnal ectotherms in a warming world.     

Ontogenetic thermal tolerance and performance of ectotherms at variable temperatures

Early experience and environmental conditions during ontogeny may affect organismal structure, physiology and fitness. Here, we assessed the effect of developmental acclimation to environmental thermal variability on walking speed in Drosophila melanogaster adults. Our results showed a shift in the performance curve to the right. Thus, upper and lower thermal limits exhibited developmental plasticity. Additionally, in constant and variable climatic scenarios, flies shifted to the right the optimum temperature but the maximum performance decreased only in flies reared on high temperatures and high thermal variability. Overall, we showed that environmental cues during ontogeny might help to construct phenotypic variation, which supports the hypothesis of ontogenetic dependence of thermal tolerances.     

Rapid within‐ and transgenerational changes in thermal tolerance and fitness in variable thermal landscapes

Phenotypic plasticity may increase the performance and fitness and allow organisms to cope with variable environmental conditions. We studied within‐generation plasticity and transgenerational effects of thermal conditions on temperature tolerance and demographic parameters in Drosophila melanogaster. We employed a fully factorial design, in which both parental (P) and offspring generations (F1) were reared in a constant or a variable thermal environment. Thermal variability during ontogeny increased heat tolerance in P, but with demographic cost as this treatment resulted in substantially lower survival, fecundity, and net reproductive rate. The adverse effects of thermal variability (V) on demographic parameters were less drastic in flies from the F1, which exhibited higher net reproductive rates than their parents. These compensatory responses could not totally overcome the challenges of the thermally variable regime, contrasting with the offspring of flies raised in a constant temperature (C) that showed no reduction in fitness with thermal variation. Thus, the parental thermal environment had effects on thermal tolerance and demographic parameters in fruit fly. These results demonstrate how transgenerational effects of environmental conditions on heat tolerance, as well as their potential costs on other fitness components, can have a major impact on populations’ resilience to warming temperatures and more frequent thermal extremes.     

Fluctuating thermal environments and time‐dependent effects on fruit fly egg‐hatching performance

Organismal performance in a changing environment is dependent on temporal patterns and duration of exposure to thermal variability. We experimentally assessed the time‐dependent effects of thermal variability (i.e., patterns of thermal exposure) on the hatching performance of Drosophila melanogaster. Flies were collected in central Chile and maintained for four generations in laboratory conditions. Fourth generation eggs were acclimated to different thermal fluctuation cycles until hatching occurred. Our results show that the frequency of extreme thermal events has a significant effect on hatching success. Eggs exposed to 24 hr cycles of thermal fluctuation had a higher proportion of eggs that hatched than those acclimated to shorter (6 and 12 hr) and longer cycles (48 hr). Furthermore, eggs subjected to frequent thermal fluctuations hatched earlier than those acclimated to less frequent thermal fluctuations. Overall, we show that, egg‐to‐adult viability is dependent on the pattern of thermal fluctuations experienced during ontogeny; thus, the pattern of thermal fluctuation experienced by flies has a significant and until now unappreciated impact on fitness.     

Grow with the flow: a latitudinal cline in physiology is associated with more variable precipitation in Erythranthe cardinalis

Local adaptation is commonly observed in nature: organisms perform well in their natal environment, but poorly outside it. Correlations between traits and latitude, or latitudinal clines, are among the most common pieces of evidence for local adaptation, but identifying the traits under selection and the selective agents is challenging. Here, we investigated a latitudinal cline in growth and photosynthesis across 16 populations of the perennial herb Erythranthe cardinalis (Phrymaceae). Using machine learning methods, we identify interannual variation in precipitation as a likely selective agent: southern populations from more variable environments had higher photosynthetic rates and grew faster. We hypothesize that selection may favour a more annualized life history – grow now rather than save for next year – in environments where severe droughts occur more often. Thus, our study provides insight into how species may adapt if Mediterranean climates become more variable due to climate change.     

The genetic architecture of defence as resistance to and tolerance of bacterial infection in Drosophila melanogaster

Defence against pathogenic infection can take two forms: resistance and tolerance. Resistance is the ability of the host to limit a pathogen burden, whereas tolerance is the ability to limit the negative consequences of infection at a given level of infection intensity. Evolutionarily, a tolerance strategy that is independent of resistance could allow the host to avoid mounting a costly immune response and, theoretically, to avoid a co‐evolutionary arms race between pathogen virulence and host resistance. Biomedically, understanding the mechanisms of tolerance and how they relate to resistance could potentially yield treatment strategies that focus on health improvement instead of pathogen elimination. To understand the impact of tolerance on host defence and identify genetic variants that determine host tolerance, we defined genetic variation in tolerance as the residual deviation from a binomial regression of fitness under infection against infection intensity. We then performed a genomewide association study to map the genetic basis of variation in resistance to and tolerance of infection by the bacterium Providencia rettgeri. We found a positive genetic correlation between resistance and tolerance, and we demonstrated that the level of resistance is highly predictive of tolerance. We identified 30 loci that predict tolerance, many of which are in genes involved in the regulation of immunity and metabolism. We used RNAi to confirm that a subset of mapped genes have a role in defence, including putative wound repair genes grainy head and debris buster. Our results indicate that tolerance is not an independent strategy from resistance, but that defence arises from a collection of physiological processes intertwined with canonical immunity and resistance.     

Local thermal adaptation and local temperature regimes drive the performance of a parasitic helminth under climate change: The case of Marshallagia marshalli from wild ungulates

Across a species' range, populations are exposed to their local thermal environments, which on an evolutionary scale, may cause adaptative differences among populations. Helminths often have broad geographic ranges and temperature-sensitive life stages but little is known about whether and how local thermal adaptation can influence their response to climate change. We studied the thermal responses of the free-living stages of Marshallagia marshalli, a parasitic nematode of wild ungulates, along a latitudinal gradient. We first determine its distribution in wild sheep species in North America. Then we cultured M. marshalli eggs from different locations at temperatures from 5 to 38°C. We fit performance curves based on the metabolic theory of ecology to determine whether development and mortality showed evidence of local thermal adaptation. We used parameter estimates in life-cycle-based host–parasite models to understand how local thermal responses may influence parasite performance under general and location-specific climate-change projections. We found that M. marshalli has a wide latitudinal and host range, infecting wild sheep species from New Mexico to Yukon. Increases in mortality and development time at higher temperatures were most evident for isolates from northern locations. Accounting for location-specific parasite parameters primarily influenced the magnitude of climate change parasite performance, while accounting for location-specific climates primarily influenced the phenology of parasite performance. Despite differences in development and mortality among M. marshalli populations, when using site-specific climate change projections, there was a similar magnitude of impact on the relative performance of M. marshalli among populations. Climate change is predicted to decrease the expected lifetime reproductive output of M. marshalli in all populations while delaying its seasonal peak by approximately 1 month. Our research suggests that accurate projections of the impacts of climate change on broadly distributed species need to consider local adaptations of organisms together with local temperature profiles and climate projections.     

Interspecific variation and elevated CO2 influence the relationship between plant chemical resistance and regrowth tolerance

To understand how comprehensive plant defense phenotypes will respond to global change, we investigated the legacy effects of elevated CO₂ on the relationships between chemical resistance (constitutive and induced via mechanical damage) and regrowth tolerance in four milkweed species (Asclepias). We quantified potential resistance and tolerance trade‐offs at the physiological level following simulated mowing, which are relevant to milkweed ecology and conservation. We examined the legacy effects of elevated CO₂ on four hypothesized trade‐offs between the following: (a) plant growth rate and constitutive chemical resistance (foliar cardenolide concentrations), (b) plant growth rate and mechanically induced chemical resistance, (c) constitutive resistance and regrowth tolerance, and (d) regrowth tolerance and mechanically induced resistance. We observed support for one trade‐off between plant regrowth tolerance and mechanically induced resistance traits that was, surprisingly, independent of CO₂ exposure. Across milkweed species, mechanically induced resistance increased by 28% in those plants previously exposed to elevated CO_2. In contrast, constitutive resistance and the diversity of mechanically induced chemical resistance traits declined in response to elevated CO₂ in two out of four milkweed species. Finally, previous exposure to elevated CO₂ uncoupled the positive relationship between plant growth rate and regrowth tolerance following damage. Our data highlight the complex and dynamic nature of plant defense phenotypes under environmental change and question the generality of physiologically based defense trade‐offs.     

Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradient

The European gypsy moth (Lymantria dispar L.) was first introduced to Massachusetts in 1869 and within 150 years has spread throughout eastern North America. This large‐scale invasion across a heterogeneous landscape allows examination of the genetic signatures of adaptation potentially associated with rapid geographical spread. We tested the hypothesis that spatially divergent natural selection has driven observed changes in three developmental traits that were measured in a common garden for 165 adult moths sampled from six populations across a latitudinal gradient covering the entirety of the range. We generated genotype data for 91,468 single nucleotide polymorphisms based on double digest restriction‐site associated DNA sequencing and used these data to discover genome‐wide associations for each trait, as well as to test for signatures of selection on the discovered architectures. Genetic structure across the introduced range of gypsy moth was low in magnitude (F_ST = 0.069), with signatures of bottlenecks and spatial expansion apparent in the rare portion of the allele frequency spectrum. Results from applications of Bayesian sparse linear mixed models were consistent with the presumed polygenic architectures of each trait. Further analyses indicated spatially divergent natural selection acting on larval development time and pupal mass, with the linkage disequilibrium component of this test acting as the main driver of observed patterns. The populations most important for these signals were two range‐edge populations established less than 30 generations ago. We discuss the importance of rapid polygenic adaptation to the ability of non‐native species to invade novel environments.     

Sap‐feeding insects on forest trees along latitudinal gradients in northern Europe: a climate‐driven patterns

Knowledge of the latitudinal patterns in biotic interactions, and especially in herbivory, is crucial for understanding the mechanisms that govern ecosystem functioning and for predicting their responses to climate change. We used sap‐feeding insects as a model group to test the hypotheses that the strength of plant–herbivore interactions in boreal forests decreases with latitude and that this latitudinal pattern is driven primarily by midsummer temperatures. We used a replicated sampling design and quantitatively collected and identified all sap‐feeding insects from four species of forest trees along five latitudinal gradients (750–1300 km in length, ten sites in each gradient) in northern Europe (59 to 70°N and 10 to 60°E) during 2008–2011. Similar decreases in diversity of sap‐feeding insects with latitude were observed in all gradients during all study years. The sap‐feeder load (i.e. insect biomass per unit of foliar biomass) decreased with latitude in typical summers, but increased in an exceptionally hot summer and was independent of latitude during a warm summer. Analysis of combined data from all sites and years revealed dome‐shaped relationships between the loads of sap‐feeders and midsummer temperatures, peaking at 17 °C in Picea abies, at 19.5 °C in Pinus sylvestris and Betula pubescens and at 22 °C in B. pendula. From these relationships, we predict that the losses of forest trees to sap‐feeders will increase by 0–45% of the current level in southern boreal forests and by 65–210% in subarctic forests with a 1 °C increase in summer temperatures. The observed relationships between temperatures and the loads of sap‐feeders differ between the coniferous and deciduous tree species. We conclude that climate warming will not only increase plant losses to sap‐feeding insects, especially in subarctic forests, but can also alter plant‐plant interactions, thereby affecting both the productivity and the structure of future forest ecosystems.     

Increased temperatures negatively affect Juniperus communis seeds: evidence from transplant experiments along a latitudinal gradient

With a distribution range that covers most of the Northern hemisphere, common juniper (Juniperus communis) has one of the largest ranges of all vascular plant species. In several regions in Europe, however, populations are decreasing in size and number due to failing recruitment. One of the main causes for this failure is low seed viability. Observational evidence suggests that this is partly induced by climate warming, but our mechanistic understanding of this effect remains incomplete. Here, we experimentally assess the influence of temperature on two key developmental phases during sexual reproduction, i.e. gametogenesis and fertilisation (seed phase two, SP2) and embryo development (seed phase three, SP3). Along a latitudinal gradient from southern France to central Sweden, we installed a transplant experiment with shrubs originating from Belgium, a region with unusually low juniper seed viability. Seeds of both seed phases were sampled during three consecutive years, and seed viability assessed. Warming temperatures negatively affected the seed viability of both SP2 and SP3 seeds along the latitudinal gradient. Interestingly, the effect on embryo development (SP3) only occurred in the third year, i.e. when the gametogenesis and fertilisation also took place in warmer conditions. We found strong indications that this negative influence mostly acts via disrupting growth of the pollen tube, the development of the female gametophyte and fertilisation (SP2). This, in turn, can lead to failing embryo development, for example, due to nutritional problems. Our results confirm that climate warming can negatively affect seed viability of juniper.     

Some (worms) like it hot: fish parasites grow faster in warmer water,and alter host thermal preferences

Elevated environmental temperatures associated with anthropogenic warming have the potential to impact host‐parasite interactions, with consequences for population health and ecosystem functioning. One way that elevated temperatures might influence parasite prevalence and intensity is by increasing life cycle completion rates. Here, we investigate how elevated temperatures impact a critical phase of the life cycle of the bird tapeworm Schistocephalus solidus – the growth of plerocercoid larvae in host fish (three‐spined sticklebacks Gasterosteus aculeatus). By 8 weeks post‐infection, plerocercoids recovered from experimentally infected sticklebacks held at 20 °C weighed on average 104.9 mg, with all exceeding 50 mg, the mass considered consistently infective to definitive hosts. In contrast, plerocercoids from sticklebacks held at 15 °C weighed on average 26.5 mg, with none exceeding 50 mg. As small increases in plerocercoid mass affect adult fecundity disproportionately in this species, enhanced plerocercoid growth at higher temperatures predicts dramatically increased output of infective parasite stages. Subsequent screening of thermal preferences of sticklebacks from a population with endemic S. solidus infection demonstrated that fish harbouring infective plerocercoids show significant preferences for warmer temperatures. Our results therefore indicate that parasite transmission might be affected in at least two ways under anthropogenic warming; by enhancing rates of parasite growth and development, and by increasing the likelihood of hosts being able to seek out proliferating warmer microhabitats. Furthermore, our results suggest the potential for positive feedback between parasite growth and host thermal preferences, which could dramatically increase the effects of even small temperature increases. We discuss the possible mechanisms underpinning our results, their likely ecological consequences and highlight key areas for further research.     

The evolution of costly acquired immune memory

A key feature of the vertebrate adaptive immune system is acquired immune memory, whereby hosts launch a faster and heightened response when challenged by previously encountered pathogens, preventing full infection. Here, we use a mathematical model to explore the role of ecological and epidemiological processes in shaping selection for costly acquired immune memory. Applying the framework of adaptive dynamics to the classic SIR (Susceptible‐Infected‐Recovered) epidemiological model, we focus on the conditions that may lead hosts to evolve high levels of immunity. Linking our work to previous theory, we show how investment in immune memory may be greatest at long or intermediate host lifespans depending on whether immunity is long lasting. High initial costs to gain immunity are also found to be essential for a highly effective immune memory. We also find that high disease infectivity and sterility, but intermediate virulence and immune period, increase selection for immunity. Diversity in host populations through evolutionary branching is found to be possible but only for a limited range of parameter space. Our model suggests that specific ecological and epidemiological conditions have to be met for acquired immune memory to evolve.     

Cryptic species of a water hyacinth biological control agent revealed in South Africa: host specificity,impact, and thermal tolerance

The discovery that cryptic species are more abundant than previously thought has implications for weed biological control, as there is a risk that cryptic species may be inadvertently released with consequences for the safety of the practice. A cryptic species of a biological control agent released for the control of the invasive alien macrophyte, water hyacinth, Eichhornia crassipes (C. Mart.) Solms. (Pontederiaceae), was recently discovered in South Africa. The two species were considered a single species prior to genetic analysis and interbreeding experiments. The original biological control agent retains the name Eccritotarsus catarinensis (Carvalho) (Heteroptera: Miridae) whereas the new species has been described as Eccritotarsus eichhorniae Henry. In this study, we compared the host specificity, efficacy, and thermal physiologies of the two species. The host specificity of the two species within the Pontederiaceae was very similar and both are safe for release in South Africa. Comparison of the per capita impact of the two species indicated that E. eichhorniae was the more damaging species but this is likely to be influenced by temperature, with E. catarinensis being more effective under lower temperatures and E. eichhorniae being more effective under higher temperatures. Releasing the correct species for the thermal environment of each release site will improve the level of control of water hyacinth in South Africa. This example highlights the need to keep populations of biological control agents from different native range collection localities separate, and to screen for host specificity and efficacy.     

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.     

The influence of bat ecology on viral diversity and reservoir status

Repeated emergence of zoonotic viruses from bat reservoirs into human populations demands predictive approaches to preemptively identify virus‐carrying bat species. Here, we use machine learning to examine drivers of viral diversity in bats, determine whether those drivers depend on viral genome type, and predict undetected viral carriers. Our results indicate that bat species with longer life spans, broad geographic distributions in the eastern hemisphere, and large group sizes carry more viruses overall. Life span was a stronger predictor of deoxyribonucleic acid viral diversity, while group size and family were more important for predicting ribonucleic acid viruses, potentially reflecting broad differences in infection duration. Importantly, our models predict 54 bat species as likely carriers of zoonotic viruses, despite not currently being considered reservoirs. Mapping these predictions as a proportion of local bat diversity, we identify global regions where efforts to reduce disease spillover into humans by identifying viral carriers may be most productive.     

The thermal mismatch hypothesis explains host susceptibility to an emerging infectious disease

Parasites typically have broader thermal limits than hosts, so large performance gaps between pathogens and their cold‐ and warm‐adapted hosts should occur at relatively warm and cold temperatures, respectively. We tested this thermal mismatch hypothesis by quantifying the temperature‐dependent susceptibility of cold‐ and warm‐adapted amphibian species to the fungal pathogen Batrachochytrium dendrobatidis (Bd) using laboratory experiments and field prevalence estimates from 15 410 individuals in 598 populations. In both the laboratory and field, we found that the greatest susceptibility of cold‐ and warm‐adapted hosts occurred at relatively warm and cool temperatures, respectively, providing support for the thermal mismatch hypothesis. Our results suggest that as climate change shifts hosts away from their optimal temperatures, the probability of increased host susceptibility to infectious disease might increase, but the effect will depend on the host species and the direction of the climate shift. Our findings help explain the tremendous variation in species responses to Bd across climates and spatial, temporal and species‐level variation in disease outbreaks associated with extreme weather events that are becoming more common with climate change.