J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 546 http://jad.tums.ac.ir Published Online: October 04, 2016 Original Article Genetic Analysis of Aedes aegypti Using Random Amplified Polymorphic DNA (RAPD) Markers from Dengue Outbreaks in Pakistan Hafiz Muhammad Ashraf 1, *Muhammad Kashif Zahoor 1, Shabab Nasir 1, Humara Naz Majeed 2, Sarwat Zahoor 3 1Department of Zoology, Government College University Faisalabad, Islamabad, Pakistan 2Department Pharmaceutical and Biological Chemistry, School of Pharmacy, University College London, London, UK 3Centre of Agricultural Biochemistry and Biotechnology (CABB), University of Agriculture Faisalabad, Faisalabad, Pakistan (Received 14 June 2014; accepted 28 July 2015) Abstract Background: Keeping in view the havoc situation of dengue fever in Pakistan, the current study was designed to demonstrate the genetic variations, gene flow and rate of migration from Lahore and Faisalabad. Methods: The larvae were collected from both natural and artificial breeding places from each collection site. The adult mosquitoes were collected by means of sweep net and battery-operated aspirator. DNA extraction was per- formed using TNE buffer method. Ten GeneLink-A series RAPD primers were used for PCR amplification and the data was analyzed through POPGENE. Results: The number of amplification products produced per primer varied from 8-12, ranging from 200 to 2000 bp with an average of 10.0 bands per primer. The percentage of polymorphic loci amplified by each primer varied from 22.5 to 51%. The UPGMA dendrogram demonstrates two distinct groups from Faisalabad and Lahore populations. The genetic diversity ranged from 0.260 in Faisalabad to 0.294 in Lahore with a total heterozygosity of 0.379. The GST value for nine populations within Lahore was 0.131 (Nm= 3.317), whereas for nine populations in Faisalabad GST value was 0.117 (Nm= 3.773). The overall genetic variation among eighteen populations showed GST= 0.341 and Nm= 1.966. Conclusion: The genetic relatedness and Nm value show that Ae. aegypti populations exhibit intra-population gene flow both in Faisalabad and Lahore. Although, both cities show a distinct pattern of genetic structure; however, few areas from both the cities show genetic similarity. The gene flow and the genetic relatedness in few populations of Lahore and Faisalabad cities need further investigation. Keywords: Aedes aegypti, Dengue, Genetic diversity, RAPD markers Introduction Aedes aegypti and Ae. albopictus, mos- quitoes are responsible to transmit dengue fever virus which is a major health concern in tropical and subtropical areas (Pinheiro and Corber 1997, Fradin and Day 2002, Gibbons and Vaughn 2002). World Health Organiza- tion (WHO) reported that about 2.5 billion people (more than 40%) are at risk of den- gue fever around the world and about 50 mil- lion new cases of dengue would be estimated mated to occur every year (WHO 2009, 2014). Dengue disease is considered endem- ic in more than 100 countries, with South- East Asia as a major affected region (Itrat et al. 2008, Khan et al. 2010, Raheel et al. 2011, Fatima et al. 2011, Idrees and Ashfaq 2012). Aedes aegypti originated from Africa and became spread to the Southeast Asia. With the passage of time, it established here and now it has become a major vector of dengue **Corresponding author: Dr Muhammad Kashif ZZahoor, E-mail: kashif.zahoor@gcuf.edu.pk J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 547 http://jad.tums.ac.ir Published Online: October 04, 2016 in the region (WHO 1997, Fraga at al. 2003). Better road infrastructure and good trans- portation facilities caused dispersal of Ae. aegypti. It can easily find breeding places in human dwelling areas; a calamity of living style of people in the sub-continent. In addi- tion, it can withstand desiccation, which also serves as a factor for its success (Raheel et al. 2011, Sajid et al. 2011). Thus, the spread mosquito and thereby the virus, has led to a resurgence of dengue (Cologna et al. 2005). The climatic change and use of insecti- cides cause appreciable effect on the genetic population structure of Ae. aegypti (Huber et al. 2002, Paupy et al. 2005, Scarpassa et al. 2008). This vector has been considered as ho- mogeneous species, however, it shows both morphological and genetic variations (Mor- lais et al. 2003). Due to epidemiological im- portance, various molecular studies have been conducted on Ae. aegypti using different DNA markers (Ayres et al. 2002, Gorrochotegui et al. 2002, Ravel et al. 2002, Ayres et al. 2003, Fraga et al. 2003, Beebe et al. 2005, Herrera et al. 2006, Bracco et al. 2007, Rasheed et al. 2013). Random amplified polymorphic DNA (RAPD) and RFLP markers are most com- monly used molecular markers to elucidate the genetic variations in Ae. aegypti (Severson et al. 1999, 2002, Fulton et al. 2001, Ayres et al. 2003, Julio et al. 2009). Unfortunately, the genetic analysis of Ae aegypti has been ne- glected in Southeast Asia despite of its seri- ous damaging aspect of dengue outbreaks. In Pakistan, dengue fever has emerged as havoc for the last few years and has become one of the key public health priorities. From the first diagnosed case in 1994, now the vi- rus has rapidly spread across the country. The dengue fever has shaped a cyclic pattern with high infection rate in the monsoon sea- son and low in winter season. More than 2,000 people were reported with dengue in- fection in 2010. Over 14,000 cases were re- ported with more than 300 deaths from La- hore alone in November 2011, likewise, 11,500 people were affected until September 2012 (Gilani 2012, Sajid et al. 2012). The latest outbreak resulted in 18,000 cases across the country including significant number of con- firmed cases reported from Lahore and Fai- salabad (Siddiqui et al. 2009, Gilani 2012, Sajid et al. 2012). According to WHO reports, ma- jority of infected people belonged to La- hore area in Punjab, Pakistan (WHO 2013). To counter this threat and to devise ef- fective control measures, there is essentially a dire need of time to know maximum about dengue and its vector. Despite the fact of not significant ecological differences in both cit- ies, more reported cases of dengue reflect the severity of disease in Lahore. Therefore, it was hypothesized genetic variation might exist in Ae. aegypti populations in Lahore. Thus, the current study was designed to an- alyze genetically the Ae. aegypti populations to demonstrate the genetic variations, gene flow and rate of migration from Lahore and Faisalabad, Punjab using random amplified polymorphic DNA (RAPD) markers. Materials and Methods Mosquito collection Eighteen populations were collected from Faisalabad and Lahore (Fig. 1). The adult mosquitoes were collected with the help of sweep net and battery-operated aspirator (Herrel et al. 2001, Shortall et al. 2009, Flor- encio et al. 2014). The larvae and pupae of mosquitoes were collected from both natural and artificial breeding places (tree holes, stagnant water, tires, waste material, lawns, homes) by using dipper from each collection site (Naeem-Ullah et al. 2010, Nikookar et al. 2010). The collected samples were preserved in 70% alcohol in labelled vials and stored at 4 °C for DNA extraction (Zahoor et al. 2013). DNA extraction DNA was extracted from individual mos- J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 548 http://jad.tums.ac.ir Published Online: October 04, 2016 quitoes from each sample using salt extrac- tion method. The mosquitoes were homog- enized in 400 μl of TNE buffer (Tris-NaCl- EDTA), and then 100 µ l of 20 µ g/µ l of Pro- teinase-K and 40 µ l of 20% Sodium Dodecyl Sulphate (SDS) were added. The homoge- nates were incubated at 55 ˚C for one hour, 300 μl of 5M NaCl was added and vortexed. The mixture was centrifuged at 15,000 rpm for 10 min and the supernatant was shifted to separate eppendorf tube. DNA was precipi- tated by adding of 300-400 µ l isopropanol or ice-cold 100% ethanol, kept at -21 ˚C for 1 hour, and then centrifuged at 15,000 rpm for 10 min. DNA pellet was washed with chilled 70% ethanol, air-dried and re-suspended in 50 µ l of sterile water (d3H2O). The DNA concentration was estimated by measuring optical density (OD) at 260nm and its qual- ity was checked on 1% agarose gel electro- phoresis (Ashraf 2013, Zahoor et al. 2013, Bibi 2015). RAPD-PCR analysis Ten Gene Link-A series RAPD primers were used for PCR amplification (Table 4). Each PCR reaction was carried out in a final volume 25 µ l, approximately 100 mg of ge- nomic DNA, 3 mM of MgCl2, 20 pmol of primer, 2.5 µ l buffer, 1.0 units of Taq DNA polymerase and 0.3 mM of each dNTPs. PCR program comprised of 35 cycles with initial denaturing of DNA at 94 ˚C for 5 min, denaturation at 94 ˚C for 1 min, primer an- nealing at 36 ˚C for 1.5 min, extension at 72 ˚C for 1.5 min, final extension at 72 ˚C for 10 min and then hold at 4 ˚C until the tubes were removed. The PCR products were run on 1.6% agarose gel electrophoresis at 80 voltages for one hour (Ashraf 2013, Zahoor et al. 2013). Statistical Analysis of Data The fingerprints examined under ultra vi- olet Transilluminator and photographed using (SynGene) Gel documentation system. The size of amplified bands or loci was compared against the DNA marker. All the fragments scored as present (1) or absent (0) for each sample. Ambiguous bands were not scored. The bands counted by starting from top to the bottom in all lanes. The RAPD markers were analyzed using the following assumptions: (1) The RAPD alleles follow the Mendelian pattern of seg- regation, (2) monomorphic fragments are ho- mologous (co-migrate), (3) loci are independ- ent from each other and (4) and mosquito populations are in Hardy-Weinberg equilib- rium as described by Ayres et al. 2003. A dendrogram was constructed using UPGMA method (unweight pair-group mean analysis) (Nei 1978). Effective migration rates (Nm) were estimated based on inbreeding indices (GST), where Nm= 0.5(1- GST)/ GST (Ayres et al. 2002). The genetic analysis was performed through POPGENE (version 1.32) software. Results Aedes aegypti populations were analyzed by RAPD-PCR using 10 oligonucleotide pri- mers (Genelink, Table 4). Four primers were further selected which produced distinct and easily detectable bands (Table 1). The repro- ducible bands over repeated PCR reactions were used for analysis. The number of frag- ments per primer varied from 8–12, ranging from 200 to 2000 bp size, with an average of 10.0 bands per primer (Table 1 and Fig. 1). Various amplified fragments differentiated the populations of Ae. aegypti from Lahore and Faisalabad (Fig. 2). The percentage of polymorphic loci amplified by each primer varied from 22.5 to 51% (Fig. 3). Total 55 polymorphic loci were scored in all samples, with an average 5.5 polymorphic loci per primer (data not shown). Four primers were selected for further due to their reproducible and distinct banding pattern (Table 2). The heterozygosity, gene flow and the variation J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 549 http://jad.tums.ac.ir Published Online: October 04, 2016 among eighteen populations were estimated using Nei’s analysis. The genetic diversity ranged from 0.260 in Faisalabad to 0.294 in Lahore with a total heterozygosity of 0.379. The GST value for nine populations within Lahore was 0.131 (Nm= 3.317), whereas for nine populations in Faisalabad GST value was 0.117 (Nm= 3.773), respectively. The over- all genetic variation among eighteen popula- tions showed GST= 0.341 and no. of migrants Nm= 1.966 (Table 3). Cluster Analysis The UPGMA dendrogram based on sim- ilarity matrix of Nei’s genetic distances demon- strates two distinct groups from Faisalabad and Lahore populations (Fig. 4). One group (A) shows seven populations from main Fai- salabad City, ie. Nazim Abad, Gulfishan Col- ony, University of Agriculture Faisalabad (UAF) Jail road, Government College Uni- versity Faisalabad (GCUF), Awami Colony, Narwala road, National Institute for Agricul- ture and Biology (NIAB). The second group (B) represents four populations from Lahore, i.e. Govt. Muhammadan Anglo Oriental Col- lege (MAO College), Qila Gujersing, UC69 and UC 81(A). There are three subgroups in group C, comprising of two populations from Faisalabad (C2, Rescue office Jhang road and AARI Jhang road) and four populations from Lahore (C3, UC 73(A) and Ghawal man- di, and and C4, Mazang chungi and UC73 (B), respectively). The BISE Lawrence road shows that both the groups A and B were originated from Lahore. The subgroup A5 and A4 comprising of UAF Jail road, Gulfishan Colony and GCUF found to be closely related population groups in Faisalabad The subgroup A3 consists of populations from Narwala road and NIAB. The genetic similarity could be speculated due to the transportation of mosquitoes in these areas. The subgroup B2 consists of closely related populations from Gujersing area and MAO College. This subgroup shows more genetic similarity to the mosquito pop- ulation from UC69 site. The subgroup B1 represents relatedness to the subgroup B2, thereby shows more genetic similarity to Gujersing area and MAO College mosquito populations. Within the main group C, there are three subgroups converged at C1 (C2, C3 and C4). The subgroup C3 shows that UC 73(A) is closely related to the mosquitoes from Gawal mandi area. Similarly, the subgroup C4 rep- resents that Mazang chungi and UC 73(B) are more closely related to each other. Inter- estingly, the main group C also contains pop- ulations from Faisalabad. The subgroup C1 consists of areas from the main high way of Jhang road, Rescue office and AARI. The collected mosquito populations show a ge- netic similarity to Lahore populations. The RAPD primer (A-04) produced frag- ments, which showed genetic variation in Ae. aegypti populations from both Lahore and Faisalabad. The genetic relatedness in Fai- salabad populations and the Nm value for an estimate of gene flow shows that the dengue vector Ae. aegypti populations exhibit gene flow among different populations. Similarly, Ae. aegypti populations from Lahore also show intra-population gene flow (Table 3). Although, both cities present a distinct pat- tern of genetic structure, however, few areas from both the cities show that Ae. aegypti populations are genetically related to each other (Fig. 4). J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 550 http://jad.tums.ac.ir Published Online: October 04, 2016 Table 1. Localities from Lahore and Faisalabad selected for sampling of Ae. aegypti populations Sr. No. Samples from Faisalabad District Sr. No. Samples from Lahore District Sampling areas Code Sampling areas Code 1 Rescue office Jhang road F1 10 UC 73 (A) L1 2 Awami colony Tariq road F2 11 UC 81 (A) L2 3 AARI Jhang road F3 12 Ghawal mandi L3 4 Gulfishan colony F4 13 Mazang chungi area L4 5 Nazimabad F5 14 Qila Gujersing L5 6 Narwala road F6 15 BISE lawerence road L6 7 GCUF F7 16 UC69 L7 8 UAF Jail road F8 17 UC 73 (B) L8 9 NIAB F9 18 MAO college L9 Table 2. Primers used in RAPD analysis showing size and number of amplified fragments obtained (No.) from Ae. aegypti populations from Lahore and Faisalabad Sr. No. Primers Nucleotide sequence Size (bp) No. 1 A-03 5’- AGTCAGCCAC - 3’ 200–1500 12 2 A-04 5’- AATCGGGCTG - 3’ 350–1850 11 3 A-06 5’- GGTCCCTGAC - 3’ 375–1200 9 4 A-18 5’- AGGTGACCGT - 3’ 210–1300 8 Table 3. Summary of Nei’s Analysis of Gene diversity among between Aedes aegypti populations from Lahore and Faisalabad through RAPD markers Populations Genetic variation GST Gene flow Nm Heterozygosity Ht Lahore 0.131 3.317 0.294 Faisalabad 0.117 3.773 0.260 Lahore, Faisalabad 0.341 1.966 0.379 Table 4. RAPD Primer Sequence (Genelink) and their codes Sr. No. Primer name (code) Primer sequence GC contents 1 GL Decamer A-03 (A-03) AGTCAGCCAC 60% 2 GL Decamer A-04 (A-04) AATCGGGCTG 60% 3 GL Decamer A-05 (A-05) AGGGGTCTTG 60% 4 GL Decamer A-06 (A-06) GGTCCCTGAC 70% 5 GL Decamer A-08 (A-08) GTGACGTAGG 60% 6 GL Decamer A-10 (A-10) GTGATCGCAG 60% 7 GL Decamer A-11 (A-11) CAATCGCCGT 60% 8 GL Decamer A-12 (A-12) TCGGCGATAG 60% 9 GL Decamer A-18 (A-18) AGGTGACCGT 60% 10 GL Decamer A-19 (A-19) CAAACGTCGG 60% J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 551 http://jad.tums.ac.ir Published Online: October 04, 2016 Fig. 1. Map showing different collection sites from Lahore and Faisalabad (Faisalabad: 73°74 E and 30°31.5 N, Lahore: 74°01-74°39′ E and 31°15′-31°45′ N) Fig. 2. Amplification profile of Aedes aegypti populations with primers A-04 from Lahore and Faisalabad. (L1-9: Lahore samples, F1-9: Faisalabad samples, Lane M= DNA Marker (size 1kb) J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 552 http://jad.tums.ac.ir Published Online: October 04, 2016 Fig. 3. Percentage of polymorphic loci amplified by each RAPD primer Fig. 4. Dendrogram based on Nei’s genetic distances among Aedes aegypti populations from Lahore and Faisalabad. (Scale Bar: branch length, *F/L: Faisalabad/Lahore) F1: Rescue office Jhang road L1: UC 73 (A) F2: Awami colony Tariq road L2: UC 81 (A) F3: AARI Jhang road L3: Ghawal mandi F4: Gulfishan colony L4: Mazang chungi area F5: Nazimabad L5: Qila Gujersing F6: Narwala road L6: BISE lawerence road F7: GCUF L7: UC69 F8: UAF Jail road L8: UC 73 (B) F9: NIAB L9: MAO College J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 553 http://jad.tums.ac.ir Published Online: October 04, 2016 Discussion Aedes aegypti is one of the most important mosquito vectors of viral agents in world, due to its role in transmitting chikungunya virus (CHIKV) and dengue virus (DENV) (Thava- ra et al. 2001, Bosio et al. 2005). The genetic diversity has been reported in Ae. aegypti pop- ulations (Gorrochotegui-Escalante et al. 2000, Souza et al. 2001, Ayres et al. 2002, Franco et al. 2002). Many DNA techniques such as RAPD, microsatellites and mitochondrial DNA markers have been widely used to detect var- iation in geographically and genetically iso- lated populations (Ballinger-Crabtree et al. 1992, Apostol et al. 1994, Ayres et al. 2003, Paduan et al. 2006). In Pakistan, few studies have been conducted on barcoding of mos- quitoes through mitochondrial COI gene and genetic structure of Ae. aegypti populations using microsatellite markers (Ashfaq et al. 2013, Rasheed et al. 2013). However, still genetic analysis of Ae. aegypti populations has not been performed yet, particularly in the dengue outbreaks. Thus, the present study reinforced the significance of genetic analy- sis of Ae. aegypti populations in dengue af- fected areas in Pakistan. The level of genetic diversity was revealed among various mosquito species and multi- ple haplotypes were reported by Ashfaq et al. (2013) in Pakistan. However, they mainly fo- cused on biodiversity and distribution of mosquitoes in KPK and Punjab provinces. Furthermore, Rasheed et al. 2013 used 13 mi- crosatellite markers and found low genetic diversity in mosquitoes collected from 13 dif- ferent cities. Low level of genetic diversity with less differentiation was revealed and considered due to better road infrastructure. Although present results are consistent to their resultsو however, based on the ampli- fied fragments in our findings, it is suggest- ed that Ae. aegypti populations are genetical- ly more diverse as previously reported in Pakistan. In addition, the genetic similarity in populations at the border area of Faisalabad to few of Lahore populations supports the idea of transportations (Rasheed et al. 2013). The heterozygosity in Ae. aegypti popu- lations from Faisalabad and Lahore (Ht. 0.311) found similar to the results obtained from Brazil using 27 RAPD loci (Ht= 0.390) (Ayres et al. 2003) and to the average heter- ozygosity of 21 loci (Ht= 0.388) in Brazil, 57 loci in Puerto Rico (Ht= 0.354) (Apostol et al. 1996) and in 131 loci of 20 Ae. aegypti populations from Mexico through RAPD markers (Franco et al. 2002). Gorrochotegui- Escalante et al. (2000) also distinguished pop- ulations of Ae. aegypti from ten cities of the Mexican northeast coast, using 60 RAPD polymorphic loci. The gene flow and genetic variations were observed within populations separated within a range of 90 to 250 km. While, Lahore-Faisalabad Ae. aegypti pop- ulations correspond to ~100 km distance. Thus, the polymorphic pattern in our find- ings showed genetic variations by distance (55 out of 70 loci, data not shown). The genetic variation (GST = 0.341, Nm= 0.966) and polymorphism (51%) found in the present study were consistent to those reported in Argentina (GST = 0.249, Nm= 0.75) (Souza et al. 2001) and Brazil (GST = 0.317, Nm= 0.54, HS= 0.274) (Ayres et al. 2003). Santos et al. (2011) found high level of polymorphism and distinguished 52 mark- ers in the four populations through RAPD analysis, ranging in size from 300 to 2072 bp. Their percentage of polymorphic loci varied from 82.69 to 94.23. In contrast, we amplified loci ranged from 200 to 2000 bp with 22.5 to 51 percentage of polymorphism, which is closer to the genetic variability and 50% similarity in Ae. aegypti populations as demonstrated by Hiragi et al. (2009) using RAPD markers. However, current findings showed maximum similarity (0.8846) and genetic distance (0.045) among Ae. aegypti J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 554 http://jad.tums.ac.ir Published Online: October 04, 2016 populations (data not shown). A population differentiation from Faisalabad (GST = 0.131) and Lahore (GST = 0.117) with higher degree of gene flow (Nm= 3.317 Faisalabad, Nm= 3.773 Lahore) showed similarity to a low population differentiation (GST = 0.208) with higher degree of gene flow (Nm= 1.90) as described by Paduan et al. (2006) in Brazil indicating both intra- and inter-population genetic variability. Current findings agree to the study of De Sousa et al. (2001). They successfully used RAPD markers to estimate the level of poly- morphism and genetic relatedness in five populations of Ae. aegypti from Argentina and Puerto Rico. However, contrary to the present study, heterozygosity and genetic dis- tance values were found significant at all loci suggesting that the population from Puerto Rico was different from Argentina populations. Subsequently, de Patarro et al. (2013) assessed the genetic variability through RAPD mark- ers in five Ae. aegypti populations from Bra- zil and USA. A total of 165 polymorphic DNA loci were generated with six primers compared to 55 loci using ten primers in the present study. The mean value of inter-pop- ulation genetic diversity (GST) ranged from 0.044 to 0.289 suggesting that genetic varia- tion existed in Ae. aegypti populations. Simi- larly, we found genetic diversity (GST) ranged from 0.260 in Faisalabad to 0.294 in Lahore. Ocampo and Wesson (2004) found genetic variation in Ae. aegypti populations in Co- lombia. The genetic heterogeneity ranged from 19% to 60% revealed through RAPD markers closer to our findings. In addition, the genetic distance and migration rate (Nm) demonstrated a moderate genetic differentia- tion among Ae. aegypti populations which is consistent to the current findings of genetic differentiation in Ae. aegypti populations ex- isted in Faisalabad and Lahore. Despite of RAPD, other DNA markers are also being used to analyse genetically mosquitoes. Paupy et al. (2012) examined genetic structure and geographical origin of Ae. aegypti in Bolivia and found polymor- phism at nine microsatellite loci and in two mitochondrial DNA regions (COI and ND4). Similar to the present work, two genetic lin- eages were obtained after genetic analysis of mtDNA sequences. One lineage was restrict- ed to rural localities and the other was found closely related to African Ae. aegypti pop- ulations. Similarly, Caldera et al. (2013) ge- netically analyzed two populations of Ae. aegypti from dengue outbreak areas in Co- lombia using mitochondrial gene ND4. The genetic differentiation with a limited gene flow (Nm= 1.40) indicated that both the pop- ulations of Ae. aegypti are genetically dis- tinct. Hlaing et al. (2010) also found genetic heterogeneity in Ae. aegypti populations us- ing microsatellite markers in Myanmar, Cam- bodia, Thailand, Sri Lanka and Nigeria. The genetic variation even in a distance of 500 m and genetic similarity in distant locations due to dispersal were observed. In parallel, Costa-da-Silva (2005) also found a distinct genetic pattern in Ae. aegypti populations in Peru using mitochondrial DNA. Further- more, Yáñez et al. 2013 tested 51 samples of Ae. aegypti populations from Peru using mi- tochondrial ND4 gene and identified five haplotypes of Ae. aegypti grouped in two lineages. This variability occurred mainly due to migration of this vector and partly through human activity. Thus, transportation has been considered as a main source of spread of mosquito borne diseases, however, other fac- tors like use of insecticides and elimination of larval habitats in and around dwellings had a demonstrable impact on gene flow and genetic structure of Ae. aegypti populations (Wallis et al. 1984, Failloux et al. 1995, Paupy et al. 2000, Lerdthusnee and Chareonviri- yaphap 2002) which needs to be further in- vestigated on larger scale using gene specific molecular markers. J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 555 http://jad.tums.ac.ir Published Online: October 04, 2016 Conclusion The genetic relatedness and Nm (number of migrants) value show that Ae. aegypti pop- ulations exhibit intra-population gene flow between Faisalabad and Lahore. Although, both cities show a distinct pattern of genetic structure, however, few areas from both the cities show genetic similarity. It is also con- cluded that Ae. aegypti populations are ge- netically more diverse as previously reported in Pakistan. Acknowledgements The facilities provided by Department of Zoology, Government College University Fai- salabad (GCUF) and the financial support from HEC Project are highly acknowledged to con- duct this research work. The authors are also thankful to Punjab Health Department, Govt. of Punjab for mosquito samples and for their cooperation during the research period. The au- thors declare that there is no conflict of interest. References Apostol BL, Black 4th WC, Reiter P, Miller BR (1996) Population genetics with RAPD-PCR markers: the breeding structure of Aedes aegypti in Puerto Rico. Heredity. 76: 325–334. Ashfaq M, Hebert PDN, Mirza JH, Khan AM, Zafar Y, Mirza MS (2014) Ana- lyzing Mosquito (Diptera: Culicidae) Diversity in Pakistan by DNA Bar- coding. PLoS ONE. 9(5): 1–12. Ashraf HM (2013) Studies on the comparative genetic analysis of Aedes aegypti pop- ulations from Lahore and Faisalabad using molecular DNA markers. MPhil thesis, Department of Zoology, Wild- life and Fisheries, Government Col- lege University Faisalabad, Pakistan. Ayres CFJ, Melo-Santos MAV, Sole-Cava AM, Furtado AF (2003) Genetic dif- ferentiation of Aedes aegypti (Diptera: Culicidae), the major dengue vector in Brazil. J Med Entomol. 40: 430–435. Ayres CFJ, Romao TPA, Melo-Santos MAV, Furtado AF (2002) Genetic diversity in Brazilian populations of Aedes al- bopictus. Mem Inst Oswaldo Cruz. 97(6): 871–875. Ballinger-Crabtree ME, Black WC, Miller B.R (1992) Use of genetic polymor- phisms detected by the random-am- plified polymorphic DNA polymerase chain reaction (RAPD-PCR) for dif- ferentiation and identification of Ae- des aegypti subspecies and popula- tions. Am J Trop Med Hyg. 47: 893– 901. Beebe NW, Whelan PI, Hurk AVD, Ritchie SA, Cooper RD (2005) Genetic diver- sity of the dengue vector Aedes ae- gypti in Australia and implications for future surveillance and mainland in- cursion monitoring. Commun Dis In- tell. 29: 299–304. Bibi M, Zahoor MK, Zahoor MA, Ashraf HM, Majeed HN, Nasir S, Rasool B, 2015 Genetic analysis of mosquitoes from rural and urban areas of Sialkot, Paki- stan. Int J Agric Biol. 17(4): 809-814. Bosio CF, Harrington LC, Jones JW, Sithiprasasna R, Norris DE, Scott TW (2005) Genetic structure of Aedes ae- gypti populations in Thailand using mitochondrial DNA. Am J Trop Med Hyg. 72: 434–442. Bracco J, Capurro M, Lourenço-DE-Oliveira R, Sallum MA (2007) Genetic varia- bility of Aedes aegypti in the Ameri- cas using a mitochondrial gene: evi- dence of multiple introductions. Mem Inst Oswaldo Cruz. 102: 573–580. Caldera SM, Jaramillo MC, Cochero S, Pé- rez-Doria A, Bejarano EE (2013) Ge- netic differences between populations of Aedes aegypti from municipalities J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 556 http://jad.tums.ac.ir Published Online: October 04, 2016 in northern Colombia, with high and low dengue incidence. Biomédica. 33: 89-98. Cologna R, Armstrong PM, Rico-Hesse R (2005) Selection for virulent dengue viruses occurs in humans and mos- quitoes. J Virol. 79: 853–859. Costa-da-Silva, ALD, Capurro ML, Bracco JE (2005) Genetic lineages in the yel- low fever mosquito Aedes (Stegomyia) aegypti (Diptera: Culicidae) from Pe- ru. Memórias do Instituto Oswaldo Cruz, 100(6): 539-544. de Patarro TF, Guirado MM, Ravazzi LM, de Bicudo H (2013) Genetic structure of Aedes aegypti populations deter- mined using pairwise comparisons. Genet Mol Res. 12(3): 3775-3787. De Sousa GB, Blanco A, Gardenal CN (2001) Genetic relationships among Aedes aegypti (Diptera: Culicidae) popula- tions from Argentina using random amplified polymorphic DNA poly- merase chain reaction markers. J Med Entomol. 38(3): 371-375. Failloux A-B, Darius H, Pasteur N (1995) Ge- netic differentiation of Aedes aegypti, the vector of dengue virus in French Polynesia. J Am Mosq Contr Assoc. 11: 457–462. Fatima Z, Idrees M, Bajwa MA, Tahir Z, Ullah O, Zia MQ, Hussain A, Akram M, Khubaib B, Afzal S, Munir S, Saleem S, Rauff B, Badar S, Naudhani M, Butt S, Aftab M, Ali L, Ali M (2011) Serotype and genotype analysis of dengue virus by sequenc- ing followed by phylogenetic analy- sis using samples from three mini out- breaks 2007–2009 in Pakistan. BMC Microbiol. 11: 200. Florencio M, Díaz-Paniagua C, Gómez- Rodríguez C, Serrano L (2014) Bio- diversity patterns in a macroinverte- brate community of a temporary pond network. Insect Conserv. Diversity. 7: 4–21. Fradin MS, Day JF (2002) Comparative effi- cacy of insect repellents against mos- quito bite. N Eng J Med. 347(1): 13–18. Fraga EC, Santos JMM, Maia JF (2003) En- zymatic variability in Aedes aegypti (Diptera: Culicidae) populations from Manaus- AM, Brazil. Genet Mol Bi- ol. 26: 181–187. Franco FG, Munoz MDL, Fuentes SL, Salas LF, Rejon JG, Beaty BJ, Black WC (2002) Large genetic distances among Aedes aegypti populations along the south pacific coast of Mexico. AM J Trop Med Hyg. 6(5): 594–598. Fulton, RE, Salasek ML, DuTeau NM, Black WC (2001) SSCP analysis of cDNA markers provides a dense link- age map of the Aedes aegypti ge- nome. Genetics. 158: 715–726. Gibbons RV, Vaughn DW (2002) Dengue: an escalating problem. BMJ. 324: 1563–1566. Gilani I (2012) 2,340 malaria cases reported. Dawn News, June 04. Available at: http://dawn.com/news/407357/malaria- plagues-faisalabad. Gorrochotegui-Escalant N, Munoz ML, Fer- nandez-Salas I, Beaty BJ and Black BJ (2000) Genetic isolation by dis- tance among Aedes aegypti popula- tions along the northeastern coast of Mexico. Am J Trop Med Hyg. 62: 200–209. Gorrochotegui-Escalante N, Gomez-Machorro C, Lozano- Fuentes S, Fernandez- Salas I, Munoz ML, Farlan-Ale JA, Garcia-Rejon J, Beaty BJ, Black WC (2002) Breeding structure of Aedes aegypti populations in Mexico varies by region. Am J Trop Med Hyg. 66: 213–222. Herrel N, Amerasinghe FP, Ensink J, Mukhtar M, Van Der Hoek W, Kon- radsen F (2001) Breeding of Anoph- eles mosquitoes in irrigated areas of J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 557 http://jad.tums.ac.ir Published Online: October 04, 2016 South Punjab, Pakistan. Med Vet En- tomol. 15: 236–248. Herrera F, Urdaneta L, Rivero J, Zoghbi N, Ruiz J, Carrasquel G, Martínez JA, Pernalete M, Villegas P, Montoya A, Rubio-Pais Y, Rojas E (2006) Popu- lation genetic structure of the dengue mosquito Aedes aegypti in Venezuela. Mem Inst Oswaldo Cruz. 101: 625– 633. Hiragi C, Simoes K, Martins E, Queiroz P, Lima L, Monnerat R (2009) Genetic variability in Aedes aegypti (L.) (Dip- tera: Culicidae) populations using RAPD markers. Neotrop Entomol. 38: 542–547. Hlaing T, Tun-Lin W, Somboon P, Socheat D, Setha T, Min S, Thaung S, Anyaele O, Silva BD, Chang MS, Prakash A, Linton Y, Walton C (2010) Spatial genetic structure of Aedes aegypti mos- quitoes in mainland Southeast Asia. Evol Applic. 3(4): 319–339. Huber K, Le Loan L, Hoang TH, Ravel S, Rodhain F, Failloux AB (2002) Ge- netic differentiation of the dengue vector, Aedes aegypti (Ho Chi Minh City, Vietnam) using microsatellite markers. Mol Ecol. 11: 1629–1635. Idrees S, Ashfaq UA (2102) A brief review on dengue molecular virology, diag- nosis, treatment and prevalence in Pakistan. Gen Vacc Ther. 10: 6. Itrat A, Khan A, Javaid S, Kamal M, Khan H, Javed S, Kalia S, Khan AH, Sethi MI, Jehan I (2008) Knowledge, aware- ness and practices regarding dengue fever among the adult population of dengue hit cosmopolitan. PLoS ONE. 3(7): 1–6. Julio NB, Chiappero MB, Rossi HJ, Rondan Duenas JCR, Gardenal CN (2009) Genetic structure of Aedes aegypti in the city of Córdoba (Argentina), a re- cently reinfested area. Mem Inst Os- waldo Cruz. Rio de Janeiro. 104(4): 626–631. Khan E, Kisat M, Khan N, Nasir A, Ayub S, Hasan R (2010) Demographic and clinical features of dengue fever in Pakistan from 2003–2007: A retro- spective cross-sectional study. PLoS ONE. 5(9): 1–6. Lerdthusnee K, Chareonviriyaphap T (2002) Genetic differentiation of Aedes ae- gypti mainland and island popula- tions from Southern Thailand. J Am Mosq Contr Assoc. 18: 173–177. Morlais I, Severson DW (2003) Intraspecific DNA variation in nuclear genes of the mosquito Aedes aegypti. Insect Mol Biol. 12: 631–639. Naeem-Ullah U, Akram W, Suhail A, Rana SA (2010) Grouping of different mosquito species on the bases of lar- val habitats. Pak J Agri Sci. 47(2): 124–131. Nei M (1978) Estimation of average hetero- zygosity and genetic distance from a small number of individuals. Genet- ics. 89: 583–590. Nikookar SH, Moosa-Kazemi SH, Oshaghi MA, Vatandoost H, Kianinasab A (2010) Species composition and di- versity of mosquitoes in Neka County, Mazandaran Province, northern Iran. J Arthropod-Borne Dis. 4(2): 26–34. Ocampo CB, Wesson DM (2004) Population dynamics of Aedes aegypti from a dengue hyperendemic urban setting in Colombia. Am J Trop Med Hyg. 71(4): 506-513. Paduan KS, Araújo-Júnior JP, Ribolla PEM (2006) Genetic variability in geograph- ical populations of Aedes aegypti (Dip- tera, Culicidae) in Brazil elucidated by molecular markers. Genet Mol Bi- ol. 29: 391–395. Paupy C, Vazeille-Falcoz M, Mousson L, Rodhain F, Failloux AB (2000) Aedes aegypti in Tahiti and Moorea (French Polynesia): isozyme differentiation in J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 558 http://jad.tums.ac.ir Published Online: October 04, 2016 the mosquito population according to human population density. Am J Trop Med Hyg. 62: 217–224. Paupy C, Chantha N, Reynes JM, Failloux AB (2005) Factors influencing the population structure of Aedes aegypti from the main cities in Cambodia. Heredity. 95: 144–147. Paupy C, Le Goff G, Brengues C, Guerra M, Revollo J, Simon ZB, Herve JP, Fon- tenille D (2012) Genetic structure and phylogeography of Aedes aegypti, the dengue and yellow-fever mosquito vec- tor in Bolivia. Infect Genet Evol. 12(6): 1260–1269. Pinheiro FP, Corber SJ (1997) Global situa- tion of dengue and dengue haemor- rhagic fever, and its emergence in the Americas. WHO Stat Q. 50: 161–169. Raheel U, Faheem M, Riaz MN, Kanwal N, Javed F, Zaidi NS, Qadri I (2011) Dengue fever in the Indian subconti- nent: an overview. J Infect Dev Ctries. 5(4): 239–247. Rasheed SB, Boots M, Frantz AC, Butlin RK (2013) Population structure of the mosquito Aedes aegypti (Stego- myia aegypti) in Pakistan. Med Vet Entomol. 27: 430–440. Ravel S, Herve JP, Diarrassouba S, Kone A, Cuny G (2002) Microsatellite mark- ers for population genetic studies in Aedes aegypti (Diptera: Culicidae) from Côte d’Ivoire: evidence for a micro- geographic genetic differentiation of mosquitoes from Bouaké. Acta Trop. 82: 39–49. Sajid A, Ikram A, Ahmed M (2012) Dengue fever outbreak 2011: clinical profile of children presenting at Madina teach- ing hospital Faisalabad. J Univ Med Dent Coll. 3 (1): 42–47. Santos JMM, Fraga EC, Maia JF, Tadei WP (2011) Genetic Diversity in Dengue Mosquito, Aedes aegypti (Diptera: Cu- licidae) from Amazon Region: Com- parative Analysis with Isozymes and RAPD Loci. The Open Trop Med J. 4: 11–20. Scarpassa VM, Cardoza TB, Cardoso-Júnior RP (2008) Population genetics and phylogeography of Aedes aegypti (Dip- tera, Culicidae) from Brazil. Am J Trop Med Hyg. 78: 895–903. Severson DW, Zaitlin D, Kassner VA (1999) Targeted identification of markers linked to malaria and filarioid nem- atode parasite resistance genes in the mosquito Aedes aegypti. Genet Res. 73: 217–224. Severson DW, Meece JK, Lovin DD, Saha G, Morlais I (2002) Linkage map or- ganization of expressed sequence tags and sequence tagged sites in the mos- quito, Aedes aegypti. Insect Mol Biol. 11: 371–378. Shortall CR, Moore A, Smith E, Hall MJ, Woiwod IP, Harrington R (2009) Long-term changes in the abundance of flying insects. Insect Conserv Di- vers. 2: 251–260. Siddiqui FJ, Haider S, Bhutta ZA (2009) En- demic Dengue Fever: a seldom rec- ognized hazard for Pakistani children. J Infect Dev Ctrie. 3(4): 306–312. Souza GB de, Blanco A, Gardenal CN (2001) Genetic relationships among Aedes ae- gypti (Diptera: Culicidae) populations from Argentina using random ampli- fied polymorphic DNA polymerase chain reaction markers. J Med En- tomol. 38: 371–375. Thavara U, Tawatsin A, Chansang C, Kong- ngamsuk W, Paosriwong S, Boon- Long J, Rongsriyam Y, Komalamisra N (2001) Larval occurrence, ovipo- sition behavior and biting activity of potential mosquito vectors of dengue on Samui Island, Thailand. J Vector Ecol. 26: 172–180. Wallis GP, Tabachnick WJ, Powell J (1984) Genetic heterogeneity among Carib- J Arthropod-Borne Dis, December 2016, 10(4): 546–559 HM Ashraf et al.: Genetic Analysis of … 559 http://jad.tums.ac.ir Published Online: October 04, 2016 bean populations of Aedes aegypti. Am J Trop Med Hyg. 33: 492–498. WHO (1997) Dengue hemorrhagic fever, diagnosis, treatment and control, 2nd ed. World Health Organization, Ge- neva, Switzerland. WHO (2009) Den- gue guidelines for diagnosis, treat- ment, prevention and control. New edition. Geneva, Switzerland. WHO (2013) Dengue Fever in Pakistan. Weekly Epidem. Monit. Available at: http://reliefweb.int/report/pakistan/who -emro-weekly-epidemiological-moni- tor-volume-volume-6-issue-37-dengue- fever-0. WHO (2014) Dengue and severe dengue. Fact sheet No117. Available at: http:// www. who.int/mediacentre/factsheets/ fs117/ en/#. Yáñez P, Mamani E, Valle J, García MP, León W, Villaseca P, Torres D, Cab- ezas C (2013) Genetic variability of Aedes Aegypti determined by mito- chondrial gene ND4 analysis in elev- en endemic areas for dengue in Peru. Rev Peru Med Exp Salud Públi- ca. 30(2): 246–250. Zahoor MK, Suhail A, Zahoor S, Iqbal A, Awan FS (2013) Molecular charac- terization of Scarab beetles (Scarabae- idae: Coleoptera) using RAPD mark- ers. Pak J Life Soc Sci. 11: 238–243.