Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates.

BACKGROUND: Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl). RESULTS: Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups. CONCLUSION: Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Whole transcriptomic analysis reveals overexpression of salivary gland and cuticular proteins genes in insecticide-resistant Anopheles arabiensis from Western Kenya.

BACKGROUND: Effective vector control is key to malaria prevention. However, this is now compromised by increased insecticide resistance due to continued reliance on insecticide-based control interventions. In Kenya, we have observed heterogenous resistance to pyrethroids and organophosphates in Anopheles arabiensis which is one of the most widespread malaria vectors in the country. We investigated the gene expression profiles of insecticide resistant An. arabiensis populations from Migori and Siaya counties in Western Kenya using RNA-Sequencing. Centers for Disease Control and Prevention (CDC) bottle assays were conducted using deltamethrin (DELTA), alphacypermethrin (ACYP) and pirimiphos-methyl (PMM) to determine the resistance status in both sites. RESULTS: Mosquitoes from Migori had average mortalities of 91%, 92% and 58% while those from Siaya had 85%, 86%, and 30% when exposed to DELTA, ACYP and PMM, respectively. RNA-Seq analysis was done on pools of mosquitoes which survived exposure ('resistant'), mosquitoes that were not exposed, and the insecticide-susceptible An. arabiensis Dongola strain. Gene expression profiles of resistant mosquitoes from both Migori and Siaya showed an overexpression mainly of salivary gland proteins belonging to both the short and long form D7 genes, and cuticular proteins (including CPR9, CPR10, CPR15, CPR16). Additionally, the overexpression of detoxification genes including cytochrome P450s (CYP9M1, CYP325H1, CYP4C27, CYP9L1 and CYP307A1), 2 carboxylesterases and a glutathione-S-transferase (GSTE4) were also shared between DELTA, ACYP, and PMM survivors, pointing to potential contribution to cross resistance to both pyrethroid and organophosphate insecticides. CONCLUSION: This study provides novel insights into the molecular basis of insecticide resistance in An. arabiensis in Western Kenya and suggests that salivary gland proteins and cuticular proteins are associated with resistance to multiple classes of insecticides.

Key gene modules and hub genes associated with pyrethroid and organophosphate resistance in Anopheles mosquitoes: a systems biology approach.

Indoor residual spraying (IRS) and insecticide-treated nets (ITNs) are the main methods used to control mosquito populations for malaria prevention. The efficacy of these strategies is threatened by the spread of insecticide resistance (IR), limiting the success of malaria control. Studies of the genetic evolution leading to insecticide resistance could enable the identification of molecular markers that can be used for IR surveillance and an improved understanding of the molecular mechanisms associated with IR. This study used a weighted gene co-expression network analysis (WGCNA) algorithm, a systems biology approach, to identify genes with similar co-expression patterns (modules) and hub genes that are potential molecular markers for insecticide resistance surveillance in Kenya and Benin. A total of 20 and 26 gene co-expression modules were identified via average linkage hierarchical clustering from Anopheles arabiensis and An. gambiae, respectively, and hub genes (highly connected genes) were identified within each module. Three specific genes stood out: serine protease, E3 ubiquitin-protein ligase, and cuticular proteins, which were top hub genes in both species and could serve as potential markers and targets for monitoring IR in these malaria vectors. In addition to the identified markers, we explored molecular mechanisms using enrichment maps that revealed a complex process involving multiple steps, from odorant binding and neuronal signaling to cellular responses, immune modulation, cellular metabolism, and gene regulation. Incorporation of these dynamics into the development of new insecticides and the tracking of insecticide resistance could improve the sustainable and cost-effective deployment of interventions.

The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa.

BACKGROUND: Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. The current study sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well as describe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzii populations. METHODS: Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing (WGS) data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Comparisons were made with An. coluzzii cohorts from West and Central Africa. RESULTS: This study reports the detection of An. coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points from which samples were analyzed and its presence confirmed through taxonomic analysis. Additionally, there was a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies up to 64%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya. CONCLUSIONS: These findings emphasize the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

Detection of Anopheles stephensi Mosquitoes by Molecular Surveillance, Kenya.

