|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Since the 1950s, the assumption that mosquitoes do not senesce (i.e., that mosquito mortality remains constant with age) has been used to assess the role of mosquitoes in pathogen transmission and predict public health consequences of vector control strategies. This assumption was first articulated by MacDonald,3 who reasoned that environmental insults, disease, and predation would kill mosquitoes before they had an opportunity to die of old age. Although evidence from a few studies4,5 seemed to support this assumption, many others suggest that senescence occurs in both captive and free-ranging mosquito populations.610 For purposes of this discussion, we define senescence as an increased probability of mortality with advancing age. In one notable study, published reports of mosquito mortality were reanalyzed, showing that many species exhibit age-dependent mortality consistent with a Gompertz mortality function, and that substantial error was introduced into pathogen transmission models when age-dependent mortality was not applied.2 In addition to increased mortality rates, adverse age-related changes in mosquito flight,11 immune function,12 salivary gland structure,13 and detoxification mechanisms14,15 have been reported, suggesting that mosquitoes are no different from most other organisms in experiencing deleterious functional and structural changes with age.
Over a decade ago, Gillies16 advocated that the assumption of age-independent mortality be "buried" because resulting estimates were, "at best, approximations." Despite this admonition and growing evidence that mosquitoes senesce, non-senescence has continued to be assumed in even the most recent disease transmission models.1720 This reluctance to change can be attributed to two factors. First, there are no published reports that provide the detailed mortality trajectory necessary to properly test the assumption. Laboratory studies have generally used populations too small (N < 1,000) to detect differences, and data from field mortality studies are imprecise. Large populations are essential to characterize mortality early in life when few mosquitoes are dying and later in life when few mosquitoes remain alive.21 Second, the assumption of age-independent mortality greatly simplifies the mathematics of models and reduces the number and complexity of variables that need to be considered.22
The major limitation of assuming non-senescence is that it leads to the unrealistic, simplified view that the vector potential of all mosquitoes is the same, regardless of their age. This age-independent assumption applies not only to mortality but also to other determinants of pathogen transmission; for example, biting rate, host preference, vector competence, dispersal, resistance to insecticides, and pathogen-induced adverse effects.22 The pervasiveness of the assumption is shown by the fact that, in the last 50 years, only a few publications have considered the effect of mosquito age on transmission factors.11,12,15,2327 Moreover, because vector mortality is a particularly sensitive component of pathogen transmission, quantitative models that assume non-senescence produce results with substantial errors.2
To rigorously test the hypothesis that mosquitoes do senesce, we carried out two mortality studies, including the largest mosquito mortality study to date (N > 100,000) and a mortality study of blood-fed and reproductively active mosquitoes. To determine the epidemiologic implications of our results, we modified the vectorial capacity28 model to include age-dependent survival and life expectancy values. Our study subject was Aedes aegypti, the most important vector worldwide of dengue virus and an urban vector of yellow fever.29
| MATERIALS AND METHODS |
|
|
|---|
Fitting data to mortality models. We used the software WinModest 1.033 to determine the mortality model that best fit observed mortality patterns and to estimate model parameters. Standard mortality models are based on age-specific hazard, ux, an instantaneous measure of mortality that can be estimated from empirical data with the approximation
![]() |
where px is the probability that an individual alive at age x survives to age x + 1.32 Briefly, the software fit a series of semi-hierarchical models to the data and calculated the log-likelihood for each model (Table 1
). A large log-likelihood indicates a better fit. A test statistic of twice the difference in the log-likelihoods of the nested models was used to determine if the increase was significant.
|
Mortality study of blood-fed mosquitoes. We conducted a second mortality study to determine the effect of blood feeding and reproduction on mortality. The F2F4 generations of an Ae. aegypti colony obtained from pupae collected in Mae Sot Province, Thailand, were used. Immature mosquitoes were raised as described above. Pupae were collected, and adults were allowed to emerge into a 30 x 30 x 30-cm screened cage for 24 hours. Approximately 100 2-day-old adult males and 100 2-day-old females were placed into each of six 15 x 60 x 90-cm screened cages containing an oviposition cup and maintained under semi-controlled climatic conditions (2530°C; ~1°C daily fluctuations, 5°C drift in mean daily temperature; 2550% RH, 12:12 L:D). Oviposition cups consisted of a 60-mm petri dish containing a water-soaked sponge covered with two layers of filter paper. Each cage was assigned to one of three treatments (two cages per treatment): 10% sucrose, human blood, or human blood plus 10% sucrose. Mosquitoes in the blood treatments were given an opportunity to feed on a human forearm introduced into the appropriate cages for 15 minutes every day. The use of human subjects for feeding mosquitoes was approved and conducted in accordance with the University of California at Davis Human Subjects Protocol 991781. Water and 10% sucrose solution (if provided) were continuously available from cotton wicks introduced through ports in each cage. Each cage was inspected daily, dead mosquitoes were removed, censored deaths were recorded, oviposition cups were collected, and life tables were constructed. Oviposition data and its relation to mortality are discussed in another publication.35 We used the software WinModest 1.033 to determine the mortality model that best fit observed mortality patterns and to estimate model parameters.
Age-dependent vectorial capacity model.
To determine the epidemiologic implications of mosquito senescence for pathogen transmission, we derived an age-dependent vectorial capacity model and used it to compare the vectorial capacity28 (predicted number of new infectious mosquito bites that result from one infectious host entering a naïve population) of simulated mosquito populations exhibiting age-independent versus age-dependent mortality with two age structures (stable age distribution and synchronous adult emergence). Specifically, the model considers the effect of mosquito age when first imbibing an infectious blood meal on subsequent vectorial capacity of that age cohort. The model defines three vectorial capacity concepts: 1) Cx = age-specific vectorial capacity: the daily number of potentially infective bites resulting from a mosquito of age x biting one infectious host; 2) Cx x
x = age class-specific vectorial capacity: the daily number of potentially infective bites resulting from all mosquitoes of age class x biting one infectious host; and 3) Ct = total population vectorial capacity: the sum of age class-specific vectorial capacities over all biting age classes. Cx (age-specific vectorial capacity) is similar to the concept of individual vectorial capacity (IC) defined by Saul and others36 as a function of age.
The modified model is described by the equation,
![]() |
where
= age at which mosquitoes begin biting hosts;
= oldest biting age class of mosquitoes in the population;
x = fraction of the population in age class x; m = mosquito density; a = mosquito daily biting rate; n = extrinsic incubation period;
i=xx+n pi = product of mosquito daily survival rates from age x when first biting an infectious host to age x + n, which is the end of the extrinsic incubation period; and ex+n = expectation of remaining infectious life at age x + n. Model simulations were performed in Excel (Microsoft, Redmond, WA). The following values were used in the model and assumed to remain constant:
= 3 days,29 m = 1.5 mosquitoes/person,37 a = 0.75 bites/d,38 and N = 10 days.37 The mean female mortality pattern from the large-scale study plus the best fitting logistic (a = 0.0018, b = 0.1416, s = 1.0730), Gompertz (a = 0.0066, b = 0.0623), and exponential (ux = a, where a = 0.0313) mortality models were used in the modified vectorial capacity model; all mortality patterns had a mean lifespan of 32 days. We used two population age structures in our simulations: stable age distribution (derived from the Lotka equation32) and synchronous emergence (total population was distributed equally among age classes 3, 4, 5, and 6 days). We calculated the stable age distribution separately for each mortality model, using the same birth rate (estimated from a laboratory reproductive study) and the appropriate mortality schedule.
We examined sensitivity of the vectorial capacity model to changes in overall mortality levels by adjusting parameters of the exponential, Gompertz, and logistic mortality patterns to achieve mean life spans of 5, 10, 16, 20, and 25 days and using this adjusted mortality data in the model.
| RESULTS |
|
|
|---|
Although mortality patterns differed among cohorts, all cohorts exhibited age-dependent mortality. In general, mortality was low for young mosquitoes (< 10 days old), increased with age, and decelerated or declined at the oldest ages (Figure 1
). Young mosquitoes are defined by very low mortality, which, in this laboratory study, was observed in mosquitoes < 10 days old.
|
Fitting data to mortality models.
Logistic or logistic-Makeham models provided the best fit for mortality patterns of all female cohorts and 28 of 29 male cohorts (1 was best fit by the Gompertz model). Best-fit models for the combined data of all cohorts were similar: logistic-Makeham for males and logistic for females. Mortality model parameter values varied significantly between the sexes. Mean rate of exponential mortality increase (b) for female cohorts was less than one half that of male cohorts (0.163 ± 0.017 versus 0.357 ± 0.071, respectively; Wilcoxon signed rank = 206.5; P < 0.0001; Figure 2B
). Similarly, mean rate of mortality deceleration (s) was significantly lower for females than males (1.417 ± 0.234 versus 2.001 ± 0.392, respectively; Wilcoxon signed rank = 123.5; P = 0.0053; Figure 2C
). Other parameter values (a and c) did not vary significantly by sex (a: Wilcoxon signed rank = 16, P = 0.74; c: Wilcoxon signed rank = 25, P = 0.37; Figure 2A and D
).
|
|
|
|
Our sensitivity analysis indicated that as overall mortality increases (average lifespan [E0] decreases), the potential for mosquito populations to transmit pathogens declines (Table 2
). At all mortality levels, synchronous emergence populations have higher total vectorial capacity (Ct) values than stable age distribution populations because of the greater numbers of young mosquitoes. When lifespan is held constant (16 days), predicted vectorial capacity (Ct) can vary 2-fold depending on specific mortality parameters used. In general, total vectorial capacity is more sensitive to changes in b (exponential mortality increase with age) than a (initial mortality rate), which is understandable because b is exponentially, whereas a is multiplicatively related to overall mortality (Table 1
).
|
| DISCUSSION |
|
|
|---|
Importantly, the companion paper to this one39 shows senescence in free-ranging Ae. aegypti in Thailand. That study used a novel synthetic cohort technique in which mosquitoes of different ages were released simultaneously into homes in a Thai village and subjected to recapture over a short period of time (2 days). Recapture rates decreased with mosquito age, which was consistent with the hypothesis that mortality increases with age. Age-dependent recapture rates could not be explained by differences in mosquito dispersal by age.
Logistic or logistic-Makeham models provided the best fit to female mortality data from both laboratory mortality studies. This result is similar to other large-scale studies on a variety of organisms including yeast, nematodes, flies, and humans, in which mortality deceleration was seen at the oldest ages.40 Previous studies reported that female mortality of several mosquito species follows a Gompertz function.2,6 Those earlier studies, however, generally used small sample sizes (< 1,000) and did not fit results to mortality models other than the Gompertz model.
Our results indicate that the age at which a mosquito bites an infectious host is critical for determining the probability that it will transmit a pathogen. The importance of old mosquitoes has long been recognized because they are the only ones that could have survived long enough to become infectious (i.e., survived to an epidemiologically dangerous age). Our simulations show that young mosquitoes are more important than may have been generally appreciated. Because of increasing risk of mortality with advancing age, a mosquito that first imbibes pathogen-infected blood when it is young has a greater potential to transmit pathogens to more hosts during its lifetime than one that first takes an infected meal when it is old. Younger mosquitoes have lower mortality and thus are more likely to survive through the pathogen incubation period to become infectious than old conspecifics.
Because of the impact of mosquito age on pathogen transmission rates, understanding mosquito population age structure is important. In our simulations, mosquito populations of the same size but with different age structures were predicted to produce significantly different levels of pathogen transmission. For example, when a pathogen, such as West Nile virus or Plasmodium falciparum, first enters an area because of migration of infectious vertebrate hosts or initiation of seasonal transmission, a young mosquito population is epidemiologically most "dangerous" because it consists entirely of individuals with low mortality rates. They have the highest probability of a few surviving long enough to become infectious after biting an infectious host. This type of age structure can be created if large numbers of adult mosquitoes emerge synchronously, an event that can be triggered by seasonally heavy rainfall or temperature increases. Several mosquito species including Anopheles dirus,41 a vector of malaria in Thailand, and Ochlerotatus provocans,42 a vector of James-town Canyon virus, have been shown to exhibit large synchronous emergences caused by heavy seasonal rainfall and the seasonal appearance of temporary aquatic development sites. All else being equal, vectorial capacity of populations composed of synchronously emerged adults will be high initially, decrease as the number of young mosquitoes declines, and increase again as the next generation emerges. Populations with a stable age distributionwhich is our general expectation for the subject of our study, Ae. aegypti31,43will tend, over time, to have more constant vectorial capacity because of continuous adult emergence. Other factors, such as temperature, which can affect pathogen transmission by expanding or contracting extrinsic incubation, may have more of an impact on intra-annual or seasonal fluctuations in dengue transmission than age structure of Ae. aegypti populations. Without perturbations in Ae. aegypti density caused by vector control, over longer periods of time, inter-annual or between-year fluctuations in endemic dengue transmission are most likely associated with shifts in herd immunity.44
Our results have three important implications. First, departure from the paradigm of constant mortality frees us to view mosquitoes and the diseases they transmit as complex, dynamic systems, not ones that can be easily described by static age-independent models. At the organismal level, this opens up a plethora of largely unexplored topics like those regarding physiologic processes that cause senescence, the effect of senescence on mosquito behavior, and the large mortality differential between the mosquito sexes. At a larger scale, mosquito senescence leads to unanswered questions concerning the effect of variable death rates on life history and evolution of insect vectors and the pathogens they transmit. At a molecular level, previously unconsidered age-dependent changes in the genetic underpinnings of a mosquitos role in pathogen transmission (i.e., gene expression) merit study.
Second, within the paradigm of age-dependent mortality, existing methods to study mosquito populations are no longer adequate. New innovative methods must be developed to better understand mosquito population dynamics. Currently the most widely used age-grading technique for Ae. aegypti can only separate mosquitoes into two age classes (parous, > 5 days old; non-parous, < 5 days old).45 New age-grading techniques with greater discriminatory power need to be developed, as well as novel means to study aging in wild mosquitoes. Techniques such as cuticular hydrocarbon age-grading46,47 and the synthetic cohort method39,48 are examples of such new methods.
Third, our results indicate that age-dependent factors should be included in vector-borne disease transmission models to accurately describe and predict the dynamics of pathogen transmission. Many vector-borne diseases have no clinical cure or vaccine; vector control is the primary option available to contain and prevent pathogen transmission. A more complete and accurate description of pathogen transmission will improve the efficiency by which control measures are developed, administered, and evaluated.
We expect that the phenomenon of mosquito age-dependent mortality is not limited to Ae. aegypti. Because of the important effects of mortality on pathogen transmission, patterns of mortality for virtually all arthropod vectors merit review within the context of their significance to disease in humans and other animals.
Received May 31, 2006. Accepted for publication September 6, 2006.
Acknowledgments: The authors thank Miriam Khaliqi, Beverly King, Sharon Minnick, Susie Turold, and Anthony Vo for assistance on this project.
Financial support: This work was supported by grants from the University of California Mosquito Research Program to T.W.S. and National Institutes of Health (AI22119 to T.W.S. and AG022500-01 and AG08761-10 to J.R.C.).
* Address correspondence to Thomas W. Scott, Department of Entomology, One Shields Avenue, Davis, CA 95616. E-mail: twscott{at}ucdavis.edu ![]()
Authors addresses: Linda M. Styer, Arbovirus Laboratories, Wadsworth Center, New York State Department of Health, 5668 State Farm Road, Slingerlands, NY 12159, Telephone: 518-862-4306, E-mail: lstyer{at}wadsworth.org. James R. Carey, Department of Entomology, One Shields Avenue, Davis, CA 95616, Telephone: 530-752-6217, E-mail: jrcarey{at}ucdavis.edu. Jane-Ling Wang, Department of Statistics, One Shields Avenue, Davis, CA 95616, Telephone: 530-752-7624, E-mail: janelwang{at}ucdavis.edu. Thomas W. Scott, Department of Entomology, One Shields Avenue, Davis, CA 95616, Telephone: 530-754-4196, Fax: 530-752-1537, E-mail: twscott{at}ucdavis.edu.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
P.A Hancock, M.B Thomas, and H.C.J Godfray An age-structured model to evaluate the potential of novel malaria-control interventions: a case study of fungal biopesticide sprays Proc R Soc B, January 7, 2009; 276(1654): 71 - 80. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |