Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s

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.

A new species in the major malaria vector complex sheds light on reticulated species evolution

Complexes of closely related species provide key insights into the rapid and independent evolution of adaptive traits. Here, we described and studied Anopheles fontenillei sp.n., a new species in the Anopheles gambiae complex that we recently discovered in the forested areas of Gabon, Central Africa. Our analysis placed the new taxon in the phylogenetic tree of the An. gambiae complex, revealing important introgression events with other members of the complex. Particularly, we detected recent introgression, with Anopheles gambiae and Anopheles coluzzii, of genes directly involved in vectorial capacity. Moreover, genome analysis of the new species allowed us to clarify the evolutionary history of the 3La inversion. Overall, An. fontenillei sp.n. analysis improved our understanding of the relationship between species within the An. gambiae complex, and provided insight into the evolution of vectorial capacity traits that are relevant for the successful control of malaria in Africa.

Phylogenetic signal from rearrangements in 18 Anopheles species by joint scaffolding extant and ancestral genomes

Background

Genomes rearrangements carry valuable information for phylogenetic inference or the elucidation of molecular mechanisms of adaptation. However, the detection of genome rearrangements is often hampered by current deficiencies in data and methods: Genomes obtained from short sequence reads have generally very fragmented assemblies, and comparing multiple gene orders generally leads to computationally intractable algorithmic questions.

Results

We present a computational method, ADseq, which, by combining ancestral gene order reconstruction, comparative scaffolding and de novo scaffolding methods, overcomes these two caveats. ADseq provides simultaneously improved assemblies and ancestral genomes, with statistical supports on all local features. Compared to previous comparative methods, it runs in polynomial time, it samples solutions in a probabilistic space, and it can handle a significantly larger gene complement from the considered extant genomes, with complex histories including gene duplications and losses. We use ADseq to provide improved assemblies and a genome history made of duplications, losses, gene translocations, rearrangements, of 18 complete Anopheles genomes, including several important malaria vectors. We also provide additional support for a differentiated mode of evolution of the sex chromosome and of the autosomes in these mosquito genomes.

Conclusions

We demonstrate the method’s ability to improve extant assemblies accurately through a procedure simulating realistic assembly fragmentation. We study a debated issue regarding the phylogeny of the Gambiae complex group of Anopheles genomes in the light of the evolution of chromosomal rearrangements, suggesting that the phylogenetic signal they carry can differ from the phylogenetic signal carried by gene sequences, more prone to introgression.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4466-7) contains supplementary material, which is available to authorized users.

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

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

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.

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.

Diversity, differentiation, and linkage disequilibrium: prospects for association mapping in the malaria vector Anopheles arabiensis

Association mapping is a widely applied method for elucidating the genetic basis of phenotypic traits. However, factors such as linkage disequilibrium and levels of genetic diversity influence the power and resolution of this approach. Moreover, the presence of population subdivision among samples can result in spurious associations if not accounted for. As such, it is useful to have a detailed understanding of these factors before conducting association mapping experiments. Here we conducted whole-genome sequencing on 24 specimens of the malaria mosquito vector, Anopheles arabiensis, to further understanding of patterns of genetic diversity, population subdivision and linkage disequilibrium in this species. We found high levels of genetic diversity within the An. arabiensis genome, with ~800,000 high-confidence, single- nucleotide polymorphisms detected. However, levels of nucleotide diversity varied significantly both within and between chromosomes. We observed lower diversity on the X chromosome, within some inversions, and near centromeres. Population structure was absent at the local scale (Kilombero Valley, Tanzania) but detected between distant populations (Cameroon vs. Tanzania) where differentiation was largely restricted to certain autosomal chromosomal inversions such as 2Rb. Overall, linkage disequilibrium within An. arabiensis decayed very rapidly (within 200 bp) across all chromosomes. However, elevated linkage disequilibrium was observed within some inversions, suggesting that recombination is reduced in those regions. The overall low levels of linkage disequilibrium suggests that association studies in this taxon will be very challenging for all but variants of large effect, and will require large sample sizes.

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.

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.

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.

Genetic differentiation between sandfly populations of Phlebotomus chinensis and Phlebotomus sichuanensis (Diptera: Psychodidae) in China inferred by microsatellites

Background

Phlebotomus chinensis is a primary vector of visceral leishmaniasis; it occurs in various biotopes with a large geographical distribution, ranging from Yangtze River to northeast China. Phlebotomus sichuanensis, a species closely related to P. chinensis in high altitude regions, has a long term disputation on its taxonomic status. Both species occur in the current epidemic regions and are responsible for the transmission of leishmaniasis. Population genetic analysis will help to understand the population structure and infer the relationship for morphologically indistinguishable cryptic species. In this study, microsatellite markers were used for studying the genetic differentiation between P. chinensis and P. sichuanensis.

Methods

Sandflies were collected in 6 representative localities in China in 2005-2009. Ten microsatellite loci were used to estimate population genetic diversity. The intra-population genetic diversity, genetic differentiation and effective population size were estimated.

Results

All 10 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy-Weinberg disequilibrium was found in 23 out of 60 (38.33%) comparisons associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed three clusters. The cluster I included almost all specimens in the sample SCD collected at high altitude habitats in Sichuan. The other two clusters were shared by the remaining 5 populations, SCJ in Sichuan, GSZ in Gansu, SXL and SXX in Shaanxi and HNS in Henan. The diversity among these 5 populations was low (F_ST = -0.003-0.090) and no isolation by distance was detected. AMOVA analysis suggested that the variations were largely derived from individuals within populations and among individuals. Consistently, the analysis of ribosomal DNA second internal transcribed spacer (ITS2) sequence uncovered three types of variants, which corresponded with the three gene pools revealed by microsatellites.

Conclusions

The data suggested that the SCD population carried a distinct gene pool, which was differentiated from the other populations. The high altitude ecological habitats, distinctive ITS2 and herein divergence inferred by microsatellite loci support the species status of P. sichuanensis. The P. chinensis populations did not have a significant divergence from each another.

Application of a qPCR assay in the investigation of susceptibility to malaria infection of the M and S molecular forms of An. gambiae s.s. in Cameroon

