| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
 |
ABSTRACT |
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
|
INTRODUCTION |
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
Traditional Ae. aegypti control measures include elimination (source reduction) or treatment of larval habitats to prevent production of adults and insecticidal space spraying to reduce adult population densities.5,6 Contemporary programs emphasize reducing Ae. aegypti populations to levels that prevent or slow virus transmission with the ultimate objective of decreasing the incidence of disease, especially severe, life-threatening illness. However, traditional entomologic surveillance techniques are based on a series of indices that were designed to detect the presence or absence of Ae. aegypti larvae. Those methods assume a strong positive correlation between the presence of larvae and adult females in a household: only adult females transmit virus to humans. There are, however, three important reasons to question the strength of the larvae-adult association. First, because larval mortality can be high, adults may not emerge from a container holding immature mosquitoes. Alternative entomologic surveillance methods, especially pupal surveys, were developed to circumvent this shortcoming.7 Second, because adults are capable of flight, they can move away and become spatially disassociated from their development sites. Third, independent of the surveillance technique (larvae, pupae, or adult collections) citywide surveys are often carried out in such a way that the number and location of households selected are derived from standard parametric sample size calculations. The assumption that there is no spatial structure among infested houses must be validated.
The purpose of this study was to characterize the spatial distribution of Ae. aegypti populations in two representative neighborhoods in the Amazonian city of Iquitos, Peru over two time periods. Specifically, from complete samples of households in two areas of Iquitos we examined the 1) underlying spatial structure of Ae. aegypti infestations (larvae, pupae, and adult), 2) temporal stability of that structure, and 3) correlation between clusters at different life stages of the mosquito. We conclude by discussing the implications of our findings on estimation of entomologic risk to epidemiologic studies of dengue and routine dengue surveillance.
|
MATERIALS AND METHODS |
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
). The major industries in Iquitos are small commercial enterprises, fishing, oil, lumber, and to some extent agriculture.
|
). We selected these two neighborhoods because they were characterized as areas of high (Maynas) and low (Tupac Amaru) prevalence of human dengue infection in previous informal studies (Morrison AC, Scott TW, unpublished data). Although Maynas could be characterized as the wealthier and older of the two neighborhoods, households within both areas vary greatly in socioeconomic status so that well constructed households with piped water and poorly constructed households with no water or sewer services exist in both neighborhoods in a patchwork. Nevertheless, there are some distinct differences between the two neighborhoods. Maynas has a higher proportion than Tupac Amaru of permanent houses constructed with bricks and concrete. Conversely, Tupac Amaru is a community in transition from predominantly temporary wood houses with palm roofs to houses constructed with brick and concrete. Even though Maynas has a better-developed sewer system than Tupac Amaru, the Maynas water supply is inconsistent. Consequently, Maynas residents are more likely than those in Tupac Amaru to store water in containers that are potential development sites for immature Ae. aegypti. In contrast, Tupac Amaru has many open sewers but because of close proximity to the city water plant most houses have a stable water supply and are less likely to store water than in Maynas. Study design. A unique-house code was painted on the front of each of the 550 houses located on 20 blocks in Maynas and the 510 houses located on 14 blocks in Tupac Amaru. Almost all houses have at least one wall in common with a neighboring house. Beginning in mid-November 1998, five two-person entomology collection teams were provided a map of a block to be surveyed with a designated start house. Households were surveyed in sequence daily along the block from the start house between 7:00 AM and 1:00 PM. Unoccupied or closed houses and houses where residents did not provide permission for the survey, businesses, offices, and schools were not sampled. Thus, we were able to survey 95% of the houses in both surveys: 528 in Maynas and 481 in Tupac Amaru. Collecting teams were rotated among blocks each day in an attempt to limit temporal and collector biases. Each day, prior to continuing surveys of unsampled households, an attempt was made to inspect houses that were previously closed or where access had been refused. Access to houses was attempted a minimum of three times. Maynas and Tupac Amaru were surveyed on alternating days. This process was carried out until all the houses in each neighborhood had been surveyed or repeated attempts to gain access failed. In mid-December 1998, immediately after termination of the first survey, the sampling procedure was repeated. The second survey was completed on January 18, 1999. To differentiate data associated with the four different collections, the two surveys will be referred to as a (November–December) and b (December–January).
Entomologic surveys. Our survey methodology was based on techniques suggested by Focks and others.9 Briefly, after asking permission to survey the household, one member of the team administered a demographic survey designed to determine the number of occupants, dimensions of the property, house construction materials, method of cooking, water use patterns, type of sewage disposal, and insecticide use. Simultaneously, the other team member began collecting adult mosquitoes using a backpack aspirator (John W. Hock Company, Gainesville, FL).10 Aspiration collections were attempted in all rooms of the house (when permitted) including walls, under furniture, and inside closets and other likely adult mosquito resting sites. Aspiration collections were similarly attempted outside the house from outside walls, under eaves, vegetation, and in and around outdoor stored materials.
In our field laboratory, larvae were identified as Ae. aegypti by the relative size of the siphon and their movement compared with the other most commonly found Culex species.11 Limatus larvae were differentiated by the characteristics on the eighth tergite.11 All larval samples were cross-checked with the entomology collection sheets provided by the field team. Pupae were counted and placed in plastic emergence vials, 
30 per vial and labeled with the house, container code, and date. Each subsequent day, emerged adults were collected and placed in a -20°C freezer. After 30 minutes to 1 hour, their species was identified, counted by sex, and data were recorded on the entomology collection sheet.
Data management. A geographic information system (GIS), using ARC/INFO and ArcView software (Environmental Systems Research Institute, Inc., Redlands, CA), was developed for the city of Iquitos. A base map of city blocks in the form of AutoCAD files was obtained from the Peruvian Navy, which they created by digitizing ortho-corrected 1995 aerial photographs. The coordinate system and datum used were Universal Transverse Mercator and WGS-84, respectively. The AutoCAD files were converted to ARC/INFO export files and all polygons (city blocks) were closed using standard ARCEdit procedures. Files were then imported into ArcView and converted to shape files.
We then divided city blocks into individual housing lots that were identified by painted codes. The front end of each house lot was measured and recorded along with the house code and street address on a rough sketch of each block. Based on maps constructed in the field, each digital block in the GIS was split into lots of appropriate width using the measuring tool in ArcView. Lot length was estimated. Lot geometric centroids were then added to each individual lot and assigned a unique project code that was included on all subsequent survey forms. Construction of maps with resolution to the level of household lots allowed all entomologic data from the four surveys to be joined to geographic coordinates via house codes. Centroids allowed for spatial analysis to be performed from the level of the individual household upwards.
Analysis of the data. Spatial patterns of four variables were examined (adult Ae. aegypti, pupae, all water-holding containers, and water-holding containers positive for larvae and pupae). Variables were explored by identifying the spatial distribution of each of the variables for each of the two time periods. Our study focused on 1) each of the two neighborhoods as a whole, 2) the magnitude of each variable in each household for each neighborhood, and 3) the presence or absence of a variable in a household for each neighborhood. Global K-functions, point and weighted, were used to identify clustering for (1) and the local statistic, Gi*, was used for (2). These statistics are some of the suite of spatial statistical programs available as part of the Point Pattern Analysis (PPA) program. The program was developed by Arthur Getis with assistance from Laura Hungerford, Dong-Mei Chen, and Jared Aldstadt. An online version is available at http://zappa.nku.edu/~longa/cgi-bin/cgi-tcl-examples/generic/ppa/ppa.cgi. For (3), we used chi-square tests to compare similarities and differences among the various patterns.
K-functions.
Pattern models are based on the K-function work of Ripley12 and Getis.13 The K-function describes the number of pairs of observations between a point, which is the center of a disk and other points that are distance d away. For a stationary, isotropic process, 
(d) is the expected number of points within distance d of an arbitrary point. The estimator of is N/A where N is the number of points in the study area A.
The estimator of K(d) is
 |
