Wing size and parity as markers of cohort demography for potential Anopheles (Culicidae: Diptera) malaria vectors in the Republic of Korea.

Wing lengths of parous (P) and nulliparous (NP) PCR-identified female Anopheles belenrae, An. kleini, An. pullus, and An. sinensis were determined from weekly trap collections at Camp Humphreys (CH), Ganghwa Island (GH), and Warrior Base (WB), Republic of Korea (ROK) during Jun-Oct, 2009. Wing length was greatest at the beginning and end of the study period. Wing length of NPs tended to be less than that of Ps before the period of maximum greening (Jul-Aug) but greater thereafter. Larger specimens tended to be Ps, and weekly wing length of Ps appeared less variable than NPs, possibly due to selection. A bimodal wing length frequency distribution of An. sinensis suggested two forms comprising small- (≤4.5 mm, SW) and large-winged females (>4.5 mm, LW). LW comprised the majority of peaks in abundance, however %SW, while still a minority, often increased during these times suggesting a density-dependent effect. At WB and GH, a two to three-week periodicity in %SW was obvious for An. sinensis and An. kleini. Analyses of weather station and satellite data showed that smaller-winged An. sinensis were associated with warmer, more humid, and greener times of the year. SW and LW specimens possibly result from agricultural practices that are common across large areas; regular synchronous peaks of SW and LW were observed from different sites. Peaks in SW Ps followed peaks in NPs in a 'ripple effect' one to two weeks apart, suggesting that wing length combined with parity could be used to follow the emergence and survival of mosquito cohorts.

Assessing the current and future potential geographic distribution of the American dog tick, Dermacentor variabilis (Say) (Acari: Ixodidae) in North America.

The American dog tick, Dermacentor variabilis, is a veterinary- and medically- significant tick species that is known to transmit several diseases to animal and human hosts. The spatial distribution of this species in North America is not well understood, however; and knowledge of likely changes to its future geographic distribution owing to ongoing climate change is needed for proper public health planning and messaging. Two recent studies have evaluated these topics for D. variabilis; however, less-rigorous modeling approaches in those studies may have led to erroneous predictions. We evaluated the present and future distribution of this species using a correlative maximum entropy approach, using publicly available occurrence information. Future potential distributions were predicted under two representative concentration pathway (RCP) scenarios; RCP 4.5 for low-emissions and RCP 8.5 for high-emissions. Our results indicated a broader current distribution of this species in all directions relative to its currently known extent, and dramatic potential for westward and northward expansion of suitable areas under both climate change scenarios. Implications for disease ecology and public health are discussed.

Forest disturbance and vector transmitted diseases in the lowland tropical rainforest of central Panama.

OBJECTIVE: To explore possible changes in the community attributes of haematophagous insects as a function of forest disturbance. We compare the patterns of diversity and abundance, plus the behavioural responses of three epidemiologically distinct vector assemblages across sites depicting various levels of forest cover. METHODS: Over a 3-year period, we sampled mosquitoes, sandflies and biting-midges in forested habitats of central Panama. We placed CDC light traps in the forest canopy and in the understorey to gather blood-seeking females. RESULTS: We collected 168 405 adult haematophagous dipterans in total, including 26 genera and 86 species. Pristine forest settings were always more taxonomically diverse than the disturbed forest sites, confirming that disturbance has a negative impact on species richness. Species of Phlebotominae and Culicoides were mainly classified as climax (i.e. forest specialist) or disturbance-generalist, which tend to decrease in abundance along with rising levels of disturbance. In contrast, a significant portion of mosquito species, including primary and secondary disease vectors, was classified as colonists (i.e. disturbed-areas specialists), which tend to increase in numbers towards more disturbed forest habitats. At pristine forest, the most prevalent species of Phlebotominae and Culicoides partitioned the vertical niche by being active at the forest canopy or in the understorey; yet this pattern was less clear in disturbed habitats. Most mosquito species were not vertically stratified in their habitat preference. CONCLUSION: We posit that entomological risk and related pathogen exposure to humans is higher in pristine forest scenarios for Culicoides and Phlebotominae transmitted diseases, whereas forest disturbance poses a higher entomological risk for mosquito-borne infections. This suggests that the Dilution Effect Hypothesis (DEH) does not apply in tropical rainforests where highly abundant, yet unrecognised insect vectors and neglected zoonotic diseases occur. Comprehensive, community level entomological surveillance is, therefore, the key for predicting potential disease spill over in scenarios of pristine forest intermixed with anthropogenic habitats. We suggest that changes in forest quality should also be considered when assessing arthropod-borne disease transmission risk.

