A novel tetra-primer ARMS-PCR approach for the molecular karyotyping of chromosomal inversion 2Ru in the main malaria vectors Anopheles gambiae and Anopheles coluzzii.

BACKGROUND: Chromosomal inversion polymorphisms have been associated with adaptive behavioral, physiological, morphological and life history traits in the two main Afrotropical malaria vectors, Anopheles coluzzii and Anopheles gambiae. The understanding of the adaptive value of chromosomal inversion systems is constrained by the feasibility of cytological karyotyping. In recent years in silico and molecular approaches have been developed for the genotyping of most widespread inversions (2La, 2Rb and 2Rc). The 2Ru inversion, spanning roughly 8% of chromosome 2R, is commonly polymorphic in West African populations of An. coluzzii and An. gambiae and shows clear increases in frequency with increasing rainfall seasonally and geographically. The aim of this work was to overcome the constraints of currently available cytological and high-throughput molecular assays by developing a simple PCR assay for genotyping the 2Ru inversion in individual specimens of both mosquito species. METHODS: We designed tetra-primer amplification refractory mutation system (ARMS)-PCR assays based on five tag single-nucleotide polymorphisms (SNPs) previously shown to be strongly correlated with 2Ru inversion orientation. The most promising assay was validated against laboratory and field samples of An. coluzzii and An. gambiae karyotyped either cytogenetically or molecularly using a genotyping-in-thousands by sequencing (GT-seq) high-throughput approach that employs targeted sequencing of multiplexed PCR amplicons. RESULTS: A successful assay was designed based on the tag SNP at position 2R, 31710303, which is highly predictive of the 2Ru genotype. The assay, which requires only one PCR, and no additional post-PCR processing other than electrophoresis, produced a clear banding pattern for 98.5% of the 454 specimens tested, which is a 96.7% agreement with established karyotyping methods. Sequences were obtained for nine of the An. coluzzii specimens manifesting 2Ru genotype discrepancies with GT-seq. Possible sources of these discordances are discussed. CONCLUSIONS: The tetra-primer ARMS-PCR assay represents an accurate, streamlined and cost-effective method for the molecular karyotyping of the 2Ru inversion in An. coluzzii and An. gambiae. Together with approaches already available for the other common polymorphic inversions, 2La, 2Rb and 2Rc, this assay will allow investigations of the adaptive value of the complex set of inversion systems observed in the two major malaria vectors in the Afrotropical region.

Mapping Malaria Vector Habitats in West Africa: Drone Imagery and Deep Learning Analysis for Targeted Vector Surveillance.

Disease control programs are needed to identify the breeding sites of mosquitoes, which transmit malaria and other diseases, in order to target interventions and identify environmental risk factors. The increasing availability of very-high-resolution drone data provides new opportunities to find and characterize these vector breeding sites. Within this study, drone images from two malaria-endemic regions in Burkina Faso and Côte d'Ivoire were assembled and labeled using open-source tools. We developed and applied a workflow using region-of-interest-based and deep learning methods to identify land cover types associated with vector breeding sites from very-high-resolution natural color imagery. Analysis methods were assessed using cross-validation and achieved maximum Dice coefficients of 0.68 and 0.75 for vegetated and non-vegetated water bodies, respectively. This classifier consistently identified the presence of other land cover types associated with the breeding sites, obtaining Dice coefficients of 0.88 for tillage and crops, 0.87 for buildings and 0.71 for roads. This study establishes a framework for developing deep learning approaches to identify vector breeding sites and highlights the need to evaluate how results will be used by control programs.

Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito.

Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.

Inferring the epidemiological benefit of indoor vector control interventions against malaria from mosquito data.

