Polymorphism at the defensin gene in the Anopheles gambiae complex: testing different selection hypotheses

Phylogeny of Anopheles darlingi (Diptera:Culicidae) based on the antimicrobial peptide genes cecropin and defensin

Cecropins and defensins are the main classes of antimicrobial peptides in the mosquito innate immune system, acting against bacteria, fungi and protozoa. There is a knowledge gap concerning these peptide genes in anopheline mosquitoes from the Brazilian Amazon. Thus, this work aimed to describe molecular techniques for detecting the genes encoding the antimicrobial peptides cecropin A (CecA) and defensin in Anopheles darlingi mosquitoes and to perform molecular phylogeny of the sequenced genes using the maximum likelihood method and Bayesian inference with other species from different geographic areas. Our results show, for the first time, a molecular biology method for detecting CecA and defensin in Anopheles darlingi that allows for the use of these molecular markers for phylogenetic analysis in anopheline species, separating the species into single and monophyletic clades.

Molecular evolution and population genetics of a Gram-negative binding protein gene in the malaria vector Anopheles gambiae (sensu lato)

BACKGROUND

Clarifying the role of the innate immune system of the malaria vector Anopheles gambiae is a potential way to block the development of the Plasmodium parasites. Pathogen recognition is the first step of innate immune response, where pattern recognition proteins like GNBPs play a central role.

RESULTS

We analysed 70 sequences of the protein coding gene GNBPB2 from two species, Anopheles gambiae (s.s.) and An. coluzzii, collected in six African countries. We detected 135 segregating sites defining 63 distinct haplotypes and 30 proteins. Mean nucleotide diversity (π) was 0.014 for both species. We found no significant genetic differentiation between species, but a significant positive correlation between genetic differentiation and geographical distance among populations.

CONCLUSIONS

Species status seems to contribute less for the molecular differentiation in GNBPB2 than geographical region in the African continent (West and East). Purifying selection was found to be the most common form of selection, as in many other immunity-related genes. Diversifying selection may be also operating in the GNBPB2 gene.

The potential for adaptive maintenance of diversity in insect antimicrobial peptides

Genes involved in immune defence are among the fastest evolving in the genomes of many species. Interestingly, however, genes encoding antimicrobial peptides (AMPs) have shown little evidence for adaptive divergence in arthropods, despite the centrality of these peptides in direct killing of microbial pathogens. This observation, coupled with a failure to detect phenotypic consequence of genetic variation in AMPs, has led to the hypothesis that individual AMPs make minor contributions to overall immune defence and that AMPs instead act as a collective cocktail. Recent data, however, have suggested an alternative explanation for the apparent lack of adaptive divergence in AMP genes. Molecular evolutionary and phenotypic data have begun to suggest that variant AMP alleles may be maintained through balancing selection in invertebrates, a pattern similar to that observed in several vertebrate AMPs. Signatures of balancing selection include high rates of non-synonymous polymorphism, trans-species amino acid polymorphisms, and convergence of amino acid states across the phylogeny. In this review, we revisit published literature on insect AMP genes and analyse newly available population genomic datasets in Drosophila, finding enrichment for patterns consistent with adaptive maintenance of polymorphism.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

A genomic comparison of two termites with different social complexity

The termites evolved eusociality and complex societies before the ants, but have been studied much less. The recent publication of the first two termite genomes provides a unique comparative opportunity, particularly because the sequenced termites represent opposite ends of the social complexity spectrum. Zootermopsis nevadensis has simple colonies with totipotent workers that can develop into all castes (dispersing reproductives, nest-inheriting replacement reproductives, and soldiers). In contrast, the fungus-growing termite Macrotermes natalensis belongs to the higher termites and has very large and complex societies with morphologically distinct castes that are life-time sterile. Here we compare key characteristics of genomic architecture, focusing on genes involved in communication, immune defenses, mating biology and symbiosis that were likely important in termite social evolution. We discuss these in relation to what is known about these genes in the ants and outline hypothesis for further testing.

Genetic polymorphism at an odorant receptor gene (Or39) among mosquitoes of the Anopheles gambiae complex in Senegal (West Africa)

BACKGROUND

Olfaction plays a significant role in insect behavior during critical steps of their life-cycle, such as host-seeking during foraging or the search for a mate. Here, we explored genetic polymorphism within and divergence between sibling species of the African malaria mosquito, Anopheles gambiae sensu lato in the gene sequence and encoded peptides of an odorant receptor, Or39. This study included sympatric specimens of An. gambiae sensu stricto, An. coluzzii and An. arabiensis sampled together in the village of Dielmo, Senegal.

