Skip to main content

Genetic variability of Aedes aegypti (Diptera: Culicidae) in El Salvador and Honduras: presence of a widespread haplotype and implications for mosquito control

Abstract

Background

This study examined population genetics of Aedes aegypti in El Salvador and Honduras, two adjacent countries in Central America. Aedes aegypti is associated with yellow fever, dengue, chikungunya, and Zika. Each year, thousands of cases of dengue are typically reported in El Salvador and Honduras.

Methods

In El Salvador, collections were obtained from five Departments. In Honduras, samples were obtained from six municipalities in four Departments. Mitochondrial DNA cytochrome oxidase I (COI) was sequenced, and consensus sequences were combined with available sequences from El Salvador to determine haplotype number, haplotype diversity, nucleotide diversity, and Tajima’s D. A haplotype network was produced to examine the relationship between genotypes.

Results

In El Salvador, there were 17 haplotypes, while in Honduras there were 4 haplotypes. In both El Salvador and Honduras, Haplotype 1 is most abundant and widespread. In El Salvador, haplotype H2 was also widespread in 10 of 11 sampled municipalities, but it was not present in Honduras. The capital of El Salvador (San Salvador) and the eastern region of ES had the highest haplotype diversity of regions sampled.

Conclusions

Haplotype 1 and H2 each belong to different phylogenetic lineages of Ae. aegypti. The most geographically widespread haplotype (H1) may have been present the longest and could be a remnant from previous eradication programs. These data may contribute to future control programs for Ae. aegypti in the two countries.

Graphical Abstract

Background

Aedes aegypti is one of the primary vectors of dengue, Zika, chikungunya, and yellow fever viruses. The mosquito originates from Africa, has been transported globally, and is now established in many tropical, subtropical, and even temperature regions [1, 2]. In addition to vectoring infectious disease, this mosquito is associated with human habitation. Areas with artificial water storage around homes and buildings provide ideal breeding sites, making this species a pest globally [1].

Aedes aegypti populations from Africa group into two mitochondrial lineages, and globally the two lineages are widely distributed with genetically distinct populations [3]. Population genetic studies have found that insects including mosquitoes have divergent populations or strains, which vary in their vector competence [4,5,6]. Populations of Ae. aegypti can survive a range of environmental conditions [7]. Determining the population genetics of a species, and the distribution of biotypes and strains and their biological attributes, can contribute to management of vectors and disease prevention. Currently, novel control approaches for Ae. aegypti such as sterile insect technique (SIT) and Wolbachia-infected mosquitoes used for incompatible insect technique (IIT) are being developed, implemented and evaluated worldwide [8, 9]. These innovative approaches, along with insecticide control of mosquitoes, may benefit from knowledge of the genetic types or strains of Ae. aegypti present in a region.

While Ae. aegypti was previously eradicated in much of the Western Hemisphere in the 1950s–1960s, including in Central America in El Salvador and Honduras, it has since reemerged as a public health threat [10,11,12,13,14]. Aedes aegypti was reintroduced and periodic outbreaks of dengue have been occurring since 1980 [11]. The mosquito is currently widespread in El Salvador and Honduras [15, 16]. In Central America, 4 to 5 million people annually are estimated to be infected with dengue [17]. In 2011–2013, over 29,000 cases of dengue were reported in El Salvador [18]. From 2005 to 2014, El Salvador reported 238,232 cases of dengue and Honduras 268,411, without considering the additional public health burden of Zika and chikungunya which emerged in 2015 [19]. Aedes aegypti-transmitted infectious disease cases have continued to increase in recent years in both countries as well as in South America [20, 21], and effective vector control interventions for this mosquito in this region are limited yet remain of the utmost importance [15,16,17,18].

Population genetics of Ae. aegypti have been investigated globally, within the Western Hemisphere, and a few regional studies have been conducted in Central America. In Central America, a study in Panama found that Ae. aegypti had high genetic diversity and primarily grouped into two genetically divergent clusters associated with wet and dry climates [22, 23]. In El Salvador, populations of Ae. aegypti from six regions were similarly found to consist of two lineages and a high genetic diversity [15]. In contrast, a study in Honduras of Ae. aegypti uncovered few haplotypes, in contrast to other studies in the Western Hemisphere [16].

