Speciation within Anopheles gambiae--the glass is half full

Genomic Analysis Reveals a New Cryptic Taxon Within the Anopheles gambiae Complex With a Distinct Insecticide Resistance Profile in the Coast of East Africa

Anopheles mosquitoes are major malaria vectors, encompassing several species complexes with diverse life histories, transmission risks and insecticide resistance profiles that challenge malaria control efforts. This study investigated the genetic structure and insecticide resistance profiles of Anopheles gambiae complex mosquitoes in Tanzania. We analysed whole-genome sequence data of 300 mosquitoes collected between 2012 and 2015 across four regions in northern Tanzania and identified An. gambiae s.s., An. arabiensis and a distinct taxonomic group that was previously unknown. This distinct taxon has a unique profile of genetic diversity and appears restricted to the coastal region, and we refer to it as the Pwani molecular form. Analysis of insecticide resistance based on target-site mutations and copy number variations (CNV) showed that these markers were strikingly absent from the Pwani molecular form in contrast to other taxa. Our analysis also revealed a pattern of geographical isolation in the An. gambiae s.s. populations, with samples from the north-western site (Muleba) clustering separately from those collected in the north-eastern site (Muheza). These geographically isolated subpopulations also had differing resistance and selection profiles, with An. gambiae s.s. from the north-western site showing genomic evidence of higher resistance to pyrethroids compared with the north-eastern population. Conversely, An. arabiensis showed no geographical population structuring, with a similar insecticide resistance profile across all sampling locations, suggesting unrestricted gene flow. Our findings underscore the need to incorporate genetic data into malaria vector surveillance and control decisions and could inform the development and deployment of new interventions.

Spatial modelling for population replacement of mosquito vectors at continental scale

Malaria is one of the deadliest vector-borne diseases in the world. Researchers are developing new genetic and conventional vector control strategies to attempt to limit its burden. Novel control strategies require detailed safety assessment to ensure responsible and successful deployments. Anopheles gambiae sensu stricto (s.s.) and Anopheles coluzzii, two closely related subspecies within the species complex Anopheles gambiae sensu lato (s.l.), are among the dominant malaria vectors in sub-Saharan Africa. These two subspecies readily hybridise and compete in the wild and are also known to have distinct niches, each with spatially and temporally varying carrying capacities driven by precipitation and land use factors.

We model the spread and persistence of a population-modifying gene drive system in these subspecies across sub-Saharan Africa by simulating introductions of genetically modified mosquitoes across the African mainland and its offshore islands. We explore transmission of the gene drive between the two subspecies that arise from different hybridisation mechanisms, the effects of both local dispersal and potential wind-aided migration to the spread, and the development of resistance to the gene drive. Given the best current available knowledge on the subspecies’ life histories, we find that an introduced gene drive system with typical characteristics can plausibly spread from even distant offshore islands to the African mainland with the aid of wind-driven migration, with resistance beginning to take over within a decade. Our model accounts for regional to continental scale mechanisms, and demonstrates a range of realistic dynamics including the effect of prevailing wind on spread and spatio-temporally varying carrying capacities for subspecies. As a result, it is well-placed to answer future questions relating to mosquito gene drives as important life history parameters become better understood.

Conventional control methods have dramatically reduced malaria, but it still kills over 300,000 children in Africa each year, and this number could increase as their effectiveness wanes. Novel control methods using gene drives rapidly reduce or modify malaria vector populations in laboratory settings, and hence are now being considered for field applications. We use modelling to assess how a gene drive might spread and persist in the malaria-carrying subspecies Anopheles gambiae sensu stricto (s.s.) and Anopheles coluzzii. These two subspecies interbreed and compete, so we model how these interactions affect the spread of the drive at a continental scale. In scenarios that allow mosquitoes to travel on prevailing wind currents, we find that a gene drive can potentially spread across national borders—and jump from offshore islands to the African mainland—but spread is eventually arrested when the drive allele is ousted by a resistant allele. As we learn more about the population dynamics of both genetically modified and wild mosquitoes, and as gene drive systems are further developed to allow local containment and evade resistance, our model will be able to answer more detailed questions about how they can be applied in the field effectively and safely.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Evidence for Divergent Selection on Immune Genes between the African Malaria Vectors, Anopheles coluzzii and A. gambiae

Simple Summary

A comparison of the genomes of the African malaria vectors, Anopheles gambiae and A. coluzzii, revealed that immune genes are highly diverged. Although these two species frequently co-occur within a single site, they occur in distinct larval habitats. Our results taken in the context of known differences in the larval habitats occupied by these taxa support the hypothesis that observed genetic divergence may be driven by immune response to microbial agents specific to these habitats. Strict within species mating may have subsequently evolved in part to maintain immunocompetence which might be compromised by dysregulation of immune pathways in hybrids. We conclude that the evolution of immune gene divergence among this important group of species may serve as a useful model to explore ecological speciation in general.

Abstract

During their life cycles, microbes infecting mosquitoes encounter components of the mosquito anti-microbial innate immune defenses. Many of these immune responses also mediate susceptibility to malaria parasite infection. In West Africa, the primary malaria vectors are Anopheles coluzzii and A. gambiae sensu stricto, which is subdivided into the Bamako and Savanna sub-taxa. Here, we performed whole genome comparisons of the three taxa as well as genotyping of 333 putatively functional SNPs located in 58 immune signaling genes. Genome data support significantly higher differentiation in immune genes compared with a randomly selected set of non-immune genes among the three taxa (permutation test p < 0.001). Among the 58 genes studied, the majority had one or more segregating mutations (72.9%) that were significantly diverged among the three taxa. Genes detected to be under selection include MAP2K4 and Raf. Despite the genome-wide distribution of immune genes, a high level of linkage disequilibrium (r² > 0.8) was detected in over 27% of SNP pairs. We discuss the potential role of immune gene divergence as adaptations to the different larval habitats associated with A. gambiae taxa and as a potential force driving ecological speciation in this group of mosquitoes.

Diversity, dynamics, direction, and magnitude of high-altitude migrating insects in the Sahel

Long-distance migration of insects impacts food security, public health, and conservation–issues that are especially significant in Africa. Windborne migration is a key strategy enabling exploitation of ephemeral havens such as the Sahel, however, its knowledge remains sparse. In this first cross-season investigation (3 years) of the aerial fauna over Africa, we sampled insects flying 40–290 m above ground in Mali, using nets mounted on tethered helium-filled balloons. Nearly half a million insects were caught, representing at least 100 families from thirteen orders. Control nets confirmed that the insects were captured at altitude. Thirteen ecologically and phylogenetically diverse species were studied in detail. Migration of all species peaked during the wet season every year across localities, suggesting regular migrations. Species differed in flight altitude, seasonality, and associated weather conditions. All taxa exhibited frequent flights on southerly winds, accounting for the recolonization of the Sahel from southern source populations. “Return” southward movement occurred in most taxa. Estimates of the seasonal number of migrants per species crossing Mali at latitude 14°N were in the trillions, and the nightly distances traversed reached hundreds of kilometers. The magnitude and diversity of windborne insect migration highlight its importance and impacts on Sahelian and neighboring ecosystems.

DNA barcoding and species delimitation of butterflies (Lepidoptera) from Nigeria

