Gambiae Species Complex Research Articles

Variation in spatial population structure in the Anopheles gambiae species complex.

Anopheles gambiae, Anopheles coluzzii, and Anopheles arabiensis are three of the most widespread vectors of malaria parasites, with geographical ranges stretching across wide swaths of Africa. Understanding the population structure of these closely related species, including the extent to which populations are connected by gene flow, is essential for understanding how vector control implemented in one location might indirectly affect vector populations in other locations. Here, we assessed the population structure of each species based on a combined data set of publicly available and newly processed whole-genome sequences. The data set included single nucleotide polymorphisms from whole genomes of 2,410 individual mosquitoes sampled from 128 locations across 19 African countries. We found that A. gambiae sampled from several countries in West and Central Africa showed low genetic differentiation from each other according to principal components analysis (PCA) and ADMIXTURE modeling. Using Estimated Effective Migration Surfaces (EEMS), we showed that this low genetic differentiation indicates high effective migration rates for A. gambiae across this region. Outside of this region, we found eight groups of sampling locations from Central, East, and Southern Africa for which A. gambiae showed higher genetic differentiation, and lower effective migration rates, between each other and the West/Central Africa group. These results indicate that the barriers to and corridors for migration between populations of A. gambiae differ across the geographical range of this malaria vector species. Using the same methods, we found higher genetic differentiation and lower migration rates between populations of A. coluzzii in West and Central Africa than for A. gambiae in the same region. In contrast, we found lower genetic differentiation and higher migration rates between populations of A. arabiensis in Tanzania, compared to A. gambiae in the same region. These differences between A. gambiae, A. coluzzii, and A. arabiensis indicate that migration barriers and corridors may vary, even between very closely related species. Overall, our results demonstrate that migration rates vary both within and between species of Anopheles mosquitoes, presumably based on species-specific responses to the ecological or environmental conditions that may impede or facilitate migration, and the geographical patterns of these conditions across the landscape. Together with previous findings, this study provides robust evidence that migration rates between populations of malaria vectors depend on the ecological context, which should be considered when planning surveillance of vector populations, monitoring for insecticide resistance, and evaluating interventions.

Measuring the Impact of Genetic Heterogeneity and Chromosomal Inversions on the Efficacy of CRISPR-Cas9 Gene Drives in Different Strains of Anopheles gambiae.

The human malaria vector Anopheles gambiae is becoming increasingly resistant to insecticides, spurring the development of genetic control strategies. CRISPR-Cas9 gene drives can modify a population by creating double-stranded breaks at highly specific targets, triggering copying of the gene drive into the cut site ("homing"), ensuring its inheritance. The DNA repair mechanism responsible requires homology between the donor and recipient chromosomes, presenting challenges for the invasion of laboratory-developed gene drives into wild populations of target species An. gambiae species complex, which show high levels of genome variation. Two gene drives (vas2-5958 and zpg-7280) were introduced into three An. gambiae strains collected across Africa with 5.3-6.6% variation around the target sites, and the effect of this variation on homing was measured. Gene drive homing across different karyotypes of the 2La chromosomal inversion was also assessed. No decrease in gene drive homing was seen despite target site heterology, demonstrating the applicability of gene drives to wild populations.

Wolbachia endosymbionts in two Anopheles species indicates independent acquisitions and lack of prophage elements

Wolbachia is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species, Wolbachia pipientis , divided into several ‘supergroups’ based on multilocus sequence typing. Wolbachia strains in mosquitoes have been shown to inhibit the transmission of human pathogens, including Plasmodium malaria parasites and arboviruses. Despite their large host range, Wolbachia strains within the major malaria vectors of the Anopheles gambiae and Anopheles funestus complexes appear at low density, established solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B Wolbachia within Anopheles demeilloni and Anopheles moucheti has opened exciting possibilities to explore naturally occurring Wolbachia endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here, we utilize genomic analyses to demonstrate that both Wolbachia strains have retained all key metabolic and transport pathways despite their smaller genome size, with this reduction potentially attributable to degenerated prophage regions. Even with this reduction, we confirmed the presence of cytoplasmic incompatibility (CI) factor genes within both strains, with wAnD maintaining intact copies of these genes while the cifB gene was interrupted in wAnM, so functional analysis is required to determine whether wAnM can induce CI. Additionally, phylogenetic analysis indicates that these Wolbachia strains may have been introduced into these two Anopheles species via horizontal transmission events, rather than by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first Wolbachia genomes, to our knowledge, that enable us to study the relationship between natural strain Plasmodium malaria parasites and their anopheline hosts.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii.

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences. Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Characterisation of the epigenetic architecture of the major malaria vector Anopheles arabiensis (Diptera: Culicidae) after treatment with epigenetic modulators and heavy metals

