Plasmodium Infections in Anopheles Mosquitoes in Ubon Ratchathani Province, Northeastern Thailand During a Malaria Outbreak

Assessing receptivity to malaria using case surveillance and forest data in a near-elimination setting in northeast Thailand

BACKGROUND

Thailand aimed to eliminate malaria by 2024, and as such is planning for future prevention of re-establishment in malaria free provinces. Understanding the receptivity of local areas to malaria allows the appropriate targeting of interventions. Current approaches to assessing receptivity involve collecting entomological data. Forest coverage is known to be associated with malaria risk, as an environment conducive to both vector breeding and high-risk human behaviours.

METHODS

Geolocated, anonymized, individual-level surveillance data from 2011 to 2021 from the Thai Division of Vector-Borne Disease (DVBD) was used to calculate incidence and estimated Rc at village level. Forest cover was calculated using raster maps of tree crown cover density and year of forest loss from the publicly available Hansen dataset. Incidence and forest cover were compared graphically and using Spearman's rho. The current foci classification system was applied to data from the last 5 years (2017-2021) and forest cover for 2021 compared between the classifications. A simple risk score was developed to identify villages with high receptivity.

RESULTS

There was a non-linear decrease in annual cases by 96.6% (1061 to 36) across the two provinces from 2011 to 2021. Indigenous Annual Parasite Index (API) and approximated Rc were higher in villages in highly forested subdistricts, and with higher forest cover within 5 km. Forest cover was also higher in malaria foci which consistently reported malaria cases each year than those which did not. An Rc > 1 was only reported in villages in subdistricts with > 25% forest cover. When applying a simple risk score using forest cover and recent case history, the classifications were comparable to those of the risk stratification system currently used by the DVBD.

CONCLUSIONS

There was a positive association between forest coverage around a village and indigenous malaria cases. Most local transmission was observed in the heavily forested subdistricts on the international borders with Laos and Cambodia, which are where the most receptive villages are located. These areas are at greater risk of importation of malaria due to population mobility and forest-going activities. Combining forest cover and recent case surveillance data with measures of vulnerability may be useful for prediction of malaria recurrence risk.

Mapping malaria transmission foci in Northeast Thailand from 2011 to 2021: approaching elimination in a hypoendemic area

BACKGROUND

Thailand is approaching local elimination of malaria in the eastern provinces. It has successfully reduced the number of cases over the past decade, but there are persistent transmission hot spots in and around forests. This study aimed to use data from the malaria surveillance system to describe the spatiotemporal trends of malaria in Northeast Thailand and fine-scale patterns in locally transmitted cases between 2011 and 2021.

METHODS

Case data was stratified based on likely location of infection and parasite species. Annual Parasite Index per 1000 population (API) was calculated for different categories. Time series decomposition was performed to identify trends and seasonal patterns. Statistically significant clusters of high (hot spots) and low (cold spots) API were identified using the Getis-Ord Gi* statistic. The stability of those hot spots and the absolute change in the proportion of API density from baseline were compared by case type.

RESULTS

The total number of confirmed cases experienced a non-linear decline by 96.6%, from 1061 in 2011 to 36 in 2021. There has been a decline in both Plasmodium vivax and Plasmodium falciparum case numbers, with only four confirmed P. falciparum cases over the last two years-a 98.89% drop from 180 in 2011. API was generally higher in Si Sa Ket province, which had peaks every 2-3 years. There was a large outbreak in Ubon Ratchathani in 2014-2016 which had a high proportion of P. falciparum reported. The proportion of cases classified increased over the study period, and the proportion of cases classed as indigenous to the village of residence increased from 0.2% to 33.3%. There were stable hot spots of indigenous and imported cases in the south of Si Sa Ket and southeast of Ubon Ratchathani. Plasmodium vivax hot spots were observed into recent years, while those of P. falciparum decreased to zero in Ubon in 2020 and emerged in the eastern part in 2021, the same year that P. falciparum hot spots in Si Sa Ket reached zero.

CONCLUSIONS

There has been a large, non-linear decline in the number of malaria cases reported and an increasing proportion of cases are classed as indigenous to the patient's village of residence. Stable hot spots of ongoing transmission in the forested border areas were identified, with transmission likely persisting because of remote location and high-risk forest-going behaviours. Future efforts should include cross-border collaboration and continued targeting of high-risk behaviours to reduce the risk of imported cases seeding local transmission.

Bacterial community structure of Anopheles hyrcanus group, Anopheles nivipes, Anopheles philippinensis, and Anopheles vagus from a malaria-endemic area in Thailand

Bacterial content of mosquitoes has given rise to the development of innovative tools that influence and seek to control malaria transmission. This study identified the bacterial microbiota in field-collected female adults of the Anopheles hyrcanus group and three Anopheles species, Anopheles nivipes, Anopheles philippinensis, and Anopheles vagus, from an endemic area in the southeastern part of Ubon Ratchathani Province, northeastern Thailand, near the Lao PDR-Cambodia-Thailand border. A total of 17 DNA libraries were generated from pooled female Anopheles abdomen samples (10 abdomens/ sample). The mosquito microbiota was characterized through the analysis of DNA sequences from the V3-V4 regions of the 16S rRNA gene, and data were analyzed in QIIME2. A total of 3,442 bacterial ASVs were obtained, revealing differences in the microbiota both within the same species/group and between different species/group. Statistical difference in alpha diversity was observed between An. hyrcanus group and An. vagus and between An. nivipes and An. vagus, and beta diversity analyses showed that the bacterial community of An. vagus was the most dissimilar from other species. The most abundant bacteria belonged to the Proteobacteria phylum (48%-75%) in which Pseudomonas, Serratia, and Pantoea were predominant genera among four Anopheles species/group. However, the most significantly abundant genus observed in each Anopheles species/group was as follows: Staphylococcus in the An. hyrcanus group, Pantoea in the An. nivipes, Rosenbergiella in An. philippinensis, and Pseudomonas in An. vagus. Particularly, Pseudomonas sp. was highly abundant in all Anopheles species except An. nivipes. The present study provides the first study on the microbiota of four potential malaria vectors as a starting step towards understanding the role of the microbiota on mosquito biology and ultimately the development of potential tools for malaria control.