RESULTS

A 1,601 bp genomic sequence composed of 6 exons and 5 introns was obtained for Or39 from 6-8 mosquitoes in each of the 3 species. DNA sequence analysis revealed a high level of molecular polymorphism (π = 0.0154; Haplotype diversity = 0.867) and high overall genetic differentiation between taxa (Fst > 0.92, P < 0.01). In total, 50 parsimony informative sites were recorded. Throughout the whole dataset, there were 13 non-synonymous mutations resulting in aminoacid changes in the encoded protein. Each of the 6 different identified peptides was species-specific and none was shared across species. Most aminoacid changes were located on the intracellular domains of the protein. However, intraspecific polymorphisms in An. gambiae and An. arabiensis as well as species-specific mutations also occurred in the first extracellular domain.

CONCLUSIONS

Although obtained from a limited number of specimens, our results point towards genetic differences between cryptic species within the An. gambiae complex in a gene of biological relevance that might be of evolutionary significance when exposed to disruptive selective forces.

Phylogeography of Aedes aegypti (yellow fever mosquito) in South Florida: mtDNA evidence for human-aided dispersal

The invasive dengue vector Aedes aegypti has persisted for > 200 years in South Florida in the United States. We tested the hypotheses that Florida's landscape creates dispersal barriers and corridors and that long-distance human-aided dispersal structures populations of Ae. aegypti. We evaluated the phylogeography of 362 individuals from Florida's East and West Coasts with a 760-bp (418- and 342-bp fragments of ND5 and ND4, respectively) mitochondrial sequence. Populations from these two coasts were not significantly differentiated, suggesting that limited urbanization in central Florida is not a strong barrier to gene flow. Evidence for long-distance dispersal between Ft. Lauderdale and the West and Ft. Myers and the East indicates the importance of human-aided dispersal. West Coast populations showed no genetic differentiation, indicating that West Coast rivers and bays did not significantly impede gene flow. Phylogeographic analysis of haplotypes showed two distinct matrilines with no geographic patterns, suggesting multiple introductions or balancing selection.

Evolution of Anopheles gambiae in Relation to Humans and Malaria

The closely related and morphologically indistinguishable mosquito species in the Afrotropical Anopheles gambiae complex differ dramatically in their contribution to malaria transmission, ranging from major vectors through minor or locally important vectors and nonvectors. Radiation of the A. gambiae complex and ongoing diversification within its nominal species appears to be a product of recent and rapid adaptation to environmental heterogeneities, notably those of anthropogenic origin. Polytene chromosome and genomic analyses suggest that paracentric chromosomal inversions and possibly other low-recombination regions have played instrumental roles in this process by facilitating ecotypic differentiation both within and across semipermeable species boundaries. Forthcoming complete genome sequences from several members of the A. gambiae complex will provide powerful tools to accelerate ongoing investigation of how genetic diversification of populations and species has shaped behavioral and physiological traits, such as vector competence, that bear on vectorial importance.

Molecular evolution of a gene cluster of serine proteases expressed in the Anopheles gambiae female reproductive tract

Background

Genes involved in post-mating processes of multiple mating organisms are known to evolve rapidly due to coevolution driven by sexual conflict among male-female interacting proteins. In the malaria mosquito Anopheles gambiae - a monandrous species in which sexual conflict is expected to be absent or minimal - recent data strongly suggest that proteolytic enzymes specifically expressed in the female lower reproductive tissues are involved in the processing of male products transferred to females during mating. In order to better understand the role of selective forces underlying the evolution of proteins involved in post-mating responses, we analysed a cluster of genes encoding for three serine proteases that are down-regulated after mating, two of which specifically expressed in the atrium and one in the spermatheca of A. gambiae females.

Results

The analysis of polymorphisms and divergence of these female-expressed proteases in closely related species of the A. gambiae complex revealed a high level of replacement polymorphisms consistent with relaxed evolutionary constraints of duplicated genes, allowing to rapidly fix novel replacements to perform new or more specific functions. Adaptive evolution was detected in several codons of the 3 genes and hints of episodic selection were also found. In addition, the structural modelling of these proteases highlighted some important differences in their substrate specificity, and provided evidence that a number of sites evolving under selective pressures lie relatively close to the catalytic triad and/or on the edge of the specificity pocket, known to be involved in substrate recognition or binding. The observed patterns suggest that these proteases may interact with factors transferred by males during mating (e.g. substrates, inhibitors or pathogens) and that they may have differently evolved in independent A. gambiae lineages.

Conclusions

Our results - also examined in light of constraints in the application of selection-inference methods to the closely related species of the A. gambiae complex - reveal an unexpectedly intricate evolutionary scenario. Further experimental analyses are needed to investigate the biological functions of these genes in order to better interpret their molecular evolution and to assess whether they represent possible targets for limiting the fertility of Anopheles mosquitoes in malaria vector control strategies.

Exceptional diversity, maintenance of polymorphism, and recent directional selection on the APL1 malaria resistance genes of Anopheles gambiae

The three-gene APL1 locus encodes essential components of the mosquito immune defense against malaria parasites. APL1 was originally identified because it lies within a mapped QTL conferring the vector mosquito Anopheles gambiae natural resistance to the human malaria parasite, Plasmodium falciparum, and APL1 genes have subsequently been shown to be involved in defense against several species of Plasmodium. Here, we examine molecular population genetic variation at the APL1 gene cluster in spatially and temporally diverse West African collections of A. gambiae. The locus is extremely polymorphic, showing evidence of adaptive evolutionary maintenance of genetic variation. We hypothesize that this variability aids in defense against genetically diverse pathogens, including Plasmodium. Variation at APL1 is highly structured across geographic and temporal subpopulations. In particular, diversity is exceptionally high during the rainy season, when malaria transmission rates are at their peak. Much less allelic diversity is observed during the dry season when mosquito population sizes and malaria transmission rates are low. APL1 diversity is weakly stratified by the polymorphic 2La chromosomal inversion but is very strongly subdivided between the M and S “molecular forms.” We find evidence that a recent selective sweep has occurred at the APL1 locus in M form mosquitoes only. The independently reported observation of a similar M-form restricted sweep at the Tep1 locus, whose product physically interacts with APL1C, suggests that epistatic selection may act on these two loci causing them to sweep coordinately.

Immune defense genes are sometimes highly variable in host populations, reflecting selective pressure to combat diverse pathogens. In other instances, where there are only a few dominant pathogens, natural selection may favor only one or a few defense alleles. Here, we show that both adaptive strategies can occur in the same genes under different circumstances. We examined diversity in the APL1 genes of the human malaria vector mosquito Anophleles gambiae, which play a role in defense against malaria parasites. We found that the APL1 genes are exceptionally polymorphic, being 10-fold more diverse than typical A. gambiae genes. The distribution of APL1 allelic diversity, however, is strongly structured depending on whether the genes are carried by the M or S “molecular forms” of the vector, which are thought to constitute newly forming species. We show that despite the evolutionary maintenance of APL1 diversity in the S form of A. gambiae, there is evidence of strong recent directional selection on APL1 genes in the M form. Independent research has shown that Tep1, a gene which encodes a protein that physically interacts with the APL1C protein, also harbors high allelic diversity in the S form and shows evidence of recent directional selection in the M form, suggesting that the evolutionary trajectories of the Tep1 and APL1 defense loci may be correlated.

Low linkage disequilibrium in wild Anopheles gambiae s.l. populations

Background

In the malaria vector Anopheles gambiae, understanding diversity in natural populations and genetic components of important phenotypes such as resistance to malaria infection is crucial for developing new malaria transmission blocking strategies. The design and interpretation of many studies here depends critically on Linkage disequilibrium (LD). For example in association studies, LD determines the density of Single Nucleotide Polymorphisms (SNPs) to be genotyped to represent the majority of the genomic information. Here, we aim to determine LD in wild An. gambiae s.l. populations in 4 genes potentially involved in mosquito immune responses against pathogens (Gambicin, NOS, REL2 and FBN9) using previously published and newly generated sequences.

Results

The level of LD between SNP pairs in cloned sequences of each gene was determined for 7 species (or incipient species) of the An. gambiae complex. In all tested genes and species, LD between SNPs was low: even at short distances (< 200 bp), most SNP pairs gave an r² < 0.3. Mean r² ranged from 0.073 to 0.766. In most genes and species LD decayed very rapidly with increasing inter-marker distance.

Conclusions

These results are of great interest for the development of large scale polymorphism studies, as LD generally falls below any useful limit. It indicates that very fine scale SNP detection will be required to give an overall view of genome-wide polymorphism. Perhaps a more feasible approach to genome wide association studies is to use targeted approaches using candidate gene selection to detect association to phenotypes of interest.

Identification of putative innate immune related genes from a cell line of the mosquito Aedes albopictus following bacterial challenge

We report identification of putative innate immune related genes from a cell line of the mosquito Aedes albopictus challenged with heat-killed bacteria. Using a subtractive hybridization and sequencing approach, we analyzed a total 309 expressed sequence tags (ESTs) which clustered in 40 contigs. Thirty-five percent of genes yielded homology to known immune genes corresponding to antimicrobial peptides (AMPs), pathogen-associated molecular patterns, protease and immune signaling cascades. Interestingly, most of the genes have not been previously described from this mosquito and thus represent a class of novel immune genes. Further, 25% sequences did not match to any known species in the non-redundant databases, appear to be specific to the mosquito A. albopictus and merit further study.

Anopheles immune genes and amino acid sites evolving under the effect of positive selection

BACKGROUND

It has long been the goal of vector biology to generate genetic knowledge that can be used to "manipulate" natural populations of vectors to eliminate or lessen disease burden. While long in coming, progress towards reaching this goal has been made. Aiming to increase our understanding regarding the interactions between Plasmodium and the Anopheles immune genes, we investigated the patterns of genetic diversity of four anti-Plasmodium genes in the Anopheles gambiae complex of species.

METHODOLOGY/PRINCIPAL FINDINGS

Within a comparative phylogenetic and population genetics framework, the evolutionary history of four innate immunity genes within the An. gambiae complex (including the two most important human malaria vectors, An. gambiae and An. arabiensis) is reconstructed. The effect of natural selection in shaping the genes' diversity is examined. Introgression and retention of ancestral polymorphisms are relatively rare at all loci. Despite the potential confounding effects of these processes, we could identify sites that exhibited dN/dS ratios greater than 1.

CONCLUSIONS/SIGNIFICANCE

In two of the studied genes, CLIPB14 and FBN8, several sites indicated evolution under positive selection, with CLIPB14 exhibiting the most consistent evidence. Considering only the sites that were consistently identified by all methods, two sites in CLIPB14 are adaptively driven. However, the analysis inferring the lineage -specific evolution of each gene was not in favor of any of the Anopheles lineages evolving under the constraints imposed by positive selection. Nevertheless, the loci and the specific amino acids that were identified as evolving under strong evolutionary pressure merit further investigation for their involvement in the Anopheles defense against microbes in general.

Molecular evolution of the three short PGRPs of the malaria vectors Anopheles gambiae and Anopheles arabiensis in East Africa

Background

Immune responses to parasites, which start with pathogen recognition, play a decisive role in the control of the infection in mosquitoes. Peptidoglycan recognition proteins (PGRPs) are an important family of pattern recognition receptors that are involved in the activation of these immune reactions. Pathogen pressure can exert adaptive changes in host genes that are crucial components of the vector's defence. The aim of this study was to determine the molecular evolution of the three short PGRPs (PGRP-S1, PGRP-S2 and PGRP-S3) in the two main African malaria vectors - Anopheles gambiae and Anopheles arabiensis.

Results

Genetic diversity of An. gambiae and An. arabiensis PGRP-S1, PGRP-S2 and PGRP-S3 was investigated in samples collected from Mozambique and Tanzania. PGRP-S1 diversity was lower than for PGRP-S2 and PGRP-S3. PGRP-S1 was the only gene differentiated between the two species. All the comparisons made for PGRP-S1 showed significant P-values for Fst estimates and AMOVA confirming a clear separation between species. For PGRP-S2 and PGRP-S3 genes it was not possible to group populations either by species or by geographic region. Phylogenetic networks reinforced the results obtained by the AMOVA and Fst values. The ratio of nonsynonymous substitutions (Ka)/synonymous substitutions (Ks) for the duplicate pair PGRP-S2 and PGRP-S3 was very similar and lower than 1. The 3D model of the different proteins coded by these genes showed that amino acid substitutions were concentrated at the periphery of the protein rather than at the peptidoglycan recognition site.

Conclusions

PGRP-S1 is less diverse and showed higher divergence between An. gambiae and An. arabiensis regardless of geographic location. This probably relates to its location in the chromosome-X, while PGRP-S2 and PGRP-S3, located in chromosome-2L, showed signs of autosomal introgression. The two short PGRP genes located in the chromosome-2L were under purifying selection, which suggests functional constraints. Different types of selection acting on PGRP-S1 and PGRP-S2 and S3 might be related to their different function and catalytic activity.

Two novel defensin-encoding genes of the Chagas disease vector Triatoma brasiliensis (Reduviidae, Triatominae): gene expression and peptide-structure modeling