where dij is the distance between the ith and jth observed points and Id(dij d) is an indicator function that is 1 if dij is less than or equal to d and 0 otherwise. For a circle centered on i passing through point j, uij is the proportion of the circumference of the circle that lies within A. When dij is less than the distance from i to one or more borders of the study area, uij is 1. The "border correction" makes 
(d) an approximately unbiased estimator of K(d) provided that d is less than the circumference of A. A square-root scale makes the function linear and stabilizes the variance. Thus, we have
 |
which is the estimator of L(d) 
[K(d)/
]. The mean of L(d) is d and the approximate variance is 
(N2).14 The expectation of L(d) given the hypothesis of complete spatial randomness (CSR) is d. CSR is a homogenous planar Poisson process where all points are independent of all other points and all locations are equally likely to contain a point. For CSR, a plot of 
(d) against d on similarly scaled axes yields a 45-degree line beginning at the natural origin. A clustered pattern occurs when (d) is greater than d and a dispersed pattern can be identified when (d) is less than d. In the spirit of a exploratory diagnostic tool, statistical significance at the P 0.05 level is assumed to exist when the observed (d) function falls outside of an envelope containing 19 permutations of the location of the N objects where each permutation is based on CSR. (d) is usually calculated for a series of distances d.
Instead of considering each point as a nominal scale variable, points can be weighted according to some measure of size or intensity,13
 |
