Anopheles gambiae on remote islands in the Indian Ocean: origins and prospects for malaria elimination by genetic modification of extant populations

The mosquito Anopheles gambiae s.s. is a primary malaria vector throughout sub-Saharan Africa including the islands of the Comoros archipelago (Anjouan, Grande Comore, Mayotte and Mohéli). These islands are located at the northern end of the Mozambique Channel in eastern Africa. Previous studies have shown a relatively high degree of genetic isolation between the Comoros islands and mainland populations of A. gambiae, but the origin of the island populations remains unclear. Here, we analyzed phylogenetic relationships among island and mainland populations using complete mitochondrial genome sequences of individual A. gambiae specimens. This work augments earlier studies based on analysis of the nuclear genome. We investigated the source population of A. gambiae for each island, estimated the number of introductions, when they occurred and explored evidence for contemporary gene flow between island and mainland populations. These studies are relevant to understanding historical patterns in the dispersal of this important malaria vector and provide information critical to assessing their potential for the exploration of genetic-based vector control methods to eliminate this disease. Phylogenetic analysis and haplotype networks were constructed from mitogenome sequences of 258 A. gambiae from the four islands. In addition, 112 individuals from seven countries across sub-Saharan Africa and Madagascar were included to identify potential source populations. Our results suggest that introduction events of A. gambiae into the Comoros archipelago were rare and recent events and support earlier claims that gene flow between the mainland and these islands is limited. This study is concordant with earlier work suggesting the suitability of these oceanic islands as appropriate sites for conducting field trial releases of genetically engineered mosquitoes (GEMs).

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

Genomic signatures of local adaptation in recent invasive Aedes aegypti populations in California

BACKGROUND

Rapid adaptation to new environments can facilitate species invasions and range expansions. Understanding the mechanisms of adaptation used by invasive disease vectors in new regions has key implications for mitigating the prevalence and spread of vector-borne disease, although they remain relatively unexplored.

RESULTS

Here, we integrate whole-genome sequencing data from 96 Aedes aegypti mosquitoes collected from various sites in southern and central California with 25 annual topo-climate variables to investigate genome-wide signals of local adaptation among populations. Patterns of population structure, as inferred using principal components and admixture analysis, were consistent with three genetic clusters. Using various landscape genomics approaches, which all remove the confounding effects of shared ancestry on correlations between genetic and environmental variation, we identified 112 genes showing strong signals of local environmental adaptation associated with one or more topo-climate factors. Some of them have known effects in climate adaptation, such as heat-shock proteins, which shows selective sweep and recent positive selection acting on these genomic regions.

CONCLUSIONS

Our results provide a genome wide perspective on the distribution of adaptive loci and lay the foundation for future work to understand how environmental adaptation in Ae. aegypti impacts the arboviral disease landscape and how such adaptation could help or hinder efforts at population control.

Population Genetics of Anopheles pretoriensis in Grande Comore Island

Anopheles pretoriensis is widely distributed across Africa, including on oceanic islands such as Grande Comore in the Comoros. This species is known to be mostly zoophylic and therefore considered to have low impact on the transmission of human malaria. However, A. pretoriensis has been found infected with Plasmodium, suggesting that it may be epidemiologically important. In the present study, we sequenced and assembled the complete mitogenome of A. pretoriensis and inferred its phylogenetic relationship among other species in the subgenus Cellia. We also investigated the genetic structure of A. pretoriensis populations on Grande Comore Island, and between this island population and sites in continental Africa, using partial sequence of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Seven haplotypes were found on the island, one of which was ubiquitous. There was no clear divergence between island haplotypes and those found on the continent. The present work contributes knowledge on this understudied, yet abundant, Anopheles species.

Mitogenome Analyses Reveal Limited Introduction of Anopheles coluzzii Into the Central African Islands of São Tomé and Príncipe

Islands possess physical characteristics that make them uniquely well-suited for initial field trials of new genetic-based technologies applied to African malaria vectors. This has led to efforts to characterize the degree of isolation of island mosquito populations. São Tomé and Príncipe (STP) is a country composed of two small islands in the Gulf of Guinea (Central Africa) where Anopheles coluzzii is the primary malaria vector. Several studies have shown a relatively high degree of genetic isolation between A. coluzzii populations in STP and the mainland compared with pairs of mainland populations separated by equivalent distances. Here, we analyzed complete mitochondrial genomes of individual A. coluzzii specimens from STP and neighboring mainland countries. The objectives are to describe the history of A. coluzzii establishment in STP, specifically to address several questions germane to their suitability as sites for a field trial release of genetically engineered mosquitoes (GEMs). These questions include: (i) What are the origins of A. coluzzii populations in STP?; (ii) How many introductions occurred?; (iii) When was A. coluzzii introduced into STP? and (iv) Is there ongoing, contemporary gene flow into STP from mainland populations? Phylogenetic analysis and haplotype networks were constructed from sequences of 345 A. coluzzii from STP, and 107 individuals from 10 countries on or near the west coast of Africa. Analysis of these data suggest that there have been two introductions of A. coluzzii onto the island of São Tomé that occurred roughly 500 years ago and that these originated from mainland West Africa. It appears that A. coluzzii has never been introduced into Príncipe Island directly from mainland Africa, but there have been at least four introductions originating from São Tomé. Our findings provide further support for the notion that contemporary populations of A. coluzzii on São Tomé and Príncipe are genetically isolated from mainland populations of this mosquito species.

Ethical Considerations for Gene Drive: Challenges of Balancing Inclusion, Power and Perspectives

