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.

Wild populations of malaria vectors can mate both inside and outside human dwellings

BACKGROUND

Wild populations of Anopheles mosquitoes are generally thought to mate outdoors in swarms, although once colonized, they also mate readily inside laboratory cages. This study investigated whether the malaria vectors Anopheles funestus and Anopheles arabiensis can also naturally mate inside human dwellings.

METHOD

Mosquitoes were sampled from three volunteer-occupied experimental huts in a rural Tanzanian village at 6:00 p.m. each evening, after which the huts were completely sealed and sampling was repeated at 11:00 p.m and 6 a.m. the next morning to compare the proportions of inseminated females. Similarly timed collections were done inside local unsealed village houses. Lastly, wild-caught larvae and pupae were introduced inside or outside experimental huts constructed inside two semi-field screened chambers. The huts were then sealed and fitted with exit traps, allowing mosquito egress but not entry. Mating was assessed in subsequent days by sampling and dissecting emergent adults caught indoors, outdoors and in exit traps.

RESULTS

Proportions of inseminated females inside the experimental huts in the village increased from approximately  60% at 6 p.m. to approximately 90% the following morning despite no new mosquitoes entering the huts after 6 p.m. Insemination in the local homes increased from approximately 78% to approximately 93% over the same time points. In the semi-field observations of wild-caught captive mosquitoes, the proportions of inseminated An. funestus were 20.9% (95% confidence interval [CI]: ± 2.8) outdoors, 25.2% (95% CI: ± 3.4) indoors and 16.8% (± 8.3) in exit traps, while the proportions of inseminated An. arabiensis were 42.3% (95% CI: ± 5.5) outdoors, 47.4% (95% CI: ± 4.7) indoors and 37.1% (CI: ± 6.8) in exit traps.

CONCLUSION

Wild populations of An. funestus and An. arabiensis in these study villages can mate both inside and outside human dwellings. Most of the mating clearly happens before the mosquitoes enter houses, but additional mating happens indoors. The ecological significance of such indoor mating remains to be determined. The observed insemination inside the experimental huts fitted with exit traps and in the unsealed village houses suggests that the indoor mating happens voluntarily even under unrestricted egress. These findings may inspire improved vector control, such as by targeting males indoors, and potentially inform alternative methods for colonizing strongly eurygamic Anopheles species (e.g. An. funestus) inside laboratories or semi-field chambers.

Using a miniaturized double-net trap (DN-Mini) to assess relationships between indoor-outdoor biting preferences and physiological ages of two malaria vectors, Anopheles arabiensis and Anopheles funestus

BACKGROUND

Effective malaria surveillance requires detailed assessments of mosquitoes biting indoors, where interventions such as insecticide-treated nets work best, and outdoors, where other interventions may be required. Such assessments often involve volunteers exposing their legs to attract mosquitoes [i.e., human landing catches (HLC)], a procedure with significant safety and ethical concerns. Here, an exposure-free, miniaturized, double-net trap (DN-Mini) is used to assess relationships between indoor-outdoor biting preferences of malaria vectors, Anopheles arabiensis and Anopheles funestus, and their physiological ages (approximated by parity and insemination states).

METHODS

The DN-Mini is made of UV-resistant netting on a wooden frame and PVC base. At 100 cm × 60 cm × 180 cm, it fits indoors and outdoors. It has a protective inner chamber where a volunteer sits and collects host-seeking mosquitoes entrapped in an outer chamber. Experiments were conducted in eight Tanzanian villages using DN-Mini to: (a) estimate nightly biting and hourly biting proportions of mosquitoes indoors and outdoors; (b) compare these proportions to previous estimates by HLC in same villages; and, (c) compare distribution of parous (proxy for potentially infectious) and inseminated mosquitoes indoors and outdoors.

RESULTS

More than twice as many An. arabiensis were caught outdoors as indoors (p < 0.001), while An. funestus catches were marginally higher indoors than outdoors (p = 0.201). Anopheles arabiensis caught outdoors also had higher parity and insemination proportions than those indoors (p < 0.001), while An. funestus indoors had higher parity and insemination than those outdoors (p = 0.04). Observations of indoor-biting and outdoor-biting proportions, hourly biting patterns and overall species diversities as measured by DN-Mini, matched previous HLC estimates.

CONCLUSIONS

