Identification and analysis of Single Nucleotide Polymorphisms (SNPs) in the mosquito Anopheles funestus, malaria vector

BACKGROUND: Single nucleotide polymorphisms (SNPs) are the most common source of genetic variation in eukaryotic species and have become an important marker for genetic studies. The mosquito Anopheles funestus is one of the major malaria vectors in Africa and yet, prior to this study, no SNPs have been described for this species. Here we report a genome-wide set of SNP markers for use in genetic studies on this important human disease vector. RESULTS: DNA fragments from 50 genes were amplified and sequenced from 21 specimens of An. funestus. A third of specimens were field collected in Malawi, a third from a colony of Mozambican origin and a third form a colony of Angolan origin. A total of 494 SNPs including 303 within the coding regions of genes and 5 indels were identified. The physical positions of these SNPs in the genome are known. There were on average 7 SNPs per kilobase similar to that observed in An. gambiae and Drosophila melanogaster. Transitions outnumbered transversions, at a ratio of 2:1. The increased frequency of transition substitutions in coding regions is likely due to the structure of the genetic code and selective constraints. Synonymous sites within coding regions showed a higher polymorphism rate than non-coding introns or 3' and 5'flanking DNA with most of the substitutions in coding regions being observed at the 3(rd )codon position. A positive correlation in the level of polymorphism was observed between coding and non-coding regions within a gene. By genotyping a subset of 30 SNPs, we confirmed the validity of the SNPs identified during this study. CONCLUSION: This set of SNP markers represents a useful tool for genetic studies in An. funestus, and will be useful in identifying candidate genes that affect diverse ranges of phenotypes that impact on vector control, such as resistance insecticide, mosquito behavior and vector competence

Genome-Wide Transcription and Functional Analyses Reveal Heterogeneous Molecular Mechanisms Driving Pyrethroids Resistance in the Major Malaria Vector Anopheles funestus Across Africa.

Pyrethroid resistance in malaria vector, An. funestus is increasingly reported across Africa, threatening the sustainability of pyrethroid-based control interventions, including long lasting insecticidal nets (LLINs). Managing this problem requires understanding of the molecular basis of the resistance from different regions of the continent, to establish whether it is being driven by a single or independent selective events. Here, using a genome-wide transcription profiling of pyrethroid resistant populations from southern (Malawi), East (Uganda), and West Africa (Benin), we investigated the molecular basis of resistance, revealing strong differences between the different African regions. The duplicated cytochrome P450 genes (CYP6P9a and CYP6P9b) which were highly overexpressed in southern Africa are not the most upregulated in other regions, where other genes are more overexpressed, including GSTe2 in West (Benin) and CYP9K1 in East (Uganda). The lack of directional selection on both CYP6P9a and CYP6P9b in Uganda in contrast to southern Africa further supports the limited role of these genes outside southern Africa. However, other genes such as the P450 CYP9J11 are commonly overexpressed in all countries across Africa. Here, CYP9J11 is functionally characterized and shown to confer resistance to pyrethroids and moderate cross-resistance to carbamates (bendiocarb). The consistent overexpression of GSTe2 in Benin is coupled with a role of allelic variation at this gene as GAL4-UAS transgenic expression in Drosophila flies showed that the resistant 119F allele is highly efficient in conferring both DDT and permethrin resistance than the L119. The heterogeneity in the molecular basis of resistance and cross-resistance to insecticides in An. funestus populations throughout sub-Saharan African should be taken into account in designing resistance management strategies

Insecticide Resistance in Malaria Vectors: An Update at a Global Scale

Malaria remains the deadliest vector-borne disease in the world. With nearly half of the world’s population at risk, 216 million people suffered from malaria in 2016, with over 400,000 deaths, mainly in sub-Saharan Africa. Important global efforts have been made to eliminate malaria leading to significant reduction in malaria cases and mortality in Africa by 42% and 66%, respectively. Early diagnosis, improved drug therapies and better health infrastructure are key components, but this extraordinary success is mainly due the use of long-lasting insecticidal nets (LLINs) and indoor residual sprayings (IRS) of insecticide. Unfortunately, the emergence and spread of resistance in mosquito populations against insecticides is jeopardising the effectiveness of the most efficient malaria control interventions. To help establish suitable resistance management strategies, it is vital to better understand the distribution of resistance, its mechanisms and impact on effectiveness of control interventions and malaria transmission. In this chapter, we present the current status of insecticide resistance worldwide in main malaria vectors as well as its impact on malaria transmission, and discuss the molecular mechanisms and future perspectives

Investigation of mechanisms of bendiocarb resistance in Anopheles gambiae populations from the city of Yaoundé, Cameroon