Accurate identification of species is a prerequisite for successful biodiversity management and further genetic studies. Species identification techniques often require both morphological diagnostics and molecular tools, such as DNA barcoding, for correct identification. In particular, the use of the subunit I of the mitochondrial cytochrome c oxidase (COI) gene for DNA barcoding has proven useful in species identification for insects. However, to date, no studies have been carried out on the DNA barcoding of Nigerian butterflies. We evaluated the utility of DNA barcoding applied for the first time to 735 butterfly specimens from southern Nigeria. In total, 699 DNA barcodes, resulting in a record of 116 species belonging to 57 genera, were generated. Our study sample comprised 807 DNA barcodes based on sequences generated from our current study and 108 others retrieved from BOLD. Different molecular analyses, including genetic distance-based evaluation (Neighbor-Joining, Maximum Likelihood and Bayesian trees) and species delimitation tests (TaxonDNA, Automated Barcode Gap Discovery, General Mixed Yule-Coalescent, and Bayesian Poisson Tree Processes) were performed to accurately identify and delineate species. The genetic distance-based analyses resulted in 163 well-separated clusters consisting of 147 described and 16 unidentified species. Our findings indicate that about 90.20% of the butterfly species were explicitly discriminated using DNA barcodes. Also, our field collections reported the first country records of ten butterfly species-Acraea serena, Amauris cf. dannfelti, Aterica galena extensa, Axione tjoane rubescens, Charaxes galleyanus, Papilio lormieri lormeri, Pentila alba, Precis actia, Precis tugela, and Tagiades flesus. Further, DNA barcodes revealed a high mitochondrial intraspecific divergence of more than 3% in Bicyclus vulgaris vulgaris and Colotis evagore. Furthermore, our result revealed an overall high haplotype (gene) diversity (0.9764), suggesting that DNA barcoding can provide information at a population level for Nigerian butterflies. The present study confirms the efficiency of DNA barcoding for identifying butterflies from Nigeria. To gain a better understanding of regional variation in DNA barcodes of this biogeographically complex area, future work should expand the DNA barcode reference library to include all butterfly species from Nigeria as well as surrounding countries. Also, further studies, involving relevant genetic and eco-morphological datasets, are required to understand processes governing mitochondrial intraspecific divergences reported in some species complexes.

Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex

Heterochromatin is identified as a potential factor driving diversification of species. To understand the magnitude of heterochromatin variation within the Anopheles gambiae complex of malaria mosquitoes, we analyzed metaphase chromosomes in An. arabiensis, An. coluzzii, An. gambiae, An. merus, and An. quadriannulatus. Using fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA), a highly repetitive fraction of DNA, and heterochromatic Bacterial Artificial Chromosome (BAC) clones, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosomes of An. gambiae. We then developed chromosome idiograms and demonstrated that the X chromosomes exhibit qualitative differences in their pattern of heterochromatic bands and position of satellite DNA (satDNA) repeats among the sibling species with postzygotic isolation, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae. The identified differences in the size and structure of the X chromosome heterochromatin point to a possible role of repetitive DNA in speciation of mosquitoes. We found that An. coluzzii and An. gambiae, incipient species with prezygotic isolation, share variations in the relative positions of the satDNA repeats and the proximal heterochromatin band on the X chromosomes. This previously unknown genetic polymorphism in malaria mosquitoes may be caused by a differential amplification of DNA repeats or an inversion in the sex chromosome heterochromatin.

Chromosome and Genome Divergence between the Cryptic Eurasian Malaria Vector-Species Anopheles messeae and Anopheles daciae

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anophelesmesseae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation.

An Evolutionary Perspective on Vector-Borne Diseases

Several aspects of the biology of the three players in a vector-borne disease that affect their evolutionary interactions are outlined. A model of the origin of a human-human cycle of vector-borne diseases is presented emphasizing the narrowing of the niche experienced by the pathogen and vector. Variation in the expected rates of evolution of the three players is discussed with the rapid rate of pathogen evolution providing them with advantages. Population sizes and fluctuations also affect the three players in very different ways. The time since the origin of a vector-borne disease likely determines how stable the interactions are and thus how easily the disease might be eliminated. Stability and variation are also linked. Human technological advances are rapidly upsetting the previously relatively slow coevolutionary adjustment of the three players. Finally, it is pointed out that development of quantitative coevolutionary models specifically addressing details of vector-borne diseases is needed to identify parameters most likely to break transmission cycles and thus control or eliminate diseases.

Effect of water source and feed regime on development and phenotypic quality in Anopheles gambiae (s.l.): prospects for improved mass-rearing techniques towards release programmes

Background

In many malaria-endemic sub-Saharan countries, insecticide resistance poses a threat to existing mosquito control measures, underscoring the need for complementary control methods such as sterile and/or genetically-modified mosquito release programmes. The sibling species Anopheles gambiae and An. coluzzii are responsible for malaria transmission in most of this region. In their natural habitat, these species generally breed in clean, soft water and it is believed that divergent preference in their larval breeding sites have played a role in their speciation process. Mosquito release programmes rely on the rearing of mosquitoes at high larval densities. Current rearing protocols often make use of deionised water regardless of the strain reared. They also depend on a delicate balance between the need for adequate feeding and the negative effect of toxic ammonia and food waste build-up on mosquito development, making managing and improving water quality in the insectary imperative.

Methods

Here, we investigated the impact of water source and feed regimes on emergence rate and phenotypic quality of mosquitoes in the insectary. First-instar larvae of An. gambiae (Kisumu strain) and An. coluzzii (Mopti and VK3 strains) were reared in three water sources with varying degrees of hardness (deionised, mineral and a mix of the two), with a daily water change. Larvae were fed daily using two standardised feeding regimes, solution and powder feed.

Results

Water source had a significant impact on mosquito size and development time for all strains. Earlier emergence of significantly larger mosquitoes was observed in mineral water with the smallest mosquitoes developing later from deionised water. Wing-length was significantly longer in mineral, mixed water and in powder feed, irrespective of sex, strains or water types. Deionised water was the least favourable for mosquito quality across all strains.

Conclusions

Mineral water and powder feed should be used in rearing protocols to improve mosquito quality where the optimal quality of mosquitoes is desired. Although results obtained were not significant for improved mosquito numbers, the phenotypic quality of mosquitoes reared was significantly improved in mineral water and mix water. Further studies are recommended on the impact mineral water has on other fitness traits such as longevity, fecundity and mating competitiveness.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3465-0) contains supplementary material, which is available to authorized users.

Mosquito microevolution drives Plasmodium falciparum dynamics

Malaria, a major cause of child mortality in Africa, is engendered by Plasmodium parasites that are transmitted by anopheline mosquitoes. Fitness of Plasmodium parasites is closely linked to the ecology and evolution of its anopheline vector. However, whether the genetic structure of vector populations impacts malaria transmission remains unknown. Here, we describe a partitioning of the African malaria vectors into generalists and specialists that evolve along ecological boundaries. We next identify the contribution of mosquito species to Plasmodium abundance using Granger causality tests for time-series data collected over two rainy seasons in Mali. We find that mosquito microevolution, defined by changes in the genetic structure of a population over short ecological timescales, drives Plasmodium dynamics in nature, whereas vector abundance, infection prevalence, temperature and rain have low predictive values. Our study demonstrates the power of time-series approaches in vector biology and highlights the importance of focusing local vector control strategies on mosquito species that drive malaria dynamics.

Vector bionomics and malaria transmission in an area of sympatry of An. arabiensis, An. coluzzii and An. gambiae

Despite extensive genetic studies on their variability and differentiation, few is known about the specific and relative role of An. coluzzii, An. gambiae and An. arabiensis in areas of sympatry. Indeed, their behavioral dissimilarities and divergent population dynamics can impact on malaria transmission level and intensity. This study was undertaken in four sympatric sites belonging to two different ecosystems with differential insecticide pressure to study the bionomics of these species and their relative role in malaria transmission. Mosquitoes were collected monthly from July to December 2011 when landing on human volunteers and by pyrethrum spray catches. Specimens belonging to the An. gambiae complex were further identified using molecular tools. Plasmodium falciparum infection and blood-feeding preferences were studied using the ELISA techniques. Overall, the three species were in sympatry in each of the four sites with the predominance of An. gambiae. Mosquito populations' dynamics varied temporally depending on the rainy season for each zone. The anthropophilic rates varied between 45.7 and 78.1% for An. arabiensis, 81.8 and 100% for An. coluzzii and 80 and 96.7% for An. gambiae. Plasmodium infection rates were higher in An. gambiae (range: 2.17%-6.54%) while for An. arabiensis and An. coluzzii it varied respectively between 0-1.24% and 0-3.66%. Malaria transmission occured in each of the four sites both indoors and outdoors and was due mainly to An. gambiae. An. arabiensis and An. coluzzii played a limited role due both to a low anthropophilic rate and a lower biting rate for An. coluzzii in comparison with An. gambiae. This study showed that, while present in sympatric areas, species from the An. gambiae complex could exhibit differential involvement in malaria transmission. Even less involved in malaria transmission, the occurrence of ecological and environmental changes tending to a good adaptation of An. coluzzii could lead to a great risk for malaria transmission in time and space in human populations.

