PLoS Neglected Tropical Diseases
Public Library of Science
image
Fine-scale population genetic structure of dengue mosquito vector, Aedes aegypti, in Metropolitan Manila, Philippines
Volume: 14, Issue: 5
DOI 10.1371/journal.pntd.0008279
  • PDF   
  • XML   
  •       
Abstract

Aedes aegypti is an efficient vector of dengue due to its highly adaptive nature to the urban environment. Although it is observed to have a short dispersal (active) capability, it has been shown to be capable of traveling long distances (passive) via human-mediated transportation. This duality may expand the distribution of the mosquito vector in urbanized areas. In this study, we examined the population genetic structure of Ae. aegypti in a highly urbanized and dengue-endemic region of the Philippines, Metropolitan Manila. Our findings indicated the dual dispersal nature of Ae. aegypti. The use of microsatellites as genetic markers also allowed us to describe the potential long-distance dispersal patterns, possibly through human-aided land transportation via the existing road networks of Metropolitan Manila.

Carvajal, Ogishi, Yaegeshi, Hernandez, Viacrusis, Ho, Amalin, Watanabe, and Kittayapong: Fine-scale population genetic structure of dengue mosquito vector, Aedes aegypti, in Metropolitan Manila, Philippines

Introduction

Dengue disease is the most prevalent mosquito-borne viral infection in tropical and subtropical countries [1] with approximately 2.5 billion people at risk of contracting the disease worldwide [2]. The Philippines considers dengue as a major public health concern and has been ranked fourth in Southeast Asia with the highest number of dengue cases reported annually [3]. In particular, the National Capital Region (NCR), also known as Metropolitan Manila, is one of the main areas affected by the disease and an average of 16% of the total number of cases reported in the Philippines from 2009 to 2014 [4]. Interestingly, the first recorded dengue epidemic in Southeast Asia occurred in this region in 1954 [5]. Although a dengue vaccine is available [6], the World Health Organization [2] still recommends disease prevention and control towards the mosquito vector.

The transmission of this arboviral disease is through the bite of either the Aedes aegypti or Ae. albopictus mosquito vector. The former is identified as the primary vector of the disease, and as such, its biology has been extensively studied including its reproduction and development [79], survivability [1011], dispersal [1214], and ecology[1516]. In many population genetic studies of Ae. aegypti , microsatellites have been used as a molecular marker of choice for differentiating macro- [17,19] and micro-geographic [20,21,22] scale mosquito populations due to its high polymorphism, co-dominance, and broad genome distribution [23]. Molecular genotyping of the mosquito vector using these markers has shown additional insights towards the microevolution of the mosquito vector [1719], revealing the population’s gene flow pattern which can be interpreted as the mosquito vector’s dispersal ability [20,24].

Urbanization has been a key factor in the widespread occurrence of dengue disease worldwide and has contributed to the dispersal of the mosquito vector [16]. Ae. aegypti has been documented to thrive in highly urbanized areas, creating the ideal condition for a dengue epidemic [25,26]. As such, investigating the impact of urbanization on Ae. aegypti ’s population genetic structure has sparked research interest, particularly regarding its micro-evolutionary process [24, 26]. A growing urban area and human population density can affect the mosquito populations on a genetic level. An increase in urbanization can reduce genetic diversity and promote homogenization of the genetic structure of mosquito populations. However, when an area reaches complete and/or stable urbanization, the mosquito populations may reach peak abundance, allowing an increase in genetic diversity and thereby, establishing heterogenous genetic structures of mosquito populations [22]; a phenomenon that has been observed in urbanized cities in Brazil [22], Thailand [21], Indonesia [27], Cambodia and Myanmar [19].

Although dengue is endemic in the Philippines, there are a limited number of studies that investigate the population genetics of Ae. aegypti in the country. These studies revealed the seasonal genetic structuring and expansion of the mosquito vector [28] and the connectivity among its populations from different island seaports [29], but provide little insight into genetic structuring of Ae. aegypti in highly urbanized and dengue-endemic regions in the country. In this study, we selected Metropolitan Manila because it represents an ideal setting to investigate how urbanization can shape the present genetic structure of Ae. aegypti . Previous studies conducted in this area have only investigated the association between dengue epidemiology and climate [5,30,31] and the population structure using wing geometric morphometrics [32], offering little information about its genetic structure and relatedness. Therefore, the study collected adult Ae. aegypti samples from several sites within Metropolitan Manila in order to characterize the population genetic structure and relatedness among mosquito populations.

Methods

Study area and mosquito sampling

Metropolitan Manila, a highly urbanized area, is the National Capital Region (NCR) of the Philippines with an area of 636 km2 . It is located at the eastern shore of Manila Bay in Southwestern Luzon (14°50′ N Latitude, 121°E Longitude). It is composed of 16 cities and 1 municipality with a total population of 12,877,253 [33]. This area is the most urbanized region in the Philippines and is the center of the national government, economy, education and culture, business, manufacturing and importation and exportation [34].

In this study, Metropolitan Manila was divided into 21 study areas (Fig 1 and S1 Table). Initially, study areas were delineated per city with at least one study area per city. However, some cities comprised either a large or small land sized area, thus designating more study areas. For example, the largest city, in the region Quezon City was divided into five study areas based on its district boundaries. In order to standardize the land size covered by each study area, we merged two small neighboring cities, San Juan and Mandaluyong. All maps were prepared in ArcGIS software version 10.2.2 (ESRI, NC, USA) for further analysis [35].

