Local adaptation and colonization are potential factors affecting sexual competitiveness and mating choice in Anopheles coluzzii populations

The mating behaviour of the malaria vector Anopheles gambiae complex is an important aspect of its reproduction biology. The success of mosquito release programmes based on genetic control of malaria crucially depends on competitive mating between both laboratory-reared and wild individuals, and populations from different localities. It is known that intrinsic and extrinsic factors can influence the mating success. This study addressed some of the knowledge gaps about factors influcencing mosquito mating success. In semi-field conditions, the study compared the mating success of three laboratory-reared and wild allopatric An. coluzzii populations originating from ecologically different locations in Burkina Faso. Overall, it was found that colonization reduced the mating competitiveness of both males and females compared to that of wild type individuals. More importly, females were more likely to mate with males of their own population of origin, be it wild or colonised, suggesting that local adaptation affected mate choice. The observations of mating behaviour of colonized and local wild populations revealed that subtle differences in behaviour lead to significant levels of population-specific mating. This is the first study to highlight the importance of local adaptation in the mating success, thereby highlighting the importance of using local strains for mass-rearing and release of An. coluzzii in vector control programmes.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences.Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Spontaneous mutation rate estimates for the principal malaria vectors Anopheles coluzzii and Anopheles stephensi

Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10−9 (95% confidence interval, 2.06 × 10−10—2.91 × 10−9) and 1.36 × 10−9 (95% confidence interval, 4.42 × 10−10—3.18 × 10−9) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anophelescoluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii

Insecticide resistance is a major threat to gains in malaria control, which have been stalling and potentially reversing since 2015. Studies into the causal mechanisms of insecticide resistance are painting an increasingly complicated picture, underlining the need to design and implement targeted studies on this phenotype. In this study, we compare three populations of the major malaria vector An. coluzzii: a susceptible and two resistant colonies with the same genetic background. The original colonised resistant population rapidly lost resistance over a 6-month period, a subset of this population was reselected with pyrethroids, and a third population of this colony that did not lose resistance was also available. The original resistant, susceptible and re-selected colonies were subject to RNAseq and whole genome sequencing, which identified a number of changes across the transcriptome and genome linked with resistance. Firstly, an increase in the expression of genes within the oxidative phosphorylation pathway were seen in both resistant populations compared to the susceptible control; this translated phenotypically through an increased respiratory rate, indicating that elevated metabolism is linked directly with resistance. Genome sequencing highlighted several blocks clearly associated with resistance, including the 2Rb inversion. Finally, changes in the microbiome profile were seen, indicating that the microbial composition may play a role in the resistance phenotype. Taken together, this study reveals a highly complicated phenotype in which multiple transcriptomic, genomic and microbiome changes combine to result in insecticide resistance.

Evolution of the Pyrethroids Target-Site Resistance Mechanisms in Senegal: Early Stage of the Vgsc-1014F and Vgsc-1014S Allelic Frequencies Shift

The evolution and spread of insecticide resistance mechanisms amongst malaria vectors across the sub-Saharan Africa threaten the effectiveness and sustainability of current insecticide-based vector control interventions. However, a successful insecticide resistance management plan relies strongly on evidence of historical and contemporary mechanisms circulating. This study aims to retrospectively determine the evolution and spread of pyrethroid resistance mechanisms among natural Anopheles gambiae s.l. populations in Senegal. Samples were randomly drawn from an existing mosquito sample, collected in 2013, 2017, and 2018 from 10 sentinel sites monitored by the Senegalese National Malaria Control Programme (NMCP). Molecular species of An. gambiae s.l. and the resistance mutations at the Voltage-gated Sodium Channel 1014 (Vgsc-1014) locus were characterised using PCR-based assays. The genetic diversity of the Vgsc gene was further analyzed by sequencing. The overall species composition revealed the predominance of Anopheles arabiensis (73.08%) followed by An. gambiae s.s. (14.48%), Anopheles coluzzii (10.94%) and Anopheles gambiae–coluzii hybrids (1.48%). Both Vgsc-1014F and Vgsc-1014S mutations were found in all studied populations with a spatial variation of allele frequencies from 3% to 90%; and 7% to 41%, respectively. The two mutations have been detected since 2013 across all the selected health districts, with Vgsc-L1014S frequency increasing over the years while Vgsc-1014F decreasing. At species level, the Vgsc-1014F and Vgsc-1014S alleles were more frequent amongst An. gambiae s.s. (70%) and An. arabiensis (20%). The Vgsc gene was found to be highly diversified with eight different haplotypes shared between Vgsc-1014F and Vgsc-1014S. The observed co-occurrence of Vgsc-1014F and Vgsc-1014S mutations suggest that pyrethroid resistance is becoming a widespread phenomenon amongst malaria vector populations, and the NMCP needs to address this issue to sustain the gain made in controlling malaria.

Transcriptomic analysis of resistance and short-term induction response to pyrethroids, in Anopheles coluzzii legs