Plasmodium falciparum is the causative agent of malaria, a disease that kills almost one million persons each year, mainly in sub-Saharan Africa. P. falciparum is transmitted to the human host by the bite of an Anopheles female mosquito, and Anopheles gambiae sensus stricto is the most tremendous malaria vector in Africa, widespread throughout the afro-tropical belt. An. gambiae s.s. is subdivided into two distinct molecular forms, namely M and S forms. The two molecular forms are morphologically identical but they are distinct genetically, and differ by their distribution and their ecological preferences. The epidemiological importance of the two molecular forms in malaria transmission has been poorly investigated so far and gave distinct results in different areas. We have developed a real-time quantitative PCR (qPCR) assay, and used it to detect P. falciparum at the oocyst stage in wild An. gambiae s.s. mosquitoes experimentally infected with natural isolates of parasites. Mosquitoes were collected at immature stages in sympatric and allopatric breeding sites and further infected at the adult stage. We next measured the infection prevalence and intensity in female mosquitoes using the qPCR assay and correlated the infection success with the mosquito molecular forms. Our results revealed different prevalence of infection between the M and S molecular forms of An. gambiae s.s. in Cameroon, for both sympatric and allopatric populations of mosquitoes. However, no difference in the infection intensity was observed. Thus, the distribution of the molecular forms of An. gambiae s.s. may impact on the malaria epidemiology, and it will be important to monitor the efficiency of malaria control interventions on the two M and S forms.

Diversity and transmission competence in lymphatic filariasis vectors in West Africa, and the implications for accelerated elimination of Anopheles-transmitted filariasis

Lymphatic Filariasis (LF) is targeted for elimination by the Global Programme for the Elimination of Lymphatic Filariasis (GPELF). The strategy adopted is based on the density dependent phenomenon of Facilitation, which hypothesizes that in an area where the vector species transmitting Wuchereria bancrofti are Anopheles mosquitoes, it is feasible to eliminate LF using Mass Drug Administration (MDA) because of the inability of Anopheles species to transmit low-density microfilaraemia. Even though earlier studies have shown Anopheles species can exhibit the process of Facilitation in West Africa, observations point towards the process of Limitation in certain areas, in which case vector control is recommended. Studies on Anopheles species in West Africa have also shown genetic differentiation, cryptic taxa and speciation, insecticide resistance and the existence of molecular and chromosomal forms, all of which could influence the vectorial capacity of the mosquitoes and ultimately the elimination goal. This paper outlines the uniqueness of LF vectors in West Africa and the challenges it poses to the 2020 elimination goal, based on the current MDA strategies.

Dry season reproductive depression of Anopheles gambiae in the Sahel

The African malaria mosquito, Anopheles gambiae, is widespread south of the Sahara including in dry savannahs and semi-arid environments where no surface water exists for several months a year. Adults of the M form of An. gambiae persist through the long dry season, when no surface waters are available, by increasing their maximal survival from 4 weeks to 7 months. Dry season diapause (aestivation) presumably underlies this extended survival. Diapause in adult insects is intrinsically linked to depressed reproduction. To determine if reproduction of the Sahelian M form is depressed during the dry season, we assessed seasonal changes in oviposition, egg batch size, and egg development, as well as insemination rate and blood feeding in wild caught mosquitoes. Results from xeric Sahelian and riparian populations were compared. Oviposition response in the Sahelian M form dropped from 70% during the wet season to 20% during the dry season while the mean egg batch size among those that laid eggs fell from 173 to 101. Correspondingly, the fraction of females that exhibited gonotrophic dissociation increased over the dry season from 5% to 45%, while a similar fraction of the population retained developed eggs despite having access to water. This depression in reproduction the Sahelian M form was not caused by a reduced insemination rate. Seasonal variation in these reproductive parameters of the riparian M form population was less extreme and the duration of reproductive depression was shorter. Blood feeding responses did not change with the season in either population. Depressed reproduction during the dry season in the Sahelian M form of An. gambiae provides additional evidence for aestivation and illuminates the physiological processes involved. The differences between the Sahelian and riparian population suggest an adaptive cline in aestivation phenotypes between populations only 130 km apart.

Evidence for X-linked introgression between molecular forms of Anopheles gambiae from Angola

The M and S molecular forms of the African malaria vector Anopheles gambiae (Diptera: Culicidae) are morphologically identical incipient species in which reproductive isolation is incomplete, enabling low-level gene flow between forms. In an attempt to find differences between the M and S forms, sequence variation was studied at loci along the X chromosome in adult female An. gambiae from Angola. A high proportion of M form specimens from Angola (79% of the 456 X chromosomes sampled) were found to contain a 16-bp insertion in intron 4 of the X-linked GPRCCK1 locus, relative to the AgamP3 release of the An. gambiae PEST genome sequence. The insertion was in Hardy-Weinberg equilibrium in Angolan M form populations. The same insertion was found in all S form specimens examined, regardless of where in Africa they were sampled, but was absent from a sample of M form specimens collected in Ghana, Bioko and Mali. In M form specimens from Angola, there was an association between alleles at the GPRCCK1 locus and those at a microsatellite locus, AGXH678, close to the centromere of the X chromosome, with significant linkage disequilibrium between loci separated by 0.472 Mbp (P < 0.033). We show that the insertion results from introgression from the S form into the M form, rather than from the retention of an ancestral character. Gene flow from the S to M form could allow genes of adaptive value to be transferred, including those conferring insecticide resistance and others influencing ecology and behaviour, and thus malaria transmission and control. We discuss factors that may have led to this introgression event.

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.

Kdr-based insecticide resistance in Anopheles gambiae s.s populations in

BACKGROUND

The spread of insecticide resistance in the malaria mosquito, Anopheles gambiae is a serious threat for current vector control strategies which rely on the use of insecticides. Two mutations at position 1014 of the S6 transmembrane segment of domain II in the voltage gated sodium channel, known as kdr (knockdown resistance) mutations leading to a change of a Leucine to a Phenylalanine (L1014F) or to a Serine (L1014S) confer resistance to DDT and pyrethroid insecticides in the insect. This paper presents the current distribution of the kdr alleles in wild Anopheles gambiae populations in Cameroon.

RESULTS