Defensins are cysteine-rich peptides involved in the innate immunity of insects and many other organisms. In the present study, two novel defensin-encoding cDNAs and the respective genomic DNAs (def3 and def4) of Triatoma brasiliensis were identified and their tissue-specific and temporal expression was characterized. Both of the deduced mature peptides consisted of 43 amino acid residues and were highly similar to previously identified triatomine defensins (81.4-100%). Semi-quantitative RT-PCR data showed that def3 was constitutively expressed in the fat body and was induced in salivary glands and the small intestine at 5 and 3 days after feeding (daf), respectively. The def4 mRNA level was highly up-regulated in the stomach and fat-body tissues at 5 and 3 daf, respectively. The three-dimensional structures of these defensins were predicted using a homology modeling approach with Def-AAA, the defensin from Anopheles gambiae, as template (62-74% identity). A map of the electrostatic potential of these models revealed that, despite their similar folding patterns, mature Def2 and Def4 have a more cationic structure than is the case for Def1 and Def3. Such differences may orient the antimicrobial profile of these defensins against distinct targets in different organs of the insect.

High, clustered, nucleotide diversity in the genome of Anopheles gambiae revealed through pooled-template sequencing: implications for high-throughput genotyping protocols

Background

Association mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in Anopheles gambiae we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance.

Results

Using two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the An. gambiae genome. Comparison of SNPs identified from pooled templates with those from individual sequences revealed a very low false positive rate. False negative rates were much higher and mostly resulted from SNPs with a low minor allele frequency. Pooled-template sequencing also provided good estimates of SNP allele frequencies. Allele frequency estimation success, along with false positive and negative call rates, improved significantly when using a qualitative measure of SNP call quality. We identified a total of 7062 polymorphic features comprising 6995 SNPs and 67 indels, with, on average, a SNP every 34 bp; a high rate of polymorphism that is comparable to other studies of mosquitoes. SNPs were significantly more frequent in members of the cytochrome p450 mono-oxygenases and carboxy/cholinesterase gene-families than in glutathione-S-transferases, other detoxification genes, and control genomic regions. Polymorphic sites showed a significantly clustered distribution, but the degree of SNP clustering (independent of SNP frequency) did not vary among gene families, suggesting that clustering of polymorphisms is a general property of the An. gambiae genome.

Conclusion

The high frequency and clustering of SNPs has important ramifications for the design of high-throughput genotyping assays based on allele specific primer extension or probe hybridisation. We illustrate these issues in the context of the design of Illumina GoldenGate assays.

Inferring selection in the Anopheles gambiae species complex: an example from immune-related serine protease inhibitors

Background

Mosquitoes of the Anopheles gambiae species complex are the primary vectors of human malaria in sub-Saharan Africa. Many host genes have been shown to affect Plasmodium development in the mosquito, and so are expected to engage in an evolutionary arms race with the pathogen. However, there is little conclusive evidence that any of these mosquito genes evolve rapidly, or show other signatures of adaptive evolution.

Methods

Three serine protease inhibitors have previously been identified as candidate immune system genes mediating mosquito-Plasmodium interaction, and serine protease inhibitors have been identified as hot-spots of adaptive evolution in other taxa. Population-genetic tests for selection, including a recent multi-gene extension of the McDonald-Kreitman test, were applied to 16 serine protease inhibitors and 16 other genes sampled from the An. gambiae species complex in both East and West Africa.

Results

Serine protease inhibitors were found to show a marginally significant trend towards higher levels of amino acid diversity than other genes, and display extensive genetic structuring associated with the 2La chromosomal inversion. However, although serpins are candidate targets for strong parasite-mediated selection, no evidence was found for rapid adaptive evolution in these genes.

Conclusion

It is well known that phylogenetic and population history in the An. gambiae complex can present special problems for the application of standard population-genetic tests for selection, and this may explain the failure of this study to detect selection acting on serine protease inhibitors. The pitfalls of uncritically applying these tests in this species complex are highlighted, and the future prospects for detecting selection acting on the An. gambiae genome are discussed.

Selection on an antimicrobial peptide defensin in ants

Ants live in crowded nests with interacting individuals, which makes them particularly prone to infectious diseases. The question is, how do ants cope with the increased risk of pathogen transmission due to sociality? We have studied the molecular evolution of defensin, a gene encoding an antimicrobial protein, in ants. Defensin sequences from several ant species were analyzed with maximum likelihood models of codon substitution to infer selection. Positive selection was detected in the mature region of defensin, whereas the signal and pro regions seem to be evolving neutrally. We also found a significantly higher rate of nonsynonymous substitutions in some phylogenetic lineages, as well as dN/dS >1, suggesting varying selection pressures in different lineages. Earlier studies on the molecular evolution of insect antimicrobial peptide genes have focused on termites and dipteran species, and detected positive selection only in duplicated termicin genes in termites. These findings, together with our present results, provide an indication that the immune systems of social insects (ants and termites) and dipteran insects may have responded differently to the selection pressure caused by microbial pathogens.