Evaluating the accuracy of morphological identification of insect pests of rice crops using DNA barcoding

Accurate identification of agricultural pests is key requirement for the successful integrated pest management (IPM) program. However, due to limitations of conventional morphological methods, other molecular method like DNA barcoding is used. The current study was designed to evaluate the accuracy of morphological identification of insect pests using DNA barcoding. Morphologically, a total of 247 insect pests, representing 10 families, 18 genera, 22 species were identified. However, molecular identifications confirmed the presence of 11 families, 16 genera, and 20 species of agricultural pests. A total of 59 specimens were processed for DNA barcoding but genomic sequences of mt COI gene up to 600 bp were revived from 48 samples. Specimens that were misidentified through morphological studies were placed to their appropriate taxon, using this molecular approach. Results revealed the existence of clear barcode gap for different pest species. Moreover, the values of distance with the nearest neighbour recorded were higher than the maximum intra-sequence divergences for all species. It is concluded that DNA barcoding is a reliable technique for identification of agricultural pests, especially for immature stages when morphometric studies are ambiguous and will be helpful in the development of more effective pest management options for regulating pest species.

Emergence of knock-down resistance in the Anopheles gambiae complex in the Upper River Region, The Gambia, and its relationship with malaria infection in children

BACKGROUND

Insecticide resistance threatens malaria control in sub-Saharan Africa. Knockdown resistance to pyrethroids and organochlorines in Anopheles gambiae sensu lato (s.l.) is commonly caused by mutations in the gene encoding a voltage-gated sodium channel which is the target site for the insecticide. The study aimed to examine risk factors for knockdown resistance in An. gambiae s.l. and its relationship with malaria infection in children in rural Gambia. Point mutations at the Vgsc-1014 locus, were measured in An. gambiae s.l. during a 2-year trial. Cross-sectional surveys were conducted at the end of the transmission season to measure malaria infection in children aged 6 months-14 years.

RESULTS

Whilst few Anopheles arabiensis and Anopheles coluzzii had Vgsc-1014 mutations, the proportion of An. gambiae sensu stricto (s.s.) mosquitoes homozygous for the Vgsc-1014F mutation increased from 64.8 to 90.9% during the study. The Vgsc-1014S or 1014F mutation was 80% higher in 2011 compared to 2010, and 27% higher in the villages with indoor residual spraying compared to those without. An increase in the proportion of An. gambiae s.l. mosquitoes with homozygous Vgsc-1014F mutations and an increase in the proportion of An. gambiae s.s. in a cluster were each associated with increased childhood malaria infection. Homozygous Vgsc-1014F mutations were, however, most common in An. gambiae s.s. and almost reached saturation during the study meaning that the two variables were colinear.

CONCLUSIONS

As a result of colinearity between homozygous Vgsc-1014F mutations and An. gambiae s.s., it was not possible to determine whether insecticide resistance or species composition increased the risk of childhood malaria infection.

Rapid selection of a pyrethroid metabolic enzyme CYP9K1 by operational malaria control activities

Since 2004, indoor residual spraying (IRS) and long-lasting insecticide-impregnated bednets (LLINs) have reduced the malaria parasite prevalence in children on Bioko Island, Equatorial Guinea, from 45% to 12%. After target site-based (knockdown resistance; kdr) pyrethroid resistance was detected in 2004 in Anopheles coluzzii (formerly known as the M form of the Anopheles gambiae complex), the carbamate bendiocarb was introduced. Subsequent analysis showed that kdr alone was not operationally significant, so pyrethroid-based IRS was successfully reintroduced in 2012. In 2007 and 2014-2015, mass distribution of new pyrethroid LLINs was undertaken to increase the net coverage levels. The combined selection pressure of IRS and LLINs resulted in an increase in the frequency of pyrethroid resistance in 2015. In addition to a significant increase in kdr frequency, an additional metabolic pyrethroid resistance mechanism had been selected. Increased metabolism of the pyrethroid deltamethrin was linked with up-regulation of the cytochrome P450 CYP9K1. The increase in resistance prompted a reversion to bendiocarb IRS in 2016 to avoid a resurgence of malaria, in line with the national Malaria Control Program plan.

Implicating Cryptic and Novel Anophelines as Malaria Vectors in Africa

Entomological indices and bionomic descriptions of malaria vectors are essential to accurately describe and understand malaria transmission and for the design and evaluation of appropriate control interventions. In order to correctly assign spatio-temporal distributions, behaviors and responses to interventions to particular anopheline species, identification of mosquitoes must be accurately made. This paper reviews the current methods and their limitations in correctly identifying anopheline mosquitoes in sub-Saharan Africa, and highlights the importance of molecular methods to discriminate cryptic species and identify lesser known anophelines. The increasing number of reports of Plasmodium infections in assumed "minor", non-vector, and cryptic and novel species is reviewed. Their importance in terms of evading current control and elimination strategies and therefore maintaining malaria transmission is emphasized.

Desiccation tolerance in Anopheles coluzzii: the effects of spiracle size and cuticular hydrocarbons

The African malaria mosquitoes Anopheles gambiae and Anopheles coluzzii range over forests and arid areas, where they withstand dry spells and months-long dry seasons, suggesting variation in their desiccation tolerance. We subjected a laboratory colony (G3) and wild Sahelian mosquitoes during the rainy and dry seasons to desiccation assays. The thoracic spiracles and amount and composition of cuticular hydrocarbons (CHCs) of individual mosquitoes were measured to determine the effects of these traits on desiccation tolerance. The relative humidity of the assay, body water available, rate of water loss and water content at death accounted for 88% of the variation in desiccation tolerance. Spiracle size did not affect the rate of water loss or desiccation tolerance of the colony mosquitoes, as was the case for the total CHCs. However, six CHCs accounted for 71% of the variation in desiccation tolerance and three accounted for 72% of the variation in the rate of water loss. Wild A. coluzzii exhibited elevated desiccation tolerance during the dry season. During that time, relative thorax and spiracle sizes were smaller than during the rainy season. A smaller spiracle size appeared to increase A. coluzzii's desiccation tolerance, but was not statistically significant. Seasonal changes in CHC composition were detected in Sahelian A. coluzzii Stepwise regression models suggested the effect of particular CHCs on desiccation tolerance. In conclusion, the combination of particular CHCs along with the total amount of CHCs is a primary mechanism conferring desiccation tolerance in A. coluzzii, while variation in spiracle size might be a secondary mechanism.

The reproductive tracts of two malaria vectors are populated by a core microbiome and by gender- and swarm-enriched microbial biomarkers

Microbes play key roles in shaping the physiology of insects and can influence behavior, reproduction and susceptibility to pathogens. In Sub-Saharan Africa, two major malaria vectors, Anopheles gambiae and An. coluzzii, breed in distinct larval habitats characterized by different microorganisms that might affect their adult physiology and possibly Plasmodium transmission. We analyzed the reproductive microbiomes of male and female An. gambiae and An. coluzzii couples collected from natural mating swarms in Burkina Faso. 16S rRNA sequencing on dissected tissues revealed that the reproductive tracts harbor a complex microbiome characterized by a large core group of bacteria shared by both species and all reproductive tissues. Interestingly, we detected a significant enrichment of several gender-associated microbial biomarkers in specific tissues, and surprisingly, similar classes of bacteria in males captured from one mating swarm, suggesting that these males originated from the same larval breeding site. Finally, we identified several endosymbiotic bacteria, including Spiroplasma, which have the ability to manipulate insect reproductive success. Our study provides a comprehensive analysis of the reproductive microbiome of important human disease vectors, and identifies a panel of core and endosymbiotic bacteria that can be potentially exploited to interfere with the transmission of malaria parasites by the Anopheles mosquito.