The objective of this study was to further examine the population genetics of Ae. aegypti in El Salvador and in Honduras, as they are two adjoining countries in Central America. El Salvador is on Pacific Coast of Central America, while Honduras borders the Caribbean on the east coast. Both countries have major port cities with international commerce and movement of goods and international tourism and migration. New collection sites were obtained from El Salvador and Honduras for examining Ae. aegypti genetic diversity. The number of lineages, genetic haplotypes, and nucleotide diversity of Ae. aegypti will be compared between samples from the two countries.

Methods

Study Sites Collection sites were in El Salvador and Honduras, Central America (Figs. 1, 2, 3). Samples were collected in El Salvador from January to March 2018. Larvae were collected in the following five Departments: Santa Ana, Ahuachapan, Cabañas, San Vicente, and La Unión (Table 1, Fig. 2). These municipalities were selected to obtain data on genotypes of Ae. aegypti present in additional regions of the country which had not been previously sampled [15]. Santa Ana is a department on the border of Guatemala and the most populated city in the west of the country. Ahuachapan also borders Guatemala, and collections were made at 1240 m above sea level. Samples in the department of Cabañas came from Sensuntepeque in the northeastern part of the country, a city geographically isolated from the rest of the territory. The collections made in San Vicente were concentrated in communities near the Pan American highway and near the Lempa River. Finally, La Unión is an old port city, located on the coast of the Gulf of Fonseca, with hot arid weather conditions (Table 1, Fig. 2).

Fig. 1
figure 1

Map of Central America including El Salvador and Honduras adapted from ArcGIS 10.8 shape file for Central America

Fig. 2
figure 2

Map of El Salvador with boundaries of each Department. New samples were obtained from five departments indicated by a mosquito icon, and existing sequences were available from an additional six departments and are represented by a DNA sequence. Overall, data used came from 11 of 14 departments throughout the country. Map is modified from ArcGIS 10.8 shape file for El Salvador

Fig. 3
figure 3

Map of Honduras with boundaries of departments. Samples were included from six municipalities in four departments, Cortés, Francisco Morazán, El Paraiso, and Choluteca, indicated with a mosquito icon (see Table 1). Map is modified from ArcGIS 10.8 shape file for Honduras

Table 1 Collection sites for Aedes aegypti from El Salvador and Honduras

In each department, one or two neighborhoods were typically visited to collect larvae, with at least two or three collection sites per neighborhood (Table 1). Sample sites included barrels of standing water and outdoor sinks called pilas. The samples at each site typically consisted of a few larvae of the numerous larvae which were present to reduce the chance that the larvae collected were siblings. A previous study found that the mean number of families represented per oviposition site for Ae. aegypti was 4.7 [24]. After collection, larvae were transported to the laboratory for rearing into adults at the University of El Salvador. Adults were identified to species using taxonomic keys [25], and females were used for DNA extraction.

Additional sequences included in analyses came from a previous study in El Salvador [15], which collected samples between May and August 2014 (Table 1). Larvae of Ae. aegypti were collected in six departments: Sonsonate, San Salvador, Chalatenango, Usulután, San Miguel, and Morazán (Fig. 2). From the prior study of six departments, 82 sequences from GenBank were included in analyses [Table 3, GenBank accessions, Usulatán (MK028219-32), Morazán (MK028233-52), Chalatenango (MK028253-62), San Miguel (MK028263-78), Sonsonate (MK028279-90), and San Salvador (MK028291-300)].

From Honduras, samples came from six collection sites in four departments. In Cortés, three neighborhoods were included, while the other departments (Francisco Morazán, El Paraíso, and Choluteca) each had one neighborhood sampled (Table 1, Fig. 3). The sites were located along a central transportation corridor in Honduras which runs through several types of ecosystems. Cortés department is a tropical region in proximity to the Caribbean Sea. Tegucigalpa is the capital city at 990 m, which has dry tropical forest landscape, and Choluteca is closer to the Pacific (Fig. 3). Colonies were formed from the initial collections of eggs from ovitraps. Eggs in initial collections were reared to adults and then identified using a taxonomic key [25]. In 2018, eggs were randomly sampled from colonies to obtain a sample of 10 eggs from each population, which were individually extracted to obtain DNA.