Caspar controls resistance to Plasmodium falciparum in diverse anopheline species

Immune responses mounted by the malaria vector Anopheles gambiae are largely regulated by the Toll and Imd (immune deficiency) pathways via the NF-kappaB transcription factors Rel1 and Rel2, which are controlled by the negative regulators Cactus and Caspar, respectively. Rel1- and Rel2-dependent transcription in A. gambiae has been shown to be particularly critical to the mosquito's ability to manage infection with the rodent malaria parasite Plasmodium berghei. Using RNA interference to deplete the negative regulators of these pathways, we found that Rel2 controls resistance of A. gambiae to the human malaria parasite Plasmodium falciparum, whereas Rel 1 activation reduced infection levels. The universal relevance of this defense system across Anopheles species was established by showing that caspar silencing also prevents the development of P. falciparum in the major malaria vectors of Asia and South America, A. stephensi and A. albimanus, respectively. Parallel studies suggest that while Imd pathway activation is most effective against P. falciparum, the Toll pathway is most efficient against P. berghei, highlighting a significant discrepancy between the human pathogen and its rodent model. High throughput gene expression analyses identified a plethora of genes regulated by the activation of the two Rel factors and revealed that the Toll pathway played a more diverse role in mosquito biology than the Imd pathway, which was more immunity-specific. Further analyses of key anti-Plasmodium factors suggest they may be responsible for the Imd pathway-mediated resistance phenotype. Additionally, we found that the fitness cost caused by Rel2 activation through caspar gene silencing was undetectable in sugar-fed, blood-fed, and P. falciparum-infected female A. gambiae, while activation of the Toll pathway's Rel1 had a major impact. This study describes for the first time a single gene that influences an immune mechanism that is able to abort development of P. falciparum in Anopheline species. Further, this study addresses aspects of the molecular, evolutionary, and physiological consequences of the observed phenotype. These findings have implications for malaria control since broad-spectrum immune activation in diverse anopheline species offers a viable and strategic approach to develop novel malaria control methods worldwide.

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.

Salivary gland transcriptome analysis during Plasmodium infection in malaria vector Anopheles stephensi

BACKGROUND

Understanding the tissue-specific molecular cross-talk mechanism during the mosquito-parasite interaction is of prime importance in the design of new strategies for malaria control. Because mosquito salivary glands are the final destination for the parasite maturation and transmission of vector-borne diseases, identification and characterization of salivary genes and their products are equally important in order to access their effect on the infectivity of the parasite. During the last five years there have been several studies on the sialomes of Anopheles mosquitoes, however very limited information is available on the changes in the salivary gland transcriptome in the presence of Plasmodium, and this information is limited to the mosquito Anopheles gambiae.

METHODS

In this study we aimed to explore and identify parasite-induced transcripts from the salivary glands of Anopheles stephensi, using a subtractive hybridization protocol.

RESULTS

Ninety-four percent of expressed sequence tags (ESTs) showed close homology to previously known families of mosquito salivary gland secretary proteins, representing the induced expression of alternative splicing and/or additional new members of the protein family. The remaining 6% of ESTs did not yield significant homology to any known proteins in the non-redundant database and thus may represent a class of unknown/novel salivary proteins. Primary analysis of the ESTs also revealed identification of several novel immune-related transcripts, including defensin and cecropins, probably involved in counter-activation of the antagonistic defense system. A comprehensive description of each family of proteins has been discussed in relation to the tissue-specific mosquito-parasite interaction.

CONCLUSION

This is the first report on the identification of new putative salivary genes, presumably activated during parasite infection.

The evolution of TEP1, an exceptionally polymorphic immunity gene in Anopheles gambiae

Background

Host-parasite coevolution can result in balancing selection, which maintains genetic variation in the susceptibility of hosts to parasites. It has been suggested that variation in a thioester-containing protein called TEP1 (AGAP010815) may alter the ability of Anopheles mosquitoes to transmit Plasmodium parasites, and high divergence between alleles of this gene suggests the possible action of long-term balancing selection. We studied whether TEP1 is a case of an ancient balanced polymorphism in an animal immune system.

Results

We found evidence that the high divergence between TEP1 alleles is the product of genetic exchange between TEP1 and other TEP loci, i.e. gene conversion. Additionally, some TEP1 alleles showed unexpectedly low variability.

Conclusion