Does extreme asymmetric dominance promote hybridization between Anopheles coluzzii and Anopheles gambiae s.s. in seasonal malaria mosquito communities of West Africa?

Background

Anopheles gambiae s.s. and An. coluzzii are two of the most important malaria vector species in sub-Saharan Africa. These recently-diverged sibling species do not exhibit intrinsic post-zygotic barriers to reproduction and are thought to be separated by strong assortative mating combined with selection against hybrids. At present, little is known about the ecological conditions that determine hybridization and introgression between these cryptic taxa.

Methods

Swarm segregation and assortative mating were studied in Western Burkina Faso in the villages of Vallée du Kou (VK7) and Soumousso which differed in terms of which sibling species was much rarer than the other, and in Bana where both occurred in similar proportions. Swarms and pairs in copula were collected and genotyped, the proportion of intra and interspecific mating determined, and interspecific sperm transfer checked genetically. Females were collected through larval and adult indoor collections and genotyped or sexed-and-genotyped via a novel multiplex PCR.

Results

A total of 3,687 males and 220 females were collected and genotyped from 109 swarms. Only 3 swarms were composed of males from both species, and these were from the village of VK7 where An. gambiae s.s. was comparatively rare. Mixed-species pairs captured in copula were only detected in that area and made for 3.62 % and 100 % of mating pairs involving An. coluzzii and An. gambiae s.s. individuals, respectively. The high An. gambiae s.s. cross-mating rate was mirrored by high rates of hybridizations estimated from female larvae and adults indoor collections. This contrasted with Soumousso, where despite being much less common than An. gambiae s.s., An. coluzzii males did not form mixed swarms, females were not found in interspecific swarms or copula and hybridization rates were low in both sibling species.

Conclusions

These data suggest that ecological conditions leading to rare An. gambiae s.s. in populations dominated by An. coluzzii may promote a breakdown of spatial swarm segregation and assortative mating between the two species. The lower overall hybridization rates observed at the larval and adult indoor stages compared to cross-mating rates support the idea that post-mating selection processes acting against hybrids may occur mostly prior to and/or at the early larval stages.

Heterosis Increases Fertility, Fecundity, and Survival of Laboratory-Produced F1 Hybrid Males of the Malaria Mosquito Anopheles coluzzii

The success of vector control strategies aiming to decrease disease transmission via the release of sterile or genetically-modified male mosquitoes critically depends on mating between laboratory-reared males and wild females. Unfortunately, mosquito colonization, laboratory rearing, and genetic manipulations can all negatively affect male competitiveness. Heterosis is commonly used to produce domestic animals with enhanced vigor and homogenous genetic background and could therefore potentially improve the mating performance of mass-reared male mosquitoes. Here, we produced enhanced hybrid males of the malaria mosquito Anopheles coluzzii by crossing two strains colonized >35 and 8 years ago. We compared the amount of sperm and mating plug proteins they transferred to females, as well as their insemination rate, reproductive success and longevity under various experimental conditions. Across experiments, widespread adaptations to laboratory mating were detected in the older strain. In large-group mating experiments, no overall hybrid advantage in insemination rates and the amount of sperm and accessory gland proteins transferred to females was detected. Despite higher sperm activity, hybrid males did not appear more fecund. However, individual-male mating and laboratory-swarm experiments revealed that hybrid males, while inseminating fewer females than older inbred males, were significantly more fertile, producing larger mating plugs and drastically increasing female fecundity. Heterotic males also showed increased longevity. These results validate the use of heterosis for creating hybrid males with improved fitness from long-established inbred laboratory strains. Therefore, this simple approach could facilitate disease control strategies based on male mosquito releases with important ultimate benefits to human health.

A multi-locus approach resolves the phylogenetic relationships of the Simulium asakoae and Simulium ceylonicum species groups in Malaysia: evidence for distinct evolutionary lineages

A multi-locus approach was used to examine the DNA sequences of 10 nominal species of blackfly in the Simulium subgenus Gomphostilbia (Diptera: Simuliidae) in Malaysia. Molecular data were acquired from partial DNA sequences of the mitochondria-encoded cytochrome c oxidase subunit I (COI), 12S rRNA and 16S rRNA genes, and the nuclear-encoded 18S rRNA and 28S rRNA genes. No single gene, nor the concatenated gene set, resolved all species or all relationships. However, all morphologically established species were supported by at least one gene. The multi-locus sequence analysis revealed two distinct evolutionary lineages, conforming to the morphotaxonomically recognized Simulium asakoae and Simulium ceylonicum species groups.

Stimulating Anopheles gambiae swarms in the laboratory: application for behavioural and fitness studies

Background

Male Anopheles mosquitoes that swarm rely in part on features of the environment including visual stimuli to locate swarms. Swarming is believed to be the primary behaviour during which mating occurs in the field, but is not a common behaviour in the laboratory. Features that stimulate male Anopheles gambiae G3 strain swarming were created in novel large indoor cages.

Methods

The following visual features were tested in all combinations to determine which were important for swarm formation. Large cages and fading ceiling lights at dusk alone did not stimulate swarming while a dark foreground and contrasting illuminated background with a contrasting landmark stimulated and localized swarm formation during artificial twilight. Given the need to test transgenic strains in as natural a setting as possible, in this study it was investigated whether induced swarm behaviour and cage size would affect relative mating performance of wild-type and transgenic β2Ppo1 and β2Ppo2 A. gambiae sexually sterile males.

Results

Even using a mosquito colony that has been in laboratory culture for 39 years, swarming behaviour was induced by this novel arrangement. The presence of swarming stimuli was associated with an increase in insemination frequency from 74.3 to 97.7% in large cages. Transgenic males showed a lower competitiveness in large cages compared to small cages regardless of the presence of swarming stimuli.

Conclusions

The results of the present study are discussed in view of the progressive evaluation of genetically modified A. gambiae strains and the potential applications of reproducing swarms in controlled conditions to dissect the mating behaviour of this species and the mechanisms controlling it.

Paratransgenesis: a promising new strategy for mosquito vector control

The three main mosquito genera, Anopheles, Aedes and Culex, transmit respectively malaria, dengue and lymphatic filariasis. Current mosquito control strategies have proved unsuccessful, and there still is a substantial number of morbidity and mortality from these diseases. Genetic control methods have now arisen as promising alternative strategies, based on two approaches: the replacement of a vector population by disease-refractory mosquitoes and the release of mosquitoes carrying a lethal gene to suppress target populations. However, substantial hurdles and limitations need to be overcome if these methods are to be used successfully, the most significant being that a transgenic mosquito strain is required for every target species, making genetically modified mosquito strategies inviable when there are multiple vector mosquitoes in the same area. Genetically modified bacteria capable of colonizing a wide range of mosquito species may be a solution to this problem and another option for the control of these diseases. In the paratransgenic approach, symbiotic bacteria are genetically modified and reintroduced in mosquitoes, where they express effector molecules. For this approach to be used in practice, however, requires a better understanding of mosquito microbiota and that symbiotic bacteria and effector molecules be identified. Paratransgenesis could prove very useful in mosquito species that are inherently difficult to transform or in sibling species complexes. In this approach, a genetic modified bacteria can act by: (a) causing pathogenic effects in the host; (b) interfering with the host's reproduction; (c) reducing the vector's competence; and (d) interfering with oogenesis and embryogenesis. It is a much more flexible and adaptable approach than the use of genetically modified mosquitoes because effector molecules and symbiotic bacteria can be replaced if they do not achieve the desired result. Paratransgenesis may therefore become an important integrated pest management tool for mosquito control.

The Seasonality and Ecology of the Anopheles gambiae complex (Dipetra: Culicidae) in Liberia Using Molecular Identification