DNA extraction

In El Salvador, individual adult females (n = 45) were used for DNA extraction, while from Honduran samples, DNA was extracted from eggs (n = 45) (Table 1). DNA extraction was completed with the Qiagen DNEasy® Blood and Tissue Extraction Kit following standard protocols [26] and using an overnight incubation of samples at 65 °C. The quantity of DNA in each sample was measured using the Qubit 2.0 fluorimeter Hs DNA kit (ThermoFisher, Waltham, MA, USA).

Mitochondrial DNA COI

For each insect, DNA was used to sequence a ~ 650-bp region of mitochondrial DNA cytochrome oxidase 1 (CO1) (known as the ‘bar code’) using a universal forward primer LCO 1490 (5′-GGTCAACAAATCATAAAGATATTGG-3′) and reverse primer HCO2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA-3′) [27, 28]. A polymerase chain reaction (PCR) mix for each sample consisted of the following: 40 µl sterile ultra-pure water, 1 µl Taq polymerase (Clonetech, Mountainview, CA), 5 µl Taq 10 × buffer, 1 µl dNTPs, 1 µl forward primer, and 1 µl reverse primer, and each reaction had 2 µl template DNA added. The PCR program was the following: an initial 1 min warm-up at 95 °C; then 40 cycles of a touchdown program consisting of 92 °C for 30 s, 43–52 °C for 30 s (with a 0.3 °C temperature increase each s), and 72 °C for 60 s; after 40 cycles, a 68 °C final extension for 10 min and then a hold at 4 °C. Samples were cleaned up using the Exo-sap-it (Affymetrix, Inc., Santa Clara, CA) cleanup kit and run on a 3730 Genetic Analyzer.

Data analyses

Resulting sequences were analyzed using Geneious software v.9.1.6 to produce consensus sequences [29]. Sequences were trimmed, forward and reverse sequences were aligned, and a consensus sequence was produced. Sequences were aligned in Geneious and exported to Mega for further analysis [29,30,31].

The number of segregating sites (S), number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (Pi), and Tajima’s D were generated using DnaSP v.6.12.03 [32]. Tajima’s D was calculated to test whether there was a departure from neutrality, such as a population expansion or contraction. A median joining haplotype network tree of Ae. aegypti was constructed using PopART [33].

The pairwise Fst was calculated in Arlequin v.3.5.2 to estimate population differentiation based on differences in haplotype frequencies [34]. Analysis of molecular variance (AMOVA) was conducted to examine the distribution of genetic variation within and among populations using 1000 permutations with Arlequin v.3.5.2 software, and additional analyses were also run to examine variation within each country [34, 35].

Results

Ninety new sequences were produced, 45 from El Salvador and 45 from Honduras (Tables 1, 2). These sequences were combined with 82 existing sequences from El Salvador for an analysis of 172 sequences from 11 municipalities in El Salvador and 6 municipalities in Honduras (Table 2). A total of 19 different haplotypes were identified in the two countries (Table 3, Fig. 4). When combining all samples from both countries, the higher frequency, broadly geographically distributed haplotypes were H1 (n = 94/172) and H2 (n = 40/172) (Fig. 4). The haplotype network shows the genetic relationships among the haplotypes. Of the 19 haplotypes, 12 (63.2%) were unique to single populations, and 7 haplotypes (36.68%) were shared among populations. Of the 19 haplotypes in the study, only two occurred in Honduras, which were not found in El Salvador (H10, H11).

Table 2 Number of individuals used from each department in El Salvador or Honduras for mtDNA COI sequences
Table 3 Number of segregating sites (S), number of haplotypes (H), haplotype diversity (Hd), and nucleotide diversity (Pi) for each of the 17 regions included from El Salvador and Honduras
Fig. 4
figure 4

Median joining haplotype network tree of Aedes aegypti. The tree was constructed using 19 haplotype CO1 mtDNA sequences from 17 sites in El Salvador and Honduras. The size of each circle is proportional to total number of samples of each haplotype. Black circles represent missing or unsampled steps. The site name abbreviations are described in Table 1