A total of 1,405 anopheline mosquitoes were collected from 21 localities throughout Cameroon and identified as An. gambiae (N = 1,248; 88.8%), An. arabiensis (N = 120; 8.5%) and An. melas (N = 37; 2.6%). Both kdr alleles 1014F and 1014S were identified in the M and S molecular forms of An. gambiae s.s. The frequency of the 1014F allele ranged from 1.7 to 18% in the M-form, and from 2 to 90% in the S-form. The 1014S allele ranged from 3-15% in the S-form and in the M-form its value was below 3%. Some specimens were found to carry both resistant kdr alleles.

CONCLUSION

This study provides an updated distribution map of the kdr alleles in wild An. gambiae populations in Cameroon. The co-occurrence of both alleles in malaria mosquito vectors in diverse ecological zones of the country may be critical for the planning and implementation of malaria vector control interventions based on IRS and ITNs, as currently ongoing in Cameroon.

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.

Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiaemolecular form identification

BACKGROUND

Anopheles gambiae M and S molecular forms, the major malaria vectors in the Afro-tropical region, are ongoing a process of ecological diversification and adaptive lineage splitting, which is affecting malaria transmission and vector control strategies in West Africa. These two incipient species are defined on the basis of single nucleotide differences in the IGS and ITS regions of multicopy rDNA located on the X-chromosome. A number of PCR and PCR-RFLP approaches based on form-specific SNPs in the IGS region are used for M and S identification. Moreover, a PCR-method to detect the M-specific insertion of a short interspersed transposable element (SINE200) has recently been introduced as an alternative identification approach. However, a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification was never carried out to evaluate whether they could be used interchangeably, as commonly assumed.

RESULTS

The genotyping of more than 400 A. gambiae specimens from nine African countries, and the sequencing of the IGS-amplicon of 115 of them, highlighted discrepancies among results obtained by the different approaches due to different kinds of biases, which may result in an overestimation of MS putative hybrids, as follows: i) incorrect match of M and S specific primers used in the allele specific-PCR approach; ii) presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches; iii) incomplete cleavage during the restriction reactions; iv) presence of different copy numbers of M and S-specific IGS-arrays in single individuals in areas of secondary contact between the two forms.

CONCLUSIONS

The results reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae M and S forms are not fully interchangeable as usually assumed, and highlight limits of the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. These limits are discussed and operational suggestions on the choice of the most convenient method for large-scale M- and S-form identification are provided, also taking into consideration technical aspects related to the epidemiological characteristics of different study areas.

The "far-west" of Anopheles gambiae molecular forms

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

Malaria vectors in the Republic of Benin: distribution of species and molecular forms of the Anopheles gambiae complex

Members of the Anopheles gambiae complex are among the best malaria vectors in the world, but their vectorial capacities vary between species and populations. A large-scale sampling of An. gambiae sensu lato was carried out in 2006 and 2007 in various bioclimatic areas of Benin (West Africa). The objective of this study was to collate data on the relative frequencies of species and forms within the An. gambiae complex and to produce a map of their spatial distribution. Sampling took place at 30 sites and 2122 females were analyzed. Two species were identified through molecular methods. The overall collection showed a preponderance of An. gambiae s.s., but unexpectedly, An. arabiensis was reported in the coastal-Guinean bioclimatic area characterized by a mean annual rainfall of >1500 mm where only An. gambiae s.s. was reported previously. Our study of Benin indicates that An. arabiensis would be adapted not only to the urban areas but also to the rural humid regions. Among 1717 An. gambiae s.s., 26.5% were of the M form and 73.3% were S form. Few hybrid specimens between the M and S forms were observed (0.2%). Only the spatial distribution of the M form appears to be mainly a function of bioclimatic area. Factors that influence the distribution of these malaria vectors are discussed. This study underlines the need of further investigations of biological, ecological, and behavioral traits of these species and forms to better appreciate their vectorial capacities. Acquisition of entomological field data appears essential to better estimate the stratification of malaria risk and help improve malaria vector control interventions.

Authentication scheme for routine verification of genetically similar laboratory colonies: a trial with Anopheles gambiae

BACKGROUND

When rearing morphologically indistinguishable laboratory strains concurrently, the threat of unintentional genetic contamination is constant. Avoidance of accidental mixing of strains is difficult due to the use of common equipment, technician error, or the possibility of self relocation by adult mosquitoes ("free fliers"). In many cases, laboratory strains are difficult to distinguish because of morphological and genetic similarity, especially when laboratory colonies are isolates of certain traits from the same parental strain, such as eye color mutants, individuals with certain chromosomal arrangements or high levels of insecticide resistance. Thus, proving genetic integrity could seem incredibly time-consuming or impossible. On the other hand, lacking proof of genetically isolated laboratory strains could question the validity of research results.

RESULTS

We present a method for establishing authentication matrices to routinely distinguish and confirm that laboratory strains have not become physically or genetically mixed through contamination events in the laboratory. We show a specific example with application to Anopheles gambiae sensu stricto strains at the Malaria Research and Reference Reagent Resource Center. This authentication matrix is essentially a series of tests yielding a strain-specific combination of results.

CONCLUSION

These matrix-based methodologies are useful for several mosquito and insect populations but must be specifically tailored and altered for each laboratory based on the potential contaminants available at any given time. The desired resulting authentication plan would utilize the least amount of routine effort possible while ensuring the integrity of the strains.

Spatial swarm segregation and reproductive isolation between the molecular forms of Anopheles gambiae

Anopheles gambiae, the major malaria vector in Africa, can be divided into two subgroups based on genetic and ecological criteria. These two subgroups, termed the M and S molecular forms, are believed to be incipient species. Although they display differences in the ecological niches they occupy in the field, they are often sympatric and readily hybridize in the laboratory to produce viable and fertile offspring. Evidence for assortative mating in the field was recently reported, but the underlying mechanisms awaited discovery. We studied swarming behaviour of the molecular forms and investigated the role of swarm segregation in mediating assortative mating. Molecular identification of 1145 males collected from 68 swarms in Donéguébougou, Mali, over 2 years revealed a strict pattern of spatial segregation, resulting in almost exclusively monotypic swarms with respect to molecular form. We found evidence of clustering of swarms composed of individuals of a single molecular form within the village. Tethered M and S females were introduced into natural swarms of the M form to verify the existence of possible mate recognition operating within-swarm. Both M and S females were inseminated regardless of their form under these conditions, suggesting no within-mate recognition. We argue that our results provide evidence that swarm spatial segregation strongly contributes to reproductive isolation between the molecular forms in Mali. However this does not exclude the possibility of additional mate recognition operating across the range distribution of the forms. We discuss the importance of spatial segregation in the context of possible geographic variation in mechanisms of reproductive isolation.

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.