Progress in gene-drive research has stimulated discussion and debate on ethical issues including community engagement and consent, policy and governance, and decision-making involved in development and deployment. Many organizations, academic institutions, foundations, and individual professionals have contributed to ensuring that these issues are considered prior to the application of gene-drive technology. Central topics include co-development of the technology with local stakeholders and communities and reducing asymmetry between developers and end-users. Important questions include with whom to conduct engagement and how to define community acceptance, develop capacity-building activities, and regulate this technology. Experts, academics, and funders have suggested that global frameworks, standards, and guidelines be developed to direct research in answering these important questions. Additionally, it has been suggested that ethical principles or commitments be established to further guide research practices. The challenging and interesting contradiction that we explore here is that the vast majority of these conversations transpire with little or no input from potential end-users or stakeholders who, we contend, should ultimately determine the fate of the technology in their communities. The question arises, whose concerns regarding marginalization, disempowerment, and inequity should be included in discussions and decisions concerning how inequities are perceived and how they may be addressed? At what stage will true co-development occur and how will opinions, perspectives and knowledge held by low-income country stakeholders be applied in determining answers to the questions regarding the ethics being debated on the academic stage? Our opinion is that the time is now.

Spontaneous mutation rate estimates for the principal malaria vectors Anopheles coluzzii and Anopheles stephensi

Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10⁻⁹ (95% confidence interval, 2.06 × 10⁻¹⁰—2.91 × 10⁻⁹) and 1.36 × 10⁻⁹ (95% confidence interval, 4.42 × 10⁻¹⁰—3.18 × 10⁻⁹) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.

Selection of sites for field trials of genetically engineered mosquitoes with gene drive

Novel malaria control strategies using genetically engineered mosquitoes (GEMs) are on the horizon. Population modification is one approach wherein mosquitoes are engineered with genes rendering them refractory to the malaria parasite, Plasmodium falciparum, coupled with a low-threshold, Cas9-based gene drive. When released into a wild vector population, GEMs preferentially transmit these parasite-blocking genes to their offspring, ultimately modifying a vector population into a nonvector one. Deploying this technology awaits ecologically contained field trial evaluations. Here, we consider a process for site selection, the first critical step in designing a trial. Our goal is to identify a site that maximizes prospects for success, minimizes risk, and serves as a fair, valid, and convincing test of efficacy and impacts of a GEM product intended for large-scale deployment in Africa. We base site selection on geographic, geological, and biological, rather than social or legal, criteria. We recognize the latter as critically important but not as a first step in selecting a site. We propose physical islands as being the best candidates for a GEM field trial and present an evaluation of 22 African islands. We consider geographic and genetic isolation, biological complexity, island size, and topography and identify two island groups that satisfy key criteria for ideal GEM field trial sites.

Evidence of Local Extinction and Reintroduction of Aedes aegypti in Exeter, California

Established populations of Aedes aegypti, a mosquito vector of multiple major arthropod-borne viruses, were first found in three California (CA) cities in 2013. From 2013 to April 2021, Ae. aegypti thwarted almost all control efforts to stop its spread and expanded its range to 308 cities, including Exeter, in 22 counties in CA. Population genomic analyses have suggested that multiple genetically distinct Ae. aegypti populations were introduced into CA. However Ae. aegypti collected for the first time in 2014 in Exeter, appeared to be different from three major genetic clusters found elsewhere in CA. Due to intense control efforts by the Delta Vector Control District (DVCD), Ae. aegypti was thought to have been eliminated from Exeter in 2015. Unfortunately, it was recollected in 2018. It was not clear if the reemergence of Ae. aegypti in Exeter was derived from the bottlenecked remnants of the original 2014 Exeter population or from an independent invasion from a different population derived from surrounding areas. The goal of this work was to determine which of these scenarios occurred (recovery after bottleneck or reintroduction after elimination) and if elimination and reintroduction occurred to identify the origin of the invading population using a population genomic approach. Our results support the reintroduction after elimination hypothesis. The source of reintroduction, however, was unexpectedly from the southern CA cluster rather than from other two geographically closer central CA genetic clusters. We also conducted a knockdown resistance mutation profile, which showed Exeter 2014 had the lowest level of resistant alleles compared to the other populations, could have contributed towards DVCD’s ability to locally eliminate Ae. aegypti in 2014.

Population modification strategies for malaria vector control are uniquely resilient to observed levels of gene drive resistance alleles

Cas9/guide RNA (gRNA)-based gene drive systems are expected to play a transformative role in malaria elimination efforts., whether through population modification, in which the drive system contains parasite-refractory genes, or population suppression, in which the drive system induces a severe fitness load resulting in population decline or extinction. DNA sequence polymorphisms representing alternate alleles at gRNA target sites may confer a drive-resistant phenotype in individuals carrying them. Modeling predicts that, for observed levels of SGV at potential target sites and observed rates of de novo DRA formation, population modification strategies are uniquely resilient to DRAs. We conclude that gene drives can succeed when fitness costs incurred by drive-carrying mosquitoes are low enough to prevent strong positive selection for DRAs produced de novo or as part of the SGV and that population modification strategies are less prone to failure due to drive resistance.

Application of the Relationship-Based Model to Engagement for Field Trials of Genetically Engineered Malaria Vectors

The transition of new technologies for public health from laboratory to field is accompanied by a broadening scope of engagement challenges. Recent developments of vector control strategies involving genetically engineered mosquitoes with gene drives to assist in the eradication of malaria have drawn significant attention. Notably, questions have arisen surrounding community and regulatory engagement activities and of the need for examples of models or frameworks that can be applied to guide engagement. A relationship-based model (RBM) provides a framework that places stakeholders and community members at the center of decision-making processes, rather than as recipients of predetermined strategies, methods, and definitions. Successful RBM application in the transformation of healthcare delivery has demonstrated the importance of open dialogue and relationship development in establishing an environment where individuals are actively engaged in decision-making processes regarding their health. Although guidelines and recommendations for engagement for gene drives have recently been described, we argue here that communities and stakeholders should lead the planning, development, and implementation phases of engagement. The RBM provides a new approach to the development of ethical, transparent, and effective engagement strategies for malaria control programs.

Complete mitogenome sequence of Anopheles coustani from São Tomé island

We report the first complete mitogenome (Mt) sequence of Anopheles coustani, an understudied malaria vector in Africa. The sequence was extracted from one individual mosquito from São Tomé island. The length of the A. coustani Mt genome was 15,408 bp with 79.3% AT content. Phylogenetic analysis revealed that A. coustani is most closely related to A. sinensis (93.5% of identity); and 90.1% identical to A. gambiae complex members.

