- Research
- Open access
- Published:
Carry-over effects of urban larval environments on the transmission potential of dengue-2 virus
Parasites & Vectors volume 11, Article number: 426 (2018)
Abstract
Background
Mosquitoes are strongly influenced by environmental temperatures, both directly and indirectly via carry-over effects, a phenomenon by which adult phenotypes are shaped indirectly by the environmental conditions experienced in previous life stages. In landscapes with spatially varying microclimates, such as a city, the effects of environmental temperature can therefore lead to spatial patterns in disease dynamics. To explore the contribution of carry-over effects on the transmission of dengue-2 virus (DENV-2), we conducted a semi-field experiment comparing the demographic and transmission rates of Aedes albopictus reared on different urban land classes in the summer and autumn season. We parameterized a model of vectorial capacity using field- and literature-derived measurements to estimate the bias introduced into predictions of vectorial capacity not accounting for carry-over effects.
Results
The larval environment of different land classes and seasons significantly impacted mosquito life history traits. Larval development and survival rates were higher in the summer than the autumn, with no difference across land class. The effect of land class on adult body size differed across season, with suburban mosquitoes having the smallest wing length in the summer and the largest wing length in the autumn, when compared to other land classes. Infection and dissemination rates were higher in the autumn and on suburban and rural land classes compared to urban. Infectiousness did not differ across land class or season. We estimate that not accounting for carry-over effects can underestimate disease transmission potential in suburban and urban sites in the summer by up to 25%.
Conclusions
Our findings demonstrate the potential of the larval environment to differentially impact stages of DENV-2 infection in Ae. albopictus mosquitoes via carry-over effects. Failure to account for carry-over effects of the larval environment in mechanistic models can lead to biased estimates of disease transmission potential at fine-scales in urban environments.
Background
Climate plays an important role in the transmission of mosquito-borne pathogens, determining the geographic range of disease vectors and shaping transmission dynamics [1, 2]. Heterogeneity in environmental conditions can directly shape individual-level variation in traits relevant to mosquito population dynamics [3] and pathogen transmission [4]. In addition to these direct effects, mosquito phenotypes can be shaped indirectly by the environmental conditions experienced in previous life history stages, a phenomenon known as carry-over effects [5]. Carry-over effects have been documented in a wide-range of species with complex life-cycles, such as amphibians [6], migratory birds [7], and damselflies [8]. Similarly, the mosquito life-cycle is characterized by ontogenetic niche shifts, with a larval aquatic stage and an adult terrestrial stage. Following these studies, we reason that the thermal environment a mosquito experiences during its larval stage is likely to have lasting impacts on adult traits, and, ultimately, on transmission potential.
Although it has been previously demonstrated that larval environmental temperature can alter individual mosquito traits important for transmission [9, 10], the net effect of temperature-mediated carry-over effects on overall transmission potential is ambiguous. Current models of mosquito-borne disease typically only incorporate direct effects of temperature, despite evidence that carry-over effects can have large impacts on adult phenotypes [11,12,13]. Additionally, laboratory studies designed to estimate temperature-mediated carry-over effects are often conducted across a wider range of temperatures than mosquitoes typically experience in the field [14]. The studies are not easily “scaled-up” to explain transmission across a landscape when incorporated into temperature-dependent models of mosquito-borne disease [15]. Urban landscapes, in particular, are composed of a variety of microclimates, which can differentially impact mosquito life-history traits leading to heterogeneity in vector population dynamics across the landscape [16]. However, it is unknown if variation in microclimate across an urban area also has implications for carry-over effects of the larval environment on adult phenotypes.
We hypothesize that relevant environmental variation across an urban landscape during the larval stage will have lasting impacts on adult traits that are important for mosquito population dynamics and pathogen transmission. Further, we predict that failure to account for carry-over effects will result in a biased estimate of vectorial capacity, the rate at which future infections arise from one infectious mosquito. To estimate the effects of the larval environment in a spatially heterogeneous, urban environment, we conducted a semi-field experiment exploring population and dengue-2 virus (DENV-2) transmission relevant life-history traits from Aedes albopictus mosquitoes reared in three urban land classes across the summer and autumn. We used a mixture of field-derived and temperature-dependent parameters to construct a model of vectorial capacity. Our modeled vectorial capacity was then compared to a calculation using the experimental grand mean for parameters affected by carry-over effects in order to estimate the bias introduced by not including these indirect effects.
Methods
We conducted a semi-field experiment across an urban gradient in Athens, GA, USA, in the summer and autumn of 2016. To explore the effects of microclimate variation across an urban landscape, we used an impervious surface map (National Land Cover Database 2011 [17]) to select three replicate sites 30 × 30 m each of low (0–5%), intermediate (6–40%), and high (41–100%) impervious surface. Percent impervious surface is an accurate predictor of land surface temperature, particularly for urban landscapes [18], and allowed us to ensure our sites exhibited the full range of urban microclimates. To select our sites, we calculated the percent impervious surface of each 30 × 30 m pixel using a moving focal window of 210 × 210 m, as the surrounding impervious surface can affect the microclimate in the pixel of interest. We then classified each pixel based on the mean impervious surface within its focal window, with 0–5% representing low, 6–40% representing intermediate, and 41–100% representing high. Because impervious surface is an effective classifier of urban land classes [19], we identified the sites as rural, suburban, and urban with low, intermediate, and high impervious surface scores, respectively. Final site selection was constrained by access and permissions, however, the final distribution of sites was chosen to ensure all sites were at least 3 km from others of the same land class, and were interspersed across the study area (Fig. 1).
Within each site, we evenly distributed four plastic trays (Sterilite, 34.60 × 20.96 × 12.38 cm), each containing 100 first-instar Ae. albopictus larvae and 1l of leaf infusion. Aedes albopictus were from a laboratory colony obtained from the Centers for Disease Control (Atlanta, GA, USA) originating from Keyport, NJ, USA in 1995 (strain ATM-NJ95) [20] and maintained following standardized protocols. Leaf infusion was prepared as described in Murdock et al. [16]. Briefly, 80 g live oak (Quercus virginiana) leaves and 3 g of 1:1 yeast:albumin mixture were infused in deionized water. Trays were screened with a fine mesh, placed in a wire cage to deter wildlife, covered with clear plastic vinyl to keep rainwater from entering, and placed in full shade. We added deionized water to trays after two weeks to maintain a total water volume at 1l. We placed data loggers [Monarch Instruments, Amherst, NH, USA: Radio Frequency Identification (RFID) Temperature Track-It Logger] in vegetation next to each tray, approximately 0.9 m above the ground. Data loggers recorded instantaneous temperature and relative humidity at ten minute intervals throughout the study period. Data loggers were also placed in the trays to measure the larval, aquatic temperature, however three and 17 loggers (of 36) failed due to water damage in the summer and autumn, respectively. Of loggers that did not fail during the experiment, water temperatures were highly correlated with ambient temperatures (ρ = 0.929); thus, only ambient temperatures are used as an approximation of larval environmental temperature. Sites were visited daily to collect emerging adults until all larvae had emerged or died (Summer Replicate: August 1 to September 3, 2016, Autumn Replicate: September 26 to November 8, 2016). We quantified the total number of adults emerging per day, and recorded the sex and wing length of each emerged adult. Adult females were collected to use in vector competence assays.
Dengue virus in vitro culturing and mosquito infections
DENV-2 stock was obtained from the World Reference Center for Emerging Viruses and Arboviruses at the University of Texas Medical Branch (PRS 225 488, originally isolated from human serum in Thailand in 1974 [21]). We propagated virus by inoculating Vero (African green monkey kidney epithelial) cells with a low MOI infection. Virus-containing supernatant was harvested when the cells exhibited more than 80% cytopathic effect. Supernatant was cleared of cell debris by centrifugation (1000× g, 1 min), aliquoted into cryo-vials, and stored at -80 °C. We quantified viral titers of virus stock using TCID-50 assays, calculated by the Spearman-Karber method [22, 23]. When mixed 1:1 with the red blood cell mixture, the final concentration of virus in the blood meal was 3.540 × 106 TCID50/ml.
Adult mosquitoes were collected as they emerged from trays, aggregated by site, and stored in reach-in incubators at 27 ± 0.5 °C, 80 ± 5% relative humidity, and a 12:12 h light:dark photocycle. To ensure infected mosquitoes were of a similar age, mosquitoes were pooled into cohorts of 4–6 days-old in the summer and 4–9 days old in the autumn (due to slower and more asynchronous emergence rates). Mosquitoes were allowed to mate and fed ad libitum with a 10% sucrose solution. Forty-eight hours prior to infection, the sucrose was replaced with deionized water, which was then removed 12–14 h before infection to encourage feeding. Infectious blood meals were administered to mosquitoes through a water-jacketed membrane feeder and consisted of 47% human red blood cells washed in DMEM (v/v), 1% sucrose (w/v), 20% FBS (v/v), 5 mM ATP, and 33% DMEM medium combined with 1 ml of virus stock [24]. Blood-fed female mosquitoes were then maintained as described above for the duration of the experiment.
For a mosquito to become infectious, arboviruses must pass through multiple tissues that impose significant barriers to infection, namely the midgut and salivary glands [25]. Therefore, we assessed mosquitoes for infection, dissemination, and infectiousness through salivation assays and tissue dissections 21 days post-infection [26]. First, mosquitoes were cold-anesthetized and immobilized by removing their legs and wings. Wings were mounted on a glass slide to measure wing length from the distal end of the alula to the apex of the wing via a dissecting scope and micrometer. The proboscis of each female was then inserted into a sterile pipette tip containing 10–20 μl of FBS (with 3 mM ATP and red food coloring) and allowed to salivate on a plate kept at 27 °C for 15 min, after which the salivation media was expelled into 500 μl of DMEM and stored at -80 °C. After salivation, we removed the head of each individual and stored the body and head separately at -80 °C.
To determine variation in the proportion of mosquitoes that become infected (bodies positive for virus), disseminated (heads positive for virus), and infectious (saliva positive for virus), we used cytopathic effect (CPE) assays to test for the presence of virus in each collected tissue [23]. Individual bodies and heads were homogenized in 500 μl of DMEM and centrifuged at 2500 rcf for 5 min. Two-hundred microliters of homogenate was added to Vero cells in a solution of DMEM (1% pen-strep, 5% FBS by volume) in a 24-well plate and kept at 37 °C and 5% CO2. Salivation media was thawed and plated on Vero cells as above. After 5 days, Vero cells were assessed for presence of DENV-2 via CPE assays. Samples were identified as positive for virus if CPE was present in the well.
All infection work was conducted in an arthropod containment level 2 (ACL-2) facility at the University of Georgia in the College of Veterinary Medicine. The physical space as well as experimental protocols have been reviewed and approved by the University of Georgia Office of Biosafety (2015-0038). Briefly, all DENV-2 exposed mosquitoes were counted initially and throughout the experiment, housed in secondary containment cages, and handled in a glove box and on ice when they were removed from secondary containment for forced salivations. All virus assays were also conducted in a biosafety cabinet in our biosafety level II (BSL-2) facility. Finally, we used designated and approved secondary containment to transport virus or infected tissues between our ACL-2 and BSL-2.
Intrinsic growth rates (r’) and vectorial capacity (VC)
We calculated the per capita population growth rate per tray following Livdahl & Sugihara [27] (Eqn. 1):
Following Livdahl & Sugihara [27], we assume N0 to represent the initial number of females before accounting for mortality during the larval stage. This enables the mortality rate to be included via the summed Ax/N0 parameter. Setting N0 equal to the number of emerged mosquitoes would imply a 100% larval survival rate, which was not the case in our study. Unfortunately, we cannot identify the sex of first-instar larvae, and must assume a ratio within the initial cohort. While we do record the proportion of emerged mosquitoes, this represents those that have survived the larval environment until emergence, and, given our findings regarding the effects of the larval microclimate on larval survival, may not be representative of the initial cohort. In our study, males emerged 1–3 days earlier than females in the summer, and up to a week earlier than females in the autumn. This additional time in the larval environment could have exposed female mosquitoes to stressful temperatures and lower resource concentrations than the male mosquitoes that emerged earlier, resulting in lower emergence rates. Eggs used in the experiment were drawn from a laboratory colony which is known to have approximately a 50:50 male:female sex ratio. Therefore, we used this value in (50% of the larvae, 50) our calculations. The other parameters are defined as follows: Ax is the number of mosquitoes emerging on day x, D is the time to reproduction following emergence (assumed to be 14 days [28]), and \( f\left(\overline{w_x}\right) \) is fecundity as a function of mean wing size on day x (\( \overline{w_x} \); Eqn. 2). This relationship is assumed to be linear and calculated via Lounibos et al. [29] (Eqn. 2):
While it is possible to reason how changes in each parameter will result in carry-over effects that individually affect disease transmission, determining the overall net effect and magnitude of the change is less straightforward. Therefore, we calculated the vectorial capacity (VC; Eqn. 3) for each site and season using a modified temperature-dependent dengue calculation defined in Mordecai et al. [30] to create a quantitative estimate of the influence of carry-over effects on disease transmission. Using the experimental mean for field-derived parameters affected by carry-over effects (fecundity and vector competence), we calculated an additional site-level VC to serve as an estimate of this value when not accounting for site-specific carry-over effects.
Here, mosquito traits are a function of temperature (T), as described in Table 1. Site-level VC was calculated using a combination of traits empirically measured in this study and traits estimated from thermal response models as described in Mordecai et al. [30].
The bite rate (a(T)), adult mosquito mortality rate (μ(T)), and extrinsic incubation rate (EIR(T)), were calculated for mosquitoes at a constant 27 °C using temperature-dependent functions from Mordecai et al. [30], to match the adult environment used in the experiment. Vector competence (b(T)c(T)) was calculated as the proportion of infectious mosquitoes per site as determined by our DENV-2 infection assays. Conventionally, vector competence is the product of the proportion of mosquitoes that become infected after biting an infected human and the proportion of bites by infectious mosquitoes that infect humans. Our estimate, the proportion of infectious mosquitoes as measured by CPE assays, is the same as the product of the proportion of mosquitoes that become infected following an infectious blood meal and the proportion of infected mosquitoes that have DENV-2 virus particles in their saliva. With this formulation we are assuming that all infectious bites result in human infection, as we are not directly measuring dengue infection outcomes in humans (i.e. effects of human immunity on DENV infection). The number of eggs produced per female per day (EFD(T)) was calculated by estimating fecundity from average female wing length following Eqn. 2, and then dividing this by the expected lifespan of mosquitoes (1/μ). The egg-to-adult survival probability (pEA(T)) was defined as empirically measured egg-to-adult survival probability (the average proportion of adult female mosquitoes emerging per site). The mosquito immature development rate (MDR(T)) was calculated as the inverse of the mean time to emergence for female mosquitoes per site, resulting in a daily rate of development. To estimate bias introduced by not including carry-over effects, we compared our site-level calculated VC to one calculated using the experimental grand mean for site-level EFD and bc. All other parameters were the same across the two models.
Statistical analysis
We used linear mixed models (LMMs) to explore if microclimate (i.e. mean, minimum, maximum, and daily ranges of temperature and relative humidity), larval development rate (1/days to emergence), female body size, and per capita growth rate differed across land class and season. Egg-to-adult survival (the proportion of adult females emerging per tray) and metrics of vector competence (i.e. infection, dissemination, and infectiousness) were fit using generalized linear mixed models (GLMMs) with binomial distributions and logit links. In all models, fixed effects included land class, season, and their interaction, with site as a random effect. The effect of body size on infection dynamics was also explored at the level of the individual mosquito, fitting a binomial GLMM including wing size as a fixed effect and site as a random effect. Vectorial capacity was calculated at the site-level, and so did not require site to be included as a random effect. We therefore used a regression model to estimate the effect of land class, season, and their interaction on site-level vectorial capacity.
To confirm the relationship between the categorical variables of land use and season and temperature, we fit additional models containing mean temperature as a covariate to the residuals of the original models including season and land use as fixed effects. This test explored if there was additional variation in the response variable due to temperature that was not explained by land class and season. To explore if the effect of temperature differed across season, we fit individual models to the above response variables including mean temperature as averaged across each season (e.g. summer or autumn) as a covariate, using the same distributions and link functions. For egg-to-adult survival, larval development, body size, and the per capita growth rate, mean temperature was calculated over each individual season (e.g. summer and autumn) at the tray level, and site was included as a random effect. Because mosquitoes were pooled by site for infection assays, temperature was aggregated to the site level and no random effects were included for analyses of infection metrics and VC.
All analyses were conducted with respect to the female subset of the population, as they are the subpopulation responsible for disease transmission. In the case of data logger failure (n = 3), imputed means from the site were used to replace microclimate data. In the case of trays failing due to wildlife tampering (two urban and one suburban in the autumn replicate), collected mosquitoes were used for infection assays, but trays were excluded from demographic analyses. For all mixed-models, significance was assessed through Wald Chi-square tests (α = 0.05) and examination of 95% confidence intervals. Pearson residuals and Q-Q plots were visually inspected for normality. All mixed models were fit using the lme4 [31] package in R v. 3.5.0 [32]. Code to run analyses and create figures is deposited on figshare.
Results
Effects of land class and season on microclimate
We found that microclimate profiles differed significantly across both season and land class (Fig. 2, Table 2). In general, temperatures were warmer in the summer and on urban sites, replicating what was found in a prior study in this system [16]. We did observe a significant interaction between season and land use on the mean daily minimum temperature and diurnal temperature range, with no effect of land use on these response variables in the summer. Urban sites in the autumn were characterized by significantly higher daily average minimum temperature and smaller diurnal temperature range relative to rural sites (Table 2). Mean relative humidity was higher in the summer than the autumn: [mean (95% CI)], summer: 87.93% (86.33–89.54 %); autumn: 73.32% (71.72– 74.92%). In the summer, minimum and mean relative humidity was significantly lower on urban sites compared to rural and suburban sites (Table 2). A similar trend was seen in the autumn, with urban sites having lower mean relative humidity compared to other land classes, but no difference in minimum relative humidity (Table 2).
Direct and carry-over effects of land class and season on population growth
Of the 3600 first-instar larvae placed in each season, a total of 2595 and 1128 mosquitoes emerged in the summer and autumn, respectively. The total female egg-to-adult survival per tray was significantly higher in summer than autumn [Table 3, mean (binomial asymptotic 95% CI), summer: 0.670 (0.598–0.735); autumn: 0.297 (0.235–0.366)], but did not differ across land class (Fig. 3a, Table 3). The mean rate of larval development per tray was significantly different between summer and autumn (Fig. 3b, Table 3), with daily mean ± SE development rates of 0.074 ± 0.002 and 0.0387 ± 0.002, respectively. There were no significant differences in larval survival or development rates across land class. We did not observe a significant carry-over effect of land class or season on mosquito wing size, however there was a significant interaction between the two (Table 3). We found a significant difference in wing size across season for mosquitoes on rural sites only, with larger bodied mosquitoes in the summer (mean ± SD 2.451 ± 0.211 mm), than the autumn (2.300 ± 0.202 mm). While urban mosquitoes tended to be larger in the autumn, and suburban mosquitoes tended to be larger in the summer, these effects were not significant.
After incorporating the number of adult females emerging per day, the day of emergence, and their body size into the per capita growth rate equation (Eqn. 1), we found that the estimated per capita growth rate was higher in the summer season than the autumn season (Fig. 3c, Table 3, mean ± SE, summer: 0.135 ± 0.005; autumn: 0.068 ± 0.006) with no difference across land class. The effect of temperature within a season was only significant for egg-to-adult survival, and differed in direction across season (mean β ± SE, summer: -0.328 ± 0.148; autumn: 0.368 ± 0.135, Additional file 1: Table S1). This mirrors a trend for the effect of land class on egg-to-adult survival to differ across season (Table 3). When controlling for land class and season, temperature explained no additional variation for any response variable (Additional file 1: Table S2).
Carry-over effects of land class and season on vector competence
A total of 319 female mosquitoes were assessed for infection status, 20 per site in the summer and varying numbers per site in the autumn due to lower emergence rates (sample sizes reported in Table 4). Carry-over effects of the larval environment on infection status were limited to infection and dissemination rates. We found that land class and season did significantly impact the probability of a mosquito becoming infected and disseminating dengue infection (Table 3). Both metrics were higher in the autumn compared to the summer replicate, with urban sites having the lowest infection and dissemination rates across both seasons (Fig. 4a, b). While there was a trend for a higher proportion of mosquitoes becoming infectious in the summer (Fig. 4c), this was not significant (χ2 = 3.63, P = 0.057). The probability of becoming infectious did not differ across land class, nor season (Fig. 4c, Table 3), despite the higher probability of mosquito infection and dissemination in the autumn, and on suburban and rural sites. Similarly, there was no effect of temperature on any infection metric within a season (Additional file 1: Table S1), and temperature did not explain any additional variation after controlling for land class and seasons (Additional file 1: Table S2). This suggests that the ability of virus to escape the midgut and invade the salivary glands differs in adults reared in the summer vs the autumn and across land class, with a higher proportion of dengue infected mosquitoes becoming infectious in the summer and on urban sites (Table 4, (χ2 = 13.65, P < 0.001). We also found the probability of infection to decline with increasing body size (χ2 = 4.776, P = 0.0289, although there was no evidence for a relationship between body size and the probability of dissemination or infectiousness.
Integrating direct and carry-over effects into estimates of transmission potential
We found VC to be higher in the summer (mean ± SE: 5.847 ± 0.768) than the autumn (0.252 ± 1.097) (Fig. 5, Table 3). In the summer season, there was a trend for VC to increase with increasing urbanization (Fig. 5). This trend was not significant, however, given the small sample size (n = 9) and the disproportional impact of having no infectious mosquitoes at one site, resulting in a value of VC = 0 for one sample. There was no effect of temperature on VC within a season (Additional file 1: Table S1), and temperature did not explain any additional variation after controlling for land class and season. When comparing VC calculations using field-based or grand mean estimates of EFD and bc, we found that the effect of land class and season were not significantly different (land class: χ2 = 0.381, P = 0.826), season: χ2 = 1.408, P = 0.235), suggesting that the omission of carry-over effects in calculations did not lead to biased estimates of relative VC in different seasons or land classes. However, the use of the grand mean did lead to an underestimate of VC on some suburban and urban sites in the summer, with a two-fold decrease in predicted VC (Fig. 5, Additional file 2: Figure S1). The calculated VC for rural sites in the summer and across all land classes in the autumn more closely resembled the grand mean calculated VC.
Discussion
Mathematical models of mosquito-borne disease rarely include mosquito larval stages [15], and of those that do, few include the influence of carry-over effects on important mosquito life-history traits (but see [33]). This is likely because there are relatively few empirical studies parameterizing carry-over effects in mosquito-pathogen systems [2], and most are laboratory studies conducted across a wider range of temperatures than those seen in the field. Here, we demonstrate that fine-scale differences in larval microclimate across land class and season generate carry-over effects on adult fecundity and vector competence for DENV-2. When integrated into a model of vectorial capacity, we find that vectorial capacity differs across season, but not land class. Further, failure to account for site-specific carry-over effects across urban land classes results in biased estimates of DENV-2 transmission potential, underestimating potential disease transmission in urban areas.
The subtle heterogeneity in microclimate we observed across season resulted in significantly different predicted population growth rates through its effects on demographic traits. Daily mean temperatures (25.43 °C) across all sites in the summer were closer to the predicted thermal optimum of Ae. albopictus (24–25 °C) [30] than in the autumn (17.69 °C), leading to higher egg-to-adult survival rates. We also observed more rapid larval development rates in the summer relative to the autumn. This is likely due to the strong positive relationship observed between development rates and mean larval temperature, as the metabolic rate of mosquitoes will increase with warming temperatures [3]. Temperature explained no additional variation in any response variable after accounting for land class and season, suggesting that our coarser characterizations of land class and season contain the temperature variation necessary to predict changes in demographic and infection rates. Additionally, we only found an effect of temperature within a season for egg-to-adult survival (Additional file 1: Table S1). While we did not find a significant influence on many traits, our trends do agree with a previous study in this system that found lower egg-to-adult survival on urban sites [16]. The variation in mean temperature across land class in our study was very small (< 1 °C), and we expect these relationships would be magnified in mega-cities that can have urban heat island effects of up to 6 °C [34].
Surprisingly, we found no main effect of land class or season on female mosquito body size, despite the difference in temperatures across season. Following allometric temperature-size relationships of ectotherms, warmer larval temperatures should lead to smaller-bodied mosquitoes [35]. However, contrary to predictions generated from the allometric temperature-size relationship, we observed mosquitoes on rural sites to be larger in the summer despite the fact that all land classes were cooler in the autumn relative to the summer. Our results contrast with many laboratory studies that have found a negative relationship between rearing temperature and mosquito body size (Ae. albopictus [36], Culex tarsalis [37], Anopheles gambiae [38]). However, these studies all used a constant temperature treatment, while mosquitoes in our field-based study experienced fluctuating temperatures. Among studies using fluctuating temperatures, there is mixed evidence for a relationship between rearing temperature and mosquito body size [16, 39]. Larger temperature fluctuations at the more extreme temperatures (cool and warm) can lead to counterintuitive effects of temperature on organismal traits if these temperatures approach or cross the thermal maximum or minimum (at which trait performance is zero) and induce thermal stress [40, 41]. Rural sites in the autumn did experience a larger average diurnal range of temperatures than in the summer, suggesting this differential effect of temperature fluctuations at thermal extremes could be acting on body size. Our findings demonstrate that, while the use of fluctuating temperatures in studies of mosquito life-history traits is relatively new, these fluctuations can have significant impacts on mosquito ecology and should be integrated in laboratory-based studies of mosquito vectors to more closely approximate field conditions.
Our results agree with laboratory studies in other arboviral systems (chikungunya [42], yellow fever [42], and Rift Valley fever [43]) that found cool larval environmental temperatures to enhance arbovirus infection relative to warmer larval environments. Studies in the Ae. albopictus-dengue virus system have also found that low larval temperatures enhance mosquito susceptibility to viral infection, although this is dependent on larval nutrition [10] and the stage of the infection (i.e. midgut vs dissemination vs saliva) [9]. While we found infection and dissemination of DENV-2 to decrease with increasing temperatures across season and land class, there was no effect on viral presence in the saliva, suggesting carry over effects due to microclimate variation may alter the overall efficiency of dengue infection. Thus, even though a smaller proportion of mosquitoes reared on urban sites and in the summer became infected and disseminated infection, these mosquitoes were more likely to become infectious, resulting in no net difference in overall vector competence across land class and season. Larval environmental temperature may differentially impact later stages of viral infection (i.e. salivary gland penetration) compared to earlier stages (i.e. midgut escape) through effects on mosquito physiology and immunity, as well as on important tissue barriers to infection [4, 42, 44, 45]. Further, our study considered only DENV-2, and other arboviruses and mosquito-borne disease are likely influenced by the mosquito’s larval environment differently.
Current models of vector-borne disease focus primarily on direct effects of environmental variables on mosquito densities and disease transmission and rarely include the effects of the larval stage, either directly or via carry-over effects [15]. While we found carry-over effects due to seasonal and urban environments to have a significant impact on DENV-2 infection and dissemination, we found no net effects on saliva positivity for the virus. Therefore, when incorporating parameters into calculations of vectorial capacity, we did not find a significant difference in predicted vectorial capacity due to land class. However, we did find VC to be higher in the summer relative to the autumn, driven by differences in demographic rates such as larval survival and development rates, rather than differences in adult vector competence. Unfortunately, given the logistical limitations imposed by a field experiment setting, we were unable to measure additional life-history traits important for disease transmission in conjunction with vector competence. Laboratory studies have found that factors such as adult longevity [46], biting rate [47], and pathogen extrinsic incubation period [48, 49] are also be impacted by carry-over effects. For example, warmer larval temperatures correspond with decreased adult longevity in mosquitoes [46], and including this relationship could mediate the seasonal differences in VC found in our study, with decreased adult longevity in the summer corresponding to decreased VC. Less is known about traits specific to transmission such as biting rate and EIP, which have only been investigated in response to larval diet and competition [47,48,49]. Carry-over effects of the larval environment can act on multiple adult phenotypes, often in conflicting ways, and the net effect of this on disease transmission has yet to be fully explored.
Our study was further limited by the difficulties in obtaining appropriate sample sizes. While semi-field experiments incorporate more realistic variation in environmental temperature than laboratory experiments, they require additional space and travel time in order to distribute replicates in a manner that meets assumptions of independence across sites. Given the size of our study area, nine was the maximum number of sites that it was possible to visit daily. Unexpectedly low emergence rates of mosquitoes in the autumn further reduced the sample size of mosquitoes that could be used in infection assays. Despite this limitation, we did find significant differences in mosquito demographic rates across season and in infection and dissemination rates across land class, suggesting that site-specific characteristics can directly and indirectly impact vector-borne disease dynamics. Yet, due to the low replication across sites, these results must be interpreted conservatively.
Carry-over effects are not simply limited to microclimate, and can result due to variation in larval nutrition [47], intra- and interspecific densities [50], and predation [33] in mosquito systems. Further, abiotic and biotic factors will likely interact to influence carry over effects [10, 51], and this interaction could be scale-dependent [52]. For example, biotic processes are predicted to be more important at local geographical scales, while abiotic processes dominate at regional geographical scales in species distribution models [53]. Future exploration of the scale-dependent contribution of different environmental factors and their interactive influence on both direct and carry-over effects is needed to improve models predicting the distribution of mosquito vector species, mosquito population dynamics and disease transmission.
Conclusions
We found fine-scale variation in microclimate across season and urban land class to shape Ae. albopictus population dynamics and DENV-2 transmission potential through direct effects on larval survival and development rates, and indirectly through carry-over effects on vector competence and fecundity. Although sample sizes were limited, our study indicates the potential effects that site-specific environments can have on mosquito demographics and infection dynamics. DENV-2 infection and dissemination rates were higher in mosquitoes from rural and suburban land classes than urban ones, and were higher in the autumn compared to the summer. However, there was no difference in overall infectiousness. Therefore, the seasonal differences in VC we observed were due to the direct effects of the larval environment on egg-to-adult survival and development rates, rather than carry-over effects. When comparing VC to a calculated VC that did not account for site-specific carry-over effects, we found that not accounting for carry-over effects results in an underestimate of predicted VC in suburban and urban sites in the summer, and an overestimate in the autumn. The interaction between the larval and adult environments, mediated by carry-over effects, could have complex consequences for adult phenotypes relevant to disease transmission for mosquitoes as well as other organisms. Given the devastating impact of disease in other species with complex life histories (e.g. chytridiomycosis in amphibians), carry-over effects in disease transmission are important, though understudied, mechanisms that must be better understood to control disease spread. Incorporating relationships between carry-over effects and organismal life-history traits into statistical and mechanistic models will lead to more accurate predictions on the distributions of species, population dynamics, and the transmission of pathogens and parasites. Mosquito-borne disease incidence is spatially heterogeneous in urban areas [54], and a better understanding of both the larval and adult environments, including their interaction, could improve the accuracy of fine-scale predictions of disease incidence across a city.
References
Reiter P. Climate change and mosquito-borne disease. Environ Health Perspect. 2001;109:141–61.
Parham PE, Waldock J, Christophides GK, Hemming D, Agusto F, Evans KJ, et al. Climate, environmental and socio-economic change: weighing up the balance in vector-borne disease transmission. Proc R Soc B Biol Soc. 2015;370:20130551.
Delatte H, Gimonneau G, Triboire A, Fontenille D. Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009;46:33–41.
Murdock C, Paaijmans K, Bell A, King J, Hillyer J, F Read A, et al. Complex effects of temperature on mosquito immune function. Proc R Soc B Biol Soc. 2012;279:3357–66.
Harrison XA, Blount JD, Inger R, Norris DR, Bearhop S. Carry-over effects as drivers of fitness differences in animals. J Anim Ecol. 2011;80:4–18.
Vonesh JR. Sequential predator effects across three life stages of the African tree frog, Hyperolius spinigularis. Oecologia. 2005;143:280–90.
Norris DR, Taylor CM. Predicting the consequences of carry-over effects for migratory populations. Biol Lett. 2006;2:148–51.
De Block M, Stoks R. Fitness effects from egg to reproduction: bridging the life history transition. Ecology. 2005;86:185–97.
Alto BW, Bettinardi D. Temperature and dengue virus infection in mosquitoes: independent effects on the immature and adult stages. Am J Trop Med Hyg. 2013;88:497–505.
Buckner EA, Alto BW, Lounibos LP. Larval temperature-food effects on adult mosquito infection and vertical transmission of dengue-1 virus. J Med Entomol. 2016;53:91–8.
Muturi EJ, Lampman R, Costanzo K, Alto BW. Effect of temperature and insecticide stress on life-history traits of Culex restuans and Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2011;48:243–50.
Muturi EJ, Kim C-H, Alto BW, Berenbaum MR, Schuler MA. Larval environmental stress alters Aedes aegypti competence for Sindbis virus. Trop Med Int Health. 2011;16:955–64.
Price DP, Schilkey FD, Ulanov A, Hansen IA. Small mosquitoes, large implications: crowding and starvation affects gene expression and nutrient accumulation in Aedes aegypti. Parasit Vectors. 2015;8:252.
Cator LJ, Thomas S, Paaijmans KP, Ravishankaran S, Justin JA, Mathai MT, et al. Characterizing microclimate in urban malaria transmission settings: a case study from Chennai, India. Malaria J. 2013;12:1–1.
Reiner RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J R Soc Interface. 2013;10:20120921.
Murdock CC, Evans MV, McClanahan TD, Miazgowicz KL, Tesla B. Fine-scale variation in microclimate across an urban landscape shapes variation in mosquito population dynamics and the potential of Aedes albopictus to transmit arboviral disease. PLoS Negl Trop Dis. 2017;11:e0005640.
Xian GZ, Homer CG, Dewitz J, Fry J, Hossain N, Wickham J. Change of impervious surface area between 2001 and 2006 in the conterminous United States. Photogramm Eng Remote Sensing. 2011;77:5.
Yuan F, Bauer ME. Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery. Remote Sens Environ. 2007;106:375–86.
Lu D, Weng Q. Use of impervious surface in urban land-use classification. Remote Sens Environ. 2006;102:146–60.
Marcombe S, Farajollahi A, Healy SP, Clark GG, Fonseca DM. Insecticide resistance status of United States populations of Aedes albopictus and mechanisms involved. PLoS One. 2014;9:e101992.
Vazeille-Falcoz M, Mousson L, Rodhain F, Chungue E, Failloux AB. Variation in oral susceptibility to dengue type 2 virus of populations of Aedes aegypti from the islands of Tahiti and Moorea, French Polynesia. Am J Trop Med Hyg. 1999;60:292–9.
Shao Q, Herrlinger S, Yang S-L, Lai F, Moore JM, Brindley MA, et al. Zika virus infection disrupts neurovascular development and results in postnatal microcephaly with brain damage. Development. 2016;143:4127–36.
Willard KA, Demakovsky L, Tesla B, Goodfellow FT, Stice SL, Murdock CC, et al. Zika virus exhibits lineage-specific phenotypes in cell culture, in Aedes aegypti mosquitoes, and in an embryo model. Viruses. 2017;9:383.
Shan C, Xie X, Muruato AE, Rossi SL, Roundy CM, Azar SR, et al. An infectious cDNA clone of Zika virus to study viral virulence, mosquito transmission, and antiviral inhibitors. Cell Host Microbe. 2016;19:1–23.
Cheng G, Liu Y, Wang P, Xiao X. Mosquito defense strategies against viral infection. Trends Parasitol. 2016;32:177–86.
Anderson SL, Richards SL, Smartt CT. A simple method for determining arbovirus transmission in mosquitoes. J Am Mosq Control Assoc. 2010;26:108–11.
Livdahl TP, Sugihara G. Non-linear interactions of populations and the importance of estimating per capita rates of change. J Anim Ecol. 1984;53:573–80.
Livdahl TP, Willey MS. Prospects for an invasion: competition between Aedes albopictus and native Aedes triseriatus. Science. 1991;253:189–91.
Lounibos LP, Suarez S, Menendez Z, Nishimura N, Escher RL, O’Connell SM, et al. Does temperature affect the outcome of larval competition between Aedes aegypti and Aedes albopictus? J Vector Ecol. 2002;27:86–95.
Mordecai EA, Cohen JM, Evans MV, Gudapati P, Johnson LR, Lippi CA, et al. Detecting the impact of temperature on transmission of Zika, dengue and chikungunya using mechanistic models. PLoS Negl Trop Dis. 2017;11:e0005568.
Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2018. https://www.R-project.org/
Roux O, Vantaux A, Roche B, Yameogo KB, Dabiré KR, Diabaté A, et al. Evidence for carry-over effects of predator exposure on pathogen transmission potential. Proc R Soc B. 2015;282:20152430.
Peng S, Piao S, Ciais P, Friedlingstein P, Ottle C, Bréon F-M, et al. Surface urban heat island across 419 global big cities. Environ Sci Technol. 2012;46:696–703.
Angilleta MJ, Steury TD, Sears MW. Temperature, growth rate, and body size in ectotherms: fitting piece of a life-history puzzle. Integr Comp Biol. 2004;44:498–509.
Reiskind MH, Zarrabi AA. Is bigger really bigger? Differential responses to temperature in measures of body size of the mosquito, Aedes albopictus. J Insect Physiol. 2012;58:911–7.
Dodson BL, Kramer LD, Rasgon JL. Effects of larval rearing temperature on immature development and West Nile virus vector competence of Culex tarsalis. Parasit Vectors. 2012;5:199.
Koella JC, Lyimo EO. Variability in the relationship between weight and wing length of Anopheles gambiae (Diptera: Culicidae). J Med Entomol. 1996;33:261–4.
Mohammed A, Chadee DD. Effects of different temperature regimens on the development of Aedes aegypti (L.) (Diptera: Culicidae) mosquitoes. Acta Trop. 2011;119:38–43.
Carrington LB, Armijos MV, Lambrechts L, Barker CM, Scott TW. Effects of fluctuating daily temperatures at critical thermal extremes on Aedes aegypti life-history traits. PLoS One. 2013;8:e58824.
Colinet H, Sinclair BJ, Vernon P, Renault D. Insects in fluctuating thermal environments. Annu Rev Entomol. 2015;60:123–40.
Adelman ZN, Anderson MAE, Wiley MR, Murreddu MG, Samuel GH, Morazzani EM, et al. Cooler temperatures destabilize RNA interference and increase susceptibility of disease vector mosquitoes to viral infection. PLoS Negl Tropl Dis. 2013;7:e2239.
Turell M. Effect of environmental temperature on the vector competence of Aedes taeniorhynchus for Rift Valley fever and Venezuelan equine encephalitis viruses. Am J Trop Med Hyg. 1993;49:672–6.
Franz AWE, Kantor AM, Passarelli AL, Clem RJ. Tissue barriers to arbovirus infection in mosquitoes. Viruses. 2015;7:3741–67.
Houk E, Obie F, Hardy J. Peritrophic membrane formation and the midgut barrier to arboviral infection in the mosquito, Culex tarsalis Coquillett (Insecta, Diptera). Acta Trop. 1979;36:39–45.
Christiansen-Jucht C, Parham PE, Saddler A, Koella JC, Basáñez M-G. Temperature during larval development and adult maintenance influences the survival of Anopheles gambiae s.s. Parasit Vectors. 2014;7:489.
Moller-Jacobs LL, Murdock CC, Thomas MB. Capacity of mosquitoes to transmit malaria depends on larval environment. Parasit Vectors. 2014;7:593.
Shapiro LLM, Murdock CC, Jacobs GR, Thomas RJ, Thomas MB. Larval food quantity affects the capacity of adult mosquitoes to transmit human malaria. Proc R Soc B Biol Soc. 2016;283:20160298-8
Bara J, Rapti Z, Cáceres CE, Muturi EJ. Effect of larval competition on extrinsic incubation period and vectorial capacity of Aedes albopictus for dengue virus. PLoS One. 2015;10:e0126703r.
Alto BW, Lounibos LP, Higgs S, Juliano SA. Larval competition differentially affects arbovirus infection in Aedes mosquitoes. Ecology. 2005;86:3279–88.
Muturi EJ, Blackshear M, Montgomery A. Temperature and density-dependent effects of larval environment on Aedes aegypti competence for an alphavirus. J Vector Ecol. 2012;37:154–61.
Leisnham PT, LaDeau SL, Juliano SA. Spatial and temporal habitat segregation of mosquitoes in urban Florida. PLoS One. 2014;9:391655.
Cohen JM, Civitello DJ, Brace AJ, Feichtinger EM, Ortega CN, Richardson JC, et al. Spatial scale modulates the strength of ecological processes driving disease distributions. Proc Natl Acad Sci USA. 2016;113:E3359–64.
Santos-Vega M, Bouma MJ, Kohli V, Pascual M. Population density, climate variables and poverty synergistically structure spatial risk in urban malaria in India. PLoS Negl Trop Dis. 2016;10:e0005155.
Acknowledgements
We thank members of the Murdock and Brindley labs for discussion and technical support conducting viral assays. We thank Diana Diaz, Abigail Lecroy, and Marco Notarangelo for assistance in the field and lab. We also thank two anonymous reviewers for their insightful comments throughout the review process.
Funding
This work was supported by the University of Georgia (Presidential Fellowship, College of Veterinary Medicine, Department of Infectious Diseases) the National Science Foundation Graduate Research Fellowship, and the National Science Foundation Research Experiences for Undergraduates (Grant No. 1156707). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets and code used in during the current study are available in the figshare repository, doi: 10.6084/m9.figshare.5558128.
Author information
Authors and Affiliations
Contributions
MVE, JMD and CCM designed the experiment. MVE, JCS and NS conducted the field-work portion of the experiment. MVE and MAB designed and conducted the infection portion of the experiment. MVE, JMD and CCM conducted statistical analyses. MVE, JMD and CCM were involved in original draft preparation and all authors were involved in reviewing and editing. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Table S1. Effect of temperature within a season. Model results from GZLMs estimating the effect of temperature within a season. Binomial models were fit with a logit-link function. Except for those models predicting infection metrics and vector competence, site was included as a random effect. Table S2. Additional variation in residuals explained by temperature. Model results from fitting temperature to residuals of original models (land class × season) for each response variable. In all models, temperature did not explain any additional variation, as evidenced by low mean sum of squares and F-statistics. (PDF 18 kb)
Additional file 2:
Figure S1. Bias in VC due to not accounting for site-level carry-over effects across land class and season. (PDF 7 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Evans, M.V., Shiau, J.C., Solano, N. et al. Carry-over effects of urban larval environments on the transmission potential of dengue-2 virus. Parasites Vectors 11, 426 (2018). https://doi.org/10.1186/s13071-018-3013-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-018-3013-3