High levels of hybridization between molecular forms of Anopheles gambiae from Guinea Bissau

In the malaria vector Anopheles gambiae Giles sensu stricto, two molecular forms denoted M and S are considered units of incipient speciation within this species. Very low hybrid frequencies and significant genetic differentiation have been found in sympatric M- and S-form populations. We studied the molecular form composition and the degree of genetic differentiation at 15 microsatellites in two samples of An. gambiae collected in two consecutive years from Bissau, Guinea Bissau. High frequencies of M/S hybrids (19-24%) were found in this area. Coincidently, very low levels of genetic differentiation were detected between forms when analysis involved microsatellites mapped at chromosome-3 (mean Fst, 0.000-0.002). The single exception was the X-linked AGXH678, for which high differentiation was measured (Fst, 0.158-0.301). This locus maps near the centromere of chromosome X, a low recombination region in which selection is likely to promote divergence between M and S forms. These results strongly suggest that the degree of isolation between M and S forms, considered the units of incipient speciation within An. gambiae, is not homogenous throughout the species distribution range.

Exploring the origin and degree of genetic isolation of Anopheles gambiae from the islands of São Tomé and Príncipe, potential sites for testing transgenic-based vector control

The evolutionary processes at play between island and mainland populations of the malaria mosquito vector Anopheles gambiae sensu stricto are of great interest as islands may be suitable sites for preliminary application of transgenic-based vector control strategies. São Tomé and Príncipe, located off the West African coast, have received such attention in recent years. This study investigates the degree of isolation of An. gambiae s.s. populations between these islands and the mainland based on mitochondrial and ribosomal DNA molecular data. We identify possible continental localities from which these island populations derived. For these purposes, we used F_ST values, haplotype networks, and nested clade analysis to estimate migration rates and patterns. Haplotypes from both markers are geographically widespread across the African continent. Results indicate that the populations from São Tomé and Príncipe are relatively isolated from continental African populations, suggesting they are promising sites for test releases of transgenic individuals. These island populations are possibly derived from two separate continental migrations. This result is discussed in the context of the history of the African slave trade with respect to São Tomé and Príncipe.

Dynamics of multiple insecticide resistance in the malaria vector Anopheles gambiae in a rice growing area in South-Western Burkina Faso

Background

Insecticide resistance of the main malaria vector, Anopheles gambiae, has been reported in south-western Burkina Faso, West Africa. Cross-resistance to DDT and pyrethroids was conferred by alterations at site of action in the sodium channel, the Leu-Phe kdr mutation; resistance to organophosphates and carbamates resulted from a single point mutation in the oxyanion hole of the acetylcholinesterase enzyme designed as ace-1^R.

Methods

An entomological survey was carried out during the rainy season of 2005 at Vallée du Kou, a rice growing area in south-western Burkina Faso. At the Vallée du Kou, both insecticide resistance mechanisms have been previously described in the M and S molecular forms of An. gambiae. This survey aimed i) to update the temporal dynamics and the circumsporozoite infection rate of the two molecular forms M and S of An. gambiae ii) to update the frequency of the Leu-Phe kdr mutation within these forms and finally iii) to investigate the occurrence of the ace-1^R mutation.

Mosquitoes collected by indoor residual collection and by human landing catches were counted and morphologically identified. Species and molecular forms of An. gambiae, ace-1^R and Leu-Phe kdr mutations were determined using PCR techniques. The presence of the circumsporozoite protein of Plasmodium falciparum was determined using ELISA.

Results

Anopheles gambiae populations were dominated by the M form. However the S form occurred in relative important proportion towards the end of the rainy season with a maximum peak in October at 51%. Sporozoite rates were similar in both forms. The frequency of the Leu-Phe kdr mutation in the S form reached a fixation level while it is still spreading in the M form. Furthermore, the ace-1^R mutation prevailed predominately in the S form and has just started spreading in the M form. The two mutations occurred concomitantly both in M and S populations.

Conclusion

These results showed that the Vallée du Kou, a rice growing area formerly occupied mainly by M susceptible populations, is progressively colonized by S resistant populations living in sympatry with the former. As a result, the distribution pattern of insecticide resistance mutations shows the occurrence of both resistance mechanisms concomitantly in the same populations. The impact of multiple resistance mechanisms in M and S populations of An. gambiae on vector control measures against malaria transmission, such as insecticide-treated nets (ITNs) and indoor residual spraying (IRS), in this area is discussed.

Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms

Background

SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms.

Methods

A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed.

Results

Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted evolution and/or mixtures of M and S IGS-sequences among the arrays of single chromatids.

Conclusion

The approach utilized allowed to develop new easy-to-use co-dominant markers for the analysis of genetic differentiation between M and S-forms and opens new perspectives in the study of the speciation process ongoing within A. gambiae.

The molecular forms of Anopheles gambiae: a phenotypic perspective

The African malaria mosquito Anopheles gambiae is undergoing speciation, being split into the M and S molecular forms. Speciation is the main process promoting biological diversity, thus, new vector species might complicate disease transmission. Genetic differentiation between the molecular forms has been extensively studied, but phenotypic differences between them, the evolutionary forces that generated divergence, and the mechanisms that maintain their genetic isolation have only recently been addressed. Here, we review recent studies suggesting that selection mediated by larval predation and competition promoted divergence between temporary and permanent freshwater habitats. These differences explain the sharp discontinuity in distribution of the molecular forms between rice fields and surrounding savanna, but they can also explain the concurrent cline between humid and arid environments due to the dependence on permanent habitats in the latter. Although less pronounced, differences in adult body size, reproductive output, and longevity also suggest that the molecular forms have adapted to distinct niches. Reproductive isolation between the molecular forms is achieved by spatial swarm segregation, although within-swarm mate recognition appears to play a role in certain locations. The implications of these results to disease transmission and control are discussed and many of the gaps in our understanding are highlighted.

Evidence of introgression of the ace-1(R) mutation and of the ace-1 duplication in West African Anopheles gambiae s. s

BACKGROUND