The cause of malaria transmission has been known for over a century but it is still unclear whether entomological measures are sufficiently reliable to inform policy decisions in human health. Decision-making on the effectiveness of new insecticide-treated nets (ITNs) and the indoor residual spraying of insecticide (IRS) have been based on epidemiological data, typically collected in cluster-randomised control trials. The number of these trials that can be conducted is limited. Here we use a systematic review to highlight that efficacy estimates of the same intervention may vary substantially between trials. Analyses indicate that mosquito data collected in experimental hut trials can be used to parameterize mechanistic models for Plasmodium falciparum malaria and reliably predict the epidemiological efficacy of quick-acting, neuro-acting ITNs and IRS. Results suggest that for certain types of ITNs and IRS using this framework instead of clinical endpoints could support policy and expedite the widespread use of novel technologies.

Multiple insecticide resistance and first evidence of V410L kdr mutation in Aedes (Stegomyia) aegypti (Linnaeus) from Burkina Faso.

The response to recent dengue outbreaks in Burkina Faso was insecticide-based, despite poor knowledge of the vector population's susceptibility to the insecticides used. Here, we report on the susceptibility to the main insecticide classes and identify important underlying mechanisms in Aedes aegypti populations in Ouagadougou and Banfora, in 2019 and 2020. Wild Ae. aegypti were tested as adults in WHO bioassays and then screened in real time melting curve qPCR analyses to genotype the F1534C, V1016I, and V410L Aedes kdr mutations. Ae. aegypti showed moderate resistance to 0.1% bendiocarb (80-95% survival post-exposure), 0.8% Malathion (60-100%), 0.21% pirimiphos-methyl (75% - 97%), and high resistance to 0.03% deltamethrin (20-70%). PBO pre-exposure partially restored pyrethroid susceptibility. Genotyping detected high frequency of 1534C allele (0.92) and moderate 1016I (0.1-0.32). The V410L mutation was detected in Burkina Faso for the first time (frequency 0.1-0.36). Mosquitoes surviving 4 h exposure to 0.03% deltamethrin had significantly higher frequencies of the F1534C mutation than dead mosquitoes (0.70 vs. 0.96, p < 0.0001) and mosquitoes surviving 2 - 4 h exposure had a significantly reduced life span. Ae. aegypti from Burkina Faso are resistant to multiple insecticide classes with multiple mechanisms involved, demonstrating the essential role of insecticide resistance monitoring within national dengue control programmes.

Optimising the deployment of vector control tools against malaria: a data-informed modelling study.

BACKGROUND: Concern that insecticide resistant mosquitoes are threatening malaria control has driven the development of new types of insecticide treated nets (ITNs) and indoor residual spraying (IRS) of insecticide. Malaria control programmes have a choice of vector control interventions although it is unclear which controls should be used to combat the disease. The study aimed at producing a framework to easily compare the public health impact and cost-effectiveness of different malaria prevention measures currently in widespread use. METHODS: We used published data from experimental hut trials conducted across Africa to characterise the entomological effect of pyrethroid-only ITNs versus ITNs combining a pyrethroid insecticide with the synergist piperonyl butoxide (PBO). We use these estimates to parameterise a dynamic mathematical model of Plasmodium falciparum malaria which is validated for two sites by comparing simulated results to empirical data from randomised control trials (RCTs) in Tanzania and Uganda. We extrapolated model simulations for a series of potential scenarios likely across the sub-Saharan African region and include results in an online tool (Malaria INtervention Tool [MINT]) that aims to identify optimum vector control intervention packages for scenarios with varying budget, price, entomological and epidemiological factors. FINDINGS: Our model indicates that switching from pyrethroid-only to pyrethroid-PBO ITNs could averted up to twice as many cases, although the additional benefit is highly variable and depends on the setting conditions. We project that annual delivery of long-lasting, non-pyrethroid IRS would prevent substantially more cases over 3-years, while pyrethroid-PBO ITNs tend to be the most cost-effective intervention per case averted. The model was able to predict prevalence and efficacy against prevalence in both RCTs for the intervention types tested. MINT is applicable to regions of sub-Saharan Africa with endemic malaria and provides users with a method of designing intervention packages given their setting and budget. INTERPRETATION: The most cost-effective vector control package will vary locally. Models able to recreate results of RCTs can be used to extrapolate outcomes elsewhere to support evidence-based decision making for investment in vector control. FUNDING: Medical Research Council, IVCC, Wellcome Trust. TRANSLATION: For the French translation of the abstract see Supplementary Materials section.