The Anopheles stephensi mosquito is an invasive malaria vector recently reported in Djibouti, Ethiopia, Sudan, Somalia, Nigeria, and Ghana. The World Health Organization has called on countries in Africa to increase surveillance efforts to detect and report this vector and institute appropriate and effective control mechanisms. In Kenya, the Division of National Malaria Program conducted entomological surveillance in counties at risk for An. stephensi mosquito invasion. In addition, the Kenya Medical Research Institute conducted molecular surveillance of all sampled Anopheles mosquitoes from other studies to identify An. stephensi mosquitoes. We report the detection and confirmation of An. stephensi mosquitoes in Marsabit and Turkana Counties by using endpoint PCR and morphological and sequence identification. We demonstrate the urgent need for intensified entomological surveillance in all areas at risk for An. stephensi mosquito invasion, to clarify its occurrence and distribution and develop tailored approaches to prevent further spread.

Screening eaves of houses reduces indoor mosquito density in rural, western Kenya.

BACKGROUND: Despite the scale-up of insecticide-treated nets and indoor residual spraying, the bulk of malaria transmission in western Kenya still occurs indoors, late at night. House improvement is a potential long-term solution to further reduce malaria transmission in the region. METHODS: The impact of eave screening on mosquito densities was evaluated in two rural villages in western Kenya. One-hundred-and-twenty pairs of structurally similar, neighbouring houses were used in the study. In each pair, one house was randomly selected to receive eave screening at the beginning of the study while the other remained unscreened until the end of the sampling period. Mosquito sampling was performed monthly by motorized aspiration method for 4 months. The collected mosquitoes were analysed for species identification. RESULTS: Compared to unscreened houses, significantly fewer female Anopheles funestus (RR = 0.40, 95% CI 0.29-0.55), Anopheles gambiae Complex (RR = 0.46, 95% CI 0.34-0.62) and Culex species (RR = 0.53, 95% CI 0.45-0.61) were collected in screened houses. No significant differences in the densities of the mosquitoes were detected in outdoor collections. Significantly fewer Anopheles funestus were collected indoors from houses with painted walls (RR = 0.05, 95% CI 0.01-0.38) while cooking in the house was associated with significantly lower numbers of Anopheles gambiae Complex indoors (RR = 0.60, 95% CI 0.45-0.79). Nearly all house owners (99.6%) wanted their houses permanently screened, including 97.7% that indicated a willingness to use their own resources. However, 99.2% required training on house screening. The cost of screening a single house was estimated at KES6,162.38 (US$61.62). CONCLUSION: Simple house modification by eave screening has the potential to reduce the indoor occurrence of both Anopheles and Culex mosquito species. Community acceptance was very high although education and mobilization may be needed for community uptake of house modification for vector control. Intersectoral collaboration and favourable government policies on housing are important links towards the adoption of house improvements for malaria control.

Western Kenyan Anopheles gambiae showing intense permethrin resistance harbour distinct microbiota.

BACKGROUND: Insecticide resistance poses a growing challenge to malaria vector control in Kenya and around the world. Following evidence of associations between the mosquito microbiota and insecticide resistance, the microbiota of Anopheles gambiae sensu stricto (s.s.) from Tulukuyi village, Bungoma, Kenya, with differing permethrin resistance profiles were comparatively characterized. METHODS: Using the CDC bottle bioassay, 133 2-3 day-old, virgin, non-blood fed female F1 progeny of field-caught An. gambiae s.s. were exposed to five times (107.5 µg/ml) the discriminating dose of permethrin. Post bioassay, 50 resistant and 50 susceptible mosquitoes were subsequently screened for kdr East and West mutations, and individually processed for microbial analysis using high throughput sequencing targeting the universal bacterial and archaeal 16S rRNA gene. RESULTS: 47 % of the samples tested (n = 133) were resistant, and of the 100 selected for further processing, 99 % were positive for kdr East and 1 % for kdr West. Overall, 84 bacterial taxa were detected across all mosquito samples, with 36 of these shared between resistant and susceptible mosquitoes. A total of 20 bacterial taxa were unique to the resistant mosquitoes and 28 were unique to the susceptible mosquitoes. There were significant differences in bacterial composition between resistant and susceptible individuals (PERMANOVA, pseudo-F = 2.33, P = 0.001), with presence of Sphingobacterium, Lysinibacillus and Streptococcus (all known pyrethroid-degrading taxa), and the radiotolerant Rubrobacter, being significantly associated with resistant mosquitoes. On the other hand, the presence of Myxococcus, was significantly associated with susceptible mosquitoes. CONCLUSIONS: This is the first report of distinct microbiota in An. gambiae s.s. associated with intense pyrethroid resistance. The findings highlight differentially abundant bacterial taxa between resistant and susceptible mosquitoes, and further suggest a microbe-mediated mechanism of insecticide resistance in mosquitoes. These results also indicate fixation of the kdr East mutation in this mosquito population, precluding further analysis of its associations with the mosquito microbiota, but presenting the hypothesis that any microbe-mediated mechanism of insecticide resistance would be likely of a metabolic nature. Overall, this study lays initial groundwork for understanding microbe-mediated mechanisms of insecticide resistance in African mosquito vectors of malaria, and potentially identifying novel microbial markers of insecticide resistance that could supplement existing vector surveillance tools.