BACKGROUND Resistance to the carbamate insecticide bendiocarb is emerging in Anopheles gambiae populations from the city of Yaoundé in Cameroon. However, the molecular basis of this resistance remains uncharacterized. The present study objective is to investigate mechanisms promoting resistance to bendiocarb in An. gambiae populations from Yaoundé. METHODS The level of susceptibility of An. gambiae s.l. to bendiocarb 0.1 % was assessed from 2010 to 2013 using bioassays. Mosquitoes resistant to bendiocarb, unexposed and susceptible mosquitoes were screened for the presence of the Ace-1(R) mutation using TaqMan assays. Microarray analyses were performed to assess the pattern of genes differentially expressed between resistant, unexposed and susceptible. RESULTS Bendiocarb resistance was more prevalent in mosquitoes originating from cultivated sites compared to those from polluted and unpolluted sites. Both An. gambiae and Anopheles coluzzii were found to display resistance to bendiocarb. No G119S mutation was detected suggesting that resistance was mainly metabolic. Microarray analysis revealed the over-expression of several cytochrome P450 s genes including cyp6z3, cyp6z1, cyp12f2, cyp6m3 and cyp6p4. Gene ontology (GO) enrichment analysis supported the detoxification role of cytochrome P450 s with several GO terms associated with P450 activity significantly enriched in resistant samples. Other detoxification genes included UDP-glucosyl transferases, glutathione-S transferases and ABC transporters. CONCLUSION The study highlights the probable implication of metabolic mechanisms in bendiocarb resistance in An. gambiae populations from Yaoundé and stresses the need for further studies leading to functional validation of detoxification genes involved in this resistance

Mapping a Quantitative Trait Locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus

BACKGROUND: Pyrethroid resistance in Anopheles funestus populations has led to an increase in malaria transmission in southern Africa. Resistance has been attributed to elevated activities of cytochrome P450s but the molecular basis underlying this metabolic resistance is unknown. Microsatellite and SNP markers were used to construct a linkage map and to detect a quantitative trait locus (QTL) associated with pyrethroid resistance in the FUMOZ-R strain of An. funestus from Mozambique. RESULTS: By genotyping 349 F(2 )individuals from 11 independent families, a single major QTL, rp1, at the telomeric end of chromosome 2R was identified. The rp1 QTL appears to present a major effect since it accounts for more than 60% of the variance in susceptibility to permethrin. This QTL has a strong additive genetic effect with respect to susceptibility. Candidate genes associated with pyrethroid resistance in other species were physically mapped to An. funestus polytene chromosomes. This showed that rp1 is genetically linked to a cluster of CYP6 cytochrome P450 genes located on division 9 of chromosome 2R and confirmed earlier reports that pyrethroid resistance in this strain is not associated with target site mutations (knockdown resistance). CONCLUSION: We hypothesize that one or more of these CYP6 P450s clustered on chromosome 2R confers pyrethroid resistance in the FUMOZ-R strain of An. funestus

High Level of Pyrethroid Resistance in an Anopheles funestus Population of the Chokwe District in Mozambique

Background Although Anopheles funestus is difficult to rear, it is crucial to analyse field populations of this malaria vector in order to successfully characterise mechanisms of insecticide resistance observed in this species in Africa. In this study we carried out a large-scale field collection and rearing of An. funestus from Mozambique in order to analyse its susceptibility status to insecticides and to broadly characterise the main resistance mechanisms involved in natural populations. Methodology/Principal Findings 3,000 F1 adults were obtained through larval rearing. WHO susceptibility assays indicated a very high resistance to pyrethroids with no mortality recorded after 1h30min exposure and less than 50% mortality at 3h30min. Resistance to the carbamate, bendiocarb was also noted, with 70% mortality after 1h exposure. In contrast, no DDT resistance was observed, indicating that no kdr-type resistance was involved. The sequencing of the acetylcholinesterase gene indicated the absence of the G119S and F455W mutations associated with carbamate and organophosphate resistance. This could explain the absence of malathion resistance in this population. Both biochemical assays and quantitative PCR implicated up-regulated P450 genes in pyrethroid resistance, with GSTs playing a secondary role. The carbamate resistance observed in this population is probably conferred by the observed altered AChE with esterases also involved. Conclusion/Significance The high level of pyrethroid resistance in this population despite the cessation of pyrethroid use for IRS in 1999 is a serious concern for resistance management strategies such as rotational use of insecticides. As DDT has now been re-introduced for IRS, susceptibility to DDT needs to be closely monitored to prevent the appearance and spread of resistance to this insecticide

Gene amplification and microsatellite polymorphism underlie a recent insect host shift

