|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
In Colombia, Ae. aegypti is the main mosquito vector transmitting dengue viruses. The incidence of DF and DHF in Colombia has increased dramatically since the re-invasion of Ae. aegypti following cessation of the eradication campaigns in the 1960s.6,7 Cali currently has hyperendemic dengue transmission,8,9 and dengue prevention is based on vector control directed principally at mosquito larvae using chemical and biologic insecticides. Adult mosquito control occurs primarily by chemical insecticides (organophosphates and pyrethroids) and is used only during epidemic situations. Additional efforts to control mosquito larvae use biologic control with copepods and fish. Educational programs are carried out principally in schools and by the local newspapers and television channels (Secretary of Health of Cali, Epidemiology Unit, Internal Annual Reports 1995, 1999, 2000, unpublished data). Although there is a standardized program in place, vector control activities change depending on the budget, availability of insecticides, and the political situation. There are no ongoing programs designed to evaluate the impact of these mosquito control activities.
The use of insecticides by mosquito control programs and by the general public has resulted in the emergence of insecticide resistance. Resistance to temephos, an organophosphate insecticide extensively used for larval control, was first identified by a bioassay in 1998 (Suarez MF and others, unpublished data). Later biochemical studies of Ae. aegypti larvae from Cali showed increasing frequencies of resistance-associated enzymes (non-specific esterases [NSE] and mixed function oxidases [MFO]), suggesting the emergence of more than one insecticide resistance mechanism (Ocampo CB and others, unpublished data). In this study, larvae from different locations in Cali showed significant differences in the types of enzyme-based resistance mechanism present and the percentage of mosquitoes showing resistance. In addition, resistance varied from one year to the next at the same location (Ocampo CB and others, unpublished data).
A potential explanation for the intersite variability observed in enzyme levels associated with insecticide resistance is the presence of genetically distinct Ae. aegypti populations at the different study sites. This could lead to the emergence of insecticide resistance at any one site that was dependent on the genetic structure of the population and the extent of local vector control activities (evolutionary selection by insecticide pressure). If a genetically subdivided population were present, an understanding of population dynamics would be essential to achieve effective vector control and could deter the emergence of insecticide resistance. Likewise, it is possible that such population level changes could affect other genetically controlled traits such as vector competence for dengue transmission. Genetic variation and variation in vector competence between closely situated populations of Ae. aegypti has been observed in other parts of the world.1021 These studies showed genetic variability and differences in vector competence in Ae. aegypti populations at a single point in time, but did not explain the dynamics of how these differences could occur. These studies suggested that insecticide pressure is a major selective force contributing to genetic variability. However, insecticide susceptibility, dengue vector competence, and population genetic variation were not examined concomitantly or in relation to temporal change.
The objective of this study was to examine the dynamics of Ae. aegypti populations in an urban environment (Cali) by investigating epidemiologically relevant biologic variables, insecticide resistance, and vector competence in conjunction with population genetic analyses. This information provides a basis for understanding patterns of transmission of dengue viruses, and for the interpretation and use of information obtained through vector surveillance.
| MATERIALS AND METHODS |
|
|
|---|
The objective of the study was to characterize the population dynamics of Ae. aegypti in a dengue endemic urban environment using epidemiologically relevant phenotypic characters (insecticide resistance and vector competence) in association with analysis of genotypic markers.
Collections of Ae. aegypti were carried out in 2000 and 2001 at five different study sites separated by a minimum distance of 4 km, which exceeds the documented flight distances of this species.2224 The rationale for this collection strategy was to increase the probability of obtaining genetically distinct populations of Ae. aegypti. The sites included the Cali residential areas of Salomia, Siloe, Agua Blanca, Puerto Mallarino, and Melendez (Figure 1
). Collections were performed during the rainy season (March-April) when mosquito abundance was highest and when the greatest number of genotypes should have been represented. A minimum of 10 larval habitats was sampled per site. These samples were obtained from different types of containers and on different weeks at each site to avoid collections of descendents from small numbers of females and to avoid the direct effects of any vector control activities in the area. For the genetic and biochemical analyses, field-collected larvae reared to five-day-old adults (F0) were used. A maximum of 10 adults were analyzed per larval habitat sampled. The bodies of the adult mosquitoes were used to evaluate insecticide resistance attributes and the legs were used to evaluate genotypic variation. Extraction of DNA from the legs rather than the whole insect decreased potential contamination by bacteria, gut contents, and symbionts, and allowed for the biochemical analyses and genetic characterization of each individual mosquito. For the vector competence studies, F1 and F2 generation progeny of F0 mosquitoes were reared under identical conditions and exposed to identical dengue infective blood meals.
|
5 days of age) were stored at 70°C until the assays were performed. These were used for detection of biochemical mechanisms of resistance to the principal insecticides applied for vector control (organophosphates, carbamates, and pyrethroids). The mechanisms evaluated included NSE, MFO, and acetylcholinesterases as previously described.2530 These biochemical techniques have been compared with bioassay methods and shown to identify mechanisms underlying resistance detected with bioassays.3135 The biochemical tests can detect resistance, determine resistance mechanism(s), and be performed on single field-collected specimens.36 Protein concentration was determined for each specimen using the Bradford method to adjust for differences in body mass of individuals that might require correction factors for the enzyme assays.25,26,37 Approximately 80100 mosquitoes were analyzed from each study site (this sample size is sufficient to identify differences
20% in the frequency of resistance between populations with a power of [1 - ß]
0.80). Individual adult mosquitoes were homogenized in 100 µL of 0.01 M potassium phosphate buffer, pH 7.2, and then suspended in 1 mL of the same buffer. Aliquots of 100 µL were transferred to microtiter plate (Dynatech Laboratories, Alexandria, VA) wells. Thirty samples were analyzed in triplicate per plate for each assay. Absorbance was measured using an enzyme-linked immunosorbent assay reader (MRX; Dynatech Laboratories). Replicates with a variability coefficient (SD/mean) greater than 0.15 were discarded to avoid differences produced due to handling errors. Three positive controls were present on all test microtiter plates. Because of the lack of highly resistant Ae. aegypti strains, the assay control was developed by constructing standard curves for each enzyme using the corresponding commercially available enzymes. The reference value for susceptibility was defined using the average absorbance value obtained from known susceptible mosquitoes (Rockefeller strain). The average/mean absorbance value of the blank control (reagents without mosquitoes) was subtracted from the values of all sampled mosquitoes to correct for background absorbance. Finally, each absorbance value was standardized as activity per unit protein using the results obtained in the protein assays.
Statistical analysis To compare Ae. aegypti populations for biochemical mechanisms of resistance, the mean absorbance values obtained for each individual population were compared with those of other populations using analysis of variance (ANOVA). Multivariate analyses were performed to compare variation in resistance values by site, sex, and year of collection using SPSS version 8.0 (SPSS, Inc., Chicago, IL). Results of the Ae. aegypti enzyme activity from the field samples from Cali were compared with those of the insecticide-susceptible Rockefeller strain and the most susceptible strain from Cali.
Vector competence The oral susceptibility to dengue 2 virus was evaluated in collaboration with the Arthropod Infectious Disease Laboratory at Colorado State University (Fort Collins, CO). Two dengue 2 (D2) strains (New Guinea C [New Guinea, 1944] and JAM1409 [Jamaica, 1983]) were used to evaluate vector competence. The JAM1409 D2 strain was used in the evaluations performed at Colorado State University in 2000 and the New Guinea C D2 strain was used for those conducted at Tulane University (New Orleans, LA) in 2001. Both strains were passaged in C6/36 mosquito cells prior to oral infection. For this study, the assumption was made that both strains had similar infectivity and would otherwise behave similarly in Ae. aegypti. All infections were conducted in a biosafety level-3 insectary at Colorado State University and Tulane University.
Adult mosquitoes collected from the Cali field sites were fed on live guinea pigs to obtain eggs (F1). Individuals from F1-F2 generations of Ae. aegypti from each study site were later exposed to D2 virus by artificial membrane feeding of virus in rhesus monkey blood.38 Within a given experiment, all mosquito colonies originating from different sites were exposed to the same D2 virus-blood mixture. Samples from all collection sites were included in each experiment. Blood-fed females were maintained in an insectary maintained at 85°F and a relative humidity of 80%, with a daily light:dark cycle of 16 hours:8 hours for 14 days (the extrinsic incubation period of the virus). Heads of the mosquitoes were tested for viral antigen by immunofluorescent antibody (IFA). Viral titers in the blood meal were quantified using the endpoint method described by Schoepp and Beaty39 to verify that the virus concentration was similar among blood meals. All specimens used in this analysis were collected and reared, and the progeny were obtained under similar conditions at similar times. The New Orleans strain of Ae. aegypti was used to standardize the methodology and as the control strain for the analysis of vector competence. Mosquitoes were infected by intrathoracic inoculations to provide positive controls for the IFA.40
Crosstabs analyses were performed within and among the populations collected during 2000 and 2001. The chi-square test was used for the statistical analysis.
Population genetics Aedes aegypti from four different sites in Cali were analyzed using a random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR).41 This technique was chosen for these studies because it detects large numbers of polymorphisms, and because it is possible to work with small amounts of DNA that can be obtained from parts of the mosquitoes (e.g., legs). From 33 to 40 mosquitoes (principally females) per site per year were analyzed. The mosquitoes were randomly selected from numerous larval habitats sampled at each site to minimize the possibility of collecting descendents of just a few females.
Random amplification of polymorphic DNA The DNA was isolated by homogenizing the legs of the mosquitoes using a Kontes pestle (Kimble/Kontes, Vineland, NJ) following the procedure described by Black and DuTeau.41 The samples were then stored at 70°C until amplification.
The PCR was set up as follows. Each 40-µL PCR mixture contained 5 µL of sample DNA, 200 µM of each dNTP, 1.5 mM MgCl2, 40 pmol of primer (final concentration = 1 µM), and 1.5 units of Taq polymerase (Platinum Taq DNA polymerase; Gibco-BRL, Gaithersburg, MD) in 1x buffer (10x buffer = 200 mM Tris-HCl, pH 8.4, 50 mM KCl). The four primers used included R07: 5'-AAG CCT CGT C-3', R32: 5'-CTT TCT ACC C-3', R22: 5'-GAT CAT AGC C-3', and R08: 5'-CAC ACT CCA G-3'. The DNA was then amplified in a thermal cycler (GeneAmp R, PDR System 2700; Applied Biosystems, Foster City, CA) using the following conditions: 95°C for two minutes (to activate the Platinum Taq DNA polymerase); then two cycles at 94°C for one minute, 30°C for two minutes, and 72°C for one minute; followed by 40 cycles at 94°C for 30 seconds, 40°C for two minutes, and 72°C for one minute; and one final cycle of 95°C for 30 seconds, 40°C for two minutes, and 72°C for five minutes. Samples were then held at 4°C. Negative controls were included in all runs. These constituted all components listed except the DNA template, which was replaced by 5 µL of sterile water.
Randomly amplified polymorphic DNA was evaluated by agarose and polyacrylamide gel electrophoresis as described by Black and DuTeau.41 The molecular weights of the scored bands were confirmed using molecular weight standards, and comparability between experiments was assured by analyzing products from different PCRs on the same polyacrylamide gels.
Analysis by RAPD Analysis of the results was based on the following assumptions: 1) the genomic regions amplified by RAPD-PCR segregate as dominant alleles; 2) each sub-population is in Hardy-Weinberg equilibrium; 3) alleles in homozygous recessive individuals are identical; and 4) dominant amplified alleles are similar. Since it is assumed that each subpopulation is in Hardy-Weinberg equilibrium, the frequency of recessive alleles (q) was estimated from the square root of the frequency of homozygous recessive (blank) individuals (q2). Population genetic analyses were carried out using the Wrights Fst estimate assuming the island model.14,42 The migration rate (Nm) was determined from Wrights Fst [Fst = 1/(4Nm + 1)]. For the analysis, gels of the RAPD products were scored for the presence or absence of the chosen bands (loci) and all data were stored in a Microsoft (Redmond, WA) Excel® spreadsheet. From here the data were imported into the FORTRAN programs RAPDBIOS, BIOSYS-2, RAPDDIST, and RAPDFST14 (Black WC IV, 1997, unpublished computer programs). RAPDBIOS and BIOSIS-2 were used to estimate the allelic frequencies, and those frequencies were used to determine which alleles would be included in further population analysis. RAPDDIST was used to compute the genetic distance among the populations. The diploid option was selected along with the Lynch and Milligan correction for small sample sizes and Neis genetic distance was calculated.42,43 RAPDFST was used with the original presence/absence matrix (excluding low variation alleles) to compute the inbreeding coefficient (Fst) and Nm.
| RESULTS |
|
|
|---|
Resistance to insecticides
Variation in enzymes that confer resistance to commonly used insecticides provided evidence of selective pressure for NSE and MFO among the Ae. aegypti populations studied (Figures 2
and 3
). These enzymes are associated with resistance to organophosphates, pyrethroids, and insect growth regulators, all of which were being used by vector control authorities and the general public in the municipality of Cali. In contrast, there was no evidence of increases in acetylcholinesterase (which acts at nerve synapses and leads to resistance to organophosphates and carbamates) as a mechanism of insecticide resistance.
|
|
Enzyme levels associated with the two biochemical mechanisms of resistance that had increased varied among populations by site and year of collection (P < 0.001, by ANOVA). Female Ae. aegypti showed significantly higher levels of enzyme activity than males (P < 0.001). However, the overall pattern of enzyme activity was similar in males and females from a given site.
Vector competence
The Ae. aegypti infection rate with D2 virus ranged from 19% to 60% infectious, as determined by IFA detection of virus antigen in head tissues (Table 1
). Vector competence varied among populations collected from different sites and by year of collection. In 2001, significant differences (P = 0.01) in vector competence were found among the populations from the five sites. Within-site variation in vector competence between 2000 and 2001 was significant only for Siloe (Table 1
).
|
Genetic markers Genetic variability analyses included the Ae. aegypti populations from Agua Blanca (east), Salomia (north), Melendez (south), and Siloe (west). The population obtained in Puerto Mallarino (northeast) was excluded because it was derived from a single yard with a large number of tires serving as breeding sites. Mosquitoes used in this analysis came from multiple larval habitats per site (as described in the Materials and Methods).
Only alleles with an average population frequency between 0.1 < P > 0.6 were used for the population analysis.44 Eleven of 21 alleles fulfilled this criterion. To determine differences between populations in different years, the following parameters were calculated: Neis genetic distance, Wrights Fst (estimates the degree to which Ae. aegypti is divided into sub-populations), and Nm. Two levels of organization were considered in conducting the analysis. First, spatial and temporal population variation was evaluated considering that all collections were made of independent populations (Table 2
, Analysis Parameter [AP] A and B). Second, the populations were grouped to analyze spatial and temporal variation in the overall population of Ae. aegypti sampled (Table 2
, AP C and D). Differences by location and time of collection are evident when the samples were analyzed as independent populations (AP A and B), but not when they are grouped by year or location (AP C and D). The values obtained for Neis genetic distance and the inbreeding coefficient (Fst) demonstrated a greater degree of differentiation among populations when they were considered independently. Similarly, gene flow was lower between populations (Table 2
). The population dynamics shown by this analysis is schematically represented in Figure 4
, which illustrates the spatial and temporal changes that could be occurring in Ae. aegypti populations from Cali or any other urban setting.
|
|
| DISCUSSION |
|
|
|---|
In this study of the population dynamics of Ae. aegypti from an urban setting, enzyme levels associated with insecticide resistance and vector competence varied independently at the different sites and times evaluated. As suggested in other studies, the variability in levels of MFO and NSE observed in this study is likely to be the result of selective pressure from insecticides used by vector control programs and the general public. In 2001, government-sponsored vector control programs were almost non-existent in Cali because of a lack of resources. Insecticide use by the general public should be considered as one of the major factors inducing selection pressures on Ae. aegypti populations during that year. Most commercially available insecticides have pyrethroids as their active ingredients. These compounds induce selection primarily on MFO enzymes4547; the enzymatic mechanism expressed at higher levels in Ae. aegypti from Cali in 2001. The observed variation in MFO levels in the Ae. aegypti from Cali could be explained by differing levels of insecticide selection pressure at the locations where collections were made. This suggests that in mosquito populations in Cali, and in other urban areas, insecticide resistance can vary from site to site at different times of the year. These results demonstrate that the Ae. aegypti population in an urban area is in continuous flux and can be genetically responsive to focal selection pressures. The population structure in turn can affect local vector control success and potential for pathogen transmission. However, high gene flow in association with decreased or varying selection pressures can also affect population genetic structure. Insecticide pressure is likely to constitute the major selective force affecting these populations, and as a result induces genetic changes through bottleneck and genetic drift effects.4850
Although this study found random changes in all the RAPD alleles evaluated (genetic drift), this does not exclude that hitchhiking events can occur on specific regions of the genome where resistance genes are present.11 This investigation did not identify specific alleles that influence insecticide susceptibility and vector competence. However, the findings provide a basis for understanding spatial and temporal variation in populations. Ideally, studies to understand the selection processes that are occurring in the populations should include specific molecular markers for insecticide resistance and vector competence. Unfortunately, such markers have not yet been identified, but availability of the genome of Ae. aegypti promises to make this possible in the near future.
The biochemical assays used in this study may not always accurately represent biologic resistance to a given insecticide.50 However, their use did reveal temporal and spatial changes in enzymes that have been linked to insecticide resistance. Biochemical assays were selected to measure the potential for development or loss of resistance because they can be linked to genetic change, and they allow individual mosquitoes to be examined, which is impossible using most bioassay methods. Ideally, in a vector control setting, the biochemical tests initially would be performed along with bioassays to assess accuracy for predicting real biologic resistance. After the initial comparison, biochemical tests could be used alone to monitor ongoing effectiveness of vector control efforts, with bioassays used to confirm any overt changes in susceptibility.51
The vector competence analysis is limited because it is only one of several variables that influence vectorial capacity (overall ability of a vector to transmit a pathogen). Spatial and temporal fluctuations in Ae. aegypti abundance and life span, human anti-mosquito behavior, and vector control activities can also influence virus transmission. However, vector competence itself is important in transmission efficiency and can influence epidemiologic situations.52 Another limitation of the vector competence evaluation was that two different strains of dengue 2 virus were used (New Guinea C in 2001 and JAM1409 in 2000). However, since only one (Siloe) of the five sites sampled showed significant differences in vector competence between populations collected in 2000 and 2001 (Table 1
), the assumption that both strains behaved in a similar manner in the mosquitoes appears reasonable. If they behaved differently, then all the sites should have shown population differences between collections made in 2000 and 2001.
The use of RAPDs for population genetics allows analysis of small amounts of DNA and the detection of a large number of polymorphisms. Nevertheless, this methodology may be limited by the difficulty to precisely replicate results because of bands associated with contaminants, and by the assumptions made in conducting the data analyses. To decrease potential methodologic limitations, the DNA was extracted from legs to avoid gut contaminants, and at least one sample per PCR set was repeated. In addition, the bands were selected from a representative mosquito population sample and only well defined bands were scored.
The present study demonstrated that Ae. aegypti from highly populated urban areas could be considered as a population with multiple local populations that can be monitored to control the emergence and spread of insecticide resistance, especially when focal outbreaks of dengue are detected. Frequent monitoring of insecticide resistance or the biochemical surrogate allows mosquito control programs to determine when insecticides need to be alternated to mitigate the onset of resistance. Evidence of high gene flow between populations within the urban area also indicates that focal insecticide resistance must be rapidly managed by eliminating the insecticide pressure completely, or by changing the class of insecticide used. Because changes in vector competence and enzyme levels associated with resistance can change over relatively short periods of time, the success of vector control activities and the ability to combat dengue epidemics in urban areas depend on frequent surveillance for insecticide resistance in appropriately selected localities. While currently this strategy of intensive surveillance would be difficult for many countries to achieve on a national basis, there are municipalities within countries that do have the resources to initiate this type of surveillance. The relevance of this type of surveillance is highlighted by the fact that insecticide resistance in mosquitoes continues to be reported from many areas around the world. Once the laboratory tools are in place, it should be possible for municipalities to initiate and continue periodic surveillance for insecticide resistance. The results presented in this paper provide a basis for surveillance strategies to take into account the complexity of urban Ae. aegypti populations and the dynamics of resistance.
Received June 16, 2003. Accepted for publication October 22, 2003.
Acknowledgments: We are grateful to Dr. Nancy Gore Saravia for her valuable comments and advice during the preparation of this manuscript. We thank Dr. Ken Olson for reviewing the manuscript and, together with Dr. Kristine Bennett (Colorado State University Fort Collins, CO), for advice and technical assistance in the assessment of vector competence studies. We also thank Dr. William Brogdon (Centers for Disease Control and Prevention, Atlanta, GA) for his advice on biochemical tests and for reviewing the manuscript, Dr. Patricia Dorn (Loyola University, New Orleans, LA) for her advice on population genetic analysis, and the personnel of the Centro Internacional de Entrenamiento e Investigaciones Medicas for logistic support.
Financial support: This work was partially supported by a Doctoral Fellowship from the Instituto Colombiano para el Desarrollo de la Ciencia y la Tecnologia Francisco José de Caldas (COLCIENCIAS), Colombia government, 19962000, the World Health Organization Special Program for Research and Training in Tropical Diseases (TDR), 2000, a TDR Research Training Grant no. M8/181/4/0.176, a Tulane Department of Tropical Medicine Scholarship, 20002002, the Lewis T. Graham Student Research Award, the Louisiana Mosquito Control Association, 2001, and the Coypu Foundation of Louisiana, 2002.
Authors addresses: Clara B. Ocampo, Centro Internacional de Entrenamiento e Investigaciones Medicas, Av. 1a N #3 - 03, Cali, Colombia, Fax: 57-2-667-2989, E-mail: claraocampo{at}cideim.org.co. Dawn M. Wesson, Department of Tropical Medicine, Tulane University, 1430 Tulane Avenue, SL17, New Orleans, LA 70112.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
V. M. Scarpassa, T. B. Cardoza, and R. P. Cardoso Junior Population Genetics and Phylogeography of Aedes aegypti (Diptera: Culicidae) from Brazil Am J Trop Med Hyg, June 1, 2008; 78(6): 895 - 903. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. MENDEZ, M. BARRETO, J. F. ARIAS, G. RENGIFO, J. MUNOZ, M. E. BURBANO, and B. PARRA HUMAN AND MOSQUITO INFECTIONS BY DENGUE VIRUSES DURING AND AFTER EPIDEMICS IN A DENGUE-ENDEMIC REGION OF COLOMBIA Am J Trop Med Hyg, April 1, 2006; 74(4): 678 - 683. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |