Mosquito immune responses and compatibility between Plasmodium parasites and anopheline mosquitoes
© Jaramillo-Gutierrez et al. 2009
Received: 18 December 2008
Accepted: 30 July 2009
Published: 30 July 2009
Skip to main content
© Jaramillo-Gutierrez et al. 2009
Received: 18 December 2008
Accepted: 30 July 2009
Published: 30 July 2009
Functional screens based on dsRNA-mediated gene silencing identified several Anopheles gambiae genes that limit Plasmodium berghei infection. However, some of the genes identified in these screens have no effect on the human malaria parasite Plasmodium falciparum; raising the question of whether different mosquito effector genes mediate anti-parasitic responses to different Plasmodium species.
Four new An. gambiae (G3) genes were identified that, when silenced, have a different effect on P. berghei (Anka 2.34) and P. falciparum (3D7) infections. Orthologs of these genes, as well as LRIM1 and CTL4, were also silenced in An. stephensi (Nijmegen Sda500) females infected with P. yoelii (17XNL). For five of the six genes tested, silencing had the same effect on infection in the P. falciparum-An. gambiae and P. yoelii-An. stephensi parasite-vector combinations. Although silencing LRIM1 or CTL4 has no effect in An. stephensi females infected with P. yoelii, when An. gambiae is infected with the same parasite, silencing these genes has a dramatic effect. In An. gambiae (G3), TEP1, LRIM1 or LRIM2 silencing reverts lysis and melanization of P. yoelii, while CTL4 silencing enhances melanization.
There is a broad spectrum of compatibility, the extent to which the mosquito immune system limits infection, between different Plasmodium strains and particular mosquito strains that is mediated by TEP1/LRIM1 activation. The interactions between highly compatible animal models of malaria, such as P. yoelii (17XNL)-An. stephensi (Nijmegen Sda500), is more similar to that of P. falciparum (3D7)-An. gambiae (G3).
Mosquitoes transmit many infectious diseases, including malaria, lymphatic filariasis, yellow fever, and dengue. Among these diseases, malaria is by far the most costly in terms of human health. It is endemic to more than 100 countries and causes 550 million cases per year, with the highest mortality in children from sub-Saharan Africa. Malaria transmission to humans requires a competent mosquito species, as Plasmodium parasites must undergo a complex developmental cycle and survive the defense responses of their insect host. In Africa, Anopheles gambiae is the major vector of Plasmodium falciparum infection, which causes the most aggressive form of human malaria.
The Plasmodium berghei (murine malaria) model is one of the most widely used experimental systems to study malaria transmission. Gene silencing by systemic injection of double-stranded RNA (dsRNA) has proven to be a very useful tool to carry out functional genomic screens aimed at identifying mosquito genes that mediate anti-parasitic responses. In general, Anopheles gambiae is considered to be susceptible to P. berghei infection, because a high prevalence of infection can be achieved and parasites are only rarely melanized; however, silencing of either thioester-containing protein 1 (TEP1) , leucine-rich repeat immune protein 1 (LRIM1) , or LRIM2 (also called APL1, ), enhances P. berghei infection by 4–5 fold; indicating that, when these effector molecules are present, about 80% of parasites are eliminated by a lytic mechanism. It is well documented that An. gambiae mosquitoes have a different transcriptional response to infection with P. berghei and P. falciparum [4, 5] and genes such as LRIM1 and C-type lectin 4 (CTL4) , which limit or enhance P. berghei infection, respectively, do not affect P. falciparum infection in An. gambiae . This raises the possibility that some antiplasmodial genes identified using the P. berghei malaria model may not be relevant to human malaria transmission.
More than 400 species of anopheline mosquitoes have been identified, but only 40 of them are considered to be important disease vectors . Different anopheline species and even particular strains of mosquitoes vary widely in their susceptibility to infection with a given Plasmodium parasite species. For example, twelve different strains of Anopheles stephensi have been shown to have very different susceptibility to P. falciparum (Welch strain) infection . Furthermore, susceptibility had a strong genetic component, which allowed selection of a An. stephensi strain (Nijmegen Sda500) that is highly susceptible to P. falciparum infection . A strain of An. gambiae (L35) was selected to be highly refractory to infection with Plasmodium cynomolgy (primate malaria). The L35 strain melanizes P. cynomolgy, as well as several other Plasmodium species such as P. berghei (murine malaria), Plasmodium gallinaceum (avian malaria), and other primate malaria parasites such as Plasmodium gonderi, Plasmodium inui, and Plasmodium knowlesi. Interestingly, P. falciparum strains from the New World are also melanized effectively, but not those of African origin, suggesting that there are genetic differences between P. falciparum strains that affect their ability to infect An. gambiae . The African strains of P. falciparum tested appeared to be better adapted to their natural mosquito vector. However, great differences in the level of resistance to P. falciparum infection have been documented in families derived from individual An. gambiae females collected in the field [3, 10], and a small region of chromosome 2L is a major determinant of genetic resistance to infection .
Drosophila melanogaster can support the development of Plasmodium gallinaceum oocysts when cultured ookinetes are injected into the hemocele . This observation opened the possibility of using a genetic approach to screen for Drosophila genes that affect Plasmodium P. gallinaceum infection. Furthermore, silencing of orthologs (or family members) of five of these candidate genes in An. gambiae (G3 strain) demonstrated that four of them also affected P. berghei infection in the mosquito .
In this study we compare how silencing a set of genes identified in the Drosophila screen affects Plasmodium infection in different vector-parasite combinations. We conclude that there is a broad range of compatibility between different Plasmodium strains and particular mosquito strains that is determined by the interaction between the parasite and the mosquito's immune system. We define compatibility as the extent to which the immune system of the mosquito is actively limiting Plasmodium infection. For example, the P. yoelii-An. stephensi and P. falciparum-An. gambiae strains used in this study are highly compatible vector-parasite combinations, as silencing several genes involved in oxidative response or immunity has no significant effect on infection. In contrast, silencing the same genes has a strong effect in less compatible vector-parasite combinations such as P. yoelii-An. gambiae or P. berghei-An. gambiae.
Effect of silencing seven An. gambiae genes or their orthologs in An. stephensi on the intensity of P. berghei, P. falciparum or P. yoelii infection.
An. gambiae Gene ID
An. gambiae P. berghei (21°C)
An. gambiae P. falciparum (26°C)
An. stephensi P. yoelii (24°C)
Silencing ArgK, Sol. Trsp., and tetraspanin genes has a similar effect on P. berghei and P. falciparum infection. ArgK is a key enzyme in cellular energy homeostasis in arthropods, with a function similar to that of creatine kinase in mammals. This enzyme catalyzes the synthesis of phosphoarginine, which serves as an energy reserve. The high-energy phosphate in phosphoarginine can be transferred to ADP to renew ATP during periods of high energy demand . Apparently, silencing this enzyme results in a physiologic state in the mosquito that does not foster the development of either P. berghei or P. falciparum. Silencing of the solute transporter has no effect, while knockdown of tetraspanin enhances infection with both parasites. Tetraspanins are proteins with four transmembrane (TM) domains that are associated extensively with one another and with other membrane proteins to form specific microdomains distinct from lipid rafts. They are expressed on the surface of numerous cell types and are involved in diverse processes from cell adhesion to signal transduction and some of them inhibit the function of other members of the same family of proteins . CD81 is a tetraspanin that has been shown to be required for hepatocyte invasion by P. falciparum and P. yoelii sporozoites . Silencing of the An. gambiae tetraspanin gene may enhance parasite invasion and/or prevent the activation of an immune cascade that limits infection with P. berghei and P. falciparum.
OXR1, GSTT1, GSTT2 and Hsc-3 silencing has a different effect on P. berghei and P. falciparum infection. In yeast and mammals, OXR1 is induced by heat and oxidative stress and prevents oxidative damage by an unknown mechanism . In An. gambiae, OXR1 silencing decreases resistance to oxidative challenge and prevents the induction of genes involved in ROS detoxification, such as catalase, following a blood meal (G. Jaramillo-Gutierrez and C. Barillas-Mury, unpublished). We have previously shown that higher ROS levels in An. gambiae reduce P. berghei infection . Thus, it is likely that the decrease in P. berghei infectivity following OXR1 silencing is due to an increase in ROS. The unexpected observation that OXR1 silencing does not affect P. falciparum infection suggests that either this parasite species is less susceptible to oxidative stress or that the ingestion of human blood results in less accumulation of ROS in the mosquito.
GSTs play an important role as antioxidants and are involved in the detoxification of xenobiotics. GSTs of the epsilon and delta class have been extensively studied for their role in insecticide resistance in mosquitoes . The GST-Theta1 (GSTT1) null genotype in human males is highly associated to increased risk of basal cell carcinoma of the skin . Furthermore, in diabetics, the deletion of one copy of the GSTT1 gene is associated with elevated markers of inflammation and lipid peroxidation . Therefore, silencing of GSTT1 and GSTT2 could result in increased lipid peroxidation, which is expected to be deleterious to P. berghei; however, it is not clear why reducing GSTT2 expression enhances P. falciparum infection.
The observed differences in the effect of silencing specific An. gambiae (G3 strain) genes on P. berghei and P. falciparum infection may reflect the degree of compatibility between these two parasite species and the mosquito strain used. Alternatively, mosquitoes may trigger different sets of effector genes in response to different Plasmodium species. To explore these possibilities, we evaluated the responses of two mosquito species that differ in their susceptibility to the same Plasmodium parasite.
An. gambiae (G3) and An. stephensi (Nijmegen Sda500) infections with P. yoelii.
Prevalence of infection
Median live oocyst number
% of midguts with melanized parasites
% of midguts with live and melanized parasites
n = 59
n = 47
The effect of silencing multiple mosquito genes in the highly compatible P. yoelii (17XNL)-An. stephensi (Nijmegen Sda500)system was very similar to that observed when P. falciparum (3D7) was used to infect An. gambiae (G3), its natural vector; suggesting that P. yoelii-An. stephensi is a representative animal model to study P. falciparum interactions with compatible vectors. Furthermore, P. yoelii-infected females can be kept at 24°C, a temperature that is more physiological for mosquitoes and closer to that used for P. falciparum infections (26°C).
Using less compatible parasite-mosquito combinations, such as the P. berghei-An. gambiae or P. yoelii-An. gambiae strains described in this study, may be particularly useful to identify and characterize immune pathways in the mosquito that could potentially limit human malaria transmission. Once a potential pathway is defined, it is possible to investigate if certain parasite strains avoid activating them, or if the effector genes are inefficient. It may also be possible to use alternative strategies (such as chemicals or fungal infections) to activate these potential antiplasmodial responses and test their effectiveness in limiting malaria transmission in natural vector-parasite combinations.
There is a broad spectrum of compatibility between different strains of Plasmodium and particular mosquito strains; for example, An. gambiae (G3) is highly compatible with P. falciparum (3D7) parasites, but has low compatibility with P. yoelii 17XNL. A given strain of Plasmodium can also be more compatible with certain mosquitoes. For example, P. yoelii 17XNL is much more compatible with An. stephensi (Nijmegen Sda500 strain) than with An. gambiae (G3). TEP1 silencing in An. gambiae (Keele strain) mosquitoes enhances infection with P. falciparum (NK54 strain), doubling the median number of oocysts . Silencing TEP1 in An. gambiae has a more dramatic effect (4–5 fold increase) on P. berghei infection . Furthermore, silencing TEP1 in An. gambiae (G3 strain) does not enhance infection with P. falciparum (NF54 strain), indicating that there are differences in compatibility between particular strains of An. gambiae and P. falciparum (M. Povelones and A. Molina-Cruz, unpublished).
Over activation of the Rel2 pathway by silencing Caspar, a critical suppressor of this cascade, drastically reduces P. falciparum (NK54 strain) infection in An. gambiae (Keele strain), An. albimanus (Santa Tecla strain) and An. stephensi mosquitoes . Double silencing experiments in An. gambiae (Keele strain) females, in which Caspar and TEP1 (or other effectors of the Rel2 pathway) were co-silenced, rescues the effect of Caspar, indicating that TEP1 is an important effector of this response. The fact that strong activation of the Rel2 pathway can very effectively prevent infection in several mosquito species that are natural vectors of P. falciparum , begs the question of why this immune response is not effective preventing disease transmission under natural field conditions.
It has been proposed that P. falciparum parasites have evolved specific mechanisms to modulate activation of the An. gambiae immune system as they adapted to their natural mosquito vector [23, 24]. The observation that P. falciparum strains from the New World, such as the Brazilian P. falciparum 7G8 strain, are melanized very effectively by the An. gambiae L35 strain but not those of African origin  adds support to the adaptation hypothesis. Recent experiments revealed that LRIM1 can also mediate immune responses against P. falciparum, because silencing this gene in An. gambiae L35 females infected with the Brazilian P. falciparum 7G8 strain completely reverts the melanization phenotype and results in live oocysts (A. Molina-Cruz, A and C. Barillas-Mury, unpublished). Selection for refractoriness to P. cynomolgy resulted in a strain of An. gambiae that is also refractory to multiple Plasmodium species. LRIM1 also mediates the antiparasitic responses of Anopheles quadriannulatus to P. berghei infection . These findings indicate that some genes, such as TEP1/LRIM1, are broad mediators of antiparasitic responses against several different Plasmodium parasites in different mosquito strains.
Under natural conditions, P. falciparum parasites must avoid or withstand the antiparasitic responses of An. gambiae to complete their life cycle and this is likely to exert selective pressure on parasite populations. For example, in Southern Mexico, three genetically distinct P. vivax populations have been identified, and experimental infections indicate that parasites are most compatible with sympatric mosquito species . The authors propose that reciprocal selection between malaria parasites and mosquito vectors has led to local adaptation of parasites to their vectors . Thus, it is likely that in well-adapted systems there is some level of immune evasion and/or suppression, and this could explain why silencing some genes involved in immunity (LRIM1, CTL4) or oxidative stress (OXR1, GSTT1 and GSTT2) in An. gambiae (G3) females, has little effect on P. falciparum (3D7 strain) infection.
There is also increasing evidence from many different studies that the interaction between Plasmodium parasites and the mosquito immune system it is a strong determinant of vectorial capacity. Nevertheless, the extent to which the mosquito immune system is effectively reducing Plasmodium infection is very variable, even between particular parasite and mosquito strains. These differences in compatibility need to be evaluated and carefully considered when selecting an experimental animal model to study malaria transmission.
An. gambiae (G3 strain) and An. stephensi (Nijmegen Sda500) mosquitoes were raised at 28°C, 75% humidity under a 12-hour light/dark cycle and maintained on a 10% sucrose solution during adult stages.
Either wild-type or GFP-P. berghei (ANKA 2.34 strain)  and the GFP-P. yoelii yoelii 17X nonlethal transgenic strain  were maintained by serial passage in 3- to 4-week-old female BALB/c mice or as frozen stocks. Mice parasitemias were monitored by light microscopy using air-dried blood smears that were methanol fixed and stained with 10% Giemsa. Female mosquitoes (4–5 days old) were fed on gametocytemic mice 2–3 days after blood inoculation from infected donor mice when parasitemias were between 5–10%. Mosquitoes infected with P. berghei or P. yoelii were kept at 21°C or 24°C, respectively, and midguts dissected 6–7 days post infection. Infection levels were determined by fluorescent (live oocyst) and light (melanized parasites) microscopy. The distribution of oocyst numbers in the different experimental groups was compared using the nonparametric Kolmogorov-Smirnov statistical test.
Individual midguts (without blood) were placed into microcentrifuge tubes containing 10 μl of HotSHOT alkaline lysis reagent (25 mM NaOH, 0.2 mM EDTA, pH 12.0) . The tubes were boiled for 5 min and immediately placed on ice; 10 μl of HotSHOT neutralizing reagent (40 mM Tris-HCl, pH 5.0) was added to each tube. The samples were centrifuged and stored at -20°C.
For the GSTT1 silencing experiment, mice were infected wild-type P. berghei (non-GFP parasites, Anka 2.34 parasites), and the level of infection in mosquitoes was determined by qPCR 6 days post infection. Genomic DNA was obtained from infected midguts, and the abundance of P. berghei 28S RNA relative to An. gambiae S7 ribosomal protein was determined. The DyNAmo SYBR Green qPCR Master mix (Finnzymes, Espoo, Finland) was used to amplify the genomic DNA, and samples were run on a MJ Research Detection system according to the manufacturer's instructions (Bio-Rad, Hercules, CA). P. berghei 28S RNA primer sequence (5/ to 3/), Fw-GTGGCCTATCGATCCTTTA and Rev: 5/GCGTCCCAATGA TAGGAAGA). Two μl of midgut genomic DNA was used to detect the number P. berghei 28S gene copies and 1 μl to determine the copies of An. gambiae ribosomal protein S7 gene in a 20-μl PCR reaction. All P. berghei 28S values shown were then normalized relative to the number of copies of S7 in the sample. The distribution of parasite/midgut genome in control (dsLacZ injected) and dsGSTT2 silenced were compared using the Kolmogorov-Smirnov test.
An. gambiae (G3) female mosquitoes were infected with P. falciparum by feeding them gametocyte cultures using an artificial membrane feeding system. The P. falciparum (3D7 strain) was maintained in O+ human erythrocytes using RPMI 1640 medium supplemented with 25 mM HEPES, 50 mg/L hypoxanthine, 25 mM NaHCO3, and 10%(v/v) heat-inactivated type O+ human serum [30, 31]. Gametocytogenesis was induced following the procedure of Ifediba and Vanderberg . Mature gametocyte cultures (stages IV and V) that were 14–16 days old were used to feed mosquitoes in 37°C warmed membrane feeders for 30 minutes. To determine the level of infection, the midguts were dissected and stained with 0.05% (w/v) mercurochrome in water and oocysts counted by light microscopy 7–9 days post blood feeding. Distribution of oocyst numbers per midgut was analyzed using the Kolmogorov-Smirnov test.
cDNA fragments of 500–600 bp were amplified for each gene using the primers shown in Additional File 1 and cDNA from 4-day-old An. gambiae females as template. The cDNA fragments were cloned into the pCR II-TOPO® vector (Invitrogen, Carlsbad, CA) and T7 sites introduced at both ends using the following vector primers (5' to 3') to amplify the cDNA insert; M13-Fw: GTAAAACGACGGCCAGT and T7-M13Rev: CTCGAGTAATACGACTCACTA TAGGGCAGGAAACAGCTATGAC. dsRNA was synthesized and purified using the MEGAscript kit (Ambion, Austin, TX). The eluted dsRNA was further cleaned and concentrated to 3 μg/μl using a Microcon YM-100 filter (Millipore, Bedford, MA).
dsRNA (207 ng in 69 nl) for each of the genes tested was injected into the thorax of cold-anesthetized 1- to 2-day-old female mosquitoes using a nano-injector (Nanoject; Drummond Scientific, Broomall, PA). In each experiment, a control group was injected with dsLacZ or dsGFP to serve as reference for intensity of infection. Gene silencing was confirmed 4 days after dsRNA injection by RT-qPCR using the ribosomal S7 gene for normalization. Poly(A) mRNA was isolated from groups of 10 adult females using Oligotex-dT beads (Qiagen, Valencia, CA) following the manufacturer's instructions. First-strand cDNA was synthesized using random hexamers and Superscript II reverse transcriptase (Invitrogen). The primers used for each gene are shown in Additional File 2. Gene expression was assessed by SYBR green qPCR (DyNAmo HS; New England Biolabs, Beverly, MA) in a Chromo4 system (Bio-Rad). PCR involved an initial denaturation at 95°C for 15 minutes, 44 cycles of 10 seconds at 94°C, 20 seconds at 58°C, and 30 seconds at 72°C. Fluorescence readings were taken at 72°C after each cycle. A final extension at 72°C for 5 minutes was completed before deriving a melting curve (70°C–95°C) to confirm the identity of the PCR product. qPCR measurements were made in duplicate.
Because all the genes tested are highly conserved across species, we tested whether it was possible to silence An. stephensi genes by injecting them with dsRNA from orthologous genes of An. gambiae. An. stephensi female mosquitoes (1–2 days old) were injected with dsRNA from An. gambiae cDNAs following the same procedure described above. Silencing efficiency was determined using qPCR 4 days after mosquitoes were injected with dsRNA. For the initial evaluation, the same primers and conditions as for An. gambiae were used, except for a lower annealing temperature (52°C instead of 58°C). For OXR1, a strong peak was obtained using the same primers as for An. gambiae, but for all other genes, several primer combinations from well conserved regions had to be designed to obtain efficient amplification that generated a single band of the expected molecular weight. For GSTT1, in was necessary to clone a fragment of An stephensi cDNA using the following degenerate primers (5/ to 3/), Fwd: CTGGCGGAAAGT GTKGCCAT and Rev: GGCCGCAGCCASACGTACTGGAA. A 180-bp fragment was amplified, sequenced, and used to generate a primer combination that would efficiently amplify AsGSTT1. Sequences of all primer sets used for qRT-PCR analysis with An. stephensi templates are shown in Additional File 3. Silencing efficiency in An. gambiae and An. stephensi, shown in Additional File 4, ranged from 55–98% and from 56–84%, respectively.
Anopheles Plasmodium-responsive leucine-rich repeat 1
C-type lectin 4
differential interference contrast
C-type lectin 4 dsRNA-injected mosquitoes
glutathione-S-transferase theta 1 dsRNA-injected mosquitoes
glutathione-S-transferase theta 2 dsRNA-injected mosquitoes
heat-shock cognate-3dsRNA-injected mosquitoes
leucine-rich repeat immune protein 1 dsRNA-injected mosquitoes
leucine-rich repeat immune protein 2 dsRNA-injected mosquitoes dsOXR1
solute transporter dsRNA-injected mosquitoes
thioester-containing protein 1 dsRNA-injected mosquitoes
tetraspanin dsRNA-injected mosquitoes
Plasmodium yoelii yoelii 17X nonlethal transgenic strain constitutively expressing green fluorescent protein
gene family, glutathione-S-transferase of the theta class gene family
glutathione-S-transferase theta 1
glutathione-S-transferase theta 2
leucine-rich repeat immune protein 1
leucine-rich repeat immune protein 2
oxidation resistance 1
polymerase chain reaction
quantitative real-time PCR
quantitative real-time reverse-transcriptase PCR
reactive oxygen species
Royal Park Memorial Institute
protein from the small ribosomal subunit S7
thioester-containing protein 1
We thank André Laughinghouse, Kevin Lee, Tovi Lehman, and Robert Gwadz for insectary support and NIAID Intramural editor Brenda Rae Marshall. This research was supported by the Intramural Research Program of the Division of Intramural Research National Institute of Allergy and Infectious Diseases, National Institutes of Health.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.