Molecular Evolution of Immune Genes in the Malaria Mosquito Anopheles gambiae

Background

As pathogens that circumvent the host immune response are favoured by selection, so are host alleles that reduce parasite load. Such evolutionary processes leave their signature on the genes involved. Deciphering modes of selection operating on immune genes might reveal the nature of host-pathogen interactions and factors that govern susceptibility in host populations. Such understanding would have important public health implications.

Methodology/Findings

We analyzed polymorphisms in four mosquito immune genes (SP14D1, GNBP, defensin, and gambicin) to decipher selection effects, presumably mediated by pathogens. Using samples of Anopheles arabiensis, An. quadriannulatus and four An. gambiae populations, as well as published sequences from other Culicidae, we contrasted patterns of polymorphisms between different functional units of the same gene within and between populations. Our results revealed selection signatures operating on different time scales. At the most recent time scale, within-population diversity revealed purifying selection. Between populations and between species variation revealed reduced differentiation (GNBP and gambicin) at coding vs. noncoding- regions, consistent with balancing selection. McDonald-Kreitman tests between An. quadriannulatus and both sibling species revealed higher fixation rate of synonymous than nonsynonymous substitutions (GNBP) in accordance with frequency dependent balancing selection. At the longest time scale (>100 my), PAML analysis using distant Culicid taxa revealed positive selection at one codon in gambicin. Patterns of genetic variation were independent of exposure to human pathogens.

Significance and Conclusions

Purifying selection is the most common form of selection operating on immune genes as it was detected on a contemporary time scale on all genes. Selection for “hypervariability” was not detected, but negative balancing selection, detected at a recent evolutionary time scale between sibling species may be rather common. Detection of positive selection at the deepest evolutionary time scale suggests that it occurs infrequently, possibly in association with speciation events. Our results provided no evidence to support the hypothesis that selection was mediated by pathogens that are transmitted to humans.     

Distinct Olfactory Signaling Mechanisms in the Malaria Vector Mosquito Anopheles gambiae

Anopheles gambiae is the principal Afrotropical vector for human malaria, in which olfaction mediates a wide range of both adult and larval behaviors. Indeed, mosquitoes depend on the ability to respond to chemical cues for feeding, host preference, and mate location/selection. Building upon previous work that has characterized a large family of An. gambiae odorant receptors (AgORs), we now use behavioral analyses and gene silencing to examine directly the role of AgORs, as well as a newly identified family of candidate chemosensory genes, the An. gambiae variant ionotropic receptors (AgIRs), in the larval olfactory system. Our results validate previous studies that directly implicate specific AgORs in behavioral responses to DEET as well as other odorants and reveal the existence of at least two distinct olfactory signaling pathways that are active in An. gambiae. One system depends directly on AgORs; the other is AgOR-independent and requires the expression and activity of AgIRs. In addition to clarifying the mechanistic basis for olfaction in this system, these advances may ultimately enhance the development of vector control strategies, targeting olfactory pathways in mosquitoes to reduce the catastrophic effects of malaria and other mosquito-borne diseases.     

The Major Yolk Protein Vitellogenin Interferes with the Anti-Plasmodium Response in the Malaria Mosquito Anopheles gambiae

When taking a blood meal on a person infected with malaria, female Anopheles gambiae mosquitoes, the major vector of human malaria, acquire nutrients that will activate egg development (oogenesis) in their ovaries. Simultaneously, they infect themselves with the malaria parasite. On traversing the mosquito midgut epithelium, invading Plasmodium ookinetes are met with a potent innate immune response predominantly controlled by mosquito blood cells. Whether the concomitant processes of mosquito reproduction and immunity affect each other remains controversial. Here, we show that proteins that deliver nutrients to maturing mosquito oocytes interfere with the antiparasitic response. Lipophorin (Lp) and vitellogenin (Vg), two nutrient transport proteins, reduce the parasite-killing efficiency of the antiparasitic factor TEP1. In the absence of either nutrient transport protein, TEP1 binding to the ookinete surface becomes more efficient. We also show that Lp is required for the normal expression of Vg, and for later Plasmodium development at the oocyst stage. Furthermore, our results uncover an inhibitory role of the Cactus/REL1/REL2 signaling cassette in the expression of Vg, but not of Lp. We reveal molecular links that connect reproduction and immunity at several levels and provide a molecular basis for a long-suspected trade-off between these two processes.     

Chromosome Inversions,Genomic Differentiation and Speciation in the African Malaria Mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.     

Molecular Characterization of Larval Peripheral Thermosensory Responses of the Malaria Vector Mosquito Anopheles gambiae

Thermosensation provides vital inputs for the malaria vector mosquito, Anopheles gambiae which utilizes heat-sensitivity within a broad spectrum of behaviors, most notably, the localization of human hosts for blood feeding. In this study, we examine thermosensory behaviors in larval-stage An. gambiae, which as a result of their obligate aquatic habitats and importance for vectorial capacity, represents an opportunistic target for vector control as part of the global campaign to eliminate malaria. As is the case for adults, immature mosquitoes respond differentially to a diverse array of external heat stimuli. In addition, larvae exhibit a striking phenotypic plasticity in thermal-driven behaviors that are established by temperature at which embryonic development occurs. Within this spectrum, RNAi-directed gene-silencing studies provide evidence for the essential role of the Transient Receptor Potential sub-family A1 (TRPA1) channel in mediating larval thermal-induced locomotion and thermal preference within a discrete upper range of ambient temperatures.     

The CLIP-Domain Serine Protease Homolog SPCLIP1 Regulates Complement Recruitment to Microbial Surfaces in the Malaria Mosquito Anopheles gambiae

The complement C3-like protein TEP1 of the mosquito Anopheles gambiae is required for defense against malaria parasites and bacteria. Two forms of TEP1 are present in the mosquito hemolymph, the full-length TEP1-F and the proteolytically processed TEP1_cut that is part of a complex including the leucine-rich repeat proteins LRIM1 and APL1C. Here we show that the non-catalytic serine protease SPCLIP1 is a key regulator of the complement-like pathway. SPCLIP1 is required for accumulation of TEP1 on microbial surfaces, a reaction that leads to lysis of malaria parasites or triggers activation of a cascade culminating with melanization of malaria parasites and bacteria. We also demonstrate that the two forms of TEP1 have distinct roles in the complement-like pathway and provide the first evidence for a complement convertase-like cascade in insects analogous to that in vertebrates. Our findings establish that core principles of complement activation are conserved throughout the evolution of animals.     

A New Role of the Mosquito Complement-like Cascade in Male Fertility in Anopheles gambiae

Thioester-containing protein 1 (TEP1) is a key immune factor that determines mosquito resistance to a wide range of pathogens, including malaria parasites. Here we report a new allele-specific function of TEP1 in male fertility. We demonstrate that during spermatogenesis TEP1 binds to and removes damaged cells through the same complement-like cascade that kills malaria parasites in the mosquito midgut. Further, higher fertility rates are mediated by an allele that renders the mosquito susceptible to Plasmodium. By elucidating the molecular and genetic mechanisms underlying TEP1 function in spermatogenesis, our study suggests that pleiotropic antagonism between reproduction and immunity may shape resistance of mosquito populations to malaria parasites.     

The Interaction between a Sexually Transferred Steroid Hormone and a Female Protein Regulates Oogenesis in the Malaria Mosquito Anopheles gambiae

Molecular interactions between male and female factors during mating profoundly affect the reproductive behavior and physiology of female insects. In natural populations of the malaria mosquito Anopheles gambiae, blood-fed females direct nutritional resources towards oogenesis only when inseminated. Here we show that the mating-dependent pathway of egg development in these mosquitoes is regulated by the interaction between the steroid hormone 20-hydroxy-ecdysone (20E) transferred by males during copulation and a female Mating-Induced Stimulator of Oogenesis (MISO) protein. RNAi silencing of MISO abolishes the increase in oogenesis caused by mating in blood-fed females, causes a delay in oocyte development, and impairs the function of male-transferred 20E. Co-immunoprecipitation experiments show that MISO and 20E interact in the female reproductive tract. Moreover MISO expression after mating is induced by 20E via the Ecdysone Receptor, demonstrating a close cooperation between the two factors. Male-transferred 20E therefore acts as a mating signal that females translate into an increased investment in egg development via a MISO-dependent pathway. The identification of this male–female reproductive interaction offers novel opportunities for the control of mosquito populations that transmit malaria.     

CYP6 P450 Enzymes and ACE-1 Duplication Produce Extreme and Multiple Insecticide Resistance in the Malaria Mosquito Anopheles gambiae

Malaria control relies heavily on pyrethroid insecticides, to which susceptibility is declining in Anopheles mosquitoes. To combat pyrethroid resistance, application of alternative insecticides is advocated for indoor residual spraying (IRS), and carbamates are increasingly important. Emergence of a very strong carbamate resistance phenotype in Anopheles gambiae from Tiassalé, Côte d''Ivoire, West Africa, is therefore a potentially major operational challenge, particularly because these malaria vectors now exhibit resistance to multiple insecticide classes. We investigated the genetic basis of resistance to the most commonly-applied carbamate, bendiocarb, in An. gambiae from Tiassalé. Geographically-replicated whole genome microarray experiments identified elevated P450 enzyme expression as associated with bendiocarb resistance, most notably genes from the CYP6 subfamily. P450s were further implicated in resistance phenotypes by induction of significantly elevated mortality to bendiocarb by the synergist piperonyl butoxide (PBO), which also enhanced the action of pyrethroids and an organophosphate. CYP6P3 and especially CYP6M2 produced bendiocarb resistance via transgenic expression in Drosophila in addition to pyrethroid resistance for both genes, and DDT resistance for CYP6M2 expression. CYP6M2 can thus cause resistance to three distinct classes of insecticide although the biochemical mechanism for carbamates is unclear because, in contrast to CYP6P3, recombinant CYP6M2 did not metabolise bendiocarb in vitro. Strongly bendiocarb resistant mosquitoes also displayed elevated expression of the acetylcholinesterase ACE-1 gene, arising at least in part from gene duplication, which confers a survival advantage to carriers of additional copies of resistant ACE-1 G119S alleles. Our results are alarming for vector-based malaria control. Extreme carbamate resistance in Tiassalé An. gambiae results from coupling of over-expressed target site allelic variants with heightened CYP6 P450 expression, which also provides resistance across contrasting insecticides. Mosquito populations displaying such a diverse basis of extreme and cross-resistance are likely to be unresponsive to standard insecticide resistance management practices.     

Behavioral Cost & Overdominance in Anopheles gambiae

In response to the widespread use of control strategies such as Insecticide Treated Nets (ITN), Anopheles mosquitoes have evolved various resistance mechanisms. Kdr is a mutation that provides physiological resistance to the pyrethroid insecticides family (PYR). In the present study, we investigated the effect of the Kdr mutation on the ability of female An. gambiae to locate and penetrate a 1cm-diameter hole in a piece of netting, either treated with insecticide or untreated, to reach a bait in a wind tunnel. Kdr homozygous, PYR-resistant mosquitoes were the least efficient at penetrating an untreated damaged net, with about 51% [39-63] success rate compared to 80% [70-90] and 78% [65-91] for homozygous susceptible and heterozygous respectively. This reduced efficiency, likely due to reduced host-seeking activity, as revealed by mosquito video-tracking, is evidence of a recessive behavioral cost of the mutation. Kdr heterozygous mosquitoes were the most efficient at penetrating nets treated with PYR insecticide, thus providing evidence for overdominance, the rarely-described case of heterozygote advantage conveyed by a single locus. The study also highlights the remarkable capacity of female mosquitoes, whether PYR-resistant or not, to locate holes in bed-nets.     

Melanotic Pathology and Vertical Transmission of the Gut Commensal Elizabethkingia meningoseptica in the Major Malaria Vector Anopheles gambiae

Background

The resident gut flora is known to have significant impacts on the life history of the host organism. Endosymbiotic bacterial species in the Anopheles mosquito gut are potent modulators of sexual development of the malaria parasite, Plasmodium, and thus proposed as potential control agents of malaria transmission.

Results

Here we report a melanotic pathology in the major African malaria vector Anopheles gambiae, caused by the dominant mosquito endosymbiont Elizabethkingia meningoseptica . Transfer of melanised tissues into the haemolymph of healthy adult mosquitoes or direct haemolymph inoculation with isolated E . meningoseptica bacteria were the only means for transmission and de novo formation of melanotic lesions, specifically in the fat body tissues of recipient individuals. We show that E . meningoseptica can be vertically transmitted from eggs to larvae and that E . meningoseptica -mono-associated mosquitoes display significant mortality, which is further enhanced upon Plasmodium infection, suggesting a synergistic impact of E . meningoseptica and Plasmodium on mosquito survival.

