Figure 1

Figure 1. Locations in Yemen where invasive Anopheles stephensi mosquitoes were detected in 2021 and 2022 (yellow diamonds). Map was created by using MapChart (https://www.mapchart.net).

We analyzed immature mosquitoes collected from 2 semiurban locations, Ad Dahi and Zabid districts (Figure 1; Appendix). We collected the Ad Dahi district specimens over the course of a single day in December 2021 while conducting Aedes surveillance during a dengue fever outbreak. The Zabid district specimens were collected during monthlong Anopheles surveillance in March 2022. Potential breeding containers surveyed included open concrete ponds and cement water tanks near block factories, washing basins in mosques, and car washing sites. Mosquitos used in this study were collected as immature specimens by using the dipping method. We reared them to adults in field insectaries and then identified them by using updated morphological keys (12). The specimens we morphologically identified as An. stephensi mosquitoes were preserved with silica gel. We sent a subsample of specimens to Baylor University (Waco, Texas, USA) for molecular analysis.

For species identification, we analyzed 2 loci, cytochrome c oxidase subunit 1 (COI) and internal transcribed spacer 2 (ITS2), by using previously described protocols (13). We used an An. stephensi–specific endpoint PCR as previously described (14). We conducted ITS2- and COI-targeted DNA sequencing, then used BLAST (https://blast.ncbi.nlm.nih.gov) and conducted phylogenetic analysis to determine species identification and evaluate the level of genetic diversity (Appendix).

Most of the Anopheles larvae were collected from cement water tanks found outside of homes, which were like other invasive settings found in eastern Ethiopia and Aden (4,10). We identified 41 mosquitoes morphologically as An. stephensi, 7 from Ad Dahi and 34 from Zabid City. Of the 41 specimens, all 7 Ad Dahi and 32 of the Zabid City specimens were confirmed to be An. stephensi mosquitoes by ITS2 endpoint assay and by COI and ITS2 targeted DNA sequence analysis (Appendix). We observed no discrepancy between the endpoint assay and sequencing assay for species identification. We did not report any discrepancies with the ITS2 endpoint assay results when compared with the DNA sequence results, but incorrect identification with the endpoint assay for other populations with unreported ITS2 variation is possible. Sequencing should accompany initial implementation of an endpoint assay in newly surveyed populations to evaluate both false negatives and false positives.

Figure 2

Figure 2. Maximum-likelihood phylogenetic analysis of cytochrome c oxidase subunit 1 gene sequences for Anopheles stephensimosquitoes collected in Yemen. Pink indicates sequences from Yemen, blue indicates sequences from...

Figure 3

Figure 3. Maximum-likelihood phylogenetic analysis of An. internal transcribed spacer 2 DNA sequences for Anopheles stephensimosquitoes collected in Yemen. Only 1 sequence (bold) is included as a representative...

Two non–An. stephensi mosquitoes were recovered; 1 was identified as Anopheles culicifacies and 1 as Aedes aegypti on the basis of ITS2 and COI sequence BLAST analysis. This change from the initial morphological identification could be the result of sorting the Anopheles and Aedes mosquitoes from the same habitat or from the selection of specimens for molecular analysis in the laboratory. Phylogenetic analysis confirmed the species identification of the An. stephensi mosquito specimens (bootstrap = 100 for COI, bootstrap = 100 for ITS2) (Figures 2, 3).

ITS2 sequences were all identical for the An. stephensi mosquito specimens and BLAST analysis revealed a 100% sequence identification in GenBank. Sequence analysis of COI revealed 2 haplotypes. One haplotype (n=35) was previously reported throughout the Horn of Africa (HoA) (COIHap 3) with the highest frequency in northeastern Ethiopia, Somaliland, and Djibouti (15). This haplotype was also noted in a recent report on the detection of An. stephensi mosquitoes in Aden, Yemen (10). The second haplotype (n=4), differing by a single nucleotide from COIHap3, has not previously been reported (designated COIHapYem1).

These findings demonstrate the successful use of ongoing vector surveillance activities for the initial detection of An. stephensi mosquitoes. The detection of a common HoA COI haplotype raises questions about the relationship between the invasive An. stephensi mosquitoes in northern HoA and Yemen. It is possible they share a common origin or that movement of An. stephensi mosquitoes has occurred between the 2 regions. Sequencing of additional loci could assist with delineating their relationship. Our findings provide information about the relative timing of the Yemen introduction of An. stephensi mosquitoes relative to the HoA detection on the basis of several observations. Fewer haplotypes were observed in the Al Hudaydah region in comparison to northeastern Ethiopia, Somalia, and Djibouti, which may indicate a recent introduction into the Al Hudaydah region. COIHap3 has a limited geographic range in comparison to the HoA-wide Hap2, which suggests the haplotype from our study could be associated with a later introduction in HoA and potentially Yemen. As for the HoA specimens and the study in Aden, no Saudi Arabian haplotypes were detected. Unlike the HoA specimens, no South Asia haplotypes were detected. Further genomic analysis and extensive An. stephensi mosquito sampling in Yemen, Saudi Arabia, and other parts of the Arabian Peninsula are needed to evaluate the hypothesis of a recent introduction in Yemen relative to the HoA.

Our findings provide insight into the genetic diversity of An. stephensi mosquito populations in Al Hudaydah governorate, Yemen. The findings also provide support for the need of long-term entomological and epidemiologic surveillance of An. stephensi mosquitoes on malaria transmission in the region. With the detection of An. stephensi mosquitoes in Yemen, concerns remain related to the status of these mosquitoes in other areas of Yemen and their relationship with the HoA An. stephensi mosquitoes. Additional genomic analysis should be conducted to further examine the relative timing of the introduction of An. stephensi mosquitoes to Yemen.