Risk factors associated with house entry of malaria vectors in an area of Burkina Faso with high, persistent malaria transmission and high insecticide resistance.

BACKGROUND: In rural Burkina Faso, the primary malaria vector Anopheles gambiae sensu lato (s.l.) primarily feeds indoors at night. Identification of factors which influence mosquito house entry could lead to development of novel malaria vector control interventions. A study was therefore carried out to identify risk factors associated with house entry of An. gambiae s.l. in south-west Burkina Faso, an area of high insecticide resistance. METHODS: Mosquitoes were sampled monthly during the malaria transmission season using CDC light traps in 252 houses from 10 villages, each house sleeping at least one child aged five to 15 years old. Potential risk factors for house entry of An. gambiae s.l. were measured, including socio-economic status, caregiver's education and occupation, number of people sleeping in the same part of the house as the child, use of anti-mosquito measures, house construction and fittings, proximity of anopheline aquatic habitats and presence of animals near the house. Mosquito counts were compared using a generalized linear mixed-effect model with negative binomial and log link function, adjusting for repeated collections. RESULTS: 20,929 mosquitoes were caught, of which 16,270 (77.7%) were An. gambiae s.l. Of the 6691 An. gambiae s.l. identified to species, 4101 (61.3%) were An. gambiae sensu stricto and 2590 (38.7%) Anopheles coluzzii. Having a metal-roof on the child's sleeping space (IRR = 0.55, 95% CI 0.32-0.95, p = 0.03) was associated with fewer malaria vectors inside the home. CONCLUSION: This study demonstrated that the rate of An. gambiae s.l. was 45% lower in sleeping spaces with a metal roof, compared to those with thatch roofs. Improvements in house construction, including installation of metal roofs, should be considered in endemic areas of Africa to reduce the burden of malaria.

Higher gametocyte production and mosquito infectivity in chronic compared to incident Plasmodium falciparum infections.

Plasmodium falciparum gametocyte kinetics and infectivity may differ between chronic and incident infections. In the current study, we assess parasite kinetics and infectivity to mosquitoes among children (aged 5-10 years) from Burkina Faso with (a) incident infections following parasite clearance (n = 48) and (b) chronic asymptomatic infections (n = 60). In the incident infection cohort, 92% (44/48) of children develop symptoms within 35 days, compared to 23% (14/60) in the chronic cohort. All individuals with chronic infection carried gametocytes or developed them during follow-up, whereas only 35% (17/48) in the incident cohort produce gametocytes before becoming symptomatic and receiving treatment. Parasite multiplication rate (PMR) and the relative abundance of ap2-g and gexp-5 transcripts are positively associated with gametocyte production. Antibody responses are higher and PMR lower in chronic infections. The presence of symptoms and sexual stage immune responses are associated with reductions in gametocyte infectivity to mosquitoes. We observe that most incident infections require treatment before the density of mature gametocytes is sufficient to infect mosquitoes. In contrast, chronic, asymptomatic infections represent a significant source of mosquito infections. Our observations support the notion that malaria transmission reduction may be expedited by enhanced case management, involving both symptom-screening and infection detection.

A PCR-RFLP method for genotyping of inversion 2Rc in Anopheles coluzzii.

