Introgression between Anopheles gambiae and Anopheles coluzzii in Burkina Faso and its associations with kdr resistance and Plasmodium infection
Malaria Journal volume 18, Article number: 127 (2019)
Insecticide resistance in Anopheles coluzzii mosquitoes has become widespread throughout West Africa including in Burkina Faso. The insecticide resistance allele (kdr or L1014F) is a prime indicator that is highly correlated with phenotypic resistance in West Africa. Studies from Benin, Ghana and Mali have suggested that the source of the L1014F is introgression of the 2L divergence island via interspecific hybridization with Anopheles gambiae. The goal of this study was to characterize local mosquito populations in the Nouna Department, Burkina Faso with respect to: (i) the extent of introgression between An. coluzzii and An. gambiae, (ii) the frequency of the L1014F mutation and (iii) Plasmodium infection rates.
A total of 95 mosquitoes were collected from ten sites surrounding Nouna town in Kossi Province, Burkina Faso in 2012. The species composition, the extent of introgression in An. coluzzii mosquitoes and their Plasmodium infection rates were identified with a modified version of the “Divergence Island SNP” (DIS) genotyping assay.
The mosquito collection contained 70.5% An. coluzzii, 89.3% of which carried a 3 Mb genomic region on the 2L chromosome with L1014F insecticide resistance mutation that was introgressed from An. gambiae. In addition, 22.4% in the introgressed An. coluzzii specimens were infected with Plasmodium falciparum, whereas none of the non-introgressed (“pure”) An. coluzzii were infected.
This paper is the first report providing divergence island SNP genotypes for natural population of Burkina Faso and corresponding Plasmodium infection rates. These observations warrant further study and could have a major impact on future malaria control strategies in Burkina Faso.
The malaria parasite, Plasmodium falciparum, continues to be a significant cause of illness and death in Burkina Faso (BF) . In 2016, it was estimated that there were 7,890,000 malaria cases due to P. falciparum in Burkina Faso, 21,300 of which resulted in death . Two of the main vector species implicated in the spread of malaria across BF are Anopheles gambiae and Anopheles coluzzii . One of the malaria intervention strategies in the country has been the distribution of insecticide-treated bed nets (ITNs): ITN coverage increased from less than 40% in 2010 to over 60% in 2014 . However, malaria cases in BF increased sharply in the period from 2010 to 2016 . This apparent paradox may be explained by the high frequency of 2L chromosome introgression, which includes the knockdown resistance mutation variants (kdr), mutations (designated L1014F and L1014S) which have been associated with increased susceptibility to Plasmodium infection in individuals carrying this mutation .
Samples collected from the Nouna Department in 2010 showed that 64% and 50% of An. gambiae and An. coluzzii, respectively, possessed the L1014F insecticide resistance mutation . Furthermore, the levels of L1014F mutation have remained relatively stable and high in the An. gambiae complex across BF, whereas the distribution and frequency of the L1014S mutation has significantly increased from < 10 to ~ 40% between 2008 and 2014 . Previous studies show that An. coluzzii collected in Benin , Mali , and Ghana  acquired the L1014F mutation by introgression of a divergence island on the 2L chromosome. Nouna Department in BF is geographically located between Benin, Ghana, and Mali; therefore, similar patterns of introgression on the 2L chromosome in An. coluzzii are expected.
Previous studies determined introgression rates in BF based on 2 markers located close to the L1014F mutation . However, whole genome sequencing of introgressed An. coluzzii reveal that the genomic island of divergence, a genomic region of 3 Mb near the centromere on the 2L chromosome, which is highly diverged between An. coluzzii and An. gambiae, was introgressed from An. gambiae to An. coluzzii [6, 7]. Recently, it has also been shown that introgressed An. coluzzii mosquitoes have higher Plasmodium infection rates [3, 9, 10]. Mitri and his co-workers suggested that this phenotype is not caused by the L1014S or L1014F mutation, but most likely by the serine protease ClipC9 gene that is located on the 3 Mb genomic region that is introgressed . The ‘Divergence Islands SNP’ assay  was improved to simultaneously identify species, introgression of the chromosome 2L 3 Mb genomic region, insecticide resistance, and Plasmodium infection. The goal of this study is to determine the introgression status and corresponding Plasmodium infection rates in the Nouna Department, BF.
Sample collection and DNA extraction
Anopheles gambiae and An. coluzzii were collected from ten villages within a 30 km radius in the Nouna Department, BF in 2012 (Table 1). DNA was extracted from head/thorax tissue using the QIAGEN Biosprint 96 system with QIAGEN blood tissue reagents following established protocols [12, 13]. A modified version of the “Divergence Island SNP” (DIS) assay [11, 14] was used to distinguish An. gambiae, An. coluzzii and introgressed An. coluzzii. Anopheles gambiae and An. coluzzii are differentiated in the DIS assay using SNP markers on the X, and 2L and 3L chromosomes . All SNPs are located within three unlinked pericentromeric genome regions known as the islands of divergence [15, 16]. Anopheles coluzzii samples were considered “pure” when the majority (> 15/18) of DIS markers had An. coluzzii specific genotypes. An An. coluzzii sample with more than two gambiae-specific markers is considered “introgressed”. This is consistent with previous studies [11, 14]. The L1014F mutation was not considered a species-specific marker. Plasmodium infection rates were assessed by including mitochondrial DNA makers of Plasmodium species that can distinguish Plasmodium ovale, Plasmodium malariae, Plasmodium vivax, and P. falciparum (see Additional file 1). Some of these markers are from the Anopheles multi-detection assay . To accommodate all these in a single multiplex assay, some of previous DIS markers (1039-358, 0407-337 SNPs) were replace with 00819-1180 SNP. The final set of markers used for the new DIS assay and the related metadata are provided (see Additional file 1).
Results and discussion
Insecticide resistant Anopheles coluzzii show patterns of introgression
A total of 95 Anopheles mosquitoes were collected from villages within 30 km radius in Nouna Department, BF. The species composition based on genotyping is listed in Table 1 and illustrated on the map in Fig. 1a. The raw genotyping results are listed (see Additional file 2). Anopheles coluzzii is the most abundant vector species in the majority of villages. In the village of Bourasso equal numbers of An. arabiensis and An. coluzzii were found. In Kamadena, the location furthest away from Nouna town, one An. gambiae sensu stricto (s.s.) specimen was found, together with An. arabiensis and An. coluzzii. The relative abundance of An. arabiensis (20%) observed near Nouna differs from previous reports that An. arabiensis was absent in that location . This observation suggests that An. arabiensis can increase in abundance after successful introduction of insecticide-treated bed nets, as observed in many locations throughout sub-Saharan Africa [18,19,20], possibly because it is less affected by insecticide-treated bed nets compared with An. gambiae . Furthermore, the relative abundance of An. gambiae observed (1.1%) is significantly lower than that reported by Namountougou and co-workers, who found An. gambiae at a frequency of 23.3% in their survey . The overall pattern is consistent with a collapse of the local An. gambiae s.s. population, apparently being replaced by introgressed An. coluzzii carrying the L1014F (kdr) mutation, as observed in Mali .
The majority of specimens we collected were An. coluzzii (75 out of 95). Of those, 89.3% were individuals that showed introgression of the chromosome 2L genomic island of divergence . The 2L introgressed An. coluzzii includes the L1014F mutation originating from An. gambiae (Table 2). One heterozygous L1014F mutation in a “pure” An. coluzzii individual was observed, which based on the genotyping assay, showed no further history of introgression with An. gambiae (all An. coluzzii specific SNPs present on the 2L chromosome). Similar patterns of introgression of the 2L chromosome were observed in Mali [6, 21], Ghana  and Guinea Bissau . In addition, 3 out of 19 An. arabiensis mosquitoes (16.8%) were heterozygous for the L1014S mutation (kdr-e) (Table 2). This variant, originally found in East Africa, is now increasingly found in West Africa [22, 23].
Introgressed Anopheles coluzzii are more likely to be infected with Plasmodium falciparum
From 95 Anopheles specimens sampled, 18 were positive for P. falciparum (19%, Table 3). This infection rate is comparable with those observed in Cameroon , Kenya  and the Democratic Republic of Congo . Infections by P. ovale, P. vivax, or P. malariae were not observed among the samples.
The infection rate in introgressed An. coluzzii (N = 67) was 22.4%, whereas none of the non-introgressed, “pure” An. coluzzii (N = 8) were infected with Plasmodium (Table 3). The SNP genotyping assay can detect P. falciparum, P. vivax, P. malariae and P. ovale, but all infections detected were with P. falciparum. The difference in Plasmodium infection rates between “pure” and introgressed An. coluzzii were not statistically significance using Fisher’s Exact Test (p = 0.345) due to the low number of pure An. coluzzii mosquitoes. Power calculations based on this data indicate that further studies with increased sample size (N > 190) would likely be required to provide sufficient statistical rigor.
Although differences in the Plasmodium infection rate between introgressed and non-introgressed, “pure” An. coluzzii is not statistically significant, there is evidence in the literature that introgressed An. coluzzii may be more susceptible to Plasmodium infection . Silencing the serine protease ClipC9 gene, which is located within the 2L introgressed genomic region, does significantly increase parasite numbers in An. coluzzii laboratory colony mosquitoes . One may also suspect that increased longevity due to increased insecticide tolerance could affect the rate of Plasmodium infection in introgressed An. coluzzii. Moreover, whether this pattern would hold over multiple collections across different seasons remain to be determined. Further research is required to understand the true impact of this gene in natural An. coluzzii populations.
The results indicate that introgressed An. coluzzii is the primary malaria vector in the Nouna Department of BF. These populations carry the L1014F and L1014S insecticide resistance mutations and also have higher Plasmodium infection rates compared to the other Anopheles populations observed in the region. These observations warrant further study and could have a major impact on future malaria control strategies in BF.
WHO. World malaria report 2017. Geneva: World Health Organization; 2017.
Namountougou M, Simard F, Baldet T, Diabaté A, Ouédraogo JB, Martin T, et al. Multiple insecticide resistance in Anopheles gambiae s.l. populations from Burkina Faso, West Africa. PLoS ONE. 2012;7:e48412.
Mitri C, Markianos K, Guelbeogo WM, Bischoff E, Gneme A, Eiglmeier K, et al. The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing. Malar J. 2015;14:391.
Dabiré RK, Namountougou M, Diabaté A, Soma DD, Bado J, Toé HK, et al. Distribution and frequency of kdr mutations within Anopheles gambiae s.l. populations and first report of the Ace. 1G119S mutation in Anopheles arabiensis from Burkina Faso (West Africa). PLoS ONE. 2014;9:e101484.
Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, et al. The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression. Insect Mol Biol. 2000;9:451–5.
Norris LC, Main BJ, Lee Y, Collier TC, Fofana A, Cornel AJ, et al. Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets. Proc Natl Acad Sci USA. 2015;12:815–20.
Clarkson CS, Weetman D, Essandoh J, Yawson AE, Maslen G, Manske M, et al. Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation. Nat Commun. 2014;5:4248.
Diabate A, Brengues C, Baldet T, Dabiré K, Hougard J-M, Akogbeto M, et al. The spread of the Leu-Phe kdr mutation through Anopheles gambiae complex in Burkina Faso: genetic introgression and de novo phenomena. Trop Med Int Health. 2004;9:1267–73.
Ndiath MO, Cailleau A, Diedhiou SM, Gaye A, Boudin C, Richard V, et al. Effects of the kdr resistance mutation on the susceptibility of wild Anopheles gambiae populations to Plasmodium falciparum: a hindrance for vector control. Malar J. 2014;13:340.
Alout H, Dabiré RK, Djogbénou LS, Abate L, Corbel V, Chandre F, et al. Interactive cost of Plasmodium infection and insecticide resistance in the malaria vector Anopheles gambiae. Sci Rep. 2016;6:29755.
Lee Y, Marsden CD, Nieman C, Lanzaro GC. A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae. Mol Ecol Resour. 2014;14:297–305.
Nieman CC, Yamasaki Y, Collier TC, Lee Y. A DNA extraction protocol for improved DNA yield from individual mosquitoes. F1000Res. 2015;4:1314.
Yamasaki YK, Nieman CC, Chang AN, Collier TC, Main BJ, Lee Y. Improved tools for genomic DNA library construction of small insects. F1000Res. 2016. https://doi.org/10.1080/23802359.2018.1495117.
Lee Y, Olson N, Yamasaki Y, Chang A, Marsden C, Ouledi A, et al. Absence of kdr resistance alleles in the Union of the Comoros, East Africa. F1000Res. 2015;4:146.
Turner TL, Hahn MW, Nuzhdin SV. Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005;3:e285.
White BJ, Cheng C, Simard F, Costantini C, Besansky NJ. Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae. Mol Ecol. 2010;19:925–39.
Lee Y, Weakley AM, Nieman CC, Malvick J, Lanzaro GC. A multi-detection assay for malaria transmitting mosquitoes. J Vis Exp. 2015;96:e52385.
Mwangangi JM, Mbogo CM, Orindi BO, Muturi EJ, Midega JT, Nzovu J, et al. Shifts in malaria vector species composition and transmission dynamics along the Kenyan coast over the past 20 years. Malar J. 2013;12:13.
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, et al. Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010;9:62.
Russell TL, Lwetoijera DW, Maliti D, Chipwaza B, Kihonda J, Charlwood JD, et al. Impact of promoting longer-lasting insecticide treatment of bed nets upon malaria transmission in a rural Tanzanian setting with pre-existing high coverage of untreated nets. Malar J. 2010;9:187.
Main BJ, Lee Y, Collier TC, Norris LC, Brisco K, Fofana A, et al. Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali. Mol Ecol. 2015;24:5145–57.
Ndiath MO, Cailleau A, Orlandi-Pradines E, Bessell P, Pagès F, Trape J-F, et al. Emerging knock-down resistance in Anopheles arabiensis populations of Dakar, Senegal: first evidence of a high prevalence of kdr-e mutation in West African urban area. Malar J. 2015;14:364.
Djègbè I, Boussari O, Sidick A, Martin T, Ranson H, Chandre F, et al. Dynamics of insecticide resistance in malaria vectors in Benin: first evidence of the presence of L1014S kdr mutation in Anopheles gambiae from West Africa. Malar J. 2011;10:261.
Ndo C, Kopya E, Donbou MA, Njiokou F, Awono-Ambene P, Wondji C. Elevated Plasmodium infection rates and high pyrethroid resistance in major malaria vectors in a forested area of Cameroon highlight challenges of malaria control. Parasit Vectors. 2018;11:157.
Ogola E, Villinger J, Mabuka D, Omondi D, Orindi B, Mutunga J, et al. Composition of Anopheles mosquitoes, their blood-meal hosts, and Plasmodium falciparum infection rates in three islands with disparate bed net coverage in Lake Victoria, Kenya. Malar J. 2017;16:360.
Riveron JM, Watsenga F, Irving H, Irish SR, Wondji CS. High Plasmodium infection rate and reduced bed net efficacy in multiple insecticide-resistant malaria vectors in Kinshasa, Democratic Republic of Congo. J Infect Dis. 2017;217:320–8.
MH, NB, GL, YL conceived and designed the study; NB collected samples; MH, YY, YL collected and analysed the data; OK mapped the data; MH, HH, IE wrote the paper. All authors read and approved the final manuscript.
We would like to thank Allison Chang for her help with the literature search and sample processing.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files.
Consent for publication
Ethics approval and consent to participate
We thank the University of California—Irvine, Malaria Initiative (UCIMI) for their support. This study was supported by University of California, Davis (Bridge Funding).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hanemaaijer, M.J., Higgins, H., Eralp, I. et al. Introgression between Anopheles gambiae and Anopheles coluzzii in Burkina Faso and its associations with kdr resistance and Plasmodium infection. Malar J 18, 127 (2019). https://doi.org/10.1186/s12936-019-2759-1
- Insecticide resistance
- Gene flow
- Malaria vector
- Burkina Faso