- Open Access
Prevalence of Plasmodium spp. in Anopheles mosquitoes in Thailand: a systematic review and meta-analysis
Parasites & Vectors volume 15, Article number: 285 (2022)
The entomological inoculation rate (EIR) is one of the key indices used to evaluate malaria transmission and vector control interventions. One of the components of the EIR is the sporozoite rate in Anopheles vectors. A systematic review and meta-analysis was performed to identify the prevalence of Plasmodium spp. in field-collected Anopheles species across Thailand.
This systematic review was registered under the PROSPERO number CRD42021297255. Studies that focused on the identification of Plasmodium spp. in Anopheles mosquitoes were identified from the electronic databases PubMed, Web of Science, and Scopus. The quality of the identified studies was determined using the Strengthening the Reporting of Observational Studies in Epidemiology approach. The proportion of Anopheles mosquitoes collected, Anopheles vectors for Plasmodium species, and specificity of Anopheles vectors for Plasmodium species were analyzed. The pooled prevalence of Plasmodium species among the primary vectors (Anopheles dirus, Anopheles minimus, and Anopheles maculatus) was estimated using the random-effects model.
Of the 1113 studies identified, 31 were included in the syntheses. Of the 100,910 Anopheles mosquitoes identified for species and sibling species, An. minimus (40.16%), An. maculatus (16.59%), and Anopheles epiroticus (9.18%) were the most prevalent Anopheles species. Of the 123,286 Anopheles mosquitoes identified, 566 (0.46%) were positive for Plasmodium species. The highest proportions of Plasmodium species were identified in Anopheles hodgkini (2/6, 33.3%), Anopheles nigerrimus (2/24, 8.33%), Anopheles balabacensis (4/84, 4.76%), An. dirus (114/4956, 2.3%), Anopheles annularis (16/852, 1.88%), Anopheles kochi (8/519, 1.54%), Anopheles vagus (3/215, 1.4%), and Anopheles baimaii (1/86, 1.16%). The pooled prevalence of Plasmodium species identified in the main Anopheles vectors was 0.4% of that of Plasmodium species identified in An. dirus was 2.1%, that of Plasmodium species identified in An. minimus was 0.4%, and that of Plasmodium species identified in An. maculatus was 0.4%.
We found a low prevalence of Plasmodium infection in Anopheles mosquitoes across Thailand. Therefore, the use of EIR to determine the impact of vector control intervention on malaria parasite transmission and elimination in Thailand must be undertaken with caution, as a large number of Anopheles specimens may be required.
Malaria is one of the major causes of morbidity and mortality in the Greater Mekong Subregion countries, including Thailand [1, 2]. Thailand is located at the center of the Indochinese Peninsula in Southeast Asia. It shares a border with Myanmar, Laos, Cambodia, and Malaysia, leading to a constant migration of foreign workers and refugees across the borders . The presence of efficient vectors and long rainy season in Thailand indicates that there are several areas that can serve as mosquito habitats and that mosquito-borne diseases present a life-threatening public health challenge from both indigenous and imported cases . Despite ongoing efforts to combat malaria, approximately 4473 confirmed cases were reported in Thailand in 2020. Plasmodium vivax is the dominant cause of infection (4098 cases, 92%), and the number of cases involving Plasmodium falciparum is much lower (256 cases, 5.7%) but still significant . In Thailand, malaria cases are reported to have a high prevalence in vulnerable forest and forest fringes along rural stretches of border areas [3, 4]. Tak Province, which is a neighbor to Myanmar, had one of the highest incidences of malaria in Thailand in 2020 (1241 cases), followed by Yala Province in peninsular Thailand, on the border with Malaysia (1075 cases), and Kanchanaburi Province, another neighbor of Myanmar, with 539 confirmed cases . Military personnel, forest workers, refugees, and local migrants are the individuals at the highest risk . Although the malaria cases in Thailand in 2020 decreased by 24% compared to those in 2019 (5859 cases) , considerable effort is needed to achieve malaria elimination by 2024 . In addition to techniques for the diagnosis of malaria and effective therapy, mosquito control is one of the most effective interventions against malaria. Distribution of insecticide-treated nets (ITNs) and indoor residual spraying (IRS) are currently used [5, 7]. With the implementation of the National Malaria Elimination Strategy 2017–2026, ITNs have started being distributed at a ratio of one ITN per two people, with the aim of at least 90% coverage in each transmission focus. Retreatment of the net is performed regularly every 6–12 months. However, if the retreatment process is inaccessible, long-lasting insecticidal nets (LLINs) are distributed. IRS is implemented in cases where both ITNs and LLINs cannot be allocated . In 2020, over 102,150 LLINs were distributed with an average of 75% coverage across Thailand . Unfortunately, these methods tackle only indoor- and late-biting vectors, and indoor-resting Anopheles mosquitoes are not suitable for reducing outdoor transmission. Thus, additional approaches are needed to protect people from outdoor- and early-biting vectors .
The malaria parasite infects humans through the bite of infected female Anopheles mosquitoes. Of the 540 Anopheles species described, 79 species are found in Thailand . These species are normally classified into species complexes to which they are closely related and phenotypically indistinguishable . However, due to their heterogeneity in ecology, bionomic pattern, and vectorial capacity across their distribution, different members of the species complex present different epidemiological roles in malaria transmission [9, 10]. Therefore, accurate species identification is important for evaluating transmission dynamics and applying the most appropriate vector control interventions . In Thailand, seven species have been implicated as primary malaria vectors, in the Minimus and Dirus complexes and the Maculatus Group [11,12,13,14]. These species include Anopheles dirus Peyton & Harrison and Anopheles baimaii Sallum & Peyton of the Dirus Complex [15, 16]; Anopheles minimus Theobald of the Minimus Complex [15, 17]; Anopheles aconitus Dönitz of the Aconitus Subgroup (Funestus Group) [16, 18]; and Anopheles maculatus Theobald, Anopheles pseudowillmori Theobald, and Anopheles sawadwongporni Rattanarithikul & Green of the Maculatus Group [15, 19, 20]. Additionally, Anopheles campestris Reid, Anopheles barbirostris van der Wulp (Barbirostris Group) [21, 22], and Anopheles epiroticus Linton & Harbach (Sundaicus Complex)  have also been identified as potential malaria vectors in Thailand. Using serological or molecular assays, the rates of natural Plasmodium infection in the Dirus Complex, Minimus Complex, Maculatus Group, Barbirostris Group, and Sundaicus Complex were 0.8–6.4%, 0.09–5%, 0.1–3.1%, 0.42–1.9%, and 0.97%, respectively . This information is crucial in determining the vector capacity of Anopheles species , as well as in planning vector control strategies.
To assess the transmission dynamics as well as the effectiveness of the vector control methods, it has been suggested that the entomological inoculation rate (EIR) should be evaluated annually . The EIR measures the frequency of infectious bites by an Anopheles mosquito per person over time, combining the human biting rate and the sporozoite rate (SR) . However, not all infectious bites result in blood-stage malaria in human hosts because Plasmodium parasites possibly experience significant bottlenecks during their sporogony cycle in mosquitoes . The SR remains one of the important entomological indicators not only in assessing EIR but also in identifying Anopheles vectors that contribute to Plasmodium transmission in endemic settings . Indicators other than EIR could also influence the transmission dynamics. As EIR can be estimated from the human biting rate, it has been demonstrated that very high vector density and high SR could sustain malaria transmission over a large part of the year . In view of the evaluation of pre-oocyst formation blocking interventions, (e.g., gametocytocidal drugs), a previous study suggested that the oocyst formation rate is another highly reliable entomological indicator of mosquito infectiveness . Although research into naturally infected Anopheles mosquito has been conducted for decades, there is a need for a comprehensive systematic review and meta-analysis focusing on the prevalence of malaria parasites in field-collected Anopheles species across Thailand. We herein mainly focused on the combined prevalence of both sporozoite and oocyst infection rates, hereby termed “Plasmodium infection.” The information collected and synthesized in this study improves our understanding of the local transmission dynamics of malaria vector species, particularly primary vectors, and may offer useful data for the evaluation of vector control interventions and malaria transmission.
The systematic review followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines . The systematic review was registered at PROSPERO with the number CRD42021297255.
Relevant studies were searched from three research databases, namely, PubMed, Web of Science, and Scopus, from inception to March 30, 2021. The reference lists of the included studies were also examined to ensure that relevant studies were not missed. We combined relevant search terms with Boolean operators, and the terms “(malaria OR Plasmodium) and (anopheles OR anopheline) and (Thailand OR Thai OR Siam)” were used to identify relevant studies in each database (Additional file 1: Table S1). Studies that focused on the identification of Plasmodium spp. in Anopheles mosquitoes were selected. The reference lists of the included studies were also screened to ensure that relevant studies were not missed. The search included only studies that were published in English between 1945 and 2021.
Eligibility criteria and study selection
Studies were selected using the PICO method. The elements of this method are as follows: P: Participants. The participants included in the study were Anopheles mosquitoes in Thailand. I: Intervention. No intervention was applied in the present study. C: Comparator. No comparator was used in the present study. O: Outcome. The outcome of interest was the presence of Plasmodium spp. in any stage that was identified in Anopheles mosquitoes. Thus, the inclusion criteria were composed of cross-sectional studies that identified Plasmodium spp. among Anopheles mosquitoes collected in Thailand. The exclusion criteria were studies with incomplete data for extraction, studies for which the full text was unavailable, reviews or systematic reviews, in vitro studies, papers describing the development of assays, and letters to the editor/comments/editorials. Two authors (CS and MK) independently screened the titles and abstracts and selected studies based on the eligibility criteria. First, the titles and abstracts generated by the electronic search were checked. Second, the full texts were examined, and studies that did not meet the eligibility criteria were excluded, with the reasons recorded. Any differences in study selection between the two authors were resolved by mutual consensus.
Data extraction and quality assessment
Pilot data extraction tables were used to collect information from each included study. The following information was collected: name of the first author, year of publication, study sites, season and time at which mosquitoes were collected, mosquito collection methods, mosquito species, number of mosquitoes, methods of Plasmodium detection, Plasmodium identified, and stage of Plasmodium spp. The quality of the studies included was determined using the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement, which includes 22 parameters . The quality of each study was assessed as high quality (> 75%), moderate quality (50–75%), or low quality (< 50%). High-quality and moderate-quality studies were included in the systematic review. Low-quality studies were excluded.
The proportion of Anopheles mosquitoes collected in the studies, Anopheles vectors positive for Plasmodium species, and Plasmodium spp. identified in Anopheles mosquitoes are presented as frequencies and percentages. The pooled prevalence of Plasmodium species among Anopheles mosquitoes and primary vectors (An. dirus, An. minimus, and An. maculatus) was estimated using a random-effects meta-analysis using the DerSimonian and Laird method . The proportions from each study were pooled using logit transformation, and back transformation to a proportion was performed using generalized linear mixed models (GLMMs). The individual study weights were not calculated by the GLMMs. The meta-analyses of proportion studies were conducted using the command “metaprop_one” in Stata version 14.0 software (StataCorp LLC, College Station, TX, USA) as described previously . Forest plots were used to depict the study-specific proportions with 95% exact confidence intervals and overall pooled estimates with 95% Wald confidence intervals. The Chi-square statistic of the likelihood ratio test was used to identify the presence of significant heterogeneity when the P-value was less than 0.05. As the individual study weights were not calculated by the GLMMs, publication bias assessment was not performed in the present study. For the meta-analyses of proportion studies with low proportion outcomes, funnel plots were not used, as they are ineffective at detecting potential publication bias .
A total of 1113 candidate studies were identified through PubMed (379 studies), Web of Science (448 studies), and Scopus (286 studies). After 483 duplicates were removed, 630 studies were screened for titles and abstracts, and 207 studies remained for further full-text examination. The examination of the full text of the studies identified 19 studies [15,16,17, 19, 20, 23, 36,37,38,39,40,41,42,43,44,45,46,47,48] that met the eligibility criteria. One hundred and eighty-eight studies were excluded for specific reasons: 46 reviews, 31 mosquito identifications, 29 parasite studies, 28 assay developments, and 18 malaria biology in Anopheles mosquitoes; 11 for incomplete data; eight letters to the editor/comments/editorials; six genetic studies of malaria; three in vitro studies; three studies for which the full text was unavailable; two systematic reviews; one animal study; one case report; and one clinical trial. An additional 12 relevant studies [49,50,51,52,53,54,55,56,57,58,59,60] from the reference lists of the included studies and Google Scholar were examined for full texts. All of them met the eligibility criteria and were included in the systematic review. Overall, 31 studies [15,16,17, 19, 20, 23, 36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60] were included in the present systematic review (Fig. 1).
Characteristics of the studies
The characteristics of the included studies are shown in Table 1. The studies were published in 1970–2000 (13 studies, 41.9%), 2001–2010 (seve studies, 22.6%), and 2011–2021 (11 studies, 35.5%). The studies were conducted in western Thailand (14 studies, 45.2%), eastern Thailand (4 studies, 12.9%), northeastern Thailand (four studies, 12.9%), and other parts of Thailand (nine studies, 29.0%). Most of the studies were conducted in Tak (nine studies, 29%) and Kanchanaburi (three studies, 9.68%). For mosquito collection, the studies used human landing collection (17 studies, 54.8%), human landing collection/animal-baited collection (six studies, 19.4%), light trap (two studies, 6.5%), human landing collection/animal-baited collection/light trap (two studies, 6.5%), human landing collection/light trap (one study, 3.22%), animal-baited collection/light trap (one study, 3.22%), and prokopack aspiration (one study, 3.22%), and one study did not specify the mosquito collection method (one study, 3.22%). Most of the studies collected mosquitoes between 18:00 and 06:00 (48.4%) and between 18:00 and 24:00 (12.9%). The details of the studies are shown in Additional file 2: Table S2. The distributions of Plasmodium-positive mosquitoes are shown in Fig. 2.
Quality of the studies
The quality of the included studies was assessed using the STROBE checklist (Additional file 3: Table S3). For the included studies identified from the three databases, most of the studies (14/19, 73.7%) [20, 23, 37,38,39,40,41,42, 44,45,46,47,48, 51] were of moderate quality, whereas four (21.1%) were of high quality [15, 16, 36, 43] and one study (5.26%)  was of low quality. For the included studies identified from the reference lists and Google Scholar, most of the studies (9/12, 75%) [49, 51, 52, 54, 56,57,58,59,60] were of moderate quality, whereas three of them (25%) [50, 53, 55] were of high quality. Overall, seven studies (22.58%) were of high quality. Twenty-three studies (74.19%) were of moderate quality. Only one low-quality study , a short report, was included in the systematic review because it contained essential data on Plasmodium infection in primary Anopheles vectors.
Anopheles mosquitoes collected in the studies
The list of the Anopheles mosquitoes collected from the studies is shown in Table 2. The total number of Anopheles mosquitoes collected from all studies was 126,025. Of those Anopheles mosquitoes, 100,910 were identified to species. Most of the Anopheles species identified were An. minimus (40.16%), An. maculatus (16.59%), An. epiroticus (9.18%), An. sawadwongporni (6.35%), An. barbirostris (5.30%), An. dirus (4.91%), An. peditaeniatus (Leicester) (3.66%), An. aconitus (2.69%), An. philippinensis Ludlow (1.34%), An. nivipes (Theobald) (1.27%), An. campestris (1.26%), and An. karwari (James) (1.03%). The studies by Nosten et al.  and Brusich et al.  did not specify the Anopheles species collected in their study, whereas a study by Carrara et al.  collected only An. minimus, An. maculatus, and An. dirus. Another study by Sriwichai et al.  collected only An. minimus and An. maculatus.
Anopheles vectors for Plasmodium species
The prevalence of the Plasmodium species identified in all Anopheles vectors was estimated using the data from 31 studies [15,16,17, 19, 20, 23, 36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60]. This meta-analysis showed that the pooled prevalence of Plasmodium species identified in all Anopheles vectors among studies conducted during 1970–2000 was 1.7% (95% CI: 0.5–6.3%, Chi-square: 347.91, P < 0.001, 12 studies). The pooled prevalence of Plasmodium species identified in all Anopheles vectors among studies conducted during 2001–2010 was 0.11% (95% CI: 0.1–0.3%, Chi-square: 56, P < 0.001, 7 studies), and that among studies conducted during 2011–2021 was 0.2% (95% CI: 0.1–0.7%, Chi-square: 131, P < 0.001, 12 studies). Overall, the pooled prevalence of Plasmodium species identified in all Anopheles vectors was 0.5% (95% CI: 0.2–1.1%, Chi-square: 738.6, P < 0.001, 31 studies; Fig. 3).
The prevalence of Plasmodium species identified in the main Anopheles vectors (An. dirus, An. maculatus, and An. minimus) was estimated using the data from 22 studies [15,16,17, 19, 20, 38, 41,42,43,44,45, 47,48,49,50,51,52,53, 56,57,58,59]. The pooled prevalence of Plasmodium species identified in the main Anopheles vectors was 0.4% (95% CI: 0.3–7%, Chi-square: 299.2, P < 0.001, 22 studies; Fig. 4). The pooled prevalence of Plasmodium species identified in An. dirus was 2.1% (95% CI: 1.3–3.2%, Chi-square: 20.3, P < 0.001, nine studies; Fig. 5). The pooled prevalence of Plasmodium species identified in An. minimus was 0.4% (95% CI: 0.3–0.7%, Chi-square: 109.1, P < 0.001, 15 studies; Fig. 6). The pooled prevalence of Plasmodium species identified in An. maculatus was 0.4% (95% CI: 0.3–0.5%, Chi-square: 0.1, P: 0.385, 11 studies; Fig. 7).
Among the particular Anopheles species identified (n = 123,286), 566 Anopheles mosquitoes (0.46%) were positive for Plasmodium species The Anopheles species that tested positive for Plasmodium were Anopheles hodgkini Reid, Anopheles nigerrimus Giles, Anopheles balabacensis Baisas, An. dirus, Anopheles annularis van der Wulp, Anopheles kochi Dönitz, Anopheles vagus Dönitz, An. baimaii, Anopheles hyrcanus (Pallas), An. karwari, An. pseudowillmori, An. aconitus, An. sawadwongporni, An. barbirostris, An. minimus, Anopheles varuna Iyengar, An. maculatus, An. nivipes, An. peditaeniatus, An. campestris, An. philippinensis, and An. epiroticus (Table 3). The highest proportions of Plasmodium-positive mosquitoes were identified as An. hodgkini (2/6, 33.3%), An. nigerrimus (2/24, 8.33%), An. balabacensis (4/84, 4.76%), An. dirus (114/4956, 2.3%), An. annularis (16/852, 1.88%), An. kochi (8/519, 1.54%), An. vagus (3/215, 1.40%), and An. baimaii (1/86, 1.16%). Other Anopheles vectors for Plasmodium species present in a smaller proportion included An. karwari (0.67%), An. pseudowillmori (0.65%), An. aconitus (0.55%), An. sawadwongporni (0.47%), An. barbirostris (0.45%), An. minimus (0.38%), An. varuna (0.37%), An. maculatus (0.3%), An. hyrcanus (0.23%), An. nivipes (0.16%), An. peditaeniatus (0.16%), An. campestris (0.16%), An. philippinensis (0.15%), and An. epiroticus (0.1%). Some of the Anopheles species (n = 2739) were negative for Plasmodium species.
Proportion of Plasmodium species identified in each Anopheles vector
Table 4 shows the proportion of Plasmodium species identified in Anopheles vectors. High proportions of P. falciparum were identified in An. annularis (8.50%), An. kochi (6.03%), and An. barbirostris (5.03%). High proportions of P. vivax were identified in An. hodgkini (33.3%), An. aconitus (28.6%), and An. nivipes (25.0%). High proportions of P. falciparum/P. vivax were identified in An. nigerrimus (14.3%), An. vagus (2.88%), and An. peditaeniatus (2.47%). High proportions of P. falciparum/P. vivax/mixed infection were identified in An. dirus (7.27%). High proportions of Plasmodium spp. (species not identified) were An. balabacansis (4.76%) and An. dirus (4.07%).
The present systematic review showed that An. minimus (40.16%), An. maculatus (16.59%), and An. epiroticus (9.18%) were the most common Anopheles mosquitoes identified in the studies included. The differences in the proportions of mosquito species might be due to the variations in the local environment, such as differences in the study site, seasonal, biology, and behavior of each species. Although An. minimus, An. maculatus, and An. epiroticus were the main mosquitoes identified in the studies, these mosquitoes harbored Plasmodium spp. at only 0.38%, 0.3%, and 0.1%, respectively. However, it is well recognized that several of the important malaria vectors, as well as other Anopheles species in Thailand, are members of closely related sibling species. Thus, entomologists classify them into complexes or groups . Although illustrated morphological keys have been published for the identification of both female adult and larval stages of Anopheles mosquitoes in Thailand , the identification of Anopheles specimens based exclusively on morphological characteristics is questionable and potentially leads to misidentification .
The three main malaria vectors in Thailand are An. minimus, An. maculatus, and An. dirus, as previously reported by the Division of Vector Borne Disease, Ministry of Public Health, Thailand . The present pooled analyses showed that 0.4%, 0.4%, and 2.1%, respectively, of these Anopheles mosquitoes harbored Plasmodium parasites. Therefore, meta-analysis in the present study suggested that there was a very low prevalence of Plasmodium spp. in the main Anopheles mosquitoes in Thailand. The success of the implementation of vector control tools, such as LLINs and IRS, across Thailand, has significantly reduced malaria transmission from 149,586 cases in 2000  to 4473 cases in 2020 , thus possibly contributing to the low prevalence of Plasmodium in field populations. Other factors such as a change in behavior of Anopheles mosquitoes due to elimination efforts , climate change, and urbanization  might also have affected the EIR in terms of both human biting rate and SRs. To determine the annual EIR and to consequently evaluate vector control interventions, we recommend that a relatively large number of mosquito specimens (more than 1000 mosquitoes) are collected for an accurate analysis of SRs in each locality, preferably as a part of a long-term study. To achieve malaria eradication in Thailand by 2024 , further Anopheles studies should be performed in provinces along the international borders, such as the Tak, Kanchanaburi, Chanthaburi, Ubon Ratchathani, and Yala provinces, in which the disease is endemic (Fig. 2) [5, 65].
To identify where malaria transmission occurs, reliable estimations of the proportions of infective Anopheles mosquitoes, as reflected by the presence of sporozoites in the salivary gland, are needed [66, 67]. Several techniques have been used to quantify sporozoites in mosquitoes. Dissection and microscopic examination of the salivary glands is considered the “gold standard” method, but this approach is labor intensive and impractical in the field . Enzyme-linked immunosorbent assay to detect circumsporozoite proteins (CSP-ELISA) is another widely used technique . However, CSP-ELISA has been shown to give false positive results, thus overestimating the real SR . Molecular-based methods have been developed to improve sensitivity and specificity . In our analysis, 19 studies (61.30%) used CSP-ELISA as a means of Plasmodium detection. Six studies (19.35%) used dissection of the salivary gland and gut. It is possible that the limitations of the dissection and CSP-ELISA methods used in the studies could have affected the estimation of sporozoite infection. Six studies (19.35%) used polymerase chain reaction (PCR)-based techniques. Sumruayphol et al.  performed nested PCR and real-time PCR on 9260 An. epiroticus specimens and found only six mosquitoes infected with P. falciparum and three with P. vivax. In another study, Tainchum et al.  also used real-time PCR for Plasmodium detection in 1090 An. minimus specimens and found only one positive sample. These studies suggested that even using methods with high sensitivity, the prevalence of Plasmodium species in Anopheles mosquitoes was very low, suggesting that the prevalence of parasite vectors in Thailand was genuinely low. Other techniques such as rapid dipstick immunochromatographic assays (Vec-Test™ Malaria)  and near-infrared spectroscopy  have also been developed for Plasmodium detection in Anopheles mosquitoes. Overall, to avoid the overestimation of SR and the EIR, it is highly recommended that all positive CSP-ELISA samples be reanalyzed or the results confirmed by performing Plasmodium-specific PCR .
Anopheles dirus had the highest pooled prevalence of Plasmodium species identified (2.1%) in our analysis. However, it should also be noted that the high prevalence of infection does not necessarily translate to the species being the main vector. Other factors also play a crucial role in the importance of primary vectors, for instance, the species must often be abundantly present and prefer to feed on humans . Previous studies also reported relatively low numbers of collected An. dirus specimens (ranging from 10–78 mosquitoes/location) recently [41, 42, 45,46,47, 53]. Therefore, it is vital to assess transmission indicators (e.g., EIR) to determine the importance of each vector species. The prevalence of Plasmodium species identified in An. dirus decreased from 5% in 1987  to 1% in 1990 . However, the prevalence of Plasmodium species identified in An. dirus increased to 4% in 1990  and decreased to 1–3% during 1991–2017 [16, 19, 20, 41, 42]. The yearly trend results were heterogeneous, and differences in study sites had to be considered. The pooled prevalence of Plasmodium species identified in An. minimus and An. maculatus was the same at 0.4%. These results suggest that the likelihood of finding an infected wild An. minimus or An. maculatus is lower than that of finding an infected An. dirus. The reasons for this observation remain to be investigated. However, several factors influencing vectorial capacity and competence have been documented, including mosquito longevity, the duration of sporogonic development, and the susceptibility or resistance of the vector to Plasmodium . A previous study of three laboratory strains of An. dirus, An. minimus, and An. sawadwongporni showed similar susceptibility to P. vivax infection using an artificial feeding system . These laboratory-raised mosquitoes are highly inbred and may be genetically dissimilar to the originally sampled population . In our analysis, natural Plasmodium infection in wild Anopheles mosquitoes was also considerably different among populations and species. In another study, large differences in P. falciparum infection were observed in a wild population of Anopheles gambiae Giles in West Africa . Therefore, it is not surprising to find differences in Plasmodium prevalence in the diverse field populations in Thailand.
The present study had some limitations. First, the studies that were included for systematic review were not performed in all areas in which malaria cases have been reported. There are missing mosquito data from the Thailand-Malaysia border, where high levels of malaria have been reported. Hence, the systematic review did not represent the overall prevalence in Thailand. Second, the majority of the studies (14 studies, 45.16%) used in our analysis used only morphological keys for species identification, and an additional 12 studies (38.71%) did not indicate which identification method was used. Therefore, some of the Anopheles specimens in these studies were only classified into complexes or groups. We therefore simply reported and analyzed the species data using the information presented in the studies. Only five studies (16.13%) used molecular techniques to confirm the species of the Anopheles specimens after the initial morphological identification into complexes or groups. To reflect the true prevalence of Plasmodium in each Anopheles species, particularly the primary and secondary malaria vectors, species confirmation using molecular techniques should also be performed. Third, it has been demonstrated that a positive SR may coincide with peak mosquito populations [16, 20, 48]. Thus, the prevalence of Plasmodium infection in each mosquito species could also be varied depending on season and mosquito abundance. In our analysis, only cross-sectional studies were included and we did not attempt to factor seasonal variation, as some included studies did not report seasonal information of positive Plasmodium infection specimens. However, is it possible that in any of the included studies, data for one species might have been collected during high transmission, whereas data for other species might have been collected during low transmission seasons, which might be a source of overestimation and underestimation of the importance of different vectors. Fourth, we included studies that identified Anopheles-harbored sporozoites and also oocysts of Plasmodium spp. As these stages are different indicators and may have different interpretations, we used the oocyst formation rate to study infection in the vector, that is, to show the susceptibility of the vector to infection; however, it does not indicate the importance of the vector in transmission. Therefore, the prevalence of Plasmodium spp. in Anopheles mosquitoes indicates the infection rates rather than the transmission capability. Finally, there are few synthesized studies, which are then subdivided into smaller subgroups for the purpose of comparing differences in mosquito species.
This systematic review confirmed the relatively low prevalence of Plasmodium species in wild Anopheles mosquitoes in Thailand. Anopheles dirus was likely to be the predominant species harboring Plasmodium species. However, the measurement of the number of sporozoites must be performed with caution to avoid overestimating the extent of Plasmodium infection. An accurate estimation of the EIR using a standardized parasite detection technique would also require the use of a relatively large number of Anopheles specimens to assess the impact of vector control interventions. The results of the present study also serve to identify potential vectors for malaria as a basis for further detailed studies. With this information, more intensive mosquito studies should be undertaken in several areas of Thailand to explore the prevalence of Plasmodium species in Anopheles mosquitoes as a basis for the development of vector control strategies.
Availability of data and materials
All data relating to the present study are available in this manuscript.
Entomological inoculation rate
Generalized linear mixed models
WHO. World malaria report 2021. Geneva: World Health Organization; 2021.
Cui L, Cao Y, Kaewkungwal J, Khamsiriwatchara A, Lawpoolsri S, Soe TN, et al. Malaria elimination in the Greater Mekong Subregion: challenges and prospects. In: Manguin S, Dev V, editors., et al., Towards malaria elimination - A leap forward. London: IntechOpen; 2018.
Corbel V, Nosten F, Thanispong K, Luxemburger C, Kongmee M, Chareonviriyaphap T. Challenges and prospects for dengue and malaria control in Thailand. Southeast Asia Trends Parasitol. 2013;29:623–33.
Suwonkerd W, Ritthison W, Ngo CT, Tainchum K, Bangs MJ, Chareonviriyaphap T. Vector biology and malaria transmission in Southeast Asia, Anopheles mosquitoes. In: Manguin S, editor. New insights into malaria vectors. London: IntechOpen; 2013.
Bureau of Vector-Borne Diseases. Annual Report 2020. Nonthaburi: Department of Disease Control, Ministry of Public Health; 2020.
WHO. World Malaria Report 2018. Geneva: World Health Organization; 2018.
Bureau of Vector-Borne Diseases. Guide to malaria elimination for Thailand’s local administrative organizations and the health network. Nonthaburi: Department of Disease Control Ministry of Public Health; 2019.
Durnez L, Coosemans M. Residual transmission of malaria: an old issue for new approaches. London: InTechOpen; 2013.
Tainchum K, Kongmee M, Manguin S, Bangs MJ, Chareonviriyaphap T. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 2015;31:109–19.
Tananchai C, Manguin S, Bangs MJ, Chareonviriyaphap T. Malaria vectors and species complexes in Thailand: Implications for vector control. Trends Parasitol. 2019;35:544–58.
Rattanarithikul R, Harrison BA, Harbach RE, Panthusiri P, Coleman RE. Illustrated keys to the mosquitoes of Thailand. IV. Anopheles. Southeast Asian J Trop Med Public Health. 2006;37:1–128.
Sinka ME, Bangs MJ, Manguin S, Chareonviriyaphap T, Patil AP, Temperley WH, et al. The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2011;4:89.
Harbach RE. The phylogeny and classification of Anopheles. In: Manguin S, editor. Anopheles mosquitoes - New insights into malaria vectors. London: IntechOpen; 2013.
Harbach RE, Rattanarithikul R, Harrison BA. Anopheles prachongae, a new species of the Gigas Complex of subgenus Anopheles (Diptera: Culicidae) in Thailand, contrasted with known forms of the complex. Zootaxa. 2017;4236:65.
Rosenberg R, Andre RG, Somchit L. Highly efficient dry season transmission of malaria in Thailand. Trans R Soc Trop Med Hyg. 1990;84:22–8.
Rattanarithikul R, Konishi E, Linthicum KJ. Detection of Plasmodium vivax and Plasmodium falciparum circumsporozoite antigen in anopheline mosquitoes collected in southern Thailand. Am J Trop Med Hyg. 1996;54:114–21.
Coleman RE, Sithiprasasna R, Kankaew P, Kiaattiut C, Ratanawong S, Khuntirat B, et al. Naturally occurring mixed infection of Plasmodium vivax VK210 and P. vivax VK247 in Anopheles mosquitoes (Diptera: Culicidae) in western Thailand. J Med Entomol. 2002;39:556–9.
Gould DJ, Esah S, Pranith U. Relation of Anopheles aconitus to malaria transmission in the central plain of Thailand. Trans R Soc Trop Med Hyg. 1967;61:441–2.
Green CA, Rattanarithikul R, Pongparit S, Sawadwongporn P, Baimai V. A newly-recognized vector of human malarial parasites in the Oriental region, Anopheles (Cellia) pseudowillmori (Theobald, 1910). Trans R Soc Trop Med Hyg. 1991;85:35–6.
Somboon P, Aramrattana A, Lines J, Webber R. Entomological and epidemiological investigations of malaria transmission in relation to population movements in forest areas of north-west Thailand. Southeast Asian J Trop Med Public Health. 1998;29:3–9.
Thongsahuan S, Baimai V, Junkum A, Saeung A, Min GS, Joshi D, et al. Susceptibility of Anopheles campestris-like and Anopheles barbirostris species complexes to Plasmodium falciparum and Plasmodium vivax in Thailand. Mem Inst Oswaldo Cruz. 2011;106:105–12.
Taai K, Harbach RE. Systematics of the Anopheles barbirostris species complex (Diptera: Culicidae: Anophelinae) in Thailand. Zool J Linn Soc. 2015;174:244–64.
Sumruayphol S, Apiwathnasorn C, Komalamisra N, Ruangsittichai J, Samung Y, Chavalitshewinkoon-Petmitr P. Bionomic status of Anopheles epiroticus Linton & Harbach, a coastal malaria vector, in Rayong Province, Thailand. Southeast Asian J Trop Med Public Health. 2010;41:541–7.
Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53.
Tusting LS, Bousema T, Smith DL, Drakeley C. Measuring changes in Plasmodium falciparum transmission: precision, accuracy and costs of metrics. Adv Parasitol. 2014;84:151–208.
Shaukat AM, Breman JG, McKenzie FE. Using the entomological inoculation rate to assess the impact of vector control on malaria parasite transmission and elimination. Malar J. 2010;9:122.
Graumans W, Jacobs E, Bousema T, Sinnis P. When is a Plasmodium-infected mosquito an infectious mosquito? Trends Parasitol. 2020;36:705–16.
Stevenson JC, Norris DE. Implicating cryptic and novel Anophelines as malaria vectors in Africa. Insects. 2016;8:1.
Epopa PS, Collins CM, North A, Millogo AA, Benedict MQ, Tripet F, et al. Seasonal malaria vector and transmission dynamics in Western Burkina Faso. Malar J. 2019;18:113.
Stone WJ, Eldering M, van Gemert GJ, Lanke KH, Grignard L, van de Vegte-Bolmer MG, et al. The relevance and applicability of oocyst prevalence as a read-out for mosquito feeding assays. Sci Rep. 2013;3:3418.
Moher D, Liberati A, Tetzlaff J, Altman DG, Group P. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med. 2009;6:e1000097.
von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP, et al. The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies. PLoS Med. 2007;4:e296.
DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7:177–88.
Nyaga VN, Arbyn M, Aerts M. Metaprop: a Stata command to perform meta-analysis of binomial data. Arch Public Health. 2014;72:39.
Hunter JP, Saratzis A, Sutton AJ, Boucher RH, Sayers RD, Bown MJ. In meta-analyses of proportion studies, funnel plots were found to be an inaccurate method of assessing publication bias. J Clin Epidemiol. 2014;67:897–903.
Brusich M, Grieco J, Penney N, Tisgratog R, Ritthison W, Chareonviriyaphap T, et al. Targeting educational campaigns for prevention of malaria and dengue fever: an assessment in Thailand. Parasit Vectors. 2015;8:43.
Frances SP, Klein TA, Wirtz RA, Eamsila C, Pilakasiri C, Linthicum KJ. Plasmodium falciparum and P. vivax circumsporozoite proteins in anophelines (Diptera: Culicidae) collected in eastern Thailand. J Med Entomol. 1996;33:990–1.
Gingrich JB, Weatherhead A, Sattabongkot J, Pilakasiri C, Wirtz RA. Hyperendemic malaria in a Thai village: dependence of year-round transmission on focal and seasonally circumscribed mosquito (Diptera: Culicidae) habitats. J Med Entomol. 1990;27:1016–26.
Harbach RE, Gingrich JB, Pang LW. Some entomological observations on malaria transmission in a remote village in Northwestern Thailand. J Eur Mosq Control Assoc. 1987;3:296–301.
Limrat D, Rojruthai B, Apiwathnasorn C, Samung Y, Prommongkol S. Anopheles barbirostris/campestris as a probable vector of malaria in Aranyaprathet, Sa Kaeo Province. Southeast Asian J Trop Med Public Health. 2001;32:739–44.
Poolphol P, Harbach RE, Sriwichai P, Aupalee K, Sattabongkot J, Kumpitak C, et al. Natural Plasmodium vivax infections in Anopheles mosquitoes in a malaria endemic area of Northeastern Thailand. Parasitol Res. 2017;116:3349–59.
Rattaprasert P, Chaksangchaichot P, Wihokhoen B, Suparach N, Sorosjinda-Nunthawarasilp P. Detection of putative antimalarial-resistant Plasmodium vivax in Anopheles vectors at Thailand-Cambodia and Thailand-Myanmar borders. Southeast Asian J Trop Med Public Health. 2016;47:182–93.
Sattabongkot J, Kiattibut C, Kumpitak C, Ponlawat A, Ryan JR, Chan AST, et al. Evaluation of the VecTest Malaria Antigen Panel assay for the detection of Plasmodium falciparum and P. vivax circumsporozoite protein in anopheline mosquitoes in Thailand. J Med Entomol. 2004;41:209–14.
Sriwichai P, Karl S, Samung Y, Kiattibutr K, Sirichaisinthop J, Mueller I, et al. Imported Plasmodium falciparum and locally transmitted Plasmodium vivax: cross-border malaria transmission scenario in Northwestern Thailand. Malar J. 2017;16:258.
Sriwichai P, Samung Y, Sumruayphol S, Kiattibutr K, Kumpitak C, Payakkapol A, et al. Natural human Plasmodium infections in major Anopheles mosquitoes in Western Thailand. Parasit Vectors. 2016;9:17.
Sumarnrote A, Corbel V, Overgaard HJ, Celhay O, Marasri N, Fustec B, et al. Plasmodium infections in Anopheles mosquitoes in Ubon Ratchathani Province, Northeastern Thailand during a malaria outbreak. J Am Mosq Control Assoc. 2018;34:11–7.
Tainchum K, Ritthison W, Chuaycharoensuk T, Bangs MJ, Manguin S, Chareonviriyaphap T. Diversity of Anopheles species and trophic behavior of putative malaria vectors in two malaria endemic areas of Northwestern Thailand. J Vector Ecol. 2014;39:424–36.
Zollner G, Sattabongkot J, Vaughan JA, Kankaew P, Robert LL, Thimasarn K, et al. Longitudinal evaluation of malaria epidemiology in an isolated village in Western Thailand: I. Study site and adult anopheline bionomics. Southeast Asian J Trop Med Public Health. 2016;47:341–65.
Baker EZ, Beier JC, Meek SR, Wirtz RA. Detection and quantification of Plasmodium falciparum and P. vivax infections in Thai-Kampuchean Anopheles (Diptera: Culicidae) by enzyme-linked immunosorbent assay. J Med Entomol. 1987;24:536–41.
Carrara VI, Sirilak S, Thonglairuam J, Rojanawatsirivet C, Proux S, Gilbos V, et al. Deployment of early diagnosis and mefloquine-artesunate treatment of falciparum malaria in Thailand: the Tak Malaria Initiative. PLoS Med. 2006;3:e183.
Coleman RE, Barth JF, Turell MJ, Gordon SW, Sattabongkot J, Copeland R, et al. Development and evaluation of a dipstick assay for detection of Plasmodium falciparum and P. vivax sporozoites in mosquitoes (Diptera: Culicidae). J Med Entomol. 2000;37:581–7.
Eamkum P, Sungvornyothin S, Kritpetcharat O, Daduang J, Lek-Uthai U, Charerntanyarak L, et al. A single-round multiplex PCR assay for the identification of Anopheles minimus related species infected with Plasmodium falciparum and Plasmodium vivax. Parasitol Int. 2014;63:442–9.
Edwards HM, Sriwichai P, Kirabittir K, Prachumsri J, Chavez IF, Hii J. Transmission risk beyond the village: entomological and human factors contributing to residual malaria transmission in an area approaching malaria elimination on the Thailand-Myanmar border. Malar J. 2019;18:221.
Junkum A, Pitasawat B, Tuetun B, Saeung A, Rattanachanpichai E, Jariyapan N, et al. Seasonal abundance and biting activity of Anopheles Aconitus (Diptera: Culicidae) in Chiang Mai, northern Thailand. Southeast Asian J Trop Med Public Health. 2007;38:215–23.
Nosten F, van Vugt M, Price R, Luxemburger C, Thway KL, Brockman A, et al. Effects of artesunate-mefloquine combination on incidence of Plasmodium falciparum malaria and mefloquine resistance in western Thailand: a prospective study. Lancet. 2000;356:297–302.
Ratanatham S, Upatham ES, Prasittisuk C, Rojanasunan W, Theerasilp N, Tremongkol A, et al. Bionomics of Anopheles minimus and its role in malaria transmission in Thailand. Southeast Asian J Trop Med Public Health. 1988;19:283–9.
Sithiprasasna R, Linthicum KJ, Liu GJ, Jones JW, Singhasivanon P. Some entomological observations on temporal and spatial distribution of malaria vectors in three villages in northwestern Thailand using a geographic information system. Southeast Asian J Trop Med Public Health. 2003;34:505–16.
Sumruayphol S, Chaiphongpachara T, Samung Y, Ruangsittichai J, Cui L, Zhong D, et al. Seasonal dynamics and molecular differentiation of three natural Anopheles species (Diptera: Culicidae) of the Maculatus group (Neocellia series) in malaria hotspot villages of Thailand. Parasit Vectors. 2020;13:574.
Upatham ES, Prasittisuk C, Ratanatham S, Green CA, Rojanasunan W, Setakana P, et al. Bionomics of Anopheles maculatus complex and their role in malaria transmission in Thailand. Southeast Asian J Trop Med Public Health. 1988;19:259–69.
Wilkinson R, Miller TA, Esah S. Anthropophilic mosquitoes in central Thailand, with notes on Anopheles balabacensis Baisas and malaria. Mosq News. 1970;30:146–8.
Somboon P, Rattanarithikul R. Mosquito surveys, rearing, preservation of mosquito specimens and identification of Anopheles in Thailand 2013. Bureau of Vector Borne Disease: Ministry of Public Health; 2013.
Bureau of Vector-Borne Diseases. Annual Report 2005. Nonthaburi: Department of Disease Control, Ministry of Public Health; 2005.
Manguin S, Garros C, Dusfour I, Harbach RE, Coosemans M. Bionomics, taxonomy, and distribution of the major malaria vector taxa of Anopheles subgenus Cellia in Southeast Asia: an updated review. Infect Genet Evol. 2008;8:489–503.
Warrell DA, Gilles HM. Essential Malariology. London: CRC Press; 2002.
Sudathip P, Kitchakarn S, Shah JA, Bisanzio D, Young F, Gopinath D, et al. A foci cohort analysis to monitor successful and persistent foci under Thailand’s Malaria elimination strategy. Malar J. 2021;20:118.
Hendershot AL, Esayas E, Sutcliffe AC, Irish SR, Gadisa E, Tadesse FG, et al. A comparison of PCR and ELISA methods to detect different stages of Plasmodium vivax in Anopheles arabiensis. Parasit Vectors. 2021;14:473.
Lardeux F, Tejerina R, Aliaga C, Ursic-Bedoya R, Lowenberger C, Chavez T. Optimization of a semi-nested multiplex PCR to identify Plasmodium parasites in wild-caught Anopheles in Bolivia, and its application to field epidemiological studies. Trans R Soc Trop Med Hyg. 2008;102:485–92.
Beier JC, Perkins PV, Koros JK, Onyango FK, Gargan TP, Wirtz RA, et al. Malaria sporozoite detection by dissection and ELISA to assess infectivity of Afrotropical Anopheles (Diptera: Culicidae). J Med Entomol. 1990;27:377–84.
Wirtz RA, Zavala F, Charoenvit Y, Campbell GH, Burkot TR, Schneider I, et al. Comparative testing of monoclonal antibodies against Plasmodium falciparum sporozoites for ELISA development. Bull World Health Organ. 1987;65:39–45.
Durnez L, Van Bortel W, Denis L, Roelants P, Veracx A, Trung HD, et al. False positive circumsporozoite protein ELISA: a challenge for the estimation of the entomological inoculation rate of malaria and for vector incrimination. Malar J. 2011;10:195.
Bass C, Nikou D, Blagborough AM, Vontas J, Sinden RE, Williamson MS, et al. PCR-based detection of Plasmodium in Anopheles mosquitoes: a comparison of a new high-throughput assay with existing methods. Malar J. 2008;7:177.
Appawu MA, Bosompem KM, Dadzie S, McKakpo US, Anim-Baidoo I, Dykstra E, et al. Detection of malaria sporozoites by standard ELISA and VecTestTM dipstick assay in field-collected anopheline mosquitoes from a malaria endemic site in Ghana. Trop Med Int Health. 2003;8:1012–7.
Maia MF, Kapulu M, Muthui M, Wagah MG, Ferguson HM, Dowell FE, et al. Detection of Plasmodium falciparum infected Anopheles gambiae using near-infrared spectroscopy. Malar J. 2019;18:85.
Cohuet A, Harris C, Robert V, Fontenille D. Evolutionary forces on Anopheles: what makes a malaria vector? Trends Parasitol. 2010;26:130–6.
Zollner GE, Ponsa N, Garman GW, Poudel S, Bell JA, Sattabongkot J, et al. Population dynamics of sporogony for Plasmodium vivax parasites from western Thailand developing within three species of colonized Anopheles mosquitoes. Malar J. 2006;5:68.
Norris DE, Shurtleff AC, Touré YT, Lanzaro GC. Microsatellite DNA polymorphism and heterozygosity among field and laboratory populations of Anopheles gambiae ss (Diptera: Culicidae). J Med Entomol. 2001;38:336–40.
Riehle MM, Markianos K, Niaré O, Xu J, Li J, Touré AM, et al. Natural malaria infection in Anopheles gambiae is regulated by a single genomic control region. Science. 2006;312:577–9.
This research was financially supported by the new strategic research project (P2P) fiscal year 2022, Walailak University, Thailand.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Additional file 2: Table S2.
Details of the included studies
Additional file 3: Table S3.
Quality of the included studies
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sukkanon, C., Masangkay, F.R., Mala, W. et al. Prevalence of Plasmodium spp. in Anopheles mosquitoes in Thailand: a systematic review and meta-analysis. Parasites Vectors 15, 285 (2022). https://doi.org/10.1186/s13071-022-05397-2