The effective population size of malaria mosquitoes: large impact of vector control.

Malaria vectors in sub-Saharan Africa have proven themselves very difficult adversaries in the global struggle against malaria. Decades of anti-vector interventions have yielded mixed results--with successful reductions in transmission in some areas and limited impacts in others. These varying successes can be ascribed to a lack of universally effective vector control tools, as well as the development of insecticide resistance in mosquito populations. Understanding the impact of vector control on mosquito populations is crucial for planning new interventions and evaluating existing ones. However, estimates of population size changes in response to control efforts are often inaccurate because of limitations and biases in collection methods. Attempts to evaluate the impact of vector control on mosquito effective population size (N(e)) have produced inconclusive results thus far. Therefore, we obtained data for 13-15 microsatellite markers for more than 1,500 mosquitoes representing multiple time points for seven populations of three important vector species--Anopheles gambiae, An. melas, and An. moucheti--in Equatorial Guinea. These populations were exposed to indoor residual spraying or long-lasting insecticidal nets in recent years. For comparison, we also analyzed data from two populations that have no history of organized vector control. We used Approximate Bayesian Computation to reconstruct their demographic history, allowing us to evaluate the impact of these interventions on the effective population size. In six of the seven study populations, vector control had a dramatic impact on the effective population size, reducing N(e) between 55%-87%, the exception being a single An. melas population. In contrast, the two negative control populations did not experience a reduction in effective population size. This study is the first to conclusively link anti-vector intervention programs in Africa to sharply reduced effective population sizes of malaria vectors.

A metagenomic analysis of the phase 2 Anopheles gambiae 1000 genomes dataset reveals a wide diversity of cobionts associated with field collected mosquitoes.

The Anopheles gambiae 1000 Genomes (Ag1000G) Consortium previously utilized deep sequencing methods to catalogue genetic diversity across African An. gambiae populations. We analyzed the complete datasets of 1142 individually sequenced mosquitoes through Microsoft Premonition's Bayesian mixture model based (BMM) metagenomics pipeline. All specimens were confirmed as either An. gambiae sensu stricto (s.s.) or An. coluzzii with a high degree of confidence ( > 98% identity to reference). Homo sapiens DNA was identified in all specimens indicating contamination may have occurred either at the time of specimen collection, preparation and/or sequencing. We found evidence of vertebrate hosts in 162 specimens. 59 specimens contained validated Plasmodium falciparum reads. Human hepatitis B and primate erythroparvovirus-1 viral sequences were identified in fifteen and three mosquito specimens, respectively. 478 of the 1,142 specimens were found to contain bacterial reads and bacteriophage-related contigs were detected in 27 specimens. This analysis demonstrates the capacity of metagenomic approaches to elucidate important vector-host-pathogen interactions of epidemiological significance.

Five years of malaria control in the continental region, Equatorial Guinea.

BACKGROUND: A successful malaria control programme began in 2004 on Bioko Island, Equatorial Guinea. From 2007, the same multiple malaria interventions, though reduced in scope for funding reasons, were introduced to the four mainland provinces of Equatorial Guinea (the continental region) aiming to recreate Bioko's success. Two provinces received long-lasting insecticidal nets (LLINs) and two provinces received biannual indoor residual spraying (IRS). Enhanced case management and communications were introduced throughout. METHODS: Estimates of intervention coverage and indicators of malaria transmission for 2007 to 2011 were derived from annual malaria indicator surveys (MIS). Results were complemented by health information system (HIS) and entomological data. The personal protection offered by LLINs and IRS against Plasmodium falciparum infection was estimated with logistic regression. RESULTS: The estimated proportion of children aged 1-4 using either an LLIN the previous night or living in a house sprayed in the last six months was 23% in 2007 and 42% in 2011. The estimated prevalence of P. falciparum in children aged 1-4 was 68% (N=1,770; 95% confidence interval [CI]: 58-76%) in 2007 and 52% (N=1,602; 95% CI: 44-61%) in 2011. Children 1-4 years had lower prevalence if they used an LLIN the previous night (N=1,124, 56%; adjusted odds ratio [aOR] 0.64, 95% CI: 0.55-0.74) or if they lived in a sprayed house (N=1,150, 57%; aOR 0.80, 95% CI: 0.62-1.03) compared to children with neither intervention (N=4,131, 66%, reference group). The minority of children who both used an LLIN and lived in a sprayed house had the lowest prevalence of infection (N=171, 45%; aOR 0.52, 95% CI: 0.35-0.78). High site-level intervention coverage did not always correlate with lower site-level P. falciparum prevalence. The malaria season peaked in either June or July, not necessarily coinciding with MIS data collection. CONCLUSIONS: Though moderate impact was achieved after five years of vector control, case management, and communications, prevalence remained high due to an inability to sufficiently scale-up coverage with either IRS or LLINs. Both LLINs and IRS provided individual protection, but greater protection was afforded to children who benefitted from both.

Genetic isolation within the malaria mosquito Anopheles melas.

Anopheles melas is a brackish water-breeding member of the Anopheles gambiae complex that is distributed along the coast of West Africa and is a major malaria vector within its range. Because little is known about the population structure of this species, we analysed 15 microsatellite markers and 1161 bp of mtDNA in 11 A. melas populations collected throughout its range. Compared with its sibling species A. gambiae, A. melas populations have a high level of genetic differentiation between them, representing its patchy distribution due to its fragmented larval habitat that is associated with mangroves and salt marsh grass. Populations clustered into three distinct groups representing Western Africa, Southern Africa and Bioko Island populations that appear to be mostly isolated. Fixed differences in the mtDNA are present between all three clusters, and a Bayesian clustering analysis of the microsatellite data found no evidence for migration from mainland to Bioko Island populations, and little migration was evident between the Southern to the Western cluster. Surprisingly, mtDNA divergence between the three A. melas clusters is on par with levels of divergence between other species of the A. gambiae complex, and no support for monophyly was observed in a maximum-likelihood phylogenetic analysis. Finally, an approximate Bayesian analysis of microsatellite data indicates that Bioko Island A. melas populations were connected to the mainland populations in the past, but became isolated, presumably when sea levels rose after the last glaciation period (≥10 000-11 000 bp). This study has exposed species-level genetic divergence within A. melas and also has implications for control of this malaria vector.

Light traps fail to estimate reliable malaria mosquito biting rates on Bioko Island, Equatorial Guinea.

BACKGROUND: The human biting rate (HBR), an important parameter for assessing malaria transmission and evaluating vector control interventions, is commonly estimated by human landing collections (HLC). Although intense efforts have been made to find alternative non-exposure mosquito collection methods, HLC remains the standard for providing reliable and consistent HBRs. The aim of this study was to assess the relationship between human landing and light trap collections (LTC), in an attempt to estimate operationally feasible conversion factors between the two. The study was conducted as part of the operational research component of the Bioko Island Malaria Control Project (BIMCP), Equatorial Guinea. METHODS: Malaria mosquitoes were collected indoors and outdoors by HLCs and LTCs in three villages on Bioko Island, Equatorial Guinea during five bimonthly collections in 2009. Indoor light traps were suspended adjacent to occupied long-lasting, insecticide-treated bed nets. Outdoor light traps were placed close to the outer wall under the roof of the collection house. Collected specimens were subjected to DNA extraction and diagnostic PCR to identify species within the Anopheles gambiae complex. Data were analysed by simple regression of log-transformed values and by Bayesian regression analysis. RESULTS: There was a poor correlation between the two collection methods. Results varied by location, venue, month, house, but also by the statistical method used. The more robust Bayesian analyses indicated non-linear relationships and relative sampling efficiencies being density dependent for the indoor collections, implying that straight-forward and simple conversion factors could not be calculated for any of the locations. Outdoor LTC:HLC relationships were weak, but could be estimated at 0.10 and 0.07 for each of two locations. CONCLUSIONS: Light trap collections in combination with bed nets are not recommended as a reliable method to assess human biting rates on Bioko Island. Different statistical analyses methods give variable and inconsistent results. Substantial variation in collection methods prevents the determination of reliable and operationally feasible conversion factors for both indoor and outdoor data. Until improved mosquito collection methods are developed that can provide reliable and unbiased HBR estimates, HLCs should continue to serve as the reference method for HBR estimation.

Limited usefulness of microsatellite markers from the malaria vector Anopheles gambiae when applied to the closely related species Anopheles melas.

Anopheles melas is a brackish water mosquito found in coastal West Africa where it is a dominant malaria vector locally. In order to facilitate genetic studies of this species, 45 microsatellite loci originally developed for Anopheles gambiae were sequenced in An. melas. Those that were suitable based on repeat number and flanking regions were examined in 2 natural populations from Equatorial Guinea. Only 15 loci were eventually deemed suitable as polymorphic markers in An. melas populations. These loci were screened in 4 populations from a wider geographic range. Heterozygosity estimates ranged from 0.18 to 0.79, and 2.5-15 average alleles were observed per locus, yielding 13 highly polymorphic markers and 2 loci with lower variability. To examine the usefulness of microsatellite markers when applied in a sibling species, the original An. gambiae specific markers were used to amplify 5 loci in An. melas. Null alleles were found for 1 An. gambiae marker. We discuss the pitfalls of using microsatellite loci across closely related species and conclude that in addition to the problem of null alleles associated with this practice, many loci may prove to be of very limited use as polymorphic markers even when used in a sibling species.

Outdoor host seeking behaviour of Anopheles gambiae mosquitoes following initiation of malaria vector control on Bioko Island, Equatorial Guinea.

BACKGROUND: Indoor-based anti-vector interventions remain the preferred means of reducing risk of malaria transmission in malaria endemic areas around the world. Despite demonstrated success in reducing human-mosquito interactions, these methods are effective solely against endophilic vectors. It may be that outdoor locations serve as an important venue of host seeking by Anopheles gambiae sensu lato (s.l.) mosquitoes where indoor vector suppression measures are employed. This paper describes the host seeking activity of anopheline mosquito vectors in the Punta Europa region of Bioko Island, Equatorial Guinea. In this area, An. gambiae sensu stricto (s.s.) is the primary malaria vector. The goal of the paper is to evaluate the importance of An gambiae s.l. outdoor host seeking behaviour and discuss its implications for anti-vector interventions. METHODS: The venue and temporal characteristics of host seeking by anopheline vectors in a hyperendemic setting was evaluated using human landing collections conducted inside and outside homes in three villages during both the wet and dry seasons in 2007 and 2008. Additionally, five bi-monthly human landing collections were conducted throughout 2009. Collections were segregated hourly to provide a time distribution of host-seeking behaviour. RESULTS: Surprisingly high levels of outdoor biting by An. gambiae senso stricto and An. melas vectors were observed throughout the night, including during the early evening and morning hours when human hosts are often outdoors. As reported previously, An. gambiae s.s. is the primary malaria vector in the Punta Europa region, where it seeks hosts outdoors at least as much as it does indoors. Further, approximately 40% of An. gambiae s.l. are feeding at times when people are often outdoors, where they are not protected by IRS or LLINs. Repeated sampling over two consecutive dry-wet season cycles indicates that this result is independent of seasonality. CONCLUSIONS: An. gambiae s.l. mosquitoes currently seek hosts in outdoor venues as much as indoors in the Punta Europa region of Bioko Island. This contrasts with an earlier pre-intervention observation of exclusive endophagy of An. gambiae in this region. In light of this finding, it is proposed that the long term indoor application of insecticides may have resulted in an adaptive shift toward outdoor host seeking in An. gambiae s.s. on Bioko Island.

Stability and detection of nucleic acid from viruses and hosts in controlled mosquito blood feeds.

Monitoring the presence and spread of pathogens in the environment is of critical importance. Rapid detection of infectious disease outbreaks and prediction of their spread can facilitate early responses of health agencies and reduce the severity of outbreaks. Current sampling methods are sorely limited by available personnel and throughput. For instance, xenosurveillance utilizes captured arthropod vectors, such as mosquitoes, as sampling tools to access blood from a wide variety of vertebrate hosts. Next generation sequencing (NGS) of nucleic acid from individual blooded mosquitoes can be used to identify mosquito and host species, and microorganisms including pathogens circulating within either host. However, there are practical challenges to collecting and processing mosquitoes for xenosurveillance, such as the rapid metabolization or decay of microorganisms within the mosquito midgut. This particularly affects pathogens that do not replicate in mosquitoes, preventing their detection by NGS or other methods. Accordingly, we performed a series of experiments to establish the windows of detection for DNA or RNA from human blood and/or viruses present in mosquito blood meals. Our results will contribute to the development of xenosurveillance techniques with respect to optimal timing of sample collection and NGS processing and will also aid trap design by demonstrating the stabilizing effect of temperature control on viral genome recovery from blood-fed mosquitoes.

Early evening questing and oviposition activity by the Culex (Diptera: Culicidae) vectors of West Nile virus in northeastern North America.

To determine whether the Culex (Diptera: Culicidae) mosquitoes that transmit West Nile virus (family Flaviviridae, genus Flavivirus, WNV) in the northeastern United States seek hosts and oviposit contemporaneously, we recorded when these mosquitoes attacked caged birds and when they deposited eggs. They traversed oviposition sites most frequently approximately 2 h after astronomical sunset, and eggs generally were deposited at that time. Although they most frequently approached avian hosts approximately 2 h after sunset during midsummer, they are more opportunistic during mid- to late fall. Because the Culex mosquitoes that serve as the main vectors of West Nile virus in the northeastern United States quest for hosts and seek to oviposit well after sunset, insecticidal aerosols would be most effective when applied at that time.

Efficacy of resmethrin aerosols applied from the road for suppressing Culex vectors of West Nile virus.