Prevalence of Plasmodium spp. in Anopheles mosquitoes in Thailand: a systematic review and meta-analysis

Background

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.

Methods

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.

Results

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%.

Conclusions

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.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05397-2.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

Trophic Behavior and Species Diversity of the Anopheles barbirostris Complex (Diptera: Culicidae) in Thailand

Species of the Anopheles barbirostris complex (Myzorhynchus Series of the subgenus Anopheles) are potential vectors of malaria and filariasis parasites. Owing to the lack of reliable identification methods, the biting activity and host preference of the species within this complex have not been previous described. In this study, the trophic behavior and host preferences of the species in the complex were determined in Thailand, and a map of their geographical distributions constructed. Adult female mosquitoes in Chanthaburi and Prachuap Khiri Khan provinces were collected using human landings (indoor and outdoor) and cow-bait. Morphological characters were used to initially identify the mosquitoes, and then, species were further separated by multiplex PCR. Four species were identified in the An. barbirostris complex: An. wejchoochotei Taai & Harbach (Diptera: Culicidae), An. saeungae Taai & Harbach (Diptera: Culicidae), An. dissidens Taai & Harbach (Diptera: Culicidae), and An. campestris Reid (Diptera: Culicidae). The species showed exophagic and zoophilic behavior, except An. wejchoochotei in Chanthaburi, which displayed slight endophilic behavior. Indoor and outdoor blood feeding behavior of An. wejchoochotei and An. saeungae peaked at 2000 h. An updated distribution map of the An. barbirostris complex in Thailand is presented. With a better understanding of the sibling vector species and their spatial distributions, more specific and effective control strategies for vectors of malaria and filariasis in Thailand can be achieved.

Anopheles ecology, genetics and malaria transmission in northern Cambodia

In the Greater Mekong Subregion, malaria cases have significantly decreased but little is known about the vectors or mechanisms responsible for residual malaria transmission. We analysed a total of 3920 Anopheles mosquitoes collected during the rainy and dry seasons from four ecological settings in Cambodia (villages, forested areas near villages, rubber tree plantations and forest sites). Using odor-baited traps, 81% of the total samples across all sites were collected in cow baited traps, although 67% of the samples attracted by human baited traps were collected in forest sites. Overall, 20% of collected Anopheles were active during the day, with increased day biting during the dry season. 3131 samples were identified morphologically as 14 different species, and a subset was also identified by DNA amplicon sequencing allowing determination of 29 Anopheles species. The investigation of well characterized insecticide mutations (ace-1, kdr, and rdl genes) indicated that individuals carried mutations associated with response to all the different classes of insecticides. There also appeared to be a non-random association between mosquito species and insecticide resistance with Anopheles peditaeniatus exhibiting nearly fixed mutations. Molecular screening for Plasmodium sp. presence indicated that 3.6% of collected Anopheles were positive, most for P. vivax followed by P. falciparum. These results highlight some of the key mechanisms driving residual human malaria transmission in Cambodia, and illustrate the importance of diverse collection methods, sites and seasons to avoid bias and better characterize Anopheles mosquito ecology in Southeast Asia.

New insights into malaria vector bionomics in Lao PDR: a nationwide entomology survey

BACKGROUND

In Laos, the malaria burden remains high despite a significant reduction of cases during the last decade. In the context of the disease elimination by 2030, a nationwide entomological survey was conducted to better understand the distribution, abundance and behaviour of major malaria vectors (Anopheles spp.) in the country.

METHODS

Mosquito collections were implemented in ten villages from ten provinces during the rainy and dry seasons of 2014 and 2015 by using human landing catch (HLC) and cow bait collection (CBC) methods. After morphological identification in the field, molecular identification of the sibling species of Anopheles mosquitoes from the Funestus, Leucosphyrus, and Maculatus groups were determined using PCR specific alleles. A screening of Plasmodium falciparum and Plasmodium vivax infections in the vectors was carried out by quantitative PCR assays.

RESULTS

A total of 14,146 adult mosquitoes representing 25 different Anopheles species were collected and morphologically identified. Molecular identification revealed the presence of 12 sibling species within the main primary vector groups, including Anopheles maculatus, Anopheles rampae, Anopheles sawadwongporni, Anopheles pseudowillmori, Anopheles dravidicus, Anopheles minimus, Anopheles aconitus, Anopheles pampanai, Anopheles harrisoni, Anopheles dirus, Anopheles baimaii, Anopheles nemophilous. Anopheles maculatus and An. minimus were predominant during both the dry and rainy seasons, but showed highly zoophilic preferences (Zoophilic index of 98% and 95%, respectively). Overall, 22% of the total malaria vectors were collected between 10:00 PM and 5:00 AM indoors when people are sleeping. Twenty-seven percent of primary and secondary vectors were collected outdoors before 10:00 PM or after 5:00 AM, times when people are usually awake and outdoors. Only two specimens were positive for P. falciparum, one An. aconitus from Phongsaly and one An. minimus from Vientiane Province CONCLUSIONS: The results indicate that people living in rural areas in Laos are constantly exposed to malaria vectors throughout the year and specifically outdoors. The use of LLINs/IRS remains important but innovative tools and new strategies are needed to address locally, the early and outdoor malaria transmission. Lack of expertise in general entomological methods may further exacerbate the situation.