Abundance of conserved CRISPR-Cas9 target sites within the highly polymorphic genomes of Anopheles and Aedes mosquitoes

A number of recent papers report that standing genetic variation in natural populations includes ubiquitous polymorphisms within target sites for Cas9-based gene drive (CGD) and that these "drive resistant alleles" (DRA) preclude the successful application of CGD for managing these populations. Here we report the results of a survey of 1280 genomes of the mosquitoes Anopheles gambiae, An. coluzzii, and Aedes aegypti in which we determine that ~90% of all protein-encoding CGD target genes in natural populations include at least one target site with no DRAs at a frequency of ≥1.0%. We conclude that the abundance of conserved target sites in mosquito genomes and the inherent flexibility in CGD design obviates the concern that DRAs present in the standing genetic variation of mosquito populations will be detrimental to the deployment of this technology for population modification strategies.

Development of a confinable gene drive system in the human disease vector Aedes aegypti

Aedes aegypti is the principal mosquito vector for many arboviruses that increasingly infect millions of people every year. With an escalating burden of infections and the relative failure of traditional control methods, the development of innovative control measures has become of paramount importance. The use of gene drives has sparked significant enthusiasm for genetic control of mosquitoes; however, no such system has been developed in Ae. aegypti. To fill this void, here we develop several CRISPR-based split gene drives for use in this vector. With cleavage rates up to 100% and transmission rates as high as 94%, mathematical models predict that these systems could spread anti-pathogen effector genes into wild populations in a safe, confinable and reversible manner appropriate for field trials and effective for controlling disease. These findings could expedite the development of effector-linked gene drives that could safely control wild populations of Ae. aegypti to combat local pathogen transmission.

Transcontinental dispersal of Anopheles gambiae occurred from West African origin via serial founder events

The mosquito Anopheles gambiae s.s. is distributed across most of sub-Saharan Africa and is of major scientific and public health interest for being an African malaria vector. Here we present population genomic analyses of 111 specimens sampled from west to east Africa, including the first whole genome sequences from oceanic islands, the Comoros. Genetic distances between populations of A. gambiae are discordant with geographic distances but are consistent with a stepwise migration scenario in which the species increases its range from west to east Africa through consecutive founder events over the last ~200,000 years. Geological barriers like the Congo River basin and the East African rift seem to play an important role in shaping this process. Moreover, we find a high degree of genetic isolation of populations on the Comoros, confirming the potential of these islands as candidate sites for potential field trials of genetically engineered mosquitoes for malaria control.

Characterization of the complete mitogenome of Anopheles aquasalis, and phylogenetic divergences among Anopheles from diverse geographic zones

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes.

Ethanol as a potential mosquito sample storage medium for RNA preservation

Sample storage for downstream RNA analysis can be challenging in some field settings, especially where access to cryogenic materials or refrigeration/freezer facilities are limited. This has limited RNA-based studies on African malaria vectors collected in the field. We evaluated RNA quality after storing mosquito samples in three different sample preservation media over a 4-week period. Storing mosquito specimens in cold (4°C) media significantly improved yields of intact RNA. Our results indicate commercially available products perform well in keeping RNA integrity as advertised. Moreover, absolute ethanol may be an economical alternative for sample preservation that can be utilized in some resource-limited settings.

Genome-wide divergence among invasive populations of Aedes aegypti in California

BACKGROUND

In the summer of 2013, Aedes aegypti Linnaeus was first detected in three cities in central California (Clovis, Madera and Menlo Park). It has now been detected in multiple locations in central and southern CA as far south as San Diego and Imperial Counties. A number of published reports suggest that CA populations have been established from multiple independent introductions.

RESULTS

Here we report the first population genomics analyses of Ae. aegypti based on individual, field collected whole genome sequences. We analyzed 46 Ae. aegypti genomes to establish genetic relationships among populations from sites in California, Florida and South Africa. Based on 4.65 million high quality biallelic SNPs, we identified 3 major genetic clusters within California; one that includes all sample sites in the southern part of the state (South of Tehachapi mountain range) plus the town of Exeter in central California and two additional clusters in central California.

CONCLUSIONS

A lack of concordance between mitochondrial and nuclear genealogies suggests that the three founding populations were polymorphic for two main mitochondrial haplotypes prior to being introduced to California. One of these has been lost in the Clovis populations, possibly by a founder effect. Genome-wide comparisons indicate extensive differentiation between genetic clusters. Our observations support recent introductions of Ae. aegypti into California from multiple, genetically diverged source populations. Our data reveal signs of hybridization among diverged populations within CA. Genetic markers identified in this study will be of great value in pursuing classical population genetic studies which require larger sample sizes.

Complete mitogenome sequence of Aedes (Stegomyia) aegypti derived from field isolates from California and South Africa

The Aedes aegypti mitogenome (Mt) sequences of field isolates from California and South Africa revealed a deletion between position 14,522 and 14,659 of the Mt contig of the AaegL5 reference genome. The length of the mitogenome of the California isolate was 16,659 bp and had 99.0% similarity with the AaegL5 Mt contig. The South African isolate sequence was 16,600 bp long and had 97.9% similarity with the reference. The region between 1496 and 1664 bp is similar to a nuclear pseudogene that might be a copy of a portion of the mitochondrial genome.

Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii

BACKGROUND

Malaria mortality rates in sub-Saharan Africa have declined significantly in recent years as a result of increased insecticide-treated bed net (ITN) usage. A major challenge to further progress is the emergence and spread of insecticide resistance alleles in the Anopheles mosquito vectors, like An. coluzzii. A non-synonymous mutation in the para voltage-gated sodium channel gene reduces pyrethroid-binding affinity, resulting in knockdown resistance (kdr). Metabolic mechanisms of insecticide resistance involving detoxification genes like cytochrome P450 genes, carboxylesterases, and glutathione S-transferases are also important. As some gene activity is tissue-specific and/or environmentally induced, gene regulatory variation may be overlooked when comparing expression from whole mosquito bodies under standard rearing conditions.