The TEP1 gene appears to be a chimera produced from at least two other TEP loci, and the divergence between TEP1 alleles is probably not caused by long-term balancing selection, but is instead due to two independent gene conversion events from one of these other genes. Nevertheless, TEP1 still shows evidence of natural selection, in particular there appears to have been recent changes in the frequency of alleles that has diminished polymorphism within each allelic class. Although the selective force driving this dynamic was not identified, given that susceptibility to Plasmodium parasites is known to be associated with allelic variation in TEP1, these changes in allele frequencies could alter the vectoring capacity of populations.

Population analysis using the nuclear white gene detects Pliocene/Pleistocene lineage divergence within Anopheles nuneztovari in South America

Anopheles (Nyssorhynchus) nuneztovari Gabaldón (Diptera: Culicidae), a locally important malaria vector in some regions of South America, has been hypothesized to consist of at least two cryptic incipient species. We investigated its phylogeographic structure in several South American localities to determine the number of lineages and levels of divergence using the nuclear white gene, a marker that detected two recently diverged genotypes in the primary Neotropical malaria vector Anopheles darlingi Root. In An. nuneztovari, five distinct lineages (1-5) were elucidated: (1) populations from northeastern and central Amazonia; (2) populations from Venezuela east and west of the Andes; (3) populations from Colombia and Venezuela west of the Andes; (4) southeastern and western Amazonian Brazil populations, and (5) southeastern and western Amazonian Brazil and Bolivian populations. There was a large amount of genetic differentiation among these lineages. The deepest and earliest divergence was found between lineage 3 and lineages 1, 2 and 4, which probably accounts for the detection of lineage 3 in some earlier studies. The multiple lineages within Amazonia are partially congruent with previous mtDNA and ITS2 data, but were undetected in many earlier studies, probably because of their recent (Pleistocene) divergence and the differential mutation rates of the markers. The estimates for the five lineages, interpreted as recently evolved or incipient species, date to the Pleistocene and Pliocene. We hypothesize that the diversification in An. nuneztovari is the result of an interaction between the Miocene/Pliocene marine incursion and Pleistocene climatic changes leading to refugial isolation. The identification of cryptic lineages in An. nuneztovari could have a significant impact on local vector control measures.

SNP discovery and molecular evolution in Anopheles gambiae, with special emphasis on innate immune system

Background

Anopheles innate immunity affects Plasmodium development and is a potential target of innovative malaria control strategies. The extent and distribution of nucleotide diversity in immunity genes might provide insights into the evolutionary forces that condition pathogen-vector interactions. The discovery of polymorphisms is an essential step towards association studies of susceptibility to infection.

Results

We sequenced coding fragments of 72 immune related genes in natural populations of Anopheles gambiae and of 37 randomly chosen genes to provide a background measure of genetic diversity across the genome. Mean nucleotide diversity (π) was 0.0092 in the A. gambiae S form, 0.0076 in the M form and 0.0064 in A. arabiensis. Within each species, no statistically significant differences in mean nucleotide diversity were detected between immune related and non immune related genes. Strong purifying selection was detected in genes of both categories, presumably reflecting strong functional constraints.

Conclusion

Our results suggest similar patterns and rates of molecular evolution in immune and non-immune genes in A. gambiae. The 3,214 Single Nucleotide Polymorphisms (SNPs) that we identified are the first large set of Anopheles SNPs from fresh, field-collected material and are relevant markers for future phenotype-association studies.

Molecular and phylogenetic analysis of a novel salivary defensin cDNA from malaria vector Anopheles stephensi

Manipulating the endogenous immune responses of the mosquito such as temporal and spatial expression of antimicrobial peptides may help in the development of a refractory mosquito, unable to transmit malaria. In mosquito several small antimicrobial peptides are activated locally in the midgut and salivary glands upon Plasmodium infection. Anopheles stephensi, the major urban malaria vector in India, has been considered as an important insect model to study vector-parasite interactions; however, so far no reports are available on the antimicrobial peptides from this mosquito species. In the present study, we report identification and molecular characterization of a novel cDNA encoding defensin like peptide, isolated from the salivary gland subtractive hybridization cDNA library of mosquito A. stephensi. Defensin cDNA is 396 base pair long, bearing an open reading frame of 96 amino acids. Deduced amino acid sequence of A. stephensi defensin (Astp_def) contains a signal peptide sequence of 24 amino acids followed by 32-amino acids long putative propeptide domain and a 40-amino acid mature peptide domain carrying 23-amino acid long consensuses sequence signature of insect defensin. Mature peptide of Astp_def carries six conserved cysteine residues, with a predicted molecular weight of 4.20kDa, and isoelectric point of 8.30, characteristic features of cationic defensins. Amino acid sequence similarity and phylogenetic analysis indicated a higher variation in the pre-propeptide region, as compared to the mature defensin peptide, assuring the presence of finely tuned immune responses to counter pathogens.