Conclusion

The high pathogenicity and permanent association of E . meningoseptica with An. Gambiae through vertical transmission constitute attractive characteristics towards the potential design of novel mosquito/malaria biocontrol strategies.     

Evidence for a Lectin Specific for Sulfated Glycans in the Salivary Gland of the Malaria Vector,Anopheles gambiae

Salivary gland homogenate (SGH) from the female mosquitoes Anopheles gambiae, An. stephensi, An. freeborni, An. dirus and An. albimanus were found to exhibit hemagglutinating (lectin) activity. Lectin activity was not found for male An. gambiae, or female Ae aegypti, Culex quinquefasciatus, Phlebotomus duboscqi, and Lutzomyia longipalpis. With respect to species-specificity, An. gambiae SGH agglutinates red blood cells (RBC) from humans, horse, sheep, goat, pig, and cow; it is less active for rats RBC, and not detectable for guinea-pigs or chicken RBC. Notably, lectin activity was inhibited by low concentrations of dextran sulfate 50–500 K, fucoidan, heparin, laminin, heparin sulfate proteoglycan, sialyl-containing glycans (e.g. 3′-sialyl Lewis X, and 6′-sialyl lactose), and gangliosides (e.g. GM3, GD1, GD1b, GTB1, GM1, GQ1B), but not by simple sugars. These results imply that molecule(s) in the salivary gland target sulfated glycans. SGH from An. gambiae was also found to promote agglutination of HL-60 cells which are rich in sialyl Lewis X, a glycan that decorates PSGL-1, the neutrophils receptor that interacts with endothelial cell P-selectin. Accordingly, SGH interferes with HL-60 cells adhesion to immobilized P-selectin. Because An. gambiae SGH expresses galectins, one member of this family (herein named Agalectin) was expressed in E. coli. Recombinant Agalectin behaves as a non-covalent homodimer. It does not display lectin activity, and does not interact with 500 candidates tested in a Glycan microarray. Gel-filtration chromatography of the SGH of An. gambiae identified a fraction with hemagglutinating activity, which was analyzed by 1D PAGE followed by in-gel tryptic digestion, and nano-LC MS/MS. This approach identified several genes which emerge as candidates for a lectin targeting sulfated glycans, the first with this selectivity to be reported in the SGH of a blood-sucking arthropod. The role of salivary molecules (sialogenins) with lectin activity is discussed in the context of inflammation, and parasite-vector-host interactions.     

A Proteomic Investigation of Soluble Olfactory Proteins in Anopheles gambiae

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are small soluble polypeptides that bind semiochemicals in the lymph of insect chemosensilla. In the genome of Anopheles gambiae, 66 genes encode OBPs and 8 encode CSPs. Here we monitored their expression through classical proteomics (2D gel-MS analysis) and a shotgun approach. The latter method proved much more sensitive and therefore more suitable for tiny biological samples as mosquitoes antennae and eggs. Females express a larger number and higher quantities of OBPs in their antennae than males (24 vs 19). OBP9 is the most abundant in the antennae of both sexes, as well as in larvae, pupae and eggs. Of the 8 CSPs, 4 were detected in antennae, while SAP3 was the only one expressed in larvae. Our proteomic results are in fairly good agreement with data of RNA expression reported in the literature, except for OBP4 and OBP5, that we could not identify in our analysis, nor could we detect in Western Blot experiments. The relatively limited number of soluble olfactory proteins expressed at relatively high levels in mosquitoes makes further studies on the coding of chemical messages at the OBP level more accessible, providing for few specific targets. Identification of such proteins in Anopheles gambiae might facilitate future studies on host finding behavior in this important disease vector.     

An Integrated Genetic Map of the African Human Malaria Vector Mosquito, Anopheles Gambiae

We present a genetic map based on microsatellite polymorphisms for the African human malaria vector, Anopheles gambiae. Polymorphisms in laboratory strains were detected for 89% of the tested microsatellite markers. Genotyping was performed for individual mosquitoes from 13 backcross families that included 679 progeny. Three linkage groups were identified, corresponding to the three chromosomes. We added 22 new markers to the existing X chromosome map, for a total of 46 microsatellite markers spanning a distance of 48.9 cM. The second chromosome has 57 and the third 28 microsatellite markers spanning a distance of 72.4 and 93.7 cM, respectively. The overall average distance between markers is 1.6 cM (or 1.1, 1.2, and 3.2 cM for the X, second, and third chromosomes, respectively). In addition to the 131 microsatellite markers, the current map also includes a biochemical selectable marker, Dieldrin resistance (Dl), on the second chromosome and five visible markers, pink-eye (p) and white (w) on the X, collarless (c) and lunate (lu) on the second, and red-eye (r) on the third. The cytogenetic locations on the nurse cell polytene chromosomes have been determined for 47 markers, making this map an integrated tool for cytogenetic, genetic, and molecular analysis.     

Genetic Dissection of Anopheles gambiae Gut Epithelial Responses to Serratia marcescens

Genetic variation in the mosquito Anopheles gambiae profoundly influences its ability to transmit malaria. Mosquito gut bacteria are shown to influence the outcome of infections with Plasmodium parasites and are also thought to exert a strong drive on genetic variation through natural selection; however, a link between antibacterial effects and genetic variation is yet to emerge. Here, we combined SNP genotyping and expression profiling with phenotypic analyses of candidate genes by RNAi-mediated silencing and 454 pyrosequencing to investigate this intricate biological system. We identified 138 An. gambiae genes to be genetically associated with the outcome of Serratia marcescens infection, including the peptidoglycan recognition receptor PGRPLC that triggers activation of the antibacterial IMD/REL2 pathway and the epidermal growth factor receptor EGFR. Silencing of three genes encoding type III fibronectin domain proteins (FN3Ds) increased the Serratia load and altered the gut microbiota composition in favor of Enterobacteriaceae. These data suggest that natural genetic variation in immune-related genes can shape the bacterial population structure of the mosquito gut with high specificity. Importantly, FN3D2 encodes a homolog of the hypervariable pattern recognition receptor Dscam, suggesting that pathogen-specific recognition may involve a broader family of immune factors. Additionally, we showed that silencing the gene encoding the gustatory receptor Gr9 that is also associated with the Serratia infection phenotype drastically increased Serratia levels. The Gr9 antibacterial activity appears to be related to mosquito feeding behavior and to mostly rely on changes of neuropeptide F expression, together suggesting a behavioral immune response following Serratia infection. Our findings reveal that the mosquito response to oral Serratia infection comprises both an epithelial and a behavioral immune component.     

Embryonic Development and Rates of Metabolic Activity in Early and Late Hatching Eggs of the Major Malaria Vector Anopheles gambiae

Anopheles gambiae eggs generally hatch at the completion of embryo development; two-three days post oviposition. However, staggered or delayed hatching has been observed whereby a single batch of eggs shows marked variation in time-to-hatch, with some eggs hatching 18 days post oviposition or later. The mechanism enabling delayed hatch has not been clearly elucidated but is likely mediated by environmental and genetic factors that either induce diapause or slow embryo development. This study aimed to compare metabolic activity and embryonic development between eggs collected from sub-colonies of the baseline Anopheles gambiae GAH colony previously selected for early or late time-to-hatch. Egg batches from early and late hatch sub-colonies as well as from the baseline colony were monitored for hatching. For both time-to-hatch selected sub-colonies and the baseline colony the majority of eggs hatched on day two post oviposition. Nevertheless, eggs produced by the late hatch sub-colony showed a significantly longer mean time to hatch than those produced by the early hatch sub-colony. The overall proportions that hatched were similar for all egg batches. CO₂ output between eggs from early and late hatch sub-colonies showed significant differences only at 3 and 7 days post oviposition where eggs from the early hatch and the late hatch sub-colony were more metabolically active, respectively. No qualitative differences were observed in embryo development between the sub-colonies. It is concluded that all viable embryos develop to maturity at the same rate and that a small proportion then enter a state of diapause enabling them to hatch later. As it has previously been shown that it is possible to at least partially select for late hatch, this characteristic is likely to involve genetic as well as environmental factors. Delayed hatching in An. gambiae is likely an adaptation to maximise reproductive output despite the increased risk of desiccation in an unstable aquatic environment.