The role of inter-specific hybridisation is of particular importance in mosquito disease vectors for predicting the evolution of insecticide resistance. Two molecular forms of Anopheles gambiae s.s., currently recognized as S and M taxa, are considered to be incipient sibling species. Hybrid scarcity in the field was suggested that differentiation of M and S taxa is maintained by limited or absent gene flow. However, recent studies have revealed shared polymorphisms within the M and S forms, and a better understanding of the occurrence of gene flow is needed. One such shared polymorphism is the G119S mutation in the ace-1 gene (which is responsible for insecticide resistance); this mutation has been described in both the M and S forms of A. gambiae s.s.

METHODS AND RESULTS

To establish whether the G119S mutation has arisen independently in each form or by genetic introgression, we analysed coding and non-coding sequences of ace-1 alleles in M and S mosquitoes from representative field populations. Our data revealed many polymorphic sites shared by S and M forms, but no diversity was associated with the G119S mutation. These results indicate that the G119S mutation was a unique event and that genetic introgression explains the observed distribution of the G119S mutation within the two forms. However, it was impossible to determine from our data whether the mutation occurred first in the S form or in the M form. Unexpectedly, sequence analysis of some resistant individuals revealed a duplication of the ace-1 gene that was observed in both A. gambiae s.s. M and S forms. Again, the distribution of this duplication in the two forms most likely occurred through introgression.

CONCLUSIONS

These results highlight the need for more research to understand the forces driving the evolution of insecticide resistance in malaria vectors and to regularly monitor resistance in mosquito populations of Africa.

High genetic differentiation between the M and S molecular forms of Anopheles gambiae in Africa

Background

Anopheles gambiae, a major vector of malaria, is widely distributed throughout sub-Saharan Africa. In an attempt to eliminate infective mosquitoes, researchers are trying to develop transgenic strains that are refractory to the Plasmodium parasite. Before any release of transgenic mosquitoes can be envisaged, we need an accurate picture of the differentiation between the two molecular forms of An. gambiae, termed M and S, which are of uncertain taxonomic status.

Methodology/Principal Findings

Insertion patterns of three transposable elements (TEs) were determined in populations from Benin, Burkina Faso, Cameroon, Ghana, Ivory Coast, Madagascar, Mali, Mozambique, Niger, and Tanzania, using Transposon Display, a TE-anchored strategy based on Amplified Fragment Length Polymorphism. The results reveal a clear differentiation between the M and S forms, whatever their geographical origin, suggesting an incipient speciation process.

Conclusions/Significance

Any attempt to control the transmission of malaria by An. gambiae using either conventional or novel technologies must take the M/S genetic differentiation into account. In addition, we localized three TE insertion sites that were present either in every individual or at a high frequency in the M molecular form. These sites were found to be located outside the chromosomal regions that are suspected of involvement in the speciation event between the two forms. This suggests that these chromosomal regions are either larger than previously thought, or there are additional differentiated genomic regions interspersed with undifferentiated regions.

Evidence for divergent selection between the molecular forms of Anopheles gambiae: role of predation

BACKGROUND

The molecular forms of Anopheles gambiae are undergoing speciation. They are characterized by a strong assortative mating and they display partial habitat segregation. The M form is mostly found in flooded/irrigated areas whereas the S form dominates in the surrounding areas, but the ecological factors that shape this habitat segregation are not known. Resource competition has been demonstrated between species undergoing divergent selection, but resource competition is not the only factor that can lead to divergence.

RESULTS

In a field experiment using transplantation of first instar larvae, we evaluated the role of larval predators in mediating habitat segregation between the forms. We found a significant difference in the ability of the molecular forms to exploit the different larval sites conditioned on the presence of predators. In absence of predation, the molecular forms outcompeted each other in their respective natural habitats however, the developmental success of the M form was significantly higher than that of the S form in both habitats under predator pressure.

CONCLUSION

Our results provide the first empirical evidence for specific adaptive differences between the molecular forms and stress the role of larval predation as one of the mechanisms contributing to their divergence.

X and Y chromosome inheritance and mixtures of rDNA intergenic spacer regions in Anopheles gambiae

We observed Anopheles gambiae sensu stricto stocks that contained both Mopti (M) and Savanna (S) rDNA intergenic spacers (IGS). ASEMBO1 male IGS sequences consistently had a mixture. A diagnostic M and S Hha I restriction enzyme site in these fragments was concordant with two SNPs associated with M and S. Standard M and S stocks demonstrated X-chromosome-only inheritance of the rDNA form, but the ASEMBO1 males showed X and Y chromosome linkage of mixed rDNA. The metaphase Y chromosomes of ASEMBO1 contained a significantly larger amount of DNA relative to the X than a standard S stock. Analysis of wild A. gambiae males from the putative location of origin of the ASEMBO1 stock did not detect the same pattern of polymorphism observed in the laboratory stock but did detect heterogeneous arrays including some missing a diagnostic Hha I restriction site. These results demonstrate that M and S IGS types can occur within the rDNA arrays of a chromatid in laboratory A. gambiae stocks, and some A. gambiae s.s. have rDNA on the Y chromosome.

Mosaic genome architecture of the Anopheles gambiae species complex

Background

Attempts over the last three decades to reconstruct the phylogenetic history of the Anopheles gambiae species complex have been important for developing better strategies to control malaria transmission.

Methodology

We used fingerprint genotyping data from 414 field-collected female mosquitoes at 42 microsatellite loci to infer the evolutionary relationships of four species in the A. gambiae complex, the two major malaria vectors A. gambiae sensu stricto (A. gambiae s.s.) and A. arabiensis, as well as two minor vectors, A. merus and A. melas.

Principal Findings

We identify six taxonomic units, including a clear separation of West and East Africa A. gambiae s.s. S molecular forms. We show that the phylogenetic relationships vary widely between different genomic regions, thus demonstrating the mosaic nature of the genome of these species. The two major malaria vectors are closely related and closer to A. merus than to A. melas at the genome-wide level, which is also true if only autosomes are considered. However, within the Xag inversion region of the X chromosome, the M and two S molecular forms are most similar to A. merus. Near the X centromere, outside the Xag region, the two S forms are highly dissimilar to the other taxa. Furthermore, our data suggest that the centromeric region of chromosome 3 is a strong discriminator between the major and minor malaria vectors.

Conclusions

