- Open Access
The effect of silencing immunity related genes on longevity in a naturally occurring Anopheles arabiensis mosquito population from southwest Ethiopia
Parasites & Vectorsvolume 12, Article number: 174 (2019)
Vector control remains the most important tool to prevent malaria transmission. However, it is now severely constrained by the appearance of physiological and behavioral insecticide resistance. Therefore, the development of new vector control tools is warranted. Such tools could include immunization of blood hosts of vector mosquitoes with mosquito proteins involved in midgut homeostasis (anti-mosquito vaccines) or genetic engineering of mosquitoes that can drive population-wide knockout of genes producing such proteins to reduce mosquito lifespan and malaria transmission probability.
To achieve this, candidate genes related to midgut homeostasis regulation need to be assessed for their effect on mosquito survival. Here, different such candidate genes were silenced through dsRNA injection in the naturally occurring Anopheles arabiensis mosquitoes and the effect on mosquito survival was evaluated.
Significantly higher mortality rates were observed in the mosquitoes silenced for FN3D1 (AARA003032), FN3D3 (AARA007751) and GPRGr9 (AARA003963) genes as compared to the control group injected with dsRNA against a non-related bacterial gene (LacZ). This observed difference in mortality rate between the candidate genes and the control disappeared when gene-silenced mosquitoes were treated with antibiotic mixtures, suggesting that gut microbiota play a key role in the observed reduction of mosquito survival.
We demonstrated that interference with the expression of the FN3D1, FN3D3 or GPRGr9 genes causes a significant reduction of the longevity of An. arabiensis mosquito in the wild.
Sub-Saharan Africa hosts some of the most efficient malaria vectors, An. gambiae, An. arabiensis and An. funestus, and carries the heaviest malaria burden worldwide. A substantial reduction in malaria related cases and deaths have been recorded in the past decade . This progress is largely attributable to the scaling-up of vector control interventions, such as long-lasting insecticide-treated nets (LLINs) and indoor residual spraying (IRS), as well as improved diagnostics and effective treatment using artemisinin-based combination therapies .
LLINs and IRS impact malaria transmission largely by reducing the daily survival rate of mosquitoes that are mostly active at night and display strong endophagic (seeking blood meals indoors) and endophilic (rest indoors following a blood meal) behavior [3, 4]. These tools have been more efficient for the endophagic/endophilic An. gambiae and An. funestus mosquitoes compared to the exophagic/exophilic An. arabiensis mosquito [5,6,7,8,9]. In addition, field studies have reported evidence of behavioral adaptation of An. arabiensis to the LLINs and/or ITNs as mosquitoes feed earlier and more frequently outdoors, and rest more frequently outdoors [7, 10, 11].
Due at least in part to this difference in vector control efficacy, a marked shift in vector composition has been observed. In east African countries where these vectors coexist, An. arabiensis is gradually replacing An. gambiae and An. funestus, e.g. in Kenya [12,13,14] and Tanzania [7, 15]. Although An. arabiensis is known to be a less efficient vector compared to An. gambiae and An. funestus, the inherent resilience of the mosquito to LLINs and IRS has been linked to reports of resurgence or stagnation in rates of malaria cases and deaths in African countries . Therefore, there is a need to develop new complementary vector control technologies targeting vectors that are resilient to the current vector control tools.
Malaria transmission blocking vaccines and mosquito population replacement via genetic modification have recently become attractive technologies to complement the existing vector control interventions [16,17,18,19,20,21]. Another novel approach is to reduce the mosquito lifespan by introducing a lethal gene or a pathogen in the mosquito population [22, 23]. Along the same lines, the mosquito lifespan could be reduced by immunizing primary mosquito blood hosts (i.e. humans and domestic bovids) with mosquito proteins involved in the midgut homeostasis. For this approach, candidate genes need to be found, which is the topic of this paper. Different An. arabiensis genes related to midgut homeostasis were screened for their potential to reduce mosquito longevity. Our work was prompted by data showing that shortly after a blood meal the number of microbiota in the mosquito midgut increases drastically, up to 1000 times [24,25,26], which in normal circumstances triggers immune reactions that soon reduce the microbiota number to the basal level [27,28,29,30,31]. We hypothesized that by compromising the immune system the mosquito would no longer be able to control the microbiota in the gut, which would lead to a shorter lifespan. Previous studies have demonstrated in a laboratory colony of An. gambiae, that when genes encoding putative bacterial receptors such as PGRPLC, type III fibronectin domain proteins (FN3Ds including FN3D1, FN3D2, FN3D3) and the gustatory receptor GPRGr9 were silenced, gut homeostasis was disrupted [28, 32]. Here, we investigated the effect of silencing some of these genes by RNA interference on longevity in naturally occurring An. arabiensis populations.
Adult An. arabiensis mosquitoes were reared from larvae and pupae collected from natural breeding sites around Jimma, south-west Ethiopia (07°40′00″N, 36°50′00″E). Larvae and/or pupae were collected using a 350 ml mosquito dipper following the standard larvae collection procedure [33, 34]. The collected larvae and pupae were transported to an on-site mud-house where the pupae were separated from the larvae and transferred into a 10 ml beaker with water and kept in a mosquito cage with a dimension of 24.5 × 24.5 × 24.5 cm3 (Bugdorm-41515; Watkins & Doncaster, Leominster, UK) until emergence. The remaining larvae were transferred to a plastic tray containing water obtained from their natural habitat and fed on yeast and tropical fish food. The water was changed every 2 days and pupae were collected daily and transferred to the adult cage. Emerged adult mosquitoes were maintained on 10% sugar solution. Adult females of zero to two day old were transported to the experimental insectary at Jimma University for double stranded RNA (dsRNA) injection.
Total RNA was extracted from ten field-collected and laboratory reared An. arabiensis mosquitoes using TRIzol (Invitrogen, Inchinnan, UK) and cleaned with Turbo DNase I (Ambion, Huntington, UK). Complementary DNA (cDNA) was synthesized by reverse transcribing 1 μg of the total RNA using a Prime-Script™ 1st-strand cDNA Synthesis Kit (TaKaRa, Kusatsu, Shangai). Fragments of the five target genes were amplified by PCR using specific gene primers tailed with the short T7 promoter sequence TAA TAC GAC TCA CTA TAG G and the cDNA as a template. The targeted genes were FN3D1 (AARA003032), FN3D2 (AARA007751), FN3D3 (AARA007751), GPRGR9 (AARA003963) and PGRPLC3 (AARA002982). PCR fragment for the LacZ gene that served as a control was synthesized using a plasmid template containing the LacZ gene (for full primer sequences see Additional file 1: Table S1). DsRNA was synthesised from purified PCR products using the TranscriptAid T7 High Yield Transcription Kit (Thermo Fisher Scientific, Waltham, USA). The dsRNA was then purified using an RNeasy Mini Kit (Qiagen, Manchester, UK), following the manufacturer’s protocol. The concentration of dsRNA was determined spectrophotometrically by a Nanodrop 1000 spectrometer (Thermo Fisher Scientific, Wilmington, USA) at 260 nm and adjusted to 3 μg/μl using ultra-pure water. Gel electrophoresis (1% TBE agarose) was performed on a sub-sample of the PCR products to confirm that products of the expected size were detected for each gene. Zero- to two-days-old An. arabiensis mosquitoes were injected with 69 nl of dsRNA specific to a target gene or the LacZ control gene following a published RNA interference technique .
Gene silencing efficiency was measured for each of the 5 silenced genes using qrtPCR. Quantification of transcript abundance was performed on cDNA synthesized from total RNA extracted from mosquitoes injected with dsRNA 3 days earlier and maintained on 10% sugar solution. Fast SYBR® Green Master Mix Real-Time PCR Master Mix (Applied Biosystems, Warrington, UK) was used in the PCR reaction and amplification was detected by a 7500 Fast Real-Time PCR system (Applied Biosystems). Each target gene was quantified in duplicate. The AgS7 gene was used as an internal control. Primer sequences are given in Additional file 1: Table S1.
Monitoring of mosquito survival
The dsRNA injected mosquitoes were put into their respectively labelled cups. For each gene between 20 and 30 mosquitoes were injected per replicate. Cups with the mosquitoes were placed inside a 2 × 1 × 0.75 m3 (length × width × height) microclimate regulatory box constructed from chipboard. It has a window of 25 cm2 on each side and the top cover, and was covered with metal mesh to allow airflow (Additional file 2: Figure S1). The box was placed in a typical rural house. To prevent ant attacks, the box was placed on a 50 cm raised stand that was dipped halfway into water. It was lined with sawdust, 30 cm deep, which was sprinkled daily with water and kept closed. The box maintained a humidity of 60–70% RH and temperature of 25–28 °C throughout the whole study period, hence providing a microclimatic condition that resembled the natural resting habitat of An. arabiensis mosquitoes. The mosquitoes were offered a 10% sugar solution daily and a blood meal every 4 days by direct-feeding on a goat. Survival was monitored daily for 25 days, starting 24 h post-injection. Six independent replicates of mosquito injection were performed per gene.
Midgut microbiota analysis
A separate experiment was performed to determine microbiota load in both blood-fed and sugar-fed mosquitoes injected with the specified dsRNA. For each gene, between 20 and 30 mosquitoes were injected and transferred to their corresponding labelled cups. The mosquitoes were then kept in a microclimate regulatory box. On day 4 post-injection, the mosquitoes were either fed on blood or kept on a sugar meal. Twenty-four hours post-feeding, 5 blood-fed and 5 sugar-fed mosquitoes were sampled and their midguts were dissected. Individual midguts were homogenized in 100 µl of 4% paraformaldehyde (PFA) in phosphate-buffered saline (PBS). The number of bacteria was quantified using flow cytometry on the 5 pooled midgut samples per replicate . A total of 3 replicates of dsRNA injections were carried out per gene.
We tested the effect of reducing/eliminating the midgut microbiota load on the survival of the gene knocked down mosquitoes. For this purpose, dsRNA-injected mosquitoes were placed in 6 different cups, each cup containing 20–30 mosquitoes. The cups were kept in the microclimate regulatory box. On the day of dsRNA injection, the mosquitoes were given a cotton ball soaked in an antibiotic cocktail of streptomycin and norfloxacin, both at a dose rate of 10 μg/ml in a 10% sugar solution. On day 4 post-injection, the mosquitoes were blood-fed on a goat. The cotton balls were re-soaked with the antibiotic cocktail every fourth day for a period of 12 days (on days 4, 8 and 12). The antibiotics/sugar feed was carried out in the morning for about 1 h, and then cotton ball was removed from the cups for the next 5 h to starve the mosquitoes. The starved mosquitoes were offered a blood meal and the cotton ball was replaced on the cups after the feed. Mosquito survival was monitored daily for 25 days starting from 24 h post-injection. For this experiment 5 replicates of dsRNA injection were performed per gene.
Data analysis was performed using the statistical software package R v.3.3.2. For the gene silencing efficiency test, the relative expression of mRNA was calculated. The standard curve method was used for real-time qPCR quantification analysis. For each test sample, the PCR cycle number at which the fluorescent intensity of the reaction curve intersects the threshold line, i.e. level of detection or the point at which a reaction reaches a fluorescent intensity above background levels, known as threshold crossing values (Ct-values), was determined for the target and reference gene. Additionally, a standard curve was generated for both the target gene and the reference gene in each assay run using serial dilution of same template. Then, for each sample, the Ct-value of the target or reference gene was standardized using corresponding a standard curve and then the expression level of the target gene was normalized to the reference gene (AgS7 gene). The relative expression of mRNA of a target gene was compared between mosquitoes for which the target gene was silenced, and mosquitoes injected with control dsLacZ by a paired t-test using the replicate as a block factor. Gene silencing efficiency was expressed as the ratio of the relative expression of the target gene in the dsLacZ injected mosquitoes and the target gene silenced mosquitoes. Survival of the target gene and LacZds RNA-injected mosquitoes was depicted by Kaplan-Meier survival curves. The effect of silencing the different target genes on survival was modeled by the Cox proportional hazards frailty model, with replicate used as frailty term . The hazard ratio of a target gene over the control dsLacZ injection was used as summary statistic, together with the median time to death for the different gene silenced mosquitoes. To investigate the effect of the antibiotics cocktail, the same Cox proportional hazards model was fitted using the hazard ratio again as summary statistic. The bacterial counts were first log-transformed and then compared by a mixed model with replicate as random effect, and treatment, feed and the two-way interaction as fixed effects factor. The F-test was used to compare the silencing of the different target genes and the injection of control dsLacZ. The ratios of bacterial counts in the control dsLacZ injections and the target gene silencing were used as summary statistics.
Effect of gene silencing on survival
A significant reduction in expression level of mRNA compared to the dsLacZ-injected mosquitoes was achieved for the targeted candidate genes FN3D1 (t = -8.36, df = 2, P = 0.007), FN3D2 (t = -7.09, df = 2, P = 0.010), FN3D3 (t = -3.82, df = 2, P = 0.031) and GPRGr9 (t = -6.60, df = 2, P = 0.011) but not for PGRPLC3 (t = -1.42, df = 2, P = 0.145). The data are presented in Fig. 1.
Mosquito survival following silencing of the target genes compared to the control dsLacZ injected mosquitoes is depicted as a function of time for the each of the above genes in Fig. 2. Significantly higher mortality rates were observed for the FN3D1, FN3D3 and GPRGr9 knocked-down mosquitoes as compared to the control group, but not for the FN3D2 and GRPLC3 knocked-down mosquitoes (Table 1).
A marked reduction in average time to death was observed in mosquito groups where target genes were silenced compared to dsLacZ-injected controls, e.g. for GRPGr9 silenced mosquitoes the median time to death was almost halved from 20 (LacZ) to 11 days (Table 2).
Effect of gene silencing on the midgut bacterial count
The global analysis demonstrated that there was no overall significant difference between the target gene silenced and the dsLacZ-injected mosquitoes with respect to bacterial count (F(5,24) = 1.148, P = 0.363). However, there was an overall effect of the feed source (F(5,24) = 20.287, P < 0.001), i.e. blood-feeding versus sugar-feeding, whereas the interaction between the factors was not significant (F(5,24) = 1.902, P = 0.131). The ratio of bacteria count in blood-fed mosquitoes as compared to sugar-fed mosquitoes was 2.13 (95% CI: 1.49–3.06).
When considering the blood-fed mosquitoes alone, no statistically significant differences in the bacterial count of the target gene silenced and the dsLacZ-injected mosquitoes were observed. The ratio of midgut bacterial count of target gene silenced to the dsLacZ-injected mosquitoes is presented in Table 3; the ratios for all target gene silenced mosquitoes compared to the control groups (LacZ injected mosquitoes) were above 1, with the highest value for FN3D1 equal to 2.66 (95% CI: 0.94–7.57) but not significantly different from 1 (P = 0.085).
Treatment of mosquitoes with an antibiotics cocktail eliminated the gene silencing effect on survival for all 5 target genes. No statistically significant differences were noted in terms of mortality between the target genes silenced and the dsLacZ-injected mosquitoes (Table 1).
We also observed that for best performing genes, FN3D1 and FN3D3, dsRNA treatment had no effect on mosquito fecundity (number of eggs laid/female) or fertility (egg hatchability) of the mosquitoes (Additional file 3: Table S2).
The scale-up of vector control interventions in conjunction with early patient diagnosis and therapy has resulted in a substantial reduction in malaria-related cases and deaths since 2000. However, the newest data suggest that this progress is progressively coming to a halt or, in fact, being reversed in some countries, indicating that the current malaria intervention tools and strategies may have reached their maximum capacity . This highlights the urgency of developing new tools to complement LLINs to achieve the malaria elimination agenda in Africa, especially against the exophagic and exophilic An. arabiensis mosquito. A viable option against this opportunistic feeder species is to challenge different blood hosts with whole protein or antigenic peptide generated from mosquito itself (i.e. mass vaccination). The resulting host antibody is ingested by the mosquito during blood-feeding on an immunised host. This anti-mosquito antibody can be designed to target mosquito molecules involved in midgut homeostasis to reduce the longevity of the mosquito to such an extent that insufficient time remains to transmit the Plasmodium parasite.
In order for a mosquito to transmit a pathogen, it must acquire a blood meal from an infected person, support pathogen replication, dissemination to the salivary glands and take a subsequent blood meal from a susceptible host; this period is commonly called extrinsic incubation period (EIP). The EIP ranges from 10 to 14 days (2 to 6 gonotrophic cycles in areas of high malaria transmission) and a female mosquito must survive longer than the EIP to transmit the parasite , hence a very small fraction of mosquitoes (< 10%) can survive long enough to transmit the pathogen. This underlines the fact that mosquito longevity is a critical factor affecting malaria transmission and provides a basis for guiding vector control strategies [38, 39]. A small reduction in the mean mosquito survival period would have a significant impact on the disease transmission [40,41,42,43].
In the present study, we assessed five midgut-expressed genes FN3D1, FN3D2, FN3D3, GPRGr9 and PGRPLC for their potential as lifespan-limiting targets in a wild An. arabiensis population. Previous studies with An. gambiae under laboratory settings have demonstrated that the above genes regulate midgut microbiota [28, 32]. Unlike these studies, our work focused on the effect of those genes on the survival of field-caught An. arabiensis. Our experimental mosquitoes were captured as larvae or pupae and maintained in an environment simulating the mosquito natural resting habitat. They were also allowed to blood-feed on a goat every fourth day, mimicking their natural habits. We have demonstrated that in these semi-natural conditions, the reduction of FN3D1, FN3D3 or GPRGr9 expression significantly reduces the longevity of the An. arabiensis mosquitoes to an average of 12, 13 and 11 days, respectively, compared to a 20-day average longevity of control mosquitoes.
The observed effect on the mosquito survival is probably linked to the disruption of the mosquito midgut homeostasis as the observed effect of gene silencing (particularly of FN3D1, FN3D3 and GPRGr9) on their survival was no longer seen when mosquitoes were treated with an antibiotic cocktail to eliminate their gut microbiota. Our hypothesis that the reduced survival is due to the inability of mosquitoes to control their gut microbiota is also supported by other studies [44, 45]. Furthermore, it has been previously demonstrated that FN3Ds and GPRGr9 have a specific effect on microbiota of the family Enterobacteriaceae .
Our results demonstrate that interfering with the expression and/or function of these genes reduces the mosquito lifespan which in turn will significantly reduce malaria transmission. A permanent gene inactivation/knock-out can be achieved in various ways. A new and most promising tool is the recently developed CRISPR/Cas9-based genome editing methodology and gene-drive systems for Anopheles mosquitoes [21, 46, 47]. Recently a CRISPR/Cas9-induced somatic gene disruption technique has been established in An. gambiae [47, 48]. Such a knock-out line can be crossed with a germline-Cas9 strain as described in to generate a germ-line gene-knockout line which might be released to introgress the life-shortening trait into the wild malaria vector population . A second approach could be through the immunization of the blood-providing host, whether human or domestic animals, with molecules derived from target vector proteins that play a role in midgut homeostasis. The purified form of such molecules (i.e. antigens) can be inoculated into the vertebrate host to induce host immune reaction, and ultimately produce specific antibodies. Mosquitoes feeding on the immunized hosts would then ingest antibodies that neutralize the function of the protein leading to disrupted midgut homeostasis and shortened lifespan. It has been previously demonstrated that An. gambiae-derived anti-midgut monoclonal antibodies significantly reduce vector survivorship . This approach is particularly attractive for zoophilic mosquitoes such as An. arabiensis. Mosquitoes that have taken an infectious blood meal will typically take three to four additional blood meals before the completion of a sporogonic period as they normally blood-feed every 2–3 days. This could ensure repeated ingestion of anti-mosquito antibodies with consequential disruption of gut bacterial homeostasis to ultimately induce reduced lifespan. This technology can also impact secondary malaria vectors including Anopheles rivolorum, Anopheles pharoensis, Anopheles coustani, Anopheles ziemanni and Anopheles squamosus that are responsible for about 5% of total malaria transmission in Africa and are often highly zoophilic .
As the anti-mosquito vaccine technologies are expected to only reduce the long-term survival of vectors of the mosquitoes, female mosquitoes will complete some of their gonotrophic cycle. Thus, the technologies are expected to have minimum selection pressure inducing resistance .
This approach has been successfully used for an anti-tick vaccine targeting the midgut antigens Bm86 and an anti-tick and anti-mosquito vaccine targeting the subolesin/akirin (SUB/AKR) antigens [53,54,55,56]. The functional model for the SUB/AKR vaccine involves the nuclear factor-kappa B (NF-kB) of vector insects to inhibit the Immunodeficiency (Imd) pathway, which is important for regulation of the gut microbiota [57, 58].
A major challenge with such protein antigens is that they can result in autoimmunity due to molecular mimicry as such molecules can possess ‘mimotopes’ that are peptides mimicking the antigenic conformation structures that are recognized by the paratope antibody, leading to autoimmunity .
In the present case, amino acid sequence analysis of the An. arabiensis FN3D1 and FN3D3 genes has shown the presence of 25–37% sequence identity with proteins of some human genes such as the protein tyrosine phosphatase receptor type F gene. This level of homology is often referred as the ‘twilight zone’ where the structural similarity between the target mosquito genes and the human genes cannot be ruled out, suggesting a potential risk of induction of autoimmune diseases in individuals upon immunization with the mosquito genes . For instance, an analysis involving over a million sequences with known structures showed that at the top cut-off of the twilight zone, about 90% of protein pairs were structurally homologous, but the level of homology reduced drastically (to < 10%) when the sequence similarity between protein pairs was below 25% . On the other hand, both mosquito genes have no significant sequence similarity in the domestic Bovidae, hence induction of an autoimmune reaction in the vaccinated animals is negligible and the protein molecules can be used to develop vaccines to immunize the animals with minimum risk.
Eliminating the expression of the midgut proteins FN3D1, FN3D3 or GPRGr9 significantly reduces the lifespan of naturally occurring An. arabiensis mosquitoes reared in field conditions. The effect is probably caused by disruption of the mosquito midgut homeostasis through interference with the midgut microbiota, eventually hampering the mosquito immuno-metabolic functions. Therefore, these proteins can be good targets of mosquito life-shortening interventions, such as anti-mosquito vaccines or mosquito genetic modification, resulting in mosquitoes that can survive long enough to complete a gonotrophic cycle but not long enough to transmit malaria parasites to a new host. Thus, the technologies are expected to have minimum selection pressure inducing resistance.
threshold crossing values
double stranded RNA
extrinsic incubation period
- FN3D1 :
fibronectin type III domain-protein 1
- FN3D2 :
fibronectin type III domain-protein 2
- FN3D3 :
fibronectin type III domain-protein 3
- GPRGR9 :
G protein coupled receptor protein 9
indoor residual spraying
long-lasting insecticide-treated nets
- PGRPLC3 :
peptidoglycan recognition protein LC3
polymerase chain reaction
quantitative real-time PCR
transmission blocking vaccine
World Health Organization. World malaria report 2018. http://www.who.int/iris/handle/10665/275867. Accessed 16 Jan 2019.
MacDonald G. The epidemiology and control of malaria. London: Oxford University Press; 1957.
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, et al. Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010;9:62.
Gratz NR. The future of vector biology and control in the World Health Organization. J Am Mosq Control Assoc. 1985;1:273–8.
Kitau J, Oxborough RM, Tungu PK, Matowo J, Malima RC, Magesa SM, et al. Species shifts in the Anopheles gambiae complex: do LLINs successfully control Anopheles arabiensis? PLoS One. 2012;7:e31481.
Eckhoff PA. A malaria transmission-directed model of mosquito life cycle and ecology. Malar J. 2011;10:303.
Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF. Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011;10:80.
Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, et al. Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010;7:e1000303.
Soko W, Chimbari MJ, Mukaratirwa S. Insecticide resistance in malaria-transmitting mosquitoes in Zimbabwe: a review. Infect Dis Poverty. 2015;4:46.
Mwangangi JM, Mbogo CM, Orindi BO, Muturi EJ, Midega JT, Nzovu J, et al. Shifts in malaria vector species composition and transmission dynamics along the Kenyan coast over the past 20 years. Malar J. 2013;12:13.
Yohannes M, Boelee E. Early biting rhythm in the afro-tropical vector of malaria, Anopheles arabiensis, and challenges for its control in Ethiopia. Med Vet Entomol. 2012;26:103–5.
Derua YA, Alifrangis M, Magesa SM, Kisinza WN, Simonsen PE. Sibling species of the Anopheles funestus group, and their infection with malaria and lymphatic filarial parasites, in archived and newly collected specimens from north-eastern Tanzania. Malar J. 2015;14:104.
Mutuku FM, King CH, Mungai P, Mbogo C, Mwangangi J, Muchiri EM, et al. Impact of insecticide-treated bed nets on malaria transmission indices on the south coast of Kenya. Malar J. 2011;10:356.
Zhou G, Afrane YA, Vardo-Zalik AM, Atieli H, Zhong D, Wamae P, et al. Changing patterns of malaria epidemiology between 2002 and 2010 in western Kenya: the fall and rise of malaria. PLoS One. 2011;6:e20318.
Trape JF, Tall A, Diagne N, Ndiath O, Ly AB, Faye J, et al. Malaria morbidity and pyrethroid resistance after the introduction of insecticide-treated bed nets and artemisinin-based combination therapies: a longitudinal study. Lancet Infect Dis. 2011;11:925–32.
Eckhoff PA, Wenger EA, Godfray HCJ, Burt A. Impact of mosquito gene drive on malaria elimination in a computational model with explicit spatial and temporal dynamics. Proc Natl Acad Sci USA. 2017;114:e255.
Alphey L, Beard CB, Billingsley P, Coetzee M, Crisanti A, Curtis C, et al. Malaria control with genetically manipulated insect vectors. Science. 2002;298:119–21.
Deredec A, Godfray HCJ, Burt A. Requirements for effective malaria control with homing endonuclease genes. Proc Natl Acad Sci USA. 2011;108:e874.
Burt A. Site-specific selfish genes as tools for the control and genetic engineering of natural populations. Proc Biol Sci. 2003;270:921–8.
James AA. Gene drive systems in mosquitoes: rules of the road. Trends Parasitol. 2005;21:64–7.
Gantz VM, Jasinskiene N, Tatarenkova O, Fazekas A, Macias VM, Bier E, et al. Highly efficient Cas9-mediated gene drive for population modification of the malaria vector mosquito Anopheles stephensi. Proc Natl Acad Sci USA. 2015;112:e6736.
Wilke ABB, Marrelli MT. Paratransgenesis: a promising new strategy for mosquito vector control. Parasit Vectors. 2015;8:342.
Ricci I, Valzano M, Ulissi U, Epis S, Cappelli A, Favia G. Symbiotic control of mosquito borne disease. Pathog Glob Health. 2012;106:380–5.
Pumpuni CB, Demaio J, Kent M, Davis JR, Beier JC. Bacterial population dynamics in three anopheline species: the impact on Plasmodium sporogonic development. Am J Trop Med Hyg. 1996;54:214–8.
Kumar S, Molina-Cruz A, Gupta L, Rodrigues J, Barillas-Mury C. A peroxidase/dual oxidase system modulates midgut epithelial immunity in Anopheles gambiae. Science. 2010;327:1644–8.
de O Gaio A, Gusmão DS, Santos AV, Berbert-Molina MA, Pimenta PF, Lemos FJ. Contribution of midgut bacteria to blood digestion and egg production in Aedes aegypti (Diptera: Culicidae) (L.). Parasit Vectors. 2011;4:105.
Gendrin M, Christophides GK. The Anopheles mosquito microbiota and their impact on pathogen transmission. In: Manguin S, editor. Anopheles mosquitoes - new insights into malaria vectors. Rijeka: InTech; 2013.
Meister S, Agianian B, Turlure F, Relogio A, Morlais I, Kafatos FC, et al. Anopheles gambiae PGRPLC-mediated defense against bacteria modulates infections with malaria parasites. PLoS Pathog. 2009;5:e1000542.
Cirimotich CM, Dong Y, Clayton AM, Sandiford SL, Souza-Neto JA, Mulenga M, et al. Natural microbe-mediated refractoriness to Plasmodium infection in Anopheles gambiae. Science. 2011;332:855–8.
Rodrigues J, Brayner FA, Alves LC, Dixit R, Barillas-Mury C. Hemocyte differentiation mediates innate immune memory in Anopheles gambiae mosquitoes. Science. 2010;329:1353–5.
Dong Y, Manfredini F, Dimopoulos G. Implication of the mosquito midgut microbiota in the defense against malaria parasites. PLoS Pathog. 2009;5:e1000423.
Stathopoulos S, Neafsey DE, Lawniczak MKN, Muskavitch MAT, Christophides GK. Genetic dissection of Anopheles gambiae gut epithelial responses to Serratia marcescens. PLoS Pathog. 2014;10:e1003897.
World Health Organization. Manual of practical entomology in malaria. Part II. Methods and techniques. Geneva: WHO Offset Publication; 1975.
Gerberg EJ. Manual for mosquito rearing and experimental techniques. Bulletin No. 5. Selma, CA: American Mosquito Control Association, Inc.; 1970.
Habtewold T, Duchateau L, Christophides GK. Flow cytometry analysis of the microbiota associated with the midguts of vector mosquitoes. Parasit Vectors. 2016;9:167.
Duchateau L, Janssen P. The frailty model. New York: Springer Verlag; 2008.
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.
Bruce-Chwatt LJ. Essential malariology. London: William Heinemann Medical Books Ltd; 1980.
Bara J, Rapti Z, Cáceres CE, Muturi EJ. Effect of larval competition on extrinsic incubation period and vectorial capacity of Aedes albopictus for dengue virus. PLoS One. 2015;10:e0126703.
Charlwood JD, Smith T, Billingsley PF, Takken W, Lyimo EOK, Meuwissen JHET. Survival and infection probabilities of anthropophagic anophelines from an area of high prevalence of Plasmodium falciparum in humans. Bull Entomol Res. 1997;87:445–53.
Scholte EJ, Ng’habi K, Kihonda J, Takken W, Paaijmans K, Abdulla S, et al. An entomopathogenic fungus for control of adult African malaria mosquitoes. Science. 2005;308:1641–2.
Dye C. Vectorial capacity: must we measure all its components? Parasitol Today. 1986;2:203–9.
Garrett-Jones C, Shidrawi GR. Malaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomology. Bull World Health Organ. 1969;40:531–45.
Asgari S. Host-microbe interactions: a case for Wolbachia dialogue. In: Stephen KW, Serap A, Dimopoulos G, editors. Arthropod vector: controller of disease transmission, volume 1: Vector microbiome and innate immunity of arthropods. London: Academic Press; 2017. p. 173–83.
Gendrin M, Rodgers FH, Yerbanga RS, Ouedraogo JB, Basanez MG, Cohuet A, Christophides GK. Antibiotics in ingested human blood affect the mosquito microbiota and capacity to transmit malaria. Nat Commun. 2015;6:5921.
Zhang H, McCarty N. CRISPR editing in biological and biochemical investigation. J Cell Biochem. 2017;118:4152–62.
Dong Y, Simões ML, Marois E, Dimopoulos G. CRISPR/Cas9-mediated gene knockout of Anopheles gambiae FREP1 suppresses malaria parasite infection. PLoS Pathog. 2018;14:e1006898.
Hammond A, Galizi R, Kyrou K, Simoni A, Siniscalchi C, Katsanos D, et al. A CRISPR-Cas9 gene drive system targeting female reproduction in the malaria mosquito vector Anopheles gambiae. Nat Biotechnol. 2016;34:78–83.
Macias VM, Ohm JR, Rasgon JL. Gene drive for mosquito control: where did it come from and where are we headed? Int J Environ Res Public Health. 2017;14:1006.
Lal AA, Patterson PS, Sacci JB, Vaughan JA, Paul C, Collins WE, Wirtz RA, Azad AF. Anti-mosquito midgut antibodies block development of Plasmodium falciparum and Plasmodium vivax in multiple species of Anopheles mosquitoes and reduce vector fecundity and survivorship. Proc Natl Acad Sci USA. 2001;98:5228–33.
Afrane YA, Bonizzoni M, Yan G. Secondary malaria vectors of sub-Saharan Africa: threat to malaria elimination on the continent? In: Rodriguez-Morale Alfonso J, editor. Current topics in malaria. Rijeka: InTech; 2016. https://doi.org/10.5772/65359
Read AF, Lynch PA, Thomas MB. How to make evolution-proof insecticides for malaria control. PLoS Biol. 2009;7:e1000058.
Jonsson NN, Matschoss AL, Pepper P, Green PE, Albrecht MS, Hungerford J, et al. Evaluation of tick GARD(PLUS), a novel vaccine against Boophilu smicroplus, in lactating Holstein-Friesian cows. Vet Parasitol. 2000;88:275–85.
da Costa M, Pinheiro-Silva R, Antunes S, Moreno-Cid JA, Custódio A, Villar M, et al. Mosquito Akirin as a potential antigen for malaria control. Malar J. 2014;13:470.
de la Fuente J, Moreno-Cid JA, Canales M, Villar M, de la Lastra JMP, Kocan KM, et al. Targeting arthropod subolesin/akirin for the development of a universal vaccine for control of vector infestations and pathogen transmission. Vet Parasitol. 2011;181:17–22.
de la Fuente J, Moreno-Cid JA, Galindo RC, Almazan C, Kocan KM, Merino O, et al. Subolesin/Akirin vaccines for the control of arthropod vectors and vector borne pathogens. Transbound Emerg Dis. 2013;60(Suppl. 2):172–8.
Naranjo V, Ayllón N, Pérez de la Lastra JM, Galindo RC, Kocan KM, Blouin EF, et al. Reciprocal regulation of NF-kB (Relish) and Subolesinin the tick vector, Ixodes scapularis. PLoS One. 2013;8:e65915.
Wan F, Lenardo MJ. Specification of DNA binding activity of NF-kappaB proteins. Cold Spring Harb Perspect Biol. 2009;1:a000067.
Huang J, Ru B, Zhu P, Nie F, Yang J, et al. MimoDB 2.0: a mimotope database and beyond. Nucleic Acids Res. 2011;40:D271–7.
Rost B. Twilight zone of protein sequence alignments. Protein Eng. 1999;12:85–94.
We thank Professor Delenasaw Yewhalaw for kindly allowing to use laboratory space and materials at the Tropical and Infectious Diseases Research Centre, Jimma University, Ethiopia, and Mr Kassahun Ebba for logistic support.
This study was funded by the Institutional University Cooperation IUC-JU project in the framework of the Flemish Interuniversity Council (VLIR-UOS). The project was additionally supported by the Welcome Trust Investigator Award to GKC, and at later stages by the Bill and Melinda Gates Foundation grant OPP1158151 to GKC.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The datasets are available from the corresponding author upon request.
TH, LD and GKC developed the research idea. SD and TH performed the mosquito field studies. TH and SD performed laboratory work. LD, TH and SD analysed and interpreted the data. All authors were involved in writing. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Ethical clearance was obtained from the research and ethical review board of the College of Agriculture and Veterinary Medicine, Jimma University.
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.