BACKGROUND: Genotyping of polymorphic chromosomal inversions in malaria vectors such as An. coluzzii Coetzee & Wilkerson is important, both because they cause cryptic population structure that can mislead vector analysis and control and because they influence epidemiologically relevant eco-phenotypes. The conventional cytogenetic method of genotyping is an impediment because it is labor intensive, requires specialized training, and can be applied only to one gender and developmental stage. Here, we circumvent these limitations by developing a simple and rapid molecular method of genotyping inversion 2Rc in An. coluzzii that is both economical and field-friendly. This inversion is strongly implicated in temporal and spatial adaptations to climatic and ecological variation, particularly aridity. METHODS: Using a set of tag single-nucleotide polymorphisms (SNPs) strongly correlated with inversion orientation, we identified those that overlapped restriction enzyme recognition sites and developed four polymerase chain reaction (PCR) restriction fragment length polymorphism (RFLP) assays that distinguish alternative allelic states at the tag SNPs. We assessed the performance of these assays using mosquito population samples from Burkina Faso that had been cytogenetically karyotyped as well as genotyped, using two complementary high-throughput molecular methods based on tag SNPs. Further validation was performed using mosquito population samples from additional West African (Benin, Mali, Senegal) and Central African (Cameroon) countries. RESULTS: Of four assays tested, two were concordant with the 2Rc cytogenetic karyotype > 90% of the time in all samples. We recommend that these two assays be employed in tandem for reliable genotyping. By accepting only those genotypic assignments where both assays agree, > 99% of assignments are expected to be accurate. CONCLUSIONS: We have developed tandem PCR-RFLP assays for the accurate genotyping of inversion 2Rc in An. coluzzii. Because this approach is simple, inexpensive, and requires only basic molecular biology equipment, it is widely accessible. These provide a crucial tool for probing the molecular basis of eco-phenotypes relevant to malaria epidemiology and vector control.

High-Throughput Genotyping of Common Chromosomal Inversions in the Afrotropical Malaria Mosquito Anopheles Funestus.

Polymorphic chromosomal inversions have been implicated in local adaptation. In anopheline mosquitoes, inversions also contribute to epidemiologically relevant phenotypes such as resting behavior. Progress in understanding these phenotypes and their mechanistic basis has been hindered because the only available method for inversion genotyping relies on traditional cytogenetic karyotyping, a rate-limiting and technically difficult approach that is possible only for the fraction of the adult female population at the correct gonotrophic stage. Here, we focus on an understudied malaria vector of major importance in sub-Saharan Africa, Anopheles funestus. We ascertain and validate tag single nucleotide polymorphisms (SNPs) using high throughput molecular assays that allow rapid inversion genotyping of the three most common An. funestus inversions at scale, overcoming the cytogenetic karyotyping barrier. These same inversions are the only available markers for distinguishing two An. funestus ecotypes that differ in indoor resting behavior, Folonzo and Kiribina. Our new inversion genotyping tools will facilitate studies of ecotypic differentiation in An. funestus and provide a means to improve our understanding of the roles of Folonzo and Kiribina in malaria transmission.

Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms.

Chromosomal inversion polymorphisms have special importance in the Anopheles gambiae complex of malaria vector mosquitoes, due to their role in local adaptation and range expansion. The study of inversions in natural populations is reliant on polytene chromosome analysis by expert cytogeneticists, a process that is limited by the rarity of trained specialists, low throughput, and restrictive sampling requirements. To overcome this barrier, we ascertained tag single nucleotide polymorphisms (SNPs) that are highly correlated with inversion status (inverted or standard orientation). We compared the performance of the tag SNPs using two alternative high throughput molecular genotyping approaches vs. traditional cytogenetic karyotyping of the same 960 individual An. gambiae and An. coluzzii mosquitoes sampled from Burkina Faso, West Africa. We show that both molecular approaches yield comparable results, and that either one performs as well or better than cytogenetics in terms of genotyping accuracy. Given the ability of molecular genotyping approaches to be conducted at scale and at relatively low cost without restriction on mosquito sex or developmental stage, molecular genotyping via tag SNPs has the potential to revitalize research into the role of chromosomal inversions in the behavior and ongoing adaptation of An. gambiae and An. coluzzii to environmental heterogeneities.

Anopheline species composition and the 1014F-genotype in different ecological settings of Burkina Faso in relation to malaria transmission.