Although further studies are needed to elucidate the basis of the phylogenetic variation among the different regions of the genome, the preponderance of sympatric admixtures among taxa strongly favor introgression of different genomic regions between species, rather than lineage sorting of ancestral polymorphism, as a possible mechanism.

Genetic population structure of Anopheles gambiae in Equatorial Guinea

BACKGROUND

Patterns of genetic structure among mosquito vector populations in islands have received particular attention as these are considered potentially suitable sites for experimental trials on transgenic-based malaria control strategies. In this study, levels of genetic differentiation have been estimated between populations of Anopheles gambiae s.s. from the islands of Bioko and Annobón, and from continental Equatorial Guinea (EG) and Gabon.

METHODS

Genotyping of 11 microsatellite loci located in chromosome 3 was performed in three island samples (two in Bioko and one in Annobón) and three mainland samples (two in EG and one in Gabon). Four samples belonged to the M molecular form and two to the S-form. Microsatellite data was used to estimate genetic diversity parameters, perform demographic equilibrium tests and analyse population differentiation.

RESULTS

High levels of genetic differentiation were found between the more geographically remote island of Annobón and the continent, contrasting with the shallow differentiation between Bioko island, closest to mainland, and continental localities. In Bioko, differentiation between M and S forms was higher than that observed between island and mainland samples of the same molecular form.

CONCLUSION

The observed patterns of population structure seem to be governed by the presence of both physical (the ocean) and biological (the M-S form discontinuity) barriers to gene flow. The significant degree of genetic isolation between M and S forms detected by microsatellite loci located outside the "genomic islands" of speciation identified in A. gambiae s.s. further supports the hypothesis of on-going incipient speciation within this species. The implications of these findings regarding vector control strategies are discussed.

Locus- and Population-Specific Selection and Differentiation between Incipient Species of Anopheles gambiae

Anopheles gambiae, the primary mosquito vector of malaria in sub-Saharan Africa, is divided into 2 sympatric incipient species known as M form and S form. Recent genomic analysis of each form revealed that differentiation between forms is clustered into 3 unlinked regions of the genome. Here, we expand the investigation of these "genomic islands of speciation" to multiple populations, including all of the genes across one of the islands. Differentiation between the M and S forms in 2 of the islands is complete across all individuals in all populations, confirming that the M and S forms are reproductively isolated taxa. Differentiation at the third island (on chromosome 2R) is limited to Cameroon populations. There is reduced variation in the M form in Cameroon at this location and increased divergence to the outgroup Anopheles arabiensis, supporting an association of adaptation with reproductive isolation.

Comparative analysis of epicuticular lipid profiles of sympatric and allopatric field populations of Anopheles gambiae s.s. molecular forms and An. arabiensis from Burkina Faso (West Africa)

We analysed by gas chromatography-mass spectrometry (GC-MS) and Gas Chromatography-Flame Ionization Detector (GC-FID) the epicuticular lipid profiles of field females of the major Afro-tropical malaria vector, Anopheles gambiae. The samples were collected in three villages in Burkina Faso (West Africa), where An. gambiae M and S molecular forms and An. arabiensis live sympatrically. The aim was to compare the cuticular hydrocarbon (CHC) composition of individual field specimens of these three taxa, to highlight possible differences among them. All the samples analysed by GC-MS (55 individuals and eight pools) were characterized by the same 48 CHCs and 10 oxygenated compounds. The 19 most abundant CHCs were quantified in 174 specimens by GC-FID: quantitative intra-taxon differences were found between allopatric populations of both An. arabiensis and S-form. Inter-taxa quantitative differences in the relative abundances of some hydrocarbons between pairs of sympatric taxa were also found, which appear to be mainly linked to local situations, with the possible exception of diMeC(35) between An. arabiensis and S-form. Moreover, MeC(29) shows some degree of differentiation between S- and M-form in all three villages. Possible causes of these differences are discussed.

Insertion polymorphism of transposable elements and population structure of Anopheles gambiae M and S molecular forms in Cameroon

The insertion polymorphism of five transposable element (TE) families was studied by Southern blots in several populations of the M and S molecular forms of the mosquito Anopheles gambiae sensu stricto from southern Cameroon. We showed that the mean TE insertion site number and the within-population insertion site polymorphism globally differed between the M and S molecular forms. The comparison of the TE insertion profiles of the populations revealed a significant differentiation between these two molecular forms (0.163 < Phi(ST) < 0.371). We cloned several insertions of a non-LTR retrotransposon (Aara8) that were fixed in one form and absent in the other one. The only insertion that could be clearly located on a chromosome arm mapped to cytological division 6 of chromosome X, confirming the importance of this region in the ongoing speciation between the M and S molecular forms.

IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis

BACKGROUND

Polymerase chain reactions to distinguish single-nucleotide polymorphisms are commonly used for mosquito identification and identifying insecticide resistance alleles. However, the existing methods used for primer design often result in analyses that are not robust or require additional steps.

METHODS

Utilizing oligonucleotides that are unique in having an intentional mismatch to both templates three bases from the SNP at the 3-prime end, three new PCR assays that distinguish SNP targets using standard gel electrophoresis of undigested DNA fragments were developed and tested. These were applied to: (1) an alternative ribosomal DNA PCR assay to distinguish five members of the Anopheles gambiae complex; (2) detection of the Mopti and Savanna rDNA types; and (3) an assay to distinguish resistance to dieldrin (Rdl) alleles in Anopheles arabiensis.

RESULTS

Reproducible specific amplification of the target alleles was observed in all three assays. The results were consistent with existing analyses but proved simpler and the results more distinct in our hands.

CONCLUSION

The simplicity and effectiveness of the method should be utilized in these and other PCR analyses to increase their specificity and simplicity. These results have the potential to be extended not only to mosquito analyses but also to parasite and human polymorphisms.

Reduced recombination rate and genetic differentiation between the M and S forms of Anopheles gambiae s.s

