Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 14572 (2023)
2
Metrics details
The voltage-gated sodium channel, para, is a target of DDT and pyrethroid class insecticides. Single nucleotide mutations in para, called knockdown resistant or kdr, which contribute to resistance against DDT and pyrethroid insecticides, have been correlated with increased susceptibility of Anopheles to the human malaria parasite Plasmodium falciparum. However, a direct role of para activity on Plasmodium infection has not yet been established. Here, using RNA-mediated silencing, we provide in vivo direct evidence for the requirement of wild-type (wt) para function for insecticide activity of deltamethrin. Depletion of wt para, which is susceptible to insecticide, causes deltamethrin tolerance, indicating that insecticide-resistant kdr alleles are likely phenocopies of loss of para function. We then show that normal para activity in An. coluzzii limits Plasmodium infection prevalence for both P. falciparum and P. berghei. A transcriptomic analysis revealed that para activity does not modulate the expression of immune genes. However, loss of para function led to enteric dysbiosis with a significant increase in the total bacterial abundance, and we show that para function limiting Plasmodium infection is microbiota dependent. In the context of the bidirectional “enteric microbiota-brain” axis studied in mammals, these results pave the way for studying whether the activity of the nervous system could control Anopheles vector competence.
Voltage-gated sodium channels (VGSC) are transmembrane proteins mediating electrical potential (voltage) in excitable cells such as neurons and muscle cells and participate in activities such as muscle contraction1,2,3, cell proliferation4 and cognition5,6. Pyrethroids represent one of the most useful insecticide classes and have been used intensively in agricultural and public health settings7,8,9. Deltamethrin is a pyrethroid that binds to the VGSC, preventing its transition from an activated to an inactivated state, leading to persistent depolarized cell membranes that paralyze and kill the insects10,11.
Malaria is a disease caused by Plasmodium parasites, which are transmitted to humans by the bite of infectious Anopheles mosquito females and has killed a predicted 619,000 persons in 202112. The use of indoor residual spraying and long-lasting insecticide-treated nets have been the foundation of mosquito vector control interventions and have provided major gains in reducing malaria transmission. However, this gain has been hampered by the development and spread of insecticide resistance in mosquito vector populations12. Anopheles resistance to the pyrethroid insecticides is conferred by metabolic detoxification through the action of the P450 monooxygenases or glutathione-S-transferases (GSTs)13,14,15,16,17, by point mutations of the VGSC that inhibit pyrethroid binding at the target site18,19 and by cuticular resistance20. This VGSC is also called para due to the paralytic phenotype observed in Drosophila VGSC mutants21,22. Several mutations identified in the VGSC lead to pyrethroid resistance known as knock-down resistance (kdr) alleles: the L1014F and L1014S – a leucine to phenylalanine substitution19 and a leucine to serine substitution18. More recently, a novel mutation I1527T, adjacent to a predicted pyrethroid-binding site, has been found An. coluzzii from Ghana, Burkina Faso, Cote d’Ivoire and Guinea23. An additional kdr mutation, formerly known as N1575Y, the so called “super kdr”, occurs with and intensifies the effect of the L995F-mediated pyrethroid resistance in An. gambiae from southern Cote d’Ivoire10,24,25.
The relationship between the insecticide resistance phenotype and Anopheles vector competence to Plasmodium falciparum infection is not well understood. In laboratory populations of An. gambiae, the presence of kdr has been correlated with increased infection susceptibility of Plasmodium26,27. In wild An. gambiae sensu lato mosquitoes, the kdr resistance allele was associated with greater susceptibility to Plasmodium28,29, while in another study no relationship was found30. Despite the genetic association, a functional link between activity of para and Anopheles vector competence for Plasmodium has not been shown. We previously performed functional genomics studies of the para gene where we did not detect a phenotype with P. falciparum infection31. As a transmembrane protein, the half-life of the VGSC is known to be long32,33, and lacking a specific antibody against the para protein, one limitation of the study in31 was that we were not able to measure para protein levels after para gene silencing. In addition, a new annotation of the para gene predicted additional splice variants and new transcripts that were not targeted by the silencing procedure used in the 2015 study.
Here, to clarify the relationship between para activity and Plasmodium infection in Anopheles mosquitoes, we designed a silencing strategy that targeted all of the annotated transcripts. We showed that disruption of para activity by RNAi increased Anopheles vector competence for infection for both P. falciparum and P. berghei. We also showed that para silencing led to dysbiosis of the enteric bacterial microbiome and that para activity in limiting Plasmodium infection is microbiota dependent. To our knowledge, this is the first functional study linking para activity with Anopheles vector competence for Plasmodium.
In conducting these experiments, we adhered to the Guide for the Care and Use of Laboratory Animals of the European Union and the French Government. The methods and experiments were approved by the French Ministry of Higher Education, Research and Innovation under the agreement reference (APAFIS #30954-2021040810159180 v1). The study is reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) 34.
The Ngousso strain (An. coluzzii) was colonized in 2006 from mosquitoes collected in Yaoundé, Cameroon. The mosquitoes are insecticide susceptible, and the colony was reared under standard insectary conditions (26 °C, 12:12 L:D, and 70% relative humidity) by the Centre for the Production and infection of Anopheles (CEPIA) of the Institute Pasteur in Paris, France.
The para gene (AGAP004707), located on chromosome arm 2L, is composed of 39 exons, which transcribe 13 different messenger RNAs (mRNA) through alternate splicing. To maximize the efficacy of the gene silencing, two double stranded (ds) RNA for para were designed and synthesized for targeting each two distinct regions along the para gene with the aim to recognize all 13 transcripts. Each of the dsPara and the dsGFP used as control were prepared using the MEGAscript™ RNAi Kit (Invitrogen by Thermo Fisher Scientific) as described in35. Primers bearing the T7 RNA polymerase recognition sequence were designed to amplify ~ 604 (Para1) and ~ 730 (Para2) base pairs regions. All the primer sequences used for the synthesis of the T7-matrix are listed in Table 1. For gene silencing procedure, one to three day old cold-anaesthetized mosquitoes were injected using the NanoJect II (Drummond) either with 1 μg of an equal quantity of the two dsPara (500 ng each per mosquito) or with 1 μg of dsGFP for the control group. dsPara1 and dsPara2 each at 10 μg/μL were used to perform the mix of dsPara used to inject the mosquitoes and 100 nL of this dsRNAmix were injected in each mosquito female. The mosquitoes were injected into the thorax and maintained on sugar water until required for downstream experiments. Efficiency of silencing was determined from cDNA by RT-qPCR, using Total RNA extracted from mosquitoes, 4 days post-injection as described previously in31.
In order to determine if silencing the para gene had any impact on mosquito survival, which might confound data from subsequent experiments, mosquitoes were injected with dsGFP or dsPara as described above, whereafter mosquitoes were maintained on sugar water (10%), ad libitum, for 14 days post injection. Three biological replicates were performed and for each replicate, mosquito mortality was monitored daily with in total 126 mosquitoes in dsGFP and 181 in dsPara. Kaplan–Meier survival curves were prepared and analysed in GraphPad Prism (Version 8.0.0 for Windows, graphpad.com) and the p-values of each replicate were combined using Fisher’s method (Fisher, 1925).
In the absence of antibodies against the hydrophobic membrane channel, and in order to establish whether the gene silencing was effective at the protein level, we prepared modified WHO bioassays using deltamethrin (0.05%) obtained from University Sains Malaysia (Malaysia). The bioassays were carried out 4, 6 and 10 days post injection (~ 1 μg dsGFP or dsPara) and Ngousso mosquitoes were exposed for a period of 7 min, minimal exposure time required for the dsGFP control group to reach close to 100% mortality 24 h after deltamethrin exposure and limit binding saturation of the insecticide on the remaining para receptor in the dsPara background. Depending on mosquito availability and the number injected, between 10 and 25 mosquitoes were exposed per assay, and a minimum of 4 biological repeats were prepared. The mortality rate of the mosquitoes exposed at d-4, d-6 and d-10 was compared between the dsPara and dsGFP treated cohorts 24 h post-exposure, and the results were analysed using the Chi-Square test where p-values were empirically determined using 105 Monte-Carlo permutations. Following independent statistical tests, the p-values from independent tests of significance were combined using the meta-analytical approach of Fisher36 when the direction of change of each independent replicate was concordant (e.g., for each independent replicate the dsPara displayed lower mortality rate than their paired dsGFP control ).
P. falciparum isolate NF54 was cultured using an automated tipper-table implemented in the CEPIA mosquito infection core facility of the Institut Pasteur, as previously described31,35. Briefly, fourteen days after initiating the subculture, gametocyte maturity was tested by exflagellation of microgametes, and parasitemia and numbers of mature male and female gametocytes were counted on Giemsa-stained slides. For experimental infection of mosquitos, 10 ml of medium containing mature gametocytes was centrifuged at 2000 rpm, and the cell pellet was resuspended in an equal volume of normal type AB human serum. The infected erythrocytes were added to fresh erythrocytes in AB human serum, mixed gently, and transferred to a membrane feeder warmed to 37 °C. At 4 days post-dsRNA injection, An. coluzzii mosquitoes were allowed to feed for 15 min, unfed females were discarded, and only fully engorged females were used for further analysis. Blood-fed mosquitoes were maintained at 26 °C and at 70% relative humidity on 10% sucrose solution.
One to three-days old female mosquitoes were injected with dsGFP or dsPara and 4 days post -injection An. coluzzii mosquitoes were fed on mice infected with P. berghei strain PbGFPCON37, which constitutively expresses green fluorescent protein (GFP). The mosquitoes were allowed to feed for 15 min. Unfed females were discarded, and only fully engorged females were maintained at 21 °C (P. berghei) and 70% relative humidity on 10% sucrose.
The same experiments were performed with mosquitoes treated with antibiotic. Briefly, immediately following adult emergence, mosquitoes were maintained on a 10% sucrose solution complemented with Penicillin 62.5 μg/mL, Streptomycin 100 μg/mL and gentamicin 50 μg/mL, and this solution was changed every day. The dsRNA injections and infectious feedings were performed as described above and only fully engorged females were maintained at 21 °C (P. berghei) and 70% relative humidity on 10% sucrose supplemented with the antibiotics.
Phenotypes were calculated from biological replicates of ≥ 30 dissected mosquitoes each and at least three independent replicates were performed for each species of Plasmodium. Mosquito midguts were dissected at 8 d post-infection. For P. falciparum, midguts were stained with 0.4% mercury dibromofluorescein (Sigma) and the number of oocysts was counted by light microscopy. For P. berghei, oocysts were counted by fluorescence microscopy. Phenotypes measured oocyst infection prevalence, which is the proportion of mosquitoes carrying ≥ 1 oocyst among the total number of dissected mosquitos, and oocyst intensity, which is the number of oocysts counted in mosquitoes with ≥ 1 oocyst. Differences in infection prevalence were statistically tested using the Chi-Square test, and analysis of oocyst intensity used the Wilcoxon signed rank non-parametric test. Statistical differences in prevalence and intensity were first tested independently for each independent replicate as described above and p-values were empirically determined using 105 Monte-Carlo permutations. Following independent statistical tests, the p-values from independent tests of significance were combined using the meta-analytical approach of Fisher36 when the direction of change of each independent replicate was concordant (e.g., each independent replicate displayed higher infection prevalence than their paired GFP controls).
Total 16S and Enterobacteriaceae 16S were measured in order to establish whether para transcript silencing had an impact and bacterial microbiome titers. Three-day old female mosquitoes were injected with dsPara or dsGFP, and after 4 days, the midguts were dissected and collected for DNA extraction. For midguts dissection, whole mosquitoes were submerged in 70% ethanol for 3 min, followed by 2 subsequent baths in sterile phosphate buffered saline. Intact (i.e. non-burst or damaged) midguts were collected in sterile Eppendorf tubes on dry ice, and midguts were retained. For each biological replicate, pools of ~ 18 midguts were extracted per treatment group (dsGFP and dsPara) and stored at -80 °C until processing. Three independent biological replicates were prepared. DNA was extracted using the DNeasy PowerSoil Kit (Qiagen). Universal primers (16S_V4q_F: 5′-GTG CCA GCM GCC GCG GTA A-3′ and 16S_V4q_R: 5′-GGA CTA CHV GGG TWT CTA AT-3′) for the V4 region of 16S rDNA were used in real-time qPCR reactions prepared as follows: 10 μl SYBR Green PCR Mix 2X (KAPA Sybr Fast 2X [Sigma-Aldrich]), 0.4 μl forward primer (10 μM), 0.4 μl reverse primer (10 μM), 8.2 μl nuclease-free H2O and 1ul template. Cycling conditions were: 95 °C/1 min followed by 40 cycles of [95 °C/3 s; 60 °C for 30 s; “plate-read”]. Enterobacteriaceae specific PCR was prepared in the same way but using the following primer pair—Entero16S_Fwd: 5′ CGT CGC AAG MMC AAA GAG 3′ and Entero16S_Rev: 5′ TTA CCG CGG CTG CTG GCA C 3′. All reactions were run in triplicate on the CFX Touch (Bio-Rad) and the relative expression was calculated by the 2-delta, delta Ct method after normalization against mosquito ribosomal protein S7 (rpS7). Difference in delta Ct distribution across the independent biological replicates between (dsPara) and (dsGFP) samples was statistically tested using Student t-test.
In order to investigate if there are pathways that might be influenced by loss of function in the voltage gated sodium channel, RNA-seq was carried out using material extracted from dsRNA-treated mosquitoes. Three-day old female females Anopheles coluzzii mosquitoes were injected with dsPara or dsGFP, and after 4 days, total RNAs from three biological replicates were extracted and sequenced. For each biological replicate, pools of ~ 18 mosquitoes were extracted per treatment group. RNA-seq was performed at the University of Minnesota Genomics Center (genomics.umn.edu) as previously described in38. Briefly, using Illumina’s Truseq RNA Sample Preparation Kit (Cat. # RS-122-2001), 1 µg of total RNA was oligo-dT purified using oligo-dT coated magnetic beads, fragmented and then reverse transcribed into cDNA. The cDNA was fragmented, blunt-ended, and ligated to indexed (barcoded) adaptors and amplified using 15 cycles of PCR. Final library size distribution was validated using capillary electrophoresis and quantified using fluorimetry (PicoGreen) and via q-PCR. Indexed libraries were then normalized, pooled and then size selected to 320 bp + /− 5% using Caliper’s XT instrument. Truseq libraries were hybridized to a paired end flow cell and individual fragments were clonally amplified by bridge amplification on the Illumina cBot. Once clustering was complete, the flow cell was loaded on the HiSeq 2000 and sequenced using Illumina’s SBS chemistry. Primary analysis of sequence reads and demultiplexing were done using CASAVA 1.8.2 and de-multiplexed FASTQ files were used for downstream analyses.
Raw read quality was assessed using FastQC version 0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) and MultiQC version 0.739. BWA-MEM version 0.7.7-r441 (https://arxiv.org/abs/1303.3997) was used for alignment against the reference genome of Anopheles gambiae str. PEST version AgamP4 (VectorBase), with default parameters. Genes were counted using featureCounts40 version 1.4.6-p3 with the annotation version AgamP4.7 and the parameters–t mRNA–g ID. Counts data were analyzed using R version 3.3.1 (R Core Team. R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing (2016)) and the Bioconductor package DESeq2 version 1.14.141 with default settings. Genes with adjusted p-values below 0.05 were considered differentially expressed. For each pairwise comparison, raw p-values were adjusted for multiple testing using the Benjamini and Hochberg procedure42. Genes with adjusted p-values below 0.05 were considered differentially expressed.
1-day old female An. coluzzii Ngousso were injected with a mix of two dsRNA (dsPara) to ensure that all 13 alternative splices are targeted and for the control group, dsGFP were injected in mosquitoes. Silencing efficiency was first confirmed by RT-qPCR at 4 days post-injection (Fig. S1). We then measured and compared the survival between dsPara and dsGFP mosquito groups for 14 days post-injection. Mortality curves revealed the absence of significant fitness difference between dsPara and dsGFP treated mosquitoes (p-value = 0.7) (Fig. 1A).
Depletion of para expression reduces deltamethrin sensitivity. (A) Daily mortality rates between dsGFP and dsPara mosquito groups were monitored for 14 days after dsRNA injection. No difference in mortality was found between dsGFP and dsPara. (B) Sensitivity to the pyrethroid deltamethrin (0.05%) was measured in dsGFP and dsPara by exposing the mosquitoes to the impregnated papers for 7 min, and 24h post-exposure, the mortality rate was counted in each group. The 7 min time exposure choice was made to avoid saturation of the para channels by the insecticide, and for which the control group is reaching close to 100% death. Mortality rate was significantly reduced in dsPara background at 4 and 6-days post-exposure, but not at 10 days.
To detect the effect of para depletion mediated by RNAi, and lacking protein-specific antibody, we tested mosquito sensitivity to the insecticide deltamethrin, known to target para receptors. In order to avoid saturation of the remaining receptor after para depletion, mosquitoes were exposed to deltamethrin (0.05%) for 7 min, which was the minimal exposure time that leads to almost 100% mortality in the dsGFP control group. Mosquitoes treated with dsPara showed significantly reduced mortality as compared to the dsGFP group, at d-4 (Chi2 = 46.86, DF = 8, p = 1.63 × 10–07; Table 2 and Fig. 1B) and d-6 (Chi2 = 33.50, DF = 8, p = 4.99 × 10–05; Table 2 and Fig. 1B) post-injection, but not at d-10 post-injection (Chi2 = 5.24, DF = 8, p = 0.73; Table 2 and Fig. 1B), suggesting that at this time point the transcript and protein abundance returned to normal levels (Fig. 1B). This reduction of deltamethrin sensitivity indicated that loss of para function can phenocopy kdr allele resistance.
We then queried the role of the VGSC on Anopheles vector competence for Plasmodium infection. An. coluzzii females were injected with dsPara or dsGFP for the control group and infected with the rodent malaria parasite P. berghei HSP70_GFP. P. berghei infection prevalence (proportion of infected mosquitoes) was significantly higher in the dsPara-treated mosquitoes as compared to the dsGFP controls (Chi2 = 19.252; P-value = 0.0038) (Fig. 2A). The same experiments performed with the human malaria parasite P. falciparum also resulted in an increase of infection prevalence at 8-d post-infection (Chi2 = 18.383; P-value = 0.0053) (Fig. 2B). There was no significant effect of dsPara treatment on infection intensity (the number of parasites in infected mosquitoes) for either P. berghei and P. falciparum (Fig.S2).
Para normal activity limits Plasmodium infection prevalence. Panel (A) and panel (B) show results of infection prevalence (proportion of infected mosquitoes) in An. coluzzii infected with P. berghei and P. falciparum, respectively. Grey color shows proportion of infected mosquitoes. Infection prevalence between dsPara and dsGFP was measured at d-8 post-infection and showed statistically significant increase in dsPara as compared to dsGFP control, for both Plasmodium species. Combined p-value (Fisher method) from the 3 independent biological replicates was indicated in each graph represented in both A and B panels. n = Total number of dissected mosquitoes. Values obtained in every biological replicate are depicted by dots and dotted lines.
We showed above that normal para levels are required for wild-type sensitivity to deltamethrin, as well as resistance to Plasmodium infection, because depletion of para caused decreased mortality due to insecticide and increased infection prevalence, respectively. To understand the mechanisms underlying these functions, we first determined the transcriptional effects of para depletion on mosquitoes, under the same conditions that displayed the functional effects. Transcriptomic analysis using RNAseq analysis of dsPara and dsGFP-treated mosquitoes, at 4 days post-dsRNA injection, indicated that only three genes were significantly differentially regulated between the two conditions (up: 20S proteasome subunit beta 5, NADH dehydrogenase subunit 2, down: annexin B10B; Sup data Table1, Sup Fig.S3). This result indicates that para activity displays a very small transcriptional footprint. This is particularly notable given the influence of para function limiting Plasmodium infection. Many immune genes are known to act against Plasmodium in a transcription-dependent fashion, but a transcriptional effect of para function on immune genes was not detected, indicating the absence of a direct link between para activity and immune gene expression. Thus, neither the effect of para on insecticide mortality nor on Plasmodium infection appear to be based strongly on transcriptional regulation and are likely have other explanations.
Bidirectional communication between gut bacteria and the brain, also referred as “gut-microbiota-brain axis” has been described to play a crucial role in the homeostasis of the gastrointestinal and microbial system of animals43. In addition, in mammals, VGSCs are expressed in primary neurons innervating visceral organs44. The enteric bacterial flora of mosquitoes is known to influence Plasmodium infection outcome45,46,47,48, and the abundance of Enterobacteriaceae is positively correlated with Anopheles vector competence for Plasmodium45. Therefore, we tested whether para depletion could alter the bacterial abundance and composition in the mosquito midgut. The enteric bacterial load was compared between dsGFP and dsPara backgrounds via quantification of 16S ribosomal RNA (16S rRNA) by qPCR. We observed that silencing of para led to increased overall bacterial abundance in the midgut, and, although this was not statistically significant, it specifically caused elevated abundance of taxa in the Enterobacteriacae (Fig. 3A, Fig S4).
Para activity limiting Plasmodium infection is microbiota dependent. (A) 16srDNA relative fold change was measured by qPCR in samples collected from antibiotic treated and non-treated backgrounds. In the non-treated background, Para silencing leads to increased level of the total abundance of An. coluzzii bacteria as compared to the dsGFP control group (dotted line), whereas in antibiotic treated background this effect is abolished. The ratio of the normalized 16srDNA in dsPara/dsGFP was calculated using triplicates. *: Statistical p-value related to the deltaCt distribution between dsPara and dsGFP on 3 independent biological replicates are mentioned in the graph. (B) Infection prevalence between dsPara and dsGFP in antibiotic (ATB)-treated background was measured at d-8post infection and showed no statistical difference between the two conditions dsPara and dsGFP. Combined p-value (Fisher method) from the 3 independent biological replicates was obtained for the infection prevalence. n = Total number of dissected mosquitoes. Values obtained in every biological replicate are depicted by dots and dotted lines.
We then queried the potential role of dsPara-dependent dysbiosis on Plasmodium infection outcome by treating mosquitoes with antibiotics. Para was depleted in An. coluzzii mosquitoes treated with antibiotic after adult emergence, and 4-d post-injection dsGFP and dsPara were infected with P. berghei. We found that antibiotic treatment reversed the effect of para depletion on Anopheles vector competence on both the enteric bacteria microbiota (Fig. 3A) and on P. berghei infection prevalence (Fig. 3B), as compared to the antibiotic-treated dsGFP controls (Fig. 3B). Antibiotic treatment did not affect P. berghei infection intensity (oocyst load) as previously observed in the absence of the treatment (Fig. S5). This result indicates that para activity controlling Plasmodium infection is microbiota dependent.
Here, we show that normal activity of the VGSC para is required for An. coluzzii sensitivity to killing by a pyrethroid insecticide, as well as to limit Plasmodium infection levels. Para activity appears to exert the influence on Plasmodium infection by altering the abundance and composition of the enteric bacterial microbiota. In mammals, excitation and contraction of gastrointestinal system muscles are regulated by membrane depolarization mediated by receptor channels49. Microbiota and gut motility are associated, and gastrointestinal motility can affect the microbiota in terms of amount, location and diversity2,50,51,52. In Aedes mosquitoes, the crop, a sugar-reach meal storage tissue, transfer the sugar into the digestive tract for digestion via peristaltic contractions, which are controlled by stretch-activated channels53,54,55. Stretch-activated channels include VGSC, therefore, we could hypothesize that Para participates in the contraction of the crop, which regulate carbohydrate resource distribution used for the mosquito itself, but these resources are also probably used by the enteric bacteria as an energy resource. As VGSCs participate in muscle contraction, the simplest explanation for our results, including the absence of direct transcriptional effect by para, is that para activity influences the abundance and composition of the enteric bacterial community by acting directly on the contractile motility of the Anopheles digestive tube or on the contractibility of the crop. 16S rRNA gene pyrosequencing of the microbiome between dsPara and dsGFP treated mosquitoes would help identifying bacteria taxa controlled by the VGSC activity, and the ones that are linked with the Plasmodium phenotype outcome. For example, higher titter of Enterobacteriaceae were significantly correlated with P. falciparum infection45.
Binding of the pyrethroid deltamethrin to the VGSC channel fixes the receptor in its active conformation, which leads to constitutive depolarization of the cell membrane and insect death. In Drosophila, para transgenic deletion led to paralysis and decreased lifespan21,22. However, in Anopheles a lethality effect was not detected after transient partial depletion of para transcript by RNAi, likely because residual para protein maintained partial physiological function of the VGSC. We also showed that reduced activity of para by RNAi-mediated silencing increased Plasmodium infection prevalence. Therefore, it could be possible that augmentation of para activity could decrease Anopheles competence for Plasmodium infection. Consistent with their implication in immune-related functions in mammals, VGSC blockers has been proposed as immune modulators for clinical applications56. Therefore, designing an insecticide structural analog or an allosteric ligand that could act as a selective para agonist, but without blocking the receptor channels in a lethal conformation, could potentially generate a tool to control malaria infection of mosquitoes, and therefore limit transmission.
In Anopheles, some observations have associated the occurrence of kdr mutations with increased P. falciparum infection27,29,31, while no relationship was observed in another study. If we assume, that kdr mutations may alter the function and activity of para, and the mosquito vector competence, the presence/absence of correlation found in these studies could be attributed to the composition of the enteric bacteria flora between the mosquitoes from these distinct studies. Previously we identified a genetic locus, named Pfin6, associated with mosquito infection in nature by the human malaria parasite P. falciparum31. In this locus, the haplotype carrying the insecticide resistant kdr allele of para was linked to increased parasite infection prevalence. However, we do not know whether allelic variations in para could alter its physiological function, and thus its role in Plasmodium development, which may explain these distinct results found in the association studies mentioned above and performed from different mosquito genetic backgrounds. In the current study, depletion of para phenocopied the kdr effect by decreasing mortality after deltamethrin exposure and increases Plasmodium infection prevalence. The 1014F mutation underlying the kdr allele was introduced using CRISPR/Cas9 genome editing into an insecticide susceptible Anopheles colony57, but the influence on Plasmodium infection was not tested, which would be interesting. It could also be informative to test the effect of RNAi-mediated para depletion, as done here, but in a mosquito colony fixed for different kdr and related mutations, and test their susceptibility to Plasmodium and effect on enteric bacterial flora.
Transcriptome analysis indicated that para activity influences transcript abundance of only three other genes. Annexin B10B, which was downregulated after para silencing, belongs to a family that can bind to Plasmodium ookinetes, and annexin B10B expression is significantly decreased 48 h after a bloodmeal58,59, but it is not known whether annexin binding influences Plasmodium infection. The two genes induced by para depletion code for the subunit beta 5 of the 20S proteasome, which is responsible for protein degradation, and NADH dehydrogenase subunit 2. Functional analyses of these genes by RNAi silencing and Plasmodium infection challenge will be required to clarify their potential role in Plasmodium infection.
Finally, in the context of the brain-gut-microbiome axis well studied in mammals, identifying new mosquito receptors from the insect nervous system known to control cell excitability such as neurons and muscles, could be tested for their effects on the mosquito microbiome. For example, testing whether others receptor channels or G-protein coupled receptors known to control brain activity could modulate the composition of the mosquito microbiome, and indirectly Anopheles vector competence, is completely unexplored. Functional study through RNAi combined with pharmacological tool using known agonists or antagonists of these receptors could be used to explore this research field and represent molecule candidate for controlling malaria transmission.
Sequence data related to analysis are available under the ArrayExpress accession E-MTAB-12896. URL link: https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-12896?key=8617f493-4d74-4aff-987f-66470012a884.
Gage, P. W., Lamb, G. D. & Wakefield, B. T. Transient and persistent sodium currents in normal and denervated mammalian skeletal muscle. J. Physiol. 418, 427–439. https://doi.org/10.1113/jphysiol.1989.sp017850 (1989).
Article CAS PubMed PubMed Central Google Scholar
Adrian, R. H. & Marshall, M. W. Sodium currents in mammalian muscle. J. Physiol. 268, 223–250. https://doi.org/10.1113/jphysiol.1977.sp011855 (1977).
Article CAS PubMed PubMed Central Google Scholar
Hernandez-Ochoa, E. O. & Schneider, M. F. Voltage sensing mechanism in skeletal muscle excitation-contraction coupling: Coming of age or midlife crisis?. Skelet. Muscle 8, 22. https://doi.org/10.1186/s13395-018-0167-9 (2018).
Article CAS PubMed PubMed Central Google Scholar
Besson, P. et al. How do voltage-gated sodium channels enhance migration and invasiveness in cancer cells?. Biochim. Biophys. Acta 1848, 2493–2501. https://doi.org/10.1016/j.bbamem.2015.04.013 (2015).
Article CAS PubMed Google Scholar
Wang, Z. et al. Voltage-gated sodium channels are involved in cognitive impairments in Parkinson’s disease-like rats. Neuroscience 418, 231–243. https://doi.org/10.1016/j.neuroscience.2019.08.024 (2019).
Article CAS PubMed Google Scholar
de Lera Ruiz, M. & Kraus, R. L. Voltage-gated sodium channels: Structure, function, pharmacology, and clinical indications. J. Med. Chem. 58, 7093–7118. https://doi.org/10.1021/jm501981g (2015).
Article CAS PubMed Google Scholar
Stehle, S. & Schulz, R. Agricultural insecticides threaten surface waters at the global scale. Proc. Natl. Acad. Sci. USA 112, 5750–5755. https://doi.org/10.1073/pnas.1500232112 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Kuivila, K. M. et al. Occurrence and potential sources of pyrethroid insecticides in stream sediments from seven U.S. metropolitan areas. Environ. Sci. Technol. 46, 4297–4303. https://doi.org/10.1021/es2044882 (2012).
Article ADS CAS PubMed Google Scholar
Farag, M. R. et al. An overview on the potential hazards of pyrethroid insecticides in fish, with special emphasis on cypermethrin toxicity. Animals (Basel) https://doi.org/10.3390/ani11071880 (2021).
Article PubMed PubMed Central Google Scholar
Davies, T. G., Field, L. M., Usherwood, P. N. & Williamson, M. S. DDT, pyrethrins, pyrethroids and insect sodium channels. IUBMB Life 59, 151–162. https://doi.org/10.1080/15216540701352042 (2007).
Article CAS PubMed Google Scholar
Field, L. M., Emyr Davies, T. G., O’Reilly, A. O., Williamson, M. S. & Wallace, B. A. Voltage-gated sodium channels as targets for pyrethroid insecticides. Eur. Biophys. J. 46, 675–679. https://doi.org/10.1007/s00249-016-1195-1 (2017).
Article CAS PubMed PubMed Central Google Scholar
World Health Organization. World Malaria Report 2021, <https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2021> (2021).
Enayati, A. A., Lengeler, C., Erlanger, T. & Hemingway, J. Field evaluation of a recombinant glutathione S-transferase-based pyrethroid quantification assay. Trans. R Soc. Trop. Med. Hyg. 99, 369–378. https://doi.org/10.1016/j.trstmh.2004.09.009 (2005).
Article CAS PubMed Google Scholar
Enayati, A. A., Ranson, H. & Hemingway, J. Insect glutathione transferases and insecticide resistance. Insect Mol. Biol. 14, 3–8. https://doi.org/10.1111/j.1365-2583.2004.00529.x (2005).
Article CAS PubMed Google Scholar
Liu, N. Insecticide resistance in mosquitoes: Impact, mechanisms, and research directions. Annu. Rev. Entomol. 60, 537–559. https://doi.org/10.1146/annurev-ento-010814-020828 (2015).
Article CAS PubMed Google Scholar
Liu, N., Li, M., Gong, Y., Liu, F. & Li, T. Cytochrome P450s–Their expression, regulation, and role in insecticide resistance. Pestic Biochem. Physiol. 120, 77–81. https://doi.org/10.1016/j.pestbp.2015.01.006 (2015).
Article CAS PubMed Google Scholar
Muller, P. et al. Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 4, e1000286. https://doi.org/10.1371/journal.pgen.1000286 (2008).
Article CAS PubMed PubMed Central Google Scholar
Ranson, H. et al. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect. Mol. Biol. 9, 491–497. https://doi.org/10.1046/j.1365-2583.2000.00209.x (2000).
Article CAS PubMed Google Scholar
Martinez-Torres, D. et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect. Mol. Biol. 7, 179–184. https://doi.org/10.1046/j.1365-2583.1998.72062.x (1998).
Article CAS PubMed Google Scholar
Balabanidou, V., Grigoraki, L. & Vontas, J. Insect cuticle: A critical determinant of insecticide resistance. Curr. Opin. Insect. Sci. 27, 68–74. https://doi.org/10.1016/j.cois.2018.03.001 (2018).
Article PubMed Google Scholar
Loughney, K., Kreber, R. & Ganetzky, B. Molecular analysis of the para locus, a sodium channel gene in Drosophila. Cell 58, 1143–1154. https://doi.org/10.1016/0092-8674(89)90512-6 (1989).
Article CAS PubMed Google Scholar
Suzuki, D. T., Grigliatti, T. & Williamson, R. Temperature-sensitive mutations in Drosophila melanogaster VII. A mutation (para-ts) causing reversible adult paralysis. Proc. Natl. Acad. Sci. USA 68, 890–893. https://doi.org/10.1073/pnas.68.5.890 (1971).
Article ADS CAS PubMed PubMed Central Google Scholar
Clarkson, C. S. et al. The genetic architecture of target-site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii. Mol. Ecol. 30, 5303–5317. https://doi.org/10.1111/mec.15845 (2021).
Article CAS PubMed PubMed Central Google Scholar
Edi, A. V. C. et al. First detection of N1575Y mutation in pyrethroid resistant Anopheles gambiae in Southern Cote d’Ivoire. Wellcome Open Res. 2, 71. https://doi.org/10.12688/wellcomeopenres.12246.1 (2017).
Article CAS PubMed PubMed Central Google Scholar
Jones, C. M. et al. Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles gambiae. Proc. Natl. Acad. Sci. USA 109, 6614–6619. https://doi.org/10.1073/pnas.1201475109 (2012).
Article ADS PubMed PubMed Central Google Scholar
Adams, K. L. et al. Selection for insecticide resistance can promote Plasmodium falciparum infection in Anopheles. PLoS Pathog. 19, e1011448. https://doi.org/10.1371/journal.ppat.1011448 (2023).
Article CAS PubMed PubMed Central Google Scholar
Alout, H. et al. Insecticide resistance alleles affect vector competence of Anopheles gambiae s.s. for Plasmodium falciparum field isolates. PLoS ONE 8, e63849. https://doi.org/10.1371/journal.pone.0063849 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Hanemaaijer, M. J. 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. https://doi.org/10.1186/s12936-019-2759-1 (2019).
Article PubMed PubMed Central Google Scholar
Ndiath, M. O. 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. 13, 340. https://doi.org/10.1186/1475-2875-13-340 (2014).
Article CAS PubMed PubMed Central Google Scholar
Traore, A. et al. Anopheline species composition and the 1014F-genotype in different ecological settings of Burkina Faso in relation to malaria transmission. Malar. J. 18, 165. https://doi.org/10.1186/s12936-019-2789-8 (2019).
Article ADS PubMed PubMed Central Google Scholar
Mitri, C. et al. The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: Genetic correlation and functional testing. Malar. J. 14, 391. https://doi.org/10.1186/s12936-015-0924-8 (2015).
Article CAS PubMed PubMed Central Google Scholar
Sherman, S. J., Chrivia, J. & Catterall, W. A. Cyclic adenosine 3’:5’-monophosphate and cytosolic calcium exert opposing effects on biosynthesis of tetrodotoxin-sensitive sodium channels in rat muscle cells. J. Neurosci. 5, 1570–1576 (1985).
Article CAS PubMed PubMed Central Google Scholar
Waechter, C. J., Schmidt, J. W. & Catterall, W. A. Glycosylation is required for maintenance of functional sodium channels in neuroblastoma cells. J. Biol. Chem. 258, 5117–5123 (1983).
Article CAS PubMed Google Scholar
du Percie Sert, N. et al. The ARRIVE guidelines 2.0: Updated guidelines for reporting animal research. PLoS Biol 18, e3000410. https://doi.org/10.1371/journal.pbio.3000410 (2020).
Article CAS Google Scholar
Mitri, C. et al. Fine pathogen discrimination within the APL1 gene family protects Anopheles gambiae against human and rodent malaria species. PLoS Pathog. 5, e1000576. https://doi.org/10.1371/journal.ppat.1000576 (2009).
Article CAS PubMed PubMed Central Google Scholar
Fisher, R. A. Statistical methods for research workers. (Oliver and Boyd, 1925).
Franke-Fayard, B. et al. A Plasmodium berghei reference line that constitutively expresses GFP at a high level throughout the complete life cycle. Mol. Biochem. Parasitol. 137, 23–33. https://doi.org/10.1016/j.molbiopara.2004.04.007 (2004).
Article CAS PubMed Google Scholar
Dieme, C. et al. Exposure of Anopheles mosquitoes to trypanosomes reduces reproductive fitness and enhances susceptibility to Plasmodium. PLoS Negl. Trop Dis. 14, e0008059. https://doi.org/10.1371/journal.pntd.0008059 (2020).
Article CAS PubMed PubMed Central Google Scholar
Ewels, P., Magnusson, M., Lundin, S. & Kaller, M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048. https://doi.org/10.1093/bioinformatics/btw354 (2016).
Article CAS PubMed PubMed Central Google Scholar
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. https://doi.org/10.1093/bioinformatics/btt656 (2014).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome. Biol. 15, 550. https://doi.org/10.1186/s13059-014-0550-8 (2014).
Article CAS PubMed PubMed Central Google Scholar
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate–a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x (1995).
Article MathSciNet MATH Google Scholar
Morais, L. H., Schreiber, H. L. T. & Mazmanian, S. K. The gut microbiota-brain axis in behaviour and brain disorders. Nat. Rev. Microbiol. 19, 241–255. https://doi.org/10.1038/s41579-020-00460-0 (2021).
Article CAS PubMed Google Scholar
Qi, F. H., Zhou, Y. L. & Xu, G. Y. Targeting voltage-gated sodium channels for treatment for chronic visceral pain. World J. Gastroenterol. 17, 2357–2364. https://doi.org/10.3748/wjg.v17.i19.2357 (2011).
Article PubMed PubMed Central Google Scholar
Boissiere, A. et al. Midgut microbiota of the malaria mosquito vector Anopheles gambiae and interactions with Plasmodium falciparum infection. PLoS Pathog. 8, e1002742. https://doi.org/10.1371/journal.ppat.1002742 (2012).
Article CAS PubMed PubMed Central Google Scholar
Dennison, N. J., Jupatanakul, N. & Dimopoulos, G. The mosquito microbiota influences vector competence for human pathogens. Curr. Opin. Insect. Sci. 3, 6–13. https://doi.org/10.1016/j.cois.2014.07.004 (2014).
Article PubMed PubMed Central Google Scholar
Gendrin, M. et al. Antibiotics in ingested human blood affect the mosquito microbiota and capacity to transmit malaria. Nat. Commun. 6, 5921. https://doi.org/10.1038/ncomms6921 (2015).
Article ADS PubMed Google Scholar
Tchioffo, M. T. et al. Modulation of malaria infection in Anopheles gambiae mosquitoes exposed to natural midgut bacteria. PLoS ONE 8, e81663. https://doi.org/10.1371/journal.pone.0081663 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Waclawikova, B., Codutti, A., Alim, K. & El Aidy, S. Gut microbiota-motility interregulation: Insights from in vivo, ex vivo and in silico studies. Gut. Microbes. 14, 1997296. https://doi.org/10.1080/19490976.2021.1997296 (2022).
Article CAS PubMed PubMed Central Google Scholar
Vantrappen, G., Janssens, J., Hellemans, J. & Ghoos, Y. The interdigestive motor complex of normal subjects and patients with bacterial overgrowth of the small intestine. J. Clin. Invest. 59, 1158–1166. https://doi.org/10.1172/JCI108740 (1977).
Article CAS PubMed PubMed Central Google Scholar
Dukowicz, A. C., Lacy, B. E. & Levine, G. M. Small intestinal bacterial overgrowth: A comprehensive review. Gastroenterol. Hepatol. (N. Y.) 3, 112–122 (2007).
PubMed Google Scholar
Guarino, M. P., Cicala, M., Putignani, L. & Severi, C. Gastrointestinal neuromuscular apparatus: An underestimated target of gut microbiota. World J. Gastroenterol. 22, 9871–9879. https://doi.org/10.3748/wjg.v22.i45.9871 (2016).
Article CAS PubMed PubMed Central Google Scholar
Stoffolano, J. G. Jr. & Haselton, A. T. The adult Dipteran crop: A unique and overlooked organ. Annu. Rev. Entomol. 58, 205–225. https://doi.org/10.1146/annurev-ento-120811-153653 (2013).
Article CAS PubMed Google Scholar
Thomson, A. J. & Holling, C. S. Experimental component analysis of feeding rate of blowfly phormia-regina (Diptera-Calliphoridae). Can. Entomol. 107, 167–173. https://doi.org/10.4039/Ent107167-2 (1975).
Article Google Scholar
Villegas, L. E. M., Radl, J., Dimopoulos, G. & Short, S. M. Bacterial communities of Aedes aegypti mosquitoes differ between crop and midgut tissues. PLoS Negl. Trop. Dis. 17, e0011218. https://doi.org/10.1371/journal.pntd.0011218 (2023).
Article CAS PubMed PubMed Central Google Scholar
Roselli, F., Livrea, P. & Jirillo, E. Voltage-gated sodium channel blockers as immunomodulators. Recent Pat. CNS Drug Discov. 1, 83–91. https://doi.org/10.2174/157488906775245255 (2006).
Article CAS PubMed Google Scholar
Grigoraki, L. et al. CRISPR/Cas9 modified. An gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes. PLoS Genet. 17, e1009556. https://doi.org/10.1371/journal.pgen.1009556 (2021).
Article CAS PubMed PubMed Central Google Scholar
Kotsyfakis, M. et al. Plasmodium berghei ookinetes bind to Anopheles gambiae and Drosophila melanogaster annexins. Mol. Microbiol. 57, 171–179. https://doi.org/10.1111/j.1365-2958.2005.04664.x (2005).
Article CAS PubMed Google Scholar
Kotsyfakis, M., Vontas, J., Siden-Kiamos, I. & Louis, C. The annexin gene family in the malaria mosquito Anopheles gambiae. Insect. Mol. Biol. 14, 555–562. https://doi.org/10.1111/j.1365-2583.2005.00586.x (2005).
Article CAS PubMed Google Scholar
Download references
This study received financial support to KV from the European Commission Horizon 2020 Infrastructures #731060 Infravec2; the Agence Nationale de la Recherche, #ANR-19-CE35-0004 ArboVec; National Institutes of Health, NIAID #AI145999; and French Laboratoire d’Excellence “Integrative Biology of Emerging Infectious Diseases” #ANR-10-LABX-62-IBEID.
Genetics and Genomics of Insect Vectors Unit, Department of Parasites and Insect Vectors, CNRS, Institut Pasteur, UMR2000, Université de Paris, 75015, Paris, France
Luisa Nardini, Emma Brito-Fravallo, Adrien Pain, Corinne Genève, Kenneth D. Vernick & Christian Mitri
Center of Bioinformatics, Biostatistics and Integrative Biology, Institut Pasteur, Université de Paris, 75015, Paris, France
Pascal Campagne
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
L.N., E.B.F., C.G., K.V. and C.M. designed the research. L.N., E.B.F., C.G., K.V. and C.M. performed the research. L.N., E.B.F., C.G., A.P., K.V. and C.M. analyzed the data. L.N., K.V. and C.M. wrote the manuscript.
Correspondence to Christian Mitri.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Nardini, L., Brito-Fravallo, E., Campagne, P. et al. The voltage-gated sodium channel, para, limits Anopheles coluzzii vector competence in a microbiota dependent manner. Sci Rep 13, 14572 (2023). https://doi.org/10.1038/s41598-023-40432-x
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-40432-x
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.