- 1Department of Biological Science, College of Science, University of Jeddah, Jeddah, Saudi Arabia
- 2Department of Biology Sciences, Faculty of Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
- 3Department of Biology, College of Science, Jazan University, Jazan, Saudi Arabia
- 4Department of Biology, College of Science, Qassim University, Buraydah, Saudi Arabia
Introduction: The Aedes aegypti constitutes the primary vector for dengue fever, yellow fever, chikungunya, Zika, West Nile, and encephalitis viruses, all of which have impacted One Health (human, animal, and environmental health) significantly. It has been distributed widely in urban settings in Saudi Arabia, particularly in cities like Jeddah and Jizan, a situation that underscores the urgent need for innovative and sustainable vector control strategies. Molecular tools, such as DNA barcoding using mitochondrial markers, have become essential for identifying mosquito species accurately and understanding their role in disease transmission. Such knowledge is vital for informing effective, climate-resilient public health interventions.
Methods: This research focuses on identifying Aedes species in various regions of Saudi Arabia using polymerase chain reaction (PCR) and sequencing techniques, in order to evaluate the molecular diversity of these dengue vectors in Jeddah and Jizan. The study utilizes the cytochrome one oxidase (COI) gene as a molecular marker for phylogenetic analysis to compare the populations of Aedes species.
Results: The findings reveal the presence of significant genetic variation among mosquito populations. In the Jeddah region, the Ae. aegypti types MN299016.1 and KU495081.1 match completely (100%) the respective populations found in Argentina and Australia, with 93.1% (27/29) and 6.9% (2/29) respectively. Meanwhile, the samples from the Jizan region are 100% and 99.6% similar to the Ae. aegypti types MN298998.1, MK300226.1, PP892777.1, and MF043259.1 found in Canada, Kenya, India, and England.
Conclusion: This study underscores the necessity of using molecular techniques in vector surveillance to mitigate the spread of Zoonotic and vector borne diseases in Saudi Arabia. Moreover, these efforts align with the objectives of the Saudi Vision 2030 by promoting environmentally responsive vector surveillance in Jeddah and Jizan, thereby supporting integrated approaches to public health and ecological sustainability.
1 Introduction
The noxious Ae. aegypti inflicts a plethora of vector-borne diseases, with yellow fever, chikungunya, Zika, West Nile, encephalitis viruses, and dengue fever responsible for significant morbidity and mortality in developing countries (1–3). These viral haemorrhagic fevers are prevalent where the urban-dwelling mosquito Ae. aegypti is incriminated as a major culprit (4, 5). There are several factors that contribute to the proliferation of mosquitoes and the spread of diseases in Jeddah and Jazan, including temperature, coastal fog, and rainfall (6, 7). In Saudi Arabia, some studies indicate a relationship between climatic conditions and the increase in the number of mosquitoes and their ability to transmit diseases (8), recorded a high density of mosquitoes in Jeddah during the winter at temperatures below 28°C and a lower density in the summer at temperatures of 40°C.
In recent years, some Saudi Arabian cities have seen an epidemic of dengue fever, which are endemic in certain cities of Saudi Arabia, such as Jeddah and Jizan (9). However, there is currently no comprehensive and authentic record of the Aedes mosquito species of Jeddah and Jizan, Saudi Arabia.
Detailed reports on mosquito species in Saudi Arabian are available in the literature (10–14). Understanding the diversity and distribution of these species is crucial for research on biodiversity and disease management. It is important to comprehend their habitats and the pathogens they transmit to develop effective vector control strategies aimed at preventing disease transmission (15, 16). Globally, only 10% of mosquito species have been identified, primarily due to a lack of taxonomic expertise (17, 18). The mitochondrial cytochrome c oxidase subunit 1 (COI) gene has demonstrated effective in the rapid and accurate identification of Aedes aegypti mosquitoes and their taxonomic features (19). The COI gene is particularly valuable for evolutionary studies (20). Molecular techniques are preferred for identifying species, assessing genetic diversity, and studying molecular evolution because they provide high resolution. These methods do not depend on the sex or age of mosquitoes to produce reliable results (21–23).
According to research conducted recently on insect taxonomy and phylogenetic analysis, DNA barcoding is an effective and convenient method for identifying various genera and species (24–27). In order to identify mosquito species accurately, taxonomically molecular biology using polymerase chain reaction (PCR) was employed by recent researchers (28). Interest in studying mosquitoes from a genetic perspective began in the mid-1950s, especially concerning the species that transmit epidemic diseases, notably those that belong to the genera Aedes, Anopheles, and Culex (29). The present study investigates the genetic diversity of Ae. aegypti, and the vector of dengue fever in the cities of Jazan and Jeddah in Saudi Arabia, using the cytochrome one oxidase (COI) gene as a significant marker. The aim of the study is to examine the COI sequences from random samples of adult Ae. aegypti collected from Jazan and Jeddah, in order to determine the similarities/diversity in their DNA sequences, and to compare them with those of the same species from other regions around the world.
2 Methods and materials
2.1 DNA extraction
All genomic DNA extraction experiments were conducted in molecular identification experimental unit, in the Dengue Mosquito Experimental Station (DMES), belonging to the Department of Biological Sciences, Faculty of Sciences, King Abdul-Aziz University, Jeddah, Saudi Arabia. As (30) explained, DNA can be extracted from preserved mosquito specimens by immersing the entire body of each specimen in absolute ethanol. The insect lysis buffer involved is prepared using the following chemical components: 16.5 grams of GuSCN (Sigma), 12ml of 0.5 M EDTA pH 8.0, 6ml of 1 M Tris-HCL pH 8.0, 1ml of Triton X-100 (Sigma), and 10ml Tween-20 (Fluka). In order to conduct the experiment in the present study, mosquito samples were placed in Eppendorf tubes, and a combination of proteinase K and insect lysis buffer was added to each tube. After vortexing and incubating the samples at 56 °C for 10 minutes, the samples were centrifuged before the supernatant was transferred to a tube for further analysis. The next step was to gently mix the supernatant with 100 litres of ethanol, and then to add the ethanol to the supernatant. Thereafter, the mixture was placed onto a silica gel-based spin column for further separation. In order to remove any possible remaining liquid, we spun the solution at 12,000 RPM for two minutes at the highest speed. The liquid was removed from the collection tube by hand.
In the same experiment, 500 litres of 70% ethanol were used in the process, and the results were the same. Afterwards, the samples were centrifuged again at 12,000 RPM on the spin column for two more minutes before they were removed from the centrifuge for further processing. It was then rewashed with ethanol to remove any remnants of the liquid that remained. To ensure that the spin column had dried completely when the washing was finished, we centrifuged it for three minutes at 13,000 RPM. In order to move the upper part of the spin column to the Eppendorf tube, the upper part of the column was removed. During the release of the DNA that was trapped within the silica gel, it was first necessary to add 20 litres of nuclease water to the centre of the gel and to leave it for about 10 minutes, so that the DNA could be released. A one-minute spin at 12,000 RPM was then applied to the sample, and it was then dried. This was repeated after another 50 litres of nuclease water was added to the solution. The DNA was then eluted from the silica gel and transferred to the Eppendorf tube, after the gel had been turned (30). Using Thermo Fisher Scientific’s NanoDrop 7000 instrument from Waltham, MA, USA, we assessed the quantity and quality of the DNA samples using the A260/Ratio, which was calculated using the NanoDrop 7000. In addition, an electrophoresis of DNA in 1% gel of agarose was conducted to assure the integrity of the DNA. Following the extraction and preservation of the DNA at a temperature of 20 °C for a short period of time, it was then stored at a temperature of -20 °C for a period, after which it could be used for experiments following preservation.
2.2 Agarose gel preparation
A microwaveable flask was filled with 70mL 1xTBE buffer, into which was dissolved 0.5% of agarose. The agarose was boiled for 1–2 minutes by swirling the flask occasionally until the agarose was completely dissolved. In order to accomplish this process, ethidium bromide (EtBr) was added to the agarose solution at a final concentration of approximately 0.2-5g/mL after the solvent had warmed up to around 50 °C for about 2–3 minutes. As the gel melted, it was poured into the tray of mini-gel, with the well comb positioned slowly in the middle, in order to avoid causing bubbles in the gel. After 20–30 minutes, once the gel had completely solidified, the comb was removed carefully, and the gel was prepared for use. In addition to covering the gel in the gel box, 1xTBE buffer was also used to load the DNA into each well. One Kbp DNA ladder was also added to each well.
2.3 Gene amplification
The mitochondrial cytochrome oxidase subunit I (COI) gene was utilized to identify polymorphism among the different populations used in this study. Gene amplification was preformed using a PCR with the universal barcoding primer set, LCO1490 (5’-GGTCAACAAATCATAAAGATATTGG-3’) and HCO2198 (5’-TAAACTTCAGGGTGACCAAAAAATC-3’) (31), generating amplicons ranging between 600–700 base pairs. Each 25 μl reaction contained 15.3 μl of 1× bovine serum albumin (BSA), 2.5 μl of 10× Reaction Buffer, 2 μl of dNTPs (2.5 μM), 1.25μl of each primer (10 μM/L), 0.2 μl Taq DNA Polymerase (1.0 U), and 2.5 μl of template DNA (>50 ng). The PCR conditions consisted of an initial denaturation at 94 °C for two minutes, followed by 40 cycles of denaturation at 94 °C for 30 seconds, annealing at 49 °C for 45 seconds and extension at 72 °C for 45 seconds, with a final extraction step at 72 °C for one minute. The PCR products were verified by electrophoresis on 1.5% agarose gel and the gel was then analysed under UV light. The amplified bands were subsequently purified from the gel and submitted for sequencing to Macrogen Company, Korea, using an applied Biosystems 3730Xl DNA Analyzer.
2.4 Sequencing
For sequencing, the samples were sent to a Macrogen company in Seoul, South Korea. A QIA quick PCR purification kit (Promega, Madison, WI, USA) was used to purify the amplified DNA for sequencing, and a Big Dye terminator cycle sequencing Ready Reaction kit (Applied Biosystems, Forster City, CA, USA) was used to perform the sequencing on an ABI Prism® 310 Genetic Analyzer (Applied Biosystems) to measure the sequence of the amplified DNA.
2.5 Sequence and phylogenetic analysis
In order to compare the nucleotide information that was extracted using the sequences discussed above within the Gene Bank database, we used the Basic Local Alignment Search Tool (BLAST) (https://blast.ncbi.nlm.nih.gov/Blast.cgi), which can be found at http://blast.ncbi.nlm.nih.gov/Blast.cgi. A nucleotide sequence alignment was conducted by utilizing the Clustal W tool in MEGA to align the raw nucleotide sequences. In order to construct the phylogenetic tree of each locus, a multiple sequence alignment was conducted using MEGA11 software18, and the sequences were then matched with those in the GenBank database of NCBI (http://www.ncbi.nlm.nih.gov), and the reference sequences were obtained to perform a phylogenetic analysis on the sequences retrieved. A phylogenetic tree of the populations was constructed based on the dissimilarity values computed using the Maximum-Likelihood method with the Tamura-Nei Model2. The distance matrix determined via this approach was utilized to infer the evolutionary relationships across the populations. In addition to the bootstrap analysis, a maximum of 1,000 replicates of the machine learning analysis was also conducted.
3 Results and discussion
3.1 Molecular identification
The DNA samples employed for this study were collected from two different regions, Jeddah and Jazan, and consisted of 29 samples from each region, with a total of 58 samples for both regions. By using the same protocol and procedure as that employed by (30), namely lysis buffers, DNA was extracted from the insects and analysed on a 1% agarose gel (Figure 1). As a preliminary step towards identification of the gene encoding, Cytochrome c oxidase subunit, I (COX1), genomic DNA was extracted from 58 randomly selected mosquito specimens taken from the two geographical regions in question. As a result of the PCR amplification using universal primers targeting the mTCOX-1 gene, 600–700 bp fragments were generated, which were then subjected to 1-5% agarose gel electrophoresis and then analysed (Figures 1, 2).

Figure 2. PCR amplification of the mT-COI gene was achieved in 48 mosquito samples. JE, Jeddah; JA, Jazan.
3.2 Gene sequencing and phylogenetic tree analysis
Phylogenetic trees (Tables 1, 2) show the convergence and divergence of the samples (Figure 2). Correlations were made between the samples and the gene bank samples (Table 1). In the results obtained from the Jeddah region, the Ae. aegypti types MN299016.1 and KU495081.1 demonstrated a complete match (100%) with the respective populations found in Argentina and Australia, with 93.1% (27/29) and 6.9% (2/29), respectively (Table 1; Figure 3). Meanwhile, in the samples from the Jazan region, the samples were 100% and 99.6% like the Ae. aegypti types MN298998.1, MK300226.1, PP892777.1, and MF043259.1 found in Canada, Kenya, India, and England (Table 1; Figure 3), with proportions of 58.6% (17/29), 20.7% (6/29), 13.8% (4/29), and 6.9% (2/29), respectively. As shown in Figure 3, the evolutionary distances between the taxa examined were compared with each other, and with the sequences related to the NCBI, in order to determine where the samples evolved from each other. A phylogenetic tree created with MEGA11 represented the convergence and divergence among the samples analysed. The phylogenetic analysis of the mTCOX-1 gene revealed that the species from Jeddah and Jazan did not match each other; instead, they formed two different clusters, according to the results of the analysis, which showed that the Ae. aegypti species formed two major groups.

Figure 3. Phylogenetic analysis of the study samples was performed using the Maximum Likelihood (ML) method, as implemented in the MEGA 11 program.

Table 2. The molecular identification of new collected Ae. aegypti from the Jeddah and Jazan regions and their accession numbers.
This study found that samples from Jazan demonstrated genetic similarity with those from Canada, Kenya, England, and India, which constituted the first cluster of samples. As a result of a further subdivision within this cluster, it was observed that there were two distinct sections, the first of which contained only the samples from Canada and Kenya, and the second contained two elements. In the first clade, the samples that matched those found in England were exclusively located, while the second clade contained samples that matched those found in India. As a result of the differentiation between the two clades, subclades were identified. The samples collected from the city of Jeddah formed part of the second cluster, as they demonstrated a genetic similarity with those from Argentina and Australia. In addition to the subdivision within this cluster, two distinct sections were identified within the cluster. Those sampled in the first section could only be compared with those from Argentina, whereas those from the second section were comparable with those from Australia. Several genetic mutations were observed in the study samples, which could be considered the cause of the differences observed between them.
This study aimed to molecularly identify Ae. aegypti mosquitoes from two distinct regions in Saudi Arabia: Jeddah and Jazan. By utilizing the mt COX1 gene as a marker, phylogenetic and sequence analyses revealed significant genetic variability between the mosquito populations from these regions. The results contributed valuable insights into the genetic structure of Ae. aegypti, providing data that may influence vector control programmes, given the importance of genetic diversity in disease transmission and control strategies.
3.3 Molecular identification of Ae. aegypti
The molecular identification of mosquito samples using the COX1 gene marker conducted by this study proved effective for distinguishing Ae. aegypti populations in Jeddah and Jazan. In both regions, all of the samples were identified as Ae. aegypti, confirming the species’ presence in these areas, as expected, due to its widespread distribution in tropical and subtropical regions (32, 33). The genetic variation observed between the two regions underscored the species’ adaptability to different environmental conditions.
The DNA extraction and PCR amplification conducted yielded fragments of 600–700 bp, corresponding to the COX1 gene, which was sequenced successfully. The sequences obtained were compared with those available in the NCBI Gene Bank, confirming the samples’ species as Ae. aegypti (Tables 1, 2). These findings corroborated previous reports regarding the utility of the COX1 gene as a molecular marker for species-level identification in mosquitoes (30).
3.4 Phylogenetic analysis and genetic divergence
A phylogenetic analysis using the Maximum Likelihood (ML) method was conducted, which revealed that there were two distinct genetic clusters; these corresponded to the mosquito populations from Jeddah and Jazan (Figure 3). This genetic divergence could be attributed to geographic isolation or to differences in selective pressures between the two regions, such as environmental conditions and human activity. The samples from Jeddah showed a complete match (100%) with Ae. aegypti strains from Argentina and Australia (MN299016.1 and KU495081.1, respectively). These samples formed a separate cluster in the phylogenetic tree, suggesting that the Jeddah population shared a closer genetic relationship with the Ae. aegypti populations from these regions. This may indicate a common origin or recent introduction, potentially due to human-mediated movement, such as international trade or travel (32). In contrast, the samples from Jazan displayed a more complex pattern of genetic similarity, matching Ae. aegypti strains from Canada, Kenya, India, and England (MN298998.1, MK300226.1, PP892777.1, and MF043259.1, respectively). This population formed two distinct sub-clusters, suggesting greater genetic diversity within the Jazan population. This may be the result of genetic introgression from multiple sources, or historical population expansion and contraction in response to environmental changes (34).
3.5 Implications of genetic divergence
The genetic divergence observed between the Jeddah and Jazan populations has important implications for mosquito control programmes. Genetic variation within mosquito populations can affect their susceptibility to insecticides and their ability to transmit diseases, such as dengue, Zika, and chikungunya. For example, populations with distinct genetic backgrounds may respond differently to chemical interventions, requiring tailored approaches for effective control (33). The divergence between the Jazan and Jeddah clusters suggested that control strategies may be required to be region-specific.
Although the present study did not analyse specific insecticide resistance genes, previous studies have documented resistance associated kdr mutations in Ae. aegypti populations from Jeddah and other regions of Saudi Arabia (35–37). These mutations, especially those in the voltage-gated sodium channel (Vssc) gene, have been strongly linked to pyrethroid resistance. A broader regional context displays the emergence of “super-resistance” due to combined kdr mutations in parts of Asia, posing serious challenges for vector control (38). These findings underscore the importance of incorporating insecticide resistance markers into future genetic analyses. This will help align population dynamic with resistance patterns, enhancing the accuracy and impact of region-specific vector control efforts.
For instance, the Jazan population’s genetic similarity to Ae. aegypti from geographically distant regions, such as Canada and Kenya may reflect a broader genetic pool that could influence vector competence or insecticide resistance. This highlighted the need for ongoing surveillance and molecular monitoring to track the genetic changes in mosquito populations and the need to adjust control strategies accordingly (39).
3.6 Geographic and environmental influence
The differences in genetic clustering between the Jeddah and Jazan populations may also reflect the geographic and environmental conditions of these regions. Jeddah, which is a coastal city, experiences higher humidity and temperature levels than Jazan, which may influence selective pressures on mosquito populations. Environmental factors, such as temperature, humidity, and vegetation, have been shown to affect mosquito breeding, survival, and behaviour, which in turn can shape genetic diversity (40). Furthermore, the human population density and urbanization in Jeddah may play a role in the genetic similarity of its mosquito population to those in Argentina and Australia, as urbanization tends to homogenize genetic diversity, due to the increased movement of people and goods. In contrast, the more rural setting of Jazan, with its agricultural activities, may provide a more heterogeneous environment, contributing to the greater genetic diversity observed in the mosquito population in this region.
4 Conclusion
This study successfully utilized molecular identification and phylogenetic analysis to investigate the genetic variation of the Aedes aegypti populations collected from Jeddah and Jazan, Saudi Arabia. The results revealed substantial genetic divergence between the two populations, emphasizing the complex population structure of this important disease vector within the different geographic regions of the Kingdom. Notably, the Jeddah population demonstrated a closer genetic relationship with populations from Argentina and Australia, suggesting a possible introduction pathway or a shared ancestral lineage. In contrast, the Jazan population exhibited a higher degree of genetic diversity with genetic similarities to mosquitoes from Canada, Kenya, India, and England, indicating a potentially more diverse origin or greater gene flow in this region. These findings highlight the need for region-specific mosquito control strategies. Genetic variability, both within and between populations, can significantly influence vector competence, insecticide resistance, and the overall effectiveness of control measures. Therefore, understanding the genetic structure of local mosquito populations is essential for designing targeted and sustainable vector control programs. This study highlights the importance of continuous genetic surveillance by using molecular tools. Monitoring the genetic dynamics of Ae. aegypti populations over time will be vital for tracking evolutionary trends, identifying potential resistance developments, and adapting public health interventions accordingly. In the context of Saudi Arabia, integrating molecular data into vector control policies will enhance the precision and success of mosquito management efforts, especially in high-risk urban areas, such as Jeddah and Jazan.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
SIA: Funding acquisition, Writing – review & editing. SMA: Investigation, Methodology, Conceptualization, Writing – original draft, Formal analysis, Validation. KA: Supervision, Writing – review & editing. SS: Conceptualization, Writing – original draft, Supervision, Validation. HA-S: Formal analysis, Writing – review & editing. AA: Validation, Writing – review & editing. HA: Writing – review & editing, Methodology. HA-R: Writing – review & editing, Investigation. JM: Validation, Supervision, Writing – original draft, Conceptualization.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number MoE-IF-UJ-R2-22-1037-1.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor FA-M declared a past co-authorship with the author JM.
Generative AI statement
The author(s) declare that Generative AI was used in the creation of this manuscript. ChatGPT-5 and Wordtune were used for language editing and structuring.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Mashlawi AM, Al-Nazawi AM, Noureldin EM, Alqahtani H, Mahyoub JA, Saingamsook J, et al. Molecular analysis of knockdown resistance (kdr) mutations in the voltage-gated sodium channel gene of Aedes aEgypti populations from Saudi Arabia. Parasites Vectors. (2022) 15:375. doi: 10.1186/s13071-022-05525-y
2. Somia ES, Ullah I, Alyahya HS, and Mahyoub JA. Identification of Wolbachia new strains from Aedes aEgypti mosquitoes, the vector of dengue fever in Jeddah Province. BMC Microbiol. (2023) 23:287. doi: 10.1186/s12866-023-03010-9
3. Mahyoub JA, Algamdi AG, and Aljameeli MM. Resistance development to bioinsecticides in Aedes aEgypti (Culicidae: Diptera), the vector of dengue fever in Saudi Arabia. Saudi J Biol Sci. (2023) 30:103776. doi: 10.1016/j.sjbs.2023.103776
4. Shepard DS, Undurraga EA, Halasa YA, and Stanaway JD. The global economic burden of dengue: a systematic analysis. Lancet Infect Dis. (2016) 16:935–41. doi: 10.1016/S1473-3099(16)00146-8
5. Fakeeh M and Zaki AM. Virologic and serologic surveillance for dengue fever in Jeddah, Saudi Arabia, 1994–1999. Am J Trop Med Hyg. (2001) 65:764–7. doi: 10.4269/ajtmh.2001.65.764
6. Alyahya HS. Prevalence of dengue fever in Saudi Arabia: Jeddah as a case study. Entomol Res. (2023) 53:539–53. doi: 10.1111/1748-5967.12685
7. Valjarević A, Algarni S, Morar C, Grama V, Stupariu M, Tiba A, et al. The coastal fog and ecological balance for plants in the Jizan region, Saudi Arabia. Saudi J Biol Sci. (2023) 30:103494. doi: 10.1016/j.sjbs.2022.103494
8. Al-Ghamdi K, Alikhan M, Mahyoub J, and Afifi ZI. Studies on identification and population dynamics of Anopheline mosquitoes from Jeddah Province of Saudi Arabia. Biosci Biotechnol Res Commun. (2008) 1:19–24.
9. Alhaeli A, Bahkali S, Ali A, Househ MS, and El-Metwally AA. The epidemiology of Dengue fever in Saudi Arabia: A systematic review. J Infect Public Health. (2016) 9:117–24. doi: 10.1016/j.jiph.2015.05.006
10. Al-Ghamdi K, Alikhan M, Mahyoub J, and Afifi ZI. Studies on identification and population dynamics of anopheline mosquitoes from Jeddah Province of Saudi Arabia. Biosci Biotechnol Res Commun. (2008) 1:19–24.
11. Al-Ghamdi K, Khan MA, and Mahyoub JA. Role of climatic factors in the seasonal abundance of Ae. aEgypti and dengue fever cases in Jeddah province of Saudi Arabia. Curr World Environ. (2010) 4.
12. Alikhan M, Al Ghamdi K, and Mahyoub JA. Aedes mosquito species in western Saudi Arabia. J Insect Sci. (2014) 14:69. doi: 10.1093/jis/14.1.69
13. Dieng H, Ahmad AH, Mahyoub JA, Turkistani AM, Mesed H, Koshike S, et al. Household survey of container–breeding mosquitoes and climatic factors influencing the prevalence of Aedes aEgypti (Diptera: Culicidae) in Makkah City, Saudi Arabia. Asian Pac J Trop Biomed. (2012) 2:849–57. doi: 10.1016/S2221-1691(12)60242-1
14. Mahyoub JA. Variation in the temporal and place effects on the dengue fever mosquito vector Aedes aEgypti (L.) in Jeddah Province. Biosci Biotechnol Res Asia. (2015) 12:45–50. Available online at: https://www.biotech-asia.org/?p=4429 (Accessed August 26, 2025).
15. Walker ED, Thibault AR, Thelen AP, Bullard BA, Huang J, Odiere MR, et al. Identification of field caught Anopheles Gambiae s.s. and Anopheles arabiensis by TaqMan single nucleotide polymorphism genotyping. Malar J. (2007) 6:23. doi: 10.1186/1475-2875-6-23
16. Laurito M and Almirón WR. Morphological variation in diagnostic features for two Culex (Culex) species of the Neotropical Region (Diptera: Culicidae). Zootaxa. (2015) 4052:573–6. doi: 10.11646/zootaxa.4052.5.5
17. Besansky NJ, Severson DW, and Ferdig MT. DNA barcoding of parasites and invertebrate disease vectors: what you don’t know can hurt you. Trends Parasitol. (2003) 19:545–6. doi: 10.1016/j.pt.2003.10.003
18. Pennisi E. Modernizing the tree of life. Science. (2003) 300:1692–7. doi: 10.1126/science.300.5626.1692
19. Mousson L, Dauga C, Garrigues T, Schaffner F, Vazeille M, and Failloux A-B. Phylogeography of Aedes (Stegomyia) aEgypti (L.) and Aedes (Stegomyia) albopictus (Skuse) (Diptera: Culicidae) based on mitochondrial DNA variations. Genet Res. (2005) 86:1–11. doi: 10.1017/S0016672305007627
20. Al Thabiani A. Distribution and mitochondrial CO1-based genetic diversity of Aedes aEgypti (L.) (Culicidae: Diptera) in Saudi Arabia. Saudi J Biol Sci. (2023) 30:103566. doi: 10.1016/j.sjbs.2023.103566
21. Shepard JJ, Andreadis TG, and Vossbrinck CR. Molecular phylogeny and evolutionary relationships among mosquitoes (Diptera: Culicidae) from the northeastern United States based on small subunit ribosomal DNA (18S rDNA) sequences. J Med Entomol. (2006) 43:443–54. doi: 10.1093/jmedent/43.3.443
22. Kumar NP, Rajavel AR, Natarajan R, and Jambulingam P. DNA barcodes can distinguish species of Indian mosquitoes (Diptera: Culicidae). J Med Entomol. (2007) 44:1–7. doi: 10.1093/jmedent/41.5.01
23. Pfeiler E, Flores-López CA, Mada-Vélez JG, Escalante-Verdugo J, and Markow TA. Genetic diversity and population genetics of mosquitoes (Diptera: Culicidae: Culex spp.) from the Sonoran Desert of North America. Sci World J. (2013) 2013:724609. doi: 10.1155/2013/724609
24. Al-Otaibi WM, Alghamdi KM, and Mahyoub JA. Molecular characterization and phylogenetic relationships among Rhynchophorus sp. haplotypes in Makkah Al-Mukarramah Region-KSA. Saudi J Biol Sci. (2022) 29:103388. doi: 10.1016/j.sjbs.2022.103388
25. Sabir JSM, Rabah S, Yacoub H, Hajrah NH, Atef A, Al-Matary M, et al. Molecular evolution of cytochrome C oxidase-I protein of insects living in Saudi Arabia. PloS One. (2019) 14:e0224336. doi: 10.1371/journal.pone.0224336
26. Al-Jameeli MM, Al-Ghamdi KM, Mahyoub JA, Sharawi SE, and Al-Ghamdi A. Morphological and molecular identification of screwworm fly in Jeddah city: Genus Chrysomya (Calliphoridae). Int J Entomol Res. (2021) 6:114–23.
27. Sharawi SE, Al Dagal MM, Alghamdi TS, Alkenani NA, Shafi ME, Alyahya HS, et al. Molecular identification of forensic insects collected from Jeddah, Saudi Arabia. Indian J Anim. Res. (2024), 1–9. doi: 10.18805/IJAR.BF-1826
28. Al-Otaibi WM, Algethami RHM, Hussien NA, and Mahyoub JA. Phylogenetic analysis of Aedes aEgypti in Taif and Jeddah governorates, Saudi Arabia. J Am Mosq Control Assoc. (2025). doi: 10.2987/24-7194. Advance online publication.
29. Soghigian J, Sither C, Justi SA, Morinaga G, Cassel BK, Vitek CJ, et al. Phylogenomics reveals the history of host use in mosquitoes. Nat Commun. (2023) 14:6252. doi: 10.1038/s41467-023-41764-y
30. Kress WJ and Erickson DL. DNA barcodes: methods and protocols. In: Kress WJ and Erickson DL, editors. DNA Barcodes: Methods and Protocols. Humana Press, Totowa, NJ (2012). p. 3–8.
31. Vrijenhoek RC. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. (1994) 3:294–9.
32. Lapadula WJ, Marcet PL, Taracena ML, Lenhart A, and Ayub MJ. Characterization of horizontally acquired ribotoxin encoding genes and their transcripts in Aedes aEgypti. Gene. (2020) 754:144857. doi: 10.1016/j.gene.2020.144857
33. Batovska J, Blacket MJ, Brown K, and Lynch SE. Molecular identification of mosquitoes (Diptera: Culicidae) in southeastern Australia. Ecol Evol. (2016) 6:3001–11. doi: 10.1002/ece3.2095
34. Makanda M, Wahome L, Onyambu G, Mutunga J, and Bargul J. Direct Submission. Nairobi, Kenya: Institute of Basic Sciences Technology and Innovation, The Pan African University (2018).
35. Al Nazawi AM, Aqili J, Alzahrani M, McCall PJ, and Weetman D. Combined target site (kdr) mutations play a primary role in highly pyrethroid resistant phenotypes of Aedes aEgypti from Saudi Arabia. Parasites Vectors. (2017) 10:1–10. doi: 10.1186/s13071-017-2096-6
36. Endersby-Harshman NM, Ali A, Alhumrani B, Alkuriji MA, Al-Fageeh MB, Al-Malik A, et al. Voltage-sensitive sodium channel (Vssc) mutations associated with pyrethroid insecticide resistance in Aedes aEgypti (L.) from two districts of Jeddah, Kingdom of Saudi Arabia: baseline information for a Wolbachia release program. Parasites Vectors. (2021) 14:1–13. doi: 10.1186/s13071-021-04867-3
37. Kasai S, Itokawa K, Uemura N, Takaoka A, Furutani S, Maekawa Y, et al. Discovery of super–insecticide-resistant dengue mosquitoes in Asia: Threats of concomitant knockdown resistance mutations. Sci Adv. (2022) 8:eabq7345. doi: 10.1126/sciadv.abq7345
38. Mashlawi AM, Al-Nazawi AM, Noureldin EM, Alqahtani H, Mahyoub JA, Saingamsook J, et al. Molecular analysis of knockdown resistance (kdr) mutations in the voltage-gated sodium channel gene of Aedes aEgypti populations from Saudi Arabia. Parasites Vectors. (2022) 15:375. doi: 10.1186/s13071-022-05525-y
39. Nayak B, Sahu B, Pattanayak AK, and Barik TK. Mitochondrial COX1 sequence of Aedes aEgypti. Berhampur, India: PG Department of Zoology, Berhampur University (2024).
Keywords: Aedes aegypti, Jeddah, Jizan, molecular insights, vector borne disease, dengue fever, phylogenetic
Citation: Alqurashi SI, Alqahtani SM, Alghamdi KMS, Sharawi SE, Al-Solami HM, Alghamdi AG, Alyahya HS, Al-Rashidi HS and Mahyoub JA (2025) Molecular insights into Aedes aegypti (L.) populations and vector surveillance in the urban areas of Jeddah and Jizan, Saudi Arabia. Front. Insect Sci. 5:1638582. doi: 10.3389/finsc.2025.1638582
Received: 30 May 2025; Accepted: 12 August 2025;
Published: 11 September 2025.
Edited by:
Fahd A. Al-Mekhlafi, King Saud University, Saudi ArabiaReviewed by:
Aleksandar Valjarević, University of Belgrade, SerbiaIracy Soares, Universidade Federal da Integração Latino-Americana, Brazil
Copyright © 2025 Alqurashi, Alqahtani, Alghamdi, Sharawi, Al-Solami, Alghamdi, Alyahya, Al-Rashidi and Mahyoub. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shatha I. Alqurashi, c2FhcXVyc2hpQHVqLmVkdS5zYQ==