- Open Access
Aedes aegypti has spatially structured and seasonally stable populations in Yogyakarta, Indonesia
Parasites & Vectorsvolume 8, Article number: 610 (2015)
Dengue fever, the most prevalent global arboviral disease, represents an important public health problem in Indonesia. Control of dengue relies on the control of its main vector, the mosquito Aedes aegypti, yet nothing is known about the population history and genetic structure of this insect in Indonesia. Our aim was to assess the spatio-temporal population genetic structure of Ae. aegypti in Yogyakarta, a densely populated region on Java with common dengue outbreaks.
We used multiple marker systems (microsatellites, nuclear and mitochondrial genome-wide single nucleotide polymorphisms generated via Restriction-site Associated DNA sequencing) to analyze 979 Ae. aegypti individuals collected from the Yogyakarta city and the surrounding hamlets during the wet season in 2011 and the following dry season in 2012. We employed individual- and group-based approaches for inferring genetic structure.
We found that Ae. aegypti in Yogyakarta has spatially structured and seasonally stable populations. The spatial structuring was significant for the nuclear and mitochondrial markers, while the temporal structuring was non-significant. Nuclear markers identified three main genetic clusters, showing that hamlets have greater genetic isolation from each other and from the inner city sites. However, one hamlet experienced unrestricted mosquito interbreeding with the inner city, forming a single genetic cluster. Genetic distance was poorly correlated with the spatial distance among mosquito samples, suggesting stronger influence of human-assisted gene flow than active mosquito movement on spatial genetic structure. A star-shaped mitochondrial haplotype network and a significant R 2 test statistic (R 2 = 0.0187, P = 0.001) support the hypothesis that Ae. aegypti in Yogyakarta originated from a small or homogeneous source and has undergone a relatively recent demographic expansion.
We report the first insights into the spatio-temporal genetic structure and the underlying processes in the dengue fever mosquito from Yogyakarta, Indonesia. Our results provide valuable information on the effectiveness of local control measures as well as guidelines for the implementation of novel biocontrol strategies such as release of Wolbachia-infected mosquitoes.
Dengue fever is the most prevalent arboviral disease in the world and Indonesia boasts one of the highest numbers of dengue infections amongst the most endemic nations . Moreover, incidence of a severe disease form known as the dengue hemorrhagic fever has dramatically increased in Indonesia over the past 45 years . A surge in dengue outbreaks is attributed to a combination of factors such as rapid human population growth, migration from rural to urban areas, inadequate basic urban infrastructure and public health measures , all of which favor viral transmission by the main vector Aedes aegypti. This mosquito has a close association with human habitation, preferentially exploiting domestic environments for its feeding and breeding requirements [4–6].
Given that an effective vaccines or treatments against dengue viruses are not yet available, the control of dengue transmission is mostly dependent on the control of the vector populations . Insecticide applications remain the predominant strategy to control the dengue fever mosquito, but their extensive and inappropriate use has led to a development of insecticide resistance in Ae. aegypti populations around the world, including Indonesia [7, 8]. In the last 15 years, a significant advancement has been made in developing and implementing alternative strategies to overcome the limitation of traditional measures and to achieve a sustainable control of vector populations . One such strategy involves the use of the intracellular bacterium Wolbachia pipientis that has been artificially transferred into Ae. aegypti [10, 11]. Wolbachia infection induces a number of pathogen-resistant phenotypes that reduce the mosquito’s vectorial competency . This strategy is based on the releases of Wolbachia-infected Ae. aegypti into the field to replace the vector competent populations with mosquitoes that have a significantly diminished capacity for pathogen transmission [13–15].
Following the successful Wolbachia releases in northern Australia and Vietnam [15, 16], the implementation of this strategy is planned for the region of Yogyakarta (http://www.eliminatedengue.com/indonesia). Given that the establishment and spread of Wolbachia infection is influenced by the dynamics of local populations as well as migration among them [17, 18], it is essential to obtain knowledge of the population genetic structure and gene flow in Ae. aegypti throughout Yogyakarta.
Apart from information on the strong genome-wide differentiation between a small sample of Ae. aegypti from Indonesia and samples from Vietnam, Australia and Brazil , there is no knowledge on the population genetic structure in the dengue fever mosquito from Indonesia. Studies from other countries in southeast Asia have generally shown significant population structuring for a range of spatial scales (5 to >2000 km) . However, the level of genetic structuring was highly variable among areas with different levels of urbanization (rural vs. urban) or between seasons (rainy vs. dry) (reviewed in ).
Here, our aim was to assess the spatio-temporal population genetic structure of Ae. aegypti in the city of Yogyakarta and the surrounding hamlets. This large and densely populated urban environment on the island of Java with an approximate population of 3.4 million people has common dengue outbreaks . For example, an average annual rate of 16.8 dengue cases per 10,000 inhabitants has been reported in Yogyakarta city from 2005 to 2007 (Yogyakarta City Health Office, 2008). We used several types of genetic marker (microsatellites, nuclear and mitochondrial genome-wide single nucleotide polymorphisms) to genotype individuals of Ae. aegypti collected during one wet and one dry season. We considered the following questions: 1) what is the level of structuring at the nuclear and mitochondrial genomes among sites in Yogyakarta? 2) Is there a seasonal variation in genetic structuring? 3) Is the level of gene flow different among urban and suburban areas (i.e. Yogyakarta city and hamlets)?
Thirteen collection sites within three districts were established in the Special Region of Yogyakarta for this study (Fig. 1, Table 1). Samples were collected from the districts of Bantul and Sleman between October to December of 2011 as a representation of the wet season and in July of 2012 as a representation of the dry season. The average precipitation in Yogyakarta in July is 40 mm and increases to 90–340 mm in October to December. The sites are separated by areas of agricultural land with many large roads/highways and rivers/canals which are likely to impede mosquito dispersal. Dry season samples also included a collection made in November 2012 from the City of Yogyakarta (Fig. 1, sites 11–13), allowing for the study of mosquito population structure between continuously populated sites without the agricultural lands. The area covered by districts Sleman, Yogyakarta City and Bantul is approximately 1,115.45 km2 with each collection site spanning between 0.1 km2 and 1.0 km2.
Within each collection site, Aedes aegypti larvae and adults were collected from ovitraps and backpack aspirators (BPA) positioned inside or very close to the houses (i.e. front- or backyards). Eggs collected from ovitraps were reared to the third larval instar to facilitate identification of species. A single larva from each ovitrap was used for analysis to minimise sampling of family groups [19, 23, 24], while multiple adults from each house were sampled and analysed.
Individuals were genotyped at nine microsatellite markers previously described (AC1, AG5, 88AT1, BbA10, BbH08, M201, 12ACG1, 470AG1 and 69TGA1 [25–27]). The loci were separated into three multiplex PCR assays as described in . Microsatellite fragment separation analysis was conducted using the Applied Biosystems 96-capillary 3730xl DNA Analyser and GeneScan™ 500 LIZ® size standard and scored using GeneMarker® Software Version 1.91 . The final microsatellite dataset included: 152 adults and 351 larvae sampled across hamlets in wet season, and 320 larvae sampled from hamlets and inner city during dry season (Table 1).
Genome-wide single nucleotide polymorphism (SNP) genotyping was done using a customized double-digest RADseq method  as described in . 100 ng of genomic DNA from each larva was digested with 10 units of restriction enzymes NlaIII and MluCI (New England Biolabs, Beverly MA, USA). Illumina adapters with customized barcode sequences (100 pM of P1, 300 pM of P2) were ligated to the genomic fragments using 100 units of T4 ligase at 16 °C overnight (New England Biolabs, Beverly, MA, USA). Size selection of 300–450 bp long fragments was completed with the Pippin Prep protocol (Sage Sciences, Beverly, MA, USA). For the final library enrichment, five repeated PCR reactions with 1 nM standard Illumina primers included the following cycling conditions: 98 °C for 30 s, 12 cycles of 98 °C for 10 s, 60 °C for 30 s, 72 °C for 90 s, and the final elongation at 72 °C for 5 min. The independent reactions were pooled and purified to create a final library. Four libraries were sequenced in four lanes of Illumina HiSeq2000 using the 100 bp paired-end chemistry (SRA Accession Numbers: SRP040064, SRX1425749).
Raw Illumina reads were filtered and trimmed based on a minimum phred score of 20. We used the short read aligner program Bowtie  to align the filtered reads to the Ae. aegypti nuclear and mitochondrial reference genomes [32, 33]. Because we expected our populations to show substantial genome divergence from the reference strain (i.e. Liverpool strain that originated from a West African population ), we allowed for the maximum number of mismatches in the alignment seed. The search algorithm was set to try as hard as possible to find valid alignments when they existed [Bowtie option—tryhard], and only uniquely aligned reads were then analysed with the Stacks pipeline . Default program parameters were used for variant and likelihood-based genotype calling at RAD stacks that had a depth of at least five reads. The final RADseq dataset included 46 larvae collected across three hamlets in the wet season and 117 larvae collected from five hamlets and three inner city locations in the dry season (Table 1, Additional file 1: Table S1). 27 larvae from the wet season had more than 30 % missing data for nuclear SNPs and were excluded from downstream analyses for this marker set.
Analyses of genetic diversity and structure
Basic population genetic analyses such as observed and expected heterozygosity and tests for departure from Hardy-Weinberg equilibrium were carried out using GenAlEx v6.5  and Fstat 184.108.40.206 . We used the program FreeNA to test for the presence of null alleles at microsatellite loci that can bias the estimates of within-population diversity and population structure .
To estimate the level of spatio-temporal genetic structuring, we performed the hierarchical analysis of molecular variance (AMOVA) and post-hoc pairwise F ST comparisons, treating each sampling location as an independent sample nested within season. Statistical significance was determined from 999 permutations using the program GenAlEx v6.5 .
Genetic structure was also inferred using the individual-based analysis implemented in the program STRUCTURE v2.3.4 [38, 39], that uses multi-locus genotypes to determine the number of genetic groups (K) and assign individuals to those groups. The parameter set included: admixture model with the information on the sampling location and correlated allele frequencies among populations. The number of burn-ins was 250,000 for the microsatellite datasets and 50,000 for the SNP datasets, followed by as many MCMC steps. Estimations were made for K = 1 to N, where N represents the number of sampling sites in a given season. The most likely number of clusters was determined using the maximum likelihood guideline by Pritchard  and the deltaK method by Evanno et al.  implemented with the program Structure Harvester . We also used a multivariate method called discriminant analysis of principal components (DAPC)  implemented within the R package adegenet , to infer the number of genetic groups and the individual assignment to those groups.
Presence of isolation by distance (IBD)  between sampling groups was estimated using the matrices of pair-wise linearized F ST values and the natural log of pairwise geographic distance between sites. IBD among individuals was tested using the matrix of unweighted squared genetic distances between individuals  and the natural log of pairwise geographic distance between individuals. Significance of the correlation between the matrices of genetic and geographic distances was determined using a Mantel test with 999 permutations in GenAlEx v6.5 . The spatial autocorrelation analysis  with 999 permutations was executed in the same program. An estimate of effective population size, Ne, was obtained using the two-sample temporal analysis method  as applied in NeEstimator 2.01  for both microsatellite and SNP data at sites 6, 7 and 8. Generations sampled were set at 0 and 10 based on the length of time between the two field collections and estimated 15 generations per year in Yogyakarta. Waples’ method  uses moments-based F statistics and, in this case, the estimator Fs  was applied as it has been developed to provide unbiased estimates of Ne for multiple types of data sets including large panels of SNPs.
The analyses of the number of mitochondrial haplotypes (Nh), nucleotide diversity (π), haplotype diversity (Hd), as well as the spatio-temporal AMOVA were done in the R package pegas . We used R 2 test  for detecting population growth, determining the statistical significance with 1000 simulations of populations under the drift-mutation equilibrium in pegas . The R 2 test has greater power than a classical mismatch distribution test for smaller sample sizes . Haplotype network analysis with the statistical parsimony method was completed with the R package TempNet .
Allelic richness over nine microsatellite loci within each site ranged between 3.81 to 5.31 in wet season and between 5.14 and 5.81 in dry season (Table 1). Site 2 exhibited the highest allelic richness in both seasons. In general, observed heterozygosity at each sampling location was lower than expected heterozygosity (Table 1) and markers did not exhibit significant linkage disequilibrium.
After removing nuclear RAD tags with more than 30 % of missing data, the final dataset included 3,178 polymorphic tags distributed across 830 scaffolds on all three chromosomes of Ae. aegypti. Overall nucleotide diversity in these RAD tags was 0.0023. Because multiple SNPs within the same RAD tag are expected to be in nearly complete linkage disequilibrium, we retained only one SNP per tag. Expected heterozygosity across 3,178 SNPs was similar for all sites, ranging between 0.205 and 0.230, and was somewhat higher than the observed heterozygosity (F IS from 0 to 0.157, Table 1).
Nuclear genetic structure
For both marker systems, hierarchical analysis of molecular variance (AMOVA) indicated strong spatial genetic structure that was temporally stable (Table 2). Three percent of total variation at microsatellite loci was attributed to differences in allele frequencies among sites within a season, while no significant variation could be attributed to changes across seasons. Spatial difference in frequencies of nuclear genome-wide SNPs contributed to nearly seven percent of total variation (Table 2). Pairwise F ST comparisons at SNPs revealed that mosquitoes from sites 2 and 6 showed the highest differentiation at within the region (Table 3). Also, adult and larval samples collected within the same hamlet did not show a significant difference in allele frequencies, except at site 6 (Additional file 1: Table S1). Genome-wide SNPs revealed much lower differentiation among city sites than among hamlets, while microsatellite markers did not reveal this pattern (Table 3).
STRUCTURE analyses of adult mosquitoes suggested four genetic groups based on the highest likelihood of data, and two groups based on the deltaK method (Additional file 1: Table S2). For the larval samples, both methods indicated that the number of genetic groups was three in the wet season and two in the following dry season (Additional file 1: Table S2).
Genome-wide SNPs showed much greater power in distinguishing Ae. aegypti individuals from different sites in Yogyakarta than microsatellite markers (Fig. 2) and DAPC without prior information on the sampling location (‘blind DAPC’) revealed three genetic clusters in the SNP dataset. Again, sites 2 and 6 showed greater separation from other sites in the region, while site 8 and inner city cites could not be distinguished. Moreover, temporal stability of genetic structure was demonstrated by the overlap of temporal samples in the PC space (Fig. 2).
Presence of IBD was highly supported between individual adult mosquitoes (microsatellites: Mantel r = 0.106, P < 0.01), and individual larvae in the season (microsatellites: Mantel r = 0.106, P < 0.01; SNPs: Mantel r = 0.247, P < 0.01). Spatial autocorrelation analysis indicated that greater relatedness between pairs of larvae can be found at distances of 8–9 km, even though it is low for larvae from different sites (i.e. autocorrelation coefficient is significantly negative at distances 3–6 km) (Fig. 3).
Site 6 had the lowest estimated Ne out of the three sites tested, using either microsatellite or SNP data (Table 4). Tighter estimates of Ne for the SNP data are in line with the observation that current methods tend to produce intervals that are too narrow when a large number of markers is used . 95 % confidence intervals for Ne generated from SNP data were 128–144, 265–298 and 689–774, for Sites 6, 7 and 8 respectively (Table 4). For microsatellite data, estimates from the same sites were 212–564, infinite and 449–1303, respectively. Infinite value for the unbiased Ne estimator means that temporal genetic variation can be explained by the sampling error alone, and not by the effects of a finite number of parents .
Mitochondrial diversity and structure
We found 16 RAD tags in eight mitochondrial genes (ND2, COXI, COX3, ATP6, ND4-ND6, CYTB) that were present in more than 80 % of individuals. We first checked for the presence of premature stop codons or heterozygous loci that would suggest paralogous/nuclear copies of mitochondrial sequences [53, 54] and then we concatenated all tags into sequences that were 990 bp long. Additionally, we confirmed that in 14 individuals ddRAD tags had sequences identical to the longer amplicons from the corresponding mitochondrial regions . Nevertheless, it is possible to have some spurious signals originating from the nuclear copies despite the control steps we employed in our study. There were 22 unique haplotypes in the full dataset (combined across sites and seasons), with haplotype 2 being shared among all eight sites, and haplotype 9 among five sites (Fig. 4). Haplotype diversity ranged between 0.250 and 0.972 and nucleotide diversity ranged between 0.0004 and 0.0044 (Table 1, Fig. 4). AMOVA revealed significant spatial structure after accounting for small seasonal changes of the mitochondrial haplotype frequencies (Table 2). The haplotype network for the entire dataset was star-shaped, with one predominant haplotype (haplotype 2) and pair-wise haplotype differences ranging between one and eight base pairs (Fig. 5). The star-shaped haplotype network suggests exponential population growth , and this hypothesis was supported with the R 2 test for recent population expansion (R 2 = 0.0187, P = 0.001).
Multiple marker systems and temporal sampling of Aedes aegypti across Yogyakarta revealed that the dengue fever mosquito has spatially structured and seasonally stable populations in this region of Indonesia. Nuclear markers (genome-wide SNPs and microsatellites) showed that mosquitoes from some hamlets exhibit greater genetic isolation, regardless of their geographic position (e.g. sites 2 and 6). Also, mosquitoes from inner city sites and one hamlet (site 8) seem to interbreed freely, constituting a single genetic cluster. Mitochondrial variation suggests relatively recent population expansion of Ae. aegypti in this region.
Aedes aegypti is thought to have invaded southeast Asia from the New World as recently as the 20th century . Mitochondrial data available in our study do not allow for identification of the origin of Ae. aegypti in Yogyakarta, but the presence of a single haplogroup supports the hypothesis of a small or genetically homogeneous source that underwent population expansion. This is in sharp contrast to Ae. aegypti from Rio de Janeiro, Brazil, where several highly divergent mitochondrial lineages have been found in one continuously distributed large population . Our results are consistent with previous reports of much lower genetic diversity in Ae. aegypti populations from southeast Asia when compared with those in the New World. For example, Bosio et al.  found that the dengue fever mosquito in Thailand harbors mitochondrial diversity that is 2–3 times lower than that in Mexico.
Analyses of nuclear markers suggest stronger influence of human-assisted gene flow than active dispersal in shaping the spatial pattern of genetic structure in Ae. aegypti from Yogyakarta. Genetic clustering did not clearly correspond to the spatial distance among sampling locations, resulting in a very weak isolation-by-distance. Generally, hamlets showed greater genetic isolation from each other and from the inner city sites. However, one hamlet (site 8) exhibited unrestricted interbreeding with mosquitoes from the inner city, forming a single genetic cluster. The spatial genetic structure of Ae. aegypti in Indonesia parallels the examples from Thailand, where mosquitoes show overall weaker differentiation within urban areas such as Bangkok or Hat Yai than in suburban or rural areas [60, 61]. Exceptions to these common patterns, like the case of high genetic connectivity between some hamlets and the urban center in Yogyakarta, could help identify specific human activities/means of transportation that facilitate extensive long distance gene flow in this mosquito.
A recent entomological study in Yogyakarta found that numbers of adult mosquitoes tended to fluctuate by season in sites 2 and 6, but were stable throughout a year in sites 7 and 8 . We did not detect any significant seasonal changes in genetic structure at site 6, but the estimated effective population size was several times lower than in site 8. Some of the local variation in mosquito numbers may be attributed to differences in rainfall patterns that varied across hamlets . Other factors could include locally varying mosquito control measures, such as removal of breeding habitats and application of insecticides which has been widespread in Indonesia since the 1970s . However, our data suggest that local control measures have had little overall success in reducing mosquito numbers to levels that would cause significant genetic bottlenecks, as seen from a temporally stable genetic structure in populations of Ae. aegypti across Yogyakarta. Diminishing efficacy of insecticide-based mosquito control is seen globally due to rapidly evolving resistance to all classes of insecticides . Recent results from bioassays and genetic screening of putative resistance alleles in Ae. aegypti demonstrated that pyrethroid insecticides are likely to be losing efficacy in Yogyakarta, urging the need for the employment of resistance management tactics . Our insights into the genetic connectivity of mosquito populations should assist in predicting the spread of the resistance alleles in this region.
Knowledge of the spatio-temporal genetic diversity and structuring in Ae. aegypti from Yogyakarta obtained through our study is also useful for implementation of Wolbachia–based strategies in this region. Replacement of local mosquito populations with Wolbachia¬infected mosquitoes could be more challenging in hamlet 8 and within the city sites due to high migration of uninfected mosquitoes into the release area. This could require a process of “adaptive management”, where an increased rate of releases is undertaken in areas with slower Wolbachia invasion .
Our study provides the first insight into the genetic diversity and structure of the dengue fever mosquito, Aedes aegypti, in Indonesia where dengue is endemic. Seasonal stability of population genetic structure and diversity suggests limited effectiveness of traditional control measures to reduce numbers of breeding mosquitoes. Finally, knowledge of the genetic connectivity/isolation of mosquitoes in some parts of Yogyakarta will be useful for optimizing the strategy that involves the releases of Wolbachia-infected mosquitoes to suppress dengue transmission.
World Health Organization. Dengue and severe dengue. WHO Fact Sheet 2012:1–4.
Karyanti MR, Uiterwaal CS, Kusriastuti R, Hadinegoro SR, Rovers MM, Heesterbeek H, et al. The changing incidence of Dengue Haemorrhagic Fever in Indonesia: a 45-year registry-based analysis. BMC Infect Dis. 2014;14:412.
Guzman MG, Halstead SB, Artsob H, Buchy P, Farrar J, Gubler DJ, et al. Dengue: a continuing global threat. Nat Rev Microbiol. 2010;8(12 Suppl):S7–16.
Harrington LC, Edman JD, Scott TW. Why do female Aedes aegypti (Diptera: Culicidae) feed preferentially and frequently on human blood? J Med Entomol. 2001;38:411–22.
Scott TW, Naksathit A, Day JF, Kittayapong P, Edman JD. A fitness advantage for Aedes aegypti and the viruses it transmits when females feed only on human blood. Am J Trop Med Hyg. 1997;57:235–9.
Service MW. Importance of ecology in Aedes aegypti control. Southeast Asian J Trop Med Public Health. 1992;23(4):681–90.
Mulyatno KC, Yamanaka A, Ngadino, Konishi E. Resistance of Aedes aegypti (L.) larvae to temephos in Surabaya, Indonesia. Southeast Asian J Trop Med Public Health. 2012;43:29–33.
Ahmad I, Astari S, Tan M. Resistance of Aedes aegypti (Diptera: Culicidae) in 2006 to pyrethroid insecticides in Indonesia and its association with oxidase and esterase levels. Pakistan J Biol Sci. 2007;10:3688–92.
McGraw EA, O’Neill SL. Beyond insecticides: new thinking on an ancient problem. Nat Rev Microbiol. 2013;11:181–93.
McMeniman CJ, Lane RV, Cass BN, Fong AWC, Sidhu M, Wang Y-F, et al. Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science. 2009;323:141–4.
Xi Z, Khoo CCH, Dobson SL. Wolbachia establishment and invasion in an Aedes aegypti laboratory population. Science. 2005;310:326–8.
Moreira LA, Iturbe-Ormaetxe I, Jeffery JA, Lu G, Pyke AT, Hedges LM, et al. A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya, and Plasmodium. Cell. 2009;139:1268–78.
Walker T, Johnson PH, Moreira LA, Iturbe-Ormaetxe I, Frentiu FD, McMeniman CJ, et al. The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations. Nature. 2011;476:450–3.
Frentiu FD, Zakir T, Walker T, Popovici J, Pyke AT, van den Hurk A, et al. Limited Dengue Virus Replication in Field-Collected Aedes aegypti Mosquitoes Infected with Wolbachia. PLoS Negl Trop Dis. 2014;8:1–10.
Hoffmann AA, Montgomery BL, Popovici J, Iturbe-Ormaetxe I, Johnson PH, Muzzi F, et al. Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature. 2011;476:454–7.
Hoffmann AA, Iturbe-Ormaetxe I, Callahan AG, Phillips BL, Billington K, Axford JK, et al. Stability of the wMel Wolbachia Infection following invasion into Aedes aegypti populations. PLoS Negl Trop Dis. 2014;8, e3115.
Hancock PA, Sinkins SP, Godfray HCJ. Population dynamic models of the spread of Wolbachia. Am Nat. 2011;177:323–33.
Barton NH, Turelli M. Spatial waves of advance with bistable dynamics: cytoplasmic and genetic analogues of allee effects. Am Nat. 2011;178:E48–75.
Rašić G, Filipović I, Weeks AR, Hoffmann AA. Genome-wide SNPs lead to strong signals of geographic structure and relatedness patterns in the major arbovirus vector, Aedes aegypti. BMC Genomics. 2014;15:275.
Hlaing T, Tun-Lin W, Somboon P, Socheat D, Setha T, Min S, et al. Spatial genetic structure of Aedes aegypti mosquitoes in mainland Southeast Asia. Evol Appl. 2010;3:319–39.
Urdaneta-Marquez L, Failloux A-B. Population genetic structure of Aedes aegypti, the principal vector of dengue viruses. Infect Genet Evol. 2011;11:253–61.
Tana S, Umniyati S, Petzold M, Kroeger A, Sommerfeld J. Building and analyzing an innovative community-centered dengue-ecosystem management intervention in Yogyakarta, Indonesia. Pathog Glob Health. 2012;106:469–78.
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.
Hoffmann AA, Goundar AA, Long SA, Johnson PH, Ritchie SA. Invasion of Wolbachia at the residential block level is associated with local abundance of Stegomyia aegypti, yellow fever mosquito, populations and property attributes. Med Vet Entomol. 2014;28 Suppl 1:90–7.
Chambers EW, Meece JK, McGowan JA, Lovin DD, Hemme RR, Chadee DD, et al. Microsatellite isolation and linkage group identification in the yellow fever mosquito Aedes aegypti. J Hered. 2007;98:202–10.
Slotman MA, Kelly NB, Harrington LC, Kitthawee S, Jones JW, Scott TW, et al. Polymorphic microsatellite markers for studies of Aedes aegypti (Diptera: Culicidae), the vector of dengue and yellow fever. Mol Ecol Notes. 2006;7:168–71.
Lovin DD, Washington KO, Hemme RR, Mori A, Epstein SR, Harker BW, et al. Genome-based polymorphic microsatellite development and validation in the mosquito Aedes aegypti and application to population genetics in Haiti. BMC Genomics. 2009;9:1–9.
Rašić G, Endersby NM, Williams C, Hoffmann AA. Using Wolbachia-based release for suppression of Aedes mosquitoes: Insights from genetic data and population simulations. Ecol Appl. 2014;24:1226–34.
Holland MM, Parson W. GeneMarker® HID: A Reliable Software Tool for the Analysis of Forensic STR Data. J Forensic Sci. 2011;56:29–35.
Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One. 2012;7, e37135.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, et al. Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007;316:1718–23.
Behura SK, Lobo NF, Haas B, DeBruyn B, Lovin DD, Shumway MF, et al. Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial dna fragments inserted in the nuclear genomes. Insect Biochem Mol Biol. 2011;41:770–7.
Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA. Stacks: An analysis tool set for population genomics. Mol Ecol. 2013;22:3124–40.
Peakall R, Smouse PE. GenALEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.
Goudet J. FSTAT: a computer program to calculate F-Statistics. J Hered. 2013;104:586–90.
Chapuis MP, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24:621–31.
Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9:1322–32.
Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87.
Pritchard JK. Documentation for structure software: Version 2.3. In Pract. 2010;6:321–6.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol Ecol. 2005;14:2611–20.
Earl DA, vonHoldt BM. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components : a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.
Jombart T, Ahmed I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011;27:3070–1.
Wright S. Isolation by Distance. Genetics. 1943;28:114–38.
Smouse PE, Peakall R. Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure. Heredity (Edinb). 1999;82(Pt5):561–73.
Waples RS. A generalized approach for estimating effective population size from temporal changes in allele frequency. Genetics. 1989;121:379–91.
Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. NeEstimator v2: Re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol Ecol Resour. 2014;14:209–14.
Jorde PE, Ryman N. Unbiased estimator for genetic drift and effective population size. Genetics. 2007;177:927–35.
Paradis E. Pegas: An R package for population genetics with an integrated-modular approach. Bioinformatics. 2010;26(3):419–20.
Ramos-Onsins SE, Rozas J. Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002;19:2092–100.
Prost S, Anderson CNK. TempNet: A method to display statistical parsimony networks for heterochronous DNA sequence data. Methods Ecol Evol. 2011;2:663–7.
Black Iv WC, Bernhardt SA. Abundant nuclear copies of mitochondrial origin (NUMTs) in the Aedes aegypti genome. Insect Mol Biol. 2009;18:705–13.
Hlaing T, Tun-Lin W, Somboon P, Socheat D, Setha T, Min S, et al. Mitochondrial pseudogenes in the nuclear genome of Aedes aegypti mosquitoes: implications for past and future population genetic studies. BMC Genet. 2009;10:11.
Yeap HL, Rašić G, Endersby-Harshman NM, Lee SF, Arguni E, Le Nguyen H, et al. Mitochondrial DNA variants help monitor the dynamics of Wolbachia invasion into host populations. Heredity (Edinb). 2015. doi:10.1038/hdy.2015.97.
Slatkin M, Hudson RR, Marjoram P, Donnelly P. Pairwise comparisons of mitochondrial DNA sequences in subdivided populations. Genetics. 1994;136:673–83.
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 (NY). 2013;68:514–25.
Rašić G, Schama R, Powell R, Maciel-de Freitas R, Endersby-Harshman NM, Filipović I, et al. Contrasting genetic structure between mitochondrial and nuclear markers in the dengue fever mosquito from Rio de Janeiro: implications for vector control. Evol Appl. 2015;8:901–15.
Bosio CF, Harrington LC, Jones JW, Sithiprasasna R, Norris DE, Scott TW. Genetic structure of Aedes aegypti populations in Thailand using mitochondrial DNA. Am J Trop Med Hyg. 2005;72:434–42.
Olanratmanee P, Kittayapong P, Chansang C, Hoffmann AA, Weeks AR, Endersby NM. Population genetic structure of Aedes (Stegomyia) aegypti (L.) at a micro-spatial scale in Thailand: implications for a dengue suppression strategy. PLoS Negl Trop Dis. 2013;7:e1913.
Sukonthabhirom S, Rongnoparut P, Saengtharatip S, Jirakanjanakit N, Chareonviriyaphap T. Genetic structure and gene flow among Aedes aegypti(Diptera: Culicidae) populations from central Thailand. J Med Entomol. 2005;42:604–9.
Tantowijoyo W, Arguni E, Johnson P, Budiwati N, Nurhayati PI, Fitriana I, et al. Spatial and temporal variation in Aedes aegypti and Aedes albopictus numbers in the Yogyakarta area of Java, Indonesia, with implications for Wolbachia releases. J Med Entomol 2015 (in press).
Kusriastuti R, Sutomo S. Evolution of dengue prevention and control program in Indonesia. Dengue Bull. 2005;29:1–7.
Liu N. Insecticide Resistance in Mosquitoes: Impact, Mechanisms, and Research Directions. Annu Rev Entomol. 2015;60:537–59.
Rochmijati Wuliandari J, Lee SF, White VL, Tantowijoyo W, Hoffmann AA, Endersby-Harshman NM. Association between Three Mutations, F1565C, V1023G and S996P, in the Voltage-Sensitive Sodium Channel Gene and Knockdown Resistance in Aedes aegypti from Yogyakarta, Indonesia. Insects. 2015;6(3):658–85.
This work is funded by the National Health and Medical Research Council, the Tahija Foundation and the Foundation for the National Institutes of Health through the Grand Challenges in Global Health Initiative of the Bill and Melinda Gates Foundation. Computational resources were provided by the NeCTAR Research Cloud at the University of Melbourne. The funders had no role in the study design, data collection, analyses and interpretation, or the manuscript writing.
The authors declare no competing interests.
GR processed and analysed the RADseq and microsatellite data, wrote the manuscript. NEH analysed the microsatellite data and wrote the manuscript. WT designed the field work and collected samples. AG and VW performed the microsatellite laboratory work and part of the data analysis. QY performed the RADseq laboratory work and data pre-processing. IF performed the RADseq data processing. PJ designed the field work and collected samples. AAH conceived the study and drafted the manuscript. EA designed and lead the field work and collections. All authors have read and approved the final manuscript.
Gordana Rašić and Nancy Endersby-Harshman are joint first authors. Ary A. Hoffmann and Eggi Arguni are joint senior authors.
Pairwise F ST comparisons between Aedes aegypti adults and larvae. Individuals collected in wet season were scored at nine microsatellite loci. Table S2. STRUCTURE Harvester results for adult and larval mosquitoes that were scored at nine microsatellite loci. (DOCX 19 kb)