Anopheles arabiensis (a member of the An. gambiae species complex) is a major vector of malaria in sub-Saharan Africa. Despite its disease vector status, there is currently a paucity of epigenetic information for this species. The aim this study was therefore to analyse global epigenetic markers and their response to metal exposure in insecticide susceptible and resistant laboratory strains of An. arabiensis. This was done using commercially available epigenetic marker quantification kits. In order to validate the efficacy of the kits, several kits were assessed to determine whether changes induced by known epigenetic modulators were detectable using these platforms. The efficacy of the dosages used were determined by examining the effect of the dosages used on insecticide resistant phenotypes. Upon confirmation that the dosages used were sufficient to induce a phenotypic change, the effect on epigenetic markers was assessed. Commercial kits were used to quantify 5-methylcysteine (5-mC) and 5-hydroxymethylcysteine (5-hmC) methylation in DNA, m6A methylation in mRNA as well as Histone Acetyl Transferase (HAT) activity. There was a marked difference in the phenotypic response in adult mosquitoes of the insecticide susceptible strain compared to that of its’ resistant counterpart. For males and females of the resistant strain, exposure to nucleic acid modifying drugs typically increased their tolerance to insecticides. The patterns of changes in 5-mC methylation by epigenetic modulators was congruent with previous studies which quantified by mass spectrometry. The two strains differed in methylation patterns under control conditions and responded differentially to larval metal exposure. In the resistant strain, which previously was demonstrated to show increased detoxification enzyme activity and insecticide tolerance after the same treatment, the potential increase in transcriptional activity appeared to be modulated by reduced methylation and increased HAT activity. This study suggests that the commercial epigenetic quantification kits can be used to characterise phenotypic changes in An. arabiensis, and also shows that epigenetic regulation of the response to metal exposure is regulated at the DNA as opposed to the RNA level.

Horizontal Transmission of the Symbiont Microsporidia MB in Anopheles arabiensis.

The recently discovered Anopheles symbiont, Microsporidia MB, has a strong malaria transmission-blocking phenotype in Anopheles arabiensis, the predominant Anopheles gambiae species complex member in many active transmission areas in eastern Africa. The ability of Microsporidia MB to block Plasmodium transmission together with vertical transmission and avirulence makes it a candidate for the development of a symbiont-based malaria transmission blocking strategy. We investigate the characteristics and efficiencies of Microsporidia MB transmission between An. arabiensis mosquitoes. We show that Microsporidia MB is not transmitted between larvae but is effectively transmitted horizontally between adult mosquitoes. Notably, Microsporidia MB was only found to be transmitted between male and female An. arabiensis, suggesting sexual horizontal transmission. In addition, Microsporidia MB cells were observed infecting the An. arabiensis ejaculatory duct. Female An. arabiensis that acquire Microsporidia MB horizontally are able to transmit the symbiont vertically to their offspring. We also investigate the possibility that Microsporidia MB can infect alternate hosts that live in the same habitats as their An. arabiensis hosts, but find no other non-anopheline hosts. Notably, Microsporidia MB infections were found in another primary malaria African vector, Anopheles funestus s.s. The finding that Microsporidia MB can be transmitted horizontally is relevant for the development of dissemination strategies to control malaria that are based on the targeted release of Microsporidia MB infected Anopheles mosquitoes.

The Genetic Architecture of Post-Zygotic Reproductive Isolation Between Anopheles coluzzii and An. quadriannulatus.