In El Salvador, 10 haplotypes (H5–7, H9, H14–19) were found only in one population each, and 7 (H1–4, H8, H12–13) were shared among populations. The highest haplotype numbers and diversities were observed in three municipalities in the eastern portion of El Salvador, in San Miguel, Morazán, La Unión, and in the capital San Salvador (Table 3). Interestingly, lower haplotype diversity was observed in the Western portion of the country, in Santa Ana, Sonsonate and Ahuachapan, all having only the two most abundant haplotypes, H1 and H2. For Santa Ana, two third of samples were H2 (11/16), while Sonsonate was almost entirely H1 (11/12), and Ahuachapan was a mix of both H1 and H2 (Table 3, Fig. 4).

In Honduras, there were a total of four haplotypes; one haplotype (H1) was shared across all six populations, and three were only found in one population each. In Honduras, nearly all the samples (93%, 42/45 samples) were haplotype H1 (Fig. 4). This same haplotype was widespread as well in El Salvador in 9 of 11 municipalities and comprised 51% (52/127) of samples in El Salvador.

In El Salvador, a second haplotype H2 was also widespread in 10 of 11 municipalities and comprised 31.5% (40/127) of samples. Haplotype 2 was more frequent in individuals from Santa Ana and San Vicente and was also present at the rest of the localities of El Salvador, except inland in San Miguel (Fig. 4). Haplotypes 4, 8, and 13 were the next most widespread haplotypes and were shared by individuals from the SM, JI, LU, SV, SS, and SFG localities (Fig. 4). At the periphery of the network were haplotypes H3, H9, H10, H12, H13, H15, H18, and H19.

Assessment of population expansion based on the neutrality test found most values were not statistically significant. The overall Tajima’s D for El Salvador populations was 1.85, which was positive but nonsignificant. In El Salvador, only Sonsonate (ES) had a negative significant Tajima’s D (−2.016). In Honduras, Tajima’s D overall was −2.0 (Table 3). In addition, Tajima’s D values were negative and statistically significant in Villanueva (HON) (D = −1.701) and Choluteca (HON) (D = −1.823) (Table 3).

The Global AMOVA test was significant with 71% of genetic variation found within the 17 combined populations and 29% among populations (Table 4). AMOVA for El Salvador populations found a similar result of significant variation among populations, while one of the populations for Honduras found no significant variation among populations. The Fst values were highest and significant primarily between both Sensuntepeque (Cabanas-ES) and San Ignacio (ES) and other populations, as well as between most collections, compared to those from Honduras (Additional File, Supplemental Table 1).

Table 4 Analysis of molecular variation (AMOVA) among seventeen populations from El Salvador and Honduras

Discussion

Aedes aegypti were collected in El Salvador and Honduras, two adjoining countries in Central America. In El Salvador, a high haplotype number and diversity were found, while in contrast Honduras had few haplotypes. The haplotype H1 was most abundant and widespread in both countries. The most widespread haplotype (H1) in both countries may have been present the longest in both countries and perhaps is a remnant of eradication programs from the 1950s and 1960s [10]. Understanding the distribution of genotypes of Ae. aegypti in the region may contribute to management of this species.

Two predominant haplotypes, H1 and H2, were identified from the combined 17 populations. The current study included GenBank accessions obtained from samples in a previous study in El Salvador [15]; the study included a phylogenetic analysis which found that H1 and H2 belonged to two distinct Ae. aegypti mitochondrial lineages [15]. One lineage contains the widespread H1 haplotype, and the second lineage contains the H2 haplotype [2, 15]. Both lineages have been dispersed into many locations worldwide. The two most abundant haplotypes, H1 and H2 (and their associated lineages), could indicate separate invasions at different times. A number of other studies have similarly found two lineages of Ae. aegypti in the Western Hemisphere. In Panama, Eskildsen et al. (2018) recovered two deeply divergent mitochondrial clades, and similarly two genetic lineages were found in Bolivia [36]. Others studying South America similarly suggest the introduction of at least two Ae. aegypti lineages. In Peru, Venezuela, and in the Amazon, two Ae. aegypti mitochondrial lineages were detected [37,38,39]. These are all consistent with our results from El Salvador and Honduras.

