Relationships between biological age, distance from aquatic habitats and pyrethroid resistance status of Anopheles funestus mosquitoes in south-eastern Tanzania

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

Genetic markers associated with the widespread insecticide resistance in malaria vector Anopheles funestus populations across Tanzania

BACKGROUND

Anopheles funestus is a leading vector of malaria in most parts of East and Southern Africa, yet its ecology and responses to vector control remain poorly understood compared with other vectors such as Anopheles gambiae and Anopheles arabiensis. This study presents the first large-scale survey of the genetic and phenotypic expression of insecticide resistance in An. funestus populations in Tanzania.

METHODS

We performed insecticide susceptibility bioassays on An. funestus mosquitoes in nine regions with moderate-to-high malaria prevalence in Tanzania, followed by genotyping for resistance-associated mutations (CYP6P9a, CYP6P9b, L119F-GSTe2) and structural variants (SV4.3 kb, SV6.5 kb). Generalized linear models were used to assess relationships between genetic markers and phenotypic resistance. An interactive R Shiny tool was created to visualize the data and support evidence-based interventions.

RESULTS

Pyrethroid resistance was universal but reversible by piperonyl-butoxide (PBO). However, carbamate resistance was observed in only five of the nine districts, and dichloro-diphenyl-trichloroethane (DDT) resistance was found only in the Kilombero valley, south-eastern Tanzania. Conversely, there was universal susceptibility to the organophosphate pirimiphos-methyl in all sites. Genetic markers of resistance had distinct geographical patterns, with CYP6P9a-R and CYP6P9b-R alleles, and the SV6.5 kb structural variant absent or undetectable in the north-west but prevalent in all other sites, while SV4.3 kb was prevalent in the north-western and western regions but absent elsewhere. Emergent L119F-GSTe2, associated with deltamethrin resistance, was detected in heterozygous form in districts bordering Mozambique, Malawi and the Democratic Republic of Congo. The resistance landscape was most complex in western Tanzania, in Tanganyika district, where all five genetic markers were detected. There was a notable south-to-north spread of resistance genes, especially CYP6P9a-R, though this appears to be interrupted, possibly by the Rift Valley.

CONCLUSIONS

This study underscores the need to expand resistance monitoring to include An. funestus alongside other vector species, and to screen for both the genetic and phenotypic signatures of resistance. The findings can be visualized online via an interactive user interface and could inform data-driven decision-making for resistance management and vector control. Since this was the first large-scale survey of resistance in Tanzania's An. funestus, we recommend regular updates with greater geographical and temporal coverage.

Fine-scale spatial distribution of deltamethrin resistance and population structure of Anopheles funestus and Anopheles arabiensis populations in Southern Mozambique

BACKGROUND

Insecticide resistance in malaria vectors can be spatially highly heterogeneous, yet population structure analyses frequently find relatively high levels of gene flow among mosquito populations. Few studies have contemporaneously assessed phenotypic, genotypic and population structure analysis on mosquito populations and none at fine geographical scales. In this study, genetic diversity, population structure, and insecticide resistance profiles of Anopheles funestus and Anopheles arabiensis were examined across mosquito populations from and within neighbouring villages.

METHODS

Mosquitoes were collected from 11 towns in southern Mozambique, as well as from different neighbourhoods within the town of Palmeira, during the peak malaria transmission season in 2016. CDC bottle bioassay and PCR assays were performed with Anopheles mosquitoes at each site to determine phenotypic and molecular insecticide resistance profiles, respectively. Microsatellite analysis was conducted on a subsample of mosquitoes to estimate genetic diversity and population structure.

RESULTS

Phenotypic insecticide resistance to deltamethrin was observed in An. funestus sensu stricto (s.s.) throughout the area, though a high level of mortality variation was seen. However, 98% of An. funestus s.s. were CYP6P9a homozygous resistant. An. arabiensis was phenotypically susceptible to deltamethrin and 99% were kdr homozygous susceptible. Both Anopheles species exhibited high allelic richness and heterozygosity. Significant deviations from Hardy-Weinberg equilibrium were observed, and high linkage disequilibrium was seen for An. funestus s.s., supporting population subdivision. However, the FST values were low for both anophelines (- 0.00457 to 0.04213), Nm values were high (9.4-71.8 migrants per generation), AMOVA results showed almost 100% genetic variation among and within individuals, and Structure analysis showed no clustering of An. funestus s.s. and An. arabiensis populations. These results suggest high gene flow among mosquito populations.

CONCLUSION

Despite a relatively high level of phenotypic variation in the An. funestus population, molecular analysis shows the population is admixed. These data indicate that CYP6P9a resistance markers do not capture all phenotypic variation in the area, but also that resistance genes of high impact are likely to easily spread in the area. Conversely, other strategies, such as transgenic mosquito release programmes will likely not face challenges in this locality.