RESULTS

We detected complex insecticide resistance in a 2014 An. coluzzii colony from southern Mali using bottle bioassays. Additional bioassays involving recombinant genotypes from a cross with a relatively susceptible 1995 An. coluzzii colony from Mali confirmed the importance of kdr and associated increased permethrin resistance to the CYP9K1 locus on the X chromosome. Significant differential expression of CYP9K1 was not observed among these colonies in Malpighian tubules. However, the P450 gene CYP6Z1 was overexpressed in resistant individuals following sublethal permethrin exposure and the carboxylesterase gene COEAE5G was constitutively overexpressed.

CONCLUSIONS

The significant P450-related insecticide resistance observed in the 2014 An. coluzzii colony indicates that ITNs treated with the P450 inhibitor piperonyl butoxide (PBO) would be more effective in this region. The known insecticide resistance gene CYP6Z1 was differentially expressed exclusively in the context of sublethal permethrin exposure, highlighting the importance of tissue-specificity and environmental conditions in gene expression studies. The increased activity of the carboxylesterase COEAE5G in the resistant An. coluzzii colony suggests resistance to other insecticides like organophosphates. Additional gene expression studies involving other tissues (e.g. fat body) would provide a more comprehensive view of genes underlying metabolic insecticide resistance in An. coluzzii from Mali. Identifying genetic markers linked to these regulatory alleles is an important next step that would substantially improve insecticide resistance surveillance and population genetic studies in this important vector species.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.

Anopheles darlingi polytene chromosomes: revised maps including newly described inversions and evidence for population structure in Manaus

Salivary gland polytene chromosomes of 4th instar Anopheles darlingi Root were examined from multiple locations in the Brazilian Amazon. Minor modifications were made to existing polytene photomaps. These included changes to the breakpoint positions of several previously described paracentric inversions and descriptions of four new paracentric inversions, two on the right arm of chromosome 3 and two on the left arm of chromosome 3 that were found in multiple locations. A total of 18 inversions on the X (n = 1) chromosome, chromosome 2 (n = 7) and 3 (n = 11) were scored for 83 individuals from Manaus, Macapá and Porto Velho municipalities. The frequency of 2Ra inversion karyotypes in Manaus shows significant deficiency of heterozygotes (p < 0.0009). No significant linkage disequilibrium was found between inversions on chromosome 2 and 3. We hypothesize that at least two sympatric subpopulations exist within the An. darlingi population at Manaus based on inversion frequencies.

The Genetic Basis of Host Preference and Resting Behavior in the Major African Malaria Vector, Anopheles arabiensis

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of "SNP heritability" for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ2, p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.

Surveillance, insecticide resistance and control of an invasive Aedes aegypti (Diptera: Culicidae) population in California

The invasion and subsequent establishment in California of Aedes aegypti in 2013 has created new challenges for local mosquito abatement and vector control districts. Studies were undertaken to identify effective and economical strategies to monitor the abundance and spread of this mosquito species as well as for its control. Overall, BG Sentinel (BGS) traps were found to be the most sensitive trap type to measure abundance and spread into new locations. Autocidal-Gravid-Ovitraps (AGO-B), when placed at a site for a week, performed equally to BGS in detecting the presence of female Ae. aegypti. Considering operational cost and our findings, we recommend use of BGS traps for surveillance in response to service requests especially in locations outside the known infestation area. We recommend AGO-Bs be placed at fixed sites, cleared and processed once a week to monitor mosquito abundance within a known infestation area. Long-term high density placements of AGO-Bs were found to show promise as an environmentally friendly trap-kill control strategy. California Ae. aegypti were found to be homozygous for the V1016I mutation in the voltage gated sodium channel gene, which is implicated to be involved in insecticide resistance. This strain originating from Clovis, California was resistant to some pyrethroids but not to deltamethrin in bottle bio-assays. Sentinel cage ultra-low-volume (ULV) trials using a new formulation of deltamethrin (DeltaGard®) demonstrated that it provided some control (average of 56% death in sentinel cages in a 91.4 m spray swath) after a single truck mounted aerial ULV application in residential areas.

Investigating associations between biting time in the malaria vector Anopheles arabiensis Patton and single nucleotide polymorphisms in circadian clock genes: support for sub-structure among An. arabiensis in the Kilombero valley of Tanzania

Background

There is growing evidence that the widespread use of Long-Lasting Insecticidal Nets (LLINs) is prompting malaria vectors to shift their biting towards times and places where people are not protected, such as earlier in the evening and/or outdoors. It is uncertain whether these behavioural shifts are due to phenotypic plasticity and/or ecological changes within vector communities that favour more exophilic species, or involve genetic factors within vector species to limit their contact with LLINs. Possibly variation in the time and location of mosquito biting has a genetic basis, but as yet this phenomenon has received little investigation. Here we used a candidate gene approach to investigate whether polymorphisms in selected circadian clock genes could explain variation in the time and location of feeding (indoors versus outside) within a natural population of the major African malaria vector Anopheles arabiensis.

Methods

Host-seeking An. arabiensis were collected from two villages (Lupiro and Sagamaganga) in Tanzania by Human Landing Catch (HLC) technique. Mosquitoes were classified into phenotypes of “early” (7 pm–10 pm) or “late” biting (4 am –7 am), and host-seeking indoors or outdoors. In these samples we genotyped 34 coding SNPs in 8 clock genes (PER, TIM, CLK, CYC, PDP1, VRI, CRY1, and CRY2), and tested for associations between these SNPs and biting phenotypes. SNPs in 8 mitochondrial genes (ATP6, ATP8, COX1, COX2, COX3, ND3, ND5 and CYTB) were also genotyped to test population subdivision within An. arabiensis.

Results

The candidate clock genes exhibited polymorphism within An. arabiensis, but it was unrelated to variation in the timing and location of their biting activity. However, there was evidence of strong genetic structure within An. arabiensis populations in association with the TIM, which was unrelated to geographic distance. Substructure within An. arabiensis was also detected using mitochondrial markers.