Genetic differentiation between the largely sympatric molecular forms M and S of Anopheles gambiae appears mostly limited to division 6 and part of division 5 of the X chromosome. This region is adjacent to the centromere and includes the rDNA that was used to define these forms. This localized differentiation between populations that experience gene flow strongly suggests that this region contains genes responsible for reproductive isolation. Regions adjacent to centromeres are known to experience less recombination in several species and it has recently been suggested that low recombination rates can facilitate the accumulation and maintenance of isolation genes in partially isolated populations. Therefore, we measured the recombination rate in division 5D/6 directly and estimate that it is at least 16-fold reduced across this region compared to the remainder of the X chromosome. Additionally, sequence data from four loci from field-collected mosquitoes from several West African countries show very strong differentiation between the molecular forms in division 5D/6, whereas none was observed in two loci elsewhere on the X chromosome. Furthermore, genetic variation was substantially lower in division 5D/6 compared to the two reference loci, and the inferred genealogies of the division 5D/6 genes show patterns consistent with selective sweeps. This suggests that the reduced recombination rate has increased the effect of selection on this region and that our data are consistent with the hypothesis that reduced recombination rates can play a role in the accumulation of isolation genes in the face of gene flow.

Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana

The malaria mosquito Anopheles gambiae s.s. is rapidly becoming a model for studies on the evolution of reproductive isolation. Debate has centered on the taxonomic status of two forms (denoted M and S) within the nominal taxon identified by point mutations in the X-linked rDNA region. Evidence is accumulating that there are significant barriers to gene flow between these forms, but that the barriers are not complete throughout the entire range of their distribution. We sampled populations from across Ghana and southern Burkina Faso, West Africa, from areas where the molecular forms occurred in both sympatry and allopatry. Neither Bayesian clustering methods nor F(ST)-based analysis of microsatellite data found differentiation between the M and S molecular forms, but revealed strong differentiation among different ecological zones, irrespective of M/S status and with no detectable effect of geographical distance. Although no M/S hybrids were found in the samples, admixture analysis detected evidence of contemporary interform gene flow, arguably most pronounced in southern Ghana where forms occur sympatrically. Thus, in the sampled area of West Africa, lack of differentiation between M and S forms likely reflects substantial introgression, and ecological barriers appear to be of greater importance in restricting gene flow.

Reproductive output of female Anopheles gambiae (Diptera: Culicidae): comparison of molecular forms

Knowledge of ecological differences between the molecular forms of Anopheles gambiae Giles (Diptera: Culicidae) might lead to understanding of their unique contribution to disease transmission, to better vector control, and to identification of the forces that have separated them. We compared female fecundity measured as egg batch size in relation to body size between the molecular forms in Mali and contrasted them with their sibling species, Anopheles arabiensis Patton. To determine whether eggs of different egg batches are of similar "quality," we compared the total protein content of first-stage larvae (L1s), collected < 2 h after hatching in deionized water. Egg batch size significantly varied between An. gambiae and An. arabiensis and between the molecular forms of An. gambiae (mean batch size was 186.3, 182.5, and 162.0 eggs in An. arabiensis and the M and the S molecular form of An. gambiae, respectively). After accommodating female body size, however, the difference in batch size was not significant. In the S molecular form, egg protein content was not correlated with egg batch size (r = -0.08, P > 0.7) nor with female body size (r = -0.18, P > 0.4), suggesting that females with more resources invest in more eggs rather than in higher quality eggs. The mean total protein in eggs of the M form (0.407 microg per L1) was 6% higher than that of the S form (0.384 microg per L1), indicating that the M form invests a greater portion of her resources into current (rather than future) reproduction. A greater investment per offspring coupled with larger egg batch size may reflect an adaptation of the M form to low productivity larval sites as independent evidence suggests.

P elements and MITE relatives in the whole genome sequence of Anopheles gambiae

Background

Miniature Inverted-repeat Terminal Elements (MITEs), which are particular class-II transposable elements (TEs), play an important role in genome evolution, because they have very high copy numbers and display recurrent bursts of transposition. The 5' and 3' subterminal regions of a given MITE family often show a high sequence similarity with the corresponding regions of an autonomous Class-II TE family. However, the sustained presence over a prolonged evolutionary time of MITEs and TE master copies able to promote their mobility has been rarely reported within the same genome, and this raises fascinating evolutionary questions.

Results

We report here the presence of P transposable elements with related MITE families in the Anopheles gambiae genome. Using a TE annotation pipeline we have identified and analyzed all the P sequences in the sequenced A. gambiae PEST strain genome. More than 0.49% of the genome consists of P elements and derivates. P elements can be divided into 9 different subfamilies, separated by more than 30% of nucleotide divergence. Seven of them present full length copies. Ten MITE families are associated with 6 out of the 9 Psubfamilies. Comparing their intra-element nucleotide diversities and their structures allows us to propose the putative dynamics of their emergence. In particular, one MITE family which has a hybrid structure, with ends each of which is related to a different P-subfamily, suggests a new mechanism for their emergence and their mobility.

Conclusion

This work contributes to a greater understanding of the relationship between full-length class-II TEs and MITEs, in this case P elements and their derivatives in the genome of A. gambiae. Moreover, it provides the most comprehensive catalogue to date of P-like transposons in this genome and provides convincing yet indirect evidence that some of the subfamilies have been recently active.

The indoor use of plastic sheeting pre-impregnated with insecticide for control of malaria vectors

OBJECTIVE

To evaluate the efficacy of permethrin-treated plastic sheeting (ITPS) applied as a lining to the ceiling or walls of rooms against pyrethroid-susceptible and pyrethroid-resistant Anopheles gambiae.

METHOD

Overnight tests were carried out in veranda-trap experimental huts in Vallée du Kou, where two molecular forms of A. gambiae, S and M, occur. The S form is mostly pyrethroid resistant due to the kdr mechanism, and the M form is mostly kdr susceptible. A variety of ITPS covered surfaces were tested, ranging from ceiling only to all walls plus ceiling covered.

RESULTS

ITPS had a major effect on the mortality of mosquitoes, the proportion killed being dependent upon the surface area covered. Homozygotes for kdr resistance showed lower rates of mortality than did heterozygotes or homozygotes for susceptibility. Deterred entry of mosquitoes and inhibition of blood feeding were also correlated with surface area covered. The mode of action and efficacy of ITPS seems to bear closer resemblance to that induced by indoor residual spraying (IRS) than to that induced by insecticide-treated nets.

CONCLUSIONS

ITPS might be conceived as being equivalent to long-lasting or permanent IRS but without some of the operational constraints normally associated with spraying. High coverage of ITPS could potentially have a mass population effect on mosquitoes and give rise to long-term community protection against malaria. A phase III trial is justified to assess the acceptability of ITPS and its efficacy against malaria.