Members of the Anopheles gambiae sensu lato (Giles) complex define a group of seven morphologically indistinguishable species, including the principal malaria vectors in Sub-Saharan Africa. Members of this complex differ in behavior and ability to transmit malaria; hence, precise identification of member species is critical to monitoring and evaluating malaria threat levels. We collected mosquitoes from five counties in Liberia every other month from May 2011 until May 2012, using various trapping techniques. A. gambiae complex members were identified using molecular techniques based on differences in the ribosomal DNA (rDNA) region between species and the molecular forms (S and M) of A. gambiae sensu stricto (s.s) specimens. In total, 1,696 A. gambiae mosquitoes were collected and identified. DNA was extracted from legs of each specimen with species identification determined by multiplex polymerase chain reaction using specific primers. The molecular forms (M or S) of A. gambiae s.s were determined by restriction fragment length polymorphism. Bivariate and multivariate logistic regression models identified environmental variables associated with genomic differentiation. Our results indicate widespread occurrence of A. gambiae s.s., the principal malaria vector in the complex, although two Anopheles melas Theobald/A. merus Donitz mosquitoes were detected. We found 72.6, 25.5, and 1.9% of A. gambiae s.s specimens were S, M, and hybrid forms, respectively. Statistical analysis indicates that the S form was more likely to be found in rural areas during rainy seasons and indoor catchments. This information will enhance vector control efforts in Liberia.

Experimental swap of Anopheles gambiae's assortative mating preferences demonstrates key role of X-chromosome divergence island in incipient sympatric speciation

Although many theoretical models of sympatric speciation propose that genes responsible for assortative mating amongst incipient species should be associated with genomic regions protected from recombination, there are few data to support this theory. The malaria mosquito, Anopheles gambiae, is known for its sympatric cryptic species maintained by pre-mating reproductive isolation and its putative genomic islands of speciation, and is therefore an ideal model system for studying the genomic signature associated with incipient sympatric speciation. Here we selectively introgressed the island of divergence located in the pericentric region of the X chromosome of An. gambiae s.s. into its sister taxon An. coluzzii through 5 generations of backcrossing followed by two generations of crosses within the introgressed strains that resulted in An. coluzzii-like recombinant strains fixed for the M and S marker in the X chromosome island. The mating preference of recombinant strains was then tested by giving virgin recombinant individuals a choice of mates with X-islands matching and non-matching their own island type. We show through genetic analyses of transferred sperm that recombinant females consistently mated with matching island-type males thereby associating assortative mating genes with the X-island of divergence. Furthermore, full-genome sequencing confirmed that protein-coding differences between recombinant strains were limited to the experimentally swapped pericentromeric region. Finally, targeted-genome comparisons showed that a number of these unique differences were conserved in sympatric field populations, thereby revealing candidate speciation genes. The functional demonstration of a close association between speciation genes and the X-island of differentiation lends unprecedented support to island-of-speciation models of sympatric speciation facilitated by pericentric recombination suppression.

Anopheles gambiae is the most important vector of malaria in Africa. This species is undergoing speciation and a number of subpopulations have been identified which can produce viable hybrid offspring but are reproductively isolated through assortative mating and ecological adaptation. This complex structure provides an ideal system for studying the unique genetic and behavioural processes required for speciation. Anopheles gambiae’s subpopulations differ genetically in limited regions of their genomes called islands of speciation. Theoretical studies predict that these islands, characterized by restricted genetic rearrangements, may protect genes of assortative mating between emerging species, and are fundamental to the speciation process. We set out to test this prediction by performing complex genetic crosses between the sister species Anopheles coluzzii and Anopheles gambiae s.s. and creating recombinant strains differing only at their X-chromosome island of speciation. We show through behavioural studies that recombinant females consistently mated with matching island-type males thereby associating assortative mating genes with the X-island of divergence. By sequencing the genetic code of the recombinant strains and natural populations, we could confirm these findings and identify candidate assortative mating genes. These findings suggest an important role of divergence islands for the genetic and behavioural processes associated with speciation.

Remarkable diversity of intron-1 of the para voltage-gated sodium channel gene in an Anopheles gambiae/Anopheles coluzzii hybrid zone

Background

Genomic differentiation between Anopheles gambiae and Anopheles coluzzii - the major malaria vectors in sub-Saharan Africa - is localized into large “islands” toward the centromeres of chromosome-X and the two autosomes. Linkage disequilibrium between these genomic islands was first detected between species-specific polymorphisms within ribosomal DNA genes (IGS-rDNA) on the X-chromosome and a single variant at position 702 of intron 1 (Int-1⁷⁰²) of the para Voltage-Gated Sodium Channel (VGSC) gene on chromosome arm 2 L. Intron-1 sequence data from West and Central Africa revealed two clearly distinct and species-specific haplogroups, each characterized by very low polymorphism, which has been attributed to a selective sweep. The aim of this study was to analyse Int-1 sequence diversity in A. gambiae and A. coluzzii populations from the Far-West of their range, in order to assess whether this selective-sweep signature could persist in a zone of high interspecific hybridization.

Methods

A 531 bp region of VGSC Int-1 was sequenced in 21 A. coluzzii, 31 A. gambiae, and 12 hybrids from The Gambia and Guinea Bissau, located within the Far-West geographical region, and in 53 A. gambiae s.l. samples from the rest of the range.

Results

Far-West samples exhibit dramatic Int-1 polymorphism, far higher within each country than observed throughout the rest of the species range. Moreover, patterning of haplotypes within A. coluzzii confirms previous evidence of a macro-geographic subdivision into a West and a Central African genetic cluster, and reveals a possible genetic distinction of A. coluzzii populations from the Far-West.

Conclusions

The results suggest a relaxation of selective pressures acting across the VGSC gene region in the hybrid zone. Genetic differentiation in the Far-West could be attributable to a founder effect within A. coluzzii, with subsequent extensive gene flow with secondarily-colonizing A. gambiae, potentially yielding a novel insight on the dynamic processes impacting genetic divergence of these key malaria vectors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-014-0522-1) contains supplementary material, which is available to authorized users.

Reproductive isolation among sympatric molecular forms of An. gambiae from inland areas of south-eastern Senegal

The Anopheles gambiae species complex includes at least seven morphologically indistinguishable species, one of which, Anopheles gambiae sensu stricto, is the primary mosquito vector responsible for the transmission of malaria across sub-Saharan Africa. Sympatric ecological diversification of An. gambiae s.s. is in progress within this complex, leading to the emergence of at least two incipient species (the M and S molecular forms now recognized as good species and named An. coluzzii and An. gambiae respectively) that show heterogeneous levels of divergence in most parts of Africa. However, this process seems to have broken down in coastal areas of West Africa at the extreme edge of the distribution. We undertook a longitudinal study to describe An. gambiae s.s. populations collected from two inland transects with different ecological characteristics in south-eastern Senegal. Analysis of samples collected from 20 sites across these two transects showed the M and S molecular forms coexisted at almost all sampled sites. Overall, similar hybridization rates (2.16% and 1.86%) were recorded in the two transects; sites with relatively high frequencies of M/S hybrids (up to 7%) were clustered toward the north-western part of both transects, often near urban settings. Estimated inbreeding indices for this putative speciation event varied spatially (range: 0.52-1), with hybridization rates being generally lower than expected under panmictic conditions. Such observations suggest substantial reproductive isolation between the M and S molecular forms, and further support the ongoing process of speciation in these inland areas. According to a recent reclassification of the An. gambiae complex, the M and S molecular forms from this zone correspond to An. coluzzii and An. gambiae, respectively. There is considerable evidence that these molecular forms differ in their behavioural and ecological characteristics. Detailed study of these characteristics will allow the development and implementation of better insect control strategies for combating malaria.

Evidence of natural Wolbachia infections in field populations of Anopheles gambiae

Wolbachia are maternally transmitted intracellular bacteria that invade insect populations by manipulating their reproduction and immunity and thus limiting the spread of numerous human pathogens. Experimental Wolbachia infections can reduce Plasmodium numbers in Anopheles mosquitoes in the laboratory, however, natural Wolbachia infections in field anophelines have never been reported. Here we show evidence of Wolbachia infections in Anopheles gambiae in Burkina Faso, West Africa. Sequencing of the 16S rRNA gene identified Wolbachia sequences in both female and male germlines across two seasons, and determined that these sequences are vertically transmitted from mother to offspring. Whole-genome sequencing of positive samples suggests that the genetic material identified in An. gambiae belongs to a novel Wolbachia strain, related to but distinct from strains infecting other arthropods. The evidence of Wolbachia infections in natural Anopheles populations promotes further investigations on the possible use of natural Wolbachia-Anopheles associations to limit malaria transmission.