where X is a random variable having values x for adult mosquitoes in houses at sites i. Equation (3) is the estimator for Lw(d), which is equal to E[w(d)]. In the cases discussed in this paper, the weights are in turn numbers of adult mosquitoes, pupae, water-holding containers, and positive containers. For each xi, there are (N - 1) values xj. In this case, the numerator of w(d) represents the product of the pairs of values xi xj within distance d of each x. The denominator is scaled such that if all x are of equal value, then (d) will be approximately equal to w(d). Thus, equation (3) represents a measure of clustering or dispersion identified in equation (2). If the number of adult mosquitoes, for example, is independently distributed within the plots of houses, (d) will be approximately equal to w(d). Upper and lower significance boundaries for w(d) can be determined by a permutation procedure in which the various observed values for number of adult mosquitoes, xi, are permuted among the house locations a specified number of times.
We also explored the increments to (d) and w(d) observed for each equal increase of distance. In a CSR pattern of adult mosquitoes, these successive values will be the same for each equal increase of d. The focus is on the non-cumulative properties of these pattern indicators. When the change in (d) is greater or less than the change in w(d) for a given distance band, the adult mosquitoes are less concentrated or more concentrated, respectively, than that expected in the observed pattern, no matter how clustered the pattern of houses. That is, the number of adult mosquitoes is not randomly distributed among the houses. In essence, we compare 
(d) with w(d) for a given small change in d.
Gi*(d) statistic. In addition to L(d), we used the local statistic, Gi*,15 to identify individual members of clusters. For Gi* we take each house as a center, one at a time, and search the nearby area for occurrences of more or fewer adult mosquitoes than expected. In this way, specific houses are identified as members or non-members of clusters. This statistic is written as
 |