The molecular evolution of four anti-malarial immune genes in the Anopheles gambiae species complex

Background

If the insect innate immune system is to be used as a potential blocking step in transmission of malaria, then it will require targeting one or a few genes with highest relevance and ease of manipulation. The problem is to identify and manipulate those of most importance to malaria infection without the risk of decreasing the mosquito's ability to stave off infections by microbes in general. Molecular evolution methodologies and concepts can help identify such genes. Within the setting of a comparative molecular population genetic and phylogenetic framework, involving six species of the Anopheles gambiae complex, we investigated whether a set of four pre-selected immunity genes (gambicin, NOS, Rel2 and FBN9) might have evolved under selection pressure imposed by the malaria parasite.

Results

We document varying levels of polymorphism within and divergence between the species, in all four genes. Introgression and the sharing of ancestral polymorphisms, two processes that have been documented in the past, were verified in this study in all four studied genes. These processes appear to affect each gene in different ways and to different degrees. However, there is no evidence of positive selection acting on these genes.

Conclusion

Considering the results presented here in concert with previous studies, genes that interact directly with the Plasmodium parasite, and play little or no role in defense against other microbes, are probably the most likely candidates for a specific adaptive response against P. falciparum. Furthermore, since it is hard to establish direct evidence linking the adaptation of any candidate gene to P. falciparum infection, a comparative framework allowing at least an indirect link should be provided. Such a framework could be achieved, if a similar approach like the one involved here, was applied to all other anopheline complexes that transmit P. falciparum malaria.

Patterns of selection in anti-malarial immune genes in malaria vectors: evidence for adaptive evolution in LRIM1 in Anopheles arabiensis

Background

Co-evolution between Plasmodium species and its vectors may result in adaptive changes in genes that are crucial components of the vector's defense against the pathogen. By analyzing which genes show evidence of positive selection in malaria vectors, but not in closely related non-vectors, we can identify genes that are crucial for the mosquito's resistance against Plasmodium.

Methodology/Principle Findings

We investigated genetic variation of three anti-malarial genes; CEC1, GNBP-B1 and LRIM1, in both vector and non-vector species of the Anopheles gambiae complex. Whereas little protein differentiation was observed between species in CEC1 and GNBP-B1, McDonald-Kreitman and maximum likelihood tests of positive selection show that LRIM1 underwent adaptive evolution in a primary malaria vector; An. arabiensis. In particular, two adjacent codons show clear signs of adaptation by having accumulated three out of four replacement substitutions. Furthermore, our data indicate that this LRIM1 allele has introgressed from An. arabiensis into the other main malaria vector An. gambiae.

Conclusions/Significance

Although no evidence exists to link the adaptation of LRIM1 to P. falciparum infection, an adaptive response of a known anti-malarial gene in a primary malaria vector is intriguing, and may suggest that this gene could play a role in Plasmodium resistance in An. arabiensis. If so, our data also predicts that LRIM1 alleles in An. gambiae vary in their level of resistance against P. falciparum.

Reduced recombination rate and genetic differentiation between the M and S forms of Anopheles gambiae s.s

Genetic differentiation between the largely sympatric molecular forms M and S of Anopheles gambiae appears mostly limited to division 6 and part of division 5 of the X chromosome. This region is adjacent to the centromere and includes the rDNA that was used to define these forms. This localized differentiation between populations that experience gene flow strongly suggests that this region contains genes responsible for reproductive isolation. Regions adjacent to centromeres are known to experience less recombination in several species and it has recently been suggested that low recombination rates can facilitate the accumulation and maintenance of isolation genes in partially isolated populations. Therefore, we measured the recombination rate in division 5D/6 directly and estimate that it is at least 16-fold reduced across this region compared to the remainder of the X chromosome. Additionally, sequence data from four loci from field-collected mosquitoes from several West African countries show very strong differentiation between the molecular forms in division 5D/6, whereas none was observed in two loci elsewhere on the X chromosome. Furthermore, genetic variation was substantially lower in division 5D/6 compared to the two reference loci, and the inferred genealogies of the division 5D/6 genes show patterns consistent with selective sweeps. This suggests that the reduced recombination rate has increased the effect of selection on this region and that our data are consistent with the hypothesis that reduced recombination rates can play a role in the accumulation of isolation genes in the face of gene flow.