OBJECTIF: Explorer les changements possibles dans les attributs communautaires des insectes hématophages en fonction de la perturbation des forêts. Nous comparons les modèles de diversité et d'abondance, ainsi que les réponses comportementales de trois assemblages de vecteurs épidémiologiquement distincts sur des sites illustrant divers niveaux de couverture forestière. MÉTHODES: Au cours d'une période de trois ans, nous avons échantillonné des moustiques, des phlébotomes et des moucherons piqueurs dans les habitats forestiers du centre de Panama. Nous avons placé des pièges à lumière CDC dans la canopée de la forêt et dans le sous-étage pour recueillir les femelles en quête de sang. RÉSULTATS: Nous avons collecté un total de 168.405 diptères hématophages adultes, dont 26 genres et 86 espèces. Les environnements de forêt intacts étaient toujours plus diversifiés du point de vue taxonomique que les sites forestiers perturbés, confirmant que les perturbations avaient un impact négatif sur la richesse en espèces. Les espèces de phlébotome et Culicoïdes étaient principalement classés comme climax (spécialiste de la forêt) ou généralistes de perturbation, qui ont tendance à diminuer en abondance parallèlement aux niveaux croissants de perturbation. En revanche, une partie importante des espèces de moustiques, y compris les vecteurs primaires et secondaires de maladies, a été classée dans la catégorie des colons (c'est-à-dire spécialistes des zones perturbées), qui ont tendance à se multiplier vers des habitats forestiers plus perturbés. Dans la forêt vierge, les espèces de phlébotomes et Culicoïdes les plus répandues cloisonnaient la niche verticale en étant actives dans la canopée de la forêt ou dans le sous-étage; pourtant, cette tendance était moins nette dans les habitats perturbés. La plupart des espèces de moustiques n'étaient pas stratifiées verticalement dans leur préférence d'habitat. CONCLUSION: Nous estimons que le risque entomologique et l'exposition associée des agents pathogènes à l'homme est plus élevé dans les scénarios de forêt vierge pour les maladies transmises par les phlébotomes et Culicoïdes, alors que la perturbation des forêts pose un risque entomologique plus élevé pour les infections transmises par les moustiques. Cela suggère que l'hypothèse de l'effet de dilution ne s'applique pas dans les forêts tropicales humides où se reproduisent très abondamment les insectes vecteurs, mais non reconnus, et où des maladies zoonotiques négligées surviennent. Une surveillance entomologique approfondie au niveau de la communauté est donc la clé pour prédire le potentiel de propagation des maladies dans des scénarios de forêt vierge mélangée à des habitats anthropiques. Nous suggérons que les changements dans la qualité des forêts soient également pris en compte lors de l'évaluation du risque de propagation de maladies transmises par les arthropodes.

An Excel Spreadsheet Tool for Exploring the Seasonality of Aedes Vector Hazard for User-Specified Administrative Regions of Brazil.

Aedes-vectored viruses are a major concern for active-duty military personnel working in South and Central America at certain times of the year. Knowledge about the seasonal changes of vector activity is important as it informs time-sensitive vector control, prophylaxis, and travel decisions. To assist in-country and extralimital efforts to anticipate when vector hazards and the risks of transmission are highest, we developed an Excel spreadsheet tool that uses published monthly habitat suitability models to display various aspects of average Aedes seasonality for user-defined second order administrative areas of Brazil. This tool expands on those previously developed by the authors for the contiguous United States, with the aim of translating global habitat suitability models into user-friendly formats to provide actionable intelligence for areas of interest.

Mosquito Surveillance Conducted by US Military Personnel in the Aftermath of the Nuclear Explosion at Nagasaki, Japan, 1945.

Mosquito surveillance data can be used to develop bionomic profiles of vector species to inform abatement plans. Thus, surveillance was conducted in the months following Allied occupation of Japan at the conclusion of World War II. Mosquito surveillance in Nagasaki, Japan, began one month after the nuclear bomb destroyed much of the city. The resulting specimens housed within the US National mosquito collection are documented here for the first time. Specimen labels were digitized and specimens were photographed to record specimen condition as part of the process for making them readily available to researchers.

Seasonal dynamics of Anopheles species at three locations in the Republic of Korea.

Weekly changes in adult Anopheles species were monitored at Camp Humphreys (CH), Ganghwa Island (GH), and Warrior Base (WB), from May-October, 2009-2010 to explore the relationship between Plasmodium vivax development and vector dynamics in the Republic of Korea (ROK). Adult females were trapped and dissected to determine parity for estimating longevity, mortality, and birthrate. A degree-day (DD) method was used to estimate the extrinsic incubation period (EIP) of P. vivax and duration of the gonotrophic cycle and other life stages. Anopheles sinensis was the predominant species, with satellite data showing peak abundance occurring after the period of maximum greenness. Abundance peaks were location dependent, comprised nulliparous and parous females, and timing could not be fully explained by DD estimation. Parity showed synchronicity between locations and years and was highest for September and lowest during maximum greenness. Mosquito longevity was predicted to exceed the EIP (when malaria transmission is possible) during weeks 29, 31, 34, and near the end of the season. Area-wide changes in parity suggest a common cause; information on local larval habitat and agricultural practices may explain location-specific effects. DD estimates of EIP and parity could be used to predict when conditions are suitable for P. vivax transmission.

Disturbance and mosquito diversity in the lowland tropical rainforest of central Panama.

The Intermediate Disturbance Hypothesis (IDH) is well-known in ecology providing an explanation for the role of disturbance in the coexistence of climax and colonist species. Here, we used the IDH as a framework to describe the role of forest disturbance in shaping the mosquito community structure, and to identify the ecological processes that increase the emergence of vector-borne disease. Mosquitoes were collected in central Panama at immature stages along linear transects in colonising, mixed and climax forest habitats, representing different levels of disturbance. Species were identified taxonomically and classified into functional categories (i.e., colonist, climax, disturbance-generalist, and rare). Using the Huisman-Olff-Fresco multi-model selection approach, IDH testing was done. We did not detect a unimodal relationship between species diversity and forest disturbance expected under the IDH; instead diversity peaked in old-growth forests. Habitat complexity and constraints are two mechanisms proposed to explain this alternative postulate. Moreover, colonist mosquito species were more likely to be involved in or capable of pathogen transmission than climax species. Vector species occurrence decreased notably in undisturbed forest settings. Old-growth forest conservation in tropical rainforests is therefore a highly-recommended solution for preventing new outbreaks of arboviral and parasitic diseases in anthropic environments.

A location-specific spreadsheet for estimating Zika risk and timing for Zika vector surveillance, using US military facilities as an example.

Local Zika virus transmission in the United States involving one or both of the known vector species, Aedes aegypti and Ae. albopictus, is of major concern. To assist efforts to anticipate the risks of transmission, we developed an Excel spreadsheet tool that uses vector and virus temperature thresholds, remotely sensed maximum temperature, and habitat suitability models to answer the questions: "is Zika transmission likely here?" and "when should we conduct vector surveillance?" An example spreadsheet, updated regularly and freely available, uses near real-time and forecast temperature data to generate guidance, based on a novel four level Zika risk code, for 733 US military facilities in the 50 states, the District of Columbia, and the territories of Guam and Puerto Rico.

A Spatiotemporal Database to Track Human Scrub Typhus Using the VectorMap Application.

Scrub typhus is a potentially fatal mite-borne febrile illness, primarily of the Asia-Pacific Rim. With an endemic area greater than 13 million km2 and millions of people at risk, scrub typhus remains an underreported, often misdiagnosed febrile illness. A comprehensive, updatable map of the true distribution of cases has been lacking, and therefore the true risk of disease within the very large endemic area remains unknown. The purpose of this study was to establish a database and map to track human scrub typhus. An online search using PubMed and the United States Armed Forces Pest Management Board Literature Retrieval System was performed to identify articles describing human scrub typhus cases both within and outside the traditionally accepted endemic regions. Using World Health Organization guidelines, stringent criteria were used to establish diagnoses for inclusion in the database. The preliminary screening of 181 scrub typhus publications yielded 145 publications that met the case criterion, 267 case records, and 13 serosurvey records that could be georeferenced, describing 13,739 probable or confirmed human cases in 28 countries. A map service has been established within VectorMap (www.vectormap.org) to explore the role that relative location of vectors, hosts, and the pathogen play in the transmission of mite-borne scrub typhus. The online display of scrub typhus cases in VectorMap illustrates their presence and provides an up-to-date geographic distribution of proven scrub typhus cases.

Records and Distribution of New World Phlebotomine Sand Flies (Psychodidae, Diptera), With Special Emphasis on Primary Types and Species Diversity.

This article includes the records and distribution of Phlebotomine sand flies (Psychodidae, Diptera) in the New World based on the specimen collections housed in 2 repositories, the US National Museum of Natural History and the Museum of Entomology, Florida State Collection of Arthropods. Approximately 128 species have primary types housed in the 2 repositories, including holotypes (47 species, 3 subspecies), "types" (7 species), allotypes (52 species, 6 subspecies), lectotypes (4 species), paratypes (93 species, 10 subspecies), and neoallotype (1 species), mounted on slides, with a total of 1,107 type slides. For species diversity, collection data from 24 countries in the sand fly database were analyzed according to the number of species present, specimen records, decade of collections, and countries where collections were conducted.

Quantitative factors proposed to influence the prevalence of canine tick-borne disease agents in the United States.

The Companion Animal Parasite Council hosted a meeting to identify quantifiable factors that can influence the prevalence of tick-borne disease agents among dogs in North America. This report summarizes the approach used and the factors identified for further analysis with mathematical models of canine exposure to tick-borne pathogens.

Geographic distribution, evolution, and disease importance of species within the Neotropical Anopheles albitarsis Group (Diptera, Culicidae).

The Anopheles albitarsis group of mosquitoes comprises eight recognized species and one mitochondrial lineage. Our knowledge of malaria vectorial importance and the distribution and evolution of these taxa is incomplete. We constructed ecological niche models (ENMs) for these taxa and used hypothesized phylogenetic relationships and ENMs to investigate environmental and ecological divergence associated with speciation events. Two major clades were identified, one north (Clade 1) and one south (Clade 2) of the Amazon River that likely is or was a barrier to mosquito movement. Clade 1 species occur more often in higher average temperature locations than Clade 2 species, and taxon splits within Clade 1 corresponded with a greater divergence of variables related to precipitation than was the case within Clade 2. Comparison of the ecological profiles of sympatric species and sister species support the idea that phylogenetic proximity is related to ecological similarity. Anopheles albitarsis I, An. janconnae, and An. marajoara ENMs had the highest percentage of their predicted suitable habitat overlapping distribution models of Plasmodium falciparum and P. vivax, and warrant additional studies of the transmission potential of these species. Phylogenetic proximity may be related to malaria vectorial importance within the Albitarsis Group.

Plasmodium-specific molecular assays produce uninterpretable results and non-Plasmodium spp. sequences in field-collected Anopheles vectors.

The Malaria Research and Reference Reagent Resource-recommended PLF/UNR/VIR polymerase chain reaction (PCR) was used to detect Plasmodium vivax in Anopheles spp. mosquitoes collected in South Korea. Samples that were amplified were sequenced and compared with known Plasmodium spp. by using the PlasmoDB.org Basic Local Alignment Search Tool/n and the National Center for Biotechnology Information Basic Local Alignment Search Tool/n tools. Results show that the primers PLF/UNR/VIR used in this PCR can produce uninterpretable results and non-specific sequences in field-collected mosquitoes. Three additional PCRs (PLU/VIV, specific for 18S small subunit ribosomal DNA; Pvr47, specific for a nuclear repeat; and GDCW/PLAS, specific for the mitochondrial marker, cytB) were then used to find a more accurate and interpretable assay. Samples that were amplified were again sequenced. The PLU/VIV and Pvr47 assays showed cross-reactivity with non-Plasmodium spp. and an arthropod fungus (Zoophthora lanceolata). The GDCW/PLAS assay amplified only Plasmodium spp. but also amplified the non-human specific parasite P. berghei from an Anopheles belenrae mosquito. Detection of P. berghei in South Korea is a new finding.

Amblyomma tapirellum (Acari: Ixodidae) collected from tropical forest canopy.

Free-ranging ticks are widely known to be restricted to the ground level of vegetation. Here, we document the capture of the tick species Amblyomma tapirellum in light traps placed in the forest canopy of Barro Colorado Island, central Panama. A total of forty eight adults and three nymphs were removed from carbon dioxide-octenol baited CDC light traps suspended 20 meters above the ground during surveys for forest canopy mosquitoes. To our knowledge, this represents the first report of questing ticks from the canopy of tropical forests. Our finding suggests a novel ecological relationship between A. tapirellum and arboreal mammals, perhaps monkeys that come to the ground to drink or to feed on fallen fruits.

Mosquitoes of eastern Amazonian Ecuador: biodiversity, bionomics and barcodes.

Two snapshot surveys to establish the diversity and ecological preferences of mosquitoes (Diptera: Culicidae) in the terra firme primary rain forest surrounding the Tiputini Biodiversity Station in the UNESCO Yasuní Biosphere Reserve of eastern Amazonian Ecuador were carried out in November 1998 and May 1999. The mosquito fauna of this region is poorly known; the focus of this study was to obtain high quality link-reared specimens that could be used to unequivocally confirm species level diversity through integrated systematic study of all life stages and DNA sequences. A total of 2,284 specimens were preserved; 1,671 specimens were link-reared with associated immature exuviae, all but 108 of which are slide mounted. This study identified 68 unique taxa belonging to 17 genera and 27 subgenera. Of these, 12 are new to science and 37 comprise new country records. DNA barcodes [658-bp of the mtDNA cytochrome c oxidase (COI) I gene] are presented for 58 individuals representing 20 species and nine genera. DNA barcoding proved useful in uncovering and confirming new species and we advocate an integrated systematics approach to biodiversity studies in future. Associated bionomics of all species collected are discussed. An updated systematic checklist of the mosquitoes of Ecuador (n=179) is presented for the first time in 60 years.

Mosquitoes of eastern Amazonian Ecuador: biodiversity, bionomics and barcodes

Two snapshot surveys to establish the diversity and ecological preferences of mosquitoes (Diptera: Culicidae) in the terra firme primary rain forest surrounding the Tiputini Biodiversity Station in the UNESCO Yasuní Biosphere Reserve of eastern Amazonian Ecuador were carried out in November 1998 and May 1999. The mosquito fauna of this region is poorly known; the focus of this study was to obtain high quality link-reared specimens that could be used to unequivocally confirm species level diversity through integrated systematic study of all life stages and DNA sequences. A total of 2,284 specimens were preserved; 1,671 specimens were link-reared with associated immature exuviae, all but 108 of which are slide mounted. This study identified 68 unique taxa belonging to 17 genera and 27 subgenera. Of these, 12 are new to science and 37 comprise new country records. DNA barcodes [658-bp of the mtDNA cytochrome c oxidase ( COI ) I gene] are presented for 58 individuals representing 20 species and nine genera. DNA barcoding proved useful in uncovering and confirming new species and we advocate an integrated systematics approach to biodiversity studies in future. Associated bionomics of all species collected are discussed. An updated systematic checklist of the mosquitoes of Ecuador (n = 179) is presented for the first time in 60 years.

SandflyMap: leveraging spatial data on sand fly vector distribution for disease risk assessments.

We feature SandflyMap (www.sandflymap.org), a new map service within VectorMap (www.vectormap.org) that allows free public online access to global sand fly, tick and mosquito collection records and habitat suitability models. Given the short home range of sand flies, combining remote sensing and collection point data give a powerful insight into the environmental determinants of sand fly distribution. SandflyMap is aimed at medical entomologists, vector disease control workers, public health officials and health planners. Data are checked for geographical and taxonomic errors, and are comprised of vouchered specimen information, and both published and unpublished observation data. SandflyMap uses Microsoft Silverlight and ESRI's ArcGIS Server 10 software platform to present disease vector data and relevant remote sensing layers in an online geographical information system format. Users can view the locations of past vector collections and the results of models that predict the geographic extent of individual species. Collection records are searchable and downloadable, and Excel collection forms with drop down lists, and Excel charts to country, are available for data contributors to map and quality control their data. SandflyMap makes accessible, and adds value to, the results of past sand fly collecting efforts. We detail the workflow for entering occurrence data from the literature to SandflyMap, using an example for sand flies from South America. We discuss the utility of SandflyMap as a focal point to increase collaboration and to explore the nexus between geography and vector-borne disease transmission.

Synchronous peaks in trap catches of malaria-infected mosquito species at Daeseongdong, a border village between North and South Korea.

Malaria continues to be a major health threat near the Demilitarized Zone (DMZ) that separates North and South Korea. Adult mosquitoes were collected from 20 July through 21 October, 2010 at Daeseongdong, a small village within the DMZ. Molecular techniques were used to identify Anopheles to species and for detection of Plasmodium vivax sporozoites in their head and thorax. Trap catches showed concordant peaks of Anopheles belenrae and An. kleini early in the study period and concordant peaks of An. pullus and An. sinensis later in the season. Three well defined peaks of the 107 sporozoite positive mosquitoes were observed: 34.6% were An. kleini, 23.4% were An. belenrae, 21.5% were An. sinensis, 19.6% were An. pullus, and 0.9% were An. lesteri. Estimation of the extrinsic incubation period from daily temperatures did not help identify preceding biting peaks of An. pullus and An. sinensis, when infection should have been acquired. We explore possible reasons for the sudden appearance and disappearance of sporozoite-infected mosquitoes, including the influx of infected mosquitoes from adjoining areas, and weather patterns. Regular surveillance for infected mosquitoes near border areas of the Republic of Korea may provide advance warning of increased malaria risk potential.

Mosquito bisection as a variable in estimates of PCR-derived malaria sporozoite rates.

BACKGROUND: Highly sensitive polymerase chain reaction (PCR) methods offer an alternative to the light microscopy examination of mosquito salivary glands for the determination of malaria sporozoite rates in wild caught female Anopheles. Removal of mosquito abdomens is assumed to eliminate false positives caused by malaria oocyst DNA in the midgut. This assumption has not been tested with current gold standard PCR assays, and for the variety of conditions that specimens could encounter in the laboratory and field. METHODS: Laboratory Anopheles stephensi were used that had been infected with Plasmodium falciparum 6-7 days and 14 days post infection (p.i.), when oocysts only and oocysts + sporozoites, respectively, are developed. Mosquitoes were killed and immediately frozen, air dried before being frozen, or stored under humid conditions overnight before being frozen, to simulate a range of conditions in the field. Additionally, abdomens were removed anterior to, at, or posterior to the junction of the abdomen and thorax, and both portions were processed using a standard nested PCR of the small sub-unit nuclear ribosomal genes (ssrDNA) with products visualized on agarose gels. RESULTS: Overall, 4.1 % (4/97) of head + thorax samples that were 6-7 days p.i. gave apparent false positives for sporozoites, compared to 9.3 % (9/97) that were positive for abdomens. No positives (0/52) were obtained when similar specimens were bisected anterior to the junction of the thorax and abdomen, compared to 21.2 % (11/52) that were positive for posterior portions. Multiple bands were noted for positives from the 'Frozen' treatment and the rate of false negatives due to DNA degradation appears higher under the 'Humid' treatment. Reproducibility of results for the 'Frozen' treatment was 90 %. CONCLUSIONS: Despite the importance of specimen condition and the bisection step in determining sporozoite rates, little attention has been paid to them in the literature. Recommendations from this study are that: 1) care needs to be taken to reduce DNA degradation in the field; 2) mosquito abdomens be separated anterior to the junction of the thorax and abdomen; and 3) DNA sequencing of a subsample of positive results should be undertaken if possible.

Mosquito species composition and Plasmodium vivax infection rates on Baengnyeong-do (island), Republic of Korea.

Vivax malaria is a significant military and civilian health threat in the north of the Republic of Korea (ROK). The island of Baengnyeong-do is the westernmost point of the ROK and is located close to the southwestern coast of the Democratic People's Republic of Korea (DPRK). Mosquitoes were collected using a black light trap on Baengnyeong-do, and Anopheles spp. were assayed by PCR, to identify the species, and screened for sporozoites of Plasmodium vivax. Of a subsample of 257 mosquitoes, Anopheles lesteri was the most frequently collected (49.8%), followed by Anopheles sinensis (22.6%), Anopheles pullus (18.7%), Anopheles kleini (7.8%), and Anopheles belenrae (1.2%). The overall sporozoite rate was 3.1%, with the highest rates observed in An. kleini (15.0%), An. sinensis (5.2%), and An. lesteri (1.6%). No sporozoite positive An. pullus or An. belenrae were observed. The results extend our knowledge of the distribution and potential role in malaria transmission of An. kleini, An. lesteri, and An. sinensis, for an area previously considered to be at a low risk for contracting vivax malaria.