In Honduras, nearly all individuals were haplotype 1 (H1). Haplotype 10 (H10), found in Northern Honduras in the department of Cortés, showed the greatest genetic distance from populations in El Salvador and southern Honduras. Geographically, this locality is characterized by its proximity to Puerto Cortés on the Honduran coast and Puerto Barrios on the Guatemalan Atlantic coast. There may be passive transport of eggs through ports into the department of Cortés Honduras, which would incorporate new alleles in this locality, causing an increase in its divergence from other populations. Results of the current study agree with the previous study in Honduras of Ae. aegypti [16], but the sampling region to date in Honduras is limited. Additional collections from other regions of the country, such as the north and eastern region, would provide more information about the genetic diversity of Ae. aegypti.

A high genetic diversity of Ae. aegypti might be expected in both countries, given the international trade between each country and exterior countries, which would provide the means for potential introduction of additional Ae. aegypti genotypes. The widespread distribution of the haplotype H1 through both countries suggests that the haplotype may have been present longer in the region than other haplotypes. One potential explanation for this could be that H1 is a remnant of eradication programs.

Populations with the highest genetic diversity indices were mostly in the east of El Salvador, in La Unión, San Francisco Gotera, San Miguel, and the capital San Salvador. La Unión (ES) has an old port and is characterized by the commercial flow with two neighboring countries, Honduras and Nicaragua, which could facilitate introductions of eggs via the Gulf of Fonseca. San Miguel and San Salvador (capital of ES) also have high commercial interconnectivity with the rest of the country. In addition, ecological factors may contribute to the observed genetic diversity. The private haplotypes observed in the network could have originated from factors such as local radiation, the adaptive processes associated with colonizing new habitats. The eastern region of the country including San Miguel and La Unión has high temperatures which facilitates rapid generation time for this mosquito species. Along with high temperatures, abundant mosquitoes and insecticide control, both anthropogenic and ecological factors could facilitate the production and movement of distinct haplotypes in this region [40, 41].

The Tajima’s D for El Salvador was positive but not significant (1.85), suggesting some genetic differentiation is occurring [42]. In addition, many populations in ES had unique haplotypes. Only Sonsonate (ES) had a significant negative Tajima’s D value. In Honduras, Tajima’s D was −2.0 overall, and values were negative and significant for Villanueva and Choluteca, indicating a recent bottleneck or selective sweep [42]. The AMOVA among all populations combined was significant, indicating gene flow among populations. However, the AMOVA for Honduras found most variation was within populations rather than between them.

This study has some limits and would benefit from additional related research. Mosquito samples could be obtained throughout additional regions of Honduras to better understand the genetic variability of Ae. aegypti through a wider region of the country. Samples from Honduras came mostly from an interior region of the country. Results from the current study suggest that sampling method may impact detection of genotypes. In El Salvador, larval samples were obtained for all populations, and a large number of haplotypes were found. In Honduras, collections originated from ovitraps and less genetic diversity were detected. Haplotype H1 was the most abundant haplotype found in Honduras, whereas El Salvador had an abundance of H1 and H2, along with numerous other haplotypes. A future study in Honduras might include larval collections and could examine if trap type influences the haplotype of Ae. aegypti collected.

Conclusions

It is important to note that one genotype, H1, was widespread and abundant in both Honduras and El Salvador and present at nearly all collection sites. To further investigate whether H1 is a remnant of previous eradication programs, Ae. aegypti sequences from other countries in the Americas could be reexamined to determine whether this particular haplotype is more abundant and widespread in other regions which were under eradication. In El Salvador, H2 was the second most abundant type after H1. It is also possible that these abundant and widespread genotypes, H1 and H2, have persisted because of significant insecticide resistance [43, 44]. This could be investigated and may benefit management. Currently, Wolbachia-infected Ae. aegypti mosquitoes are being implemented in El Salvador and Honduras to reduce populations of Ae. aegypti [45]. The current study provides a foundation to understand the genetic variants of Ae. aegypti present in the region before these mosquito releases occurred. Understanding the mating frequency between the two mitochondrial lineages of Ae. aegypti (H1 and H2) may also be helpful and could contribute to the success of novel mosquito management strategies that rely on cross mating among populations. There remains a constant global challenge to control Ae. aegypti and reduce its populations to prevent vector borne disease.