Conclusions

The variable timing and location of biting in An. arabiensis could not be linked to candidate clock genes that are known to influence behaviour in other Diptera. This finding does not rule out the possibility of a genetic basis to biting behaviour in this malaria vector, but suggests these are complex phenotypes that require more intensive ecological, neuronal and genomic analyses to understand.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1394-8) contains supplementary material, which is available to authorized users.

Surveillance, insecticide resistance and control of an invasive Aedes aegypti (Diptera: Culicidae) population in California

The invasion and subsequent establishment in California of Aedes aegypti in 2013 has created new challenges for local mosquito abatement and vector control districts. Studies were undertaken to identify effective and economical strategies to monitor the abundance and spread of this mosquito species as well as for its control. Overall, BG Sentinel (BGS) traps were found to be the most sensitive trap type to measure abundance and spread into new locations. Autocidal-Gravid-Ovitraps (AGO-B), when placed at a site for a week, performed equally to BGS in detecting the presence of female Ae. aegypti. Considering operational cost and our findings, we recommend use of BGS traps for surveillance in response to service requests especially in locations outside the known infestation area. We recommend AGO-Bs be placed at fixed sites, cleared and processed once a week to monitor mosquito abundance within a known infestation area. Long-term high density placements of AGO-Bs were found to show promise as an environmentally friendly trap-kill control strategy. California Ae. aegypti were found to be homozygous for the V1016I mutation in the voltage gated sodium channel gene, which is implicated to be involved in insecticide resistance. This strain originating from Clovis, California was resistant to some pyrethroids but not to deltamethrin in bottle bio-assays. Sentinel cage ultra-low-volume (ULV) trials using a new formulation of deltamethrin (DeltaGard®) demonstrated that it provided some control (average of 56% death in sentinel cages in a 91.4 m spray swath) after a single truck mounted aerial ULV application in residential areas.

Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali

In certain cases, a species may have access to important genetic variation present in a related species via adaptive introgression. These novel alleles may interact with their new genetic background, resulting in unexpected phenotypes. In this study, we describe a selective sweep on standing variation on the X chromosome in the mosquito Anopheles coluzzii, a principal malaria vector in West Africa. This event may have been influenced by the recent adaptive introgression of the insecticide resistance gene known as kdr from the sister species Anopheles gambiae. Individuals carrying both kdr and a nearly fixed X-linked haplotype, encompassing at least four genes including the P450 gene CYP9K1 and the cuticular protein CPR125, have rapidly increased in relative frequency. In parallel, a reproductively isolated insecticide-susceptible A. gambiae population (Bamako form) has been driven to local extinction, likely due to strong selection from increased insecticide-treated bed net usage.

Defining Genetic, Taxonomic, and Geographic Boundaries Among Species of the Psorophora confinnis (Diptera: Culicidae) Complex in North and South America

The Psorophora confinnis complex is currently composed of three species--Psorophora confinnis sensu stricto (Lynch Arribalzaga) in South America, Psorophora columbiae (Dyar and Knab) in North America, and Psorophora jamaicensis (Theobald) in the Caribbean. Members of the complex are of considerable importance as vectors of arboviruses, for example, Venezuelan equine encephalitis virus, and are significant biting pests throughout their range. The biological and geographic boundaries of Ps. confinnis and Ps. columbiae are unclear. In fact, the name Ps. columbiae is presently designated as "provisional." In this article, we aim to clarify the taxonomy and geographic distributions of species within the Ps. confinnis complex. A population genetics approach was employed using gene and genotypic frequency data at 26 isozyme loci. The results suggest that the Ps. confinnis complex in North and South America is composed of four species. Ps. confinnis s.s. and Ps. columbiae are distinct species in South and North America, respectively. Populations in Colombia, South America, formally designated as Ps. funiculus (Dyar) and populations in the southwestern United States and western Mexico, formally designated Ps. toltecum (Dyar and Knab), are distinct species. Psorophora toltecum and Psorophora funiculus species names should be resurrected from synonymy. In addition we identified a Ps. columbiae and Ps. toltecum hybrid zone in central Texas in a region described as being one of 13 North American suture zones, being geographical areas in which closely related species occur in sympatry and frequently hybridize.

A multi-detection assay for malaria transmitting mosquitoes

The Anopheles gambiae species complex includes the major malaria transmitting mosquitoes in Africa. Because these species are of such medical importance, several traits are typically characterized using molecular assays to aid in epidemiological studies. These traits include species identification, insecticide resistance, parasite infection status, and host preference. Since populations of the Anopheles gambiae complex are morphologically indistinguishable, a polymerase chain reaction (PCR) is traditionally used to identify species. Once the species is known, several downstream assays are routinely performed to elucidate further characteristics. For instance, mutations known as KDR in a para gene confer resistance against DDT and pyrethroid insecticides. Additionally, enzyme-linked immunosorbent assays (ELISAs) or Plasmodium parasite DNA detection PCR assays are used to detect parasites present in mosquito tissues. Lastly, a combination of PCR and restriction enzyme digests can be used to elucidate host preference (e.g., human vs. animal blood) by screening the mosquito bloodmeal for host-specific DNA. We have developed a multi-detection assay (MDA) that combines all of the aforementioned assays into a single multiplex reaction genotyping 33SNPs for 96 or 384 samples at a time. Because the MDA includes multiple markers for species, Plasmodium detection, and host blood identification, the likelihood of generating false positives or negatives is greatly reduced from previous assays that include only one marker per trait. This robust and simple assay can detect these key mosquito traits cost-effectively and in a fraction of the time of existing assays.

Colonization of malaria vectors under semi-field conditions as a strategy for maintaining genetic and phenotypic similarity with wild populations

Background

Malaria still accounts for an estimated 207 million cases and 627,000 deaths worldwide each year. One proposed approach to complement existing malaria control methods is the release of genetically-modified (GM) and/or sterile male mosquitoes. As opposed to laboratory colonization, this requires realistic semi field systems to produce males that can compete for females in nature. This study investigated whether the establishment of a colony of the vector Anopheles arabiensis under more natural semi-field conditions can maintain higher levels of genetic diversity than achieved by laboratory colonization using traditional methods.