Malaria vectors that are behaviourally adapted to bite humans outdoors also have their older, potentially infectious sub-populations concentrated outdoors, while those adapted to bite indoors have their older sub-populations concentrated indoors. Here, potentially infectious An. arabiensis more likely bite outdoors than indoors, while potentially infectious An. funestus more likely bite indoors. These observations validate previous evidence that even outdoor-biting mosquitoes regularly enter houses when young. They also demonstrate efficacy of DN-Mini for measuring indoor-outdoor biting behaviours of mosquitoes, their hourly biting patterns and epidemiologically relevant parameters, e.g., parity and insemination status, without exposure to volunteers. The trap is easy-to-use, easy-to-manufacture and affordable (prototypes cost ~ 100 US$/unit).

Swarms of the malaria vector Anopheles funestus in Tanzania

Background

Anopheles funestus mosquitoes currently contribute more than 85% of ongoing malaria transmission events in south-eastern Tanzania, even though they occur in lower densities than other vectors, such as Anopheles arabiensis. Unfortunately, the species ecology is minimally understood, partly because of difficulties in laboratory colonization. This study describes the first observations of An. funestus swarms in Tanzania, possibly heralding new opportunities for control.

Method

Using systematic searches by community-based volunteers and expert entomologists, An. funestus swarms were identified in two villages in Ulanga and Kilombero districts in south-eastern Tanzania, starting June 2018. Swarms were characterized by size, height, start- and end-times, presence of copulation and associated environmental features. Samples of male mosquitoes from the swarms were examined for sexual maturity by observing genitalia rotation, species identity using polymerase chain reaction and wing sizes.

Results

581 An. funestus (98.1% males (n = 570) and 1.9% (n = 11) females) and 9 Anopheles gambiae sensu lato (s.l.) males were sampled using sweep nets from the 81 confirmed swarms in two villages (Ikwambi in Kilombero district and Tulizamoyo in Ulanga district). Six copulation events were observed in the swarms. Mean density (95% CL) of An. funestus caught/swarm/village/evening was 6.6 (5.9–7.2) in Tulizamoyo and 10.8 (5.8–15.8) in Ikwambi. 87.7% (n = 71) of the swarms were found in Tulizamoyo, while 12.3% (n = 10) were in Ikwambi. Mean height of swarms was 1.7 m (0.9–2.5 m), while mean duration was 12.9 (7.9–17.9) minutes. The PCR analysis confirmed that 100% of all An. funestus s.l. samples processed were An. funestus sensu stricto. Mean wing length of An. funestus males was 2.47 mm (2.0–2.8 mm), but there was no difference between swarming males and indoor-resting males. Most swarms (95.0%) occurred above bare ground, sometime on front lawns near human dwellings, and repeatedly in the same locations.

Conclusion

This study has demonstrated occurrence of An. funestus swarms for the first time in Tanzania. Further investigations could identify new opportunities for improved control of this dominant malaria vector, possibly by targeting the swarms.

Housing gaps, mosquitoes and public viewpoints: a mixed methods assessment of relationships between house characteristics, malaria vector biting risk and community perspectives in rural Tanzania

BACKGROUND

House improvement and environmental management can significantly improve malaria transmission control in endemic communities. This study assessed the influence of physical characteristics of houses and surrounding environments on mosquito biting risk in rural Tanzanian villages, and examined knowledge and perceptions of residents on relationships between these factors and malaria transmission. The study further assessed whether people worried about these risks and how they coped.

METHODS

Entomological surveys of indoor mosquito densities were conducted across four villages in Ulanga district, south-eastern Tanzania. The survey involved 48 sentinel houses sampled monthly and other sets of 48 houses randomly recruited each month for one-off sampling over 12 months. Physical characteristics of the houses and surrounding environments were recorded. Questionnaire surveys were administered to 200 household heads to assess their knowledge and concerns regarding the observed housing and environmental features, and whether they considered these features when constructing houses. Focus group discussions, were conducted to clarify emergent themes on people's perceptions on relationships between housing or environmental factors and malaria transmission.

RESULTS

The entomological surveys showed statistically higher indoor densities of the malaria vectors (Anopheles arabiensis and Anopheles funestus) in houses with mud walls compared to plastered or brick walls, open eaves compared to closed eaves and unscreened windows compared to screened windows. Most respondents reported that their houses allowed mosquito entry, at least partially. Participants were aware that house structure and environmental characteristics influenced indoor mosquito densities and consequently malaria transmission. They were concerned about living in poorly-constructed houses with gaps on eaves, walls, windows and doors but were constrained by low income.

CONCLUSION