The Anopheles gambiae complex is comprised of eight morphologically indistinguishable species and has emerged as a model system for the study of speciation genetics due to the rapid radiation of its member species over the past two million years. Male hybrids between most An. gambiae complex species pairs are sterile, and some genotype combinations in hybrid males cause inviability. We investigated the genetic basis of hybrid male inviability and sterility between An. coluzzii and An. quadriannulatus by measuring segregation distortion and performing a QTL analysis of sterility in a backcross population. Hybrid males were inviable if they inherited the An. coluzzii X chromosome and were homozygous at one or more loci in 18.9 Mb region of chromosome 3. The An. coluzzii X chromosome has a disproportionately large effect on hybrid sterility when introgressed into an An. quadriannulatus genetic background. Additionally, an epistatic interaction between the An. coluzzii X and a 1.12 Mb, pericentric region of the An. quadriannulatus 3L chromosome arm has a statistically significant contribution to the hybrid sterility phenotype. This same epistatic interaction occurs when the An. coluzzii X is introgressed into the genetic background of An. arabiensis, the sister species of An. quadriannulatus, suggesting that this may represent one of the first Dobzhansky–Muller incompatibilities to evolve early in the radiation of the Anopheles gambiae species complex. We describe the additive effects of each sterility QTL, epistatic interactions between them, and genes within QTL with protein functions related to mating behavior, reproduction, spermatogenesis, and microtubule morphogenesis, whose divergence may contribute to post-zygotic reproductive isolation between An. coluzzii and An. quadriannulatus.

The Impact of Cross-Species Gene Flow on Species Tree Estimation.

Recent analyses of genomic sequence data suggest cross-species gene flow is common in both plants and animals, posing challenges to species tree estimation. We examine the levels of gene flow needed to mislead species tree estimation with three species and either episodic introgressive hybridization or continuous migration between an outgroup and one ingroup species. Several species tree estimation methods are examined, including the majority-vote method based on the most common gene tree topology (with either the true or reconstructed gene trees used), the UPGMA method based on the average sequence distances (or average coalescent times) between species, and the full-likelihood method based on multilocus sequence data. Our results suggest that the majority-vote method based on gene tree topologies is more robust to gene flow than the UPGMA method based on coalescent times and both are more robust than likelihood assuming a multispecies coalescent (MSC) model with no cross-species gene flow. Comparison of the continuous migration model with the episodic introgression model suggests that a small amount of gene flow per generation can cause drastic changes to the genetic history of the species and mislead species tree methods, especially if the species diverged through radiative speciation events. Estimates of parameters under the MSC with gene flow suggest that African mosquito species in the Anopheles gambiae species complex constitute such an example of extreme impact of gene flow on species phylogeny. [IM; introgression; migration; MSci; multispecies coalescent; species tree.].

Gene copy number and function of the APL1 immune factor changed during Anopheles evolution

BackgroundThe recent reference genome assembly and annotation of the Asian malaria vector Anopheles stephensi detected only one gene encoding the leucine-rich repeat immune factor APL1, while in the Anopheles gambiae and sibling Anopheles coluzzii, APL1 factors are encoded by a family of three paralogs. The phylogeny and biological function of the unique APL1 gene in An. stephensi have not yet been specifically examined.MethodsThe APL1 locus was manually annotated to confirm the computationally predicted single APL1 gene in An. stephensi. APL1 evolution within Anopheles was explored by phylogenomic analysis. The single or paralogous APL1 genes were silenced in An. stephensi and An. coluzzii, respectively, followed by mosquito survival analysis, experimental infection with Plasmodium and expression analysis.ResultsAPL1 is present as a single ancestral gene in most Anopheles including An. stephensi but has expanded to three paralogs in an African lineage that includes only the Anopheles gambiae species complex and Anopheles christyi. Silencing of the unique APL1 copy in An. stephensi results in significant mosquito mortality. Elevated mortality of APL1-depleted An. stephensi is rescued by antibiotic treatment, suggesting that pathology due to bacteria is the cause of mortality, and indicating that the unique APL1 gene is essential for host survival. Successful Plasmodium development in An. stephensi depends upon APL1 activity for protection from high host mortality due to bacteria. In contrast, silencing of all three APL1 paralogs in An. coluzzii does not result in elevated mortality, either with or without Plasmodium infection. Expression of the single An. stephensi APL1 gene is regulated by both the Imd and Toll immune pathways, while the two signaling pathways regulate different APL1 paralogs in the expanded APL1 locus.ConclusionsAPL1 underwent loss and gain of functions concomitant with expansion from a single ancestral gene to three paralogs in one lineage of African Anopheles. We infer that activity of the unique APL1 gene promotes longevity in An. stephensi by conferring protection from or tolerance to an effect of bacterial pathology. The evolution of an expanded APL1 gene family could be a factor contributing to the exceptional levels of malaria transmission mediated by human-feeding members of the An. gambiae species complex in Africa.