Host plant shifts of herbivorous insects may be a first step toward sympatric speciation and can create new pests of agriculturally important crops; however, the molecular mechanisms that mediate this process are poorly understood. Certain races of the polyphagous aphid Myzus persicae have recently adapted to feed on tobacco (Myzus persicae nicotianae) and show a reduced sensitivity to the plant alkaloid nicotine and cross-resistance to neonicotinoids a class of synthetic insecticides widely used for control. Here we show constitutive overexpression of a cytochrome P450 (CYP6CY3) allows tobacco-adapted races of M. persicae to efficiently detoxify nicotine and has preadapted them to resist neonicotinoid insecticides. CYP6CY3, is highly overexpressed in M. persicae nicotianae clones from three continents compared with M. persicae s.s. and expression level is significantly correlated with tolerance to nicotine. CYP6CY3 is highly efficient (compared with the primary human nicotine-metabolizing P450) at metabolizing nicotine and neonicotinoids to less toxic metabolites in vitro and generation of transgenic Drosophila expressing CYP6CY3 demonstrate that it confers resistance to both compounds in vivo. Overexpression of CYP6CY3 results from the expansion of a dinucleotide microsatellite in the promoter region and a recent gene amplification, with some aphid clones carrying up to 100 copies. We conclude that the mutations leading to overexpression of CYP6CY3 were a prerequisite for the host shift of M. persicae to tobacco and that gene amplification and microsatellite polymorphism are evolutionary drivers in insect host adaptation

Genomic Footprints of Selective Sweeps from Metabolic Resistance to Pyrethroids in African Malaria Vectors Are Driven by Scale up of Insecticide-Based Vector Control

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the `resistance curve' and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistanceassociated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised

Geographical distribution of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) and genetic diversity of invading population of Ae. albopictus in the Republic of the Congo [version 1; referees: 2 approved]

Background: The arbovirus vector, Aedes albopictus, originating from Asia, has recently invaded African countries, including the Republic of the Congo, where it was associated with a chikungunya outbreak. Up until now, little was known about its distribution in relation to the native Aedes aegypti and how the invasion will modify the epidemiology of arboviral diseases. Here, we assessed the current distribution of Ae. albopictus and Ae. aegypti in the Republic of the Congo and explored the genetic diversity of the invading species, Ae. albopictus. Methods: Immature stages of Aedes were collected in nine locations in the Republic of the Congo in 2017 following a north-south transect and reared to adult stage. Adults were morphologically identified, counted and grouped according to species and location. Genetic diversity of Ae. albopictus was assessed by analyzing the cytochrome oxidase I (COI) gene. Results: Ae. albopictus and Ae. aegypti were found together across the country in all the locations investigated. The invasive species is predominant over the native species in all locations except Brazzaville, suggesting that Ae. albopictus is displacing Ae. aegypti across Congo. When comparing the species distributions across the two largest cities, Brazzaville and Pointe Noire, Ae. albopictus was more prevalent than Ae. aegypti in the suburbs whereas the opposite situation was reported in the city centre. Mitochondrial DNA analysis revealed very low genetic diversity of Ae. albopictus with only three haplotypes recorded across the country supporting the recent introduction of this species in the Republic of the Congo. Phylogenetic tree analysis revealed that Ae. albopictus from Congo originated from other tropical Asian countries such as China, likely as a result of increasing trade links. Conclusion: These findings are important for the implementation of vector control strategies and can serve as a foundation for further research on these vectors in the country

Pyrethroid Resistance in an Anopheles funestus Population from Uganda

Background: The susceptibility status of Anopheles funestus to insecticides remains largely unknown in most parts of Africa because of the difficulty in rearing field-caught mosquitoes of this malaria vector. Here we report the susceptibility status of the An. funestus population from Tororo district in Uganda and a preliminary haracterisation of the putative resistance mechanisms involved. Methodology/Principal Findings: A new forced egg laying technique used in this study significantly increased the numbers of field-caught females laying eggs and generated more than 4000 F1 adults. WHO bioassays indicated that An. funestus in Tororo is resistant to pyrethroids (62% mortality after 1 h exposure to 0.75% permethrin and 28% mortality to 0.05% deltamethrin). Suspected DDT resistance was also observed with 82% mortality. However this population is fully susceptible to bendiocarb (carbamate), malathion (organophosphate) and dieldrin with 100% mortality observed after exposure to each of these insecticides. Sequencing of a fragment of the sodium channel gene containing the 1014 codon conferring pyrethroid/DDT resistance in An. gambiae did not detect the L1014F kdr mutation but a correlation between haplotypes and resistance phenotype was observed indicating that mutations in other exons may be conferring the knockdown resistance in this species. Biochemical assays suggest that resistance in this population is mediated by metabolic resistance with elevated level of GSTs, P450s and pNPA compared to a susceptible strain of Anopheles gambiae. RT-PCR further confirmed the involvement of P450s with a 12-fold over-expression of CYP6P9b in the Tororo population compared to the fully susceptible laboratory colony FANG. Conclusion: This study represents the first report of pyrethroid/DDT resistance in An. funestus from East Africa. With resistance already reported in southern and West Africa, this indicates that resistance in An. funestus may be more widespread than previously assumed and therefore this should be taken into account for the implementation and management of vector control programs in Africa