BACKGROUND: A three-year longitudinal study was conducted in four sentinel sites from different ecological settings in Burkina Faso, between 2008 and 2010 to identify longitudinal changes in insecticide resistance within Anopheles gambiae complex species based on larval collection. During this study, adult mosquitoes were also collected indoor and outdoor using several methods of collection. The present study reports the diversity of malaria vectors and the 1014F-genotype from this adult collection and investigates the association between this 1014F-genotype and sporozoite rate. METHODS: Adult mosquitoes were collected from July to August (corresponding to the start of rainy season) and October to November (corresponding to the end of rainy season) over 3 years (2008-2010) at four sites across the country, using pyrethrum spray catches (PSC), exit traps and pit shelters. Anopheles gambiae complex mosquitoes were identified to species and genotyped for the L1014F kdr mutation by PCR using genomic DNA. The circumsporozoite antigen of Plasmodium falciparum was detected in mosquitoes using sandwich ELISA. RESULTS: Overall 9212 anopheline mosquitoes were collected during the study period. Of those, 6767 mosquitoes were identified as Anopheles gambiae sensu lato (s.l.). Anopheles arabiensis, Anopheles coluzzii, Anopheles gambiae and or Anopheles funestus were incriminated as vectors of P. falciparum in the study area with an average sporozoite rate of 5%, (95% CI 4.14-5.99%). The kdr1014F-genotype frequencies were 11.44% (95% CI 2.5-39.85%), 19.2% (95% CI 4.53-53.73%) and 89.9 (95% CI 63.14-97.45%), respectively for An. arabiensis, An. coluzzii and An. gambiae. The proportion of the 1014F-genotype varied between sporozoite-infected and uninfected An. gambiae s.l. group. There was no significant difference in the 1014F-genotype frequency between infected and uninfected mosquitoes. CONCLUSION: The current study shows the diversity of malaria vectors and significant interaction between species composition and kdr1014F-genotype in An. gambiae complex mosquitoes from Burkina Faso. In this study, no associations were found between the 1014F-genotype and P. falciparum infection in the major malaria vector An. gambiae s.l.

Efficacy of Olyset Duo, a bednet containing pyriproxyfen and permethrin, versus a permethrin-only net against clinical malaria in an area with highly pyrethroid-resistant vectors in rural Burkina Faso: a cluster-randomised controlled trial.

BACKGROUND: Substantial reductions in malaria incidence in sub-Saharan Africa have been achieved with massive deployment of long-lasting insecticidal nets (LLINs), but pyrethroid resistance threatens control. Burkina Faso is an area with intense malaria transmission and highly pyrethroid-resistant vectors. We assessed the effectiveness of bednets containing permethrin, a pyrethroid, and pyriproxyfen, an insect growth regulator, versus permethrin-only (standard) LLINs against clinical malaria in children younger than 5 years in Banfora, Burkina Faso. METHODS: In this two-group, step-wedge, cluster-randomised, controlled, superiority trial, standard LLINs were incrementally replaced with LLINs treated with permethrin plus pyriproxyfen (PPF) in 40 rural clusters in Burkina Faso. In each cluster, 50 children (aged 6 months to 5 years) were followed up by passive case detection for clinical malaria. Cross-sectional surveys were done at the start and the end of the transmission seasons in 2014 and 2015. We did monthly collections from indoor light traps to estimate vector densities. Primary endpoints were the incidence of clinical malaria, measured by passive case detection, and the entomological inoculation rate. Analyses were adjusted for clustering and for month and health centre. This trial is registered as ISRCTN21853394. FINDINGS: 1980 children were enrolled in the cohort in 2014 and 2157 in 2015. At the end of the study, more than 99% of children slept under a bednet. The incidence of clinical malaria was 2·0 episodes per child-year in the standard LLIN group and 1·5 episodes per child-year in the PPF-treated LLIN group (incidence rate ratio 0·88 [95% CI 0·77-0·99; p=0·04]). The entomological inoculation rate was 85 (95% CI 63-108) infective bites per transmission season in the standard LLIN group versus 42 (32-52) infective bites per transmission season in the PPF-treated LLIN group (rate ratio 0·49, 95% CI 0·32-0·66; p<0·0001). INTERPRETATION: PPF-treated LLINs provide greater protection against clinical malaria than do standard LLINs and could be used as an alternative to standard LLINs in areas with intense transmission of Plasmodium falciparum malaria and highly pyrethroid-resistant vectors. FUNDING: EU Seventh Framework Programme.