The Anopheles coluzzii range extends into Kenya: Detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa.

Background: Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. In the current study, we sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well asdescribe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzi populations. Methods: Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Results: This study reports the detection of Anopheles coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points sampled and its presence confirmed through taxonomic analysis. Additionally, we found a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies of ~60%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya. Conclusions: These findings emphasise the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

Microsporidia MB in the primary malaria vector Anopheles gambiae sensu stricto is avirulent and undergoes maternal and horizontal transmission.

BACKGROUND: The demonstration that the recently discovered Anopheles symbiont Microsporidia MB blocks malaria transmission in Anopheles arabiensis and undergoes vertical and horizontal transmission suggests that it is a promising candidate for the development of a symbiont-based malaria transmission-blocking strategy. The infection prevalence and characteristics of Microsporidia MB in Anopheles gambiae sensu stricto (s.s.), another primary vector species of malaria in Kenya, were investigated. METHODS: Field-collected females were confirmed to be Microsporidia MB-positive after oviposition. Egg counts of Microsporidia MB-infected and non-infected individuals were used to infer the effects of Microsporidia MB on fecundity. The time to pupation, adult sex ratio and survival were used to determine if Microsporidia MB infection has similar characteristics in the host mosquitoes An. gambiae and An. arabiensis. The intensity of Microsporidia MB infection in tissues of the midgut and gonads, and in carcasses, was determined by quantitative polymerase chain reaction. To investigate horizontal transmission, virgin males and females that were either Microsporidia MB-infected or non-infected were placed in standard cages for 48 h and allowed to mate; transmission was confirmed by quantitative polymerase chain reaction targeting Microsporidia MB genes. RESULTS: Microsporidia MB was found to naturally occur at a low prevalence in An. gambiae s.s. collected in western Kenya. Microsporidia MB shortened the development time from larva to pupa, but other fitness parameters such as fecundity, sex ratio, and adult survival did not differ between Microsporidia MB-infected and non-infected hosts. Microsporidia MB intensities were high in the male gonadal tissues. Transmission experiments indicated that Microsporidia MB undergoes both maternal and horizontal transmission in An. gambiae s.s. CONCLUSIONS: The findings that Microsporidia MB naturally infects, undergoes maternal and horizontal transmission, and is avirulent in An. gambiae s.s. indicate that many of the characteristics of its infection in An. arabiensis hold true for the former. The results of the present study indicate that Microsporidia MB could be developed as a tool for the transmission-blocking of malaria across different Anopheles species.

Comparison of four outdoor mosquito trapping methods as potential replacements for human landing catches in western Kenya.