Availability of Data and Materials

The 90 new sequences generated during the current study are available in GenBank (El Salvador sample accession numbers PP732749-PP732793, Honduras sample numbers PP732794-PP732838). The 82 existing additional accessions (sequences) from GenBank included in the analysis were the following from El Salvador: Usulatán (MK028219-32), Morazán (MK028233-52), Chalatenango (MK028 253-62), San Miguel (MK028263-78), Sonsonate (MK028279-90), and San Salvador (MK028291-300)].

References

  1. Brown JE, Evans BR, Zheng W, Obas V, Barrera-Martinez L, Egizi A, et al. Human impacts have shaped historical and recent evolution in Aedes aegypti, the dengue and yellow fever mosquito. Evolution. 2014;68:514–25.

    Article  CAS  PubMed  Google Scholar 

  2. Gloria-Soria A, Lima A, Lovin DD, Cunningham JM, Severson DW, Powell JR. Origin of a high-latitude population of Aedes aegypti in Washington DC. Am J Trop Med Hyg. 2018;98:445–52.

    Article  PubMed  Google Scholar 

  3. Gloria-Soria A, Ayala D, Bheecarry A, Calderon-Arguedas O, Chadee DD, Chiappero M, et al. Global genetic diversity of Aedes aegypti. Mol Ecol. 2016;25:5377–95.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Bosio CF, Beaty BJ, Black WC IV. Quantitative genetics of vector competence for dengue-2 virus in Aedes aegypti. Am J Trop Med Hyg. 2018;59:965–70.

    Article  Google Scholar 

  5. Dickson LB, Sanchez-Vargas I, Sylla M, Flemming K, Black WC IV. Vector competence in West African Aedes aegypti is flavivirus species and genotype dependent. PLoS Negl Trop Dis. 2014;8:e3153.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Amoa-Bosompem M, Kobayashi D, Itokawa K, et al. Determining vector competence of Aedes aegypti from Ghana in transmitting dengue virus serotypes 1 and 2. Parasite Vector. 2021;14:228.

    Article  CAS  Google Scholar 

  7. Rheinhold JM, Lazzari CR, Lahondere C. Effects of the environmental temperature on Aedes aegypti and Aedes albopictus mosquitoes: a review. Insects. 2018;9:158.

    Article  Google Scholar 

  8. Crawford JE, Clarke DW, Criswell V, et al. Efficient production of male Wolbachia-infected Aedes aegypti mosquitoes enables large-scale suppression of wild populations. Nat Biotechnol. 2020;38:482–92.

    Article  CAS  PubMed  Google Scholar 

  9. Sánchez-Aldana-Sánchez GA, Liedo P, Bond JG, Dor A. Release of sterile Aedes aegypti mosquitoes: chilling effect on mass-reared males survival and escape ability and on irradiated males sexual competitiveness. Sci Rep. 2023;13:3797.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Soper FL. Aedes aegypti and yellow fever. Bull World Health Organ. 1967;36:521–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Gubler DJ, Trent DW. Emergence of epidemic dengue/dengue hemorrhagic fever as a public health problem in the Americas. Infect Agents Dis. 1993;2:383–93.

    CAS  PubMed  Google Scholar 

  12. Salvatella-Agrelo R. Aedes aegypti, Aedes albopictus (Diptera, Culicidae) y su papel como vectores en las Américas. La situación de Uruguay. Rev Med Uruguay. 1996;12:28–36.

    Google Scholar 

  13. Hayes JM, Garcia-Rivera E, Flores-Reyna R, Suarez-Rangel G, Rodriguez-Mata T, Coto-Portillo R, et al. Risk factors for infection during a severe dengue outbreak in El Salvador in 2000. Am J Trop Med Hyg. 2003;69:629–33.

    Article  PubMed  Google Scholar 

  14. Monteiro FA, Shama R, Martins AJ, Gloria-Soria A, Brown JE, Powell JR. Genetic diversity of Brazilian Aedes aegypti: patterns following an eradication program. PLoS Negl Trop Dis. 2014;8:e3167.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Joyce AL, Torres MM, Torres R, Moreno M. Genetic variability of the Aedes aegypti (Diptera: Culicidae) mosquito in El Salvador, vector of dengue, yellow fever, chikungunya and Zika. Parasite Vector. 2018;11:637.

    Article  Google Scholar 

  16. Escobar D, Ortiz B, Urrutia O, Fontecha G. Genetic diversity among four populations of Aedes aegypti (Diptera: Culicidae) from honduras as revealed by mitochondrial DNA cytochrome oxidase I. Pathogens. 2022;11:620.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Hotez PJ, Woc-Colburn L, Bottazzi ME. Neglected tropical diseases in Central America and panama: review of their prevalence, populations at risk and impact on regional development. Int J Parasitol. 2014;44:597–603.

    Article  PubMed  Google Scholar 

  18. Joyce AL, Alvarez FS, Hernandez E. Forest coverage and socioeconomic factors associated with dengue in El Salvador, 2011–2013. Vect Borne Zoon Dis. 2021;21:8.

    Google Scholar 

  19. Ávila-Agüero ML, Camacho-Badilla K, Brea-del-Castillo J, Cerezo L, Duenas L, Luque M, et al. Epidemiología del dengue en centroamérica y república dominicana. Rev Chilena Infectol. 2019;36:698–706.

    Article  PubMed  Google Scholar 

  20. MINSAL 2020. Boletines epidemiológicos semanales https://www.salud.gob.sv/boletines-epidemiologicos-2020/. Accessed 28 Mar 2024.

  21. OPS. Dengue 2022. https://www3.paho.org/data/index.php/es/temas/indicadores-dengue.html. Accessed 28 Mar 2024.

  22. Bennett KL, McMillan OW, Loaiza JR. The genomic signal of local environmental adaptation in Aedes aegypti mosquitoes. Evol Appl. 2021;14:1301–13.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Eskildsen GA, Rovira JR, Smith O, Miller MJ, Bennett KL, McMillan WO, et al. Maternal invasion history of Aedes aegypti and Aedes albopictus into the Isthmus of panama: implications for the control of emergent viral disease agents. PLoS ONE. 2018;13:e0194874.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Apostol BL, Black WC, Reiter P, Miller BR. Use of randomly amplified polymorphic DNA amplified by polymerase chain reaction markers to estimate the number of Aedes aegypti families at oviposition sites in San Juan, Puerto Rico. Am J Trop Med Hyg. 1994;51:89–97.

    Article  CAS  PubMed  Google Scholar 

  25. Burkett-Cadena ND. Mosquitoes of the Southeastern United States. Tuscaloosa: University of Alabama; 2013. p. 202.

    Google Scholar 

  26. Qiagen. DNeasy tissue handbook. Valencia: Qiagen; 2006.

    Google Scholar 

  27. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.

    CAS  PubMed  Google Scholar 

  28. Cywinska A, Hunter FF, Hebert PDN. Identifying canadian mosquito species through DNA barcodes. Med Vet Entomol. 2006;20:413–24.

    Article  CAS  PubMed  Google Scholar 

  29. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38:3022–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Rozas J, et al. DnaSP 6: DNA sequence polymorphism analysis of large datasets. Mol Biol Evol. 2017;2017:3299–302.

    Article  Google Scholar 

  33. Leigh JW, Bryant D. POPART: full feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.

    Article  Google Scholar 

  34. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under linux and windows. Mol Ecol Resour. 2010;10:564–7.

    Article  PubMed  Google Scholar 

  35. Peakall R, Smouse PE. GenAlEx 6: genetic analysis in excel. population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  36. Paupy C, Le Goff G, Brengues C, Guerra M, Revolla J, Barja Simon Z, et al. Genetic structure and phylogeography of Aedes aegypti, the dengue and yellow-fever mosquito vector in Bolivia. Infect Genet Evol. 2012;12:1260–9.

    Article  PubMed  Google Scholar 

  37. Costa da Silva AL, Capurro ML, Bracco JE. Genetic lineages in the yellow fever mosquito Aedes (Stegomyia) aegypti (Diptera: Culicidae) from Peru. Mem Inst Oswaldo Cruz. 2005;100:539–44.

    Article  PubMed  Google Scholar 

  38. Herrera F, Urdaneta L, Rivero J, Zoghbi J, Ruiz J, Carrasquel G, et al. Population genetic structure of the dengue mosquito Aedes aegypti in Venezuela. Mem Inst Oswaldo Cruz. 2006;101:625–33.

    Article  CAS  PubMed  Google Scholar 

  39. Lima RS Jr, Scarpassa VM. Evidencia de dos linajes del vector del dengue Aedes aegypti en la Amazonía brasileña, con base en secuencias del gen ND4 del ADN mitocondrial. Genet Mol Biol. 2009;32:414–22.

    Article  CAS  PubMed  Google Scholar 

  40. Crisci JV. The voice of historical biogeography. J Biogeogr. 2001;28:157–68.

    Article  Google Scholar 

  41. Piñero D, Barahona A, Eguiarte L, Rocha Olivares A, Salas Lizana R. La diversidad genética como instrumento para la conservación y el aprovechamiento de la biodiversidad: estudios en especies mexicanas, en Capital natural de México. Conabio, México: Vol.I: conocimiento actual de la biodiversidad; 2008. p. 437–94.

    Google Scholar 

  42. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Bisset Lazcano JL, Rodríguez MM, San Martin JL, Romero JE, Montoya R. Evaluación de la resistencia a insecticidas de una cepa de Aedes aegypti de El Salvador. Rev Panam Salud Publica. 2009;26:229–34.

    Article  Google Scholar 

  44. Franco DN, Nelly D, Cortes SP, Soledad P. Resistencia y susceptibilidad de Aedes aegypti a insecticidas y themephos en el municipio del distrito Central. Honduras Centro América. Tesis MSc en Epidemiología y Salud: Universidad Nacional Autónoma de Nicaragua, León; 2017.

    Google Scholar 

  45. World Mosquito Program. WMP El Salvador Fact Sheet. 2023. https://www.worldmosquitoprogram.org/en/global-progress/el-salvador. Accessed 26 Mar 2024.

