Figures
Abstract
Background
Aedes aegypti is the main mosquito vector of the four serotypes of dengue virus (DENV). Previous population genetic and vector competence studies have demonstrated substantial genetic structure and major differences in the ability to transmit dengue viruses in Ae. aegypti populations in Mexico.
Methodology/Principal Findings
Population genetic studies revealed that the intersection of the Neovolcanic axis (NVA) with the Gulf of Mexico coast in the state of Veracruz acts as a discrete barrier to gene flow among Ae. aegypti populations north and south of the NVA. The mosquito populations north and south of the NVA also differed in their vector competence (VC) for dengue serotype 2 virus (DENV2). The average VC rate for Ae. aegypti mosquitoes from populations from north of the NVA was 0.55; in contrast the average VC rate for mosquitoes from populations from south of the NVA was 0.20. Most of this variation was attributable to a midgut infection and escape barriers. In Ae. aegypti north of the NVA 21.5% failed to develop midgut infections and 30.3% of those with an infected midgut failed to develop a disseminated infection. In contrast, south of the NVA 45.2% failed to develop midgut infections and 62.8% of those with an infected midgut failed to develop a disseminated infection.
Conclusions
Barriers to gene flow in vector populations may also impact the frequency of genes that condition continuous and epidemiologically relevant traits such as vector competence. Further studies are warranted to determine why the NVA is a barrier to gene flow and to determine whether the differences in vector competence seen north and south of the NVA are stable and epidemiologically significant.
Author Summary
The Neovolcanic axis (NVA) traverses Mexico at the 19th parallel and is considered to be a geographic barrier to many species. We have demonstrated that the intersection of the NVA with the coast in Veracruz state is a barrier to gene flow in Ae. aegypti. This was unexpected because the intersection of the NVA with the Pacific Coast is not a barrier to gene flow. Further studies to identify the actual mechanism(s) that is(are) contributing to the lack of gene flow will provide important information on the trafficking potential of Ae. aegypti, which will be of great value to Ae. aegypti control programs. There are significant differences in vector competence for dengue virus between mosquitoes north and south of the NVA, but the epidemiological significance of these finding remains to be determined. Future studies will determine if, for example, the genes that condition midgut infection and vector competence of Ae. aegypti populations provide biomarkers for risk of dengue transmission. Such biomarkers could be of great value to control programs in resource limited environments by allowing targeting of vector control efforts to areas at most risk for epidemic dengue and dengue hemorrhagic fever.
Citation: Lozano-Fuentes S, Fernandez-Salas I, de Lourdes Munoz M, Garcia-Rejon J, Olson KE, Beaty BJ, et al. (2009) The Neovolcanic Axis Is a Barrier to Gene Flow among Aedes aegypti Populations in Mexico That Differ in Vector Competence for Dengue 2 Virus. PLoS Negl Trop Dis 3(6): e468. https://doi.org/10.1371/journal.pntd.0000468
Editor: Duane Gubler, Duke University-National University of Singapore, Singapore
Received: January 23, 2009; Accepted: May 27, 2009; Published: June 30, 2009
Copyright: © 2009 Lozano-Fuentes et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by grants AI-49256 and AI-45430 from the National Institutes of Health and in part by the Innovative Vector Control Consortium. SLF was supported by NIH Fogarty Center Training Grant AI-46753. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors declare that there are no competing interests.
Introduction
The mosquito Aedes aegypti is the main vector of the four serotypes of Dengue virus (DENV1-4). There are 50–100 million DENV infections each year [1],[2] and while most of these are mild or asymptomatic, the numbers of severe infections with shock and hemorrhage have increased dramatically in many parts of the world [3],[4]. Aedes aegypti populations exhibit a large amount of genetic variation in their ability to become infected with, propagate, and eventually transmit flaviviruses [5]–[8], including DENV1–4. Vector competence for flaviviruses is thought to be controlled by at least two physiological mechanisms, a midgut infection barrier (MIB) and a midgut escape barrier (MEB) [9],[10] with environmental factors contributing up to 60% of variation [9]. Our genetic studies suggested that infection rates among natural populations of Ae. aegypti may be due to segregation of alleles at up to 8 loci [11]–[14].
We previously conducted studies to determine the breeding structure and vector competence of Ae. aegypti populations in Mexico [15],[16]. For the population genetic studies, Ae. aegypti were collected from throughout the coastal regions of Mexico, and 25 haplotypes of the Nicotinamide Adenine dinucleotide dehydrogenase subunit 4 mitochondrial (ND4) gene were detected by SSCP analysis. These studies revealed that northeastern Mexican Ae. aegypti were genetically differentiated from the Yucatan and Pacific Coast mosquitoes. FST values revealed extensive gene flow along the Pacific Coast, but not in the Yucatan Peninsula and northeastern Mexico. These studies also revealed a barrier to gene flow somewhere along the Gulf of Mexico between Tuxpan and Moloacan/Minatitlan in northern and southern Veracruz State, respectively. Ae. aegypti collected for the population genetic studies were also phenotyped for vector competence for DENV2, which revealed considerable variation in vector competence for DENV2 in Mexico [8]. Interestingly, the Ae. aegypti collections from southern Veracruz Coastal Plain differed significantly in vector competence; mosquitoes from Merida, Chetumal, and Cancun in the Yucatan were the most vector competent and those from Nuevo Laredo and Houston, the least vector competent [8]. Unfortunately, in both the population genetic and vector competence studies, no sites were sampled in Veracruz state (∼750 km from north to south). This prevented us from identifying the specific barriers to gene flow.
This also prevented us from examining vector competence in Veracruz. This is of special interest because a 1986 serological survey conducted by the Secretaria de Salud of Mexico [17] revealed major differences in dengue seroprevalence rates in cities and towns in Veracruz state. The dengue seroprevalence rate was 58% (29/50) in people from Martinez de la Torre (in northern Veracruz) versus 0% (0/50) in samples from Moloacan (southern Veracruz).
The present study is therefore an attempt to define more precisely the geographic barrier to gene flow previously observed between the northern [16] and southern Gulf of Mexico Coastal Plain [15] and to characterize more thoroughly the vector competence of mosquitoes separated in southern Veracruz. We obtained 10 Ae. aegypti collections between Tuxpan in the north to Minatitlan in the south (Figure 1). Nine of these same 10 sites were resampled in 2004 to test the consistency of our 2003 results. These collections were analyzed with the same mitochondrial ND4 marker gene as in earlier studies [15],[16]. The same mosquitoes were assessed for VC and midgut infection and escape barriers using established protocols [8].
Pie charts indicate the proportion of mosquitoes that were vector competent (black), midgut negative (red) and head negative (green). The VC rates were interpolated by Inverse Distance Weighting and geographic areas are colored from yellow to red according to predicted vector competence rates. R2 = 0.66 and root mean square error = 9.6.
Methods
Mosquito collection
Mosquitoes were collected as larvae from the cities listed in Table 1. In each city multiple locations were visited (Figure 1) and at each location at least 3 separate breeding sites separated by at least 500 meters were sampled. These larvae were returned to the laboratory and emerged adults were individually examined to confirm that they were Ae aegypti. The 2003 collection was processed for analysis of vector competence and mtDNA markers; the 2004 collection only for mtDNA analyses. All experiments used F1–F4 mosquitoes to minimize effects of colonization and inbreeding.
Vector competence
The DENV2 strain used was dengue 2 JAM1409, which was isolated in 1983 in Jamaica and belongs to the American Asian genotype [18],[19]. Procedures for growing virus in 14 day cell culture, quantifying the virus and infecting mosquitoes with membrane feeders covered with sterile hog gut membranes are published [8]. A highly DENV2 susceptible Aedes aegypti colony called D2S3 [20] served as an internal control in each experimental feed to test for consistency in the titer and infectiousness of the DENV2 meal preparation. Undiluted virus titers ranged from 7.5–8.5 log10 infectious virus/ml, which resulted in infection of 100% of the D2S3 mosquitoes in each feeding experiment.
Fully engorged mosquitoes were removed from the feeding carton and held for 14-days at a constant 27°C and 80% relative humidity in an insectary with a 12-hour photoperiod. Mosquitoes were frozen at −70°C until processed. Heads and abdomen were assayed for infections by immunofluorescence assay (IFA) using a mouse derived primary monoclonal antibody directed against a flavivirus E gene epitope [21],[22]. DNA was then extracted from the thorax [23] for population genetic studies.
For IFA, detection of DENV antigen in head tissues revealed a disseminated infection; these mosquitoes were scored as head positive (H+). The H+ mosquitoes were considered to be vector competent (VC), because salivary glands become infected in disseminated DENV infections and the H+ mosquitoes are presumably capable of transmitting the virus [24]. If no viral antigen was detected in the head tissues, the mosquito was scored as head negative (H−) and vector incompetent (VIC). To determine the anatomic basis for VIC, the H− mosquitoes were then examined to determine if the midgut was infected. H− mosquitoes with no detectable antigen in the midgut were scored as having a midgut infection barrier (MIB). H− mosquitoes with detectable viral antigen in the midgut were scored as having a midgut escape barrier (MEB). Because mosquitoes with a MIB could not be phenotyped for a MEB, we also determined the overall head negative rate (H-R) = H−/N.
The 95% confidence interval around VC, MIB, MEB and H-R was calculated as the Wald interval:where
VC, MIB, MEB and H-R. Estimates are either adjusted by adding half of the squared Z-critical value (1.96) to the numerator and the entire squared critical value to the denominator before computing the interval [25].
Population structure
Primers used to amplify the ND4 and all the polymerase chain reaction (PCR) and Single Strand Conformation Polymorphism (SSCP) conditions were reported earlier [15],[16]. The ND4 PCR products from mosquitoes containing each of the 9 haplotypes were sequenced at least once along both strands using an ABI sequencer (Davis Sequencing, Davis, California). Products from at least two mosquitoes representing each haplotype were sequenced. These 20 sequences were compared to sequences reported previously and assigned the same numeric labels [15],[16]. Phylogenetic relationships among haplotypes have been previously described [15],[16].
Statistical analysis of mitochondrial haplotype frequencies
Variation in haplotype frequencies within and among collection sites and regions was examined using Molecular Analysis of Variance (AMOVA) [26]. Arlequin3 estimated pairwise Slatkin's “linearized FST” [FST/(1−FST)] [27] among collections and computed the significance of the variance components associated with each level of genetic structure by a nonparametric permutation test with 100,000 pseudoreplicates [26]. A distance matrix containing linearized FST values was collapsed to construct a dendrogram using unweighted pair-group method with arithmetic averaging analysis [28] in the NEIGHBOR procedure in PHYLIP3.5C [29].
Spatial analysis of vector competence
Inverse Distance Weighting interpolations are based on the assumption that the interpolating surface should be influenced most by the nearby points and less by the more distant points [30],[31]. The transformed VC values (arcsin√VC) were interpolated and the resulting surface was then back transformed. The maximum search area considered was 2.5° with no anisotropy (i.e. circular search area); the search was continued until five geographically most proximate collections (neighbors) were identified.
Results
Gene flow
The ND4 was amplified and surveyed for variation by SSCP analysis [32],[33] among 654 mosquitoes in 19 collections (Table 1). These were 10 collections obtained in 2003 and 9 obtained in 2004 (no mosquitoes were collected in Cosoleacaque). Nine different ND4 haplotypes were detected with SSCP. The ND4 gene was sequenced in 20 mosquitoes. All the sequenced haplotypes were compared to those previously reported (GenBank accession numbers AF334841–AF334865), and no novel haplotypes were detected. Accordingly all haplotypes in this study retain the same numerical designations as those in GenBank. As reported in previous studies [15],[16], sequences of mosquitoes with identical SSCP patterns were identical within each haplotype, and SSCP patterns differed among mosquitoes with one or a few nucleotide differences.
Figure 2 is a UPGMA cluster analysis of pairwise linearized FST values [27] among 46 collections including 19 from the present study, 12 from previous studies north of Panuco in 1996–1997 [16] and 15 from south and east of Minatitlan in 1998–1999 [15]. With the single exception of Nuevo Laredo, all collections in northern Veracruz fall within a single cluster. The genetic distinctness of Nuevo Laredo Ae. aegypti was previously reported for both RAPD and mtDNA markers [16]. Northern collections cluster independently of collection year. Figure 3 indicates that most mosquitoes in northern Veracruz have haplotypes 1–9 and that haplotypes 10–18 are absent.
With the exceptions of Minatitlan 1996, Moloacan 1998, and Zempoala 2003, all collections in southern Veracruz fall into a cluster that is distinct from the northern cluster (Figure 2). Figure 3 indicates that most mosquitoes in southern Veracruz have haplotypes 10–18, 19 and 24. In 2003 Zempoala (Figure 1) mosquitoes (indicated with an arrow in Figure 2) clustered with the northern collections and were most similar to Martinez de la Torre (Figure 3), the collection site just north of Zempoala. But in 2004, Zempoala mosquitoes clustered with the southern collections. In contrast to northern mosquitoes, southern collections cluster according to collection year with 1998–1999 collections clustering independently of the 2003–2004 collections. Figure 2 indicates that the Minatitlan 1996 and Moloacan 1998 collections are similar to one another but genetically very distinct from Minatitlan 2003–2004. Moloacan and Minatitlan are in close geographic proximity to one another, yet they are also close to Cosoleacaque, Acayucan, and Coatzacoalcos.
Haplotype number designations correspond to those in GenBank accessions AF334841–AF334865.
Nested analysis of haplotype frequencies
AMOVA [26] was used to compare haplotype frequencies 1) among all 46 collections in northern and southern Veracruz, 2) among the nineteen 2003 and 2004 collections in northern and southern Veracruz, and 3) in 2003 vs. 2004 collections (Table 2). When analyzing all 46 collections, a significant 16% of the variation in haplotype frequencies arose between collections in northern and southern Veracruz (Table 2) and an additional 20% arose among collections made either in northern or southern Veracruz. A similar pattern was detected when analyzing the 2003 and 2004 collections alone. However a greater percentage of the variation (24.5%) arose between collections in northern and southern Veracruz and, probably because we eliminated variation arising from the 1996–1998 collections, less variation (13%) arose among collections in northern and southern Veracruz. Whether collections were made in 2003 or 2004 made no difference. A negative and non-significant percentage of the variation arose between years.
Analysis of vector competence
Table 3 provides the results of the vector competence studies for DENV2 of mosquitoes from the 2003 collections. In each site we report the proportion of mosquitoes with virus in the head tissues (H+) or not (H−), and for H− individuals the presence of virus in the midgut (M+) or not (M−). The VC (H+/N), VIC (H−/N), MIB (M−/N) and MEB (H−/M+) rates were calculated for each population (Table 3). The VC and VIC rates as well as the MIB rate for each of the 10 populations are presented in pie charts in Figure 1.
In northern Veracruz, the VC rate ranged from 0.38–0.75 with an average of 0.55, while VC among mosquitoes in southern Veracruz ranged from 0.11–0.33 and averaged 0.20. These differences were significant (Wilcoxon tests for unpaired samples, p = 0.0095). This variation was attributable to the greater proportion of mosquitoes with both MIBs and MEBs in southern Veracruz. Mosquitoes with uninfected guts constituted 21.5% of collections in northern Veracruz while 45.2% of mosquitoes in southern Veracruz had a MIB. This 23.7% difference in MIB rate was significant (Wilcoxon tests for unpaired samples, p = 0.0191). MEB rate varied by 32.5% between northern MEB% = 30.3%) and southern (MEB% = 62.8%) collections ((Wilcoxon tests for unpaired samples, p = 0.0095).
The VC rates among the 10 Veracruz sites were interpolated by Inverse Distance Weighting (IDW) [30],[31] with (arcsin√VC)/100 using ArcInfo 9.1. The model derived by jackknifing over the 10 sites using the “leave-one-out” procedure had a R2 = 0.66 and a root mean square error of 9.6. The interpolated predicted values appear in colors from red (susceptible) to yellow (refractory) in Figure 1.
Both the original measurements of VC (pie charts) and the predicted values from IDW interpolation suggest that VC declines precipitously south of the intersection of the Neovolcanic axis (NVA) with the Gulf of Mexico coast (Figure 1). The overall pattern in VC among 34 collections of Ae. aegypti made over an 8 year period throughout Mexico and two sites in the southern United States (Figure 4) demonstrates that the VC of mosquitoes from Alvarado, Acayucan, Coatzacoalcos and Cosoleacaque is among the lowest in Mexico.
The bars represent mosquito susceptibility to DENV2 (Jam1409). The number next to the city name is the mosquito susceptibility ((H+/N)×100).
Discussion
Our results are consistent with an hypothesis that the intersection of the NVA with the Gulf of Mexico coast is the barrier to gene flow previously observed between Ae. aegypti collections north [16] and south on coastal plain along the Gulf of Mexico [15]. The Transverse Volcanic Belt of Mexico [34] divides the state of Veracruz into northern and southern Coastal Plains. This belt began to develop during the Oligocene and then later, during the Pliocene–Pleistocene, intense orogenic activity raised the Neovolcanic axis. The NVA extends from near the Pacific Coast east to the Gulf of Mexico and intersects the Atlantic coast in the state of Veracruz. The NVA favored a warm and dry climate in the south of Mexico, and promoted the establishment of tropical deciduous forests [35],[36]. Near the NVA, the onset of rainfall is earlier than in the semiarid highlands. Mexico's six highest mountains are part of the NVA, which constitutes the largest east west mountain range on the North American continent and have played an enormous role in vicariance and allopatric speciation events in a large number of plant and animal species [37]–[39]. We observed a local change in mitochondrial haplotype frequencies in Zempoala from a northern type of pattern (very similar to Martinez de la Torre) in 2003 to a more southern type of pattern (high frequency of haplotype 14) in 2004 (Figure 3). This suggests that local gene flow can occur across the NVA intersection. However, the overall pattern among collections made over an 8 year period (Figure 2) argues strongly that the narrow corridor between the NVA and the Atlantic Ocean restricts gene flow in the long term.
It isn't clear why the NVA acts as a barrier to gene flow in Ae. aegypti. We examined differences in climatic factors such as solar radiation, precipitation, and land use as potential barriers to gene flow in the Veracruz Coastal Plain, but there were no obvious consistent differences in these factors north and south of the NVA [40]. The NVA could serve as physical barrier to gene flow because the distribution of Ae. aegypti in Mexico is largely limited to elevations < 610 m (∼2,000 feet) above sea level [41]. However, elevations also exceed this limit where the NVA intersects the Pacific Ocean and mosquitoes from Tapachula north to Tucson Arizona appeared to represent a single panmictic population [15].
One major difference between the Pacific and Atlantic coasts of Mexico is the amount of movement of people and commerce. Aedes aegypti is generally considered to have low mobility via flight [42],[43], but is facilely moved about locally and globally through human transportation and commerce [14], [44]–[47]. In comparison with the Pacific Coast of Mexico, where there are major roads and railways (albeit not always on the edge of the coastal plain), the corridor between the Atlantic Ocean and the NVA contains only a single, two lane road that is used only for local travel. Most automobile and truck traffic between northern and southern Mexico goes through Mexico city [40]. The Pacific Coast also has robust maritime and cruise ship activity which may traffic Ae. aegypti along the coast. In contrast, there is little such activity between cities on the Gulf of Mexico, which could also limit gene flow [40]. Overall these considerations suggest that the principal barrier is human trafficking and commerce, but further investigations will be required to determine if this is true.
Unexpectedly the NVA is also associated with significantly different VC phenotypes. Figure 4 shows the overall pattern in vector competence of Ae. aegypti in Mexico and the southern United States and demonstrates that the low vector competence of mosquitoes from sites just south of the NVA is unusual. The reasons for this remain to be determined. The genetic mechanisms conditioning the differences in VC remain to be determined. Association mapping studies to determine if the early trypsin and late trypsin genes conditioned VC revealed no consistent associations between segregating sites in the genes and VC for DENV2 [48]. Importantly the DNA from each of the mosquitoes phenotyped for DENV2 VC has been archived and as new candidate genes for VC are identified, the potential role of the genes in VC can be rapidly tested using these materials. It is also important to note that these studies have been done with only one dengue virus serotype/genotype. It will be important to confirm these results with additional dengue serotypes and genotypes that are circulating in Mexico and Latin America [19].
The temporal stability of the VC patterns north and south of the NVA is of interest. VC for DENV2 appears to be a quantitative genetic trait with up to 60% of the variation in VC being associated with random, or uncontrolled environmental effects [9],[12]. QTL mapping of genome regions conditioning MIB and MEB have identified 8 different genome regions [11],[13],[14] three associated with a MEB, and five associated with an MIB and three of these mapping families originated from northeastern Mexico [11],[13],[14]. An ongoing reevaluation of VC in mosquitoes collected north and south of the NVA in 2005 indicates that the reduced VC south of the NVA is stable (S. Bernhardt, personal communication).
An interesting alternative hypothesis is that the patterns that we are detecting may have little to do with environmental and ecological factors and may instead represent the introduction of Ae aegypti formosus south of the NVA [46],[49]. The two subspecies are sympatric in Senegal [50] and other parts of West Africa, and could have been introduced independently and multiple times into the New World. We recently discovered chromosomal inversions in Senegalese Ae. aegypti formosus [51]. Such inversions might also act as barriers to gene flow if they condition prezygotic reproductive isolating mechanisms. If mating does occur then inversions might cause excess chromosome breakage following crossing over during meiosis in hybrids yielding aneuploid gametes. Aedes aegypti aegypti and Ae aegypti formosus differ dramatically in their vector competence for yellow fever virus [5],[7] and DENV [50]. Provocatively, the Ae. aegypti populations north and south of the NVA also differed significantly in VC for DENV2. However a focal distribution of Ae. aegypti formosus in southern Veracruz would not explain why no barriers to gene flow were detected between southern Veracruz collections and collections in the Yucatan Peninsula.
A major goal of our dengue research program in Mexico is to determine if mosquito VC is correlated with dengue incidence. If so, identification of genes that are biomarkers of VC could permit targeting of control efforts to areas at greatest risk for dengue epidemics. In this regard, it is intriguing that the Ae. aegypti south of the NVA exhibited low VC and some cities, for example, in Moloacan, the seroprevalence rate was zero [17]. This historical data, which is 20 years old, may not reflect current conditions in southern Veracruz. The epidemiology of dengue in Mexico has changed dramatically in the last two decades, and dengue is now hyperendemic in Coastal Plains of Mexico. A serosurvey for dengue antibodies conducted in Jatilpan, Veracruz, which is in the same region as Moloacan, revealed a seroprevalence rate of 80% [52]. However, in the Secretaria de Salud report there were cities/towns in southern Veracruz that had similar seroprevalence rates in 1986. Clearly this is a complex situation, and conducting prospective VC studies and seroprevalence surveys in cities and towns in this unique region. Such studies would provide important information on the importance of VC in dengue incidence.
Supporting Information
Alternative Language Abstract S1.
Translation of the Abstract into Spanish by Saul Lozano-Fuentes
https://doi.org/10.1371/journal.pntd.0000468.s001
(0.03 MB DOC)
Acknowledgments
We gratefully acknowledge the support of vector control programs in Veracruz State for assistance in collection of Ae. aegypti eggs. We also thank Cindy Meredith and staff of the Arthropod-Borne and Infectious Disease Laboratory for their assistance in Ae. aegypti population colonization and maintenance.
Author Contributions
Conceived and designed the experiments: SLF KEO BJB WCB. Performed the experiments: SLF. Analyzed the data: SLF WCB. Contributed reagents/materials/analysis tools: SLF IFS MdLM JGR KEO BJB WCB. Wrote the paper: SLF BJB WCB.
References
- 1. Gubler DJ (2002) The global emergence/resurgence of arboviral diseases as public health problems. Arch Med Res 33: 330–342.
- 2. Monath TP (1994) Dengue: the risk to developed and developing countries. Proc Natl Acad Sci U S A 91: 2395–2400.
- 3. Gubler DJ (2002) Epidemic dengue/dengue hemorrhagic fever as a public health, social and economic problem in the 21st century. Trends Microbiol 10: 100–103.
- 4.
Organization PAH (1994) Dengue and Dengue Hemorrhagic Fever in the Americas: Guidelines for Prevention and Control. Washington, D.C.
- 5. Aitken TH, Downs WG, Shope RE (1977) Aedes aegypti strain fitness for yellow fever virus transmission. Am J Trop Med Hyg 26: 985–989.
- 6. Rosen L, Roseboom LE, Gubler DJ, Lien JC, Chaniotis BN (1985) Comparative susceptibility of mosquito species and strains to oral and parenteral infection with dengue and Japanese encephalitis viruses. Am J Trop Med Hyg 34: 603–615.
- 7. Tabachnick WJ, Wallis GP, Aitken TH, Miller BR, Amato GD, et al. (1985) Oral infection of Aedes aegypti with yellow fever virus: geographic variation and genetic considerations. Am J Trop Med Hyg 34: 1219–1224.
- 8. Bennett KE, Olson KE, Munoz Mde L, Fernandez-Salas I, Farfan-Ale JA, et al. (2002) Variation in vector competence for dengue 2 virus among 24 collections of Aedes aegypti from Mexico and the United States. Am J Trop Med Hyg 67: 85–92.
- 9. Bosio CF, Beaty BJ, Black WC (1998) Quantitative genetics of vector competence for dengue-2 virus in Aedes aegypti. Am J Trop Med Hyg 59: 965–970.
- 10. Miller BR, Mitchell CJ (1991) Genetic selection of a flavivirus-refractory strain of the yellow fever mosquito Aedes aegypti. Am J Trop Med Hyg 45: 399–407.
- 11. Bennett KE, Flick D, Fleming KH, Jochim R, Beaty BJ, et al. (2005) Quantitative trait loci that control dengue-2 virus dissemination in the mosquito Aedes aegypti. Genetics 170: 185–194.
- 12. Black WC, Bennett KE, Gorrochotegui-Escalante N, Barillas-Mury CV, Fernandez-Salas I, et al. (2002) Flavivirus susceptibility in Aedes aegypti. Arch Med Res 33: 379–388.
- 13. Bosio CF, Fulton RE, Salasek ML, Beaty BJ, Black WC (2000) Quantitative trait loci that control vector competence for dengue-2 virus in the mosquito Aedes aegypti. Genetics 156: 687–698.
- 14. Gomez-Machorro C, Bennett KE, del Lourdes Munoz M, Black WC IV (2004) Quantitative trait loci affecting dengue midgut infection barriers in an advanced intercross line of Aedes aegypti. Insect Mol Biol 13: 637–648.
- 15. Gorrochotegui-Escalante N, Gomez-Machorro C, Lozano-Fuentes S, Fernandez-Salas L, De Lourdes Munoz M, et al. (2002) Breeding structure of Aedes aegypti populations in Mexico varies by region. Am J Trop Med Hyg 66: 213–222.
- 16. Gorrochotegui-Escalante N, Munoz ML, Fernandez-Salas I, Beaty BJ, Black WC (2000) Genetic isolation by distance among Aedes aegypti populations along the northeastern coast of Mexico. Am J Trop Med Hyg 62: 200–209.
- 17.
Gomez H (1992) Monografia sobre la epidemiologia del dengue [Monograph]. Distrito Federal, Mexico: Secretaria de Salud. pp. 11–59.
- 18. Deubel V, Kinney RM, Trent DW (1986) Nucleotide sequence and deduced amino acid sequence of the structural proteins of dengue type 2 virus, Jamaica genotype. Virology 155: 365–377.
- 19. Diaz FJ, Black WC, Farfan-Ale JA, Lorono-Pino MA, Olson KE, et al. (2006) Dengue virus circulation and evolution in Mexico: A phylogenetic perspective. Arch Med Res 37: 760–773.
- 20. Bennett KE, Beaty BJ, Black WC (2005) Selection of D2S3, an Aedes aegypti (Diptera: Culicidae) strain with high oral susceptibility to Dengue 2 virus and D2MEB, a strain with a midgut barrier to Dengue 2 escape. J Med Entomol 42: 110–119.
- 21. Gould EA, Buckley A, Cammack N (1985) Use of the biotin-streptavidin interaction to improve flavivirus detection by immunofluorescence and ELISA tests. J Virol Methods 11: 41–48.
- 22. Gould EA, Buckley A, Cammack N, Barrett AD, Clegg JC, et al. (1985) Examination of the immunological relationships between flaviviruses using yellow fever virus monoclonal antibodies. J Gen Virol 66(Pt 7): 1369–1382.
- 23.
Black WC, DuTeau NM (1997) RAPD-PCR and SSCP analysis for insect population genetic studies. In: Crampton J, Beard CB, Louis C, editors. The Molecular Biology of Insect Disease Vectors: A Methods Manual. New York: Chapman and Hall. pp. 361–373.
- 24. Salazar MI, Richardson J, Sanchez-Vargas I, Olson K, Beaty B (2007) Dengue virus type 2: replication and tropisms in orally infected Aedes aegypti mosquitoes. BMC Microbiol 7: 9.
- 25. Agresti A, Coull BA (1998) Approximate is better than “exact” for interval estimation of binomial proportions. Am Stat 52: 119–126.
- 26. Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131: 479–491.
- 27. Slatkin M (1993) Isolation by distance in equilibrium and nonequilibrium populations. Evolution 47: 264–279.
- 28.
Sokal RR, Sneath PHA (1963) Principles of Numerical Taxonomy. San Francisco: Freeman.
- 29.
Felsenstein J (2006) Phylogeny Inference Package. Version 3.66 ed. Seattle, WA: University of Washington.
- 30.
Shepard D (1968) A two-dimensional interpolation function for irregularly-spaced data. ACM. pp. 517–524.
- 31. Wartenberg D, Uchrin C, Coogan P (1991) Estimating exposure using kriging: a simulation study. Environ Health Perspect 94: 75–82.
- 32. Hayashi K (1991) PCR-SSCP: a simple and sensitive method for detection of mutations in the genomic DNA. PCR Methods Appl 1: 34–38.
- 33. Orita M, Iwahana H, Kanazawa H, Hayashi K, Sekiya T (1989) Detection of polymorphisms of human DNA by gel electrophoresis as single-strand conformation polymorphisms. Proc Natl Acad Sci U S A 86: 2766–2770.
- 34. Alaniz-Alvarez SA, Nieto-Samaniego AF, Ferrari L (1998) Effect of strain rate in the distribution of monogenetic and polygenetic volcanism in the Transmexican volcanic belt. Geology 26: 591–594.
- 35. Lozano-Garcia S, Sosa-Najera S, Sugiura Y, Caballero M (2005) 23,000 yr of vegetation history of the Upper Lerma, a tropical high-altitude basin in Central Mexico. Quat Res 64: 70–82.
- 36. Metcalfe SE, O'Hara SL, Caballero M, Davies SJ (2000) Records of Late Pleistocene-Holocene climatic change in Mexico - a review. Quat Sci Rev 19: 699–721.
- 37. Fa JE (1989) Conservation-motivated analysis of mammalian biogeography in the Trans-Mexican Neovolcanic Belt. Natl Geogr Res 5: 296–316.
- 38. Halffter G (1987) Biogeography of the Montane Entomofauna of Mexico and Central-America. Annu Rev Entomol 32: 95–114.
- 39. Huidobro L, Morrone JJ, Villalobos JL, Alvarez F (2006) Distributional patterns of freshwater taxa (fishes, crustaceans and plants) from the Mexican Transition Zone. J Biogeogr 33: 731–741.
- 40.
Lozano-Fuentes S (2004) Aedes aegypti vector competence and gene flow in Mexico. Association mapping software for testing candidates genes associated with a phenotype. [Dissertation]. Fort Collins, Colorado: Colorado State University.
- 41. Ibanez-Bernal S, Gomez-Dantes H (1995) Vectors of dengue in Mexico: a critical review). Salud Publica Mex 37: (Suppl)S53–S63.
- 42. Edman JD, Scott TW, Costero A, Morrison AC, Harrington LC, et al. (1998) Aedes aegypti (Diptera : Culicidae) movement influenced by availability of oviposition sites. J Med Entomol 35: 578–583.
- 43. Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, et al. (2005) Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg 72: 209–220.
- 44. Huber K, Le Loan L, Chantha N, Failloux AB (2004) Human transportation influences Aedes aegypti gene flow in Southeast Asia. Acta Trop 90: 23–29.
- 45. Lounibos LP (2002) Invasions by insect vectors of human disease. Annu Rev Entomol 47: 233–266.
- 46. Tabachnick WJ (1991) The yellow fever mosquito: evolutionary genetics and arthropod-borne disease. Am Entomol 37: 14–24.
- 47. Tabachnick WJ, Powell JR (1979) A world-wide survey of genetic variation in the yellow fever mosquito, Aedes aegypti. Genet Res 34: 215–229.
- 48. Gorrochotegui-Escalante N, Lozano-Fuentes S, Bennett KE, Molina-Cruz A, Beaty BJ, et al. (2005) Association mapping of segregating sites in the early trypsin gene and susceptibility to dengue-2 virus in the mosquito Aedes aegypti. Insect Biochem Mol Biol 35: 771–788.
- 49. Powell JR, Tabachnick WJ, Arnold J (1980) Genetics and the origin of a vector population: Aedes aegypti, a case study. Science 208: 1385–1387.
- 50. Sylla M, Bosio C, Urdaneta-Marquez L, Ndiaye M, Black WC IV (2009) Gene flow, subspecies composition, and dengue virus-2 susceptibility among Aedes aegypti collections in Senegal. PLoS Negl Trop Dis 3: e408.
- 51. Scott A, Bernhardt CB, Sylla M, Bosio C, Black WC IV (2009) Evidence of multiple chromosomal inversions in Aedes aegypti formosus from Senegal. Insect Mol Biol.. In press.
- 52. Navarrete-Espinosa J, Acevedo-Vales JA, Huerta-Hernández E, Torres-Barranca J, Gavaldón-Rosas DG (2006) Prevalencia de anticuerpos contra dengue y leptospira en la población de Jáltipan, Veracruz. Salud Publica Mex 48: 220–228.