where wij (d) is the i, jth element of a one/zero spatial weights matrix with ones if the jth house is within d of a given ith house; all other elements are zero; Wi* = 
wij(d), where wii is included, and S1i* = w2ij(d) (all j). The mean of the adult mosquitoes in houses is 
and s is the standard deviation. The value of Gi*(d) is given in normal standard deviates. Note that this statistic has as its expectation, Wi , which controls for the number of houses within d of each house. Note, too, that Gi*(d) is 0 in a pattern where adult mosquitoes are randomly distributed within d of house i. For this study, we arbitrarily define values greater than 2.575 (the 0.01 level of confidence) as representing houses which are members of clusters of adult mosquitoes. The statistics used in the analysis and the test criteria are summarized in Table 1.
|
|
RESULTS |
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
Neighborhood pattern analysis.
The results of the K-function analysis for adult Ae. aegypti in Maynas in time period a are shown in Table 2. Adult mosquito clustering occurs if values of (d) are higher not only than adult mosquitoes distributed at random in the Maynas neighborhood for a given distance (i.e., d), but also higher than the (d) value for the pattern of houses at that same distance. Clearly, it is not enough that adult mosquitoes were spatially concentrated at the same rate as the spatial concentration of houses. Note that in column 3 in Table 2, the w(d) value for adult mosquitoes at 10 meters is 22.86, which is quite a bit higher than the 10.00 (random expectation) shown in column 1. However, houses were much more clustered than random (16.33 versus 10.00 at 10 meters). Even so, adult mosquitoes were more clustered than houses. In addition, using 19 permutations to identify the range of possible values for adult mosquitoes among houses (at the 0.05 level), we find that adult mosquitoes at 22.86 fall outside of that range (low of 11.88 to high of 19.10) at 10 meters. This gives strong statistical evidence that adult mosquitoes were clustered in the Maynas neighborhood during time period a. Clustering is at the 10-meter level; thus, we can conclude that there is clustering around houses to at least 10 meters distant.
|
, as distance increases to 20 meters, 30 meters, and so on, the (d) values for houses increase at a rate that is not dissimilar from random expectation. This means that although houses are closely spaced at short distances, there is little or no increase in clustering as distance increases. The w(d) value for adult mosquitoes shown in column 3 at 20 and 30 meters, however, increases at a slightly higher rate than houses (column 5 versus column 4), indicating a continuing of the clustering identified at 10 meters to at least 30 meters. This pattern of increase changes by 40 meters (the increment is 10.55, less than the house increment of 14.15) indicating an end to the increase in clustering. That is, beyond 30 meters, any further clustering of adult mosquitoes corresponds to clustering of houses. We conclude that adult mosquitoes cluster heavily at nearest house distances and moderately to approximately 30 meters. In Maynas, the mean house width was 7 ± 3 meters; thus, adult clusters could extend to about 2 households on each side.
We altered equations (1) and (3) to include houses themselves; that is, we allowed i to equal j (Table 3; see Getis13 for an explanation of the methodology). Our focus now is on houses and their neighbors rather than neighboring houses only. In this circumstance, the clustering of houses (column 2) is inflated to include not only near neighbors at 10 meters, but also the houses themselves. The original value of 16.33 at 10 meters now increases to 21.44 for houses indicating that in this view, houses are more clustered than was indicated previously (an increase of 31%). More importantly, however, are the results when adult mosquitoes within houses are taken into account. Here the value at 10 meters increases to 39.30 from 22.86, an increase of 72%. The implication is that adult mosquitoes are heavily clustered within houses. Note also that as distance increases, the increment to houses and adult mosquitoes is approximately 10, indicating that there is a cessation of clustering beyond 10 meters. Again there is additional, albeit weak clustering up to 30 meters because the increase in the mosquito value is higher than that for the houses at 20 meters and 30 meters. These results taken together with the earlier ones unequivocally indicate that adult mosquitoes cluster heavily within or among nearest neighboring houses. In addition, there is evidence of further, albeit minor, clustering as far as 30 meters. The clustering within houses in the Maynas neighborhood quantitatively overwhelms this further clustering.
|
). If clustering was within households, the Gi* statistic will be above +2.575 at short distances, say 1 meter at the 0.01 level of statistical significance. If clustering continues to near neighbors within 10 meters of a house, the value of Gi* will be higher at 10 meters than at 1 meter. If values of Gi* do not increase with increases in distance, then whatever clustering existed at the shorter distance ceases to exist at longer distances. The houses that are members of significant clusters at 1, 10, 20, and 30 meters are shown in Figure 3. Note that of the 528 houses in Maynas during time period a, 35 (6.6%) are members of statistically significant clusters of adult mosquitoes. Of the 35, 10 exhibit clustering with near neighbors beyond the house itself. Of these 10, seven show clustering to 10 meters, two to 20 meters, and one to 30 meters. This result reinforces the notion that adult mosquitoes tend to cluster in single households with a modest spread to as far as 30 meters.
|
|
is a map of the presence of one or more mosquitoes in households. One hundred sixty-four (31.1%) of the houses had one or more adult mosquitoes present; however, only 35 of them (21.3%) were members of statistically significant clusters. This indicates that clusters were made up mainly of household concentrations, and that 79.7% of the households with mosquitoes were spread about in a random pattern among all households.
|
and 5 allow for the comparison of K-function values for water-holding containers, positive containers, and pupae with house and adult mosquito patterns in the Maynas neighborhood (equations 1 and 3). The d = 10 meters row in Table 4 shows, as before, that adult mosquitoes cluster more so than houses (22.86-16.33), but the pattern of water-holding containers and positive containers is more nearly like the pattern of houses (16.25–16.33 and 15.40–16.33). Thus, there is evidence of no clustering for these variables. In the case of pupae, however, there is a significantly lower value (12.03), indicating that pupae do not cluster beyond the household and, in fact, are dispersed rather evenly throughout the neighborhood. However, when we allow i to equal j (Table 5), pupae increase from 12.03 to 56.13, an extremely high and statistically significant value. This indicates that pupae cluster strongly within houses, but households infested with pupae are dispersed rather evenly throughout the neighborhood (Table 4).
|
|
and 5), we conclude that water-holding containers are ubiquitous in Maynas. That is, nearly all houses have water-holding containers. Conversely, containers positive for pupae and/or larvae are more concentrated in some houses than others and infested houses are dispersed evenly throughout the neighborhood.
Continuing on to 20 meters, 30 meters, and further (Tables 4 and 5), only pupae act differently than containers and positive containers. For both of the container variables, increases mirror those of houses, reinforcing our earlier results that show ubiquitous occurrences of these variables. Pupae values (Table 5), however, increase at a much slower rate than houses after 10 meters, indicating that households infested with pupae are less common than households with water-holding containers or positive containers, and that the spatial pattern of pupae is characterized by strong clustering within households.
Household pattern analysis of non-adult mosquitoes. Our Gi* statistic results show that there is a lack of statistically significant clustering beyond households for container and immature mosquito variables. In the case of pupae, there were 18 households exhibiting clustering with no clustering beyond the household. Of the 24 houses with clusters of containers, only two were clustered to a neighboring distance of 10 meters. For positive containers, 23 houses exhibit clustering, but only three of those were clustered beyond the household, two to 10 meters, and one to 20 meters.
Patterns of pupae: presence or absence in houses.
In this analysis, the concern is less with numbers of pupae in houses and more with their spatial occurrence in houses. Data in Figure 5 were derived from a Gi* analysis that assigned a 1 to houses with one or more pupae present and 0 for the absence of pupae. We found that 18 (3.4%) of the 528 houses can be considered as members of clusters at the 99% level of confidence. There are two distinct clusters: one in the middle block in the south and a smaller cluster in the north. These concentrations raise the question of the relationship of the location of pupae to adult mosquitoes.
|
Association among clusters.
In Table 6 we see, as before, that of the 528 houses in Maynas, 35 were members of clusters of adult mosquitoes and 18 were members of clusters of pupae in time period a. Only three houses were constituents of both clusters, a non-statistically significant result at the 0.05 level (
2 = 1.60, degrees of freedom = 1, Yates’ correction for small expectations). There was not a significant correlation between pupal and adult abundance within household or neighborhood clusters detected during the same survey.
|
reveals the overlap of households that have as few as one mosquito or one pupae present. Note that of the 528 houses in Maynas, 164 had at least one mosquito present and 155 had at least one pupae present in time period a. Expectation from a chi-square two-by-two contingency test indicate that the two types of occurrence come together in households 48 times. A total of 66 households were infested with both pupae and adults, demonstrating that the presence of these two variables are not independent (P < 0.01).
|
shows the results of Spearman’s rank correlation test where the number of water-holding containers per household were ranked from 1 to 14. Ranks 15 and 16 were made up of 15–19 and 20–35 containers, respectively. The final two ranks were grouped because of the few numbers of observations at these high levels. The mean number of adult mosquitoes per house was ranked for each container level. The result was a moderately high positive correlation for adults (+0.615, P < 0.05), and a modest correlation for pupae (+0.487, not significant). Our analysis indicates that elevated numbers of water-holding containers in houses increase the likelihood for elevated numbers of adult mosquitoes and/or pupae to be present.
|
1 positive container) was 45% in Maynas a, 38% in Maynas b, 29% in Tupac Amaru a, and 23% in Tupac Amaru b.
Clustering patterns of adult mosquitoes and pupae were consistent among the four surveys, but the level of clustering was greatest during the first Tupac Amaru survey. Table 9 shows the (d) values (i may equal j) for each of the four surveys for 10 meters. Houses in Tupac Amaru were slightly more clustered than in Maynas (25.00-21.44). Note also that in both neighborhoods water-holding containers are distributed much the same as were houses, but positive containers tend to cluster. Maynas with 29.05 and 31.00 in the two time periods are approximately 8–10 L units higher than the pattern of houses. Tupac Amaru with 38.56 and 44.30 are about 13–19 units higher than the pattern of houses. This implies that positive containers were more clustered in Tupac Amaru than Maynas, which may be a reflection of lower infestation rates in Tupac Amaru. Nevertheless, in both sites the level of clustering was relatively low.
|
Overall patterns of adult mosquito and pupae household clustering.
Table 6 shows the number of houses that were members of statistically significant clusters of pupae and adult Ae. aegypti. The number of houses included in clusters for pupae in Maynas decreased from 18 to 4 from time period a to b. Interestingly, the location of adult clusters changed between the two surveys. Twenty-eight households were members of adult clusters in the first Maynas survey that were not members of clusters in the second, a statistically significant finding that was not the case in Tupac Amaru. Only seven households were members of adult clusters in both Maynas surveys. Twenty Maynas households were members of clusters in the second but not first survey. The same type of result, changing cluster locations, was evident with member houses of pupae clusters. In Maynas none of the houses were members of pupae clusters in both surveys, whereas six households were part of pupae clusters during both time periods in Tupac Amaru (Table 6). This result indicates that the spatial distribution of entomologic data varies greatly within short periods of time.
Association among households having one or more of each variable present in each neighborhood over time.
Although clusters of positive containers, pupae, and adult mosquitoes identified by Gi* were not consistent with time, Ae. aegypti infestations of individual households were clearly a risk factor for future infestation. That is, there is evidence of repeat offenders. Table 7 shows the number of houses observed to be infested with either pupae or adults in survey a, survey b, or both. Pupae in a are again found in the same houses in b in both neighborhoods between 29% and 45% of the time, a statistically significant result. The implication is that for unknown reasons mosquitoes are more likely to lay eggs in containers on some house lots than others. Another risk factor for infestation is the number of water-holding containers in a household. Results in Table 8 indicate that there is a tendency for houses in both neighborhoods and both time periods to contain more pupae when more water-holding containers are present.
|
DISCUSSION |
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
The lack of spatial structure for immature forms of Ae. aegypti supports recommended vector surveillance strategies where standard sample size calculations and resource limitations are used to determine in a systematic way the number of houses to be sampled, typically every ith house. Our K-function analysis indicates that individual households are the appropriate spatial unit for entomologic surveys. From a temporal perspective because water-holding containers were ubiquitous in Iquitos, all households are at risk of infestation over any considerable period of time. Our results, however, imply that as the number of containers on a premise increases so does the risk of Ae. aegypti pupae and adult infestations. In other words, positive containers and pupae cluster within individual households, but the location of clusters changes through time. Biologically this makes sense. Infestation of a household is largely a function of container management practices by the occupants of the property and the ecology of Ae. aegypti egg-laying behavior. We did not detect larger scale structure that might have been affected by other factors (data not presented or discussed in this paper) such as the availability of piped water, local temperature, rainfall patterns, or garbage disposal.
Identification of "key premises" or households that are superproducers of Ae. aegypti has been proposed as a way to streamline surveys.18 The idea is that the presence of pupae or adults during an initial survey is a significant risk factor for observing the same life stage at the same location during subsequent surveys. If we adopt the notion of controlling key premises as a way of reducing but not eliminating Ae. aegypti populations, the fundamental need to refine our understanding of entomologic thresholds is reinforced. Until we quantitatively define the relationship between mosquito density and risk of virus transmission, we cannot predict the effect that eliminating key premises will have on the risk of human infection and disease. For example, eliminating key premises may not reduce the adult mosquito population below the threshold density and, depending on the nature of the relationship between virus transmission and vector density, the pattern of human infections could continue unabated. Interestingly, the transient pattern of immature mosquito cluster locations observed in our study indicates that even if key premises can be identified and eliminated there may still be a sufficient number of Ae. aegypti to sustain dengue virus transmission. It should be noted, however, that because Iquitos has a relatively low percentage Ae. aegypti production in permanent water holding containers, our results may be site specific. The same kind of thorough examination may need to be carried out (large sample sizes and spatial analysis) at other locations.
Although small, there was significant spatial structure of adult mosquito populations compared with pupae and positive containers. Adults cluster most to distances of approximately 10 meters and to a lesser extent out to 30 meters, which could include neighboring houses. This finding is consistent with our conclusion to use the household as the basic unit of entomologic surveillance. It also superficially supports focal insecticide treatments for dengue control, a practice in which households are treated with insecticides within a 50–100-meter radius of the residence of a detected dengue case.16 There are, however, at least three shortcomings to focal treatments that extend beyond spatial patterns of adult Ae. aegypti. The approach does not take into account 1) the time delay between when a person is infective to mosquitoes and they are detected as being clinically ill with dengue, 2) that infected people can transport virus rapidly over greater distances than flying infected mosquitoes, and 3) that viremic people can have an inapparent infection or may not seek medical assistance, the homes and surrounding areas of many people infective to mosquitoes will not be sprayed.
Our statistical approach corroborates results from mark-release-recapture experiments on the dispersal of adult Ae. aegypti. Most researchers have concluded that the typical flight range of this species is short (<100 meters). Rodhain and Rosen19 stated that spontaneous dispersal of adult Ae. aegypti averages from 30 to 50 meters per day, so that females are rarely expected to visit more than two or three houses in their lifetime. The length of an Ae. aegypti lifetime is difficult to estimate, but is generally believed to range from 8 to 16 days.9 Ordonez and others20 reported minimum and maximum daily flight distance for Ae. aegypti of 8 meters and 120 meters, respectively, with a mean of 30.5 meters. In a Kenyan village, McDonald21 found that most adult Ae. aegypti dispersed to less than 20 meters and the majority of those recaptured were collected in the same house where they were released. Edman and others22 similarly collected most of their recaptured Ae. aegypti in Puerto Rico from their release house. In Kenya, Trpis and Hausermann23 reported 57 meters as the mean daily flight distance for females, with a maximum dispersal of 154 meters. Sixty percent of their recaptured females were collected in 11 houses that were within 50 meters from their release point. Our spatial analysis agrees with the preponderance of evidence that in a place such as Iquitos most adult Ae. aegypti do not fly far from the container where they developed as larvae and pupae.
Spatial referencing of our adult survey data and application of statistical tools, such as K-function and Gi*, provided insights into adult dispersal behavior that help explain patterns of human dengue infections. We propose that over short periods of time the restricted flight range and frequent blood-feeding behavior of Ae. aegypti24 are underlying factors in the clustering patterns of human dengue infections. In addition to the studies cited above on Ae. aegypti dispersal, several researchers have reported spatial and temporal clusters of clinically ill dengue patients in the same household or adjacent houses.2,3,25,26 In the first spatial statistics analysis of this phenomenon, Morrison and others4 found that dengue cases reported within a three-day interval during an epidemic in Florida, Puerto Rico clustered up to 10 meters. With regard to blood-feeding behavior, Ae. aegypti is know to frequently and preferentially imbibe human blood meals24,27 and infected females can transmit dengue virus to as many as 20 consecutive hosts, one after another.28 It is conceivable that a single or very few infected Ae. aegypti that remain in the same general area could bite and transmit virus to several susceptible family members or their immediate neighbors within a period of a few days.
Upon further investigation, we may discover that the extent to which infected humans are clustered is influenced by house construction and distribution. For example, households in our study area were small and often located close together; most were row houses with common walls. Although features of housing in Iquitos might facilitate Ae. aegypti movement, we do not expect that the tendency for adult females to disperse will be dramatically different at other locations. In Iquitos, water-holding containers were found in all households surveyed, something that is expected to decrease the probability of female dispersal.22
Abundance of adult female mosquitoes should be the most appropriate measure of entomologic risk because they are in the life stage from which viruses are transmitted. Interestingly, in at least one previous study adult Ae. aegypti abundance was correlated with diagnosed dengue cases.29 The value of larval indices was recently challenged because their relationship with adult densities is questionable.6 Pupal indices are now being considered as alternatives to traditional larval indices.7,9 Pupal indices are attractive for three reasons. First, it is theoretically possible to make absolute counts of their abundance, something that cannot be done for flying and difficult to capture adults. Second, pupal mortality is low. The magnitude of the pupal population should, therefore, be directly and relatively easily correlated with adult densities. Third, because the pupa is the life stage that directly precedes the virus-transmitting adult, pupae should be a more direct measure of transmission risk than larvae, which are a developmental step removed from adults.
Results from our spatial analyses, however, identified some limitations of pupal indices. The transient nature and high variability of containers positive for pupae can lead to misleading survey results, especially if the goal is to identify "key premises" and if only a single survey is carried out. Examination of spatial correlations among water-holding containers, larvae, pupae, and adults reveal significant correlations between life stages that are directly linked in their developmental sequence. For example, larval clusters correlated with pupal clusters and pupal with adults, but larval clusters were not correlated with adult clusters. This indicates that many containers exhibited a cohort effect. That is to say, cohorts of mosquitoes in a given container move in synchrony through the different stages of their life cycle without overlapping other cohorts. A noteworthy observation in that regard is that we did not consistently collect all stages of mosquitoes at the same time in the same household. This indicates that containers in Iquitos are not in equilibrium with the mosquito population. Instead houses are positive for a limited period of time as mosquitoes develop, disperse, and the household reverts to being negative. Other households subsequently become positive and the process repeats itself. In locations where positive containers are ubiquitous and permanent a different pattern of cluster spatial stability may emerge.
We conclude that pattern analysis can efficiently describe local Ae. aegypti populations and substantially aid in our understanding of dengue epidemiology and the development of dengue surveillance and control strategies. We argue that development of long-term entomologic risk assessment strategies requires thorough surveys of all mosquito life stages. Our results highlight the importance of scale when investigating the dynamics of dengue transmission. In Iquitos, the appropriate scale for assessing mosquito vector density is the household level at frequent time intervals.
This work is being extended with more extensive studies in additional areas of Iquitos, including an entire city study of the affinity that Ae. aegypti may have for particular types of water-holding containers and the relationship of various measures of mosquito density to human dengue infection. In addition, related work is underway in Thailand, which will allow comparison of concepts and processes described for Iquitos to results from an ecologically and epidemiologically distinct study area.
Received March 8, 2003. Accepted for publication July 29, 2003.
Acknowledgments: We thank the residents of Maynas and Tupac Amaru, Iquitos, Peru for allowing us to work with them in their homes. We greatly appreciate support of the Loreto Regional Health Department, including Dr. Carlos Calampa, Dr. Jorge Reyes, Dr. Ruben Naupay, Dr. Carlos Vidal, Dr. Martin Casapia, and Dr. Moises Sihuincha, who have all facilitated our work in Iquitos. Helvio Astete and Gerson Perez Rodriguez supervised the collection and processing of mosquitoes. Entomologic surveys were carried out by Jimmy Maykol Castillo Pizango, Rusbel Inapi Tamani, Juan Luiz Sifuentes Rios, Nestor Jose Nonato Lancha, Federico Reategui Viena, Victor Elespuru Hidalgo, Edson Pilco Mermao, Abner Enrique Varzallo Lachi, Fernando Chota Ruiz, Angel Puertas Lozano, Guillermo Inapi Huaman, and Manuel Ruiz Rioja. Jimmy Roberto Espinoza Benevides and Fernando Espinoza Benevides entered data into our database. Jose Elespuro Bastos checked data entry for accuracy.
Financial support: This study was supported by grant AI-42332 from the National Institutes of Health.
Authors’ addresses: Arthur Getis and Kenneth Gray, Department of Geography, San Diego State University, San Diego, CA 92182-4493. Amy C. Morrison and Thomas W. Scott, Department of Entomology, University of California, 1 Shields Avenue, Davis, CA 95616-8584, Telephone: 530-754-4196, Fax: 530-752-1537, E-mail: twscott{at}ucdavis.edu.
|
REFERENCES |
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONREFERENCES |
|---|
This article has been cited by other articles:
|
|
 |
|
  J. Garcia-Rejon, M. A. Lorono-Pino, J. A. Farfan-Ale, L. Flores-Flores, E. D. P. Rosado-Paredes, N. Rivero-Cardenas, R. Najera-Vazquez, S. Gomez-Carro, V. Lira-Zumbardo, P. Gonzalez-Martinez, et al. Dengue Virus-Infected Aedes aegypti in the Home Environment Am J Trop Med Hyg, December 1, 2008; 79(6): 940 - 950. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C.-C. Kan, P.-F. Lee, T.-H. Wen, D.-Y. Chao, M.-H. Wu, N. H. Lin, S. Y.-J. Huang, C.-S. Shang, I-C. Fan, P.-Y. Shu, et al. Two Clustering Diffusion Patterns Identified from the 2001-2003 Dengue Epidemic, Kaohsiung, Taiwan Am J Trop Med Hyg, September 1, 2008; 79(3): 344 - 352. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. J. M. Koenraadt, J. Aldstadt, U. Kijchalao, R. Sithiprasasna, A. Getis, J. W. Jones, and T. W. Scott Spatial and Temporal Patterns in Pupal and Adult Production of the Dengue Vector Aedes aegypti in Kamphaeng Phet, Thailand Am J Trop Med Hyg, August 1, 2008; 79(2): 230 - 238. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. Chansang and P. Kittayapong Application of Mosquito Sampling Count and Geospatial Methods to Improve Dengue Vector Surveillance Am J Trop Med Hyg, November 1, 2007; 77(5): 897 - 902. [Abstract] [Full Text] [PDF]
|
|
|
|