Fine-Mapping Complex Inversion Breakpoints and Investigating Somatic Pairing in the Anopheles gambiae Species Complex Using Proximity-Ligation Sequencing.

Chromosomal inversions are fundamental drivers of genome evolution. In the main Afrotropical malaria vector species, belonging to the Anopheles gambiae species complex, inversions play an important role in local adaptation and have a rich history of cytological study. Despite the importance and ubiquity of some chromosomal inversions across the species complex, inversion breakpoints are often challenging to map molecularly due to the presence of large repetitive regions. Here, we develop an approach that uses Hi-C sequencing data to molecularly fine-map the breakpoints of inversions. We demonstrate that this approach is robust and likely to be widely applicable for both identification and fine-mapping inversion breakpoints in species whose inversions have heretofore been challenging to characterize. We apply our method to interrogate the previously unknown inversion breakpoints of 2Rbc and 2Rd in An. coluzzii We found that inversion breakpoints occur in large repetitive regions, and, strikingly, among three inversions analyzed, two breakpoints appear to be reused in two separate inversions. These breakpoint-adjacent regions are strongly enriched for the presence of a 30 bp satellite repeat sequence. Because low frequency inversion breakpoints are not correlated with genomic regions containing this satellite, we suggest that interrupting this particular repeat may result in arrangements with higher relative fitness. Additionally, we use heterozygous individuals to quantitatively investigate the impacts of somatic pairing in the regions immediately surrounding inversion breakpoints. Finally, we discuss important considerations for possible applications of this approach for inversion breakpoint identification in a range of organisms.

Analysis-ready datasets for insecticide resistance phenotype and genotype frequency in African malaria vectors

The impact of insecticide resistance in malaria vectors is poorly understood and quantified. Here a series of geospatial datasets for insecticide resistance in malaria vectors are provided, so that trends in resistance in time and space can be quantified, and the impact of resistance found in wild populations on malaria transmission in Africa can be assessed. Specifically, data have been collated and geopositioned for the prevalence of insecticide resistance, as measured by standard bioassays, in representative samples of individual species or species complexes. Data are provided for the Anopheles gambiae species complex, the Anopheles funestus subgroup, and for nine individual vector species. Data are also given for common genetic markers of resistance to support analyses of whether these markers can improve the ability to monitor resistance in low resource settings. Allele frequencies for known resistance-associated markers in the Voltage-gated sodium channel (Vgsc) are provided. In total, eight analysis-ready, standardised, geopositioned datasets encompassing over 20,000 African mosquito collections between 1957 and 2017 are released.

Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex.

Deep coalescence and introgression make it challenging to infer phylogenetic relationships among closely related species that arose through radiative speciation events. Despite numerous phylogenetic analyses and the availability of whole genomes, the phylogeny in the Anopheles gambiae species complex has not been confidently resolved. Here we extract over 80, 000 coding and noncoding short segments (called loci) from the genomes of six members of the species complex and use a Bayesian method under the multispecies coalescent model to infer the species tree, which takes into account genealogical heterogeneity across the genome and uncertainty in the gene trees. We obtained a robust estimate of the species tree from the distal region of the X chromosome: (A. merus, ((A. melas, (A. arabiensis, A. quadriannulatus)), (A. gambiae, A. coluzzii))), with A. merus to be the earliest branching species. This species tree agrees with the chromosome inversion phylogeny and provides a parsimonious interpretation of inversion and introgression events. Simulation informed by the real data suggest that the coalescent approach is reliable while the sliding-window analysis used in a previous phylogenomic study generates artifactual species trees. Likelihood ratio test of gene flow revealed strong evidence of autosomal introgression from A. arabiensis into A. gambiae (at the average rate of migrants per generation), but not in the opposite direction, and introgression of the 3 L chromosomal region from A. merus into A. quadriannulatus. Our results highlight the importance of accommodating incomplete lineage sorting and introgression in phylogenomic analyses of species that arose through recent radiative speciation events.