Methods

Wild females of the African malaria vector An. arabiensis were collected from a village in southern Tanzania and used to establish new colonies under different conditions at the Ifakara Health Institute. Levels of genetic diversity and inbreeding were monitored in colonies of An. arabiensis that were simultaneously established in small cage colonies in the SFS and in a large semi-field (SFS) cage and compared with that observed in the original founder population. Phenotypic traits that determine their fitness (body size and energetic reserves) were measured at 10^th generation and compared to founder wild population.

Results

In contrast to small cage colonies, the SFS population of An. arabiensis exhibited a higher degree of similarity to the founding field population through time in several ways: (i) the SFS colony maintained a significantly higher level of genetic variation than small cage colonies, (ii) the SFS colony had a lower degree of inbreeding than small cage colonies, and (iii) the mean and range of mosquito body size in the SFS colony was closer to that of the founding wild population than that of small cage colonies. Small cage colonies had significantly lower lipids and higher glycogen abundances than SFS and wild population.

Conclusions

Colonization of An. arabiensis under semi-field conditions was associated with the retention of a higher degree of genetic diversity, reduced inbreeding and greater phenotypic similarity to the founding wild population than observed in small cage colonies. Thus, mosquitoes from such semi-field populations are expected to provide more realistic representation of mosquito ecology and physiology than those from small cage colonies.

Plasmodium falciparum infection rates for some Anopheles spp. from Guinea-Bissau, West Africa