Characterization of swarming and mating behaviour between Anopheles coluzzii and Anopheles melas in a sympatry area of Benin

The swarm structure of two sibling species, Anopheles gambiae coluzzii and Anopheles melas, was characterize to explore the ecological and environmental parameters associated with the formation of swarms and their spatial distribution. Swarms and breeding sites were searched and sampled between January and December 2010, and larval and adult samples were identified by PCR. During the dry season, 456 swarms of An. gambiae s.l. were sampled from 38 swarm sites yielding a total of 23,274 males and 76 females. Of these 38 swarming sites, 18 were composed exclusively of An. gambiae coluzzii and 20 exclusively of An. melas, presenting clear evidence of reproductive swarm segregation. The species makeup of couples sampled from swarms also demonstrated assortative mating. The swarm site localization was close to human dwellings in the case of the An. gambiae coluzzii and on salt production sites for An. melas. At the peak of the rainy season, swarms of An. melas were absent. These findings offer evidence that the ecological speciation of these two sibling species of mosquitoes is associated with spatial swarm segregation and assortative mating, providing strong support for the hypothesis that mate recognition is currently maintaining adaptive differentiation and promoting ecological speciation. Further studies on the swarming and mating systems of An. gambiae, with the prospect of producing a predictive model of swarm distribution, are needed to inform any future efforts to implement strategies based on the use of GMM or SIT.

Swarming behaviour in natural populations of Anopheles gambiae and An. coluzzii: review of 4 years survey in rural areas of sympatry, Burkina Faso (West Africa)

The swarming behaviour of natural populations of Anopheles gambiae and An. coluzzii (formerly known as An. gambiae S and M forms, respectively) were investigated through longitudinal surveys conducted between July 2006 and October 2009 in two rural areas of south-western Burkina Faso where these forms are sympatric. In both sites, the majority of swarms were recorded above visual markers localised among houses. In Soumousso, a wooded area of savannah, 108 pairs caught in copula from 205 swarms were sampled; in VK7, a rice growing area, 491 couples from 250 swarms were sampled. If segregated swarms were the norm in both sites, many visual markers were shared by the two forms of An. gambiae. Furthermore, mixed swarms were collected annually in frequencies varying from one site to another, though no mixed inseminations were recorded, corroborating the low hybrid rate previously reported in the field. The occurrence of inter-specific mate-recognition mechanisms, which allow individuals to avoid hybridisation, is discussed.

Differential effects of inbreeding and selection on male reproductive phenotype associated with the colonization and laboratory maintenance of Anopheles gambiae

Background

Effective mating between laboratory-reared males and wild females is paramount to the success of vector control strategies aiming to decrease disease transmission via the release of sterile or genetically modified male mosquitoes. However mosquito colonization and laboratory maintenance have the potential to negatively affect male genotypic and phenotypic quality through inbreeding and selection, which in turn can decrease male mating competitiveness in the field. To date, very little is known about the impact of those evolutionary forces on the reproductive biology of mosquito colonies and how they ultimately affect male reproductive fitness.

Methods

Here several male reproductive physiological traits likely to be affected by inbreeding and selection following colonization and laboratory rearing were examined. Sperm length, and accessory gland and testes size were compared in male progeny from field-collected females and laboratory strains of Anopheles gambiae sensu stricto colonized from one to over 25 years ago. These traits were also compared in the parental and sequentially derived, genetically modified strains produced using a two-phase genetic transformation system. Finally, genetic crosses were performed between strains in order to distinguish the effects of inbreeding and selection on reproductive traits.

Results

Sperm length was found to steadily decrease with the age of mosquito colonies but was recovered in refreshed strains and crosses between inbred strains therefore incriminating inbreeding costs. In contrast, testes size progressively increased with colony age, whilst accessory gland size quickly decreased in males from colonies of all ages. The lack of heterosis in response to crossing and strain refreshing in the latter two reproductive traits suggests selection for insectary conditions.

Conclusions

These results show that inbreeding and selection differentially affect reproductive traits in laboratory strains overtime and that heterotic ‘supermales’ could be used to rescue some male reproductive characteristics. Further experiments are needed to establish the exact relationship between sperm length, accessory gland and testes size, and male reproductive success in the laboratory and field settings.

Seasonal variation in spatial distributions of Anopheles gambiae in a Sahelian village: evidence for aestivation

Changes in spatial distribution of mosquitoes over time in a Sahelian village were studied to understand the sources of the mosquitoes during the dry season when no larval sites are found. At that time, the sources of Anopheles gambiae Giles may be local shelters used by aestivating mosquitoes or migrants from distant populations. The mosquito distribution was more aggregated during the dry season, when few houses had densities 7- to 24-fold higher than expected. The high-density houses during the dry season differed from those of the wet season. Most high-density houses during the dry season changed between years, yet their vicinity was rather stable. Scan statistics confirmed the presence of one or two adjacent hotspots in the dry season, usually found on one edge of the village. These hotspots shifted between the early and late dry season. During the wet season, the hotspots were relatively stable near the main larval site. The locations of the hotspots in the wet season and early and late dry season were similar between years. Season-specific, stable, and focal hotspots are inconsistent with the predictions based on the arrival of migrants from distant localities during the dry season, but are consistent with the predictions based on local shelters used by aestivating mosquitoes. Targeting hotspots in Sahelian villages for vector control may not be effective because the degree of aggregation is moderate, the hotspots are not easily predicted, and they are not the sources of the population. However, targeting the dry-season shelters may be highly cost-effective, once they can be identified and predicted.

Freshwater biodiversity and aquatic insect diversification

Inland waters cover less than 1% of Earth's surface but harbor more than 6% of all insect species: Nearly 100,000 species from 12 orders spend one or more life stages in freshwater. Little is known about how this remarkable diversity arose, although allopatric speciation and ecological adaptation are thought to be primary mechanisms. Freshwater habitats are highly susceptible to environmental change and exhibit marked ecological gradients. Standing waters appear to harbor more dispersive species than running waters, but there is little understanding of how this fundamental ecological difference has affected diversification. In contrast to the lack of evolutionary studies, the ecology and habitat preferences of aquatic insects have been intensively studied, in part because of their widespread use as bioindicators. The combination of phylogenetics with the extensive ecological data provides a promising avenue for future research, making aquatic insects highly suitable models for the study of ecological diversification.

Genetic and environmental factors associated with laboratory rearing affect survival and assortative mating but not overall mating success in Anopheles gambiae sensu stricto

Anopheles gambiae sensu stricto, the main vector of malaria in Africa, is characterized by its vast geographical range and complex population structure. Assortative mating amongst the reproductively isolated cryptic forms that co-occur in many areas poses unique challenges for programs aiming to decrease malaria incidence via the release of sterile or genetically-modified mosquitoes. Importantly, whether laboratory-rearing affects the ability of An. gambiae individuals of a given cryptic taxa to successfully mate with individuals of their own form in field conditions is still unknown and yet crucial for mosquito-releases. Here, the independent effects of genetic and environmental factors associated with laboratory rearing on male and female survival, mating success and assortative mating were evaluated in the Mopti form of An. gambiae over 2010 and 2011. In semi-field enclosures experiments and despite strong variation between years, the overall survival and mating success of male and female progeny from a laboratory strain was not found to be significantly lower than those of the progeny of field females from the same population. Adult progeny from field-caught females reared at the larval stage in the laboratory and from laboratory females reared outdoors exhibited a significant decrease in survival but not in mating success. Importantly, laboratory individuals reared as larvae indoors were unable to mate assortatively as adults, whilst field progeny reared either outdoors or in the laboratory, as well as laboratory progeny reared outdoors all mated significantly assortatively. These results highlight the importance of genetic and environment interactions for the development of An. gambiae's full mating behavioral repertoire and the challenges this creates for mosquito rearing and release-based control strategies.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