In rural south-eastern Tanzania, significant proportions of people still live in houses with open eaves, unscreened windows and gaps on doors. Though they are fully aware of associated mosquito biting and pathogen transmission risks, they are constrained by low-income levels. The study proposes that community-based house improvement initiatives combined with targeted subsidies could lower the financial barriers, improve access to essential construction materials or designs, and significantly accelerate malaria transmission control in these communities.

New evidence of mating swarms of the malaria vector, Anopheles arabiensis in Tanzania

Background: Malaria mosquitoes form mating swarms around sunset, often at the same locations for months or years. Unfortunately, studies of Anopheles swarms are rare in East Africa, the last recorded field observations in Tanzania having been in 1983. Methods: Mosquito swarms were surveyed by trained volunteers between August-2016 and June-2017 in Ulanga district, Tanzania. Identified Anopheles swarms were sampled using sweep nets, and collected mosquitoes killed by refrigeration then identified by sex and taxa. Sub-samples were further identified by PCR, and spermatheca of females examined for mating status. Mosquito ages were estimated by observing female ovarian tracheoles and rotation of male genitalia. GPS locations, types of swarm markers, start/end times of swarming, heights above ground, mosquito counts/swarm, and copulation events were recorded. Results: A total of 216 Anopheles swarms were identified, characterized and mapped, from which 7,142 Anopheles gambiae s.l and 13 Anopheles funestus were sampled. The An. gambiae s.l were 99.6% males and 0.4% females, while the An. funestus were all males. Of all An. gambiae s.l analyzed by PCR, 86.7% were An. arabiensis, while 13.3% returned non-amplified DNA. Mean height (±SD) of swarms was 2.74±0.64m, and median duration was 20 (IQR; 15-25) minutes. Confirmed swarm markers included rice fields (25.5%), burned grounds (17.2%), banana trees (13%), brick piles (8.8%), garbage heaps (7.9%) and ant-hills (7.4%). Visual estimates of swarm sizes by the volunteers was strongly correlated to actual sizes by sweep nets (R=0.94; P=<0.001). All females examined were nulliparous and 95.6% [N=6787] of males had rotated genitalia, indicating sexual maturity. Conclusions: This is the first report of Anopheles swarms in Tanzania in more than three decades. The study demonstrates that the swarms can be identified and characterized by trained community-based volunteers, and highlights potential new interventions, for example targeted aerosol spraying of the swarms to improve malaria control.

Interventions that effectively target Anopheles funestus mosquitoes could significantly improve control of persistent malaria transmission in south-eastern Tanzania

Malaria is transmitted by many Anopheles species whose proportionate contributions vary across settings. We re-assessed the roles of Anopheles arabiensis and Anopheles funestus, and examined potential benefits of species-specific interventions in an area in south-eastern Tanzania, where malaria transmission persists, four years after mass distribution of long-lasting insecticide-treated nets (LLINs). Monthly mosquito sampling was done in randomly selected households in three villages using CDC light traps and back-pack aspirators, between January-2015 and January-2016, four years after the last mass distribution of LLINs in 2011. Multiplex polymerase chain reaction (PCR) was used to identify members of An. funestus and Anopheles gambiae complexes. Enzyme-linked immunosorbent assay (ELISA) was used to detect Plasmodium sporozoites in mosquito salivary glands, and to identify sources of mosquito blood meals. WHO susceptibility assays were done on wild caught female An. funestus s.l, and physiological ages approximated by examining mosquito ovaries for parity. A total of 20,135 An. arabiensis and 4,759 An. funestus were collected. The An. funestus group consisted of 76.6% An. funestus s.s, 2.9% An. rivulorum, 7.1% An. leesoni, and 13.4% unamplified samples. Of all mosquitoes positive for Plasmodium, 82.6% were An. funestus s.s, 14.0% were An. arabiensis and 3.4% were An. rivulorum. An. funestus and An. arabiensis contributed 86.21% and 13.79% respectively, of annual entomological inoculation rate (EIR). An. arabiensis fed on humans (73.4%), cattle (22.0%), dogs (3.1%) and chicken (1.5%), but An. funestus fed exclusively on humans. The An. funestus populations were 100% susceptible to organophosphates, pirimiphos methyl and malathion, but resistant to permethrin (10.5% mortality), deltamethrin (18.7%), lambda-cyhalothrin (18.7%) and DDT (26.2%), and had reduced susceptibility to bendiocarb (95%) and propoxur (90.1%). Parity rate was higher in An. funestus (65.8%) than An. arabiensis (44.1%). Though An. arabiensis is still the most abundant vector species here, the remaining malaria transmission is predominantly mediated by An. funestus, possibly due to high insecticide resistance and high survival probabilities. Interventions that effectively target An. funestus mosquitoes could therefore significantly improve control of persistent malaria transmission in south–eastern Tanzania.