Download references

Acknowledgements

We thank Mackenzie Hoyt, Michaela Hoyt, and Fabian Gonzalez, University of California Merced students, who helped with molecular work. We thank Francisco Samuel Alvarez of for assistance with constructing the maps in ArcGis.

Funding

This project was supported by the University of California Merced, which contributed to the cost of research supplies and undergraduate student research assistants.

Author information

Authors and Affiliations

Authors

Contributions

ALJ, MM, LP, ROC, and DG designed the study. LP and ROC provided samples. LP and ALJ performed the molecular work. ALJ MM LP analyzed the data; AJ, MM, LP, ROC, and DE discussed the results, wrote, edited and approved the final manuscript.

Corresponding author

Correspondence to A. L. Joyce.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

13071_2024_6312_MOESM1_ESM.docx

Additional File1 Table S1. Pairwise FST distance among populations. Populations from El Salvador: Ataco (AT), Sensuntepeque (SEN), Jiquilisco (JI), La Unión (LU), San Francisco Gotera (SFG), Santa Ana (SA), San Ignacio (SI), San Miguel (SM), Sonsonate (SON), San Salvador (SS), San Vicente (SV). Populations from Honduras: Tegucigalpa (TEG), Choluteca (CHO), Cofradia (COF), Villanueva (VIL), Danli (DAN), Los Angeles (LA). Collection locations of all populations are detailed in Table 1.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Joyce, A.L., Moreno, M., Palomo, L. et al. Genetic variability of Aedes aegypti (Diptera: Culicidae) in El Salvador and Honduras: presence of a widespread haplotype and implications for mosquito control. Parasites Vectors 17, 229 (2024). https://doi.org/10.1186/s13071-024-06312-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-024-06312-7

Keywords