The Effect of Hybridization on Dosage Compensation in Member Species of the Anopheles gambiae Species Complex.

Dosage compensation has evolved in concert with Y-chromosome degeneration in many taxa that exhibit heterogametic sex chromosomes. Dosage compensation overcomes the biological challenge of a “half dose” of X chromosome gene transcripts in the heterogametic sex. The need to equalize gene expression of a hemizygous X with that of autosomes arises from the fact that the X chromosomes retain hundreds of functional genes that are actively transcribed in both sexes and interact with genes expressed on the autosomes. Sex determination and heterogametic sex chromosomes have evolved multiple times in Diptera, and in each case the genetic control of dosage compensation is tightly linked to sex determination. In the Anopheles gambiae species complex (Culicidae), maleness is conferred by the Y-chromosome gene Yob, which despite its conserved role between species is polymorphic in its copy number between them. Previous work demonstrated that male An. gambiae s.s. males exhibit complete dosage compensation in pupal and adult stages. In the present study, we have extended this analysis to three sister species in the An. gambiae complex: An. coluzzii, An. arabiensis, and An. quadriannulatus. In addition, we analyzed dosage compensation in bi-directional F1 hybrids between these species to determine if hybridization results in the mis-regulation and disruption of dosage compensation. Our results confirm that dosage compensation operates in the An. gambiae species complex through the hypertranscription of the male X chromosome. Additionally, dosage compensation in hybrid males does not differ from parental males, indicating that hybridization does not result in the mis-regulation of dosage compensation.

A selective αvβ5 integrin antagonist hidden into the anophelin family protein cE5 from the malaria vector Anopheles gambiae

AbstractA RGD motif was identified in the N‐terminal region of cE5, a potent salivary thrombin inhibitor from the African malaria vector Anopheles gambiae. A peptide (APQ30) encompassing the first 30 amino acids residues of the protein and including the RGD tripeptide was tested in cell adhesion assays and found to inhibit αvβ3 and αvβ5 mediated adhesion. A shorter peptide (APQ16), strongly conserved among members of the A. gambiae species complex and including only the first 16 residues, retained adhesion inhibitory properties, however with enhanced specificity toward αvβ5. In addition, migration and invasion assays showed its capacity to inhibit the invasiveness of the malignant cell lines HepG2 and MDA‐MB231. Altogether our data point to APQ16 as a new promising candidate as theranostic agent.

The Anopheles gambiae 2La chromosome inversion is associated with susceptibility to Plasmodium falciparum in Africa.

Chromosome inversions suppress genetic recombination and establish co-adapted gene complexes, or supergenes. The 2La inversion is a widespread polymorphism in the Anopheles gambiae species complex, the major African mosquito vectors of human malaria. Here we show that alleles of the 2La inversion are associated with natural malaria infection levels in wild-captured vectors from West and East Africa. Mosquitoes carrying the more-susceptible allele (2L+a) are also behaviorally less likely to be found inside houses. Vector control tools that target indoor-resting mosquitoes, such as bednets and insecticides, are currently the cornerstone of malaria control in Africa. Populations with high levels of the 2L+a allele may form reservoirs of persistent outdoor malaria transmission requiring novel measures for surveillance and control. The 2La inversion is a major and previously unappreciated component of the natural malaria transmission system in Africa, influencing both malaria susceptibility and vector behavior.

Preferential suppression of Anopheles gambiae host sequences allows detection of the mosquito eukaryotic microbiome