Molecular characterization of the malaria vector Anopheles gambiae s.s. in Madagascar

Anopheles gambiae s.s. Giles (Diptera: Culicidae), the primary African malaria vector, has been characterized at the subspecies level in Madagascar, where only the molecular form S and haplotype gIA occur. The haplotype gIC proposed by other authors was not observed amongst the 35 mosquito genomes sequenced. These S/gIA characteristics are also found on the Comoros archipelago and in continental Africa.

The distribution of hatching time in Anopheles gambiae

Background

Knowledge of the ecological differences between the molecular forms of Anopheles gambiae and their sibling species, An. arabiensis might lead to understanding their unique contribution to disease transmission and to better vector control as well as to understanding the evolutionary forces that have separated them.

Methods

The distributions of hatching time of eggs of wild An. gambiae and An. arabiensis females were compared in different water types. Early and late hatchers of the S molecular form were compared with respect to their total protein content, sex ratio, development success, developmental time and adult body size.

Results

Overall, the distribution of hatching time was strongly skewed to the right, with 89% of the eggs hatching during the second and third day post oviposition, 10% hatching during the next four days and the remaining 1% hatching over the subsequent week. Slight, but significant differences were found between species and between the molecular forms in all water types. Differences in hatching time distribution were also found among water types (in each species and molecular form), suggesting that the eggs change their hatching time in response to chemical factors in the water. Early hatchers were similar to late hatchers except that they developed faster and produced smaller adults than late hatchers.

Conclusion

Differences in hatching time and speed of development among eggs of the same batch may be adaptive if catastrophic events such as larval site desiccation are not rare and the site's quality is unpredictable. The egg is not passive and its hatching time depends on water factors. Differences in hatching time between species and molecular forms were slight, probably reflecting that conditions in their larval sites are rather similar.

Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville, Gabon

Point mutations in the voltage-gated sodium channel gene involved in knockdown resistance to DDT and pyrethroid insecticides have been described in several insect species. In the malaria vector Anopheles gambiae Giles sensu stricto (Diptera: Culicidae) two mutations have been identified. The first, consisting of a leucine-phenylalanine substitution at amino acid position 1014, is widespread in West Africa. The second, a leucine-serine substitution at the same position, has to date only been detected in western Kenya. Analysis of the kdr polymorphism in a sample of 106 An. gambiae s.s. of the rDNA S-form/Type I collected in Libreville (Gabon) surprisingly revealed the presence of both East and West African kdr mutations with frequencies of 63% and 37%, respectively. No wild-type alleles were detected and there was an excess of heterozygous genotypes (P = 0.04). In addition, an inconsistency was found during the kdr genotyping procedures by polymerase chain reaction, which could have lead to an underestimation of resistance alleles. The implications of these findings are discussed.

Malaria vectors in Angola: distribution of species and molecular forms of the Anopheles gambiae complex, their pyrethroid insecticide knockdown resistance (kdr) status and Plasmodium falciparum sporozoite rates

Background

Malaria is by far the greatest cause of morbidity and mortality in Angola, being responsible for 50% of all outpatient attendance and around 22% of all hospital deaths, yet by 2003 only 2% of under-5s used insecticide-treated nets. Entomological studies are an essential foundation for rational malaria control using insecticide-treated nets and indoor residual spraying, but there have been no published studies of malaria vectors in Angola over the 27 years of the civil war, to its end in 2002. This paper describes studies arising from a WHO-sponsored visit in support of the National Malaria Control Programme.

Methods

During April 2001, mosquitoes were sampled by indoor pyrethrum spray collection from four sites in the semi-arid coastal provinces of Luanda and Benguela and two sites in Huambo province, in the humid tropical highlands. Collections took place towards the end of the rainy season and were used to determine the Anopheles species present, their sporozoite rates and the frequency of a kdr allele conferring resistance to pyrethroid insecticides.

Results

A PCR test for the Anopheles gambiae complex showed a preponderance of An. gambiae, with indoor resting densities ranging from 0.9 to 23.5 per house. Of 403 An. gambiae identified to molecular form, 93.5% were M-form and 6.5% S-form. M and S were sympatric at 4 sites but no M/S hybrids were detected. The highest proportion of S-form (20%) was in samples from Huambo, in the humid highlands. Anopheles funestus was found at one site near Luanda. The sporozoite rate of mosquitoes, determined by an ELISA test, was 1.9% for An. gambiae (n = 580) and 0.7% for An. funestus (n = 140). Of 218 An. gambiae (195 M-form and 23 S-form) genotyped for the West African kdr-resistance allele, all were homozygous susceptible.

Conclusion

An. gambiae M-form is the most important and widespread malaria vector in the areas studied but more extensive studies of malaria vectors are required to support the malaria control programme in Angola. These should include standard insecticide resistance biossays and molecular assays that can detect both metabolic resistance and target site insensitivity.

Genetic structure of Anopheles gambiae populations on islands in northwestern Lake Victoria, Uganda

Background

Alternative means of malaria control are urgently needed. Evaluating the effectiveness of measures that involve genetic manipulation of vector populations will be facilitated by identifying small, genetically isolated vector populations. The study was designed to use variation in microsatellite markers to look at genetic structure across four Lake Victoria islands and two surrounding mainland populations and for evidence of any restriction to free gene flow.

Methods

Four Islands (from 20–50 km apart) and two surrounding mainland populations (96 km apart) were studied. Samples of indoor resting adult mosquitoes, collected over two consecutive years, were genotyped at microsatellite loci distributed broadly throughout the genome and analysed for genetic structure, effective migration (Nem) and effective population size (Ne).

Results

Ne estimates showed island populations to consist of smaller demes compared to the mainland ones. Most populations were significantly differentiated geographically, and from one year to the other. Average geographic pair-wise FST ranged from 0.014–0.105 and several pairs of populations had Ne m < 3. The loci showed broad heterogeneity at capturing or estimating population differences.

Conclusion

These island populations are significantly genetically differentiated. Differences reoccurred over the study period, between the two mainland populations and between each other. This appears to be the product of their separation by water, dynamics of small populations and local adaptation. With further characterisation these islands could become possible sites for applying measures evaluating effectiveness of control by genetic manipulation.