We determined whether aerosol applications of resmethrin, delivered from the road, suppress the reproductive activity of Culex pipiens pipiens and Cx. restuans mosquitoes in suburban sites located near Boston. Oviposition implies a prior blood-feeding event and hence a potential West Nile virus (WNV) transmission-related event. Droplet size, rate of delivery and meteorological conditions were monitored. The target populations proved to be fully susceptible to the insecticide that was used. The roads in the test sites generally gave adequate opportunity for insecticidal coverage. We found that the aerosol plume may have failed to contact the target mosquitoes and conclude that such insecticidal aerosols, delivered from the road, may not effectively reduce the force of transmission of WNV in our test sites.

Insecticide resistance allele frequencies in Anopheles gambiae before and after anti-vector interventions in continental Equatorial Guinea.

Anti-malaria interventions that rely on insecticides can be compromised by insecticide-resistance alleles among malaria vectors. We examined frequency changes of resistance alleles at two loci, knockdown resistance (kdr) and acetylcholinesterase-1 (ace-1), which confer resistance to pyrethroids and DDT, and carbamates, respectively. A total of 7,059 Anopheles gambiae sensu stricto mosquitoes were analyzed from multiple sites across continental Equatorial Guinea. A subset of sites included samples collected pre-intervention (2007) and post-intervention (2009-2011). Both L1014S and L1014F resistance alleles were observed in almost all pre-intervention collections. In particular, L1014F was already at substantial frequencies in M form populations (17.6-74.6%), and at high frequencies (> 50%) in all but two S form populations. Comparison before and throughout anti-vector interventions showed drastic increases in L1014F, presumably caused by intensified selection pressure imposed by pyrethroids used in vector control efforts. In light of these findings, inclusion of other insecticide classes in any anti-vector intervention can be considered prudent.

Malaria transmission after five years of vector control on Bioko Island, Equatorial Guinea.

BACKGROUND: Malaria is endemic with year-round transmission on Bioko Island. The Bioko Island Malaria Control Project (BIMCP) started in 2004 with the aim to reduce malaria transmission and to ultimately eliminate malaria. While the project has been successful in reducing overall malaria morbidity and mortality, foci of high malaria transmission still persist on the island. Results from the 2009 entomological collections are reported here. METHODS: Human landing collections (HLC) and light trap collections (LTC) were carried out on Bioko Island, Equatorial Guinea in 2009. The HLCs were performed in three locations every second month and LTCs were carried out in 10 locations every second week. Molecular analyses were performed to identify species, detect sporozoites, and identify potential insecticide resistance alleles. RESULTS: The entomological inoculation rates (EIR) on Bioko Island ranged from 163 to 840, with the outdoor EIRs reaching > 900 infective mosquito bites per year. All three human landing collection sites on Bioko Island had an annual EIR exceeding the calculated African average of 121 infective bites per year. The highest recorded EIRs were in Punta Europa in northwestern Bioko Island with human biting rates of 92 and 66 mosquito landings per person per night, outdoors and indoors, respectively. Overall, the propensity for mosquito biting on the island was significantly higher outdoors than indoors (p < 0.001). Both Anopheles gambiae s.s. and An. melas were responsible for malaria transmission on the island, but with different geographical distribution patterns. Sporozoite rates were the highest in An. gambiae s.s. populations ranging from 3.1% in Punta Europa and 5.7% in Riaba in the southeast. Only the L1014F (kdr-west) insecticide resistance mutation was detected on the island with frequencies ranging from 22-88% in An. gambiae s.s. No insecticide resistance alleles were detected in the An. melas populations. CONCLUSIONS: In spite of five years of extensive malaria control and a generalized reduction in the force of transmission, parasite prevalence and child mortality, foci of very high transmission persist on Bioko Island, particularly in the northwestern Punta Europa area. This area is favorable for anopheline mosquito breeding; human biting rates are high, and the EIRs are among the highest ever recorded. Both vector species collected in the study have a propensity to bite outdoors more frequently than indoors. Despite current vector control efforts mosquito densities remain high in such foci of high malaria transmission. To further reduce transmission, indoor residual spraying (IRS) needs to be supplemented with additional vector control interventions.

Spatial and temporal variation in vector competence of Culex pipiens and Cx. restuans mosquitoes for West Nile virus.

Vector competence, the probability that a vector will transmit a pathogen after feeding on an infected host, is known to vary among vector species, populations, days since feeding, and temperature during the extrinsic incubation period. However, the extent of spatio-temporal variability and consistency in vector competence of populations is not known. We examined vector competence of Culex pipiens Linnaeus and Cx. restuans Theobald mosquitoes for West Nile virus collected over 3 years from 17 sites to measure spatial and temporal scales of variation in vector competence. We found extreme variation with 0-52% of mosquitoes transmitting West Nile virus at a single site between different sampling periods, and similar variation across populations. However, we also found that within a smaller geographic range, vector competence varied somewhat synchronously, suggesting that environmental and population genetic factors might influence vector competence. These results highlight the spatio-temporal variability in vector competence and the role of local processes.