Search

Abstract Background Disease-vector mosquito monitoring is an essential prerequisite to optimize control interventions and evidence-based risk predictions. However, conventional entomological monitoring methods are labor- and time-consuming and do not allow high temporal/spatial resolution. In 2022, a novel system coupling an optical sensor with machine learning technologies (VECTRACK) proved effective in counting and identifying Aedes albopictus and Culex pipiens adult females and males. Here, we carried out the first extensive field evaluation of the VECTRACK system to assess: (i) whether the catching capacity of a commercial BG-Mosquitaire trap (BGM) for adult mosquito equipped with VECTRACK (BGM + VECT) was affected by the sensor; (ii) the accuracy of the VECTRACK algorithm in correctly classifying the target mosquito species genus and sex; (iii) Ae. albopictus capture rate of BGM with or without VECTRACK. Methods The same experimental design was implemented in four areas in northern (Bergamo and Padua districts), central (Rome) and southern (Procida Island, Naples) Italy. In each area, three types of traps—one BGM, one BGM + VECT and the combination of four sticky traps (STs)—were rotated each 48 h in three different sites. Each sampling scheme was replicated three times/area. Collected mosquitoes were counted and identified by both the VECTRACK algorithm and operator-mediated morphological examination. The performance of the VECTRACK system was assessed by generalized linear mixed and linear regression models. Aedes albopictus capture rates of BGMs were calculated based on the known capture rate of ST. Results A total of 3829 mosquitoes (90.2% Ae. albopictus) were captured in 18 collection-days/trap/site. BGM and BGM + VECT showed a similar performance in collecting target mosquitoes. Results show high correlation between visual and automatic identification methods (Spearman Ae. albopictus: females = 0.97; males = 0.89; P

Abstract The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae – are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as “intermediate” between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identified a putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather may have originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, this taxon lacks insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Abstract Background Pyrethroid resistance is one of the major threats for effectiveness of insecticide-treated bed nets (ITNs) in malaria vector control. Genotyping of mutations in the voltage gated sodium channel (VGSC) gene is widely used to easily assess the evolution and spread of pyrethroid target-site resistance among malaria vectors. L1014F and L1014S substitutions are the most common and best characterized VGSC mutations in major African malaria vector species of the Anopheles gambiae complex. Recently, an additional substitution involved in pyrethroid resistance, i.e. V402L, has been detected in Anopheles coluzzii from West Africa lacking any other resistance alleles at locus 1014. The evolution of target-site resistance mutations L1014F/S and V402L was monitored in An. coluzzii and Anopheles arabiensis specimens from a Burkina Faso village over a 10-year range after the massive ITN scale-up started in 2010. Methods Anopheles coluzzii (N = 300) and An. arabiensis (N = 362) specimens collected both indoors and outdoors by different methods (pyrethrum spray catch, sticky resting box and human landing collections) in 2011, 2015 and 2020 at Goden village were genotyped by TaqMan assays and sequencing for the three target site resistance mutations; allele frequencies were statistically investigated over the years. Results A divergent trend in resistant allele frequencies was observed in the two species: 1014F decreased in An. coluzzii (from 0.76 to 0.52) but increased in An. arabiensis (from 0.18 to 0.70); 1014S occurred only in An. arabiensis and slightly decreased over time (from 0.33 to 0.23); 402L increased in An. coluzzii (from 0.15 to 0.48) and was found for the first time in one An. arabiensis specimen. In 2020 the co-occurrence of different resistance alleles reached 43% in An. coluzzii (alleles 410L and 1014F) and 32% in An. arabiensis (alleles 1014F and 1014S). Conclusions Overall, an increasing level of target-site resistance was observed among the populations with only 1% of the two malaria vector species being wild type at both loci, 1014 and 402, in 2020. This, together with the co-occurrence of different mutations in the same specimens, calls for future investigations on the possible synergism between resistance alleles and their phenotype to implement local tailored intervention strategies.

Abstract Modelling approaches play a crucial role in supporting local public health agencies by estimating and forecasting vector abundance and seasonality. However, the reliability of these models is contingent on the availability of standardized, high-quality data. Addressing this need, our study focuses on collecting and harmonizing egg count observations of the mosquito Aedes albopictus, obtained through ovitraps in monitoring and surveillance efforts across Albania, France, Italy, and Switzerland from 2010 to 2022. We processed the raw observations to obtain a continuous time series of ovitraps observations allowing for an extensive geographical and temporal coverage of Ae. albopictus population dynamics. The resulting post-processed observations are stored in the open-access database VectAbundance.This initiative addresses the critical need for accessible, high-quality data, enhancing the reliability of modelling efforts and bolstering public health preparedness.

Abstract 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. Graphical Abstract

Abstract Background Mosquito and human behaviour interaction is a key determinant of the maximum level of protection against malaria that can be provided by insecticide-treated nets (ITNs). Nevertheless, scant literature focuses on this interaction, overlooking a fundamental factor for efficient malaria control. This study aims to estimate malaria transmission risk in a Burkina Faso village by integrating vector biting rhythms with some key information about human habits. Methods Indoor/outdoor human landing catches were conducted for 16 h (16:00–08:00) during 8 nights (September 2020) in Goden village. A survey about net usage and sleeping patterns was submitted to half the households (October–December 2020). A subsample of collected specimens of Anopheles gambiae sensu lato was molecularly processed for species identification, Plasmodium detection from heads-thoraxes and L1014F pyrethroid-resistance allele genotyping. Hourly mosquito abundance was statistically assessed by GLM/GAM, and the entomological inoculation rate (EIR) was corrected for the actual ITN usage retrieved from the questionnaire. Results Malaria transmission was mainly driven by Anopheles coluzzii (68.7%) followed by A. arabiensis (26.2%). The overall sporozoite rate was 2% with L1014F estimated frequency of 0.68 (N = 1070 out of 15,201 A. gambiae s.l. collected). No major shift in mosquito biting rhythms in response to ITN or differences between indoor and outdoor catches were detected. Impressive high biting pressure (mean 30.3 mosquitoes/person/hour) was exerted from 20:00 to 06:00 with a peak at 4:00. Human survey revealed that nearly all inhabitants were awake before 20:00 and after 7:00 and at least 8.7% had no access to bednets. Adjusting for anthropological data, the EIR dropped from 6.7 to 1.2 infective bites/person/16 h. In a scenario of full net coverage and accounting only for the human sleeping patterns, the daily malaria transmission risk not targetable by ITNs was 0.69 infective bites. Conclusions The high mosquito densities and interplay between human/vector activities means that an estimated 10% of residual malaria transmission cannot be prevented by ITNs in the village. Locally tailored studies, like the current one, are essential to explore the heterogeneity of human exposure to infective bites and, consequently, to instruct the adoption of new vector control tools strengthening individual and community protection. Graphical Abstract

BackgroundEmerging arboviral diseases in Europe pose a challenge due to difficulties in detecting and diagnosing cases during the initial circulation of the pathogen. Early outbreak detection enables public health authorities to take effective actions to reduce disease transmission. Quantification of the reporting delays of cases is vital to plan and assess surveillance and control strategies. Here, we provide estimates of reporting delays during an emerging arboviral outbreak and indications on how delays may have impacted onward transmission.Methodology/principal findingsUsing descriptive statistics and Kaplan-Meyer curves we analyzed case reporting delays (the period between the date of symptom onset and the date of notification to the public health authorities) during the 2017 Italian chikungunya outbreak. We further investigated the effect of outbreak detection on reporting delays by means of a Cox proportional hazard model. We estimated that the overall median reporting delay was 15.5 days, but this was reduced to 8 days after the notification of the first case. Cases with symptom onset after outbreak detection had about a 3.5 times higher reporting rate, however only 3.6% were notified within 24h from symptom onset. Remarkably, we found that 45.9% of identified cases developed symptoms before the detection of the outbreak.Conclusions/significanceThese results suggest that efforts should be undertaken to improve the early detection and identification of arboviral cases, as well as the management of vector species to mitigate the impact of long reporting delays.

Abstract Background Colonization of large part of Europe by the Asian tiger mosquito Aedes albopictus is causing autochthonous transmission of chikungunya and dengue exotic arboviruses. While pyrethroids are recommended only to reduce/limit transmission, they are widely implemented to reduce biting nuisance and to control agricultural pests, increasing the risk of insurgence of resistance mechanisms. Worryingly, pyrethroid resistance (with mortality

Malaria is still a leading cause mortality in Côte d’Ivoire despite extensive LLINs coverage. We present the results of an entomological survey conducted in a coastal and in an inland village with the aim to estimate Anopheles gambiae sensu lato (s.l.) female’s abundance indoor/outdoor and Plasmodium falciparum infection rate and analyze the occurrence of blood-feeding in relation to LLINs use. Pyrethrum spray (PSC) and window exit traps (WT) collections were carried out to target endophagic/endophilic and endophagic/exophilic females, respectively. Data on LLINs use in sampled houses were collected. (1) high levels of malaria transmission despite LLINs coverage >70% (~1 An. gambiae s.l. predicted mean/person/night and ~5% Plasmodium falciparum infection rate); (2) 46% of females in the PSC sample were blood-fed, suggesting that they fed on an unprotected host inside the house; (3) 81% of females in WT were unfed, suggesting that they were leaving the house to find an available host. Model estimates that if everyone sleeps under LLINs the probability for a mosquito to bite decreases of 48% and 95% in the coastal and inland village, respectively. The results show a high proportion of mosquito biting and resting indoors despite extensive LLINs. The biological/epidemiological determinants of accounting for these results merit deeper investigations.

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

The Mosquito Alert dataset includes occurrence records of adult mosquitoes collected worldwide in 2014–2020 through Mosquito Alert, a citizen science system for investigating and managing disease-carrying mosquitoes. Records are linked to citizen science-submitted photographs and validated by entomologists to determine the presence of five targeted European mosquito vectors: Aedes albopictus, Ae. aegypti, Ae. japonicus, Ae. koreicus, and Culex pipiens. Most records are from Spain, reflecting Spanish national and regional funding, but since autumn 2020 substantial records from other European countries are included, thanks to volunteer entomologists coordinated by the AIM-COST Action, and to technological developments to increase scalability. Among other applications, the Mosquito Alert dataset will help develop citizen science-based early warning systems for mosquito-borne disease risk. It can also be reused for modelling vector exposure risk, or to train machine-learning detection and classification routines on the linked images, to assist with data validation and establishing automated alert systems.

Human and animal vector-borne diseases, particularly mosquito-borne diseases, are emerging or re-emerging worldwide. Six Aedes invasive mosquito (AIM) species were introduced to Europe since the 1970s: Aedes aegypti, Ae. albopictus, Ae. japonicus, Ae. koreicus, Ae. atropalpus and Ae. triseriatus. Here, we report the results of AIMSurv2020, the first pan-European surveillance effort for AIMs. Implemented by 42 volunteer teams from 24 countries. And presented in the form of a dataset named “AIMSurv Aedes Invasive Mosquito species harmonized surveillance in Europe. AIM-COST Action. Project ID: CA17108”. AIMSurv2020 harmonizes field surveillance methodologies for sampling different AIMs life stages, frequency and minimum length of sampling period, and data reporting. Data include minimum requirements for sample types and recommended requirements for those teams with more resources. Data are published as a Darwin Core archive in the Global Biodiversity Information Facility- Spain, comprising a core file with 19,130 records (EventID) and an occurrences file with 19,743 records (OccurrenceID). AIM species recorded in AIMSurv2020 were Ae. albopictus, Ae. japonicus and Ae. koreicus, as well as native mosquito species.

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.

Abstract Background The spatial spread of many mosquito-borne diseases occurs by focal spread at the scale of a few hundred meters and over longer distances due to human mobility. The relative contributions of different spatial scales for transmission of chikungunya virus require definition to improve outbreak vector control recommendations. Methods We analyzed data from a large chikungunya outbreak mediated by the mosquito Aedes albopictus in the Lazio region, Italy, consisting of 414 reported human cases between June and November 2017. Using dates of symptom onset, geographic coordinates of residence, and information from epidemiological questionnaires, we reconstructed transmission chains related to that outbreak. Results Focal spread (within 1 km) accounted for 54.9% of all cases, 15.8% were transmitted at a local scale (1–15 km) and the remaining 29.3% were exported from the main areas of chikungunya circulation in Lazio to longer distances such as Rome and other geographical areas. Seventy percent of focal infections (corresponding to 38% of the total 414 cases) were transmitted within a distance of 200 m (the buffer distance adopted by the national guidelines for insecticide spraying). Two main epidemic clusters were identified, with a radius expanding at a rate of 300–600 m per month. The majority of exported cases resulted in either sporadic or no further transmission in the region. Conclusions Evidence suggest that human mobility contributes to seeding a relevant number of secondary cases and new foci of transmission over several kilometers. Reactive vector control based on current guidelines might allow a significant number of secondary clusters in untreated areas, especially if the outbreak is not detected early. Existing policies and guidelines for control during outbreaks should recommend the prioritization of preventive measures in neighboring territories with known mobility flows to the main areas of transmission.

Abstract Background Despite the overall major impact of long-lasting insecticide treated nets (LLINs) in eliciting individual and collective protection to malaria infections, some sub-Saharan countries, including Burkina Faso, still carry a disproportionately high share of the global malaria burden. This study aims to analyse the possible entomological bases of LLIN limited impact, focusing on a LLIN-protected village in the Plateau Central region of Burkina Faso. Methods Human landing catches (HLCs) were carried out in 2015 for 12 nights both indoors and outdoors at different time windows during the highest biting activity phase for Anopheles gambiae (s.l.). Collected specimens were morphologically and molecularly identified and processed for Plasmodium detection and L1014F insecticide-resistance allele genotyping. Results Almost 2000 unfed An. gambiae (s.l.) (54% Anopheles coluzzii and 44% Anopheles arabiensis) females landing on human volunteers were collected, corresponding to a median number of 23.5 females/person/hour. No significant differences were observed in median numbers of mosquitoes collected indoors and outdoors, nor between sporozoite rates in An. coluzzii (6.1%) and An. arabiensis (5.5%). The estimated median hourly entomological inoculation rate (EIR) on volunteers was 1.4 infective bites/person/hour. Results do not show evidence of the biting peak during night hours typical for An. gambiae (s.l.) in the absence of bednet protection. The frequency of the L1014F resistant allele (n = 285) was 66% in An. coluzzii and 38% in An. arabiensis. Conclusions The observed biting rate and sporozoite rates are in line with the literature data available for An. gambiae (s.l.) in the same geographical area before LLIN implementation and highlight high levels of malaria transmission in the study village. Homogeneous biting rate throughout the night and lack of preference for indoor-biting activity, suggest the capacity of both An. coluzzii and An. arabiensis to adjust their host-seeking behaviour to bite humans despite bednet protection, accounting for the maintenance of high rates of mosquito infectivity and malaria transmission. These results, despite being limited to a local situation in Burkina Faso, represent a paradigmatic example of how high densities and behavioural plasticity in the vector populations may contribute to explaining the limited impact of LLINs on malaria transmission in holo-endemic Sudanese savannah areas in West Africa.

Abstract Background Chromosomal inversion polymorphisms play a role in adaptation to heterogeneous environments. Inversion polymorphisms are implicated in the very high ecological flexibility of the three main malaria vector species of the Afrotropical Anopheles gambiae complex, facilitating the exploitation of anthropogenic environmental modifications and promoting a strong association with humans. In addition to extending the species’ spatial and temporal distribution, inversions are associated with epidemiologically relevant mosquito behavior and physiology, underscoring their medical importance. We here present novel PCR-RFLP based assays strongly predictive of genotype for the cosmopolitan 2Rb inversion in An. coluzzii and An. gambiae, a development which overcomes the numerous constraints inherent to traditional cytological karyotyping. Methods We designed PCR-RFLP genotyping assays based on tag SNPs previously computationally identified as strongly predictive (> 95%) of 2Rb genotype. We targeted those tags whose alternative allelic states destroyed or created the recognition site of a commercially available restriction enzyme, and designed assays with distinctive cleavage profiles for each inversion genotype. The assays were validated on 251 An. coluzzii and 451 An. gambiae cytologically karyotyped specimens from nine countries across Africa and one An. coluzzii laboratory colony. Results For three tag SNPs, PCR-RFLP assays (denoted DraIII, MspAI, and TatI) reliably produced robust amplicons and clearly distinguishable electrophoretic profiles for all three inversion genotypes. Results obtained with the DraIII assay are ≥ 95% concordant with cytogenetic assignments in both species, while MspAI and TatI assays produce patterns highly concordant with cytogenetic assignments only in An. coluzzii or An. gambiae, respectively. Joint application of species-appropriate pairs of assays increased the concordance levels to > 99% in An. coluzzii and 98% in An. gambiae. Potential sources of discordance (e.g. imperfect association between tag and inversion, allelic dropout, additional polymorphisms in the restriction target site, incomplete or failed restriction digestion) are discussed. Conclusions The availability of highly specific, cost effective and accessible molecular assays for genotyping 2Rb in An. gambiae and An. coluzzii allows karyotyping of both sexes and all developmental stages. These novel tools will accelerate deeper investigations into the role of this ecologically and epidemiologically important chromosomal inversion in vector biology.

In the last decades, the colonization of Mediterranean Europe and of other temperate regions by Aedes albopictus created an unprecedented nuisance problem in highly infested areas and new public health threats due to the vector competence of the species. The Sterile Insect Technique (SIT) and the Incompatible Insect Technique (IIT) are insecticide-free mosquito-control methods, relying on mass release of irradiated/manipulated males, able to complement existing and only partially effective control tools. The validation of these approaches in the field requires appropriate experimental settings, possibly isolated to avoid mosquito immigration from other infested areas, and preliminary ecological and entomological data. We carried out a 4-year study in the island of Procida (Gulf of Naples, Italy) in strict collaboration with local administrators and citizens to estimate the temporal dynamics, spatial distribution, and population size of Ae. albopictus and the dispersal and survival of irradiated males. We applied ovitrap monitoring, geo-spatial analyses, mark-release-recapture technique, and a citizen-science approach. Results allow to predict the seasonal (from April to October, with peaks of 928-9,757 males/ha) and spatial distribution of the species, highlighting the capacity of Ae. albopictus population of Procida to colonize and maintain high frequencies in urban as well as in sylvatic inhabited environments. Irradiated males shown limited ability to disperse (mean daily distance travelled

Anopheles coluzzii and An. gambiae—the two most recently radiated species of the An. gambiae complex and the major Afrotropical malaria vector species—are identified by markers in the X-centromeric IGS rDNA region. Putative IGS-hybrids are rarely found in the field, except in restricted areas where genomic studies have led to the hypothesis that the observed IGS-patterns are due to cryptic taxa rather than to hybridization between the two species. We investigated the genome-wide levels of admixture in two villages in Côte d’Ivoire where high levels of IGS-hybrids have been detected, confirming unparalleled high frequencies in the coastal village. Genotyping of 24 Ancestry Informative Markers (AIMs) along the three chromosomes produced discordant results between the IGS-marker and the multilocus genotype obtained for AIMs across the whole genome (29%) as well as AIMs on chromosome-X (considered to be fundamental for species reproductive isolation) only (21%). Results highlight a complicated pattern of admixture that deserves deeper genomic analyses to understand better possible underlying causes (from extensive processes of hybridization to the existence of different cryptic taxa), and stress the need of developing advanced diagnostics for An. coluzzii, An. gambiae and putative new taxa, instrumental for assessing taxon-specific epidemiological characters.

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained.

Abstract The increasing number of exotic arbovirus cases imported in Europe and the 2017 chikungunya outbreak in central/southern Italy highlight the urgency of evidence-based outbreak management plans to predict, prevent or interrupt spreading of these arboviruses to non-endemic countries in temperate regions. We here present the results of three mark-release-recapture experiments conducted in a peri-urban area of North-East Italy to estimate the spatio-temporal dynamics of the dispersal of Aedes albopictus females looking for oviposition sites. The Flight Range of 90% of the mosquito population (FR90) was found to exceed 200 m, consistently with data obtained from a previous study conducted in a highly urbanised area in Rome (Central Italy). Modelling results showed that dispersal can be so rapid that insecticide spraying within a 200m-radius around a potential infected case leaves >10% probability that a potentially infected mosquito escapes the treatment, even if this is carried out after only 2–3 days since the importation of a viremic case. These data provide evidence in favour of an update of guidelines for the control of exotic autochthonous arbovirus transmission in temperate areas and highlight the need of effective surveillance approaches and rapid response to contain the risks associated to imported viremic cases.

Mosquitoes-borne diseases are major public health issues particularly in Africa. Vector control interventions and human-made environmental/climatic changes significantly affect the distribution and abundance of vector species. We carried out an entomological survey targeting host-seeking mosquitos in two different ecological contexts—coastal and inland—in Senegal, by CDC-light and BG-sentinel traps. Results show high predominance of Culex quinquefasciatus (90%) and of Anopheles arabiensis within malaria vectors (46%), with mean numbers of females/trap/nights =8 and Cx. quinquefasciatus and of CDC traps for An. gambiae s.l. and highlight that, despite both traps target the host-seeking fraction of the population, they provide different patterns of species abundance, temporal dynamics and host-seeking activity, leading to possible misinterpretation of the species bionomics. This draws attention to the need of taking into account trapping performance, in order to provide realistic quantification of the number of mosquitoes per units of space and time, the crucial parameter for evaluating vector–human contact, and estimating risk of pathogen transmission.

BackgroundIn the last four decades, the Asian tiger mosquito, Aedes albopictus, vector of several human arboviruses, has spread from its native range in South-East Asia to all over the world, largely through the transportation of its eggs via the international trade in used tires. Albania was the first country invaded in Europe in 1979, followed by Italy in 1990 and other Mediterranean countries after 2000.Methods/principal findingsWe here inferred the invasion history and migration patterns of Ae. albopictus in Italy (today the most heavily-infested country in Europe), Greece and Albania, by analyzing a panel of >100,000 single nucleotide polymorphisms (SNPs) obtained by sequencing of double-digest Restriction site-Associated DNA (ddRADseq). The obtained dataset was combined with samples previously analyzed from both the native and invasive range worldwide to interpret the results using a broader spatial and historical context. The emerging evolutionary scenario complements the results of other studies in showing that the extraordinary worldwide expansion of Ae. albopictus has occurred thanks to multiple independent invasions by large numbers of colonists from multiple geographic locations in both native and previously invaded areas, consistently with the role of used tires shipments to move large numbers of eggs worldwide. By analyzing mosquitoes from nine sites across ~1,000-km transect in Italy, we were able to detect a complex interplay of drift, isolation by distance mediated divergence, and gene flow in shaping the species very recent invasion and range expansion, suggesting overall high connectivity, likely due to passive transportation of adults via ground transportation, as well as specific adaptations to local conditions.Conclusions/significanceResults contribute to characterize one of the most successful histories of animal invasion, and could be used as a baseline for future studies to track epidemiologically relevant characters (e.g. insecticide resistance).

Polymerase chain reaction (PCR)-based genotyping of mutations in the voltage-sensitive sodium channel (vssc) associated with resistance to pyrethroid insecticides is widely used and represents a potential early warning and monitoring system for insecticide resistance arising in mosquito populations, which are vectors of different human pathogens. In the secondary vector Aedes albopictus—an Asian species that has invaded and colonized the whole world, including temperate regions—sequencing of domain II of the vssc gene is still needed to detect the V1016G mutation associated with pyrethroid resistance. In this study we developed and tested a novel allele-specific PCR (AS-PCR) assay to genotype the V1016G mutation in this species and applied it to the analysis of wild populations from Italy. The results confirm the high accuracy of the novel AS-PCR and highlight frequencies of the V1016G allele as >5% in most sampling sites, with peaks of 20–45% in coastal touristic sites where pyrethroid treatments are extensively implemented, mostly for mosquito nuisance reduction. The high frequency of this mutation observed in Italian Ae. albopictus populations should serve as a warning bell, advocating for increased monitoring and management of a phenomenon which risks neutralizing the only weapon today available to counteract (risks of) arbovirus outbreaks.

The increasing spread of the Asian tiger mosquito, Aedes albopictus, in Europe and US raises public health concern due to the species competence to transmit several exotic human arboviruses, among which dengue, chikungunya and Zika, and urges the development of suitable modeling approach to forecast the spatial and temporal distribution of the mosquito. Here we developed a dynamical species distribution modeling approach forecasting Ae. albopictus eggs abundance at high spatial (0.01 degree WGS84) and temporal (weekly) resolution over 10 Balkan countries, using temperature times series of Modis data products and altitude as input predictors. The model was satisfactorily calibrated and validated over Albania based observed eggs abundance data weekly monitored during three years. For a given week of the year, eggs abundance was mainly predicted by the number of eggs and the mean temperature recorded in the preceding weeks. That is, results are in agreement with the biological cycle of the mosquito, reflecting the effect temperature on eggs spawning, maturation and hatching. The model, seeded by initial egg values derived from a second model, was then used to forecast the spatial and temporal distribution of eggs abundance over the selected Balkan countries, weekly in 2011, 2012 and 2013. The present study is a baseline to develop an easy-handling forecasting model able to provide information useful for promoting active surveillance and possibly prevention of Ae. albopictus colonization in presently non-infested areas in the Balkans as well as in other temperate regions.

Abstract Background Urban mosquitoes in temperate regions may represent a high nuisance and are associated with the risk of arbovirus transmission. Common practices to reduce this burden, at least in Italian highly infested urban areas, imply calendar-based larvicide treatments of street catch basins – which represent the main non-removable urban breeding site – and/or insecticide ground spraying. The planning of these interventions, as well as the evaluation of their effectiveness, rarely benefit of adequate monitoring of the mosquito abundance and dynamics. We propose the use of adhesive traps to monitor Aedes albopictus and Culex pipiens adults and to evaluate the efficacy of insecticide-based control strategies. Methods We designed two novel types of adhesive traps to collect adult mosquitoes visiting and/or emerging from catch basins. The Mosquito Emerging Trap (MET) was exploited to assess the efficacy of larvicide treatments. The Catch Basin Trap (CBT) was exploited together with the Sticky Trap (ST, commonly used to collect ovipositing/resting females) to monitor adults abundance in the campus of the University of Rome “Sapienza” - where catch basins were treated with Insect Growth Regulators (IGR) bi-monthly and Low-Volume insecticide spraying were carried out before sunset - and in a nearby control area. Results Results obtained by MET showed that, although all monitored diflubenzuron-treated catch basins were repeatedly visited by Ae. albopictus and Cx. pipiens, adult emergence was inhibited in most basins. Results obtained by ST and CBT showed a significant lower adult abundance in the treated area than in the untreated one after the second adulticide spraying, which was carried out during the major phase of Ae. albopictus population expansion in Rome. Spatial heterogeneities in the effect of the treatments were also revealed. Conclusions The results support the potential of the three adhesive traps tested in passively monitoring urban mosquito adult abundance and seasonal dynamics and in assessing the efficacy of control measures. ST showed higher specificity for Ae. albopictus and CBT for Cx. pipiens. The results also provide a preliminary indication on the effectiveness of common mosquito control strategies carried out against urban mosquito in European urban areas.

Abstract Background Despite the expanding worldwide distribution of Aedes albopictu s and its increasing relevance as arboviral vector, current methods to collect adult specimens are not optimal. Improved approaches are thus needed to monitor their density and pathogen infections, and to establish baseline data for control interventions. A widely used device is the BG-Sentinel (BG-trap) which mostly targets host-seeking females attracted by release of CO2 and/or a synthetic odour blend (the BG lure). We compared the attractiveness of this blend to that of the Mbita (MB5) lure, a new synthetic blend of proven efficiency in attracting Afrotropical malaria vectors, and evaluated the additional effect of CO2 to the two odour baits. Findings We carried out 6x6 Latin square experiments in two Ae. albopictus-infested areas in Rome, baiting the BG-traps as follows: CO2, BG lure, MB5 lure, BG lure + CO2, MB5 lure + CO2, no bait. CO2 was derived from yeast-fermented sugar. Overall, 949 females and 816 males were collected. Baited traps collected significantly more females than unbaited ones. Traps baited with either lures in combination with CO2 were more effective than those baited with CO2 alone. No significant differences were observed in female captures between traps baited with any of the two lures, nor between the two lures, independently from the addition of CO2. The use of BG lure + CO2 significantly increased males catches compared to unbaited traps. Conclusions The results suggest a broad significance of the MB5 lure for sampling medically important mosquito species and highlight the high efficacy of the combination of lures + CO2 for female Ae. albopictus and of BG lure + CO2 for males, leading to consider CO2 as an essential additional cue for the sampling of this species.

Background Aedes albopictus is an aggressive invasive mosquito species that represents a serious health concern not only in tropical areas, but also in temperate regions due to its role as vector of arboviruses. Estimates of mosquito biting rates are essential to account for vector-human contact in models aimed to predict the risk of arbovirus autochthonous transmission and outbreaks, as well as nuisance thresholds useful for correct planning of mosquito control interventions. Methods targeting daytime and outdoor biting Ae. albopictus females (e.g., Human Landing Collection, HLC) are expensive and difficult to implement in large scale schemes. Instead, egg-collections by ovitraps are the most widely used routine approach for large-scale monitoring of the species. The aim of this work was to assess whether ovitrap data can be exploited to estimate numbers of adult biting Ae. albopictus females and whether the resulting relationship could be used to build risk models helpful for decision-makers in charge of planning of mosquito-control activities in infested areas. Method Ovitrap collections and HLCs were carried out in hot-spots of Ae. albopictus abundance in Rome (Italy) along a whole reproductive season. The relationship between the two sets of data was assessed by generalized least square analysis, taking into account meteorological parameters. Result The mean number of mosquito females/person collected by HLC in 15′ (i.e., females/HLC) and the mean number of eggs/day were 18.9 ± 0.7 and 39.0 ± 2.0, respectively. The regression models found a significant positive relationship between the two sets of data and estimated an increase of one biting female/person every five additional eggs found in ovitraps. Both observed and fitted values indicated presence of adults in the absence of eggs in ovitraps. Notably, wide confidence intervals of estimates of biting females based on eggs were observed. The patterns of exotic arbovirus outbreak probability obtained by introducing these estimates in risk models were similar to those based on females/HLC (R0 > 1 in 86% and 40% of sampling dates for Chikungunya and Zika, respectively; R0 1 for Chikungunya is also to be expected when few/no eggs/day are collected by ovitraps. Discussion This work provides the first evidence of the possibility to predict mean number of adult biting Ae. albopictus females based on mean number of eggs and to compute the threshold of eggs/ovitrap associated to epidemiological risk of arbovirus transmission in the study area. Overall, however, the large confidence intervals in the model predictions represent a caveat regarding the reliability of monitoring schemes based exclusively on ovitrap collections to estimate numbers of biting females and plan control interventions.

Abstract Background Monitoring densities of adult mosquito populations is a major challenge in efforts to evaluate the epidemiology of mosquito-borne diseases, and their response to vector control interventions. In the case of malaria, collection of outdoor-resting Anophelines is rarely incorporated into surveillance and control, partially due to the lack of standardized collection tools. Such an approach, however, is increasingly important to investigate possible changes in mosquito behaviour in response to the scale up of Insecticide Treated Nets and Indoor Residual Spraying. In this study we evaluated the Sticky Resting Box (SRB) - i.e. a sticky variant of previously investigated mosquito Resting Box, which allows passive collection of mosquitoes entering the box – and compared its performance against traditional methods for indoor and outdoor resting mosquito sampling. Methods Daily collections were carried out in two neighbouring villages of Burkina Faso during rainy season 2011 and dry season 2012 by SRB located indoors and outdoors, and by Back-Pack aspiration inside houses (BP) and in ad hoc built outdoor pit-shelters (PIT). Results Overall, almost 20,000 Culicidae specimens belonging to 16 species were collected and morphologically identified. Malaria vectors included Anopheles coluzzii (53%), An. arabiensis (12%), An. gambiae s.s. (2.0%) and An. funestus (4.5%). The diversity of species collected in the two villages was similar for SRB and PIT collections outdoors, and significantly higher for SRB than for BP indoors. The population dynamics of An. gambiae s.l. mosquitoes, as obtained by SRB-collections was significantly correlated with those obtained by the traditional methods. The predicted mean estimates of An. gambiae s.l. specimens/sampling-unit/night-of-collections was 6- and 5-times lower for SRB than for BP indoors and PIT outdoors, respectively. Conclusions Overall, the daily performance of SRB in terms of number of malaria vectors/trap was lower than that of traditionally used approaches for in- and outdoor collections. However, unlike these methods, SRB could be set up to collect mosquitoes passively over at least a week. This makes SRB a promising tool for passively monitoring anopheline resting populations, with data presented here providing guidance for how to set up SRB-based collections to acquire information comparable to those obtained with other methods.

Aedes albopictus is a tropical invasive species which in the last decades spread worldwide, also colonizing temperate regions of Europe and US, where it has become a public health concern due to its ability to transmit exotic arboviruses, as well as severe nuisance problems due to its aggressive daytime outdoor biting behaviour. While several studies have been carried out in order to predict the potential limits of the species expansions based on eco-climatic parameters, few studies have so far focused on the specific effects of these variables in shaping its micro-geographic abundance and dynamics. The present study investigated eco-climatic factors affecting Ae. albopictus abundance and dynamics in metropolitan and sub-urban/rural sites in Rome (Italy), which was colonized in 1997 and is nowadays one of the most infested metropolitan areas in Southern Europe. To this aim, longitudinal adult monitoring was carried out along a 70 km-transect across and beyond the most urbanized and densely populated metropolitan area. Two fine scale spatiotemporal datasets (one with reference to a 20m circular buffer around sticky traps used to collect mosquitoes and the second to a 300m circular buffer within each sampling site) were exploited to analyze the effect of climatic and socio-environmental variables on Ae. albopictus abundance and dynamics along the transect. Results showed an association between highly anthropized habitats and high adult abundance both in metropolitan and sub-urban/rural areas, with "small green islands" corresponding to hot spots of abundance in the metropolitan areas only, and a bimodal seasonal dynamics with a second peak of abundance in autumn, due to heavy rains occurring in the preceding weeks in association with permissive temperatures. The results provide useful indications to prioritize public mosquito control measures in temperate urban areas where nuisance, human-mosquito contact and risk of local arbovirus transmission are likely higher, and highlight potential public health risks also after the summer months typically associated with high mosquito densities.

Seasonal-long larvicide treatments and/or outdoor space-spray applications of insecticides are frequently applied to reduce Aedes albopictus nuisance in urban areas of temperate regions, where the species has become a permanent pest affecting people's quality of life and health. However, assessments of the effectiveness of sequential interventions is a difficult task, as it requires to take into account the cumulative and combined effect of multiple treatments, as well as the mosquito seasonal dynamics (rather than mosquito abundance before and after single treatments). We here present the results of the effectiveness assessment of a seasonal-long calendar-based control intervention integrating larvicide treatments of street catch basins and night-time adulticide ground spraying in the main University hospital in Rome (Italy). Cage-experiments and an intensive monitoring of wild mosquito abundance in treated and untreated sites were carried out along an entire season. Sticky traps were used to monitor adult abundance and site-specific eco-climatic variations (by recording water left over in each trap), in order to disentangle the effect of insecticide treatments from eco-climatic drivers on mosquito seasonal dynamics. Despite the apparent limited impact of single adulticide sprayings assessed based on mortality in caged and wild mosquitoes, the results of the temporal analysis showed that mosquito seasonal patterns were initially comparable in the two sites, diverged in the absence of diverging eco-climatic conditions and remained stable afterwards. This allowed to attribute the lack of the expected Ae. albopictus population expansion in the treated site to the combined effect of multiple adulticide sprayings and larvicide treatments carried out during the whole season. The approach proposed was proved to be successful to assess effects of seasonal-long control treatments on adult mosquito population dynamics and could represent a valuable instrument to assess the effectiveness of other control interventions, to evaluate their actual cost-benefits and to possibly minimize space-spraying applications to reduce mosquito nuisance.

"Far-West" Africa is known to be a secondary contact zone between the two major malaria vectors Anopheles coluzzii and A. gambiae. We investigated gene-flow and potentially adaptive introgression between these species along a west-to-east transect in Guinea Bissau, the putative core of this hybrid zone. To evaluate the extent and direction of gene flow, we genotyped site 702 in Intron-1 of the para Voltage-Gated Sodium Channel gene, a species-diagnostic nucleotide position throughout most of A. coluzzii and A. gambiae sympatric range. We also analyzed polymorphism in the thioester-binding domain (TED) of the innate immunity-linked thioester-containing protein 1 (TEP1) to investigate whether elevated hybridization might facilitate the exchange of variants linked to adaptive immunity and Plasmodium refractoriness. Our results confirm asymmetric introgression of genetic material from A. coluzzii to A. gambiae and disruption of linkage between the centromeric "genomic islands" of inter-specific divergence. We report that A. gambiae from the Guinean hybrid zone possesses an introgressed TEP1 resistant allelic class, found exclusively in A. coluzzii elsewhere and apparently swept to fixation in West Africa (i.e. Mali and Burkina Faso). However, no detectable fixation of this allele was found in Guinea Bissau, which may suggest that ecological pressures driving segregation between the two species in larval habitats in this region may be different from those experienced in northern and more arid parts of the species' range. Finally, our results also suggest a genetic subdivision between coastal and inland A. gambiae Guinean populations and provide clues on the importance of ecological factors in intra-specific differentiation processes.

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are small soluble polypeptides that bind semiochemicals in the lymph of insect chemosensilla. In the genome of Anopheles gambiae, 66 genes encode OBPs and 8 encode CSPs. Here we monitored their expression through classical proteomics (2D gel-MS analysis) and a shotgun approach. The latter method proved much more sensitive and therefore more suitable for tiny biological samples as mosquitoes antennae and eggs. Females express a larger number and higher quantities of OBPs in their antennae than males (24 vs 19). OBP9 is the most abundant in the antennae of both sexes, as well as in larvae, pupae and eggs. Of the 8 CSPs, 4 were detected in antennae, while SAP3 was the only one expressed in larvae. Our proteomic results are in fairly good agreement with data of RNA expression reported in the literature, except for OBP4 and OBP5, that we could not identify in our analysis, nor could we detect in Western Blot experiments. The relatively limited number of soluble olfactory proteins expressed at relatively high levels in mosquitoes makes further studies on the coding of chemical messages at the OBP level more accessible, providing for few specific targets. Identification of such proteins in Anopheles gambiae might facilitate future studies on host finding behavior in this important disease vector.

The main constraint to the fight against container-breeding mosquito vectors of human arboviruses is the difficulty in targeting the multiplicity of larval sources, mostly represented by small man-made water containers. The aim of this work is to assess the feasibility of the "auto-dissemination" approach, already tested for Aedes aegypti, as a possible alternative to traditional, inefficient control tools, against Ae. albopictus in urban areas. The approach is based on the possibility that wild adult females, exposed to artificial resting sites contaminated with pyriproxyfen, can disseminate this juvenile hormone analogue to larval habitats, thus interfering with adult emergence.We carried out four field experiments in two areas of Rome that are typically highly infested with Ae. albopictus, i.e. the main cemetery and a small green area within a highly urbanised neighbourhood. In each area we used 10 pyriproxyfen "dissemination" stations, 10 "sentinel" sites and 10 covered, control sites. The sentinel and control sites each contained 25 Ae. albopictus larvae. These were monitored for development and adult emergence.When a 5% pyriproxyfen powder was used to contaminate the dissemination sites, we observed significantly higher mortality at the pupal stage in the sentinel sites (50-70%) than in the controls (

The main Afrotropical malaria vector, Anopheles gambiae sensu stricto, is undergoing a process of sympatric ecological diversification leading to at least two incipient species (the M and S molecular forms) showing heterogeneous levels of divergence across the genome. The physically unlinked centromeric regions on all three chromosomes of these closely related taxa contain fixed nucleotide differences which have been found in nearly complete linkage disequilibrium in geographic areas of no or low M-S hybridization. Assays diagnostic for SNP and structural differences between M and S forms in the three centromeric regions were applied in samples from the western extreme of their range of sympatry, the only area where high frequencies of putative M/S hybrids have been reported. The results reveal a level of admixture not observed in the rest of the range. In particular, we found: i) heterozygous genotypes at each marker, although at frequencies lower than expected under panmixia; ii) virtually all possible genotypic combinations between markers on different chromosomes, although genetic association was nevertheless detected; iii) discordant M and S genotypes at two X-linked markers near the centromere, suggestive of introgression and inter-locus recombination. These results could be indicative either of a secondary contact zone between M and S, or of the maintenance of ancestral polymorphisms. This issue and the perspectives opened by these results in the study of the M and S incipient speciation process are discussed.

MALDI profiling and imaging mass spectrometry (IMS) are novel techniques for direct analysis of peptides and small proteins in biological tissues. In this work we applied them to the study of Anopheles gambiae antennae, with the aim of analysing expression of soluble proteins involved in olfaction perireceptor events. MALDI spectra obtained by direct profiling on single antennae and by the analysis of extracts, showed similar profiles, although spectra obtained through profiling had a richer ion population and higher signal to noise ratio. Male and female antennae showed distinct protein profiles. MALDI imaging experiments were also performed and differences were observed in the localization of some proteins. Two proteins were identified through high resolution measurement and top-down MS/MS experiments. A 8 kDa protein only present in the male antennae matched with an unannotated sequence of the An. gambiae genome, while the presence of odorant binding protein 9 (OBP-9) was confirmed through experiments of 2-DE, followed by MS and MS/MS analysis of digested spots. This work shows that MALDI MS profiling is a technique suitable for the analysis of proteins of small and medium MW in insect appendices, and allows obtaining data for several specimens which can be investigated for differences between groups. Proteins of interest can be identified through other complementary MS approaches.

How often insecticide resistance mutations arise in natural insect populations is a fundamental question for understanding the evolution of resistance and also for modeling its spread. Moreover, the development of resistance is regarded as a favored model to study the molecular evolution of adaptive traits. In the malaria vector Anopheles gambiae two point mutations (L1014F and L1014S) in the voltage-gated sodium channel gene, that confer knockdown resistance (kdr) to DDT and pyrethroid insecticides, have been described. In order to determine whether resistance alleles result from single or multiple mutation events, genotyping of the kdr locus and partial sequencing of the upstream intron-1 was performed on a total of 288 A. gambiae S-form collected from 28 localities in 15 countries. Knockdown resistance alleles were found to be widespread in West Africa with co-occurrence of both 1014S and 1014F in West-Central localities. Differences in intron-1 haplotype composition suggest that kdr alleles may have arisen from at least four independent mutation events. Neutrality tests provided evidence for a selective sweep acting on this genomic region, particularly in West Africa. The frequency and distribution of these kdr haplotypes varied geographically, being influenced by an interplay between different mutational occurrences, gene flow and local selection. This has important practical implications for the management and sustainability of malaria vector control programs.