Antennal transcriptome profiles of anopheline mosquitoes reveal human host olfactory specialization in Anopheles gambiae

BACKGROUND

Two sibling members of the Anopheles gambiae species complex display notable differences in female blood meal preferences. An. gambiae s.s. has a well-documented preference for feeding upon human hosts, whereas An. quadriannulatus feeds on vertebrate/mammalian hosts, with only opportunistic feeding upon humans. Because mosquito host-seeking behaviors are largely driven by the sensory modality of olfaction, we hypothesized that hallmarks of these divergent host seeking phenotypes will be in evidence within the transcriptome profiles of the antennae, the mosquito's principal chemosensory appendage.

RESULTS

To test this hypothesis, we have sequenced antennal mRNA of non-bloodfed females from each species and observed a number of distinct quantitative and qualitative differences in their chemosensory gene repertoires. In both species, these gene families show higher rates of sequence polymorphisms than the overall rates in their respective transcriptomes, with potentially important divergences between the two species. Moreover, quantitative differences in odorant receptor transcript abundances have been used to model potential distinctions in volatile odor receptivity between the two sibling species of anophelines.

CONCLUSION

This analysis suggests that the anthropophagic behavior of An. gambiae s.s. reflects the differential distribution of olfactory receptors in the antenna, likely resulting from a co-option and refinement of molecular components common to both species. This study improves our understanding of the molecular evolution of chemoreceptors in closely related anophelines and suggests possible mechanisms that underlie the behavioral distinctions in host seeking that, in part, account for the differential vectorial capacity of these mosquitoes.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

BACKGROUND

Anopheles (Kerteszia) cruzii (Diptera: Culicidae) is a primary vector of human and simian malaria parasites in southern and southeastern Brazil. Earlier studies using chromosome inversions, isoenzymes and a number of molecular markers have suggested that An. cruzii is a species complex.

RESULTS

In this study, a multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was carried out to investigate in more detail the genetic differentiation between the An. cruzii populations from Florianópolis-Santa Catarina (southern Brazil) and Itatiaia-Rio de Janeiro States (southeastern Brazil). The analyses were performed first comparing Florianópolis and Itatiaia, and then comparing the two putative sympatric incipient species from Itatiaia (Itatiaia A and Itatiaia B). The analysis revealed high FST values between Florianópolis and Itatiaia (considering Itatiaia A and B together) and also between the sympatric Itatiaia A and Itatiaia B, irrespective of their function. Also, using the IM program, no strong indication of migration was found between Florianópolis and Itatiaia (considering Itatiaia A and B together) using all loci together, but between Itatiaia A and Itatiaia B, the results show evidence of migration only in the direction of Itatiaia B.

CONCLUSIONS

The results of the multilocus analysis indicate that Florianópolis and Itatiaia represent different species of the An. cruzii complex that diverged around 0.6 Mya, and also that the Itatiaia sample is composed of two sympatric incipient species A and B, which diverged around 0.2 Mya. Asymmetric introgression was found between the latter two species despite strong divergence in some loci.

Insecticide resistance monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance

Insecticide resistance in the malaria vector Anopheles gambiae s.l. (Diptera: Culicidae) threatens insecticide-based control efforts, necessitating regular monitoring. We assessed resistance in field-collected An. gambiae s.l. from Jinja, Uganda using World Health Organization (WHO) bioassays. Only An. gambiae s.s. and An. arabiensis (≈70%) were present. Female An. gambiae exhibited extremely high pyrethroid resistance (permethrin LT50 > 2 h; deltamethrin LT50 > 5 h). Female An. arabiensis were resistant to permethrin and exhibited reduced susceptibility to deltamethrin. However, while An. gambiae were DDT resistant, An. arabiensis were fully susceptible. Both species were fully susceptible to bendiocarb and fenitrothion. Kdr 1014S has increased rapidly in the Jinja population of An. gambiae s.s. and now approaches fixation (≈95%), consistent with insecticide-mediated selection, but is currently at a low frequency in An. arabiensis (0.07%). Kdr 1014F was also at a low frequency in An. gambiae. These frequencies preclude adequately-powered tests for an association with phenotypic resistance. PBO synergist bioassays resulted in near complete recovery of pyrethroid susceptibility suggesting involvement of CYP450s in resistance. A small number (0.22%) of An. gambiae s.s. ×An. arabiensis hybrids were found, suggesting the possibility of introgression of resistance alleles between species. The high levels of pyrethroid resistance encountered in Jinja threaten to reduce the efficacy of vector control programmes which rely on pyrethroid-impregnated bednets or indoor spraying of pyrethroids.

Assortative mating in mixed swarms of the mosquito Anopheles gambiae s.s. M and S molecular forms, in Burkina Faso, West Africa

The molecular form composition of Anopheles gambiae Giles s.s. (Diptera: Culicidae) mating swarms and the associated mating pairs (copulae) were investigated during two rainy seasons (July to October, 2005 and July to November, 2006) in the villages of Soumousso and Vallée du Kou (VK7). Although the habitats of these villages differ markedly, sympatric populations of M and S molecular forms of An. gambiae s.s. occur in both places periodically. The main aim was to assess the degree to which these molecular forms mate assortatively. In Soumousso, a wooded savannah habitat, the majority of swarm samples consisted of only S-form males (21/28), although a few M-form males were found in mixed M- and S-form swarms. In VK7, a rice growing area, the majority of swarm samples consisted of only M-form males (38/62), until October and November 2006, when there were nearly as many mixed-form as single-form swarms. Overall, ∼60% of M- and S-form swarms were temporally or spatially segregated; the two forms were effectively prevented from encountering each other. Of the remaining 40% of swarms, however, only about half were single-form and the rest were mixed-form. Of the 33 copulae collected from mixed-form swarms, only four were mixed-form pairs, significantly fewer than expected by random pairing between forms (χ(2) = 10.34, d.f. = 2, P < 0.01). Finally, all specimens of inseminated females were of the same form as the sperm contained within their spermatheca (n = 91), even for the four mixed-form copulae. These findings indicate that assortative mating occurs within mixed-form swarms, mediated most probably by close-range mate recognition cues.

Chromosomal inversions among insecticide-resistant strains of Anopheles stephensi Liston, a malaria mosquito

Polytene chromosomes were prepared from the ovarian nurse cells of semi-gravid females of ten insecticide-resistant strains of Anopheles stephensi. Altogether, 16 heterozygous paracentric inversions, namely b/+ (11D-16C) in alphamethrin; i/+ (14B-18A) and h/+ (27B-28A) in DDT; j/+ (14A-16B) in chlorpyrifos; k/+ (11D-16B) in cyfluthrin; l/+ (11A-16C) in deltamethrin; m/+ (14B-15C) and e/+ (32A-33B) in bifenthrin; n/+ (12D-14B), f/+ (33A-36A) and g/+ (33C-34A) in propoxur; o/+ (11A-12D), h/+ (37A-37C) and i/+ (31C-32C) in temephos; d/+ (33D-35C) in carbofuran and a/+ (41C-43B) in neem strains, were reported. No inversions were observed in X chromosome so far. The frequency of inversions in different insecticides was found to be highest in the 2R arm, followed by the 3R arm. Such inversions were not reported in the corresponding susceptible strains or in the parental stocks.

Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow

The primary Afrotropical malaria mosquito vector Anopheles gambiae sensu stricto has a complex population structure. In west Africa, this species is split into two molecular forms and displays local and regional variation in chromosomal arrangements and behaviors. To investigate patterns of macrogeographic population substructure, 25 An. gambiae samples from 12 African countries were genotyped at 13 microsatellite loci. This analysis detected the presence of additional population structuring, with the M-form being subdivided into distinct west, central, and southern African genetic clusters. These clusters are coincident with the central African rainforest belt and northern and southern savannah biomes, which suggests restrictions to gene flow associated with the transition between these biomes. By contrast, geographically patterned population substructure appears much weaker within the S-form.

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.

