How will mosquitoes adapt to climate warming?
Abstract
The potential for adaptive evolution to enable species persistence under a changing climate is one of the most important questions for understanding impacts of future climate change. Climate adaptation may be particularly likely for short-lived ectotherms, including many pest, pathogen, and vector species. For these taxa, estimating climate adaptive potential is critical for accurate predictive modeling and public health preparedness. Here, we demonstrate how a simple theoretical framework used in conservation biology—evolutionary rescue models—can be used to investigate the potential for climate adaptation in these taxa, using mosquito thermal adaptation as a focal case. Synthesizing current evidence, we find that short mosquito generation times, high population growth rates, and strong temperature-imposed selection favor thermal adaptation. However, knowledge gaps about the extent of phenotypic and genotypic variation in thermal tolerance within mosquito populations, the environmental sensitivity of selection, and the role of phenotypic plasticity constrain our ability to make more precise estimates. We describe how common garden and selection experiments can be used to fill these data gaps. Lastly, we investigate the consequences of mosquito climate adaptation on disease transmission using Aedes aegypti-transmitted dengue virus in Northern Brazil as a case study. The approach outlined here can be applied to any disease vector or pest species and type of environmental change.
Introduction
Climate change is expected to have major impacts on species distributions in coming decades, and predicting these impacts is an area of intense research interest. As their basic physiological and ecological traits depend heavily on temperature, climate impacts are expected to be particularly strong for ectotherms (Deutsch et al., 2008). This encompasses many taxa that threaten human health and well-being, including agricultural and forest pests, wildlife and plant pathogens, and disease vectors, for which accurately predicting distributions under climate change is critical for protecting human and animal health. Several prominent reviews have found that climate change is expected to increase, decrease, or, most commonly, shift the distributions of these taxa due to nonlinear and interactive effects of temperature and other climatic factors (Porter et al., 1991; Harvell et al., 2002; Deutsch et al., 2008; Rohr et al., 2011; Altizer et al., 2013; Lafferty and Mordecai, 2016; Pinsky et al., 2019; Lehmann et al., 2020; Rohr and Cohen, 2020). However, these predictions largely assume that species climate responses are fixed, ignoring the potential for adaptive responses.
Evidence of evolutionary adaptation to contemporary climate change has emerged for diverse taxa including mammals (Réale et al., 2003), fish (Kovach et al., 2012), plants (Franks et al., 2007; Exposito-Alonso et al., 2018), birds (Nussey et al., 2005; Karell et al., 2011), reptiles (Logan et al., 2018), and insects (Umina et al., 2005). However, while climate adaptation has typically been studied in the context of conservation biology, population genetics theory suggests that evolutionary adaptation is most likely for short-lived species with high population growth rates—properties of many pest, pathogen, and vector species (Lynch and Lande, 1993; Bürger and Lynch, 1995; Kingsolver, 2009). For several of these species, recent research demonstrates the potential for climate adaptation within a few decades. For example, in the European gypsy moth (Lymantria dispar)—one of the world’s most destructive forest pests (Montgomery and Wallner, 1988)—shifts in thermal tolerance were evident within 30 years of population expansion (Friedline et al., 2019). Similarly, in the Asian tiger mosquito (Aedes albopictus)—a vector of yellow fever, dengue, and chikungunya viruses—adaptive responses to novel temperature conditions were detected within 10–30 years of population expansions (Medley, 2010; Urbanski et al., 2012; Egizi et al., 2015; Medley et al., 2019). Despite this emerging evidence of the potential for rapid climate adaptation in pest and vector taxa, a general strategy for understanding and estimating their climate adaptive potential is lacking.
The potential for climate adaptation is a particularly important open question in mosquito-borne disease biology (Mordecai et al., 2019; Franklinos et al., 2019). Mosquito-borne diseases are a major public health burden, causing an estimated 500 million cases and millions of deaths globally each year (World Health Organization, 2014; World Health Organization, 2018). Environmental drivers of mosquito-borne disease transmission have been relatively well studied, and consistently highlight temperature—and by extension, climate warming—as a fundamental driver (Shragai et al., 2017; Mordecai et al., 2019; Franklinos et al., 2019; Shocket et al., 2021). Temperature influences mosquito-borne disease dynamics because it directly affects mosquito physiology, life cycles, behavior, and competence for disease transmission (Cator et al., 2020). For mosquitoes and other ectotherms, temperature has strong, nonlinear effects on traits such as survival and fecundity that lead to unimodal effects of temperature on fitness, where temperatures above and below intermediate thermal optima limit mosquito population growth (Huey and Stevenson, 1979; Huey and Berrigan, 2001; Angilletta, 2009; Amarasekare and Savage, 2012; Mordecai et al., 2019). Recent forecasts based on the unimodal relationship with temperature predict that in some areas where disease risk is currently high, future warming will decrease transmission risk as temperatures exceed mosquito thermal optima and limits (Gething et al., 2010; Ryan et al., 2015; Ryan et al., 2019; Mordecai et al., 2019; Mordecai et al., 2020). However, these predictions are likely to underestimate future disease risk if mosquitoes adapt to climate warming. As a result, estimating the potential for mosquito thermal adaptation is critical for more accurate predictive modeling.
Here, we outline a theoretical framework for investigating climate adaptive potential in ectotherm pests and pathogens, which draws from evolutionary rescue models typically used in conservation biology. We use mosquito adaptation to warming temperature as a focal case given the high global health burden of mosquito-borne disease and substantial recent research progress on mosquito thermal biology. In the following sections, we: (1) outline the theoretical framework and the specific parameters needed to estimate adaptive potential, (2) synthesize available information for mosquito thermal adaptation and identify key data gaps for predictive modeling, (3) highlight priorities and describe specific empirical approaches for filling these gaps, (4) explore the consequences of mosquito thermal adaptation on disease transmission, and (5) discuss the application of this framework to other vector and pest species. We focus here on adaptation to warming temperatures that exceed current mosquito thermal optima. We consider temperature in isolation despite its influence on relative humidity, which has strong impacts on mosquito population dynamics and host-seeking behavior and is predicted to shift with anthropogenic climate change. Several studies have investigated mosquito adaptation to desiccation (e.g. Kearney et al., 2009; Simard et al., 2009; Fouet et al., 2012), but few have investigated mosquito responses to simultaneous variation in temperature and humidity (but see Yamana and Eltahir, 2013; Yamana et al., 2016). Given this lack of empirical data, we focus on temperature specifically, but we discuss the inclusion of aridity in predictive modeling approaches (see ‘Climate factors currently limiting population persistence’). Similarly, we do not include non-climate factors such as biotic interactions, land use change, and human activities that may impact mosquito population persistence, because their effects on mosquito responses to temperature and thermal adaptive potential remain poorly understood. We discuss the adaptive potential of mosquitoes broadly, but our principal interest is in populations of major disease vector species (e.g. Aedes aegypti, Ae. albopictus, Anopheles gambiae, Culex pipiens, Cx. quinquefasciatus, which transmit dengue, chikungunya, Zika, and West Nile viruses, malaria, and other pathogens) and we discuss species-specific responses where possible. After presenting mosquito thermal adaptation as a focal case, we discuss how the approach we describe here can be applied to study the adaptive potential of any species in response to any specific environmental change.
Framework for investigating climate adaptation
Species may respond to warming temperatures through three primary mechanisms: tracking suitable temperatures through range shifts, avoiding or temporarily coping with stressful temperatures through phenotypic plasticity (e.g. shifting biting activity to cooler times), and tolerating warming through genetic evolutionary adaptation (e.g. evolved shifts in thermal tolerance resulting from selection). Here, we focus on evolutionary adaptation as it would enable in situ population persistence under sustained environmental change and is currently the least well-understood climate response (Merilä and Hendry, 2014; Urban et al., 2016; González-Tokman et al., 2020). Investigating the potential for evolutionary climate adaptation requires identifying: (1) the climate factors currently limiting population persistence, (2) the most climate-sensitive and fitness-relevant traits, and (3) the potential evolutionary rates of these traits (Figure 1). We describe these factors further below, using mosquito thermal adaptation as a focal case.
Climate factors currently limiting population persistence
Temperature fundamentally limits mosquito ranges and persistence through its influence on mosquito survival, development, and reproductive rates. However, the precise aspects of temperature that determine these limits remain unclear (Christophers, 1960; Yang et al., 2009; Brady et al., 2013; Shapiro et al., 2017; Tesla et al., 2018; Shocket et al., 2018; Mordecai et al., 2019). Temperature averages, variability, extremes, and interactions among these factors all impact ectotherm fitness (Nene et al., 2007; Lambrechts et al., 2011; Bozinovic et al., 2011; Kingsolver et al., 2013; Paaijmans et al., 2013; Blanford et al., 2013; Dowd et al., 2015; Buckley and Huey, 2016; Ma et al., 2021). The temperature variable that most strongly constrains persistence for a particular mosquito population will likely vary based on breeding and resting habitat preferences and by location, as thermal regimes can differ greatly between water sources used for oviposition, between indoor and outdoor environments, and by latitude and altitude (Paaijmans et al., 2008; Paaijmans et al., 2010; Paaijmans and Thomas, 2011). For example, changes in diurnal temperature variation may impose strong selection on East African An. arabiensis populations that primarily rest outdoors (Githeko et al., 1996), but have less impact on An. gambiae s.s and An. funestus populations that primarily rest indoors, where temperatures vary less drastically (Bødker et al., 2003; Minakawa et al., 2006). Similarly, endemic populations of a given species may experience vastly different climate regimes than invasive populations, leading to different constraints on their persistence (Moran and Alexander, 2014; Egizi et al., 2015). For example, Aedes japonicus japonicus populations in their native range in East Asia must contend with extreme cold winter temperatures, which they survive as diapausing eggs, while invasive populations in Hawaii and the southeastern U.S. experience warmer year-round temperatures (Reeves and Korecki, 2004; Larish and Savage, 2005; Kaufman and Fonseca, 2014). Based on evidence from other arthropod species, maximum annual temperature has been identified as the strongest driver of species extinctions (out of the 19 WorldClim Bioclimatic variables; Román-Palacios and Wiens, 2020), thus changes in maximum temperatures may exert the strongest selection pressure on mosquito populations near their warm range limits. As most temperature – trait responses are studied under constant temperatures (Deutsch et al., 2008; Angilletta, 2009; Paaijmans et al., 2013; Vasseur et al., 2014; Dowd et al., 2015; Buckley and Huey, 2016), we consider mosquito responses to increases in mean temperature as a focal example. Similarly, we consider the independent effects of temperature on mosquitoes despite its influence on moisture availability. Restrictions on moisture availability can readily be incorporated into predictive models by setting thresholds for annual precipitation or vegetative greenness (an indicator of local moisture availability; Suzuki et al., 2006), and masking out regions falling below these levels. For example, in estimating future temperature-based suitability of malaria transmission, Ryan et al., 2015 applied a threshold for Normalized Difference Vegetation Index (NDVI)—two consecutive months of NDVI above 0.125—and masked out all regions falling below this limit. In general, the framework we present can readily be applied to any specific measure of temperature or other environmental variable, such as temperature extremes, precipitation, wind patterns, land use change, and human activities (Reiter, 2001; Patz et al., 2008; Paaijmans and Thomas, 2011; Mordecai et al., 2019; Franklinos et al., 2019; Rocklöv and Dubrow, 2020).
Mosquito traits
As selection acts on specific phenotypes, ‘adaptation’ here refers to a change in the thermal tolerance of a specific mosquito trait—an emergent property that reflects underlying physiological changes (e.g. changes in protein thermal sensitivity; Somero, 1995, Somero, 2003, Somero, 2010, González-Tokman et al., 2020). The critical traits to examine are those with the strongest climate sensitivity and the strongest impact on overall mosquito fitness. We hypothesize that mosquito life history traits with the lowest critical thermal maxima will experience the strongest thermal selection as environmental temperatures that exceed this limit have especially strong negative impacts on ectotherm fitness (Deutsch et al., 2008; Kingsolver et al., 2013). For several mosquito vector species, these traits are adult lifespan and fecundity, suggesting thermal selection may be strongest on the adult life stage (Figure 1; reviewed in Mordecai et al., 2019). However, whether these life history traits also pose the strongest constraints on mosquito fitness and persistence at high temperatures remains poorly understood because prior work has largely focused on mosquito traits related to disease transmission during the activity season, but there may be additional traits that help mosquitoes tolerate climate extremes (e.g. diapause, aestivation). For this reason, we consider several fitness-relevant mosquito life history traits (e.g. survival, development rates, and fecundity). Additionally, we consider various metrics that describe trait thermal tolerance (e.g. time to partial paralysis, known as ‘knockdown time’ at high temperatures; trait performances at high temperatures; and temperatures causing 50% sample mortality) because they may provide differing information on species adaptive potential (Hangartner and Hoffmann, 2016). The framework outlined here can be applied to any specific trait and measurement.
Potential evolutionary rates of climate-sensitive traits
After identifying the climate factors and mosquito traits that limit population persistence, we can now compare their rates of change to predict whether populations can adapt apace with environmental change (Figure 1). To do so, we turn to evolutionary rescue models, which estimate the maximum rate of evolutionary change (i.e. adaptive genetic turnover) of a population and compare it to the projected rate of environmental change. Populations can persist only when their maximum sustainable evolutionary rate exceeds the required rate of evolution dictated by the environment (Bell and Gonzalez, 2009; Hoffmann and Sgrò, 2011; Gomulkiewicz and Shaw, 2013; Gonzalez et al., 2013; Carlson et al., 2014; Bell, 2017). Evolutionary rescue models explicitly model demographic rates and assume that populations are comprised of different genotypes with different reproductive advantages. As these models track population responses to sustained, directional environmental change, they are well-suited to estimating the potential for thermal adaptation in response to climate warming (Huey and Kingsolver, 1993; Bürger and Lynch, 1995; Chevin et al., 2010; Bay et al., 2017; Cotto et al., 2017; Diniz‐Filho et al., 2019), and have provided valuable estimates of climate adaptation potential across a range of taxa (Gienapp et al., 2013; Bush et al., 2016; Cotto et al., 2017; Diniz‐Filho et al., 2019). Even with incomplete or imprecise knowledge of all parameters, these models can place bounds on the climate response space to indicate where adaptation is highly unlikely and to inform future data collection efforts.
Here, we consider the analytic, quantitative-genetic evolutionary rescue model described by Chevin et al., 2010. This model estimates population adaptive potential under climate warming using (Box 1; Figure 1): (1) the maximum population growth rate under optimal conditions ), (2) the population generation time (T), (3) the phenotypic variance in the trait of interest , (4) the strength of selection imposed by temperature change , (5) the trait heritability (, (6) the degree of phenotypic plasticity in thermal tolerance (b), (7) how the trait optimum changes with temperature (i.e. environmental sensitivity of selection; B), and (8) the expected rate of temperature change during the time period . Although the simplicity of this analytic evolutionary rescue model may constrain the accuracy of its projections, it illustrates the basic factors likely to affect population persistence, which we consider to be the minimum information needed to make initial predictions (see Appendix 1 for additional unmodeled factors and the ‘Priorities and approaches’ section for methods to incorporate additional complexity). We present the main findings from the available information below, including information from the closely related model organism Drosophila when little information is available for mosquitoes.
Box 1.
Evolutionary rescue model formula (Chevin et al., 2010) and parameter definitions.
T: population generation time: (for populations with discrete, non-overlapping generations), the mean time between reproduction in one cohort to reproduction in the successive cohort.
: maximum population growth rate: the intrinsic rate of increase under optimal conditions (i.e. no intra- or inter-specific competition).
: phenotypic variance: the measured variance in the trait of interest.
: heritability: the proportion of phenotypic variance in a trait attributable to additive genetic effects.
: strength of selection: the impact on fitness from deviations from the optimal trait value under a given environment. As in Kearney et al., 2009, a standardized version of selection strength can be approximated from temperature-dependent survival rates by:
where s is the percentage survival under a given environmental change, and i is given as the change in phenotype (in standard deviations) between the starting and selected populations (Falconer and Mackay, 1996; Matsumura et al., 2012).
b: phenotypic plasticity: the ability of individual genotypes to produce alternative phenotypes in different environments (Via et al., 1995). Here, plasticity encompasses thermal acclimation, dormancy and behavioral thermoregulation including shifts in mosquito biting, microhabitat usage, and oviposition sites and timing.
B: environmental sensitivity of selection: the change in the optimum phenotype with environmental change.
: maximum rate of environmental change: the highest rate of sustained environmental change under which long-term population persistence is possible.
Mosquito thermal adaptation: evidence and data gaps
Generation time (T) and maximum population growth rate ()
Short generation times enable rapid evolutionary responses (Lynch and Lande, 1993; Bürger and Lynch, 1995), and high intrinsic population growth rates reduce the chance of extinction prior to adaptation (Bürger and Lynch, 1995; Orr and Unckless, 2008; Gomulkiewicz and Houle, 2009). The generally rapid life cycles and large population sizes of mosquitoes favor rapid evolution, but precise demographic estimates under natural conditions are unavailable for most species (but see Appendix 2—table 1 for growth rates and generation times for Ae. aegypti, Anopheles spp., and Cx. pipiens) and will vary with biotic and abiotic conditions. However, even high estimates of mosquito lifespans of approximately 3 months (Macdonald, 1952; Nayar and Sauerman, 1971; Papadopoulos et al., 2016; Joubert et al., 2016) are on par with or well below those of other species that have already demonstrated evolutionary responses to climate change (e.g. Drosophila subobscura, Rodríguez-Trelles and Rodríguez, 1998; Balanyá et al., 2006; Tamiaschirus hudsonicus, Réale et al., 2003; Brassica rapa, Franks et al., 2007, several bird species, Gienapp et al., 2007; Cepaea nemoralis, Ożgo and Schilthuizen, 2012; Oncorhynchus gorbuscha, Kovach et al., 2012). Further, high intrinsic population growth rates (r) of 0.19–0.38 per generation have been calculated for several major vector species (Appendix 2—table 1, Equation 1, Figure 1; Amarasekare and Savage, 2012; Johnson et al., 2015; Mordecai et al., 2017; Shocket et al., 2020) and census population size estimates on the order of 1,000–10,000 individuals have been found across studies of varying mosquito species and settings (Touré et al., 1998; Lehmann et al., 1998; Maciel-de-Freitas et al., 2008; Neira et al., 2014; Le Goff et al., 2019). Placing these mosquito results in context, a Drosophila modeling study showed that growth rates and population sizes in this range facilitated population persistence for over 300 generations under heat-knockdown selection (Willi and Hoffmann, 2009). Mosquito demographic characteristics therefore favor thermal adaptation.
Variation in thermal tolerance ()
Higher genetically based variance in a trait results in higher rates of phenotypic evolution (Lande, 1976). While no studies (to our knowledge) have measured within-population variation in mosquito thermal tolerance, several studies have investigated variation between populations (Mogi, 1992; Reisen, 1995; Dodson et al., 2012; Vorhees et al., 2013; Ruybal et al., 2016; Chu et al., 2019). Overall, these studies find genetically based, but often trait-specific variation that did not always clearly support local thermal adaptation (i.e. a correlation between trait values and local climatic conditions; Appendix 3—table 1). Some studies have found thermal tolerance varying predictably with the population’s thermal environment of origin. For example, upper thermal limits of mosquito respiration and survival after heat shock were positively correlated with the temperature of origin for Cx. tarsalis and An. gambiae, respectively (Rocca et al., 2009; Vorhees et al., 2013). However, several other studies have found the opposite pattern (Ruybal et al., 2016), found minimal or no variation in thermal responses between populations (Dodson et al., 2012; Mogi, 1992), or found that certain populations had uniformly higher or lower trait performance at all experimental temperatures independent of their climate of origin (Ruybal et al., 2016; Reisen, 1995; Chu et al., 2019; Dodson et al., 2012). Taken together, mosquito populations do sometimes vary in their thermal performance, but there is no clear evidence for existing local thermal adaptation across temperature gradients of similar magnitude to those predicted by climate change over the next several decades. This may suggest either barriers to thermal adaptation or relatively weak selection on thermal performance (see ‘Strength of selection’ section). However, the lack of within-population sampling and/or idiosyncratic, trait-specific temperature relationships may have obscured true patterns of local adaptation (Bradshaw et al., 2000).
Heritability of thermal tolerance
Higher heritability—the proportion of phenotypic variance in a population attributable to genetic effects—enables faster evolutionary rates because populations respond more efficiently to selection (Falconer and Mackay, 1996). To our knowledge, there are no estimates of the heritability of trait thermal tolerance for mosquitoes. However, evolutionary theory and empirical work in other ectotherm taxa suggest that thermal tolerance heritability is generally low. In particular, highly polygenic, complex, or environmentally-dependent traits—as expected for thermal tolerances—typically have low heritability (Bay et al., 2017). Supporting this expectation, a meta-analysis of heritability data for upper thermal limits in Drosophila resulted in an overall estimate of 0.28 (i.e. 28% of the populations phenotypic variance is due to genetic variance; Diamond, 2017), which is similar to heritability estimates for other Drosophila life history traits (average h2 = 0.26; Roff and Mousseau, 1987, Mousseau and Roff, 1987) and indicates moderately low heritability (but see Jenkins and Hoffmann, 1994, h2 = 0.5). However, more recent evidence suggests Drosophila can rapidly adapt to novel temperatures through multiple, alternative genetic pathways that lead to similar increases in thermal tolerance (i.e. ‘genetic redundancy’), challenging the notion that highly polygenic traits have low heritability (Barghi et al., 2019). In general, uncertainty surrounding the ecological relevance of laboratory measurements of insect thermal tolerance (Terblanche et al., 2007; Chown et al., 2009; Mitchell and Hoffmann, 2010) and the divergent evolutionary histories of Drosopholids and mosquitoes limit our understanding of thermal tolerance trait heritability in mosquitoes.
Strength of selection ()
For heritable traits, stronger natural selection—the differential survival or reproduction of mosquitoes with different trait values—would lead to faster adaptive responses, despite causing high initial mortality (Box 1; Lynch and Lande, 1993; Hartl and Clark, 1997). Temperature-imposed selection on mosquitoes, which can be approximated from temperature-dependent survival rates (Box 1,; Falconer and Mackay, 1996), is likely to be strong. Upper thermal limits for adult and larval survival are as low as 32–38°C (reviewed in Mordecai et al., 2019), which many mosquito populations—particularly those in the tropics—already experience and will increasingly face in a warming climate (Deutsch et al., 2008). Further, steep declines in survival between thermal optima and critical limits have been observed across mosquito species (Focks et al., 1993; Alto and Juliano, 2001; Kamimura et al., 2002; Delatte et al., 2009; Muturi et al., 2011; Mordecai et al., 2019). This high selection pressure may facilitate mosquito adaptation, provided that heritable variation in trait thermal tolerance exists.
Phenotypic plasticity (b)
Phenotypic plasticity—the ability of individual genotypes to produce varying phenotypes based on the environment (West-Eberhard, 2003)—provides an alternative mechanism for coping with climate change that is more rapid than evolutionary adaptation. However, because plasticity impedes natural selection on genetically based variation, it may ultimately inhibit population persistence under long-term directional change (Gienapp et al., 2008; Whitman and Agrawal, 2009; Chevin et al., 2010; Chevin et al., 2013; Merilä and Hendry, 2014). For mosquitoes, potentially important plastic responses include changes in activity patterns, biting behavior, or microhabitat selection, thermal acclimation, and initiation of dormancy, as reviewed below (and see Appendix 3—table 2). Phenotypic plasticity may itself vary across genotypes and thus could evolve in response to environmental change, but experimental evidence of the evolution of plasticity is lacking (Dewitt et al., 1998; Scheiner and Berrigan, 1998; Stinchcombe et al., 2004). Overall, mosquitoes possess a variety of potential plastic responses, but the capacity for these responses to increase thermal tolerance, their potential fitness costs, and how these plastic responses might interact with the process of evolutionary adaptation remain poorly understood. Below, we review current knowledge of different potential plastic responses.
Behavioral thermoregulation
Larval and adult mosquitoes could temporarily cope with warming, particularly high-temperature extremes, through behavioral avoidance. In laboratory thermal preference studies, Aedes, Anopheles, and Culex spp. have demonstrated behavioral avoidance of high temperatures when exposed to a thermal gradient (Thomson, 1938; Blanford et al., 2009; Verhulst et al., 2020). In natural settings, several studies have shown shifts in the biting time or habitat selection of mosquitoes, particularly An. gambiae, seasonally or in response to insecticide spraying (Taylor, 1975; Reisen and Aslamkhan, 1978; Voorham, 2002; Pates and Curtis, 2005; Manda et al., 2011). While such behavioral shifts have not been conclusively linked to temperature, studies have found increasing usage of underground or shaded oviposition sites that was correlated with increasing temperature, and not associated with change in habitat availability or accompanied by genetic differentiation (Ae. aegypti, Somers et al., 2011; Chadee and Martinez, 2016). Similarly, larvae in permafrost regions were observed to rest in deeper, cooler pond water when surface water temperatures became exceptionally high (Ae. communis, a non-vector species; Haufe and Burgess, 1956). Seeking out and accessing cooler microclimates may buffer mosquitoes from warm temperature extremes, reducing mortality and decreasing the strength of selection. However, evidence for mosquito behavioral thermoregulation more generally remains limited (Paaijmans and Thomas, 2011; Waldock et al., 2013), and trade-offs in resource acquisition from restricted foraging and activity time, and a lack of readily available cooler microhabitats would constrain this behavior (Angilletta, 2009; Sears et al., 2016; Huey and Kingsolver, 2019). Conversely, the absence of evidence for this phenomenon may be due to measurement challenges associated with observing mosquitoes in the field (Paaijmans and Thomas, 2011).
Thermal acclimation
Increases in thermal tolerance after exposure to warmer temperatures during development—a form of thermal acclimation—have been documented in several mosquito species (An. albimanus, Benedict et al., 1991; An. arabiensis and An. funestus, Lyons et al., 2012; Cx. pipiens, Gray, 2013; Ae. aegypti, Sivan et al., 2021). However, increases in thermal limits were typically minimal, suggesting a limited capacity for thermal acclimation to reduce mortality at high temperatures and enable population persistence. For example, the critical thermal limits of respirometry, motor function, or survival increased by less than 2°C for populations developing in 5°C warmer environments (Benedict et al., 1991; Lyons et al., 2012; Gray, 2013). Similarly, critical thermal maxima varied minimally with acclimation temperatures across a diverse range of over 200 ectotherm species (Gunderson and Stillman, 2015; Somero et al., 2016; Heerwaarden et al., 2016; Rohr et al., 2018).
Dormancy
Temporarily unfavorable environmental conditions could be overcome through dormancy—the interruption or reduction of metabolic activity through diapause or quiescence—a response that has been demonstrated in all major vector species (reviewed in Diniz et al., 2017). Dry-season dormancy (i.e. aestivation) is likely one mechanism enabling An. gambiae and An. coluzzi to persist during the 3- to 6-month long dry season in the Sahel, as evidenced by very low population sizes during the dry season followed by rapid increases after the first rain (Lehmann et al., 2010; Lehmann et al., 2014; Adamou et al., 2011; Yaro et al., 2012; Dao et al., 2014). However, there are no known examples of dormancy mechanisms in ectotherms that respond solely to high temperatures, thus this may be an unlikely response for mosquitoes, particularly tropical species facing warming temperatures in humid environments.
Environmental sensitivity of selection (B)
Environmental sensitivity of selection refers to how the optimum phenotype shifts with changes in the environment and is typically measured as the slope of the relationship between the optimal trait value and the environmental variable (e.g. the rate of change in the optimal upper thermal limit of adult life span against maximum summer temperature; Figure 1; Chevin et al., 2010, Chevin et al., 2015). A larger difference between the environmental sensitivity of selection and phenotypic plasticity (i.e. a greater deviation in the phenotype from the optimal value) necessitates faster adaptation (Chevin et al., 2010). For mosquitoes, as for nearly all other organisms, the environmental sensitivity of trait thermal tolerance has not been empirically measured (Chevin et al., 2010). However, across mosquito populations (Lyons et al., 2012; Vorhees et al., 2013) and species (Mordecai et al., 2019), upper thermal limits for most mosquito life history traits were less variable than lower thermal limits and optima. These patterns could reflect strong environmental sensitivity on lower thermal tolerance, intermediate sensitivity for the optimum, and weak sensitivity on upper thermal tolerance. However, it may more likely reflect evolutionary constraints on upper thermal tolerance (Kellermann et al., 2012; Hoffmann et al., 2013), or result from competing selection pressures, greater metabolic costs of heat versus cold tolerance, genetic constraints, or gene flow hindering thermal adaptation (Angilletta, 2009; Kristensen et al., 2016).
Expected rate of environmental change ()
Rates of environmental change will vary based on the specific temperature variable being considered (e.g. mean temperature of the hottest month or quarter, maximum temperature in the dry season, etc.), and depend on climate policy: projections of global mean annual surface temperatures in 2100 vary by over 3°C depending on the future climate scenario (Collins et al., 2013). However, while greater warming is projected for higher latitudes (IPCC, 2007), faster rates of adaptation may be necessary for tropical mosquito populations that already experience environmental temperatures close to their thermal optima and may experience large fitness costs under additional warming in the absence of adaptation (Deutsch et al., 2008; Somero, 2010; Ryan et al., 2015; Mordecai et al., 2019). For example, although Ae. aegypti and Ae. albopictus are highly adaptable and have expanded into temperate climates, they are vulnerable to climate warming in their tropical ranges as temperatures here are expected to exceed their thermal optima and upper thermal limits in coming decades (Ryan et al., 2019).
Priorities and approaches for measuring adaptive potential
Addressing several key data and knowledge gaps will improve our ability to estimate mosquito adaptive potential (Table 1). As outlined above, there are virtually no estimates of the heritability of thermal tolerance traits, environmental sensitivity of selection, and within-population variation (and few estimates of between-population variation) in thermal tolerance for mosquitoes specifically. Additionally, we have a limited understanding of the role of phenotypic plasticity, particularly behavioral thermoregulation, in mosquito thermal tolerance. Although other parameters of the evolutionary rescue model (i.e. the strength of selection imposed by temperature change, mosquito generation time, and maximum population growth rate) are often not measured directly or precisely, we have relatively more information about these parameters and they are unlikely to be the primary constraints on evolutionary adaptation (see ‘Estimating evolutionary rates’). We therefore recommend that future research focus on measuring environmental sensitivity of selection, plasticity, and within-population variation and heritability in thermal tolerance. We discuss the most promising and feasible approaches for doing so below.
Selection experiments are a powerful tool for investigating the evolution of complex traits (reviewed in Fuller et al., 2005; Garland and Rose, 2009; Swallow et al., 2009) that can be used to estimate several of the parameters in evolutionary rescue model parameters. In artificial selection experiments, where individuals are chosen to advance to the next generation based on their value for a particular trait (e.g. time to thermal knockdown), heritability can be measured as: (Falconer and Mackay, 1996). Here, R is the mean difference in the trait between control and selected lines, is the trait standard deviation in the control lines, and i, the intensity of selection, is determined based on what proportion of the population is selected each generation (see Box 1). In laboratory natural selection—in which the treatment conditions, rather than the researcher, impose the selection pressure—selection strength itself can be approximated based on the survival rates between generations held at specific temperatures (see Box 1). Both selection designs have been used extensively with model organisms such as Drosophila spp., Daphnia spp., and Escherichia coli to measure changes in upper limits of trait thermal tolerance. While no thermal selection experiments have yet been published on mosquitoes (Dennington et al., in prep), several major vector species, including Ae. aegypti and Cx. quinquefasciatus, can be readily maintained and manipulated in the lab (Munstermann, 1997; Kauffman et al., 2017) and can therefore be used in experiments to obtain estimates of the heritability of thermal tolerance and the selection strength imposed by different temperature conditions.
Common garden experiments, where traits are measured for distinct populations or genotypes exposed to the same environmental conditions, enable measurement of nearly all rescue model parameters (Clausen et al., 1941; Merilä and Hendry, 2014; de Villemereuil et al., 2016; Berend et al., 2019). In mosquitoes, common garden experiments have been used to investigate variation in thermal tolerance between populations sampled across a thermal gradient (Mogi, 1992; Reisen, 1995; Rocca et al., 2009; Vorhees et al., 2013; Ruybal et al., 2016; Chu et al., 2019), but this approach could also be used to measure within-population variation if thermal tolerance traits were measured at the individual level. For example, larval development rates or adult survival time could be measured by tracking mosquitoes housed individually under different temperature conditions (as in Bedhomme et al., 2003). To measure genetically based variation in thermal tolerance, and to avoid confounding parental effects and thermal acclimation in the original environment, collected populations should be reared for at least one generation in the lab before experimentation. However, plasticity itself can be measured by, for example, varying larval rearing temperature (Dodson et al., 2012) or measuring thermoregulatory behavior as the trait of interest (e.g. Logan et al., 2018). Common garden experiments can also be used to measure the environmental sensitivity of selection if fitness is measured in addition to thermal tolerance traits (Chevin et al., 2010). Lastly, by tracking parentage and measuring thermal tolerance traits (e.g. by housing mating pairs of mosquitoes separately), the heritability of thermal tolerance can be measured based on the slope of the trait values of parent and offspring (Falconer and Mackay, 1996).
In addition to providing estimates of rescue model parameters, common garden experiments can be combined with genomic approaches to identify genetic variants associated with climate-adaptive traits (for examples in non-mosquito species, see De Kort et al., 2014; de Villemereuil et al., 2016; Exposito-Alonso et al., 2019; Capblancq et al., 2020). In the closest example of this approach in mosquito populations, the hypothesized thermo-adaptive role of a particular genotype (chromosomal inversion 2La) associated with aridity clines in Africa in An. gambiae (Coluzzi et al., 1979) was confirmed based on thermal tolerance experiments on the two genotypes (homokaryotypic populations 2La+ and 2La) (Rocca et al., 2009). In other taxa, common garden experiments have been combined with genome scans to quantify and predict climate-driven selection along the genome of the plant Arabidopsis thaliana (Exposito-Alonso et al., 2019), and to identify 162 candidate genes underlying climate adaptation in the harlequin fly Chironomus riparius (Waldvogel et al., 2018). In these studies, whole-genome sequencing would provide greater power to detect causal loci, thus this approach would be most feasible for mosquito species with available reference genomes, namely Ae. aegypti (e.g. Nene et al., 2007; Matthews et al., 2018), Ae. albopictus (Chen et al., 2015), An. darlingi (Marinotti et al., 2013), An. gambiae (Holt et al., 2002), and An. stephensi (Jiang et al., 2014).
Selection experiments and common garden experiments provide the means to obtain critical missing information on mosquito adaptive potential, but there are several challenges to these approaches. For any experimental test of adaptive potential, regardless of the methodology used, one must identify appropriate temperature treatments and assess thermal tolerance on the mosquito life history traits most relevant for fitness. Arbitrary choices for these details make it more difficult to extrapolate from these results to natural systems. Experiments commonly use treatments with constant temperatures above mean ambient temperatures. However, temperature minima or maxima, seasonal variability, and/or accumulated thermal stress may be more relevant to adaptive potential. For example, increases in minimum temperatures affect overnight recovery from heat stress in mosquitoes (Murdock et al., 2012; Bai et al., 2019). Further, given trade-offs in isolating the effect of temperature versus incorporating realistic ecological variation and in maximizing replication between versus within populations, no single study can definitively determine a species’ adaptive potential. As a first step, controlled and replicated lab studies measuring mosquito fitness (either directly or as a composite of individual life history traits) under realistic projected thermal regimes that incorporate natural diurnal variation in temperature, combined with genomics approaches, will greatly improve our understanding on current and potential mosquito thermal adaptation (Andriamifidy et al., 2019). Such studies will inform parameters of evolutionary rescue models and, more broadly, enable investigation of the dynamics and limits of thermal adaptation.
While these empirical approaches will address data gaps that we have emphasized within the evolutionary rescue framework (Table 1), the model itself (Chevin et al., 2010) has several important limitations. Notably, these include the lack of potential genotype-by-environment interactions in the expression of phenotypes, evolution in plasticity, gene flow, genetic correlations between traits associated with thermal tolerance, and demographic or environmental stochasticity (Chevin et al., 2010). These simplifying assumptions make the model tractable but may limit the accuracy of the predictions if these factors play a large role in adaptation. Adding complexity would require additional data collection and may make predictive models too computationally intensive to solve analytically but can be implemented through simulations (Bürger and Lynch, 1995). Several studies have effectively used simulations to incorporate environmental stochasticity (Ashander et al., 2016), demographic stochasticity (Martin et al., 2013), dispersal (Schiffers et al., 2013), carrying capacity (Bridle et al., 2010), and evolution in plasticity (Scheiner et al., 2017) into an evolutionary rescue model framework. Simulation results can be used to investigate transient evolutionary dynamics and can be compared with analytic results to determine the impact of these processes on evolutionary rescue. For example, Ashander et al., 2016 estimated population extinction risk using both analytic approximations and simulations to find that evolving plasticity only facilitated evolutionary rescue when the environmental change was sufficiently predictable. Using simulation to model more realistically complex evolutionary scenarios will likely be necessary when more precise forecasting is a priority, and is becoming a more approachable method through the availability of individual-based evolutionary simulation tools such as SLiM (Haller and Messer, 2017). For example, Matz et al., 2020 used the SLiM framework to estimate the adaptive potential of a coral metapopulation under varying levels of mutation, migration, and selection efficiency, enabling them to identify the main predictors of adaptive potential and the scenarios enabling long-term coral persistence.
Consequences for disease transmission
Current modeling approaches for predicting mosquito-borne disease transmission under climate change do not incorporate evolutionary adaptation in mosquito thermal tolerance. In particular, several studies have used a temperature-dependent modeling approach (where is the average number of secondary cases that result from a single infected individual introduced into a fully susceptible population) to estimate transmission of mosquito-borne diseases including dengue, chikungunya, Zika, and malaria under projected temperature conditions (e.g. Ryan et al., 2015; Ryan et al., 2019; Mordecai et al., 2017; Tesla et al., 2018). These studies rely on relationships between temperature and mosquito life history traits previously measured in the lab and currently provide the best estimates of mosquito-borne disease transmission under climate warming. However, if mosquito thermal optima and upper thermal limits increase, these predictions would underestimate future disease risk.
To investigate the consequences of shifts in mosquito thermal limits on disease transmission predictions, we consider a case study using Aedes aegypti-transmitted dengue virus in Northern Brazil (Appendix 4). This region, which includes the North and Northeastern Brazilian macroregions, experiences approximately 250,000 dengue cases annually (National Notifiable Diseases Information System (SINAN), 2019), primarily transmitted by Aedes aegypti (Chouin-Carneiro and Barreto dos Santos, 2017). In the absence of mosquito thermal adaptation, Ryan et al., 2019 projected that year-round transmission suitability would decrease in this area by 2080 under an upper climate change scenario (representative concentration pathway (RCP) 8.5). We repeat the modeling approach used in this projection to examine the rate of evolutionary adaptation required by Aedes aegypti to maintain current levels of dengue transmission suitability (Appendix 4). We use the same temperature-dependent model to estimate the number of months per year where temperatures do not prevent dengue transmission (i.e. > 0, as defined previously in Ryan et al., 2019) under current (2021) climate conditions and in 2080 under RCP 8.5. We then estimate the amount of evolutionary change in Aedes aegypti thermal limits necessary to maintain current levels of transmission suitability. We assume that adult fecundity is the mosquito trait under thermal selection as it has the lowest critical thermal maximum (34.61°C) of the Ae. aegypti and dengue virus life history traits and thus sets the warm temperature limit for dengue transmission (Mordecai et al., 2017; Mordecai et al., 2019). As in Ryan et al., 2019, we use mean monthly temperature when estimating temperature-based suitability for transmission, although this is not necessarily the climate factor that most strongly limits Ae. aegypti persistence in this region.
We find that in the absence of thermal adaptation in Ae. aegypti fecundity, the average number of months per year with suitable temperatures for dengue transmission in Northern Brazil would decrease from 12.0 in 2021 to 10.3 in 2080 (Figure 2). To maintain current (2021) levels of dengue transmission suitability under 2080 temperatures, the critical thermal maximum of Ae. aegypti fecundity would need to increase by an average of 1.57°C in this time period, or roughly 0.03°C/year. This evolutionary rate is on par with sustainable evolutionary rates estimated for other taxa and traits in the face of climate warming (e.g. great tit breeding time: 0.03–0.10°C/year; Gienapp et al., 2013). However, determining whether this is a plausible rate of evolutionary change in fecundity for Ae. aegypti in this region will require collecting missing information on the evolutionary rescue model parameters (Table 1) through the empirical approaches described above. In this case study, estimating the thermal adaptive potential of Ae. aegypti fecundity would help determine whether or not the dengue transmission season may decrease by nearly 2 months in a region containing approximately 69 million people (IBGE, 2010).
As explored in the case study, maintaining disease transmission under climate warming may require evolved increases in mosquito upper thermal limits. However, such evolutionary shifts could maintain, increase, or decrease transmission depending on whether they are accompanied by shifts in lower thermal limits, on the strength of thermodynamic constraints, and on genetic correlations between traits. In the absence of other changes to thermal performance, upward shifts in thermal limits could maintain current levels of disease transmission under rising temperatures, particularly if lower temperatures are infrequently experienced. However, disease transmission may increase if peak performances for mosquito traits like fecundity and biting rate increase with their thermal optima. This is an expectation of the ‘hotter-is-better’ hypothesis, but how the shape of thermal performance curves evolves is a point of ongoing debate and empirical uncertainty (Angilletta et al., 2010; Latimer et al., 2011; Kontopoulos et al., 2020). Regardless, genetic correlations between mosquito traits under direct selection and other traits that may impact disease transmission (e.g., development time and immunocompetence, as observed in Ae. aegypti; Koella and Boëte, 2002) could still constrain mosquito-borne disease transmission under thermal adaptation (Lande and Arnold, 1983).
Mosquitoes, like other ectotherms, may cope with warming temperatures through a variety of other mechanisms besides shifts in thermal physiology, such as accelerated life cycles, phenological shifts, and behavioral thermoregulation, with varying consequences for disease transmission (Huey and Kingsolver, 1993; Bradshaw et al., 2000; Stearns et al., 2000; Angilletta et al., 2003; Waldvogel et al., 2020). Evolved increases in life cycle speed can mitigate increases in daily mortality rates, and were suggested to occur in Anopheles spp. in response to vector control interventions (Ferguson et al., 2012). Adult mosquito longevity is already the main limitation on transmission near upper thermal limits for many major mosquito-borne diseases (Mordecai et al., 2019). Further reductions in adult lifespan could cause large declines in transmission for pathogens with longer incubation periods. In particular, transmission of malaria parasites, which have a minimum incubation period of approximately nine days (Paaijmans et al., 2012; Blanford et al., 2013), may be more negatively impacted under shortened mosquito lifespans than viral pathogens such as dengue virus and chikungunya virus, which have generally faster incubation periods—as low as three to five days at temperatures above 30°C (Tjaden et al., 2013; Rudolph et al., 2014; Mordecai et al., 2019; Winokur et al., 2020). The implications of warming-driven life cycle adaptation therefore depend on the interaction between vector and pathogen traits, which vary across species and environments.
Behavioral thermoregulation and phenological shifts could increase, maintain, or decrease disease transmission, primarily depending on how these shifts impact mosquito – human contact rates and the effectiveness of vector control activities (Ferreira et al., 2017). For example, if rising temperatures promote shifts in biting activity towards the cooler, night-time hours when humans are more likely to be protected by bed nets, disease transmission may be reduced (Taylor, 1975; Pates and Curtis, 2005; Moiroux et al., 2012; Thomsen et al., 2017; Carrasco et al., 2019). However, in the absence of vector control, shifts towards night-time biting, as well as thermoregulatory shifts favoring indoor versus outdoor biting, could increase mosquito – human contact rates and transmission (Takken, 2002). Similarly, phenological shifts in mosquito activity could lead to changes in the length or timing of disease transmission, potentially maintaining, increasing, or decreasing disease transmission. For example, increasing monthly mean temperatures in portions of California have effectively doubled the potential transmission season of St. Louis encephalitis virus, such that elderly persons traveling to California for the winter are newly at risk (Patz and Reisen, 2001). Failing to account for phenological shifts in mosquito activity may render vector control programs less effective at reducing mosquito populations and disease transmission. In general, the impact of mosquito thermal adaptation on disease transmission will vary based on the mechanism of thermal adaptation, making identifying what adaptive strategies are most likely in different contexts a priority for future research.
Applications to other vector and pest taxa
The same properties favoring mosquito thermal adaptation—short generation times, high population growth rates, and strong climate sensitivity—apply to many insect taxa that threaten human, animal, and plant health. These include vectors of major human, wildlife, and plant disease (e.g., species of tsetse flies, biting midges, psyllids, and aphids), as well as pests of crops and forest resources (e.g., species of beetles, moths, fruit flies, and fire ants). Despite the substantial societal cost adaptation in pest and disease vector species could impose, their potential to adapt to climate change remains poorly understood. This remains challenging to predict given the many determinants of evolutionary rates, incomplete data on these determinants for most taxa, and the inability to perform a single, conclusive experiment.
Drawing from conservation biology techniques used to study climate adaptive potentials in threatened and endangered species, we have outlined a framework and empirical approaches for investigating mosquito thermal adaptation that can be applied to any vector or pest species and type of environmental change. For example, in the Eastern U.S., range retractions of the invasive European gypsy moth have been linked to the duration of exposure above the optimal temperature for larval and pupal development (28°C; Tobin et al., 2014). Further, a recent common garden experiment found moth populations at the southern range edge, which experience the strongest thermal selection, have higher thermal tolerance in egg hatching, highlighting the potential for adaptive evolution in this species (Faske et al., 2019). The climate adaptive potential of this species could thus be estimated by comparing the potential evolutionary rates of thermal limits in immature development with projected rates of warming in this region. More precisely estimating the thermal adaptive potential of this pest species would enable forest management personnel to tailor intervention and control strategies in the face of ongoing warming. Similarly, estimates of mosquito thermal adaptive potential would enable vector control personnel to better target surveillance and insecticide applications to the appropriate locations and seasons of mosquito activity. More broadly, understanding and estimating the potential for climate adaptation in taxa of concern to human health is critical for accurately predicting and preparing for their persistence or shifts in their distributions under climate change.
Conclusion: How will mosquitoes adapt to climate warming?
Our synthesis makes clear that some general aspects of mosquito demographics and strong temperature-imposed selection may facilitate rapid evolution and adaptation to climate warming. In particular, typical mosquito generation times and population growth rates are on par with those of species that have already demonstrated evolutionary responses to climate change. Further, the steep declines in survival at temperatures exceeding mosquito thermal optima, as well as evidence of some population-level variation in mosquito temperature-trait responses, indicate that the selective pressures and raw genetic material necessary for evolutionary adaptation exist. However, making more accurate predictions about mosquito persistence and adaptation under climate warming will require identifying: (1) which life history traits experience the strongest thermal selection for a particular mosquito population, (2) how the optimal trait thermal tolerance varies with environmental temperature, (3) the extent of heritability and within-population variation in trait thermal tolerance, and (4) the role of phenotypic plasticity (particularly behavioral thermoregulation) in evolutionary adaptation and persistence. Empirical approaches such as common garden or selection experiments to obtain multiple pieces of missing information at once and leveraging information from related taxa where applicable can be used to address these key knowledge gaps. This would enable better estimates of mosquito adaptive potential and its implications for the future of mosquito-borne disease in a warming climate.
Data and code accessibility statement
All data and code supporting the results are either referenced in text, or available in the Appendices or on Github (https://github.com/lcouper/MosquitoAdaptationCaseStudy; copy archived at swh:1:rev:19a0d661adb3c0079bd5631be757cc1f255a854a, Couper, 2021).
Appendix 1
Additional factors influencing population persistence
In addition to the eight parameters in the Chevin et al., 2010 evolutionary rescue model (main text, Box 1), other factors influence a population’s potential for thermal adaptation. The starting population size informs the probability of obtaining rescue variants, and thus the probability of adaptation (Orr and Unckless, 2008; Bell, 2013; Martin et al., 2013; Carlson et al., 2014). The degree of stochasticity in temperature change influences an organism's ability to respond via plastic changes (Manenti et al., 2014; Catullo et al., 2015). The breadth of thermal performance and the genetic architecture of thermal tolerance inform the strength of selection (Huey and Kingsolver, 1993; Kopp and Matuszewski, 2014). Rates of gene flow and dispersal may hinder local adaptation due to the influx of maladapted genes but may promote evolutionary rescue through increasing genetic variation and bolstering small population sizes (Garant et al., 2007; Baskett and Gomulkiewicz, 2011; Kirkpatrick and Peischl, 2013; Schiffers et al., 2013; Bourne et al., 2014; Carlson et al., 2014). The degree of density dependence determines the strength of genetic bottlenecks (Nordstrom et al., 2020). Biotic interactions such as interspecific competition or predation may alter selection pressures (Culler et al., 2015), constraining rates of adaptation (Barbour et al., 2020), or limit population persistence despite sufficient rates of adaptation (Angilletta, 2009; Carlson et al., 2014; Johnson et al., 2019; Huey and Kingsolver, 2019). Fitness costs of phenotypic plasticity, such as reduced foraging associated with behavioral thermoregulation or energetic costs associated with maintaining the physiological machinery for acclimation may constrain levels of plasticity (Dewitt et al., 1998; Angilletta, 2009; Chevin et al., 2010).
Appendix 2
Methods for calculating mosquito population growth rate
The formula for population growth rate (r) as a function of temperature (T) is derived from the Euler-Lotka equation (Amarasekare and Savage, 2012) as follows:
Population growth rate is a function of adult mortality (µ), mosquito development rate (MDR), fecundity (E), and juvenile mortality (µj). W is the upper branch of the Lambert function.
Appendix 3
Appendix 4
Methods for ‘Consequences for disease transmission’ case study: Aedes aegypti-transmitted dengue virus in Northern Brazil
Case study context
In this case study, we estimate the monthly temperature suitability for dengue virus transmission by Aedes aegypti in Northern Brazil. This includes the North and Northeast macroregions, with a combined population size of approximately 69 million people (Instituto Brasileiro de Geografia e Estatistica and população, 2016) and 250,000 probable dengue cases per year (National Notifiable Diseases Information System (SINAN), 2019). Based on model projections that do not incorporate mosquito thermal adaptation (Ryan et al., 2019), Ae. aegypti-transmitted dengue suitability is expected to decline in this region by 2080 under an upper climate change scenario (RCP 8.5) as warming temperatures exceed mosquito upper thermal limits. Thermal adaptation by Ae. aegypti in this region could enable dengue transmission suitability to be maintained.
Modeling approach
We use a temperature-dependent model of —the number of secondary infections expected from a single infected individual introduced into a fully susceptible population—using the following expression (Mordecai et al., 2013):
where (T) indicates a temperature-dependent trait, a is the mosquito biting rate, b is the proportion of infectious mosquito bites resulting in infected humans, c is the proportion of bites on infected humans resulting in infected mosquitoes, µ is the adult mosquito mortality rate (here calculated as 1/ lifespan), PDR is the parasite development rate, EFD is the number of eggs laid per female per day, is the mosquito egg-to-adulthood survival probability, MDR is the sub-adult mosquito development rate, N is the density of humans, and r is the human recovery rate.
We use the data compiled by Mordecai et al., 2017 on the temperature-dependence of Aedes aegypti and dengue virus traits involved in transmission to parameterize the model. We use this model to estimate the number of months per year in which temperatures are suitable for dengue transmission (i.e., the number of months in which ), a conservative threshold that defines the range of temperatures at which dengue transmission is not prohibited (Ryan et al., 2019).
We first estimate the number of suitable months for transmission under current (2021) and future (2080) temperatures assuming no mosquito adaptation. We then estimate the extent of mosquito thermal adaptation needed to maintain current levels of transmission suitability under future temperatures. We assume adult fecundity to be the life history trait under thermal selection as it has lowest critical thermal maximum (; 34.61°C) of the Ae. aegypti and dengue virus life history traits (MacLean et al., 2019) and is thus the trait that determines the maximum temperature for dengue transmission. Current mean monthly temperatures do not exceed the for fecundity at any location within Northern Brazil in 2021. By 2080 under RCP 8.5, mean monthly temperatures are often expected to exceed , indicating temperatures are projected to be unsuitable for dengue transmission (i.e., ). We ask: What degree of increase in the of Ae. aegypti fecundity would enable current monthly transmission suitability to be maintained in 2080? In this example, we consider fecundity as the only trait under thermal selection, although future temperatures may exceed the of other Ae. aegypti and dengue virus life history traits, and there may be genetic correlations between fecundity and mosquito survival and development.
Climate data
For consistency with Ryan et al., 2019, we use mean monthly temperature when estimating Other climate variables such as diurnal temperature variation or monthly precipitation may limit future Ae. aegypti persistence in this region, but their effects on mosquito life history traits are less well understood than those of mean temperature, and thus they are not considered here for simplicity. Regions in which extreme drought or other limiting climate conditions are projected could be masked from the analyses as in Ryan et al., 2015 for more precise estimates. For current and projected mean monthly temperature data, we used the Hadley general circulation model (GCM) (HadGEM2-ES), which has high performance in Brazil and is the most commonly used GCM (Ryan et al., 2019; Almagro et al., 2020). Data were accessed through the Earth System Grid Federation. We used temperature projections made under representative concentration pathway (RCP) 8.5 for consistency with Ryan et al., 2019. RCP 8.5 is considered a ‘business-as-usual’ fossil fuel emissions scenario and corresponds to an 8.5 increase in solar radiation by 2100 (Riahi et al., 2011).
All climate analyses and mapping were conducted in R version 4.0.2. The R code and climate data files used in the analysis are available on Github (https://github.com/lcouper/MosquitoAdaptationCaseStudy).
Data availability
All data and code supporting the results are either referenced in text, or available in the Appendices or on Github (https://github.com/lcouper/MosquitoAdaptationCaseStudy; copy archived at https://archive.softwareheritage.org/swh:1:rev:19a0d661adb3c0079bd5631be757cc1f255a854a).
References
-
A framework for elucidating the temperature dependence of fitnessThe American Naturalist 179:178–191.https://doi.org/10.1086/663677
-
Do we know how mosquito disease vectors will respond to climate change?Emerging Topics in Life Sciences 3:115–132.https://doi.org/10.1042/ETLS20180125
-
Tradeoffs and the evolution of thermal reaction normsTrends in Ecology & Evolution 18:234–240.https://doi.org/10.1016/S0169-5347(03)00087-9
-
BookThermal Adaptation: A Theoretical and Empirical SynthesisOxford: OUP Oxford Press.https://doi.org/10.1093/acprof:oso/9780198570875.001.1
-
Thermodynamic effects on organismal performance: is hotter better?Physiological and Biochemical Zoology 83:197–206.https://doi.org/10.1086/648567
-
Predicting evolutionary rescue via evolving plasticity in stochastic environmentsProceedings of the Royal Society B: Biological Sciences 283:20161690.https://doi.org/10.1098/rspb.2016.1690
-
BookDrosophila A Laboratory HandbookCold Spring Harbor, New york: Cold Spring Harbor Laboratory Press.
-
Loss of consumers constrains phenotypic evolution in the resulting food webEvolution Letters 4:266–277.https://doi.org/10.1002/evl3.170
-
Unusual productivity of Aedes aegypti in septic tanks and its implications for dengue controlMedical and Veterinary Entomology 22:62–69.https://doi.org/10.1111/j.1365-2915.2008.00720.x
-
Introgressive hybridization as a mechanism for species rescueTheoretical Ecology 4:223–239.https://doi.org/10.1007/s12080-011-0118-0
-
Predicting responses to contemporary environmental change using evolutionary response architecturesThe American Naturalist 189:463–473.https://doi.org/10.1086/691233
-
Sex-specific reaction norms to intraspecific larval competition in the mosquito aedes aegyptiJournal of Evolutionary Biology 16:721–730.https://doi.org/10.1046/j.1420-9101.2003.00576.x
-
Evolutionary rescue and the limits of adaptationPhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120080.https://doi.org/10.1098/rstb.2012.0080
-
Evolutionary rescueAnnual Review of Ecology, Evolution, and Systematics 48:605–627.https://doi.org/10.1146/annurev-ecolsys-110316-023011
-
Heat-shock mortality and induced thermotolerance in larvae of the mosquito anopheles albimanusJournal of the American Mosquito Control Association 7:547–550.
-
Correlations between measures of heat resistance and acclimation in two species of Drosophila and their hybridsBiological Journal of the Linnean Society 64:449–462.https://doi.org/10.1111/j.1095-8312.1998.tb00343.x
-
Relationship between altitude and intensity of malaria transmission in the Usambara Mountains, TanzaniaJournal of Medical Entomology 40:706–717.https://doi.org/10.1603/0022-2585-40.5.706
-
Between migration load and evolutionary rescue: dispersal, adaptation and the response of spatially structured populations to environmental changeProceedings of the Royal Society B: Biological Sciences 281:20132795.https://doi.org/10.1098/rspb.2013.2795
-
The mean and variance of environmental temperature interact to determine physiological tolerance and fitnessPhysiological and Biochemical Zoology 84:543–552.https://doi.org/10.1086/662551
-
How extreme temperatures impact organisms and the evolution of their thermal toleranceIntegrative and Comparative Biology 56:98–109.https://doi.org/10.1093/icb/icw004
-
Genomic prediction of (Mal)Adaptation across current and future climatic landscapesAnnual Review of Ecology, Evolution, and Systematics 51:245–269.https://doi.org/10.1146/annurev-ecolsys-020720-042553
-
Evolutionary rescue in a changing worldTrends in Ecology & Evolution 29:521–530.https://doi.org/10.1016/j.tree.2014.06.005
-
Behavioural adaptations of mosquito vectors to insecticide controlCurrent Opinion in Insect Science 34:48–54.https://doi.org/10.1016/j.cois.2019.03.005
-
Evolutionary potential of thermal preference and heat tolerance in Drosophila subobscuraJournal of Evolutionary Biology 32:818–824.https://doi.org/10.1111/jeb.13483
-
The role of vector trait variation in Vector-Borne disease dynamicsFrontiers in Ecology and Evolution 8:189.https://doi.org/10.3389/fevo.2020.00189
-
An analytical study of population growth in Drosophila melanogasterEcological Monographs 20:173–206.https://doi.org/10.2307/1948580
-
ReportTransmission of Major Arboviruses in Brazil: The Role of Aedes Aegypti and Aedes Albopictus Vectors. Page Biological Control of Pest and Vector InsectsInTech, Rijeka, Croatia.
-
BookAëdes Aegypti (L.), the Yellow Fever Mosquito: Its Life History, Bionomics, and StructureCambridge: Cambridge University Press.https://doi.org/10.1126/science.133.3463.1473-a
-
Regional differentiation in plant speciesThe American Naturalist 75:231–250.https://doi.org/10.1086/280955
-
ReportLong-Term Climate Change: Projections, Commitments and IrreversibilityClimate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.
-
Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complexTransactions of the Royal Society of Tropical Medicine and Hygiene 73:483–497.https://doi.org/10.1016/0035-9203(79)90036-1
-
In a warmer arctic, mosquitoes avoid increased mortality from predators by growing fasterProceedings of the Royal Society B: Biological Sciences 282:20151549.https://doi.org/10.1098/rspb.2015.1549
-
Costs and limits of phenotypic plasticityTrends in Ecology & Evolution 13:77–81.https://doi.org/10.1016/S0169-5347(97)01274-3
-
Evolutionary potential of upper thermal tolerance: biogeographic patterns and expectations under climate change: biogeography of evolutionary potentialAnnals of the New York Academy of Sciences 1389:5–19.https://doi.org/10.1111/nyas.13223
-
Review: thermal preference in DrosophilaJournal of Thermal Biology 34:109–119.https://doi.org/10.1016/j.jtherbio.2008.11.007
-
Thermal variation, thermal extremes and the physiological performance of individualsJournal of Experimental Biology 218:1956–1967.https://doi.org/10.1242/jeb.114926
-
Evidence that implicit assumptions of 'no evolution' of disease vectors in changing environments can be violated on a rapid timescalePhilosophical Transactions of the Royal Society B: Biological Sciences 370:20140136.https://doi.org/10.1098/rstb.2014.0136
-
Development, reproductive output and population growth of the fruit fly pest Drosophila suzukii (Diptera: drosophilidae) on artificial dietJournal of Economic Entomology 107:1392–1398.https://doi.org/10.1603/EC13504
-
Genomic basis and evolutionary potential for extreme drought adaptation in Arabidopsis thalianaNature Ecology & Evolution 2:352–358.https://doi.org/10.1038/s41559-017-0423-0
-
Drosophila melanogaster as a model system to study mitochondrial biologyMethods in Molecular Biology 372:33–49.https://doi.org/10.1007/978-1-59745-365-3_3
-
Dynamic life table model for aedes aegypti (Diptera: culicidae): analysis of the literature and model developmentJournal of Medical Entomology 30:1003–1017.https://doi.org/10.1093/jmedent/30.6.1003
-
The effect of global change on mosquito-borne diseaseThe Lancet Infectious Diseases 19:e302–e312.https://doi.org/10.1016/S1473-3099(19)30161-6
-
Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradientMolecular Ecology 28:2206–2223.https://doi.org/10.1111/mec.15069
-
How and when selection experiments might actually be usefulIntegrative and Comparative Biology 45:391–404.https://doi.org/10.1093/icb/45.3.391
-
BookExperimental Evolution: Concepts, Methods, and Applications of Selection ExperimentsBerkeley: University of California Press.https://doi.org/10.1093/icb/icq093
-
No place to hide: microclimates of sonoran desert DrosophilaJournal of Thermal Biology 28:353–362.https://doi.org/10.1016/S0306-4565(03)00011-1
-
Predicting demographically sustainable rates of adaptation: can great tit breeding time keep pace with climate change?Philosophical Transactions of the Royal Society B: Biological Sciences 368:20120289.https://doi.org/10.1098/rstb.2012.0289
-
Resting behaviour, ecology and genetics of malaria vectors in large scale agricultural Areas of western KenyaParassitologia 38:481–489.
-
Demographic and genetic constraints on evolutionThe American Naturalist 174:E218–E229.https://doi.org/10.1086/645086
-
Evolutionary rescue beyond the modelsPhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120093.https://doi.org/10.1098/rstb.2012.0093
-
Evolutionary rescue: an emerging focus at the intersection between ecology and evolutionPhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120404.https://doi.org/10.1098/rstb.2012.0404
-
Plasticity in thermal tolerance has limited potential to buffer ectotherms from global warmingProceedings of the Royal Society B: Biological Sciences 282:20150401.https://doi.org/10.1098/rspb.2015.0401
-
SLiM 2: flexible, interactive forward genetic simulationsMolecular Biology and Evolution 34:230–240.https://doi.org/10.1093/molbev/msw211
-
BookPrinciples of Population Genetics (Fourth Edition)Sunderland: Sinauer Associates.
-
Limited scope for plasticity to increase upper thermal limitsFunctional Ecology 30:1947–1956.https://doi.org/10.1111/1365-2435.12687
-
Geographical variation in the acclimation responses of Drosophila to temperature extremesThe American Naturalist 142 Suppl 1:S93–S113.https://doi.org/10.1086/285525
-
Temperature, demography, and ectotherm fitnessThe American Naturalist 158:204–210.https://doi.org/10.1086/321314
-
Evolution of resistance to high temperature in ectothermsThe American Naturalist 142:S21–S46.https://doi.org/10.1086/285521
-
Climate warming, resource availability, and the metabolic meltdown of ectothermsThe American Naturalist 194:E140–E150.https://doi.org/10.1086/705679
-
ReportClimate Change 2007: The Physical Science BasisWorking Group I Contribution to the Fourth Assessment Report of the IPCC.
-
An ecological framework for modeling the geography of disease transmissionTrends in Ecology & Evolution 34:655–668.https://doi.org/10.1016/j.tree.2019.03.004
-
Effect of temperature on the development of aedes aegypti and Aedes albopictusMedical Entomology and Zoology 53:53–58.https://doi.org/10.7601/mez.53.53_1
-
Climate change drives microevolution in a wild birdNature Communications 2:1–7.https://doi.org/10.1038/ncomms1213
-
Invasion biology of Aedes japonicus japonicus (Diptera: culicidae)Annual Review of Entomology 59:31–49.https://doi.org/10.1146/annurev-ento-011613-162012
-
The well-temperatured biologist. (American society of naturalists presidential address)The American Naturalist 174:755–768.https://doi.org/10.1086/648310
-
Heat stress and the fitness consequences of climate change for terrestrial ectothermsFunctional Ecology 27:1415–1423.https://doi.org/10.1111/1365-2435.12145
-
Evolutionary rescue by beneficial mutations in environments that change in space and timePhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120082.https://doi.org/10.1098/rstb.2012.0082
-
Rapid evolution of quantitative traits: theoretical perspectivesEvolutionary Applications 7:169–191.https://doi.org/10.1111/eva.12127
-
Genetic change for earlier migration timing in a pink salmon populationProceedings of the Royal Society B: Biological Sciences 279:3870–3878.https://doi.org/10.1098/rspb.2012.1158
-
Proteomic data reveals a physiological basis for costs and benefits associated with thermal acclimationJournal of Experimental Biology 16:969–976.https://doi.org/10.1242/jeb.132696
-
The measurement of selection on correlated charactersEvolution 37:1210–1226.https://doi.org/10.2307/2408842
-
Introduction and establishment of Aedes (Finlaya) Japonicus japonicus (Theobald) on the island of hawaii: implications for arbovirus transmissionJournal of the American Mosquito Control Association 21:318–321.https://doi.org/10.2987/8756-971X(2005)21[318:IAEOAF]2.0.CO;2
-
Quantitative genetic variation for thermal performance curves within and among natural populations of Drosophila serrataJournal of Evolutionary Biology 24:965–975.https://doi.org/10.1111/j.1420-9101.2011.02227.x
-
The effective population size of anopheles gambiae in Kenya: implications for population structureMolecular Biology and Evolution 15:264–276.https://doi.org/10.1093/oxfordjournals.molbev.a025923
-
Aestivation of the african malaria mosquito, anopheles gambiae in the sahelThe American Journal of Tropical Medicine and Hygiene 83:601–606.https://doi.org/10.4269/ajtmh.2010.09-0779
-
Seasonal variation in spatial distributions of anopheles gambiae in a sahelian village: evidence for aestivationJournal of Medical Entomology 51:27–38.https://doi.org/10.1603/ME13094
-
Complex responses of global insect pests to climate warmingFrontiers in Ecology and the Environment 18:141–150.https://doi.org/10.1002/fee.2160
-
Disparate patterns of thermal adaptation between life stages in temperate vs. tropical Drosophila melanogasterJournal of Evolutionary Biology 31:323–331.https://doi.org/10.1111/jeb.13234
-
Consequences of heat hardening on a field fitness component in Drosophila depend on environmental temperatureThe American Naturalist 169:175–183.https://doi.org/10.1086/510632
-
Thermal physiology and thermoregulatory behaviour exhibit low heritability despite genetic divergence between lizard populationsProceedings of the Royal Society B: Biological Sciences 285:20180697.https://doi.org/10.1098/rspb.2018.0697
-
Evolution and extinction in response to environmental changeBiotic Interactions and Global Change 49:234–250.https://doi.org/10.1111/j.1558-5646.2007.00301.x
-
Survive a warming climate: insect responses to extreme high temperaturesAnnual Review of Entomology 66:163–184.https://doi.org/10.1146/annurev-ento-041520-074454
-
The analysis of the sporozoite rateTropical Diseases Bulletin 49:569–586.https://doi.org/10.1016/j.actatropica.2016.02.016
-
Evolution and plasticity of thermal performance: an analysis of variation in thermal tolerance and fitness in 22 Drosophila speciesPhilosophical Transactions of the Royal Society B: Biological Sciences 374:20180548.https://doi.org/10.1098/rstb.2018.0548
-
Effects of irritant chemicals on aedes aegypti resting behavior: is there a simple shift to untreated "safe sites"?PLOS Neglected Tropical Diseases 5:e1243.https://doi.org/10.1371/journal.pntd.0001243
-
Predictability rather than amplitude of temperature fluctuations determines stress resistance in a natural population of Drosophila simulansJournal of Evolutionary Biology 27:2113–2122.https://doi.org/10.1111/jeb.12463
-
The genome of anopheles darlingi, the main neotropical malaria vectorNucleic Acids Research 41:7387–7400.https://doi.org/10.1093/nar/gkt484
-
The probability of evolutionary rescue: towards a quantitative comparison between theory and evolution experimentsPhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120088.https://doi.org/10.1098/rstb.2012.0088
-
Estimating the potential for coral adaptation to global warming across the Indo-West pacificGlobal Change Biology 26:3473–3481.https://doi.org/10.1111/gcb.15060
-
Rapid local adaptation to northern winters in the invasive Asian tiger mosquito Aedes albopictus : A moving targetJournal of Applied Ecology 56:2518–2527.https://doi.org/10.1111/1365-2664.13480
-
Climate change, adaptation, and phenotypic plasticity: the problem and the evidenceEvolutionary Applications 7:1–14.https://doi.org/10.1111/eva.12137
-
Malaria vector productivity in relation to the Highland environment in KenyaThe American Journal of Tropical Medicine and Hygiene 75:448–453.https://doi.org/10.4269/ajtmh.2006.75.448
-
Temperature and photoperiod effects on larval and ovarian development of New Zealand strains of Culex quinquefasciatus (Diptera: Culicidae)Annals of the Entomological Society of America 85:58–66.https://doi.org/10.1093/aesa/85.1.58
-
Changes in anopheles funestus biting behavior following universal coverage of long-lasting insecticidal nets in BeninThe Journal of Infectious Diseases 206:1622–1629.https://doi.org/10.1093/infdis/jis565
-
BookThe Gypsy MothIn: Berryman A. A, editors. Dynamics of Forest Insect Populations: Patterns, Causes, Implications. Boston: Springer. pp. 353–375.https://doi.org/10.1007/978-1-4899-0789-9
-
Evolutionary responses to global change: lessons from invasive speciesEcology Letters 17:637–649.https://doi.org/10.1111/ele.12262
-
Detecting the impact of temperature on transmission of zika, dengue, and Chikungunya using mechanistic modelsPLOS Neglected Tropical Diseases 11:e0005568.https://doi.org/10.1371/journal.pntd.0005568
-
Thermal biology of mosquito-borne diseaseEcology Letters 22:1690–1708.https://doi.org/10.1111/ele.13335
-
Climate change could shift disease burden from malaria to arboviruses in africaThe Lancet Planetary Health 4:e416–e423.https://doi.org/10.1016/S2542-5196(20)30178-9
-
Natural selection and the heritability of fitness componentsHeredity 59 (Pt 2):181–197.https://doi.org/10.1038/hdy.1987.113
-
BookCare and maintenance of Aedes mosquito coloniesIn: Crampton J. M, Beard C. B, Louis C, editors. The Molecular Biology of Insect Disease Vectors: A Methods Manual. Dordrecht: Springer Netherlands. pp. 13–20.https://doi.org/10.1007/978-94-009-1535-0
-
Complex effects of temperature on mosquito immune functionProceedings of the Royal Society B: Biological Sciences 279:3357–3366.https://doi.org/10.1098/rspb.2012.0638
-
Effect of temperature and insecticide stress on life-history traits of culex restuans and aedes albopictus (Diptera: culicidae)Journal of Medical Entomology 48:243–250.https://doi.org/10.1603/ME10017
-
The effects of diet on life-span, fecundity and flight potential of aedes taeniorhynchus adultsJournal of Medical Entomology 8:506–513.https://doi.org/10.1093/jmedent/8.5.506
-
Estimation of aedes aegypti (Diptera: Culicidae) population size and adult male survival in an urban area in PanamaMemórias Do Instituto Oswaldo Cruz 109:879–886.https://doi.org/10.1590/0074-0276140136
-
ConferenceNegative density dependence constrains evolutionary rescueESA Annual Meeting.
-
Population extinction and the genetics of adaptationThe American Naturalist 172:160–169.https://doi.org/10.1086/589460
-
Thermal tolerance in widespread and tropical Drosophila species: does phenotypic plasticity increase with latitude?The American Naturalist 178 Suppl 1:S80–S96.https://doi.org/10.1086/661780
-
Evolutionary change in Cepaea nemoralis shell colour over 43 yearsGlobal Change Biology 18:74–81.https://doi.org/10.1111/j.1365-2486.2011.02514.x
-
The effect of water turbidity on the near-surface water temperature of larval habitats of the malaria mosquito anopheles gambiaeInternational Journal of Biometeorology 52:747–753.https://doi.org/10.1007/s00484-008-0167-2
-
Temperature variation makes ectotherms more sensitive to climate changeGlobal Change Biology 19:2373–2380.https://doi.org/10.1111/gcb.12240
-
Seasonality of Post-capture longevity in a Medically-Important mosquito (Culex pipiens)Frontiers in Ecology and Evolution 4:63.https://doi.org/10.3389/fevo.2016.00063
-
Mosquito behavior and vector controlAnnual Review of Entomology 50:53–70.https://doi.org/10.1146/annurev.ento.50.071803.130439
-
Disease emergence from global climate and land use changeMedical Clinics of North America 92:1473–1491.https://doi.org/10.1016/j.mcna.2008.07.007
-
Immunology, climate change and vector-borne diseasesTrends in Immunology 22:171–172.https://doi.org/10.1016/S1471-4906(01)01867-1
-
The potential effects of climatic change on agricultural insect pestsAgricultural and Forest Meteorology 57:221–240.https://doi.org/10.1016/0168-1923(91)90088-8
-
Genetic and plastic responses of a northern mammal to climate changeProceedings of the Royal Society of London. Series B: Biological Sciences 270:591–596.https://doi.org/10.1098/rspb.2002.2224
-
Ochlerotatus japonicus japonicus (Theobald) (Diptera: culicidae), a new invasive mosquito for Georgia and south CarolinaProceedings of the Entomological Society of Washington 106:233–234.
-
Effect of temperature on culex tarsalis (Diptera: culicidae) from the Coachella and san joaquin valleys of CaliforniaJournal of Medical Entomology 32:636–645.https://doi.org/10.1093/jmedent/32.5.636
-
Biting rhythms of some Pakistan mosquitoes (Diptera: culicidae)Bulletin of Entomological Research 68:313–330.https://doi.org/10.1017/S0007485300007392
-
Climate change and mosquito-borne diseaseEnvironmental Health Perspectives 109 Suppl 1:141–161.https://doi.org/10.1289/ehp.01109s1141
-
Quantitative genetics and fitness: lessons from DrosophilaHeredity 58 (Pt 1):103–118.https://doi.org/10.1038/hdy.1987.15
-
Frontiers in climate change-disease researchTrends in Ecology & Evolution 26:270–277.https://doi.org/10.1016/j.tree.2011.03.002
-
The complex drivers of thermal acclimation and breadth in ectothermsEcology Letters 21:1425–1439.https://doi.org/10.1111/ele.13107
-
Genetic variation for tolerance to high temperatures in a population of Drosophila melanogasterEcology and Evolution 8:10374–10383.https://doi.org/10.1002/ece3.4409
-
Incubation periods of mosquito-borne viral infections: a systematic reviewThe American Journal of Tropical Medicine and Hygiene 90:882–891.https://doi.org/10.4269/ajtmh.13-0403
-
Mapping physiological suitability limits for malaria in africa under climate changeVector-Borne and Zoonotic Diseases 15:718–725.https://doi.org/10.1089/vbz.2015.1822
-
Global expansion and redistribution of Aedes-borne virus transmission risk with climate changePLOS Neglected Tropical Diseases 13:e0007213.https://doi.org/10.1371/journal.pntd.0007213
-
Limited evolutionary rescue of locally adapted populations facing climate changePhilosophical Transactions of the Royal Society B: Biological Sciences 368:20120083.https://doi.org/10.1098/rstb.2012.0083
-
BookEnvironmental drivers of vector-borne diseaseIn: Shocket M. S, editors. Population Biology of Vector-Borne Diseases. Oxford University Press. pp. 85–118.https://doi.org/10.1093/oso/9780198853244.001.0001
-
Zika and Chikungunya: mosquito-borne viruses in a changing worldAnnals of the New York Academy of Sciences 1399:61–77.https://doi.org/10.1111/nyas.13306
-
Population growth of Drosophila melanogaster (Diptera: drosophilidae) at constant and alternating Temperatures1Annals of the Entomological Society of America 65:993–1001.https://doi.org/10.1093/aesa/65.5.993
-
Thermotolerance and acclimation in the immature stages of aedes aegypti (L) (Diptera: Culicidae) to simulated thermal stressInternational Journal of Tropical Insect Science 41:333–344.https://doi.org/10.1007/s42690-020-00211-x
-
Proteins and temperatureAnnual Review of Physiology 57:43–68.https://doi.org/10.1146/annurev.ph.57.030195.000355
-
Protein adaptations to temperature and pressure: complementary roles of adaptive changes in amino acid sequence and internal milieuComparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology 136:577–591.https://doi.org/10.1016/S1096-4959(03)00215-X
-
The physiology of climate change: how potentials for acclimatization and genetic adaptation will determine ‘winners’ and ‘losers’Journal of Experimental Biology 213:912–920.https://doi.org/10.1242/jeb.037473
-
Biochemical AdaptationResponse to Environmental Challenges from Life's Origins to the Anthropocene, Biochemical Adaptation, Sunderland, Sinauer Associates, Inc.
-
Genetics and morphology of aedes aegypti (Diptera: culicidae) in septic tanks in puerto ricoJournal of Medical Entomology 48:1095–1102.https://doi.org/10.1603/ME11129
-
Flowering time plasticity in Arabidopsis thaliana: a reanalysis of Westerman & Lawrence (1970)Journal of Evolutionary Biology 17:197–207.https://doi.org/10.1046/j.1420-9101.2003.00641.x
-
Global analyses of satellite-derived vegetation index related to climatological wetness and warmthInternational Journal of Climatology 26:425–438.https://doi.org/10.1002/joc.1256
-
BookSelection Experiments and Experimental Evolution of Performance and PhysiologyIn: Garland T, Rose M. R, editors. Experimental Evolution: Concepts, Methods, and Applications of Selection Experiments. Berkeley: University of California Press. pp. 909–910.https://doi.org/10.1093/icb/icq093
-
Do insecticide-treated bednets have an effect on malaria vectors?Tropical Medicine and International Health 7:1022–1030.https://doi.org/10.1046/j.1365-3156.2002.00983.x
-
Negligible senescence during reproductive dormancy in Drosophila melanogasterThe American Naturalist 158:248–258.https://doi.org/10.1086/321320
-
Changes in the feeding behaviour of a malaria vector, anopheles farauti Lav., following use of DDT as a residual spray in houses in the British Solomon Islands ProtectorateTransactions of the Royal Entomological Society of London 127:277–292.https://doi.org/10.1111/j.1365-2311.1975.tb00576.x
-
Critical thermal limits depend on methodological contextProceedings of the Royal Society B: Biological Sciences 274:2935–2943.https://doi.org/10.1098/rspb.2007.0985
-
Temperature drives Zika virus transmission: evidence from empirical and mathematical modelsProceedings of the Royal Society B: Biological Sciences 285:20180795.https://doi.org/10.1098/rspb.2018.0795
-
Mosquito behavior change after distribution of bednets results in decreased protection against malaria exposureThe Journal of Infectious Diseases 215:790–797.https://doi.org/10.1093/infdis/jiw615
-
The reactions of mosquitoes to temperature and humidityBulletin of Entomological Research 29:125–140.https://doi.org/10.1017/S0007485300026158
-
Extrinsic incubation period of dengue: knowledge, backlog, and applications of temperature dependencePLOS Neglected Tropical Diseases 7:e2207.https://doi.org/10.1371/journal.pntd.0002207
-
Supraoptimal temperatures influence the range dynamics of a non-native insectDiversity and Distributions 20:813–823.https://doi.org/10.1111/ddi.12197
-
Increased temperature variation poses a greater risk to species than climate warmingProceedings of the Royal Society B: Biological Sciences 281:20132612.https://doi.org/10.1098/rspb.2013.2612
-
Thermal preferences of subtropical Aedes aegypti and temperate Ae. japonicus mosquitoesJournal of Thermal Biology 91:102637.https://doi.org/10.1016/j.jtherbio.2020.102637
-
Adaptive phenotypic plasticity: consensus and controversyTrends in Ecology & Evolution 10:212–217.https://doi.org/10.1016/S0169-5347(00)89061-8
-
Thermal resistance and performance correlate with climate in populations of a widespread mosquitoPhysiological and Biochemical Zoology 86:73–81.https://doi.org/10.1086/668851
-
The role of environmental variables on Aedes albopictus biology and chikungunya epidemiologyPathogens and Global Health 107:224–241.https://doi.org/10.1179/2047773213Y.0000000100
-
The genomic footprint of climate adaptation in chironomus ripariusMolecular Ecology 27:1439–1456.https://doi.org/10.1111/mec.14543
-
Maximum voluntary temperature of insect larvae reveals differences in their thermal biologyJournal of Thermal Biology 33:380–384.https://doi.org/10.1016/j.jtherbio.2008.06.002
-
BookDevelopmental Plasticity and EvolutionOxford University Press.https://doi.org/10.1038/hdy.2015.14
-
BookWhat is Phenotypic Plasticity and Why is it Important?In: Whitman D, Ananthakrishnan T, editors. Phenotypic Plasticity of Insects. Science Publishers. pp. 1–63.
-
Demographic factors and genetic variation influence population persistence under environmental changeJournal of Evolutionary Biology 22:124–133.https://doi.org/10.1111/j.1420-9101.2008.01631.x
-
Impact of temperature on the extrinsic incubation period of zika virus in aedes aegyptiPLOS Neglected Tropical Diseases 14:e0008047.https://doi.org/10.1371/journal.pntd.0008047
-
ReportA Global Brief on Vector-Borne DiseasesWorld Health Organization Technical Report.
-
Climate change unlikely to increase malaria burden in West AfricaNature Climate Change 6:1009–1013.https://doi.org/10.1038/nclimate3085
-
Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengueEpidemiology and Infection 137:1188–1202.https://doi.org/10.1017/S0950268809002040
-
Dry season reproductive depression of Anopheles gambiae in the SahelJournal of Insect Physiology 58:1050–1059.https://doi.org/10.1016/j.jinsphys.2012.04.002
Article and author information
Author details
Funding
National Institute of General Medical Sciences (R35GM133439)
- Johannah E Farner
- Devin G Kirk
- Eloise B Skinner
- Erin A Mordecai
National Science Foundation (DEB-1518681)
- Johannah E Farner
- Marta Shocket
- Erin A Mordecai
National Science Foundation (DEB-2011147)
- Erin A Mordecai
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank George Somero for enlightening discussions about the molecular mechanisms of thermal tolerance. We thank Nina Dennington for sharing insight on experimental methods for studying mosquito thermal adaptation. We are grateful to Chris Anderson, Tejas Athni, Alex Becker, Caroline Glidden, and Morgan Kain for helpful feedback on the manuscript. EAM, DGK, JEF, and EBS were funded by the NIH National Institute of General Medical Sciences R35 MIRA program (R35GM133439). EAM, JEF, and MSS were funded by the NSF Ecology and Evolution of Infectious Diseases (EEID) program (DEB-1518681). EAM was funded by the NSF EEID program (DEB-2011147, with support from the Fogarty International Center), a Terman Award, a Stanford King Center for Global Development seed grant, and a seed grant from the Stanford Center for Innovation in Global Health and Woods Institute for the Environment. MLC was supported by the Illich-Sadowsky Fellowship through the Stanford Interdisciplinary Graduate Fellowship program. NN was supported by the Stanford Data Science Scholarship and the P.E.O. Scholar Award, International Chapter of the P.E.O. Sisterhood. LIC was funded by the Stanford Graduate Fellowship. JEF was funded by the Bing-Mooney Graduate Fellowship. MJH was funded by the Knight-Hennessy Scholarship.
Copyright
© 2021, Couper et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
Climatic warming can shift community composition driven by the colonization-extinction dynamics of species with different thermal preferences; but simultaneously, habitat fragmentation can mediate species’ responses to warming. As this potential interactive effect has proven difficult to test empirically, we collected data on birds over 10 years of climate warming in a reservoir subtropical island system that was formed 65 years ago. We investigated how the mechanisms underlying climate-driven directional change in community composition were mediated by habitat fragmentation. We found thermophilization driven by increasing warm-adapted species and decreasing cold-adapted species in terms of trends in colonization rate, extinction rate, occupancy rate and population size. Critically, colonization rates of warm-adapted species increased faster temporally on smaller or less isolated islands; cold-adapted species generally were lost more quickly temporally on closer islands. This provides support for dispersal limitation and microclimate buffering as primary proxies by which habitat fragmentation mediates species range shift. Overall, this study advances our understanding of biodiversity responses to interacting global change drivers.
-
- Ecology
- Microbiology and Infectious Disease
Predicting how species diversity changes along environmental gradients is an enduring problem in ecology. In microbes, current theories tend to invoke energy availability and enzyme kinetics as the main drivers of temperature-richness relationships. Here, we derive a general empirically-grounded theory that can explain this phenomenon by linking microbial species richness in competitive communities to variation in the temperature-dependence of their interaction and growth rates. Specifically, the shape of the microbial community temperature-richness relationship depends on how rapidly the strength of effective competition between species pairs changes with temperature relative to the variance of their growth rates. Furthermore, it predicts that a thermal specialist-generalist tradeoff in growth rates alters coexistence by shifting this balance, causing richness to peak at relatively higher temperatures. Finally, we show that the observed patterns of variation in thermal performance curves of metabolic traits across extant bacterial taxa is indeed sufficient to generate the variety of community-level temperature-richness responses observed in the real world. Our results provide a new and general mechanism that can help explain temperature-diversity gradients in microbial communities, and provide a quantitative framework for interlinking variation in the thermal physiology of microbial species to their community-level diversity.