Rapid and Non-destructive Detection and Identification of Two Strains of Wolbachia in Aedes aegypti by Near-Infrared Spectroscopy

The release of Wolbachia infected mosquitoes is likely to form a key component of disease control strategies in the near future. We investigated the potential of using near-infrared spectroscopy (NIRS) to simultaneously detect and identify two strains of Wolbachia pipientis (wMelPop and wMel) in male and female laboratory-reared Aedes aegypti mosquitoes. Our aim is to find faster, cheaper alternatives for monitoring those releases than the molecular diagnostic techniques that are currently in use. Our findings indicate that NIRS can differentiate females and males infected with wMelPop from uninfected wild type samples with an accuracy of 96% (N = 299) and 87.5% (N = 377), respectively. Similarly, females and males infected with wMel were differentiated from uninfected wild type samples with accuracies of 92% (N = 352) and 89% (N = 444). NIRS could differentiate wMelPop and wMel transinfected females with an accuracy of 96.6% (N = 442) and males with an accuracy of 84.5% (N = 443). This non-destructive technique is faster than the standard polymerase chain reaction diagnostic techniques. After the purchase of a NIRS spectrometer, the technique requires little sample processing and does not consume any reagents.

Near infrared spectroscopy (NIRS) is a technique that measures specific frequencies of light absorbed by C-H, O-H, S-H and N-H functional groups. Mosquito samples are grouped based upon absorption differences between their chemical properties. In this study, we used NIRS to differentiate 1) Aedes aegypti infected with either of the two strains of intracellular bacterium Wolbachia (wMel and wMelPop) from wild type Ae. aegypti and 2) Aedes aegypti infected with wMel from those infected with wMelPoP. NIRS facilitated the differentiation of wMel and wMelPop from wild type samples and samples infected with either of the Wolbachia infected strains with high prediction accuracies over their lifespan. Predictive models were derived from initial calibration data sets and validated against independent cohorts. Prediction accuracies were high (82–98%) regardless of the cohort mosquitoes were sampled from. The results show that NIRS may have real potential as an alternative method for monitoring Wolbachia incidence in mosquitoes. A rapid, simple and cost-effective surveillance tool suitable for resource-poor areas and large urban release programs would be of great utility for evaluating Wolbachia-based interventions. The models developed during this study require further validation using field collections.

Correlations between household occupancy and malaria vector biting risk in rural Tanzanian villages: implications for high-resolution spatial targeting of control interventions

Background

Fine-scale targeting of interventions is increasingly important where epidemiological disease profiles depict high geographical stratifications. This study verified correlations between household biomass and mosquito house-entry using experimental hut studies, and then demonstrated how geographical foci of mosquito biting risk can be readily identified based on spatial distributions of household occupancies in villages.

Methods

A controlled 4 × 4 Latin square experiment was conducted in rural Tanzania, in which no, one, three or six adult male volunteers slept under intact bed nets, in experimental huts. Mosquitoes entering the huts were caught using exit interception traps on eaves and windows. Separately, monthly mosquito collections were conducted in 96 randomly selected households in three villages using CDC light traps between March-2012 and November-2013. The number of people sleeping in the houses and other household and environmental characteristics were recorded. ArcGIS 10 (ESRI-USA) spatial analyst tool, Gi* Ord Statistic was used to analyse clustering of vector densities and household occupancy.

Results

The densities of all mosquito genera increased in huts with one, three or six volunteers, relative to huts with no volunteers, and direct linear correlations within tested ranges (P < 0.001). Significant geographical clustering of indoor densities of malaria vectors, Anopheles arabiensis and Anopheles funestus, but not Culex or Mansonia species occurred in locations where households with highest occupancy were also most clustered (Gi* P ≤ 0.05, and Gi* Z-score ≥1.96).

Conclusions

This study demonstrates strong correlations between household occupancy and malaria vector densities in households, but also spatial correlations of these variables within and between villages in rural southeastern Tanzania. Fine-scale clustering of indoor densities of vectors within and between villages occurs in locations where houses with highest occupancy are also clustered. The study indicates potential for using household census data to preliminarily identify households with greatest Anopheles mosquito biting risk.