Presence of Plasmodiumfalciparum circumsporozoite protein (CSP) was detected by enzyme linked immunosorbent assay (ELISA) in a sample of Anophelesgambiae s.s., A. melas and A. pharoensis collected in Guinea-Bissau during October and November 2009. The percentage of P. falciparum infected samples (10.2% overall; confidence interval (CI): 7.45-13.6%) was comparable to earlier studies from other sites in Guinea-Bissau (9.6-12.4%). The majority of the specimens collected were identified as A. gambiae which had an individual infection rate of 12.6 % (CI: 8.88-17.6) across collection sites. A small number of specimens of A. coluzzii, A. coluzzii x A. gambiae hybrids, A. melas and A. pharoensis were collected and had infection rates of 4.3% (CI:0.98-12.4), 4.1% (CI:0.35-14.5), 11.1% (CI:1.86-34.1) and 33.3% (CI:9.25-70.4) respectively. Despite being present in low numbers in indoor collections, the exophilic feeding behaviors of A. melas (N=18) and A. pharoensis (N=6) and high infection rates observed in this survey suggest falciparum-malaria transmission potential outside of the protection of bed nets.

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR-based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single-nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so-called islands of divergence leading us to describe this method as the 'Divergence Island SNP' (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon

Background

Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance.

Methods

A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali.

Results

Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion.

Discussion

More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

An analysis of two island groups as potential sites for trials of transgenic mosquitoes for malaria control

Considerable technological advances have been made towards the generation of genetically modified mosquitoes for vector control. In contrast, less progress has been made towards field evaluations of transformed mosquitoes which are critical for evaluating the success of, and hazards associated with, genetic modification. Oceanic islands have been highlighted as potentially the best locations for such trials. However, population genetic studies are necessary to verify isolation. Here, we used a panel of genetic markers to assess for evidence of genetic isolation of two oceanic island populations of the African malaria vector, Anopheles gambiae s.s. We found no evidence of isolation between the Bijagós archipelago and mainland Guinea-Bissau, despite separation by distances beyond the known dispersal capabilities of this taxon. Conversely, the Comoros Islands appear to be genetically isolated from the East African mainland, and thus represent a location worthy of further investigation for field trials. Based on assessments of gene flow within and between the Comoros islands, the island of Grande Comore was found to be genetically isolated from adjacent islands and also exhibited local population structure, indicating that it may be the most suitable site for trials with existing genetic modification technologies.

Chromosome inversions, genomic differentiation and speciation in the African malaria mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.

Altitudinal genetic and morphometric variation among populations of Culex theileri Theobald (Diptera: Culicidae) from northeastern Turkey

Enviromental conditions, including such important climatic variables as temperature and precipitation, change with altitude; thus, elevation plays a significant role in determining population and community structure in a variety of organisms. Using single nucleotide polymorphisms (SNPs) and geometric morphometrics, nine populations of Culex theileri Theobald occurring in different ecological subregions at altitudes between 808-2,130 m in northeastern Turkey were compared. The wing size and shape data indicate that there are significant phenotypic differences among them, while Cx theileri populations are not genetically differentiated in the northeast part of Turkey. The size and shape variation analysis of wings showed that there is a positive correlation between wing (body) size/shape and altitude.

Identification and characterization of single nucleotide polymorphisms (SNPs) in Culex theileri (Diptera: Culicidae)

Culex theileri Theobald (Diptera: Culicidae) is one of the most common mosquito species in northeastern Turkey and serves as a vector for various zoonotic diseases including West Nile virus. Although there have been some studies on the ecology of Cx. theileri, very little genetic data has been made available. We successfully sequenced 11 gene fragments from Cx. theileri specimens collected from the northeastern part of Turkey. On average, we found a Single nucleotide polymorphism every 45 bp. Transitions outnumbered transversions, at a ratio of 2:1. This is the first report of genetic polymorphisms in Cx. theileri and Single nucleotide polymorphism discovered from this study can be used to investigate population structure and gene-environmental interactions.

Differential Plasmodium falciparum infection of Anopheles gambiae s.s. molecular and chromosomal forms in Mali

Background

Anopheles gambiae sensu stricto (s.s.) is a primary vector of Plasmodium falciparum in sub-Saharan Africa. Although some physiological differences among molecular and chromosomal forms of this species have been demonstrated, the relative susceptibility to malaria parasite infection among them has not been unequivocally shown. The objective of this study was to investigate P. falciparum circumsporozoite protein infection (CSP) positivity among An. gambiae s.s. chromosomal and molecular forms.

Methods

Wild An. gambiae from two sites Kela (n = 464) and Sidarebougou (n = 266) in Mali were screened for the presence of P. falciparum CSP using an enzyme-linked immunosorbent assay (ELISA). Samples were then identified to molecular form using multiple PCR diagnostics (n = 713) and chromosomal form using chromosomal karyotyping (n = 419).

Results

Of 730 An. gambiae sensu lato (s.l.) mosquitoes, 89 (12.2%) were CSP ELISA positive. The percentage of positive mosquitoes varied by site: 52 (11.2%) in Kela and 37 (13.9%) in Sidarebougou. Eighty-seven of the positive mosquitoes were identified to molecular form and they consisted of nine Anopheles arabiensis (21.4%), 46 S (10.9%), 31 M (12.8%), and one MS hybrid (14.3%). Sixty of the positive mosquitoes were identified to chromosomal form and they consisted of five An. arabiensis (20.0%), 21 Savanna (15.1%), 21 Mopti (30.4%), 11 Bamako (9.2%), and two hybrids (20.0%).

Discussion

In this collection, the prevalence of P. falciparum infection in the M form was equivalent to infection in the S form (no molecular form differential infection). There was a significant differential infection by chromosomal form such that, P. falciparum infection was more prevalent in the Mopti chromosomal forms than in the Bamako or Savanna forms; the Mopti form was also the most underrepresented in the collection. Continued research on the differential P. falciparum infection of An. gambiae s.s. chromosomal and molecular forms may suggest that Plasmodium – An. gambiae interactions play a role in malaria transmission.

High degree of single nucleotide polymorphisms in California Culex pipiens (Diptera: Culicidae) sensu lato

Resolution of systematic relationships among members of the Culex pipiens (L.) complex has important implications for public health as well as for studies on the evolution of sibling species. Currently held views contend that in California considerable genetic introgression occurs between Cx. pipiens and Cx. quinquefasciatus Say, and as such, these taxa behave as if they are a single species. Development of high throughput SNP genotyping tools for the analysis of Cx. pipiens complex population structure is therefore desirable. As a first step toward this goal, we sequenced 12 gene fragments from specimens collected in Marin and Fresno counties. On average, we found a higher single nucleotide polymorphism (SNP) density than any other mosquito species reported thus far. Coding regions contained significantly higher GC content (median 54.7%) than noncoding regions (42.4%; Wilcoxon rank sum test, P = 5.29 x 10(-5)). Differences in SNP allele frequencies observed between mosquitoes from Marin and Fresno counties indicated significant genetic divergence and suggest that SNP markers will be useful for future detailed population genetic studies of this group. The high density of SNPs highlights the difficulty in identifying species within the complex and may be associated with the large degree of phenotypic variation observed in this group of mosquitoes.

Single-nucleotide polymorphisms for high-throughput genotyping of Anopheles arabiensis in East and southern Africa

Anopheles arabiensis Patton is one of the principal vectors of malaria in sub-Saharan Africa, occupying a wide variety of ecological zones. This species is increasingly responsible for malaria transmission in Africa and is becoming the dominant vector species in some localities. Despite its growing importance, little is known about genetic polymorphisms in this species. Multiple sequences of various gene fragments from An. arabiensis isolates from Cameroon were obtained from GenBank. In total, 20 gene fragments containing single-nucleotide polymorphisms (SNPs) at moderate density were selected for direct sequencing from field collected specimens from Tanzania and Zambia. We obtained 301 SNPs in total from the 20 gene fragments, 60 of which were suitable for Illumina GoldenGate SNP genotyping. A greater number of SNPs (n = 185) was suitable for analysis using Sequenom iPLEX, an alternative high-throughput genotyping technology using mass spectrometry. An SNP was present every 59 (+/- 44.5) bases on average. Overall, An. arabiensis from Tanzania and Zambia are genetically closer (mean F(ST) = 0.075) than either is to populations in Cameroon (F(ST, TZ-CM) = 0.250, F(ST,ZA-CM) = 0.372). A fixed polymorphism between East/southern and Central Africa was identified on AGAP000574, a gene on the X chromosome. We have identified SNPs in natural populations of An. arabiensis. SNP densities in An. arabiensis were higher than Anopheles gambiae s.s., suggesting a greater challenge in the development of high-throughput SNP analysis for this species. The SNP markers provided in this study are suitable for a high-throughput genotyping analysis and can be used for population genetic studies and association mapping efforts.

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.

Morphological differentiation may mediate mate-choice between incipient species of Anopheles gambiae s.s

The M and S molecular forms of Anopheles gambiae s.s. have been considered incipient species for more than ten years, yet the mechanism underlying assortative mating of these incipient species has remained elusive. The discovery of the importance of harmonic convergence of wing beat frequency in mosquito mating and its relation to wing size have laid the foundation for exploring phenotypic divergence in wing size of wild populations of the two forms. In this study, wings from field collected mosquitoes were measured for wing length and wing width from two parts of the sympatric distribution, which differ with respect to the strength of assortative mating. In Mali, where assortative mating is strong, as evidenced by low rates of hybridization, mean wing lengths and wing widths were significantly larger than those from Guinea-Bissau. In addition, mean wing widths in Mali were significantly different between molecular forms. In Guinea-Bissau, assortative mating appears comparatively reduced and wing lengths and widths did not differ significantly between molecular forms. The data presented in this study support the hypothesis that wing beat frequency may mediate assortative mating in the incipient species of A. gambiae and represent the first documentation of a morphological difference between the M and S molecular forms.

Population genetic structure of Anopheles arabiensis and Anopheles gambiae in a malaria endemic region of southern Tanzania

BACKGROUND

Genetic diversity is a key factor that enables adaptation and persistence of natural populations towards environmental conditions. It is influenced by the interaction of a natural population's dynamics and the environment it inhabits. Anopheles gambiae s.s. and Anopheles arabiensis are the two major and widespread malaria vectors in sub-Saharan Africa. Several studies have examined the ecology and population dynamics of these vectors. Ecological conditions along the Kilombero valley in Tanzania influence the distribution and population density of these two vector species. It remains unclear whether the ecological diversity within the Kilombero valley has affected the population structure of An. gambiae s.l. populations. The goal of this study was to characterise the genetic structure of sympatric An. gambiae s.s and An. arabiensis populations along the Kilombero valley.

METHODOLOGY

Mosquitoes were collected from seven locations in Tanzania: six from the Kilombero valley and one outside the valley (-700 km away) as an out-group. To archive a genome-wide coverage, 13 microsatellite markers from chromosomes X, 2 and 3 were used.

RESULTS

High levels of genetic differentiation among An. arabiensis populations was observed, as opposed to An. gambiae s.s., which was genetically undifferentiated across the 6,650 km2 of the Kilombero valley landscape. It appears that genetic differentiation is not attributed to physical barriers or distance, but possibly by ecological diversification within the Kilombero valley. Genetic divergence among An. arabiensis populations (FST = 0.066) was higher than that of the well-known M and S forms of An. gambiae s. s. in West and Central Africa (FST = 0.035), suggesting that these populations are maintained by some level of reproductive isolation.

CONCLUSION

It was hypothesized that ecological diversification across the valley may be a driving force for observed An. arabiensis genetic divergence. The impact of the observed An. arabiensis substructure to the prospects for new vector control approaches is discussed.

Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics

Culex quinquefasciatus (the southern house mosquito) is an important mosquito vector of viruses such as West Nile virus and St. Louis encephalitis virus, as well as of nematodes that cause lymphatic filariasis. C. quinquefasciatus is one species within the Culex pipiens species complex and can be found throughout tropical and temperate climates of the world. The ability of C. quinquefasciatus to take blood meals from birds, livestock, and humans contributes to its ability to vector pathogens between species. Here, we describe the genomic sequence of C. quinquefasciatus: Its repertoire of 18,883 protein-coding genes is 22% larger than that of Aedes aegypti and 52% larger than that of Anopheles gambiae with multiple gene-family expansions, including olfactory and gustatory receptors, salivary gland genes, and genes associated with xenobiotic detoxification.

Identification of three single nucleotide polymorphisms in Anopheles gambiae immune signaling genes that are associated with natural Plasmodium falciparum infection

Background

Laboratory studies have demonstrated that a variety of immune signaling pathways regulate malaria parasite infection in Anopheles gambiae, the primary vector species in Africa.

Methods

To begin to understand the importance of these associations under natural conditions, an association mapping approach was adopted to determine whether single nucleotide polymorphisms (SNPs) in selected immune signaling genes in A. gambiae collected in Mali were associated with the phenotype of Plasmodium falciparum infection.

Results

Three SNPs were identified in field-collected mosquitoes that were associated with parasite infection in molecular form-dependent patterns: two were detected in the Toll5B gene and one was detected in the gene encoding insulin-like peptide 3 precursor. In addition, one infection-associated Toll5B SNP was in linkage disequilibrium with a SNP in sequence encoding a mitogen-activated protein kinase that has been associated with Toll signaling in mammalian cells. Both Toll5B SNPs showed divergence from Hardy-Weinberg equilibrium, suggesting that selection pressure(s) are acting on these loci.

Conclusions

Seven of these eight infection-associated and linked SNPs alter codon frequency or introduce non-synonymous changes that would be predicted to alter protein structure and, hence, function, suggesting that these SNPs could alter immune signaling and responsiveness to parasite infection.

Ecological and genetic relationships of the Forest-M form among chromosomal and molecular forms of the malaria vector Anopheles gambiae sensu stricto

BACKGROUND

Anopheles gambiae sensu stricto, one of the principal vectors of malaria, has been divided into two subspecific groups, known as the M and S molecular forms. Recent studies suggest that the M form found in Cameroon is genetically distinct from the M form found in Mali and elsewhere in West Africa, suggesting further subdivision within that form.

METHODS

Chromosomal, microsatellite and geographic/ecological evidence are synthesized to identify sources of genetic polymorphism among chromosomal and molecular forms of the malaria vector Anopheles gambiae s.s.

RESULTS

Cytogenetically the Forest M form is characterized as carrying the standard chromosome arrangement for six major chromosomal inversions, namely 2La, 2Rj, 2Rb, 2Rc, 2Rd, and 2Ru. Bayesian clustering analysis based on molecular form and chromosome inversion polymorphisms as well as microsatellites describe the Forest M form as a distinct population relative to the West African M form (Mopti-M form) and the S form. The Forest-M form was the most highly diverged of the An. gambiae s.s. groups based on microsatellite markers. The prevalence of the Forest M form was highly correlated with precipitation, suggesting that this form prefers much wetter environments than the Mopti-M form.

CONCLUSION

Chromosome inversions, microsatellite allele frequencies and habitat preference all indicate that the Forest M form of An. gambiae is genetically distinct from the other recognized forms within the taxon Anopheles gambiae sensu stricto. Since this study covers limited regions of Cameroon, the possibility of gene flow between the Forest-M form and Mopti-M form cannot be rejected. However, association studies of important phenotypes, such as insecticide resistance and refractoriness against malaria parasites, should take into consideration this complex population structure.

Desiccation resistance among subpopulations of Anopheles gambiae s.s. from Selinkenyi, Mali

Certain forms of Anopheles gambiae s.s. actively maintain malaria transmission in the driest areas and months of the year because of considerable drought tolerance. We monitored desiccation resistance of F1 offspring of both the M and S forms of field-collected An. gambiae s.s. Our results indicate that the geographic cline in the distribution of the two forms, as observed in Mali, corresponds to a physiological difference in response to arid environments. In addition, female mosquitoes survived significantly longer than males, enhancing the vector competence for the malaria parasite. Our study supports a genetic link to the drought tolerance phenotype, a phenotype with important consequences to malaria transmission in many places in Africa.