Geographic midpoints of Ae. aegypti study areas (colored diamond) with their corresponding household sites (colored circles) in Metropolitan Manila.
Fig 1
Details of each study area can be seen in S1 Table. The map was prepared using ArcGIS version 10.2.2 from LandsatLook Viewer (http://landsatlook.usgs.gov/).Geographic midpoints of Ae. aegypti study areas (colored diamond) with their corresponding household sites (colored circles) in Metropolitan Manila.

Households on each study area were selected based on voluntary informed consent for collecting adult Ae. aegypti mosquitoes inside their premises. The number of households per study area ranged from 2–14 with an average of six households in each area and a total of 151 households. The maximum distance between households within study areas ranged from 1.39–6.17 km. Since the households were widely dispersed across each study area, we calculated the geographical midpoint [36] to assign a single georeferenced location for each study area in the subsequent genetic analysis. The distance among study areas (midpoints) ranged from 2.85–39.66 km.

The majority of the previous population genetic studies of Ae. aegypti have either used only larval or reared larval-to-adult samples which could lead to potential bias in estimating population genetic parameters due to the sampling of full sibling mosquito larvae [37]. This was evidenced in a population genetic structure study of amphibian larval samples that led to inaccurate estimates of differentiation among populations when compared to adult samples [38,39]. As such, this study targets adult Ae. aegypti samples rather than conventional egg, or larval samples to prevent collecting mosquito full siblings. Collection of Ae. aegypti mosquitoes was done by installing a commercially available mosquito UV Light Trap (MosquitoTrap, Jocanima Corporation, Las Pinas City, Philippines) in each household’s outdoor premises for 3–5 days. Collected adult mosquito individuals were sorted, then identified accordingly based on the pictorial keys of Rueda et al. [40] and preserved in 99% ethanol. From May 2014 until January 2015, a total of 526 adult Ae. aegypti were collected, ranging from 10 to 42 individuals per study area (S1 Table).

DNA extraction and microsatellite genotyping

The total genomic DNA of individual mosquito sample was extracted using the QIAGEN Blood and Tissue DNEasy Kit (Qiagen, Hilden, Germany) following the modified protocol of Crane [41]. The average DNA concentration per individual mosquito sample was measured at 41.68 ± 34.50 ng/ μL using NanoDrop2000 Spectrophotometer (Thermo Scientific Inc., Waltham, Massachusetts, USA). We selected 11 microsatellite markers out of the 24 loci as listed by Slotman et al. [42] and Chambers et al. [23] by conducting trial runs (PCR and gel electrophoresis) on a subset of Ae. aegypti samples (n = 8). The selection of the markers was based on the reliability of the amplified products in all samples from different trial runs. The selected markers used in this study were congruent with previous population studies conducted in Southeast Asia [19,21]. The loci were grouped into four sets for multiplex PCR (S2 Table). Each set consisted of fluorescent labeled forward primers with different annealing temperatures during the amplification process. Each set was composed of 1.2 μL of 10X buffer (TaKaRa, Shiga, Japan), 0.8 μL of 25 mM MgCl2 , 1.6 μL of 10 mM of each dNTP, 0.6 μL of 10 μM forward and reverse primers and 0.08 μL 5.0U/ μL of Taq DNA polymerase (TaKaRa), 1.5 μL of 10% dimethyl sulfoxide (DMSO) (Sigma-Aldrich, St. Louis, Missouri, USA) and 1 μL of template DNA with a final volume of 10 μL. Thermocycling conditions were as follows: initial denaturation at 94°C for 5 minutes, denaturation at 94°C for 30 seconds, annealing (see S2 Table for temperatures and durations for each primer set) extension at 72°C for 30 seconds (35 cycles) and a final incubation step of 72°C for 5 minutes. PCR amplicons were checked by electrophoresis in 3% agarose gels and visualized under UV light using the Chemidoc XRS Chemiluminescent Gel Documentation Cabinet (BIO-RAD). Prior to fragment size analyses, multiplex PCR products were diluted in 1/15 water and then pooled together. One microliter of each diluted pool was added with 0.5 μL of GS 500 Liz Internal Size Standard (Applied Biosystems, USA) and HD formamide for a total volume of 20 μL. Fragment analysis of the amplified products were done using ABI 3500 Genetic Analyzer (Life Technologies) while genotyping was conducted in GeneMapper (Applied Biosystems). Microsatellite data were checked for error, and the presence of null alleles was confirmed with MicroChecker [43]. The dataset was deposited and made publicly available at vectorbase.org [44] with the population biology project ID: VBP0000554.

The exact Hardy-Weinberg equilibrium (HWE) test and estimations of the linkage disequilibrium (LD) among all pairs of loci were conducted using GENEPOP v4.2.1 [45,46]. Significance levels for multiple testing were corrected using the Bonferroni procedure. Measures of genetic diversity, including mean number of different alleles (MNa), mean number of effective alleles (MNe), allelic richness (AR), private allelic richness (PAR), observed heterozygosity (Ho), expected heterozygosity (He ), and inbreeding coefficients (Fis) were estimated using Genetic Analysis in Excel (GenAlEx) version 6.3 [47]. To assess the magnitude of genetic differentiation among sites, uncorrected and corrected pairwise FST values were calculated using Arlequin v3.5.1.3 [48] and FreeNA [49], respectively, with 10,000 permutations. The analysis of molecular variance (AMOVA) was performed using GenAlEx version 6.3 [47]. Since the corrected pairwise FST values were calculated by FreeNA [49], we manually calculated the number of pairwise effective migrants using the formula of Slatkin and Barton [50] in GenAlEx version 6.3 [47]. FreeNA [49] was also used to calculate the Cavalli Sforza and Edwards distance. A dendrogram was constructed based on the Cavalli Sforza and Edwards distance using the unweighted pair group with arithmetic mean (UPGMA) in the fastcluster package [51] and the optimal number of clusters were determined using the cindex index in the NbClust package [52] of the R program [53].

Genetic structure

The number of genetic groups (K) was inferred using the Bayesian approach in the software Structure version 2.3.2 [54]. The study used admixture models where the alpha value may vary, and independent allele frequencies set at lambda equals to one. Twenty independent runs were performed for each value of K (1–20) with a burn-in phase of 200,000 iterations followed by 600,000 replications. To determine the most likely number of clusters, we calculated the commonly used ΔK statistic as developed by Evanno et al (2005) [55] using the web application Structure Harvester version 0.6.93 [56]. To avoid the effects of uneven sampling [57,58] in the Bayesian analysis using the software Structure version 2.3.2 [54], we standardized the number of individuals per study area to 10 and created five different datasets. These datasets were subjected to Bayesian analysis using the same parameters and analyzed using the Evanno et al (2005) [55] method. The ‘pophelper’ package [59] in R [53] was used to summarize and visualize the bar plots for the best K statistic identified for each dataset.

Isolation by distance and spatial autocorrelation

Pairwise geographic distances (km) among study areas and households of the mosquito populations were calculated using Vincenty’s formulae [60] in Microsoft Excel 2016. To test isolation by distance (IBD), pairwise FST and geographical distance (km) were examined using Mantel’s test of correlation with 10,000 permutations. Spatial autocorrelation was performed using pairwise Provesti’s genetic distance among mosquito individuals and geographic distance (km) among households with 10,000 permutations and bootstrap replications. Pairwise Provesti’s genetic distance among mosquito individuals was calculated using the ‘poppr’ package [61] in R [53]. Results of the permutation were considered significant at the 5% level. In this analysis, a correlogram was produced with 45 distance classes at 1 km intervals. Both analyses yielded a correlation coefficient of the two data matrices ranging from –1 to +1, with a test for a significant relationship by random permutation. All analyses were performed using GenAlEx version 6.3 [47].

Results

Marker assessment

We observed a total of 113 alleles across 11 microsatellite loci in 21 mosquito populations from Metropolitan Manila (S3 Table). The number of alleles per loci ranged from 3 (F06) to 19 (B07) with an average of 10.25 alleles per loci, suggesting the high polymorphic nature of the selected microsatellites. A total of 91 out of 231 (39%) deviated significantly from the Hardy-Weinberg Equilibrium expectations after sequential Bonferroni correction. Seventy-two of these significant deviations indicated He > Ho, suggesting heterozygosity deficiency (S4 Table). Null alleles were present in four loci (M313, AC4, AG7 and H08) and the null allele frequency ranges from 0.00–0.33 in all loci (S5 Table). The LD test showed a total of 131 out of 1155 (11.30%) pairs of loci with significant LD after Bonferroni corrections, with no locus pair consistently correlated across all mosquito populations (S6 Table).

Genetic diversity and differentiation

Table 1 shows a summary of the genetic diversity of each mosquito population. The mean number of different alleles ranged from 3.82 (TGG) to 6.36 (QZC-3) while the mean number of effective alleles areas ranged from 2.74 (QZC-1) to 3.51 (LSP-2). On the other hand, the mean allelic richness ranged from 3.24 (QZC-1) to 3.85 (MRK) while the proportion of private alleles ranged from 0.02 (LSP-2) to 0.23 (PSY). All mosquito populations except for one (MKT) did not conform to Hardy-Weinberg equilibrium expectations (He > Ho), indicating heterozygosity deficits and the possibility of inbreeding within each mosquito population.

Table 1
Genetic diversity of 21 Ae. aegypti populations based on 11 microsatellites in Metropolitan Manila, Philippines.
Study AreasMNaMNeARPARHoHeFIS
CAL-N5.913.093.680.080.4930.5760.169
CAL-S4.732.933.490.110.4680.5890.215
VAL4.272.933.490.040.5080.5670.143
QZC-14.452.743.240.100.4700.5340.111
QZC-25.822.943.450.060.5200.5450.047
QZC-36.362.923.440.030.5250.5650.071
QZC-46.003.203.630.080.5410.5740.077
QZC-55.553.183.370.050.5000.5650.117
MNL-16.093.303.530.080.5240.5810.090
MNL-25.362.753.550.030.4680.5300.132
MRK5.643.093.850.070.5110.5790.112
SJ-MND5.823.123.290.030.4960.5740.191
PSG4.912.883.560.080.5240.5730.099
MKT5.183.003.590.060.5580.5510.036
PSY5.642.943.710.230.5060.5510.080
TGG3.822.793.270.050.5360.5670.071
PRNQ5.642.903.780.040.5640.5800.076
LSP-15.823.133.650.120.5360.5780.104
LSP-25.003.513.440.020.4320.5890.296
MTLP-15.733.083.790.130.4380.5800.240
MTLP-26.003.023.510.030.4490.5770.240
MNa = mean number of different alleles; MNe = mean number of effective alleles; AR = allelic richness; PAR = private allelic richness; Ho = observed heterozygosity; He = expected heterozygosity; FIS = inbreeding coefficient

The global FST was estimated to be 0.013 where pairwise corrected FST values among mosquito populations ranged from -0.009 to 0.064 (S7 Table). Among these comparisons, 87 out of 201 (41.4%) uncorrected pairwise FST values showed significant genetic differences (S7 Table). The effective number of migrants among mosquito populations ranged from 3.67 to 428.57 migrants (S8 Table) with 13 pairwise estimates were not calculated due to the negative or zero FST values. The results from the AMOVA showed that the FIS and FIT values were 0.131 and 0.137, respectively (both P < 0.001) (S9 Table), indicating that genetic differences exist among and within individuals.

The dendrogram based on Cavalli Sforza and Edwards distance revealed the spatial pattern and distribution of genetically similar mosquito populations (Fig 2). The “cindex” index in NbClust determined the optimal number of cluster groups to be four. One group (green line/circles, Fig 2) had 14 mosquito populations that were highly connected or genetically similar. An extended pattern can be observed from the north down to the southern parts of Metropolitan Manila.

Dendrogram showing the genetic relatedness of each study area based on Cavalli Sforza and Edwards distance.
Fig 2
Colored lines indicate the genetic groups (n = 4). The map showing the study areas in respect to their genetic group assignment (colored circles). Blue lines in the map show the primary road network of the region. The map was prepared using ArcGIS version 10.2.2 from LandsatLook Viewer (http://landsatlook.usgs.gov/).Dendrogram showing the genetic relatedness of each study area based on Cavalli Sforza and Edwards distance.

The Mantel test between the pairwise genetic (FST) and geographic distances of all mosquito populations showed very low and non-significant correlation (R2 = 0.006, p = 0.25), indicating no isolation by distance (S1 Fig). High genetic similarity among adult mosquito individuals was only found up to 1 km based on the spatial autocorrelation analysis (S2 Fig), suggesting limited dispersal capability for Ae. aegypti.

Genetic structure

According to Evanno et al.’s method (2005) [55], the most probable number of genetically differentiated groups using all individuals across mosquito populations is K = 4 (Fig 3, S3 Fig). Furthermore, the standardized datasets showed the most probably numbers of K ranged from 3 to 6 (Fig 3, S3 Fig). All barplots demonstrated no clear genetically-distinct clusters, suggesting high admixture among all mosquito populations.

Individual barplots from the Bayesian analysis using Structure software indicated the distribution of genetic clusters across mosquito populations using (a) all (n = 526) sample individuals and (b) five standardized datasets (n = 10 per study area). Each individual is represented by a single horizontal line. Brackets are shown to separate study areas. The most likely number of genetic clusters in (a) is K = 4 while, in (b) the number ranges from K = 3 to K = 6.
Fig 3
Individual barplots from the Bayesian analysis using Structure software indicated the distribution of genetic clusters across mosquito populations using (a) all (n = 526) sample individuals and (b) five standardized datasets (n = 10 per study area). Each individual is represented by a single horizontal line. Brackets are shown to separate study areas. The most likely number of genetic clusters in (a) is K = 4 while, in (b) the number ranges from K = 3 to K = 6.

Discussion

Our study revealed low genetic differentiation (FST) among mosquito populations across Metropolitan Manila which correlates well with the findings from other population genetic studies of Ae. aegypti in a microgeographic-scale setting [19,21,22]. In a study conducted in Brazil, the level of genetic differentiation ranged from 0.002 to 0.094 with a maximum distance among collection sites of 30 km [22]. Urban cities in Southeast Asian countries with spatial scales of 5–50 km also showed low genetic differentiation (0.026–0.032) [19]. The same outcome (average FST = 0.037) was observed among villages in Thailand with geographical distances up to 27 km [21]. Our findings, along with those of previous studies, suggest that mosquito populations within fine-scale areas consist of similar allele frequencies, exhibiting continuous and active exchange or sharing of alleles. Moreover, our study revealed the lack of isolation-by-distance (IBD) which further implies that mosquito populations are not in equilibrium of migration and genetic drift at this spatial scale.

The lack of IBD and low FST imply high gene flow and weak genetic drift among mosquito populations in Metropolitan Manila. This high gene flow may indicate the possibility of “passive” dispersal of Ae. aegypti throughout the region. Other studies claim that this mosquito vector occasionally travels long distances by taking advantage of human-aided transportation routes via land, sea, or air [26,29,62]. This is evidenced by the eggs, larvae, and adults found in commercial trucks and ships through tire importation [25,63] while larvae and pupae are littered across transportation zones such as airports [64] and docks/ports [29,63]. Such a mechanism could facilitate high gene flow, resulting in genetically-similar mosquito populations. Additionally, spatial autocorrelation revealed limited dispersal distance up to 1 km, suggesting the individual flight range or “active dispersal” of Ae. aegypti [12]. The duality of both “passive” and “active” dispersal has been a common finding among human-adapted mosquitoes, especially in Ae. aegypti found in Southeast Asia [19,21,27, 65] and is considered a key factor in the persistence and resurgence of mosquito-borne diseases [66].

Furthermore, the low genetic differences due to weak genetic drift among mosquito populations are also reported. As evidenced by studies in Vietnam [20] and in the Philippines [28], genetic differences among mosquito populations tend to be lower during the rainy season compared to the dry season. These studies surmise that mosquito populations from the dry season recolonize breeding sites, leading to the expansion and proliferation of the mosquito vector. The large population size at this period can decrease genetic drift which leads to low genetic differences among mosquito populations. A similar phenomenon was observed in areas where insecticide is applied periodically. Mosquito populations in untreated areas facilitated high mosquito abundance and weakened genetic drift, resulting to decreased genetic differentiation among mosquito populations [20].

Our Bayesian analysis inferred sympatric multiple genetic clusters (K = 3 to K = 6) similar to the results reported by other population genetic studies with micro-geographic scales [19,22,67,68] and those conducted in the Philippines [28,29]. A possible explanation for this phenomenon could be the result of divergence from a single ancestry that produced multiple genetic clusters over time. It is argued that a single house or closely situated houses could act as an assembling unit in the formation of multiple genetic clusters in fine-scale areas [19,26,69]. Another possible explanation for multiple genetic clusters could be various invasions or introductions of Ae. aegypti populations from neighboring regions or provinces surrounding Metropolitan Manila. These has been shown in studies in China [70] and the US [71]. We attempted to test this hypothesis by obtaining publicly available genetic data from previous studies in the Philippines [17, 28, 29], however only two microsatellite loci (AC2 and AC4) matched with our set of markers. Furthermore, we expanded our search to neighboring Southeast Asian countries (e.g. Thailand, Vietnam and Indonesia) using the database vectorbase.org [44], but obtained the same outcome of matched loci (AC2 and AC4). This perspective is an interesting avenue for future investigations regarding population genetic studies of Ae. aegypti in the Philippines.

Although this analysis indicated multiple genetic clusters, mixed ancestry is likewise noticeable where each individual mosquito is assigned to different genetic clusters. The bar plots showed clear patterns of genetic admixture across mosquito populations. The interpretations of our Bayesian analysis are consistent with other population genetic studies [19,22,28,29,67,68] which demonstrated low genetic differentiation with no discernable genetic structures. The constant migration of mosquito populations facilitated by passive dispersal (e.g. human-mediated transportation) could have produced the genetic admixture mosquito populations seen in our results.

The constructed dendrogram allows us to further examine the passive dispersal of Ae. aegypti by describing potential dispersal patterns and its impact on the fine-scale genetic structuring. The first group (Fig 2, green line/circles) demonstrates the spanning network of genetically similar mosquito populations from the northern up to the southern parts of the region. The long-distance genetic similarity among these mosquito populations could be due to the major roads that connect the north and south. This supports the claim that long-distance dispersal of the mosquito vector could be facilitated through human-aided transportation.

Interestingly, the mosquito populations in the eastern (Fig 2, red line/circles) and northwestern (Fig 2, yellow line/circles) areas highlight the patterns of dispersal and features that could affect fine-scale genetic structuring. These areas have a unique landscape feature of being marginally separated from neighboring cities. For example, a large highway and river split a sizeable land area in cities of Marikina (MRK) and Pasig (PSG) while large highways (e.g., expressways) divide Valenzuela (VAL) and South Caloocan (CAL-S). As a result, only a few roads and bridges connect to these cities which limits access to neighboring areas. This circumstance may limit migration between mosquito populations, resulting in formation of distinct genetic groups as shown in the dendrogram (MRK & PSG and VAL & CAL-S). Landscape fragmentation can limit gene flow among mosquito populations as evidenced in the West Indies, producing defined mosquito populations on the east and west side of a major highway [72]. This could also explain the TGG mosquito population from Taguig City since the land area is separated by two large highways and a river which limits access to neighboring cities (Fig 2, black line/circle).

The genetic grouping shown in the dendrogram can also be explained by the occurrence and magnitude of flooding in selected cities within the region. From 2009 to 2014, Metropolitan Manila had experienced several flooding events which displaced residents and damaged properties and businesses [73]. Porio [74] identified the cities of Taguig, Pasig, Marikina, Valenzuela, and Caloocan as flood-prone areas, among many others. Because of these cities marginally separated, we hypothesize that, after a flooding event, alleles within these populations could be generally preserved and its allelic combination may become distinct from other mosquito populations over time due to limited gene flow.

Our findings along with those from previous reports in the Philippines [28,29] provide a partial picture of the current population genetic status of the mosquito vector, Ae. aegypti, in the country. The low genetic differentiation observed in all the reports including ours, indicates its widespread distribution not only within cities and villages but also in different islands of the country. Human-mediated transportation (marine or road transportation) facilitate the long dispersal capability of Ae. aegypti, enabling the observed admixture within the mosquito population. This comprehensive illustration of the population biology of Ae. aegypti in the Philippines strengthens our understanding for developing better vector control strategies. If and when the country pursues innovative vector control approaches, these findings will be integral to the success of the program.

References

1 

SC Weaver, WK Reisen. . Present and future arboviral threats. Antiviral research. 201021;85(2):, pp.328–45. , doi: 10.1016/j.antiviral.2009.10.008

2 

World Health Organization: Western Pacific Region (2012) Dengue: dengue in the Western Pacific region. Available: http://www.wpro.who.int/emerging_diseases/Dengue/en/index.html. Accessed June 2017

3 

FE Edillo, YA Halasa, FM Largo, JN Erasmo, NB Amoin, MT Alera, IK Yoon, AC Alcantara, DS Shepard. . Economic cost and burden of dengue in the Philippines. The American journal of tropical medicine and hygiene. 201524;92(2):, pp.360–6. , doi: 10.4269/ajtmh.14-0139

4 

TM Carvajal, KM Viacrusis, LF Hernandez, HT Ho, DM Amalin, K Watanabe. . Machine learning methods reveal the temporal pattern of dengue incidence using meteorological factors in metropolitan Manila, Philippines.BMC infectious diseases.201812;18(1):, pp.183, doi: 10.1186/s12879-018-3066-0

5 

L Bravo, VG Roque, J Brett, R Dizon, M L'Azou. . Epidemiology of dengue disease in the Philippines (2000–2011): a systematic literature review.PLoS neglected tropical diseases.2014116;8(11):, pp.e3027, doi: 10.1371/journal.pntd.0003027

6 

Sanofi Pasteur. Sanofi Pasteur’s Dengue Vaccine Approved in the Philippines. Available at: http://www.sanofipasteur.com/en/articles/sanofi-pasteur-dengue-vaccine-approved-in-thephilippines.aspx (2016). Accessed 26 Oct 2017

7 

LC Farnesi, AJ Martins, D Valle, GL Rezende. . Embryonic development of Aedes aegypti (Diptera: Culicidae): influence of different constant temperatures.Memórias do Instituto Oswaldo Cruz.20092;104(1):, pp.124–6. , doi: 10.1590/s0074-02762009000100020

8 

TW Scott, PH Amerasinghe, AC Morrison, LH Lorenz, GG Clark, D Strickman, P Kittayapong, JD Edman. . Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: blood feeding frequency.Journal of medical entomology. 20001;37(1):, pp.89–101. , doi: 10.1603/0022-2585-37.1.89

9 

LC Harrington, JD Edman, TW Scott. . Why do female Aedes aegypti (Diptera: Culicidae) feed preferentially and frequently on human blood?.Journal of medical entomology. 20015;38(3):, pp.411–22. , doi: 10.1603/0022-2585-38.3.411

10 

LM Rueda, KJ Patel, RC Axtell, RE Stinner. . Temperature-dependent development and survival rates of Culex quinquefasciatus and Aedes aegypti (Diptera: Culicidae).Journal of medical entomology. 199091;27(5):, pp.892–8. , doi: 10.1093/jmedent/27.5.892

11 

W Tun‐Lin, TR Burkot, BH Kay. . Effects of temperature and larval diet on development rates and survival of the dengue vector Aedes aegypti in north Queensland, Australia. Medical and veterinary entomology. 200031;14(1):, pp.31–7. , doi: 10.1046/j.1365-2915.2000.00207.x

12 

P Reiter, MA Amador, RA Anderson, GG Clark. . Dispersal of Aedes aegypti in an urban area after blood feeding as demonstrated by rubidium-marked eggs. The American journal of tropical medicine and hygiene. 199521;52(2):, pp.177–9. , doi: 10.4269/ajtmh.1995.52.177

13 

LC Harrington, TW Scott, K Lerdthusnee, RC Coleman, A Costero, GG Clark, JJ Jones, S Kitthawee, P Kittayapong, R Sithiprasasna, JD Edman. . Dispersal of the dengue vector Aedes aegypti within and between rural communities. The American journal of tropical medicine and hygiene. 200521;72(2):, pp.209–20.

14 

NA Honório, WD Silva, PJ Leite, JM Gonçalves, LP Lounibos, R Lourenço-de-Oliveira. . Dispersal of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in an urban endemic dengue area in the State of Rio de Janeiro, Brazil.Memórias do Instituto Oswaldo Cruz.20033;98(2):, pp.191–8. , doi: 10.1590/s0074-02762003000200005

15 

S Christophers. . Aedes aegypti (L.) the yellow fever mosquito: its life history, bionomics and structure.Aëdes aegypti (L.) the Yellow Fever Mosquito: its Life History, Bionomics and Structure.1960.

16 

DJ Gubler. . Dengue, urbanization and globalization: the unholy trinity of the 21st century. Tropical medicine and health. 2011;39(4SUPPLEMENT):, pp.S3–11.

17 

JE Brown, CS McBride, P Johnson, S Ritchie, C Paupy, H Bossin, J Lutomiah, I Fernandez-Salas, A Ponlawat, AJ Cornel, WC Black. . Worldwide patterns of genetic differentiation imply multiple ‘domestications’ of Aedes aegypti, a major vector of human diseases.Proceedings of the Royal Society of London B: Biological Sciences. 2011112:, pp.rspb20102469.

18 

L Urdaneta-Marquez, AB Failloux. . Population genetic structure of Aedes aegypti, the principal vector of dengue viruses. Infection, Genetics and Evolution. 201131;11(2):, pp.253–61. , doi: 10.1016/j.meegid.2010.11.020

19 

T Hlaing, W Tun‐Lin, P Somboon, D Socheat, T Setha, S Min, S Thaung, O Anyaele, B De Silva, MS Chang, A Prakash. . Spatial genetic structure of Aedes aegypti mosquitoes in mainland Southeast Asia. Evolutionary applications. 201071;3(4):, pp.319–39. , doi: 10.1111/j.1752-4571.2009.00113.x

20 

K Huber, L Le Loan, TH Hoang, TK Tien, F Rodhain, AB Failloux. . Temporal genetic variation in Aedes aegypti populations in ho chi Minh City (Vietnam).Heredity. 20027;89(1):, pp.7, doi: 10.1038/sj.hdy.6800086

21 

P Olanratmanee, P Kittayapong, C Chansang, AA Hoffmann, AR Weeks, NM Endersby. . Population genetic structure of Aedes (Stegomyia) aegypti (L.) at a micro-spatial scale in Thailand: implications for a dengue suppression strategy.PLoS neglected tropical diseases.2013110;7(1):, pp.e1913, doi: 10.1371/journal.pntd.0001913

22 

AB Wilke, R Wilk-da-Silva, MT Marrelli. . Microgeographic population structuring of Aedes aegypti (Diptera: Culicidae).PloS one.2017920;12(9):, pp.e0185150, doi: 10.1371/journal.pone.0185150

23 

EW Chambers, JK Meece, JA McGowan, DD Lovin, RR Hemme, DD Chadee, K McAbee, SE Brown, DL Knudson, DW Severson. . Microsatellite isolation and linkage group identification in the yellow fever mosquito Aedes aegypti. Journal of Heredity. 200749;98(3):, pp.202–10. , doi: 10.1093/jhered/esm015

24 

S Ravel, N Monteny, DV Olmos, JE Verdugo, G Cuny. . A preliminary study of the population genetics of Aedes aegypti (Diptera: Culicidae) from Mexico using microsatellite and AFLP markers.Acta Tropica. 2001330;78(3):, pp.241–50. , doi: 10.1016/s0001-706x(01)00083-3

25 

DD Chadee. . Seasonal incidence and horizontal distribution patterns of oviposition by Aedes aegypti in an urban environment in Trinidad, West Indies. Journal of the American Mosquito Control Association. 199291;8(3):, pp.281–4.

26 

JE Brown, BR Evans, W Zheng, V Obas, L Barrera‐Martinez, A Egizi, H Zhao, A Caccone, JR Powell. . Human impacts have shaped historical and recent evolution in Aedes aegypti, the dengue and yellow fever mosquito. Evolution. 20142;68(2):, pp.514–25. , doi: 10.1111/evo.12281

27 

G Rašić, N Endersby-Harshman, W Tantowijoyo, A Goundar, V White, Q Yang, I Filipović, P Johnson, AA Hoffmann, E Arguni. . Aedes aegypti has spatially structured and seasonally stable populations in Yogyakarta, Indonesia.Parasites & vectors.2015Dec;8(1):, pp.610.

28 

SL Sayson, A Gloria-Soria, JR Powell, FE Edillo. . Seasonal genetic changes of Aedes aegypti (Diptera: Culicidae) populations in selected sites of Cebu City, Philippines.Journal of medical entomology. 2015814;52(4):, pp.638–46. , doi: 10.1093/jme/tjv056

29 

E Fonzi, Y Higa, AG Bertuso, K Futami, N Minakawa. . Human-mediated marine dispersal influences the population structure of Aedes aegypti in the Philippine Archipelago. PLoS neglected tropical diseases. 201563;9(6):, pp.e0003829, doi: 10.1371/journal.pntd.0003829

30 

TM Carvajal, KM Viacrusis, LF Hernandez, HT Ho, DM Amalin, K Watanabe. . Machine learning methods reveal the temporal pattern of dengue incidence using meteorological factors in metropolitan Manila, Philippines.BMC infectious diseases.201812;18(1):, pp.183, doi: 10.1186/s12879-018-3066-0

31 

A Sumi, EF Telan, H Chagan-Yasutan, MB Piolo, T Hattori, N Kobayashi. . Effect of temperature, relative humidity and rainfall on dengue fever and leptospirosis infections in Manila, the Philippines.Epidemiology & Infection.20171;145(1):, pp.78–86.

32 

TM Carvajal, LF Hernandez, HT Ho, MG Cuenca, BM Orantia, CR Estrada, KM Viacrusis, DM Amalin, K Watanabe. . Spatial analysis of wing geometry in dengue vector mosquito, Aedes aegypti (L.)(Diptera: Culicidae), populations in Metropolitan Manila, Philippines. Journal of vector borne diseases. 201661;53(2):, pp.127

33 

Philippine Statistics Authority: Population and Housing. http://psa.gov.ph/ (2016). Accessed on Jun 2016

34 

Asia Development Bank. Urban metabolism of six Asian cities: Metro Manila. Mandaluyong City, Philippines: Asian Development Bank2014 pp. , pp.29–34

35 

ESRI. 2014. ArcGIS Desktop: Release 10. Redlands, CA: Environmental Systems Research Institute

36 

GeoMidpoint.com. Midpoint Calculator. 2007. Accessed August 2017. Available from: http://www.geomidpoint.com/

37 

FW Allendorf, SR Phelps. . Use of allelic frequencies to describe population structure. Canadian Journal of Fisheries and Aquatic Sciences. 1981121;38(12):, pp.1507–14.

38 

W Peterman, ER Brocato, RD Semlitsch, LS Eggert. . Reducing bias in population and landscape genetic inferences: the effects of sampling related individuals and multiple life stages.PeerJ.2016314;4:, pp.e1813, doi: 10.7717/peerj.1813

39 

CS Goldberg, LP Waits. . Quantification and reduction of bias from sampling larvae to infer population and landscape genetic structure. Molecular Ecology Resources. 20103;10(2):, pp.304–13. , doi: 10.1111/j.1755-0998.2009.02755.x

40 

LM Rueda. Pictorial keys for the identification of mosquitoes (Diptera: Culicidae) associated with dengue virus transmission.Walter Reed Army Inst Of Research Washington Dc Department Of Entomology; 200483.

41 

Crane S. DNA Extraction From Archival Museum Insect Specimens. 2011. Accessed Mar 2015. Available from: https://s3-eu-west-1.amazonaws.com/pfigshare-u-files/1114092/extractionmuseum.pdf

42 

MA Slotman, NB Kelly, LC Harrington, S Kitthawee, JW Jones, TW Scott, A Caccone, JR Powell. . Polymorphic microsatellite markers for studies of Aedes aegypti (Diptera: Culicidae), the vector of dengue and yellow fever.Molecular Ecology Resources. 200711;7(1):, pp.168–71.

43 

C Van Oosterhout, WF Hutchinson, DP Wills, P Shipley. . MICRO‐CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Resources. 200491;4(3):, pp.535–8

44 

GI Giraldo-Calderón, SJ Emrich, RM MacCallum, G Maslen, E Dialynas, P Topalis, N Ho, S Gesing, VectorBase Consortium, G Madey, FH Collins. . VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic acids research. 20141215;43(D1):, pp.D707–13.

45 

F Rousset. . genepop’007: a complete re‐implementation of the genepop software for Windows and Linux. Molecular ecology resources. 200811;8(1):, pp.103–6. , doi: 10.1111/j.1471-8286.2007.01931.x

46 

F Rousset. . GENEPOP (Version 1.2): Population genetics software for exact tests and ecumenicalism.J. Hered. 1995; 83: , pp.239.

47 

R Peakall, PE Smouse. . GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and researchdan update. Bioinformatics28, , pp.2537e2539. , doi: 10.1093/bioinformatics/bts460

48 

L Excoffier, HE Lischer. . Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular ecology resources. 201051;10(3):, pp.564–7. , doi: 10.1111/j.1755-0998.2010.02847.x

49 

MP Chapuis, A Estoup. . Microsatellite null alleles and estimation of population differentiation. Molecular biology and evolution. 2006125;24(3):, pp.621–31. , doi: 10.1093/molbev/msl191

50 

M Slatkin, NH Barton. . A comparison of three indirect methods for estimating average levels of gene flow. Evolution. 1989111;43(7):, pp.1349–68. , doi: 10.1111/j.1558-5646.1989.tb02587.x

51 

D Müllner. . fastcluster: Fast hierarchical, agglomerative clustering routines for R and Python.Journal of Statistical Software. 2013529;53(9):, pp.1–8.

52 

M Charrad, N Ghazzali, V Boiteau, A Niknafs, MM Charrad. . Package ‘NbClust’. Journal of Statistical Software. 20141015;61:, pp.1–36.

53 

R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria2016.

54 

JK Pritchard, M Stephens, P Donnelly. . Inference of population structure using multilocus genotype data. Genetics. 200061;155(2):, pp.945–59.

55 

G Evanno, S Regnaut, J Goudet. . Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular ecology. 20057;14(8):, pp.2611–20. , doi: 10.1111/j.1365-294X.2005.02553.x

56 

DA Earl. . STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method.Conservation genetics resources. 201261;4(2):, pp.359–61.

57 

SJ Puechmaille. . The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem. Molecular Ecology Resources. 20165;16(3):, pp.608–27. , doi: 10.1111/1755-0998.12512

58 

J Wang. . The computer program structure for assigning individuals to populations: easy to use but easier to misuse. Molecular ecology resources. 20179;17(5):, pp.981–90. , doi: 10.1111/1755-0998.12650

59 

RM Francis. . pophelper: an R package and web app to analyse and visualize population structure. Molecular ecology resources. 20171;17(1):, pp.27–32. , doi: 10.1111/1755-0998.12509

60 

T Vincenty. . Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations.Survey review.197541;23(176):, pp.88–93.

61 

ZN Kamvar, JF Tabima, NJ Grünwald. . Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction.PeerJ.2014; 2: , pp.e281, doi: 10.7717/peerj.281

62 

A Egizi, J Kiser, C Abadam, DM Fonseca. . The hitchhiker's guide to becoming invasive: exotic mosquitoes spread across a US state by human transport not autonomous flight. Molecular ecology. 201671;25(13):, pp.3033–47. , doi: 10.1111/mec.13653

63 

M Suleman, M Arshad, K Khan. . Yellowfever mosquito (Diptera: Culicidae) introduced into Landi Kotal, Pakistan, by tire importation.Journal of medical entomology. 199671;33(4):, pp.689–93. , doi: 10.1093/jmedent/33.4.689

64 

N Sukehiro, N Kida, M Umezawa, T Murakami, N Arai, T Jinnai, S Inagaki, H Tsuchiya, H Maruyama, Y Tsuda. . First report on invasion of yellow fever mosquito, Aedes aegypti, at Narita International Airport, Japan in August 2012. Japanese journal of infectious diseases. 2013;66(3):, pp.189–94. , doi: 10.7883/yoken.66.189

65 

C Paupy, N Chantha, K Huber, N Lecoz, JM Reynes, F Rodhain, AB Failloux. . Influence of breeding sites features on genetic differentiation of Aedes aegypti populations analyzed on a local scale in Phnom Penh Municipality of Cambodia. The American journal of tropical medicine and hygiene. 200471;71(1):, pp.73–81

66 

AM Lutambi, MA Penny, T Smith, N Chitnis. . Mathematical modelling of mosquito dispersal in a heterogeneous environment. Mathematical biosciences. 201321;241(2):, pp.198–216. , doi: 10.1016/j.mbs.2012.11.013

67 

PO Vidal, L Suesdek. . Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti. Infection, Genetics and Evolution. 201241;12(3):, pp.591–6. , doi: 10.1016/j.meegid.2011.11.013

68 

SA Guagliardo, Y Lee, AA Pierce, J Wong, YY Chu, AC Morrison, H Astete, B Brosi, G Vazquez-Prokopec, TW Scott, U Kitron. . The genetic structure of Aedes aegypti populations is driven by boat traffic in the Peruvian Amazon. PLoS neglected tropical diseases. 2019918;13(9):, pp.e0007552, doi: 10.1371/journal.pntd.0007552

69 

A Getis, AC Morrison, K Gray, TW Scott. Characteristics of the spatial pattern of the dengue vector, Aedes aegypti, in Iquitos, Peru.InPerspectives on Spatial Data Analysis2010 (pp. , pp.203–225). Springer, Berlin, Heidelberg.

70 

QM Shi, HD Zhang, G Wang, XX Guo, D Xing, YD Dong, L Xiao, J Gao, QM Liu, AJ Sun, CX Li. . The genetic diversity and population structure of domestic Aedes aegypti (Diptera: Culicidae) in Yunnan Province, southwestern China.Parasites & vectors.201712;10(1):, pp.292.

71 

E Pless, A Gloria-Soria, BR Evans, V Kramer, BG Bolling, WJ Tabachnick, JR Powell. . Multiple introductions of the dengue vector, Aedes aegypti, into California.PLoS neglected tropical diseases. 2017810;11(8):, pp.e0005718, doi: 10.1371/journal.pntd.0005718

72 

RR Hemme, CL Thomas, DD Chadee, DW Severson. . Influence of urban landscapes on population dynamics in a short-distance migrant mosquito: evidence for the dengue vector Aedes aegypti.PLoS Neglected Tropical Diseases.2010316;4(3):, pp.e634, doi: 10.1371/journal.pntd.0000634

73 

BC Balgos, LB dela Cruz, VP Valenzuela. . Preparing Metro Manila toward Urban Resiliency: Prospects and Retrospect.Prepared for the Global Assessment Report on Disaster Risk Reduction.2014

74 

E Porio. . Vulnerability, adaptation, and resilience to floods and climate change-related risks among marginal, riverine communities in Metro Manila.Asian Journal of Social Science. 201111;39(4):, pp.425–45.

https://www.researchpad.co/tools/openurl?pubtype=article&doi=10.1371/journal.pntd.0008279&title=Fine-scale population genetic structure of dengue mosquito vector, <i>Aedes aegypti</i>, in Metropolitan Manila, Philippines&author=Thaddeus M. Carvajal,Kohei Ogishi,Sakiko Yaegeshi,Lara Fides T. Hernandez,Katherine M. Viacrusis,Howell T. Ho,Divina M. Amalin,Kozo Watanabe,Pattamaporn Kittayapong,&keyword=&subject=Research Article,Biology and Life Sciences,Evolutionary Biology,Population Genetics,Biology and Life Sciences,Genetics,Population Genetics,Biology and Life Sciences,Population Biology,Population Genetics,People and Places,Geographical Locations,Asia,Philippines,Medicine and Health Sciences,Infectious Diseases,Disease Vectors,Insect Vectors,Mosquitoes,Aedes Aegypti,Biology and Life Sciences,Species Interactions,Disease Vectors,Insect Vectors,Mosquitoes,Aedes Aegypti,Biology and Life Sciences,Organisms,Eukaryota,Animals,Invertebrates,Arthropoda,Insects,Mosquitoes,Aedes Aegypti,Biology and Life Sciences,Genetics,Genetic Loci,Biology and Life Sciences,Evolutionary Biology,Evolutionary Processes,Genetic Drift,Biology and Life Sciences,Evolutionary Biology,Population Genetics,Genetic Drift,Biology and Life Sciences,Genetics,Population Genetics,Genetic Drift,Biology and Life Sciences,Population Biology,Population Genetics,Genetic Drift,Biology and Life Sciences,Genetics,Heredity,Gene Flow,Biology and Life Sciences,Evolutionary Biology,Population Genetics,Gene Flow,Biology and Life Sciences,Genetics,Population Genetics,Gene Flow,Biology and Life Sciences,Population Biology,Population Genetics,Gene Flow,People and Places,Geographical Locations,Asia,Biology and Life Sciences,Genetics,Gene Types,Microsatellite Loci,