INTRODUCTION: Longitudinal monitoring of outdoor-biting malaria vector populations is becoming increasingly important in understanding the dynamics of residual malaria transmission. However, the human landing catch (HLC), the gold standard for measuring human biting rates indoors and outdoors, is costly and raises ethical concerns related to increased risk of infectious bites among collectors. Consequently, routine data on outdoor-feeding mosquito populations are usually limited because of the lack of a scalable tool with similar sensitivity to outdoor HLC. METHODOLOGY: The Anopheles trapping sensitivity of four baited proxy outdoor trapping methods-Furvela tent trap (FTT), host decoy trap (HDT), mosquito electrocuting traps (MET) and outdoor CDC light traps (OLT)-was assessed relative to HLC in a 5 × 5 replicated Latin square conducted over 25 nights in two villages of western Kenya. Indoor CDC light trap (ILT) was run in one house in each of the compounds with outdoor traps, while additional non-Latin square indoor and outdoor HLC collections were performed in one of the study villages. RESULTS: The MET, FTT, HDT and OLT sampled approximately 4.67, 7.58, 5.69 and 1.98 times more An. arabiensis compared to HLC, respectively, in Kakola Ombaka. Only FTT was more sensitive relative to HLC in sampling An. funestus in Kakola Ombaka (RR = 5.59, 95% CI 2.49-12.55, P < 0.001) and Masogo (RR = 4.38, 95% CI 1.62-11.80, P = 0.004) and in sampling An. arabiensis in Masogo (RR = 5.37, 95% CI 2.17-13.24, P < 0.001). OLT sampled significantly higher numbers of An. coustani in Kakola Ombaka (RR = 3.03, 95% CI 1.65-5.56, P < 0.001) and Masogo (RR = 2.88, 95% CI 1.15-7.22, P = 0.02) compared to HLC. OLT, HLC and MET sampled mostly An. coustani, FTT had similar proportions of An. funestus and An. arabiensis, while HDT sampled predominantly An. arabiensis in both villages. FTT showed close correlation with ILT in vector abundance for all three species at both collection sites. CONCLUSION: FTT and OLT are simple, easily scalable traps and are potential replacements for HLC in outdoor sampling of Anopheles mosquitoes. However, the FTT closely mirrored indoor CDC light trap in mosquito indices and therefore may be more of an indoor mimic than a true outdoor collection tool. HDT and MET show potential for sampling outdoor host-seeking mosquitoes. However, the traps as currently designed may not be feasible for large-scale, longitudinal entomological monitoring. Therefore, the baited outdoor CDC light trap may be the most appropriate tool currently available for assessment of outdoor-biting and malaria transmission risk.

Impact of indoor residual spraying with pirimiphos-methyl (Actellic 300CS) on entomological indicators of transmission and malaria case burden in Migori County, western Kenya.

Indoor residual spraying (IRS) of insecticides is a major vector control strategy for malaria prevention. We evaluated the impact of a single round of IRS with the organophosphate, pirimiphos-methyl (Actellic 300CS), on entomological and parasitological parameters of malaria in Migori County, western Kenya in 2017, in an area where primary vectors are resistant to pyrethroids but susceptible to the IRS compound. Entomological monitoring was conducted by indoor CDC light trap, pyrethrum spray catches (PSC) and human landing collection (HLC) before and after IRS. The residual effect of the insecticide was assessed monthly by exposing susceptible An. gambiae s.s. Kisumu strain to sprayed surfaces in cone assays and measuring mortality at 24 hours. Malaria case burden data were extracted from laboratory records of four health facilities within the sprayed area and two adjacent unsprayed areas. IRS was associated with reductions in An. funestus numbers in the intervention areas compared to non-intervention areas by 88% with light traps (risk ratio [RR] 0.12, 95% CI 0.07-0.21, p < 0.001) and 93% with PSC collections (RR = 0.07, 0.03-0.17, p < 0.001). The corresponding reductions in the numbers of An. arabiensis collected by PSC were 69% in the intervention compared to the non-intervention areas (RR = 0.31, 0.14-0.68, p = 0.006), but there was no significant difference with light traps (RR = 0.45, 0.21-0.96, p = 0.05). Before IRS, An. funestus accounted for over 80% of Anopheles mosquitoes collected by light trap and PSC in all sites. After IRS, An. arabiensis accounted for 86% of Anopheles collected by PSC and 66% by CDC light trap in the sprayed sites while the proportion in non-intervention sites remained unchanged. No sporozoite infections were detected in intervention areas after IRS and biting rates by An. funestus were reduced to near zero. Anopheles funestus and An. arabiensis were fully susceptible to pirimiphos-methyl and resistant to pyrethroids. The residual effect of Actellic 300CS lasted ten months on mud and concrete walls. Malaria case counts among febrile patients within IRS areas was lower post- compared to pre-IRS by 44%, 65% and 47% in Rongo, Uriri and Nyatike health facilities respectively. A single application of IRS with Actellic 300CS in Migori County provided ten months protection and resulted in the near elimination of the primary malaria vector An. funestus and a corresponding reduction of malaria case count among out-patients. The impact was less on An. arabiensis, most likely due to their exophilic nature.