Anopheles mosquitoes are vectors of the human malaria parasite, Plasmodium falciparum. The vector microbiota is a likely factor influencing parasite transmission. The prokaryotic microbiota of mosquitoes is efficiently surveyed by sequencing of hypervariable regions of the 16s ribosomal RNA (rRNA) gene. However, identification of the eukaryotic microbiota by targeting the 18s rRNA gene is challenging due to simultaneous amplification of the abundant 18s rRNA gene target in the mosquito host. Consequently, the eukaryotic microbial diversity of mosquitoes is vastly underexplored. An efficient methodology is needed to identify this component of the microbiota, expected to include relatives of Plasmodium. Here, we use defined panels of Anopheles samples from West Africa to test two experimental PCR clamp approaches to maximize the specific amplification of 18s rRNA gene hypervariable regions from eukaryotic microbes: anneal-inhibiting blocking primers and peptide-nucleic acid (PNA) oligonucleotide blockers. Of the two, PNA blockers were the only efficient blocking strategy, allowing a reduction of mosquito 18s rRNA gene sequences by more than 80% for the V4 hypervariable region. These PNA blockers will facilitate taxonomic profiling of the eukaryotic microbiota of the A. gambiae species complex, and contribute to a better understanding of microbial influence upon immunity and pathogen infection.

Geographical distributions of African malaria vector sibling species and evidence for insecticide resistance

BackgroundMany of the mosquito species responsible for malaria transmission belong to a sibling complex; a taxonomic group of morphologically identical, closely related species. Sibling species often differ in several important factors that have the potential to impact malaria control, including their geographical distribution, resistance to insecticides, biting and resting locations, and host preference. The aim of this study was to define the geographical distributions of dominant malaria vector sibling species in Africa so these distributions can be coupled with data on key factors such as insecticide resistance to aid more focussed, species-selective vector control.ResultsWithin the Anopheles gambiae species complex and the Anopheles funestus subgroup, predicted geographical distributions for Anopheles coluzzii, An. gambiae (as now defined) and An. funestus (distinct from the subgroup) have been produced for the first time. Improved predicted geographical distributions for Anopheles arabiensis, Anopheles melas and Anopheles merus have been generated based on records that were confirmed using molecular identification methods and a model that addresses issues of sampling bias and past changes to the environment. The data available for insecticide resistance has been evaluated and differences between sibling species are apparent although further analysis is required to elucidate trends in resistance.ConclusionsSibling species display important variability in their geographical distributions and the most important malaria vector sibling species in Africa have been mapped here for the first time. This will allow geographical occurrence data to be coupled with species-specific data on important factors for vector control including insecticide resistance. Species-specific data on insecticide resistance is available for the most important malaria vectors in Africa, namely An. arabiensis, An. coluzzii, An. gambiae and An. funestus. Future work to combine these data with the geographical distributions mapped here will allow more focussed and resource-efficient vector control and provide information to greatly improve and inform existing malaria transmission models.

Genome-Wide Divergence in the West-African Malaria Vector Anopheles melas.

Anopheles melas is a member of the recently diverged An. gambiae species complex, a model for speciation studies, and is a locally important malaria vector along the West-African coast where it breeds in brackish water. A recent population genetic study of An. melas revealed species-level genetic differentiation between three population clusters. An. melas West extends from The Gambia to the village of Tiko, Cameroon. The other mainland cluster, An. melas South, extends from the southern Cameroonian village of Ipono to Angola. Bioko Island, Equatorial Guinea An. melas populations are genetically isolated from mainland populations. To examine how genetic differentiation between these An. melas forms is distributed across their genomes, we conducted a genome-wide analysis of genetic differentiation and selection using whole genome sequencing data of pooled individuals (Pool-seq) from a representative population of each cluster. The An. melas forms exhibit high levels of genetic differentiation throughout their genomes, including the presence of numerous fixed differences between clusters. Although the level of divergence between the clusters is on a par with that of other species within the An. gambiae complex, patterns of genome-wide divergence and diversity do not provide evidence for the presence of pre- and/or postmating isolating mechanisms in the form of speciation islands. These results are consistent with an allopatric divergence process with little or no introgression.

Seeing is believing: the nocturnal malarial mosquito Anopheles coluzzii responds to visual host-cues when odour indicates a host is nearby.