Identification of cryptic Anopheles mosquito species by molecular protein profiling

Vector control is the mainstay of malaria control programmes. Successful vector control profoundly relies on accurate information on the target mosquito populations in order to choose the most appropriate intervention for a given mosquito species and to monitor its impact. An impediment to identify mosquito species is the existence of morphologically identical sibling species that play different roles in the transmission of pathogens and parasites. Currently PCR diagnostics are used to distinguish between sibling species. PCR based methods are, however, expensive, time-consuming and their development requires a priori DNA sequence information. Here, we evaluated an inexpensive molecular proteomics approach for Anopheles species: matrix assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). MALDI-TOF MS is a well developed protein profiling tool for the identification of microorganisms but so far has received little attention as a diagnostic tool in entomology. We measured MS spectra from specimens of 32 laboratory colonies and 2 field populations representing 12 Anopheles species including the A. gambiae species complex. An important step in the study was the advancement and implementation of a bioinformatics approach improving the resolution over previously applied cluster analysis. Borrowing tools for linear discriminant analysis from genomics, MALDI-TOF MS accurately identified taxonomically closely related mosquito species, including the separation between the M and S molecular forms of A. gambiae sensu stricto. The approach also classifies specimens from different laboratory colonies; hence proving also very promising for its use in colony authentication as part of quality assurance in laboratory studies. While being exceptionally accurate and robust, MALDI-TOF MS has several advantages over other typing methods, including simple sample preparation and short processing time. As the method does not require DNA sequence information, data can also be reviewed at any later stage for diagnostic or functional patterns without the need for re-designing and re-processing biological material.

Physiological correlates of ecological divergence along an urbanization gradient: differential tolerance to ammonia among molecular forms of the malaria mosquito Anopheles gambiae

BACKGROUND

Limitations in the ability of organisms to tolerate environmental stressors affect their fundamental ecological niche and constrain their distribution to specific habitats. Evolution of tolerance, therefore, can engender ecological niche dynamics. Forest populations of the afro-tropical malaria mosquito Anopheles gambiae have been shown to adapt to historically unsuitable larval habitats polluted with decaying organic matter that are found in densely populated urban agglomerates of Cameroon. This process has resulted in niche expansion from rural to urban environments that is associated with cryptic speciation and ecological divergence of two evolutionarily significant units within this taxon, the molecular forms M and S, among which reproductive isolation is significant but still incomplete. Habitat segregation between the two forms results in a mosaic distribution of clinally parapatric patches, with the M form predominating in the centre of urban agglomerates and the S form in the surrounding rural localities. We hypothesized that development of tolerance to nitrogenous pollutants derived from the decomposition of organic matter, among which ammonia is the most toxic to aquatic organisms, may affect this pattern of distribution and process of niche expansion by the M form.

RESULTS

Acute toxicity bioassays indicated that populations of the two molecular forms occurring at the extremes of an urbanization gradient in Yaounde, the capital of Cameroon, differed in their response to ammonia. The regression lines best describing the dose-mortality profile differed in the scale of the explanatory variable (ammonia concentration log-transformed for the S form and linear for the M form), and in slope (steeper for the S form and shallower for the M form). These features reflected differences in the frequency distribution of individual tolerance thresholds in the two populations as assessed by probit analysis, with the M form exhibiting a greater mean and variance compared to the S form.

CONCLUSIONS

In agreement with expectations based on the pattern of habitat partitioning and exposure to ammonia in larval habitats in Yaounde, the M form showed greater tolerance to ammonia compared to the S form. This trait may be part of the physiological machinery allowing forest populations of the M form to colonize polluted larval habitats, which is at the heart of its niche expansion in densely populated human settlements in Cameroon.

Can field-based mosquito feeding assays be used for evaluating transmission-blocking interventions?

A recent meta-analysis of mosquito feeding assays to determine the Plasmodium falciparum transmission potential of naturally infected gametocyte carriers highlighted considerable variation in transmission efficiency between assay methodologies and between laboratories. This begs the question as to whether mosquito feeding assays should be used for the evaluation of transmission-reducing interventions in the field and whether these field-based mosquito assays are currently standardized sufficiently to enable accurate evaluations. Here, we address biological and methodological reasons for the observed variations, discuss whether these preclude the use of field-based mosquito feeding assays in field evaluations of transmission-blocking interventions, and propose how we can maximize the precision of estimates. Altogether, we underscore the significant advantages of field-based mosquito feeding assays in basic malaria research and field trials.

Physiology and development of the M and S molecular forms of Anopheles gambiae in Burkina Faso (West Africa)

In West Africa, M and S molecular forms of Anopheles gambiae sensu stricto (Diptera: Culicidae) Giles, frequently occur together, although with different population bionomics. The S form typically breeds in rain-dependant water collections and is present during the rainy season only whereas the M form can thrive all year long in areas with permanent breeding opportunities. In the present study, we explored physiological and developmental trade-offs at play in laboratory colonies and field populations of the M and S forms that originated from an area of sympatry in Burkina Faso, where M and S larvae exhibit such habitat segregation. In the laboratory, larvae of the M form developed slower than the S form (mean values 9.51 and 8.85 days, respectively, Wilcoxon's test, P < 0.001). Although wing length and dry weight at emergence showed large variations, M females were on average 8% heavier than S females of similar wing length. Higher nutritional reserves (proteins and lipids) in teneral adults explained part of this weight difference, reflecting a better ability of the M form to garner resources at the larval stage. Furthermore, a higher rate of ovarian maturation was observed in the M form after a single bloodmeal. The relevance of these findings for parasite transmission is discussed.

Evidence for population-specific positive selection on immune genes of Anopheles gambiae

Host-pathogen interactions can be powerful drivers of adaptive evolution, shaping the patterns of molecular variation at the genes involved. In this study, we sequenced alleles from 28 immune-related loci in wild samples of multiple genetic subpopulations of the African malaria mosquito Anopheles gambiae, obtaining unprecedented sample sizes and providing the first opportunity to contrast patterns of molecular evolution at immune-related loci in the recently discovered GOUNDRY population to those of the indoor-resting M and S molecular forms. In contrast to previous studies that focused on immune genes identified in laboratory studies, we centered our analysis on genes that fall within a quantitative trait locus associated with resistance to Plasmodium falciparum in natural populations of A. gambiae. Analyses of haplotypic and genetic diversity at these 28 loci revealed striking differences among populations in levels of genetic diversity and allele frequencies in coding sequence. Putative signals of positive selection were identified at 11 loci, but only one was shared by two subgroups of A. gambiae. Striking patterns of linkage disequilibrium were observed at several loci. We discuss these results with respect to ecological differences among these strata as well as potential implications for disease transmission.

Anthropogenic habitat disturbance and ecological divergence between incipient species of the malaria mosquito Anopheles gambiae

BACKGROUND

Anthropogenic habitat disturbance is a prime cause in the current trend of the Earth's reduction in biodiversity. Here we show that the human footprint on the Central African rainforest, which is resulting in deforestation and growth of densely populated urban agglomerates, is associated to ecological divergence and cryptic speciation leading to adaptive radiation within the major malaria mosquito Anopheles gambiae.

METHODOLOGY/PRINCIPAL FINDINGS

In southern Cameroon, the frequency of two molecular forms--M and S--among which reproductive isolation is strong but still incomplete, was correlated to an index of urbanisation extracted from remotely sensed data, expressed as the proportion of built-up surface in each sampling unit. The two forms markedly segregated along an urbanisation gradient forming a bimodal cline of ∼6-km width: the S form was exclusive to the rural habitat, whereas only the M form was present in the core of densely urbanised settings, co-occurring at times in the same polluted larval habitats of the southern house mosquito Culex quinquefasciatus--a species association that was not historically recorded before.

CONCLUSIONS/SIGNIFICANCE

Our results indicate that when humans create novel habitats and ecological heterogeneities, they can provide evolutionary opportunities for rapid adaptive niche shifts associated with lineage divergence, whose consequences upon malaria transmission might be significant.