Dynamics of the Human Infectious Reservoir for Malaria Determined by Mosquito Feeding Assays and Ultrasensitive Malaria Diagnosis in Burkina Faso.

BACKGROUND: Plasmodium falciparum gametocytes are essential for malaria transmission. Malaria control measures that aim at reducing transmission require an accurate characterization of the human infectious reservoir. METHODS: We longitudinally determined human infectiousness to mosquitoes and P. falciparum carriage by an ultrasensitive RNA-based diagnostics in 130 randomly selected inhabitants of an endemic area. RESULTS: At least 1 mosquito was infected by 32.6% (100 of 307) of the blood samples; in total, 7.6% of mosquitoes (916 of 12 079) were infected. The proportion of infectious individuals and infected mosquitoes were negatively associated with age and positively with asexual parasites (P < .001). Human infectiousness was higher at the start of the wet season and subsequently declined at the peak of the wet season (adjusted odds ratio, 0.52; P = .06) and in the dry season (0.23; P < .001). Overall, microscopy-negative individuals were responsible for 28.7% of infectious individuals (25 of 87) and 17.0% of mosquito infections (145 of 855). CONCLUSIONS: Our study reveals that the infectious reservoir peaks at the start of the wet season, with prominent roles for infections in children and submicroscopic infections. These findings have important consequences for strategies and the timing of interventions, which need to include submicroscopic infections and be implemented in the dry season.

Dynamics of malaria transmission and susceptibility to clinical malaria episodes following treatment of Plasmodium falciparum asymptomatic carriers: results of a cluster-randomized study of community-wide screening and treatment, and a parallel entomology study.

BACKGROUND: In malaria-endemic countries, large proportions of individuals infected with Plasmodium falciparum are asymptomatic and constitute a reservoir of parasites for infection of newly hatched mosquitoes. METHODS: Two studies were run in parallel in Burkina Faso to evaluate the impact of systematic identification and treatment of asymptomatic carriers of P. falciparum, detected by rapid diagnostic test, on disease transmission and susceptibility to clinical malaria episodes. A clinical study assessed the incidence of symptomatic malaria episodes with a parasite density >5,000/µL after three screening and treatment campaigns ~1 month apart before the rainy season; and an entomological study determined the effect of these campaigns on malaria transmission as measured by entomological inoculation rate. RESULTS: The intervention arm had lower prevalence of asymptomatic carriers of asexual parasites and lower prevalence of gametocyte carriers during campaigns 2 and 3 as compared to the control arm. During the entire follow-up period, out of 13,767 at-risk subjects, 2,516 subjects (intervention arm 1,332; control arm 1,184) had symptomatic malaria. Kaplan-Meier analysis of the incidence of first symptomatic malaria episode with a parasite density >5,000/µL showed that, in the total population, the two treatment arms were similar until Week 11-12 after campaign 3, corresponding with the beginning of the malaria transmission season, after which the probability of being free of symptomatic malaria was lower in the intervention arm (logrank p < 0.0001). Similar trends were observed in infants and children <5 years and in individuals ≥5 years of age. In infants and children <5 years old who experienced symptomatic malaria episodes, the geometric mean P. falciparum density was lower in the intervention arm than the control arm. This trend was not seen in those individuals aged ≥5 years. Over the year, monthly variation in mosquito density and entomological inoculation rate was comparable in both arms, with September peaks in both indices. CONCLUSION: Community screening and targeted treatment of asymptomatic carriers of P. falciparum had no effect on the dynamics of malaria transmission, but seemed to be associated with an increase in the treated community's susceptibility to symptomatic malaria episodes after the screening campaigns had finished. These results highlight the importance of further exploratory studies to better understand the dynamics of disease transmission in the context of malaria elimination.