- Open Access
The ecdysone receptor regulates several key physiological factors in Anopheles funestus
Malaria Journal volume 21, Article number: 97 (2022)
Malaria is a devastating disease, transmitted by female Anopheles mosquitoes infected with Plasmodium parasites. Current insecticide-based strategies exist to control the spread of malaria by targeting vectors. However, the increase in insecticide resistance in vector populations hinder the efficacy of these methods. It is, therefore, essential to develop novel vector control methods that efficiently target transmission reducing factors such as vector density and competence. A possible vector control candidate gene, the ecdysone receptor, regulates longevity, reproduction, immunity and other physiological processes in several insects, including malaria vectors. Anopheles funestus is a prominent vector in sub-Saharan Africa, however, the function of the ecdysone receptor in this mosquito has not previously been studied. This study aimed to determine if the ecdysone receptor depletion impacts An. funestus longevity, reproduction and susceptibility to Plasmodium falciparum infection.
RNA interference was used to reduce ecdysone receptor expression levels in An. funestus females and investigate how the above-mentioned phenotypes are influenced. Additionally, the expression levels of the ecdysone receptor, and reproduction genes lipophorin and vitellogenin receptor as well as the immune gene, leucine rich immune molecule 9 were determined in ecdysone receptor-depleted mosquitoes using quantitative polymerase chain reaction.
Ecdysone receptor-depleted mosquitoes had a shorter lifespan, impaired oogenesis, were less fertile, and had reduced P. falciparum infection intensity.
Overall, this study provides the first experimental evidence that supports ecdysone receptor as a potential target in the development of vector control measures targeting An. funestus.
Malaria is a severe disease transmitted by female Anopheles mosquitoes infected with protozoan Plasmodium parasites. In 2019, a concerning 229 million malaria cases and 409,000 deaths were recorded . Malaria is endemic to tropical and sub-tropical regions of the world. Africa is severely affected, with more than 90% of cases and mortalities occurring on the continent annually . The main malaria vectors in sub-Saharan Africa are either members of the Anopheles gambiae complex or the Anopheles funestus group . The An. funestus group comprises eleven African species, but contains a single significant malaria vector, An. funestus [3,4,5,6]. Due to its highly anthropophilic (preference for a human host) and endophilic (preference to rest indoors) nature, An. funestus is one of the main vectors in sub-Saharan Africa [2, 7].
Efforts to eradicate malaria have been implemented by the World Health Organization (WHO) in the form of The Global Technical Strategy for Malaria (GTS) 2016–2030 . The GTS milestones aim to eliminate at least 90% of malaria by 2030 especially in countries that were most affected in the year 2015 . Unfortunately, the 2020 milestones were not achieved due to several factors , emphasizing the need for novel and effective control strategies. Malaria is controlled by adopting strategies targeting either the Plasmodium parasite or Anopheles vector.
Anopheles mosquitoes are primarily controlled using two insecticide-based control interventions: indoor residual spraying (IRS) and long-lasting insecticidal nets (LLIN). These methods are widely employed and their popularity is attributed to their efficiency, cost-effectiveness and ease of implementation . In Africa, LLINs and IRS have proved to be very effective forms of vector control, preventing malaria transmission by 68% and 13% respectively between the years 2000–2015 . Unfortunately, vectors have undergone adaptations that have hindered efforts in malaria eradication. Specifically, An. funestus mosquitoes have developed genetic insecticide resistance mechanisms [12,13,14,15,16,17], physical changes that limit insecticide uptake  and behavioural adaptations, such as diurnal and outdoor feeding [19, 20], to evade the current insecticide based interventions. Intensified efforts in vector control are, therefore, required to achieve the WHO targets and potentially eradicate malaria in endemic areas. One solution could be to explore the development of alternative control methods.
There is a growing interest in genetically-based control methods due to the availability of Anopheles genomes [21, 22] and development of molecular techniques, such as gene editing by clustered regularly interspaced short palindromic repeats (CRISPR) Cas9 technology  among others. The identification of potential genes involved in vector susceptibility to Plasmodium or other factors, such as vector longevity and reproduction would be ideal for genetically-based control interventions to reduce malaria transmission. One such gene is the ecdysone receptor (EcR). The EcR is a nuclear receptor that functions as a ligand controlled transcription factor . The EcR binds to the Ultraspiracle protein (USP) and the 20-hydroxyecdysone (20E) hormone to form the EcR-USP complex [25,26,27]. Together the functional EcR-USP complex binds to a region of DNA either in the form of an inverted or direct repeat, termed the ecdysone response element (EcRE) to activate transcription of the target gene identified as an ‘’early gene’’ [24, 28, 29]. Upon activation, these early genes subsequently regulate the expression of several late genes. This process is termed the 20E signalling pathway. Importantly, EcR via the 20E pathway regulates several physiological functions that ultimately affect malaria transmission in vector species .
Vector longevity, reproduction and susceptibility to Plasmodium are all factors that affect the density of vectors and/or their competence to transmit malaria . The 20E signaling pathway influences the survival of vectors and treatment of vectors with 20E agonists reduced their longevity [32, 33]. The role that EcR depletion plays in the 20E pathway, i.e. impairing the pathway, by regulating vector longevity however has not yet been characterized in An. funestus. Manipulation of this hormone receptor has provided invaluable insights about the 20E signaling pathway in the mosquito. Egg development is hindered and a decrease in the number of eggs developed is observed in EcR depleted Aedes aegypti and An. gambiae mosquitoes due to decreased 20E signaling [34, 35]. Two genes, Lipophorin (Lp) and Vitellogenin receptor (VgR), regulated by the 20E pathway, largely regulate the reproductive processes vitellogenesis and oogenesis [35,36,37,38,39] therefore to determine if EcR regulates these genes through 20E signaling is of interest. Furthermore, EcR depleted An. gambiae also developed fewer P. falciparum oocysts in comparison to the controls . Implying that functional EcR is required in the 20E signaling pathway, to induce a successful immune response in vectors against Plasmodium. An important gene associated with Plasmodium immunity, leucine rich immune molecule 9 (LRIM 9), regulated by EcR in An. gambiae was discovered to exhibit immunity against Plasmodium berghei in a previous study . The current study therefore aimed to (i) characterize the role of EcR in longevity, (ii) explore its function in fertility, oogenesis and regulation of Lp and VgR genes and (iii) to determine its role in LRIM 9 regulation and immunity against P. falciparum, all in an important African malaria vector, An. funestus.
A laboratory strain of An. funestus mosquitoes originating from Mozambique (FUMOZ) was used in this study. Female FUMOZ mosquitoes were used for all data collection purposes and were aged accordingly. The FUMOZ strain was reared in the Maureen Coetzee insectary at the Wits Research Institute for Malaria (WRIM) under standard insectary conditions at ± 26 °C with ± 80% humidity and a 12:12 day:night cycle .
Gametocytes were produced as per Reader et al. . Briefly, PfNF54 gametocytes were induced from a highly synchronized (> 95%) asexual ring-stage population at 0.5% parasitaemia and 6% haematocrit. Medium was aspirated daily and replaced with fresh medium pre-warmed to 37 °C. Gametocyte health was monitored with microscopy of Giemsa-stained smears. Mature PfNF54 gametocytes (> 98% stage V, 1.5−2.5% gametocytaemia) were used to infect An. funestus.
RNA extraction, DNase I treatment, and cDNA synthesis
RNA was extracted using the TRIzol™ Plus RNA Purification Kit (12183555, Invitrogen, USA) according to the manufacturer's instructions. RNA quality control consisted of spectrophotometry using the Nanodrop One (ND-ONE-W, Thermo Scientific, MA, USA) to determine RNA concentration and purity as well as agarose gel electrophoresis to confirm RNA integrity. A DNase treatment with the TURBO DNA-free™ Kit (AM1907, Ambion, TX, USA) was conducted according to manufacturer’s instructions to remove any contaminating genomic DNA carried over from RNA extraction. Once RNA was free of genomic DNA, First strand cDNA synthesis was conducted with the iScript™ cDNA Synthesis Kit (1708891, BioRad, CA, USA).
Each 1X qPCR reaction consisted of 5 µl of IQ™ SYBR super-mix (1708880, Bio-Rad, CA, USA), 300 nM of each primer, 1 µl of cDNA and nuclease free water to make up a 10 µl reaction. Cycling conditions consisted of an initial denaturation at 94 °C for 3 min, 35 cycles of denaturation, annealing and extension at 94 °C for 20 s, annealing temperature determined by optimization for 25 s and 72 °C for 30 s respectively with a final extension step at 72 °C for 10 min. A melt peak analysis and a no-template control excluding cDNA were included to assess for specificity and contamination, respectively. Reference genes RPS7 and RPL19 were used to normalize samples. Primer sequences are listed in Table 1.
The starting material for the production of dsRNA molecules consisted of a plasmid containing GFP or for EcR, cDNA prepared from RNA using the iScript™ cDNA synthesis kit and subject to PCR to obtain a 484 bp amplicon. Total RNA for this purpose was extracted from female An. funestus at 24 h post-blood meal, since this is the peak of EcR expression . Subsequently, to add T7 promoter sequences to the 5′ end of the GFP clone and EcR template, a PCR was conducted using the T7-containing RNAi primers and the GFP or EcR template. The T7-containing 544 bp GFP and 530 bp EcR templates were concentrated and purified with ethanol precipitation to a final concentration of 1 µg/µl. Double stranded RNA molecules homologous to the EcR and GFP templates were synthesised with the MEGAscript® RNAi Kit (AM1626, Thermo Scientific, MA, USA) according to the manufacturer’s instructions. Specific dsRNA products were then concentrated by precipitation to a final concentration of 10 µg/µl. All primers used are described in Table 1.
To deliver dsRNA molecules to mosquitoes, nanoinjection was conducted using the Nanoject II (3-000-204, Drummond, AL, USA). All mosquitoes were fed on a 10% (w/v) sucrose solution prior to nanoinjection. Mosquitoes were cold anaesthetized prior to nanoinjection in the thorax with approximately 500 ng dsRNA. After nanoinjection, mosquitoes were placed in recovery cages where they were provided with a 10% (w/v) sucrose solution and kept under standard insectary conditions to be used for further assays. An uninjected control was included in the biological assays to determine if the phenotypes observed occurred as a result of the nanoinjection procedure.
One day old female An. funestus mosquitoes were administered with either dsEcR or dsGFP (control) according to the nanoinjection method described previously. Post injection mosquitoes were provided with a 10% sucrose solution and maintained under standard insectary conditions. Mortality among mosquitoes was recorded each day until all the mosquitoes had died. Three biological replicates were conducted. Knockdown of EcR was confirmed with qPCR. Approximately fifteen mosquitoes from each biological replicate were randomly selected from each treatment group 48 h post injection and subject to RNA extraction, RNA quality control, DNase treatment, cDNA synthesis and qPCR as described previously. Knockdown of EcR was evaluated at the 24, 48 and 72 h time points to determine the most suitable time to confirm dsEcR treatment was successful (Additional file 1: Fig. S1).
Forty newly emerged males and forty newly emerged females were combined into 17.5 cm3 BugDorm cages (BD41515, MegaView Science, Taiwan, China) per replicate and allowed to mate for 12 days to achieve the optimal mating success rate (Additional file 2: Fig. S2). Six biological replicates were conducted. Mosquitoes were maintained on a 10% sucrose solution under standard insectary conditions for the duration of this experiment. When females were 10 days old they were isolated prior to dsRNA treatment. After injection, the females were placed back with males for the remaining two days to fulfil the 12 day mating period. On the eleventh day, the 10% sucrose solution was removed and replaced with distilled water to encourage blood feeding. On the twelfth day, fifteen females from each treatment group were removed for RNA extraction and qPCR to confirm knockdown. Simultaneously, females were blood fed for 1 h on bovine blood with the Hemotek artificial membrane feeding system. Blood fed females were subsequently removed and placed into individual 250 ml oviposition cups containing approximately 10 ml of distilled water and lined with filter paper over a black background. The number of eggs laid by mated females were counted by microscopy to determine fecundity. Females that had oviposited were dissected after oviposition and their spermathecae observed by microscopy at 200× magnification to determine mating status; virgin females were subsequently excluded from the study. If mosquitoes had not oviposited seven days post blood meal, they were dissected to determine their mating status and observe if eggs were retained in ovaries.
Plasmodium infection assay
Anopheles funestus females aged seven to ten days old (to allow for optimal mating and encourage blood feeding) were injected with 10 µg/µl of dsEcR or dsGFP as described previously. Ninety mosquitoes were included per replicate and a total of eight biological replicates were conducted. Subsequent to injection, the 10% sucrose solution was removed and replaced with distilled water to encourage blood feeding. Concurrently, mosquitoes were isolated for RNA extraction, DNase 1 treatment, cDNA synthesis and qPCR. Twenty-four hours after nanoinjection, mosquitoes were offered a PfNF54 infected blood meal (> 98% stage V gametocytes, 1.5−2.5% gametocytaemia, 50% (v/v) A + male human serum, Interstate blood bank Inc, Memphis, Tennessee, USA) for 40 min using a glass feeder. Unfed mosquitoes were discarded while fed mosquitoes were maintained on a 10% sucrose solution for eight days post feeding. After this time, mosquitoes were aspirated into ethanol to immobilize them and subsequently transferred to 1× PBS. Mosquito midguts were dissected 8 days post infection and stained on a microscope slide using 0.1% mercurochrome (M7011, Sigma, MO, USA) . Midguts were subsequently viewed under a compound microscope between 20 and 40× magnification and the intensity and prevalence of oocysts in each midgut was counted and recorded.
All qPCR data were analysed using relative gene expression analysis as per 2−∆∆Ct method . The 2−∆∆Ct values were subsequently log transformed and statistical significance was determined using an unpaired student’s t-test . The Kaplan Meier survival curve was used for longevity analysed using the Log-Rank test to calculate statistical significance (GraphPad Prism 8). Statistical analysis for the eggs oviposited, fertility and P. falciparum oocyst intensity were evaluated for normality using the Shapiro–Wilk test. If data were not normally distributed the Mann–Whitney test was used, otherwise an unpaired students t-test was used. To determine statistical significance for the number of ovaries containing mature eggs compared to those with immature follicles, a two tailed Fisher’s exact test was used. Plasmodium falciparum oocyst prevalence data was analysed using a Chi-squared test. For all analyses, p values < 0.05 were considered to be statistically significant.
Injection of dsEcR reduces EcR expression but does not influence Lp, VgR and LRIM 9 expression in An. funestus
To assess different phenotypes associated with EcR depletion, An. funestus females were injected with dsRNA targeting EcR (dsEcR) alongside dsGFP as a control. Relative expression, as determined by qPCR, was used to confirm gene silencing. EcR was significantly downregulated (t(6) = 4.3, p = 0.0051, unpaired student’s t-test) with a mean expression ± SEM of − 2.23 ± 0.32 fold in dsEcR injected samples compared to − 1.50e−006 ± 0.40 fold confirming that the dsEcR injection significantly reduced EcR transcription (Fig. 1A).
Likewise, the expression of several other transcripts Lp, VgR and LRIM 9 previously demonstrated to be involved in mosquito reproduction, or susceptibility to Plasmodium species [38, 40, 47], were measured. In dsEcR injected females, Lp expression was downregulated to a mean ± SEM of − 1.58 ± 0.12 fold compared to 1.33e−006 ± 0.66 fold in dsGFP injected females (Fig. 1B). This difference was however not statistically significant (t(4) = 2.35, p = 0.0785, unpaired student’s t-test). Expression of LRIM 9 was reduced to a mean ± SEM of − 0.56 ± 0.19 fold in dsEcR injected females compared to dsGFP injected females with a mean ± SEM of − 1.00e−006 ± 0.08 fold, however no statistical significance was observed (t(4) = 2.703, p = 0.0539, unpaired student’s t-test) (Fig. 1C). Moreover, VgR expression was depleted to a mean ± SEM of − 0.84 ± 0.03 fold in dsEcR treated females compared to − 6.67e−007 ± 0.36 fold in dsGFP treated females (Fig. 1D). This difference was not statistically significant (t(4) = 2.327, p = 0.0805, unpaired student’s t-test). No significant difference was observed in Lp, LRIM 9 nor VgR, suggesting depletion of EcR does not affect expression of these genes in this study.
EcR depletion results in decreased longevity in An. funestus
To determine the effect of EcR depletion on An. funestus longevity, mortality was recorded daily in dsEcR and dsGFP treated females until all mosquitoes had succumbed. Mosquitoes injected with dsEcR did not survive past 20 days whereas those injected with dsGFP survived up to 37 days. There was a statistically significant difference in the probability of survival between the two groups (χ2(1, N= 240) = 30.18, p < 0.0001, Log rank test). The median survival (time taken to reach a survival of 50%) of dsGFP females was 14 (95% CI of ratio: 2.96–5.41) days compared to 3.5 (95% CI of ratio: 0.19 to 0.34) days for females injected with dsEcR. This translates to a 54% decrease in longevity in dsEcR injected mosquitoes (Fig. 2). No significant difference (χ2(1, N= 240) = 2.011, p = 0.1561¸ Log rank test) was observed between the dsGFP and uninjected control groups, suggesting that the nanoinjection procedure did not influence longevity (Additional file 4: Table S1).
EcR depletion decreases reproductive success in An. funestus.
To determine EcRs effect on fecundity, the number of eggs oviposited by mated females was counted in dsEcR and dsGFP treated females. Although a total number of 157 females were induced to lay eggs, the number of females that oviposited were low. Altogether, 9 dsEcR treated females and 19 dsGFP treated females oviposited. No significant difference was observed in the blood feeding rates between dsEcR and dsGFP groups (Additional file 3: Fig. S3). The mean number of eggs per female ± SEM oviposited in dsEcR treated females was 36.11 ± 4.01 compared to 48.11 ± 4.12 in dsGFP treated females (Fig. 3A). This decrease however, was not significant (t(26) = 1.813, p = 0.0814, unpaired student’s t-test), suggesting that EcR depletion does not affect the number of eggs oviposited. Additionally no statistical difference (t(91) = 0.1522, p = 0.8793, unpaired student’s t-test) was observed between the dsGFP and uninjected control groups, suggesting that the nanoinjection procedure had no effect on the number of eggs oviposited (Additional file 4: Table S1). Those mosquitoes that had not oviposited were dissected to observe if eggs were retained in the ovaries. Surprisingly, the majority of ovaries of dsEcR injected mated An. funestus females closely resembled that of virgin An. funestus females, containing immature and undifferentiated oocytes whereas the majority of dsGFP injected females ovaries contained mature eggs (Fig. 3B). In dsEcR injected females, 32% of mated females developed eggs whereas 89% of dsGFP injected females developed eggs. This difference observed was statistically significantly different (p < 0.0001, Fisher’s exact test) (Fig. 3C). This suggests that dsEcR injected An. funestus females typically develop fewer eggs compared to controls. Together these results substantiate that the depletion of EcR had a deleterious effect on oogenesis by preventing the development of oocytes into mature eggs in An. funestus. No significant difference (p > 0.9999, Fisher’s exact test) was observed between the dsGFP and uninjected control groups, indicating the nanoinjection procedure did not affect the oviposition process (Additional file 4: Table S1). Although very few dsEcR injected females oviposited eggs (n = 9), fertility was also compared between the dsEcR and dsGFP treated mosquitoes by monitoring the hatching of the eggs oviposited. The median percentage fertility (25–75% percentile) in dsEcR treated mated females was significantly lower at 69% (55–80%) compared to 86% (76–93%) in dsGFP treated mated females (t(26) = 3.169, p = 0.0039, unpaired student’s t-test) (Fig. 3D). These results indicate that EcR depletion results in An. funestus females that are less fertile than control females. Furthermore, no significant differences (t(46) = 0.5975, p = 0.5531, unpaired student’s t-test) existed between the dsGFP and uninjected control groups, corroborating that the nanoinjection procedure did not affect fertility in An. funestus (Additional file 4: Table S1). Taken together, results showed that EcR influences reproductive processes such as oogenesis and fertility in An. funestus.
EcR depletion decreases P. falciparum infection intensity in An. funestus.
The effect of EcR depletion on An. funestus infection by P. falciparum was determined by comparing the intensity of oocysts and the prevalence of infection between dsEcR-injected and dsGFP-injected mosquitoes. Oocyst intensity was significantly decreased by 50% in dsEcR treated An. funestus females with 3.80 ± 0.54 oocysts per midgut compared to dsGFP treated An. funestus females with 9.37 ± 1.17 oocysts per midgut (Mann–Whitney U = 1751, p = 0.0013, 8 biological replicates) (Fig. 4). EcR depletion significantly decreased the number of oocysts that developed in the midgut. Additionally, the mean prevalence of infected An. funestus females was 69.45% (± 6.05) in dsEcR injected An. funestus females compared to 80.46% (± 4.28) in dsGFP injected An. funestus females (Fig. 4). No significance in infection prevalence was observed between dsEcR and dsGFP treated females (χ2(1, N= 720) = 2.205, p = 0.1376, Chi-squared test), suggesting that EcR depletion does not influence the prevalence of infection by P. falciparum. These results corroborate that EcR influences P. falciparum oocyst development in An. funestus but does not affect the incidence of infection by P. falciparum. No significant difference was observed between the dsGFP and uninjected control groups for either the P. falciparum infection intensity (Mann–Whitney U = 2300, p = 0.5337) or prevalence (χ2(1, N= 720) = 0.2584, p = 0.6112, Chi-squared test), suggesting that the nanoinjection procedure did not influence these processes (Additional file 4: Table S1).
This study provides the first experimental evidence that EcR plays a crucial biological role in the major vector species An. funestus. Vectors such as An. funestus among others require novel control strategies due to the increase in insecticide resistance which reduces the efficacy of current vector control interventions. The transcription factor EcR was identified as a potential gene to be investigated due to its pleiotropic effects discovered in mosquitoes . The current study sought out to determine how 20E signaling, using EcR as a proxy, regulates several factors that influence vector density and competence in An. funestus.
Depletion of EcR in An. funestus reduced longevity by over a half compared to the control. Similarly, in other insect species such as Drosophila melanogaster, Nilaparvata lugens and Sitobion avenae, reduction of EcR using dsEcR resulted in decreased longevity [48,49,50]. Together these findings confirm that reducing transcription of EcR results in deleterious effects on longevity, corroborating that EcR is essential in the regulation of longevity. Longevity plays a major role in malaria transmission. Wild An. funestus can survive approximately 30 days  and in this time are able to complete the extrinsic incubation period (EIP) required for Plasmodium parasites to mature and become infectious to humans. Several factors such as temperature, vector species, parasite species, vector nutrition and survival can influence the EIP . Defining the duration of the EIP is therefore difficult, however it typically ranges from 12 to 16 days [51,52,53]. Reducing EcR in An. funestus females could therefore decrease longevity to a point where the EIP cannot be completed in the field, thereby reducing malaria transmission. Additionally, decreased vector longevity can reduce the vector to human ratio, possibly reducing malaria transmission rates .
Malaria transmission and vector density are also impacted by vector reproduction. Moreover, ingestion of a blood meal (a prerequisite for anautogenous mosquito reproduction) increases the risk of Plasmodium transmission. Two reproductive processes, fecundity and fertility were assessed after EcR depletion. Impaired oogenesis was observed after reduction of EcR levels in An. funestus. Consistent with these findings, other studies reported reduced follicle length  and a reduction in the number of eggs developed  after EcR depletion in Ae. aegypti and An. gambiae. Lp and VgR are two important proteins involved in lipid transport to oocytes and the uptake of the YPP vitellogenin providing nutrients for oocytes to develop respectively [35,36,37,38,39]. In the current study, depletion of EcR did not affect the expression of the Lp and VgR genes, however 20E signaling was confirmed to regulate Lp and Vg gene expression in other studies [35, 38, 54, 55]. However, these studies are not directly comparable due to differences in methodologies used and in addition discrepancy may be due to the complex interaction of several other genes in the 20E pathway suggesting that EcR is not solely responsible for the regulation of these genes.
The profound mechanism of gene regulation by EcR is not yet fully understood but studies have discovered several key elements of this process. Ecdysone responsive genes such as Vg i.e. those containing EcR binding sites, also contain binding sites for several other transcription factors which regulate gene expression differentially depending on temporal and spatial requirements . Changes in gene regulation in these ecdysone responsive genes are brought about indirectly by EcR coupled with the action of other transcription factors . However, EcR can also directly regulate the specific temporal or spatial expression of genes . Moreover, ecdysone responsive genes such as Lp are regulated differentially in mosquito species such as Ae. aegypti and An. gambiae through a currently unknown mechanism [35, 38]. The current study has limited information on gene regulation by EcR but the above studies confirm that future studies should concentrate on this process to provide a clearer understanding of this process in An. funestus.
Depletion of EcR significantly reduced fertility in An. funestus. A source of uncertainty arises from the small sample size in this study, however it is known that An. funestus is refractory to colonization due to low blood feeding rates, mating success and oviposition rates (Koekemoer, pers comm). This experiment should be repeated to provide more clarity. Nevertheless, results from the current study are noteworthy, a haem peroxidase (HPX15) was found to control fertility in An. gambiae . This enzyme requires 20E for normal functioning  and providing a mechanistic basis for the role of EcR in fertility. Silencing HPX15 was found to decrease fertility in females over multiple gonotrophic cycles . EcRs regulation of fertility could, therefore, be due to its function as a transcription factor of HPX15, but will have to be investigated in future.
Furthermore, EcR was found to control transcription of the immune molecule LRIM 9 in a study by Upton et al. . Depletion of LRIM 9 with dsLRIM 9 in An. gambiae resulted in a threefold higher P. berghei oocyst intensity, validating the importance of LRIM 9 in providing immunity against P. berghei . However, this was not evident when EcR was depleted in An. funestus as a decrease in P. falciparum oocyst intensity was observed instead. This discrepancy could presumably be due to a variance in the activity of LRIM 9 against different Plasmodium species studied. Similarly, it could be attributed to differences in the immune systems and regulation of Plasmodium infection between An. gambiae and An. funestus. Immune genes of An. funestus are genetically more similar to the Asian malaria vector Anopheles stephensi than An. gambiae due to evolutionary divergences in these vectors [59,60,61]. The contradiction in these results suggest a different mechanism is involved in decreasing the intensity of P. falciparum oocysts when EcR is knocked down in An. funestus. Similar to the current study, Werling et al.  observed a decrease in oocyst intensity when dsEcR treated An. gambiae was infected with P. falciparum. Strikingly, this decrease in P. falciparum oocysts was accompanied by larger, faster maturing sporozoites which was brought about by Lp . Presently, limitations such as a decrease in the longevity of dsEcR injected females (in this study) limited the sample size of infected mosquitoes available 16–18 days post infection and hindered the investigation of faster maturing sporozoites in An. funestus. These findings should however be further investigated in An. funestus once limitations are overcome to determine if the reduction in oocyst intensity results in faster maturing sporozoites and a more rapid EIP for P. falciparum or if a different immune strategy is adopted by the different Anopheles species. If the latter is proven true, an RNAi based approach targeting An. funestus EcR will need to be taken into consideration. It will also be interesting to determine the length of time that that EcR remains depleted in An. funestus after dsEcR treatment as this is a key factor in the development for a potential RNAi based control method.
If it is found to be true that depleting An. funestus EcR in any way increases parasite infectivity, genetically based control methods, such as gene drives , paratransgenesis , transmission blocking vaccines  or simply any form of EcR gene modification cannot be considered. Additionally, it will be important to confirm these findings on wild An. funestus populations from different geographical areas. Fortunately, decreasing EcR levels is only one of the possibilities to target ecdysone signalling to reduce malaria transmission. Nonsteroidal 20-hydroxyecdysone agonists share structural similarities with 20E, hence they are able to competitively bind to the EcR-USP complex, amplifying 20E signalling and all of its effects . Five 20E agonists, namely chromafenozide, fufenozide, halofenozide, methoxyfenozide and tebufenozide are currently available as insecticides targeting agricultural pests [66, 67]. These 20E agonists are promising as potential malaria vector control strategies as they exhibit low toxicity to non-target organisms and penetrate the mosquito cuticle [66,67,68]. The efficacy of some of these compounds has also been demonstrated against several Anopheles vectors [32, 33, 69].
The current study provided some insight on the biological function of EcR in An. funestus as a regulator of longevity, oogenesis, fertility and susceptibility to P. falciparum. Further research is required to determine the mode of action of, and genes involved in, EcRs regulation of pathways governing longevity, reproductive success and importantly Plasmodium infectivity in An. funestus. For example, the complex interaction of EcR and the 20E signalling pathway can further be elucidated through next generation RNA sequencing. This will allow for the identification of the complete profile of early and late genes regulated by EcR, providing a holistic overview of this pathway. These genes can thereafter be examined to determine their functions in regulating mosquito physiological process that potentially target vector density and competence. Subsequently, success in both laboratory and field based aspects of this research could result in approved EcR based vector control methods that will potentially reduce malaria transmission.
Availability of data and materials
The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.
Double stranded RNA
Green fluorescent protein
Indoor residual spraying
Long-lasting insecticidal nets
- LRIM 9:
Leucine rich immune molecule 9
World Health Organization
WHO. World malaria report 2020: 20 years of global progress and challenges. Geneva, World Health Organization, 2020. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2020. Accessed 5 Jan 2021
Gillies MT, De Meillon B. The Anophelinae of Africa south of the Sahara (Ethiopian zoogeographical region). Johannesburg: South African Institute Medical Research; 1968. p. 343.
Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa South of the Sahara. Publ S Afr Inst Med Res. 1987;55:1–143.
Cohuet A, Simard F, Toto JC, Kengne P, Coetzee M, Fontenille D. Species identification within the Anopheles funestus group of malaria vectors in Cameroon and evidence for a new species. Am J Trop Med Hyg. 2003;69:200–5.
Harbach RE. The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships. Bull Entomol Res. 2004;94:537–53.
Spillings BL, Brooke BD, Koekemoer LL, Chiphwanya J, Coetzee M, Hunt RH. A new species concealed by Anopheles funestus Giles, a major malaria vector in Africa. Am J Trop Med Hyg. 2009;81:510–5.
Wiebe A, Longbottom J, Gleave K, Shearer FM, Sinka ME, Massey NC, et al. Geographical distributions of African malaria vector sibling species and evidence for insecticide resistance. Malar J. 2017;16:85.
WHO. Global technical strategy for malaria 2016–2030. Geneva: World Health Organization; 2015. https://www.who.int/publications/i/item/9789241564991. Accessed 15 Apr 2020
WHO. Global technical strategy for malaria 2016–2030, 2021 update. Geneva: World Health Organization; 2021. https://www.who.int/publications/i/item/9789240031357. Accessed 28 Sept 2021
WHO. WHO Global malaria programme, Global plan for insecticide resistance management in malaria vectors (GPIRM). Geneva: World Health Organization; 2012. https://apps.who.int/iris/handle/10665/44846. Accessed 7 Apr 2020
Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.
Wondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, Morgan JC. Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa. Insect Biochem Mol Biol. 2011;41:484–91.
Okoye PN, Brooke BD, Koekemoer LL, Hunt RH, Coetzee M. Characterisation of DDT, pyrethroid and carbamate resistance in Anopheles funestus from Obuasi, Ghana. Trans R Soc Trop Med Hyg. 2008;102:591–8.
Amenya DA, Naguran R, Lo TCM, Ranson H, Spillings BL, Wood OR, et al. Over expression of a cytochrome P450 (CYP6P9) in a major African malaria vector, Anopheles funestus, resistant to pyrethroids. Insect Biochem Mol Biol. 2008;17:19–25.
Christian RN, Strode C, Ranson H, Coetzer N, Coetzee M, Koekemoer LL. Microarray analysis of a pyrethroid resistant African malaria vector, Anopheles funestus, from southern Africa. Pestic Biochem Physiol. 2011;99:140–7.
Riveron JM, Yunta C, Ibrahim SS, Djouaka R, Irving H, Menze BD, et al. A single mutation in the GSTe2 gene allows tracking of metabolically based insecticide resistance in a major malaria vector. Genome Biol. 2014;15:R27.
Coetzee M, Koekemoer LL. Molecular systematics and insecticide resistance in the major African malaria vector Anopheles funestus. Annu Rev Entomol. 2013;58:393–412.
Wood OR, Hanrahan S, Coetzee M, Koekemoer LL, Brooke BD. Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus. Parasit Vectors. 2010;3:67.
Moiroux N, Gomez MB, Pennetier C, Elanga E, Djènontin A, Chandre F, et al. Changes in Anopheles funestus biting behavior following universal coverage of long-lasting insecticidal nets in Benin. J Infect Dis. 2012;206:1622–9.
Sougoufara S, Diédhiou SM, Doucouré S, Diagne N, Sembène PM, Harry M, et al. Biting by Anopheles funestus in broad daylight after use of long-lasting insecticidal nets: a new challenge to malaria elimination. Malar J. 2014;13:125.
Neafsey DE, Christophides GK, Collins FH, Emrich SJ, Fontaine MC, Gelbart W, et al. The evolution of the Anopheles 16 genomes project. G3 (Besthesda). 2013;3:1191–4.
Ghurye J, Koren S, Small ST, Redmond S, Howell P, Phillippy AM, et al. A chromosome-scale assembly of the major African malaria vector Anopheles funestus. GigaScience. 2019;8:giz063.
Doudna JA, Charpentier E. The new frontier of genome engineering with CRISPR-Cas9. Science. 2014;346:1258096.
Koelle MR, Talbot WS, Segraves WA, Bender MT, Cherbas P, Hogness DS. The Drosophila EcR gene encodes an ecdysone receptor, a new member of the steroid receptor superfamily. Cell. 1991;67:59–77.
Oro AE, McKeown M, Evans RM. Relationship between the product of the Drosophila ultraspiracle locus and the vertebrate retinoid X receptor. Nature. 1990;347:298–301.
Yao T-P, Forman BM, Jiang Z, Cherbas L, Chen J-D, McKeown M, et al. Functional ecdysone receptor is the product of EcR and ultraspiracle genes. Nature. 1993;366:476–9.
Wang S-F, Li C, Zhu J, Miura K, Miksicek RJ, Raikhel AS. Differential expression and regulation by 20-hydroxyecdysone of mosquito ultraspiracle isoforms. Dev Biol. 2000;218:99–113.
Cherbas L, Lee K, Cherbas P. Identification of ecdysone response elements by analysis of the Drosophila Eip28/29 gene. Genes Dev. 1991;5:120–31.
Dobens L, Rudolph K, Berger EM. Ecdysterone regulatory elements function as both transcriptional activators and repressors. Mol Cellular Biol. 1991;11:1846–53.
Ekoka E, Maharaj S, Nardini L, Dahan-Moss Y, Koekemoer LL. 20-Hydroxyecdysone (20E) signaling as a promising target for the chemical control of malaria vectors. Parasit Vectors. 2021;14:86.
Shaw WR, Catteruccia F. Vector biology meets disease control: using basic research to fight vector-borne diseases. Nat Microbiol. 2018;4:20–34.
Childs LM, Cai FY, Kakani EG, Mitchell SN, Paton D, Gabrieli P, et al. Disrupting mosquito reproduction and parasite development for malaria control. PLoS Pathog. 2016;12: e1006060.
Brown F, Paton DG, Catteruccia F, Ranson H, Ingham VA. A steroid hormone agonist reduces female fitness in insecticide-resistant Anopheles populations. Insect Biochem Mol Biol. 2020;121: 103372.
Bryant B, Raikhel AS. Programmed autophagy in the fat body of Aedes aegypti is required to maintain egg maturation cycles. PLoS ONE. 2011;6: e25502.
Werling K, Shaw WR, Itoe MA, Westervelt KA, Marcenac P, Paton DG, et al. Steroid hormone function controls non-competitive Plasmodium development in Anopheles. Cell. 2019;177:315-325.e14.
Raikhel AS. Vitellogenesis in mosquitoes. Adv Dis Vector Res. 1992;9:1–39.
Sappington TW, Hays AR, Raikhel AS. Mosquito vitellogenin receptor: purification, developmental and biochemical characterization. Insect Biochem Mol Biol. 1995;25:807–17.
Sun J, Hiraoka T, Dittmer NT, Cho K-H, Raikhel AS. Lipophorin as a yolk protein precursor in the mosquito, Aedes aegypti. Insect Biochem Mol Biol. 2000;30:1161–71.
Rono MK, Whitten MMA, Oulad-Abdelghani M, Levashina EA, Marois E. The major yolk protein vitellogenin interferes with the anti-Plasmodium response in the malaria mosquito Anopheles gambiae. PLoS Biol. 2010;8: e1000434.
Upton LM, Povelones M, Christophides GK. Anopheles gambiae blood feeding initiates an anticipatory defense response to Plasmodium berghei. J Innate Immun. 2015;7:74–86.
Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M. Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector Anopheles funestus. Med Vet Entomol. 2005;19:271–5.
Reader J, Botha M, Theron A, Lauterbach SB, Rossouw C, Engelbrecht D, et al. Nowhere to hide: interrogating different metabolic parameters of Plasmodium falciparum gametocytes in a transmission blocking drug discovery pipeline towards malaria elimination. Malar J. 2015;14:213.
Cho W-L, Kapitskaya MZ, Raikhel AS. Mosquito ecdysteroid receptor: analysis of the cDNA and expression during vitellogenesis. Insect Biochem Mol Biol. 1995;25:19–27.
BEI. Methods in Anopheles research. 2016. https://www.beiresources.org/Portals/2/VectorResources/2016-Methods-in-Anopheles-Research-full-manual.pdf. Accessed 30 Nov 2020.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25:402–8.
Bradburn S. How to perform the delta-delta ct method in Excel. 2018. https://toptipbio.com/delta-delta-ct-pcr/. Accessed 31 Aug 2021.
Dhadialla TS, Hays AR, Raikhel AS. Characterization of the solubilized mosquito vitellogenin receptor. Insect Biochem Mol Biol. 1992;22:803–16.
Tricoire H, Battisti V, Trannoy S, Lasbleiz C, Pret A-M, Monnier V. The steroid hormone receptor EcR finely modulates Drosophila lifespan during adulthood in a sex-specific manner. Mech Ageing Dev. 2009;130:547–52.
Yan T, Chen H, Sun Y, Yu X, Xia L. RNA interference of the ecdysone receptor genes EcR and USP in grain aphid (Sitobion avenae F.) affects its survival and fecundity upon feeding on wheat plants. Int J Mol Sci. 2016;17:2098.
Yu R, Xu X, Liang Y, Tian H, Pan Z, Jin S, et al. The insect ecdysone receptor is a good potential target for RNAi-based pest control. Int J Biol Sci. 2014;10:1171.
Ohm JR, Baldini F, Barreaux P, Lefevre T, Lynch PA, Suh E, et al. Rethinking the extrinsic incubation period of malaria parasites. Parasit Vectors. 2018;11:178.
WHO. Manual on practical entomology in malaria. Geneva: World Health Organization; 1995. https://apps.who.int/iris/handle/10665/42481. Accessed 25 Nov 2021.
Vaughan JA, Noden BH, Beier JC. Population dynamics of Plasmodium falciparum sporogony in laboratory-infected Anopheles gambiae. J Parasitol. 1992;78:716–24.
Baldini F, Gabrieli P, South A, Valim C, Mancini F, Catteruccia F. The interaction between a sexually transferred steroid hormone and a female protein regulates oogenesis in the malaria mosquito Anopheles gambiae. PLoS Biol. 2013;11: e1001695.
Kokoza VA, Martin D, Mienaltowski MJ, Ahmed A, Morton CM, Raikhel AS. Transcriptional regulation of the mosquito vitellogenin gene via a blood meal-triggered cascade. Gene. 2001;274:47–65.
Martı́n D, Wang S-F, Raikhel AS. The vitellogenin gene of the mosquito Aedes aegypti is a direct target of ecdysteroid receptor. Mol Cell Endocrinol. 2001;173:75–86.
Uyehara CM, McKay DJ. Direct and widespread role for the nuclear receptor EcR in mediating the response to ecdysone in Drosophila. Proc Natl Acad Sci USA. 2019;116:9893–902.
Shaw WR, Teodori E, Mitchell SN, Baldini F, Gabrieli P, Rogers DW, et al. Mating activates the heme peroxidase HPX15 in the sperm storage organ to ensure fertility in Anopheles gambiae. Proc Natl Acad Sci USA. 2014;111:5854–9.
Kamali M, Marek PE, Peery A, Antonio-Nkondjio C, Ndo C, Tu Z, et al. Multigene phylogenetics reveals temporal diversification of major African malaria vectors. PLoS ONE. 2014;9: e93580.
Couto J, Ferrolho J, de la Fuente J, Domingos A. Anopheles gambiae and A. stephensi immune response during Plasmodium berghei infection uncovered by sialotranscriptomic analysis. Int Trends Immununol. 2016;4:58–72.
Mitri C, Bischoff E, Eiglmeier K, Holm I, Dieme C, Brito-Fravallo E, et al. Gene copy number and function of the APL1 immune factor changed during Anopheles evolution. Parasit Vectors. 2020;13:18.
Hammond AM, Galizi R. Gene drives to fight malaria: current state and future directions. Pathog Glob Health. 2017;111:412–23.
Wilke ABB, Marrelli MT. Paratransgenesis: a promising new strategy for mosquito vector control. Parasit Vectors. 2015;8:342.
Carter R. Transmission blocking malaria vaccines. Vaccine. 2001;19:2309–14.
Dhadialla TS, Carlson GR, Le DP. New insecticides with ecdysteroidal and juvenile hormone activity. Annu Rev Entomol. 1998;43:545–69.
Nakagawa Y. Nonsteroidal ecdysone agonists. Vitam Horm. 2005;73:131–73.
Smagghe G, Gomez LE, Dhadialla TS. Bisacylhydrazine insecticides for selective pest control. Adv Insect Physiol. 2012;43:163–249.
Schneider MI, Smagghe G, Pineda S, Vinuela E. Action of insect growth regulator insecticides and spinosad on life history parameters and absorption in third-instar larvae of the endoparasitoid Hyposoter didymator. Biol Control. 2004;31:189–98.
Reynolds RA, Smith RC. The 20-hydroxyecdysone agonist, halofenozide, primes anti-Plasmodium immunity in Anopheles gambiae via the ecdysone receptor. bioRxiv. 2020;3:1–9.
The authors would like to thank Prof Arbuthnot from the antiviral gene therapy unit at the University of Witwatersrand for providing the GFP plasmid.
Funding was provided by the Department of Science and Innovation (DSI) and the National Research Foundation (NRF) South African Research Chairs Initiative (UID 64763 to L.L.K. and UID 84627 to L.M.B.); the NRF Communities of Practice (UID 110666 to L.M.B. and L.L.K.) and the South African Medical Research Council Strategic Health Innovation Partnerships (SHIP) with funds from DSI.
Ethics approval and consent to participate
All in vitro experiments involving human blood donors and human malaria parasites hold ethics approval from the University of Pretoria Research Health Sciences Ethics Committee (506/2018) and Natural and Agricultural Sciences Ethics Committee (NAS180000094). The work abides by the Declaration of Helsinki principles.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Relative EcR expression levels in dsEcR injected An. funestus females compared to dsGFP injected An. funestus females. The EcR gene was knocked down in dsEcR injected An. funestus females as EcR expression levels were drastically reduced compared to the GFP control. Statistically significant knockdown was evident in dsEcR injected An. funestus females 24, 48 and 72 h after injection as EcR expression in dsEcR injected An. funestus females was 0.11 ± 0.006 (p < 0.05), 0.01 ± 0.001 (p < 0.01) and 0.2 ± 0.06 (p < 0.05) respectively when compared to the GFP injected control of 1. This data confirmed EcR knockdown in An. funestus females injected with dsEcR. Data is representative of 2 biological replicates and normalised using an average of RPS7 and RPL19 reference genes. Expression levels calculated using relative quantification method (∆∆Ct). At each time point statistical significance was assessed with the unpaired student’s t-test. *p < 0.05, **p < 0.01. Error bars represent standard deviation.
The highest mating success rate was achieved when An. funestus males and females are combined for 12 days after which no further increases are observed. The percentage mating success rate increased progressively until it reached its highest value of 62.2% after 12 days of mating. After this point, the mating success rate reached a plateau until day 20. Statistical significance was calculated using one-way ANOVA with Tukey’s post hoc analysis to correct for multiple comparison. Data represents the means of 3 biological replicates. Error bars represent standard deviation of means. ns = not statistically significant p > 0.05; ** = p < 0.01; *** = p < 0.001; **** = p < 0.0001. (n) = number of females per time point across 3 biological replicates.
Blood feeding rates did not differ between treatment groups. Insignificant differences amongst treatment groups confirmed that the blood feeding rates did not influence any changes observed in the phenotypes of dsEcR treated An. funestus females (p > 0.05). Statistical significance calculated using an unpaired student’s t-test. Error bars represent standard deviation. ns = not statistically significant p > 0.05.
Statistical significance between dsGFP and uninjected controls from the various biological assays conducted.
About this article
Cite this article
Maharaj, S., Ekoka, E., Erlank, E. et al. The ecdysone receptor regulates several key physiological factors in Anopheles funestus. Malar J 21, 97 (2022). https://doi.org/10.1186/s12936-022-04123-8
- African malaria vector
- Plasmodium falciparum
- RNA interference