- Open Access
Age-structured vectorial capacity reveals timing, not magnitude of within-mosquito dynamics is critical for arbovirus fitness assessment
Parasites & Vectors volume 13, Article number: 310 (2020)
Transmission dynamics of arboviruses like Zika virus are often evaluated by vector competence (the proportion of infectious vectors given exposure) and the extrinsic incubation period (EIP, the time it takes for a vector to become infectious), but vector age is another critical driver of transmission dynamics. Vectorial capacity (VC) is a measure of transmission potential of a vector-pathogen system, but how these three components, EIP, vector competence and vector age, affect VC in concert still needs study.
The interaction of vector competence, EIP, and mosquito age at the time of infection acquisition (Ageacquisition) was experimentally measured in an Aedes aegypti-ZIKV model system, as well as the age-dependence of probability of survival and the willingness to bite. An age-structured vectorial capacity framework (VCage) was then developed using both EIPMin and EIPMax, defined as the time to first observed minimum proportion of transmitting mosquitoes and the time to observed maximum proportion of transmitting mosquitoes.
The within-mosquito dynamics of vector competence/EIP were not significant among treatments where mosquitoes were exposed at different ages. However, VCage revealed: (i) age-dependence in vector-virus interactions is important for transmission success; (ii) lower vector competence but at shorter EIPs was sufficient for transmission perpetuation; and (iii) R0 may be overestimated by using non-age-structured VC.
The results indicate that ultimately the temporal component of the virus-vector dynamics is most critical, especially when exposure occurred at advanced mosquito age. While our study is limited to a single virus-vector system, and a multitude of other factors affect both vector competence and mosquito mortality, our methods can be extrapolated to these other scenarios. Results indicate that how ‘highly’ or ‘negligibly’ competent vectors are categorized may need adjustment.
The transmission dynamics of arboviruses such as Zika virus (ZIKV) are evaluated over several characteristics, namely vector competence and the extrinsic incubation period (EIP). Aedes aegypti is the primary vector of ZIKV and several studies have evaluated its competence to transmit the virus [1,2,3,4]. Vector competence is the ability of a mosquito to acquire and ultimately transmit a virus [5, 6]. The time it takes for this process to occur is referred to as the extrinsic incubation period (EIP) . Vector competence and EIP are interrelated measures of the proportion of vectors that become infectious given exposure and the time it takes for a vector to become infectious given exposure, respectively [6, 8]. Thus, EIP can be described as the temporality of vector competence and has been used to evaluate the relative fitness of arbovirus systems [6, 8, 9]. In addition, EIP and vector competence are influenced by many parameters including vector species, mosquito-immune system, microbiota fauna, discrete populations within species, and environmental factors [2, 3, 5, 10,11,12,13,14,15].
Indeed, changes in arbovirus fitness and thus transmission dynamics have been predicated on altered vector competence, especially as a critical component of vectorial capacity [6, 16,17,18,19,20,21,22,23,24]. The composite of vector competence and EIP into a single, dynamic measure allows for a more comprehensive understanding of this process [5, 10, 25,26,27]. Not all mosquitoes that are exposed will be able to transmit (vector competence) and the time it takes for those mosquitoes that will transmit is not a constant (EIP), and so understanding this composite over several days post-infection is critical [6, 28].
Vectorial capacity (VC) was derived as a measure of transmission potential of a vector-borne pathogen by a competent vector, and incorporates both vector competence and EIP [6, 7, 17, 18, 25, 29]. VC is the vector-centric component of the basic reproduction number (R0), and VC represents the number of secondary cases resulting from the introduction of a single infectious human individual per infectious day of that human index case [30, 31]. VC is given by:
where ‘m’ is the density of mosquitoes relative to humans, ‘a’ is the biting rate, ‘b’ is the vector competence of the mosquito for a particular virus, ‘N’ is the extrinsic incubation period, and ‘p’ is the daily probability of mosquito survival . And R0 is given by:
where ‘r’ is the human recovery rate, or the reciprocal of the average human infectious period, and ‘c’ is the probability of transmission from the human to the vector given contact. R0 is a threshold parameter, whereby an R0 value < 1 is not expected to sustain transmission and an R0 value ≥ 1 is likely to result in an outbreak [11, 32].
VC uses an average daily probability of survival to calculate the probability of a mosquito living through the EIP (pN) and pN/(− ln(p) represents the expected infectious days given N and p . It is intuitive how mortality could impact transmission, but most studies address mortality based on a constant age at which the mosquito acquires the infection [9, 10, 12, 34, 35]. Mathematical modeling studies have addressed age structure with respect to EIP or bite-structure, but no studies have empirically investigated all three processes simultaneously [9, 36]. Thus, herein, we experimentally test the hypothesis that transmission dynamics of arboviruses are as or more affected by time as a function of age versus time as a function of days post-infection. We further develop an age-structured vectorial capacity equation (VCage) to quantify these potential effects.
We first wanted to determine if age and/or prior blood meals affected the within-mosquito viral dynamics, as well as various life traits of the mosquito. To that end, we designated three treatment groups: YOUNG, OLDER and S.OLDER. The YOUNG group was offered an infectious blood meal at approximately 5 days post-emergence (dpe). The OLDER group was offered a non-infectious blood meal at 5 dpe and then an infectious blood meal 1 week later (12 dpe). The S.OLDER group was not offered a prior, non-infectious blood meal, but a single infectious blood meal at approximately 12 dpe (to match the OLDER timing). S refers to sugar, which was done in order to distinguish the absence of a blood meal at 5 dpe. All non-infectious blood meals are referred to as “mock” blood meals, as they contain non-infectious cell culture media in place of infectious media.
Virus and mosquitoes
ZIKV strain PRVABC59 (Asian lineage), originally isolated from a human patient in Puerto Rico in 2015, was provided by Dr Barbara Johnson at the USA Centers for Disease Control and Prevention (GenBank: KU501215). The viral stock was passaged three times in Vero cells before passage on Vero cells for mosquito exposure. Supernatant was collected 3 days post-inoculation and titer determined as previously described . The titer of ZIKV was verified by qRT-PCR at approximately 8 × 107 plaque forming units (pfu)/ml, matched across all exposure experiments. Virus used for mosquito exposure was never frozen, as this has been shown to negatively affect vector competence when compared to freshly cultured virus [38, 39]. Colony Aedes aegypti (Rockefeller) were provided by Dr Daniel Swale of the LSU Entomology Department, Baton Rouge, Louisiana, USA. To isolate the effect(s) of age, mosquitoes were maintained at constant conditions with 16:8 h light:dark periods, approximately 80% humidity (monitored using a digital RH monitor), and at 28 °C constant temperature. Mosquitoes were supplied with a 10% sucrose water solution after emergence via soaked cotton pledgets, which were replenished every 24 h. Sucrose solution was removed 24 h before experiments and was again provided after blood feeding/mock starvation.
Blood-feeding and oral exposure of Ae. aegypti
Infectious and non-infectious blood meals were prepared with ZIKV-infected cell culture supernatant and non-infected cell culture supernatant, respectively. Whole bovine blood in Alsever’s solution (Hemostat Labs, Dixon, CA, USA) was used in a 2:1 blood to supernatant ratio . Mosquitoes were fed via the Hemotek membrane feeding system (Hemotek, Blackburn, UK) for 45 min, after which mosquitoes were cold anesthetized and blood-fed females were placed into clean cartons until further treatment. Because starvation and cold anesthetization could affect mortality rates, every treatment was subjected to the same schedule of each of these conditions. For example, all groups were starved 24 h prior to days 5 and 12 post-emergence, regardless of whether that group would receive a blood meal. In addition, not only were blood-fed mosquitoes cold anesthetized and moved to a new carton, but other groups that did not feed on that schedule were treated the same in order to mimic these conditions and control for these factors. The experimental design is depicted in Additional file 1: Figure S1.
In order to determine if the timing of the infectious blood meal affected the within-mosquito viral kinetics, mosquitoes were sampled at 5, 8 and 11 days post-infection (dpi) to test for infection and dissemination across all groups (sample sizes are provided in Additional file 1: Table S1). Mosquito legs and bodies were placed into separate tubes containing 900 μl BA-1 diluent media and BBs. Samples were then homogenized twice at 25 Hz for 3 min using a Qiagen Tissuelyzer (Qiagen, Hilden, Germany). RNA was extracted (5× MagMax 96 Viral RNA Isolation Kit; Applied Biosystems, Foster City, CA, USA) and tested for the presence of viral RNA via qRT-PCR (SuperScript III One-Step RT-PCR System; Invitrogen, Carlsbad, CA, USA) as in . Each treatment was repeated a total of three times and data was averaged over these replicates. Differences among the treatment groups for infection and dissemination rates were tested by a chi-square test for multiple proportions on 5, 8 and 11 dpi.
To directly assess transmission potential, 20 mosquitoes per treatment were force-salivated. Briefly, ZIKV-exposed mosquitoes were immobilized on ice before removing legs and wings. Mosquitoes were then placed on double-sided tape, and the proboscis of each mosquito was placed into a pipette tip containing 35 μl FBS with 3 mmol/l ATP for 30 min, as previously described . Contents of the pipette tip were ejected into 100 μl BA-1 diluent and stored at − 80 °C before testing (see below).
Mosquitoes from the two OLDER groups were sampled at 5, 8 and 11 dpi, as well as 16 dpi to represent the end of the mortality study (28 dpe). In order to characterize the EIP across the mosquito’s lifetime, mosquitoes in group YOUNG were sampled at 5, 8 and 11 dpi, as well as additional days post-infection in order to match the age at which the older groups were sampled. Thus, by adding 7 to the time points in the older groups, we additionally sampled the YOUNG group at 12, 15 and 18 dpi, as well as 23 dpi which corresponds to the oldest time point of 28 days old. Samples were tested for the presence of viral RNA in saliva using the techniques described above [40, 41]. Differences among the treatment groups’ transmission rates were tested by Chi-square test for multiple proportions at each sampling timepoint.
Mortality studies were performed for the same three treatments (YOUNG, OLDER, and S.OLDER). We added additional mock blood meal controls (that is, accounting for any infectious blood meal-associated alteration of mortality) where a mock blood meal was used in place of infectious blood meals. The three controls were: (i) a mock blood meal at 5 dpe (M.Y) to correspond to the YOUNG treatment; (ii) a mock blood meal at 5 dpe, followed by another mock blood meal at 12 dpe (M.M) to correspond to the OLDER two blood meal treatment; and (iii) a mock blood meal at 12 dpe (S.M) to correspond to the S.OLDER, one blood meal treatment. An additional negative control treatment was performed where the mosquitoes were never blood-fed (S). All treatments were cold anesthetized at 5 and 12 dpe, regardless of whether they were offered a blood meal so that all mosquitoes experienced the same treatment, and mosquito density per carton was kept relatively constant with an average of 47 mosquitoes/carton (range 36–58). Each mortality treatment was repeated a total of three times and data are averages of the three replicates.
Cartons were checked daily and, when present, dead mosquitoes were counted and removed up to 28 dpe (approximately 1 month), as this has been shown to be the upper limit of field survival of Ae. aegypti and is similar to the range used in Tesla et al. [13, 41]. Only mosquitoes that took all offered blood meals were included in the mortality analyses.
To test for differences in mortality rates among treatments, Kaplan-Meier survival analyses were conducted and the average time to death (TTD) was estimated. Daily mortality estimates relative to age were then predicted using best fits in R (version 3.5.2) .
Determining age-structured willingness to probe
One of the most influential parameters for determining transmission dynamics of vector-borne pathogens is the biting rate (because it impacts vector-host transmission in both directions). Biting rate is sometimes parameterized as the reciprocal of the number of days between feeds . This assumes that the waiting time between bites is exponentially distributed, and we make this assumption for three biting rates: 0.5 (once every two days); 1 (once a day); and 2 (twice a day) . However, we wanted to determine if the willingness of the mosquito to probe was affected by (i) timing of infectious blood meal and/or (ii) age of the mosquito . While mosquito biting is a function of many factors, it has been shown that heat cues are sufficient to initiate host-seeking behaviors . We use this to determine the willingness or probability of a mosquito probing, which can lead to transmission [46,47,48].
Ten to twelve mosquitoes per each ZIKV-infected treatment (YOUNG, OLDER and S.OLDER) were placed in individual, clear plastic canisters (Bioquip, Rancho Dominquez, CA, USA) 24 h before being provided a blood meal via a membrane feeder using 1 ml discs with 800 μl of blood (Hemotek, Blackburn, UK). This was done at the same dpi schedule as the vector competence studies above. Willingness to bite was assessed using a two-tiered approach by a single observer to control for observation bias. First, mosquitoes were observed through the clear canister for their general position in the canister and secondly, the disc was removed to determine if they were on or near the mesh at the top of the canister. In all cases, these two methods of observation matched. That is, if a mosquito was observed to be at the blood meal prior to disc removal (looking through the canister), it did not move to the bottom of the canister upon disc removal.
This observation was performed at 1, 20 and 45 min post-placement of the disc and the disc was replaced between observation time points. Thus, a mosquito was assessed as “landed” and recorded as “1” if the female was at the top of the canister at any of the observation times. The female was otherwise classified as “not landed” and coded as “0” if she was at the bottom of the canister for all three observation times. We then calculated the probability of biting (Z) as a function of age. Z(age) was determined by fitting the proportion of mosquitoes that landed or fed at least once a day using a self-starting non-linear least squares regression:
where ‘s’ is the asymptote, ‘g’ is the zero-response parameter, and ‘h’ is the logarithmic rate constant. Comparison of probing differences among treatments was assessed both per day post-infection as well as per mosquito age by Kruskal-Wallis test.
Age-structured vectorial capacity and R0
We re-formulated the vectorial capacity equation to estimate VC as a function of the age at which the mosquito acquires an infectious blood meal, redefining the parameters with respect to age at the time of acquisition of infection. We define Ageacquisition as the age at which a mosquito acquires an infection (the day of taking the infectious blood meal) and Agetransmission as the age at which the mosquito subsequently transmits (Fig. 1):
In VCage, ‘m’ is still the mosquito-to-human density and, for illustrative purposes, held constant here at 1. The parameters zacquisition and ztransmission are the probability of a mosquito biting at age of acquisition and age at time of transmission, respectively. The traditional calculation of pN represents the probability of a mosquito living through a fixed number of days N, the EIP, and pN/(− ln(p)) represents the expected infectious days given N and p . In the context of an age-dependent vectorial capacity framework, we can calculate a more precise probability of surviving based on the day the mosquito obtained an infection and the cumulative survival probability to Agetransmission where the cumulative probability of living through the EIP given Ageacquisition is given by:
where ‘pj’ is the probability of daily survival on day j post emergence. We further estimate the infectious period (L) in an age-structured way by numerically deriving L:
where TTD is the average time-to-death derived from the experimental mortality study.
Using the data from our experimental studies, we calculated Agetransmission for two scenarios: EIPmin and EIPmax from our observed transmission data , and calculated VCage as:
VC is a component of R0, which also includes the average infectious period of the human and the probability of transmission from the human to the vector. For illustrative purposes, we calculated R0 using both the traditional calculation of VC (Eq. 1), as well as VCage. We made the following assumptions: mosquito density is held constant at m = 1; the average infectious period of the human is 5 days (r−1) ; and the average probability of transmission from human to vector is parameterized experimentally as the average proportion of mosquitoes that develop a midgut infection given exposure (c). From these data, we can calculate the value of VC needed to push R0 above this threshold. We use this value as a “threshold” to compare VCage across treatment groups and with the traditional VC measures.
All statistics and subsequent graphics were performed and generated using R version 3.5.2 . All packages used are provided in Additional file 1: Text S1. All functions and function parameters used to fit the data and obtain age-dependent distributions of these parameters are given in Additional file 1: Table S2. Goodness-of-fit was assessed either through AIC (for non-linear models) or R2 for linear models.
Within-mosquito dynamics are more affected by time as a function of age versus days post-infection
We first wanted to determine if the age at which a mosquito is offered an infectious blood meal and acquires the ZIKV infection (Ageacquisition) impacted the within-vector kinetics of the mosquito. To do this, we looked at infection, dissemination, and transmission of the three treatments in the context of days post-infection. We found no significant difference in the infection (P > 0.05) or dissemination rates (P > 0.05) across all treatments (Additional file 1: Table S1). However, when direct age comparisons are made, the effect of Ageacquisition on infection becomes obvious (Fig. 2). The average rate of infection of mosquitoes given exposure was 78.2% across all treatments and days post-infection (95% confidence interval: 74.3–82.0%).
When we evaluated the proportion transmitting in terms of days post-exposure (5, 8 and 11 only) none of the treatments were transmitting on days 5 or 8 dpi. On day 11, the transmission rates were 5%, 10% and 15% in groups OLDER, S.OLDER and YOUNG, respectively. These differences were not statistically significant using the Chi-square test for proportions (P > 0.5) (Additional file 1: Table S1).
However, when we investigated mosquitoes from group YOUNG at time points that age-matched OLDER and S.OLDER treatments (Table 1), the YOUNG group achieved a maximum transmission of 45% at 28 days old (23 dpi) versus only 15% for the S.OLDER group and 10% for the OLDER group at 28 days old (16 dpi).
We define the EIPmax as the earliest EIP (days post-infection) where the maximum proportion of mosquitoes are transmitting, and we define the EIPmin as the earliest EIP where the proportion of mosquitoes transmitting is minimal, but greater than 0. EIP values at different vector competence levels have been used to evaluate the vectorial capacity for malaria . For the YOUNG group, EIPmax was 23 dpi at 45%, and EIPmin was 11 dpi at 15%. For the OLDER group, EIPmax was 16 dpi at 10% and EIPmin was 11 dpi at 5%. For the S.OLDER group, EIPmax was 16 dpi at 15% and the EIPmin was 11 dpi at 10%.
Mortality of Aedes aegypti is modestly affected by timing of infectious blood meal
When each infectious treatment was compared to a matching non-infectious control group, the only significant difference in average time to death (TTD) was between the S.OLDER and S.M treatments, with an estimated difference in TTD of two days (Additional file 1: Table S3). Of interest, the non-blood-fed sugar-only controls died significantly faster than any of the blood-fed treatments with an average TTD of 19.6 days (Additional file 1: Table S3), which has been previously shown .
Pairwise comparisons of the ZIKV-exposed treatments determined that group YOUNG had a significantly longer average time to death (TTD) when compared to groups OLDER and S.OLDER, though this difference was modest (0.6 and 1.4 days, respectively) (Fig. 3). The TTD for the YOUNG group was 25.9 days, 25.3 days for the OLDER group, and 24.5 days for the S.OLDER group. This corresponds to average daily survival probabilities of 0.962, 0.961 and 0.960 for the YOUNG, OLDER and S.OLDER groups, respectively. Predicted daily survival rates were generated for the YOUNG group using a non-linear fit and a linear model was fit to the OLDER and S.OLDER groups. The parameters of these models and goodness-of-fit assessments are given in Additional file 1: Table S2, and the observed and predicted values are shown in Additional file 1: Figure S2 for all three treatment groups. Additional comparisons were made among the treatment and control groups, detailed in Additional file 1: Text S2, Figure S3 and Table S3.
Age-dependence of willingness to probe
There was no significant difference in the willingness to probe based on treatment (YOUNG, OLDER and S.OLDER) at each day post-infection via the Kruskal-Wallis test, but there was a significant effect of age (P < 0.05). We then were able to fit a daily probability of probing based on age (Fig. 4).
VCage framework for assessing transmission as a function of age and effect on R0
To calculate VCage, we used treatment-group specific average TTD and predicted daily probabilities of survival (Additional file 1: Figure S4) and the overall daily prediction of willingness to bite in Fig. 4. Since biting rate (the number of bites per day) is a field-derived parameter, we calculated VCage over three biting rates: 0.5 (once every 2 days), 1 (once a day), and 2 (twice a day) . For comparison, we calculated the traditional VC using the average life-dependent traits determined experimentally above, EIPmin and EIPmax, and the biting rates referenced above. These VC values are given in Table 2, along with corresponding values for R0. In the absence of published data regarding the probability of human to mosquito ZIKV transmission success, we used the average infection rate from our experimental data (78.2%) to parameterize human to mosquito transmission, recognizing that this model system used a high dose of ZIKV titer, although such titers in individuals are not unheard of [52, 53].
To achieve an R0 of 1, the vectorial capacity needed to be at or above 0.256. When parameterized with EIPmax, the VCage at a biting rate of 0.5 (once every 2 days) did not achieve this minimal value, but at biting rates of once or twice per day, VCage was sufficient to push R0 over the threshold of one (Fig. 5). We next determined the window of opportunity, which encompasses the days post-emergence when a mosquito acquires an infection and which results in VCage values where R0 would be greater than one. The window of opportunity when we parameterized VCage with EIPmax indicated that a virus must be acquired within the first 2 days post-emergence when the bite rate was twice or once per day (Fig. 5). Interestingly, when parameterized with EIPmin, all bite rates reached sufficient VCage values. The window of opportunity was 14 days post-emergence for a bite rate of twice a day, 12 days post-emergence for a bite rate of once a day, and 7 days post-emergence for a bite rate of once every 2 days (Fig. 5). Even at its highest, VCage indicates that traditional VC calculations are likely overestimates.
OLDER and S.OLDER groups
For both the OLDER and S.OLDER groups, the infectious blood meal occurred at 12 days post-emergence, which was outside the window of opportunity for VCage in all but one scenario (Additional file 1: Figure S4). Only at a bite rate of twice daily did the OLDER group achieve the minimum VCage of 0.256 on or after Ageacquisition of 12 days-old, and only when parameterized with the EIPmin. The window of opportunity for this scenario was days 12–13 post-emergence.
With these vector age-dependent traits, both EIP and vector competence would need to be altered significantly to achieve a VCage ≥ 0.256 at 12 days post-emergence for other scenarios. For example, we consider a hypothetical EIPmax of EIP50 (the time it takes for 50% of exposed mosquitoes to transmit) and a hypothetical EIPmin of EIP10 to demonstrate how these group-specific mortality rates and the age structure of willingness to bite interact to drive vectorial capacity above the threshold in the context of Ageacquisition [13, 49, 54]. The S.OLDER group would need a maximum EIP50 of 12 dpi or 11 dpi for biting rates of twice and once daily, respectively, and 8 dpi for a biting rate of once every 2 days. For the OLDER group, the maximum EIP50 would need to be 13, 12, or 11 dpi for biting rates of twice or once daily, or once every 2 days, respectively. With the EIPmin, similar maximum EIPs were noted of 12 and 9 dpi for biting rates of twice or once daily, respectively. For both the OLDER and S.OLDER groups, with a hypothetical EIPmin of 10%, VCage never got above the threshold needed to get R0 at or above 1.
Prediction of vector-borne disease spread remains difficult, as transmission of vector-borne disease is a dynamic, multifaceted system. This includes life traits of the mosquito, environmental factors, and vector-virus interactions [14, 16, 23, 55, 56]. Here we demonstrate through a combination of experimental and computational methods that mosquito age at the time of pathogen acquisition is a powerful driver of transmission potential due, in large part, to the age-dependence of daily mortality and biting habits. Further, these drivers lead to differences in the estimates of R0.
Our study, which focused upon the age at the time of infectious blood meal, showed no significant impacts on the vector competence of colony Ae. aegypti for ZIKV. While a recent study did show there were significant effects of multiple blood meals on vector competence , we delivered the second blood meal at a much longer interval (7 days versus 3 dpi in ), which likely plays a role in this disparate finding and indicates that there may be age-dependence in this phenomenon as well. However, this hypothesis is outside the scope of the present study. In our mortality study, we did observe a modest difference in TTD. We also observed a very short TTD in mosquitoes only given sugar with no blood meal. We chose not to pursue these effects, as the differences in TTD were modest and the impacts of no blood meal have been previously observed.
The technology for determining the age structure of natural mosquito populations in the field is currently still in development. For example, a study using near-infrared spectroscopy was able to predict the age of female Ae. aegypti ± 2 days, indicating that determining the age structure of a mosquito population is possible, and that such technology could be refined for field studies [58, 59]. Further, mid-infrared spectroscopy had varying, but some promising results in determining the age structure of Anopheles mosquitoes [59, 60]. As these technologies are pursued and refined, there will be a need for ways to understand and quantify age-dependent interactions among vector competence, EIP, mosquito lifespan, and biting behavior [16, 61].
The methodologies in this study highlight the need to understand the quantitative interaction between vector competence, mosquito mortality, and age at time of infection acquisition. Thus, we anticipate that this methodology could be used to explore other important modifiers of vector competence and mosquito mortality and the interaction of the two, such as innate mosquito immunity response to infection, adult and larval microbiome, as well as other extrinsic factors known to affect these traits [2, 3, 5, 10,11,12,13,14,15]. Several recent studies have focused on environmental factors such as temperature, and found that temperature not only affects vector competence of many arboviruses, especially in Aedes aegypti, but also several life traits of the mosquito [41, 55, 56]. This means that transmission is ultimately impacted by the interactions of all of these temporally dependent processes, and future studies should consider the age structure when assessing these impacts. Studies have also shown the impact of vector species, in particular the difference between field-derived and laboratory strains [2, 3]. Here, we used a laboratory strain of Ae. aegypti, which has been shown to vary in its competence when compared to a field-derived strain. Future use of VCage could highlight these differences, as well as differences among various arbovirus-arthropod systems.
Our results demonstrate the importance of age structure when evaluating the fitness of a mosquito-virus system and indicate that R0 may be overestimated when it is not considered. This framework further shows that how within-mosquito arbovirus fitness is measured, often by comparing proportions of transmitting mosquitoes at arbitrary time points, is not sufficient. Here, when the Ageacquisition was advanced, the difference in hypothetical EIPs necessary for the system to succeed was not very different. For example, in the OLDER and S.OLDER groups, there was only one day difference between the EIPmin and EIPmax needed for success at biting rates of 2 and 1, though these quantities represented a difference of 40% in the proportion of mosquitoes transmitting. The same trend was demonstrated in the YOUNG group where EIPmin resulted in longer windows of opportunity compared to EIPmax. Additionally, when the minimum time to transmission is shorter, a larger portion of the mosquito population contributes to transmission, because VCage implicates older mosquitoes in this scenario given no significant reduction in mortality due to extrinsic factors . Thus, VCage reveals that the temporality of the within-mosquito arbovirus dynamics is more impactful than the ultimate magnitude of this widely used fitness measure. More investigations into the earlier and minimal transmission rates, such as at earlier time points that we did not consider, may be warranted to quantify fully the contribution of lower but faster dissemination profiles in arbovirus systems.
A study by Althouse et al.  also found that the temporality of transmission from non-human primates was sometimes more impactful than the magnitude of the viremia leading to transmission to the mosquito. They proposed a “tortoise-and-the-hare” (TatH) model to describe this relationship between arboviral viremia profiles in non-human primates and the predicted transmission success to vectors, showing that the strategy of “slow and steady” viremia, i.e. lower levels for longer periods, resulted in higher predicted transmission success of arboviruses . This same TatH model recently described macro-transmission dynamics in Colombia, where it was demonstrated that slow burn-in epidemics actually resulted in cumulatively more cases and higher R0 values than in initially explosive outbreaks . Vector competence is a continuous process over time, and as such, EIPmin and EIPmax are not independent measures. Thus, vector competence profiles with higher transmission rates at earlier times (“hare strategy”) will always be more fit than lower proportions at the same EIP. However, this study demonstrates that in systems where the maximum measured vector competence is low, but the time to minimum transmission is short (“tortoise strategy”) , there is still a good chance the system will succeed. VCage also suggests that how arbovirus phenotypes are compared and ranked, and perhaps even how the field identifies ‘highly’ or ‘negligibly’ competent vectors, may need adjustment in the context of mosquito age.
We demonstrate that mosquito age may not affect experimentally determined infection, dissemination, or transmission rates when assessed in the traditional days post-infection manner. However, when we consider mosquito mortality, EIP, and vector competence in the context of a more holistic measure of transmission potential, vectorial capacity, and further adjust VC for the effect of age at the time of mosquito infection, we demonstrate that age can be an important factor, and that how some arbovirus fitness measures are assessed may need reconsideration.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional file.
extrinsic incubation period
Costa-da-Silva AL, Ioshino RS, Araujo HR, Kojin BB, Zanotto PM, Oliveira DB, et al. Laboratory strains of Aedes aegypti are competent to Brazilian Zika virus. PLoS ONE. 2017;12:e0171951.
Roundy CM, Azar SR, Rossi SL, Huang JH, Leal G, Yun R, et al. Variation in Aedes aegypti mosquito competence for Zika virus transmission. Emerg Infect Dis. 2017;23:625–32.
Chouin-Carneiro T, Vega-Rua A, Vazeille M, Yebakima A, Girod R, Goindin D, et al. Differential susceptibilities of Aedes aegypti and Aedes albopictus from the Americas to Zika virus. PLoS Negl Trop Dis. 2016;10:e0004543.
Musso D, Gubler DJ. Zika virus. Clin Microbiol Rev. 2016;29:487–524.
Valderrama A, Diaz Y, Lopez-Verges S. Interaction of Flavivirus with their mosquito vectors and their impact on the human health in the Americas. Biochem Biophys Res Commun. 2017;492:541–7.
Kramer LD, Ciota AT. Dissecting vectorial capacity for mosquito-borne viruses. Curr Opin Virol. 2015;15:112–8.
Christofferson RC, Mores CN, Wearing HJ. Bridging the gap between experimental data and model parameterization for chikungunya virus transmission predictions. J Infect Dis. 2016;214(Suppl. 5):S466–70.
Ohm JR, Baldini F, Barreaux P, Lefevre T, Lynch PA, Suh E, et al. Rethinking the extrinsic incubation period of malaria parasites. Parasit Vectors. 2018;11:178.
Bellan SE. The importance of age dependent mortality and the extrinsic incubation period in models of mosquito-borne disease transmission and control. PLoS ONE. 2010;5:e10165.
Maciel-de-Freitas R, Koella JC, Lourenco-de-Oliveira R. Lower survival rate, longevity and fecundity of Aedes aegypti (Diptera: Culicidae) females orally challenged with dengue virus serotype 2. Trans R Soc Trop Med Hyg. 2011;105:452–8.
Brady OJ, Godfray HC, Tatem AJ, Gething PW, Cohen JM, McKenzie FE, et al. Vectorial capacity and vector control: reconsidering sensitivity to parameters for malaria elimination. Trans R Soc Trop Med Hyg. 2016;110:107–17.
Robert MA, Christofferson RC, Silva NJ, Vasquez C, Mores CN, Wearing HJ. Modeling mosquito-borne disease spread in U.S. urbanized areas: the case of dengue in Miami. PLoS ONE. 2016;11:e0161365.
Goindin D, Delannay C, Ramdini C, Gustave J, Fouque F. Parity and longevity of Aedes aegypti according to temperatures in controlled conditions and consequences on dengue transmission risks. PLoS ONE. 2015;10:e0135489.
Hardy JL, Houk EJ, Kramer LD, Reeves WC. Intrinsic factors affecting vector competence of mosquitoes for arboviruses. Annu Rev Entomol. 1983;28:229–62.
Souza-Neto JA, Powell JR, Bonizzoni M. Aedes aegypti vector competence studies: a review. Infect Genet Evol. 2019;67:191–209.
Lounibos LP, Kramer LD. Invasiveness of Aedes aegypti and Aedes albopictus and vectorial capacity for chikungunya virus. J Infect Dis. 2016;214(Suppl. 5):S453–8.
Macdonald G. The epidemiology and control of malaria. London: Oxford University Press; 1957.
Garrett-Jones C. Prognosis for interruption of malaria transmission through assessment of the mosquito’s cectorial capacity. Nature. 1964;204:1173–5.
Bagny L, Delatte H, Quilici S, Fontenille D. Progressive decrease in Aedes aegypti distribution in Reunion Island since the 1900s. J Med Entomol. 2009;46:1541–5.
Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog. 2007;3:e201.
Vazeille M, Moutailler S, Coudrier D, Rousseaux C, Khun H, Huerre M, et al. Two chikungunya isolates from the outbreak of La Reunion (Indian Ocean) exhibit different patterns of infection in the mosquito, Aedes albopictus. PLoS ONE. 2007;2:e1168.
Kilpatrick AM, Meola MA, Moudy RM, Kramer LD. Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mosquitoes. PLoS Pathog. 2008;4:e1000092.
Ciota AT, Kramer LD. Vector-virus interactions and transmission dynamics of West Nile virus. Viruses. 2013;5:3021–47.
Moudy RM, Meola MA, Morin LL, Ebel GD, Kramer LD. A newly emergent genotype of West Nile virus is transmitted earlier and more efficiently by Culex mosquitoes. Am J Trop Med Hyg. 2007;77:365–70.
Dye C. The analysis of parasite transmission by bloodsucking insects. Annu Rev Entomol. 1992;37:1–19.
Ruckert C, Ebel GD. How do virus-mosquito interactions lead to viral emergence? Trends Parasitol. 2018;34(4):310–21.
Petersen MT, Silveira IDD, Tatila-Ferreira A, David MR, Chouin-Carneiro T, Van den Wouwer L, et al. The impact of the age of first blood meal and Zika virus infection on Aedes aegypti egg production and longevity. PLoS One. 2018;13:e0200766.
Christofferson RC, Mores CN. Estimating the magnitude and direction of altered arbovirus transmission due to viral phenotype. PLoS ONE. 2011;6:e16298.
Smith DL, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004;3:13.
Massad E, Coutinho FA. Vectorial capacity, basic reproduction number, force of infection and all that: formal notation to complete and adjust their classical concepts and equations. Mem Inst Oswaldo Cruz. 2012;107:564–7.
Dye C. Vectorial capacity: must we measure all its components? Parasitol Today. 1986;2:203–9.
Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS Pathog. 2012;8:e1002588.
Garrett-Jones C, Shidrawi GR. Malaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomology. Bull World Health Organ. 1969;40:531–45.
Styer LM, Carey JR, Wang JL, Scott TW. Mosquitoes do senesce: departure from the paradigm of constant mortality. Am J Trop Med Hyg. 2007;76:111–7.
Novoseltsev VN, Michalski AI, Novoseltseva JA, Yashin AI, Carey JR, Ellis AM. An age-structured extension to the vectorial capacity model. PLoS ONE. 2012;7:e39479.
Rock KS, Wood DA, Keeling MJ. Age- and bite-structured models for vector-borne diseases. Epidemics. 2015;12:20–9.
Kawiecki AB, Mayton EH, Dutuze MF, Goupil BA, Langohr IM, Del Piero F, et al. Tissue tropisms, infection kinetics, histologic lesions, and antibody response of the MR766 strain of Zika virus in a murine model. Virol J. 2017;14:82.
Ciota AT, Bialosuknia SM, Zink SD, Brecher M, Ehrbar DJ, Morrissette MN, et al. Effects of Zika virus strain and Aedes mosquito species on vector competence. Emerg Infect Dis. 2017;23:1110–7.
Richards SL, Pesko K, Alto BW, Mores CN. Reduced infection in mosquitoes exposed to blood meals containing previously frozen flaviviruses. Virus Res. 2007;129:224–7.
Tramonte AR, Christofferson RC. Investigating the probability of establishment of Zika virus and detection through mosquito surveillance under different temperature conditions. PLoS ONE. 2019;14:e0214306.
Tesla B, Demakovsky LR, Mordecai EA, Ryan SJ, Bonds MH, Ngonghala CN, et al. Temperature drives Zika virus transmission: evidence from empirical and mathematical models. Proc Biol Sci. 2018;285:20180795.
R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. https://www.r-project.org/.
Paaijmans KP, Cator LJ, Thomas MB. Temperature-dependent pre-bloodmeal period and temperature-driven asynchrony between parasite development and mosquito biting rate reduce malaria transmission intensity. PLoS ONE. 2013;8:e55777.
Christofferson RC, Mores CN, Wearing HJ. Characterizing the likelihood of dengue emergence and detection in naive populations. Parasit Vectors. 2014;7:282.
Zermoglio PF, Robuchon E, Leonardi MS, Chandre F, Lazzari CR. What does heat tell a mosquito? Characterization of the orientation behaviour of Aedes aegypti towards heat sources. J Insect Physiol. 2017;100:9–14.
Styer LM, Kent KA, Albright RG, Bennett CJ, Kramer LD, Bernard KA. Mosquitoes inoculate high doses of West Nile virus as they probe and feed on live hosts. PLoS Pathog. 2007;3:1262–70.
Matsuoka H, Yoshida S, Hirai M, Ishii A. A rodent malaria, Plasmodium berghei, is experimentally transmitted to mice by merely probing of infective mosquito, Anopheles stephensi. Parasitol Int. 2002;51:17–23.
Dubrulle M, Mousson L, Moutailler S, Vazeille M, Failloux AB. Chikungunya virus and Aedes mosquitoes: saliva is infectious as soon as two days after oral infection. PLoS ONE. 2009;4:e5895.
Paaijmans KP, Blanford S, Chan BH, Thomas MB. Warmer temperatures reduce the vectorial capacity of malaria mosquitoes. Biol Lett. 2012;8:465–8.
Kucharski AJ, Funk S, Eggo RM, Mallet HP, Edmunds WJ, Nilles EJ. Transmission dynamics of Zika virus in island populations: a modelling analysis of the 2013–14 French Polynesia outbreak. PLoS Negl Trop Dis. 2016;10:e0004726.
Styer LM, Minnick SL, Sun AK, Scott TW. Mortality and reproductive dynamics of Aedes aegypti (Diptera: Culicidae) fed human blood. Vector Borne Zoonotic Dis. 2007;7:86–98.
Perkasa A, Yudhaputri F, Haryanto S, Hayati RF, Ma’roef CN, Antonjaya U, et al. Isolation of Zika virus from febrile patient, Indonesia. Emerg Infect Dis. 2016;22:924–5.
Waggoner JJ, Gresh L, Vargas MJ, Ballesteros G, Tellez Y, Soda KJ, et al. Viremia and clinical presentation in Nicaraguan patients infected with Zika virus, chikungunya virus, and dengue virus. Clin Infect Dis. 2016;63:1584–90.
Carrington LB, Armijos MV, Lambrechts L, Scott TW. Fluctuations at a low mean temperature accelerate dengue virus transmission by Aedes aegypti. PLoS Negl Trop Dis. 2013;7:e2190.
Alto BW, Wiggins K, Eastmond B, Ortiz S, Zirbel K, Lounibos LP. Diurnal temperature range and chikungunya virus infection in invasive mosquito vectors. J Med Entomol. 2018;55:217–24.
Muttis E, Balsalobre A, Chuchuy A, Mangudo C, Ciota AT, Kramer LD, et al. Factors related to Aedes aegypti (Diptera: Culicidae) populations and temperature determine differences on life-history traits with regional implications in disease transmission. J Med Entomol. 2018;55:1105–12.
Armstrong PM, Ehrlich HY, Magalhaes T, Miller MR, Conway PJ, Bransfield A, et al. Successive blood meals enhance virus dissemination within mosquitoes and increase transmission potential. Nat Microbiol. 2020;5:239–47.
Sikulu-Lord MT, Milali MP, Henry M, Wirtz RA, Hugo LE, Dowell FE, et al. Near-infrared spectroscopy, a rapid method for predicting the age of male and female wild-type and Wolbachia infected Aedes aegypti. PLoS Negl Trop Dis. 2016;10:e0005040.
Lambert B, Sikulu-Lord MT, Mayagaya VS, Devine G, Dowell F, Churcher TS. Monitoring the age of mosquito populations using near-infrared spectroscopy. Sci Rep. 2018;8:5274.
Gonzalez Jimenez M, Babayan SA, Khazaeli P, Doyle M, Walton F, Reedy E, et al. Prediction of mosquito species and population age structure using mid-infrared spectroscopy and supervised machine learning. Wellcome Open Res. 2019;4:76.
Johnson BJ, Hugo LE, Churcher TS, Ong OTW, Devine GJ. Mosquito age grading and vector-control programmes. Trends Parasitol. 2020;36:39–51.
Onyango MG, Bialosuknia SM, Payne AF, Mathias N, Kuo L, Vigneron A, et al. Increased temperatures reduce the vectorial capacity of Aedes mosquitoes for Zika virus. Emerg Microbes Infect. 2020;9:67–77.
Althouse BM, Hanley KA. The tortoise or the hare? Impacts of within-host dynamics on transmission success of arthropod-borne viruses. Philos Trans R Soc Lond B Biol Sci. 2015. https://doi.org/10.1098/rstb.2014.0299.
Pena-Garcia VH, Christofferson RC. Correlation of the basic reproduction number (R0) and eco-environmental variables in Colombian municipalities with chikungunya outbreaks during 2014–2016. PLoS Negl Trop Dis. 2019;13:e0007878.
Goddard LB, Roth AE, Reisen WK, Scott TW. Vector competence of California mosquitoes for West Nile virus. Emerg Infect Dis. 2002;8:1385–91.
We would like to thank Dr. Barbara Johnson for providing our viral strain, as well as Dr. Daniel Swale for providing us with mosquito eggs. We would like to thank all our laboratory colleagues for putting up with the hogging of freezer space to hold all of the samples generated for this study. And last but not least, we’d like to thank Aesop for his apparently ubiquitously applicable fable.
This study was funded through NIH/NIGMS R01GM122077.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Illustration of main treatment design for vector competence experimentation. Table S1. Infection and dissemination rates for each day post-infection (dpi) and corresponding mosquito age for each of the three ZIKV treatments. Table S2. Modeled fits of parameters, the type of model, and the parameter values. Figure S2. Observed and predicted daily probabilities of survival for the three treatment groups. Text S1. R packages used in the study. Text S2. Comparison of mortality of Aedes aegypti with respect to bloodmeals. Figure S3. Mortality curves for all treatments and associated controls. Table S3. Mean time to death and sample size for ZIKV-infection treatments and unexposed controls used in the mortality study. Figure S4. Demonstration of the VCage for the two treatments exposed at older ages.
About this article
Cite this article
Mayton, E.H., Tramonte, A.R., Wearing, H.J. et al. Age-structured vectorial capacity reveals timing, not magnitude of within-mosquito dynamics is critical for arbovirus fitness assessment. Parasites Vectors 13, 310 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s13071-020-04181-4
- Extrinsic incubation period
- Vector competence
- Aedes aegypti
- Vectorial capacity
- Biting rate