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

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

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

<p>Abstract</p> <p>Background</p> <p><it>Anopheles gambiae sensu stricto</it>, 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.</p> <p>Methods</p> <p>Chromosomal, microsatellite and geographic/ecological evidence are synthesized to identify sources of genetic polymorphism among chromosomal and molecular forms of the malaria vector <it>Anopheles gambiae s.s</it>.</p> <p>Results</p> <p>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 <it>An. gambiae s.s</it>. 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.</p> <p>Conclusion</p> <p>Chromosome inversions, microsatellite allele frequencies and habitat preference all indicate that the Forest M form of <it>An. gambiae </it>is genetically distinct from the other recognized forms within the taxon <it>Anopheles gambiae sensu stricto</it>. 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.</p

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

BACKGROUND: Anopheles gambiae sensu stricto (s.s.) is a primary vector of Plasmodium falciparum in sub-Saharan Africa. Although some physiological differences among molecular and chromosomal forms of this species have been demonstrated, the relative susceptibility to malaria parasite infection among them has not been unequivocally shown. The objective of this study was to investigate P. falciparum circumsporozoite protein infection (CSP) positivity among An. gambiae s.s. chromosomal and molecular forms. METHODS: Wild An. gambiae from two sites Kela (n = 464) and Sidarebougou (n = 266) in Mali were screened for the presence of P. falciparum CSP using an enzyme-linked immunosorbent assay (ELISA). Samples were then identified to molecular form using multiple PCR diagnostics (n = 713) and chromosomal form using chromosomal karyotyping (n = 419). RESULTS: Of 730 An. gambiae sensu lato (s.l.) mosquitoes, 89 (12.2%) were CSP ELISA positive. The percentage of positive mosquitoes varied by site: 52 (11.2%) in Kela and 37 (13.9%) in Sidarebougou. Eighty-seven of the positive mosquitoes were identified to molecular form and they consisted of nine Anopheles arabiensis (21.4%), 46 S (10.9%), 31 M (12.8%), and one MS hybrid (14.3%). Sixty of the positive mosquitoes were identified to chromosomal form and they consisted of five An. arabiensis (20.0%), 21 Savanna (15.1%), 21 Mopti (30.4%), 11 Bamako (9.2%), and two hybrids (20.0%). DISCUSSION: In this collection, the prevalence of P. falciparum infection in the M form was equivalent to infection in the S form (no molecular form differential infection). There was a significant differential infection by chromosomal form such that, P. falciparum infection was more prevalent in the Mopti chromosomal forms than in the Bamako or Savanna forms; the Mopti form was also the most underrepresented in the collection. Continued research on the differential P. falciparum infection of An. gambiae s.s. chromosomal and molecular forms may suggest that Plasmodium – An. gambiae interactions play a role in malaria transmission

PCR-based karyotyping of Anopheles gambiae inversion 2Rj identifies the BAMAKO chromosomal form.

BACKGROUND: The malaria vector Anopheles gambiae is polymorphic for chromosomal inversions on the right arm of chromosome 2 that segregate nonrandomly between assortatively mating populations in West Africa. One such inversion, 2Rj, is associated with the BAMAKO chromosomal form endemic to southern Mali and northern Guinea Conakry near the Niger River. Although it exploits a unique ecology and both molecular and chromosomal data suggest reduced gene flow between BAMAKO and other A. gambiae populations, no molecular markers exist to identify this form. METHODS: To facilitate study of the BAMAKO form, a PCR assay for molecular karyotyping of 2Rj was developed based on sequences at the breakpoint junctions. The assay was extensively validated using more than 700 field specimens whose karyotypes were determined in parallel by cytogenetic and molecular methods. As inversion 2Rj also occurs in SAVANNA populations outside the geographic range of BAMAKO, samples were tested from Senegal, Cameroon and western Guinea Conakry as well as from Mali. RESULTS: In southern Mali, where 2Rj polymorphism in SAVANNA populations was very low and most of the 2Rj homozygotes were found in BAMAKO karyotypes, the molecular and cytogenetic methods were almost perfectly congruent. Elsewhere agreement between the methods was much poorer, as the molecular assay frequently misclassified 2Rj heterozygotes as 2R+j standard homozygotes. CONCLUSION: Molecular karyotyping of 2Rj is robust and accurate on 2R+j standard and 2Rj inverted homozygotes. Therefore, the proposed approach overcomes the lack of a rapid tool for identifying the BAMAKO form across developmental stages and sexes, and opens new perspectives for the study of BAMAKO ecology and behaviour. On the other hand, the method should not be applied for molecular karyotyping of j-carriers within the SAVANNA chromosomal form

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

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

Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets

Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa

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

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

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

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

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division [version 2; peer review: 2 approved]

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