BackgroundThe immediate aim of our study was to analyse the behaviour of the malarial mosquito Anopheles coluzzii (An. gambiae species complex) near a human host with the ultimate aim of contributing to our fundamental understanding of mosquito host-seeking behaviour and the overall aim of identifying behaviours that could be exploited to enhance sampling and control strategies.ResultsBased on 3D video recordings of individual host-seeking females in a laboratory wind-tunnel, we found that despite being a nocturnal species, An. coluzzii is highly responsive to a visually conspicuous object, but only in the presence of host-odour. Female mosquitoes approached and abruptly veered away from a dark object, which suggests attraction to visual cues plays a role in bringing mosquitoes to the source of host odour. It is worth noting that the majority of our recorded flight tracks consisted of highly stereotyped ‘dipping’ sequences near the ground, which have been mentioned in the literature, but never before quantified.ConclusionsOur quantitative analysis of female mosquito flight patterns within ~1.5 m of a host has revealed highly relevant information about responsiveness to visual objects and flight height that could revolutionise the efficacy of sampling traps; the capturing device of a trap should be visually conspicuous and positioned near the ground where the density of host-seeking mosquitoes would be greatest. These characteristics are not universally present in current traps for malarial mosquitoes. The characterisation of a new type of flight pattern that is prevalent in mosquitoes suggests that there is still much that is not fully understood about mosquito flight behaviour.Electronic supplementary materialThe online version of this article (doi:10.1186/s13071-016-1609-z) contains supplementary material, which is available to authorized users.

Circumsporozoite protein rates, blood-feeding pattern and frequency of knockdown resistance mutations in Anopheles spp. in two ecological zones of Mauritania.

BackgroundMosquitoes belonging to Anopheles gambiae species complex are the main malaria vector in Mauritania but data on their vector capacities, feeding habits and insecticide susceptibility are still scanty. The objectives of this study were to fill this gap.MethodsAdult Anopheles spp. mosquitoes were collected using pyrethrum spray catch method from two ecological zones of Mauritania: Nouakchott (Saharan zone) and Hodh Elgharbi region (Sahelian zone). Circumsporozoite proteins (CSP) for P. falciparum, P. vivax VK210 and P. vivax VK247 were detected by enzyme-linked immunosorbent assay (ELISA) from the female anopheline mosquitoes. To confirm CSP-ELISA results, polymerase chain reaction (PCR) was also performed. Blood meal identification was performed in all engorged females by partial sequencing of the mitochondrial cytochrome b gene. Molecular assessments of pyrethroid knockdown resistance (kdr) and insensitive acetylcholinesterase resistance (ace-1) were conducted.ResultsIn Nouakchott, the only species of Anopheles identified during the survey was Anopheles arabiensis (356 specimens). In Hodh Elgharbi, 1016 specimens of Anopheles were collected, including 578 (56.9 %) Anopheles rufipes, 410 (40.35 %) An. arabiensis, 20 (1.96 %) An. gambiae, 5 (0.5 %) An. pharoensis and 3 (0.3 %) An. funestus. Three of 186 female An. arabiensis collected in Nouakchott and tested by ELISA were found positive for Plasmodium vivax VK210, corresponding to a sporozoite rate of 1.6 %; however PCR confirmed infection by P. vivax sporozoite in only one of these. In Hodh Elgharbi, no mosquito was found positive for Plasmodium spp. infection. There was a statistically significant difference in the percentage of human blood-fed Anopheles spp. between Nouakchott (58.7 %, 47 of 80 blood-engorged An. arabiensis females) and Hodh Elgharbi (11.1 %, 2 of 18 blood-engorged mosquitoes). Analysis of the kdr polymorphisms showed 48.2 % (70/145) of East African kdr mutation (L1014S) in Nouakchott compared to 10 % (4/40) in Hodh Elgharbi region (P < 0.001). Nevertheless, West African kdr mutation (L1014F) was found only in An. gambiae populations (4/40, 10 %) from Hodh Elgharbi region. No ace-1 mutation was found in mosquito specimens from the two study zones.ConclusionsOverall, this study confirmed the autochthonous P. vivax malaria transmission in Nouakchott, involving An. arabiensis as the main vector. It also described for the first time the absence of ace-1 mutation, the co-occurrence of both West and East African kdr mutation in An. gambiae in Mauritania, and highlighted the regional variations in the prevalence and type of kdr mutations.