Background Insecticide-treated bed nets and indoor residual spraying comprise the major control measures against Anopheles gambiae sl, the dominant vector in sub-Saharan Africa. The primary site of contact with insecticide is through the mosquitoes’ legs, which represents the first barrier insecticides have to bypass to reach their neuronal targets. Proteomic changes and leg cuticle modifications have been associated with insecticide resistance that may reduce the rate of penetration of insecticides. Here, we performed a multiple transcriptomic analyses focusing on An. coluzzii legs. Results Firstly, leg-specific enrichment analysis identified 359 genes including the pyrethroid-binder SAP2 and 2 other chemosensory proteins, along with 4 ABCG transporters previously shown to be leg enriched. Enrichment of gene families included those involved in detecting chemical stimuli, including gustatory and ionotropic receptors and genes implicated in hydrocarbon-synthesis. Subsequently, we compared transcript expression in the legs of a highly resistant strain (VK7-HR) to both a strain with very similar genetic background which has reverted to susceptibility after several generations without insecticide pressure (VK7-LR) and a lab susceptible population (NG). Two hundred thirty-two differentially expressed genes (73 up-regulated and 159 down-regulated) were identified in the resistant strain when compared to the two susceptible counterparts, indicating an over-expression of phase I detoxification enzymes and cuticular proteins, with decrease in hormone-related metabolic processes in legs from the insecticide resistant population. Finally, we analysed the short-term effect of pyrethroid exposure on An. coluzzii legs, comparing legs of 1 h-deltamethrin-exposed An. coluzzii (VK7-IN) to those of unexposed mosquitoes (VK7-HR) and identified 348 up-regulated genes including those encoding for GPCRs, ABC transporters, odorant-binding proteins and members of the divergent salivary gland protein family. Conclusions The data on An. coluzzii leg-specific transcriptome provides valuable insights into the first line of defense in pyrethroid resistant and short-term deltamethrin-exposed mosquitoes. Our results suggest that xenobiotic detoxification is likely occurring in legs, while the enrichment of sensory proteins, ABCG transporters and cuticular genes is also evident. Constitutive resistance is primarily associated with elevated levels of detoxification and cuticular genes, while short-term insecticide-induced tolerance is linked with overexpression of transporters, GPCRs and GPCR-related genes, sensory/binding and salivary gland proteins.

Evaluation of the interaction between insecticide resistance-associated genes and malaria transmission in Anopheles gambiae sensu lato in central Côte d’Ivoire

Background There is evidence that the knockdown resistance gene (Kdr) L1014F and acetylcholinesterase-1 gene (Ace-1R) G119S mutations involved in pyrethroid and carbamate resistance in Anopheles gambiae influence malaria transmission in sub-Saharan Africa. This is likely due to changes in the behaviour, life history and vector competence and capacity of An. gambiae. In the present study, performed as part of a two-arm cluster randomized controlled trial evaluating the impact of household screening plus a novel insecticide delivery system (In2Care Eave Tubes), we investigated the distribution of insecticide target site mutations and their association with infection status in wild An. gambiae sensu lato (s.l.) populations. Methods Mosquitoes were captured in 40 villages around Bouaké by human landing catch from May 2017 to April 2019. Randomly selected samples of An. gambiae s.l. that were infected or not infected with Plasmodium sp. were identified to species and then genotyped for Kdr L1014F and Ace-1R G119S mutations using quantitative polymerase chain reaction assays. The frequencies of the two alleles were compared between Anopheles coluzzii and Anopheles gambiae and then between infected and uninfected groups for each species. Results The presence of An. gambiae (49%) and An. coluzzii (51%) was confirmed in Bouaké. Individuals of both species infected with Plasmodium parasites were found. Over the study period, the average frequency of the Kdr L1014F and Ace-1R G119S mutations did not vary significantly between study arms. However, the frequencies of the Kdr L1014F and Ace-1R G119S resistance alleles were significantly higher in An. gambiae than in An. coluzzii [odds ratio (95% confidence interval): 59.64 (30.81–131.63) for Kdr, and 2.79 (2.17–3.60) for Ace-1R]. For both species, there were no significant differences in Kdr L1014F or Ace-1R G119S genotypic and allelic frequency distributions between infected and uninfected specimens (P > 0.05). Conclusions Either alone or in combination, Kdr L1014F and Ace-1R G119S showed no significant association with Plasmodium infection in wild An. gambiae and An. coluzzii, demonstrating the similar competence of these species for Plasmodium transmission in Bouaké. Additional factors including behavioural and environmental ones that influence vector competence in natural populations, and those other than allele measurements (metabolic resistance factors) that contribute to resistance, should be considered when establishing the existence of a link between insecticide resistance and vector competence. Graphical Abstract

Evaluating the contribution of aestivation to the persistence of malaria mosquitoes through the Sahelian dry season using stable isotopes

Data suggests the malaria vector Anopheles coluzzii persists in the Sahel by dry-season aestivation though evidence is scant. We have marked Anopheles mosquitoes using deuterium (2H) to assess the contribution of aestivation to persistence of mosquitoes through the seven-month dry season. If local aestivation is the only way A. coluzzii persists, the frequency of marked mosquitoes should remain stable throughout, whereas finding no marked mosquitoes would be evidence against aestivation. Larval sites were spiked with 2H at the end of the 2017 wet season in two Sahelian villages in Mali. We monitored 2H-enriched populations until the onset of rains. By the end of the enrichment period, 33% of A. coluzzii mosquitoes were clearly marked. Expectedly, 2H levels in marked mosquitoes degraded over time, resulting in a partial overlap of the marked and non-marked 2H distributions. We utilized three methods to estimate the fraction of marked mosquitoes in the population. Seven months after enrichment, 7% of the population had 2H values above the highest pre-enrichment value. An excess of 21% exceeded the 3rd quartile of the pre-enrichment population. A finite mixed population model showed 2.5% represented a subpopulation of marked mosquitoes with elevated 2H, compatible with our predictions. We provide evidence that aestivation is a major persistence mechanism of A. coluzzii in the Sahel, contributing at least 20% of the adults at the onset of rains, suggesting A. coluzzii utilizes multiple persistence strategies enabling its populations rapid buildup, facilitating subsequent malaria resurgence. These may complicate vector control and malaria elimination campaigns.