Population genetics of Plasmodium resistance genes in Anopheles gambiae: no evidence for strong selection

ITN-VIROINF: Understanding (Harmful) Virus-Host Interactions by Linking Virology and Bioinformatics

Many recent studies highlight the fundamental importance of viruses. Besides their important role as human and animal pathogens, their beneficial, commensal or harmful functions are poorly understood. By developing and applying tailored bioinformatical tools in important virological models, the Marie Skłodowska-Curie Initiative International Training Network VIROINF will provide a better understanding of viruses and the interaction with their hosts. This will open the door to validate methods of improving viral growth, morphogenesis and development, as well as to control strategies against unwanted microorganisms. The key feature of VIROINF is its interdisciplinary nature, which brings together virologists and bioinformaticians to achieve common goals.

Detecting Adaptation with Genome-Scale Molecule Evolutionary Analysis: An Educational Primer for Use with "RNA Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates"

Hosts and pathogens impose coevolutionary pressure on each other as pathogens strive to establish themselves and hosts seek to suppress infection. RNA interference (RNAi) is a mechanism by which cells repress viruses and transposable elements, thereby serving as a form of immune defense. Previous studies have shown that antiviral RNAi genes evolve extraordinarily quickly in the fruit fly Drosophila melanogaster, suggesting that they may adaptively coevolve with viruses and transposable elements. An article by Palmer and colleagues extends this observation to nematodes and multiple insects. Their article can be combined with this Primer to demonstrate the use of comparative genomics and molecular evolutionary analyses in the measurement of natural selection.

RNA-Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates

Conflict between organisms can lead to a reciprocal adaptation that manifests as an increased evolutionary rate in genes mediating the conflict. This adaptive signature has been observed in RNA-interference (RNAi) pathway genes involved in the suppression of viruses and transposable elements in Drosophila melanogaster, suggesting that a subset of Drosophila RNAi genes may be locked in an arms race with these parasites. However, it is not known whether rapid evolution of RNAi genes is a general phenomenon across invertebrates, or which RNAi genes generally evolve adaptively. Here we use population genomic data from eight invertebrate species to infer rates of adaptive sequence evolution, and to test for past and ongoing selective sweeps in RNAi genes. We assess rates of adaptive protein evolution across species using a formal meta-analytic framework to combine data across species and by implementing a multispecies generalized linear mixed model of mutation counts. Across species, we find that RNAi genes display a greater rate of adaptive protein substitution than other genes, and that this is primarily mediated by positive selection acting on the genes most likely to defend against viruses and transposable elements. In contrast, evidence for recent selective sweeps is broadly spread across functional classes of RNAi genes and differs substantially among species. Finally, we identify genes that exhibit elevated adaptive evolution across the analyzed insect species, perhaps due to concurrent parasite-mediated arms races.

The Anopheles gambiae 2La chromosome inversion is associated with susceptibility to Plasmodium falciparum in Africa

Chromosome inversions suppress genetic recombination and establish co-adapted gene complexes, or supergenes. The 2La inversion is a widespread polymorphism in the Anopheles gambiae species complex, the major African mosquito vectors of human malaria. Here we show that alleles of the 2La inversion are associated with natural malaria infection levels in wild-captured vectors from West and East Africa. Mosquitoes carrying the more-susceptible allele (2L+^a) are also behaviorally less likely to be found inside houses. Vector control tools that target indoor-resting mosquitoes, such as bednets and insecticides, are currently the cornerstone of malaria control in Africa. Populations with high levels of the 2L+^a allele may form reservoirs of persistent outdoor malaria transmission requiring novel measures for surveillance and control. The 2La inversion is a major and previously unappreciated component of the natural malaria transmission system in Africa, influencing both malaria susceptibility and vector behavior.

The Anopheles FBN9 immune factor mediates Plasmodium species-specific defense through transgenic fat body expression

Mosquitoes have a multifaceted innate immune system that is actively engaged in warding off various pathogens, including the protozoan malaria parasite Plasmodium. Various immune signaling pathways and effectors have been shown to mediate a certain degree of defense specificity against different Plasmodium species. A key pattern recognition receptor of the Anopheles gambiae immune system is the fibrinogen domain-containing immunolectin FBN9, which has been shown to be transcriptonally induced by Plasmodium infection, and to mediate defense against both rodent and human malaria parasites and bacteria. Here we have further studied the defense specificity of FBN9 using a transgenic approach, in which FBN9 is overexpressed in the fat body tissue after a blood meal through a vitellogenin promoter. Interestingly, the Vg-FBN9 transgenic mosquitoes showed increased resistance only to the rodent parasite P. berghei, and not to the human parasite P. falciparum, pointing to differences in the mosquito's defense mechanisms against the two parasite species. The Vg-FBN9 transgenic mosquitoes were also more resistant to infection with both Gram-positive and Gram-negative bacteria and showed increased longevity when infected with P. berghei. Our study points to the importance of both experimentally depleting and enriching candidate anti-Plasmodium effectors in functional studies in order to ascertain their suitability for the development of transgenic mosquito-based malaria control strategies.

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.

Malaria vector populations across ecological zones in Guinea Conakry and Mali, West Africa

Background

Malaria remains a pervasive public health problem in sub-Saharan West Africa. Here mosquito vector populations were explored across four sites in Mali and the Republic of Guinea (Guinea Conakry). The study samples the major ecological zones of malaria-endemic regions in West Africa within a relatively small distance.

Methods

Mosquito vectors were sampled from larval pools, adult indoor resting sites, and indoor and outdoor human-host seeking adults. Mosquitoes were collected at sites spanning 350 km that represented arid savannah, humid savannah, semi-forest and deep forest ecological zones, in areas where little was previously known about malaria vector populations. 1425 mosquito samples were analysed by molecular assays to determine species, genetic attributes, blood meal sources and Plasmodium infection status.

Results

Anopheles gambiae and Anopheles coluzzii were the major anophelines represented in all collections across the ecological zones, with A. coluzzii predominant in the arid savannah and A. gambiae in the more humid sites. The use of multiple collection methodologies across the sampling sites allows assessment of potential collection bias of the different methods. The L1014F kdr insecticide resistance mutation (kdr-w) is found at high frequency across all study sites. This mutation appears to have swept almost to fixation, from low frequencies 6 years earlier, despite the absence of widespread insecticide use for vector control. Rates of human feeding are very high across ecological zones, with only small fractions of animal derived blood meals in the arid and humid savannah. About 30 % of freshly blood-fed mosquitoes were positive for Plasmodium falciparum presence, while the rate of mosquitoes with established infections was an order of magnitude lower.

Conclusions

The study represents detailed vector characterization from an understudied area in West Africa with endemic malaria transmission. The deep forest study site includes the epicenter of the 2014 Ebola virus epidemic. With new malaria control interventions planned in Guinea, these data provide a baseline measure and an opportunity to assess the outcome of future interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1242-5) contains supplementary material, which is available to authorized users.

A mosquito mediator of parasite-induced immune priming

Immune memory is a central feature of the mammalian adaptive immune system. The more primitive innate immune system of insects has also been shown to comprise memory, or immune priming. A recent study has shed new light on how Plasmodium primes the mosquito immune system for greater resistance to a subsequent infection with the same pathogen.

Population genetics of Anopheles coluzzii immune pathways and genes

Natural selection is expected to drive adaptive evolution in genes involved in host–pathogen interactions. In this study, we use molecular population genetic analyses to understand how natural selection operates on the immune system of Anopheles coluzzii (formerly A. gambiae “M form”). We analyzed patterns of intraspecific and interspecific genetic variation in 20 immune-related genes and 17 nonimmune genes from a wild population of A. coluzzii and asked if patterns of genetic variation in the immune genes are consistent with pathogen-driven selection shaping the evolution of defense. We found evidence of a balanced polymorphism in CTLMA2, which encodes a C-type lectin involved in regulation of the melanization response. The two CTLMA2 haplotypes, which are distinguished by fixed amino acid differences near the predicted peptide cleavage site, are also segregating in the sister species A. gambiae (“S form”) and A. arabiensis. Comparison of the two haplotypes between species indicates that they were not shared among the species through introgression, but rather that they arose before the species divergence and have been adaptively maintained as a balanced polymorphism in all three species. We additionally found that STAT-B, a retroduplicate of STAT-A, shows strong evidence of adaptive evolution that is consistent with neofunctionalization after duplication. In contrast to the striking patterns of adaptive evolution observed in these Anopheles-specific immune genes, we found no evidence of adaptive evolution in the Toll and Imd innate immune pathways that are orthologously conserved throughout insects. Genes encoding the Imd pathway exhibit high rates of amino acid divergence between Anopheles species but also display elevated amino acid diversity that is consistent with relaxed purifying selection. These results indicate that adaptive coevolution between A. coluzzii and its pathogens is more likely to involve novel or lineage-specific molecular mechanisms than the canonical humoral immune pathways.

Genomic analyses of three malaria vectors reveals extensive shared polymorphism but contrasting population histories

Anopheles gambiae s.l. are important malaria vectors, but little is known about their genomic variation in the wild. Here, we present inter- and intraspecies analysis of genome-wide RADseq data, in three Anopheles gambiae s.l. species collected from East Africa. The mosquitoes fall into three genotypic clusters representing described species (A. gambiae, A. arabiensis, and A. merus) with no evidence of cryptic breeding units. Anopheles merus is the most divergent of the three species, supporting a recent new phylogeny based on chromosomal inversions. Even though the species clusters are well separated, there is extensive shared polymorphism, particularly between A. gambiae and A. arabiensis. Divergence between A. gambiae and A. arabiensis does not vary across the autosomes but is higher in X-linked inversions than elsewhere on X or on the autosomes, consistent with the suggestion that this inversion (or a gene within it) is important in reproductive isolation between the species. The 2La/2L+(a) inversion shows no more evidence of introgression between A. gambiae and A. arabiensis than the rest of the autosomes. Population differentiation within A. gambiae and A. arabiensis is weak over approximately 190-270 km, implying no strong barriers to dispersal. Analysis of Tajima's D and the allele frequency spectrum is consistent with modest population increases in A. arabiensis and A. merus, but a more complex demographic history of expansion followed by contraction in A. gambiae. Although they are less than 200 km apart, the two A. gambiae populations show evidence of different demographic histories.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon

Background

Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance.

Methods

A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali.

Results

Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion.

Discussion

More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

Animal mitochondria, positive selection and cyto-nuclear coevolution: insights from pulmonates

Pulmonate snails have remarkably high levels of mtDNA polymorphism within species and divergence between species, making them an interesting group for the study of mutation and selection on mitochondrial genomes. The availability of sequence data from most major lineages – collected largely for studies of phylogeography - provides an opportunity to perform several tests of selection that may provide general insights into the evolutionary forces that have produced this unusual pattern. Several protein coding mtDNA datasets of pulmonates were analyzed towards this direction. Two different methods for the detection of positive selection were used, one based on phylogeny, and the other on the McDonald-Kreitman test. The cyto-nuclear coevolution hypothesis, often implicated to account for the high levels of mtDNA divergence of some organisms, was also addressed by assessing the divergence pattern exhibited by a nuclear gene. The McDonald-Kreitman test indicated multiple signs of positive selection in the mtDNA genes, but was significantly biased when sequence divergence was high. The phylogenetic method identified five mtDNA datasets as affected by positive selection. In the nuclear gene, the McDonald-Kreitman test provided no significant results, whereas the phylogenetic method identified positive selection as likely present. Overall, our findings indicate that: 1) slim support for the cyto-nuclear coevolution hypothesis is present, 2) the elevated rates of mtDNA polymorphims and divergence in pulmonates do not appear to be due to pervasive positive selection, 3) more stringent tests show that spurious positive selection is uncovered when distant taxa are compared and 4) there are significant examples of positive selection acting in some cases, so it appears that mtDNA evolution in pulmonates can escape from strict deleterious evolution suggested by the Muller’s ratchet effect.

No evidence for positive selection at two potential targets for malaria transmission-blocking vaccines in Anopheles gambiae s.s

Human malaria causes nearly a million deaths in sub-Saharan Africa each year. The evolution of drug-resistance in the parasite and insecticide resistance in the mosquito vector has complicated control measures and made the need for new control strategies more urgent. Anopheles gambiae s.s. is one of the primary vectors of human malaria in Africa, and parasite-transmission-blocking vaccines targeting Anopheles proteins have been proposed as a possible strategy to control the spread of the disease. However, the success of these hypothetical technologies would depend on the successful ability to broadly target mosquito populations that may be genetically heterogeneous. Understanding the evolutionary pressures shaping genetic variation among candidate target molecules offers a first step towards evaluating the prospects of successfully deploying such technologies. We studied the population genetics of genes encoding two candidate target proteins, the salivary gland protein saglin and the basal lamina structural protein laminin, in wild populations of the M and S molecular forms of A. gambiae in Mali. Through analysis of intraspecific genetic variation and interspecific comparisons, we found no evidence of positive natural selection at the genes encoding these proteins. On the contrary, we found evidence for particularly strong purifying selection at the laminin gene. These results provide insight into the patterns of genetic diversity of saglin and laminin, and we discuss these findings in relation to the potential development of these molecules as vaccine targets.

Evidence for population-specific positive selection on immune genes of Anopheles gambiae

Host-pathogen interactions can be powerful drivers of adaptive evolution, shaping the patterns of molecular variation at the genes involved. In this study, we sequenced alleles from 28 immune-related loci in wild samples of multiple genetic subpopulations of the African malaria mosquito Anopheles gambiae, obtaining unprecedented sample sizes and providing the first opportunity to contrast patterns of molecular evolution at immune-related loci in the recently discovered GOUNDRY population to those of the indoor-resting M and S molecular forms. In contrast to previous studies that focused on immune genes identified in laboratory studies, we centered our analysis on genes that fall within a quantitative trait locus associated with resistance to Plasmodium falciparum in natural populations of A. gambiae. Analyses of haplotypic and genetic diversity at these 28 loci revealed striking differences among populations in levels of genetic diversity and allele frequencies in coding sequence. Putative signals of positive selection were identified at 11 loci, but only one was shared by two subgroups of A. gambiae. Striking patterns of linkage disequilibrium were observed at several loci. We discuss these results with respect to ecological differences among these strata as well as potential implications for disease transmission.

Single-nucleotide polymorphisms for high-throughput genotyping of Anopheles arabiensis in East and southern Africa

Anopheles arabiensis Patton is one of the principal vectors of malaria in sub-Saharan Africa, occupying a wide variety of ecological zones. This species is increasingly responsible for malaria transmission in Africa and is becoming the dominant vector species in some localities. Despite its growing importance, little is known about genetic polymorphisms in this species. Multiple sequences of various gene fragments from An. arabiensis isolates from Cameroon were obtained from GenBank. In total, 20 gene fragments containing single-nucleotide polymorphisms (SNPs) at moderate density were selected for direct sequencing from field collected specimens from Tanzania and Zambia. We obtained 301 SNPs in total from the 20 gene fragments, 60 of which were suitable for Illumina GoldenGate SNP genotyping. A greater number of SNPs (n = 185) was suitable for analysis using Sequenom iPLEX, an alternative high-throughput genotyping technology using mass spectrometry. An SNP was present every 59 (+/- 44.5) bases on average. Overall, An. arabiensis from Tanzania and Zambia are genetically closer (mean F(ST) = 0.075) than either is to populations in Cameroon (F(ST, TZ-CM) = 0.250, F(ST,ZA-CM) = 0.372). A fixed polymorphism between East/southern and Central Africa was identified on AGAP000574, a gene on the X chromosome. We have identified SNPs in natural populations of An. arabiensis. SNP densities in An. arabiensis were higher than Anopheles gambiae s.s., suggesting a greater challenge in the development of high-throughput SNP analysis for this species. The SNP markers provided in this study are suitable for a high-throughput genotyping analysis and can be used for population genetic studies and association mapping efforts.

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.

The demographic histories of the M and S molecular forms of Anopheles gambiae s.s

Anopheles gambiae is a primary vector of Plasmodium falciparum, a human malaria parasite that causes over a million deaths each year in sub-Saharan Africa. Population genetic tests have been employed to detect natural selection at suspected A. gambiae antimalaria genes, but these tests have generally been compromised by the lack of demographically correct null models. Here, we used a coalescent simulation approach within a maximum likelihood framework to fit population growth, bottleneck, and migration models to polymorphism data from Cameroonian A. gambiae. The best-fit models for both the "M" and the "S" molecular forms of A. gambiae included ancient population growth and a high rate of migration from an unsampled subpopulation. After correcting for differences in effective population size, our models suggest that the molecular forms expanded at different times and both expansions significantly predate the advent of agriculture. We show that correcting null models for demography increases the power to detect natural selection in A. gambiae.

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.

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.

Rodent malaria-resistant strains of the mosquito, Anopheles gambiae, have slower population growth than -susceptible strains

Background

Trade-offs between anti-parasite defence mechanisms and other life history traits limit the evolution of host resistance to parasites and have important implications for understanding diseases such as malaria. Mosquitoes have not evolved complete resistance to malaria parasites and one hypothesis is that anti-malaria defence mechanisms are costly.

Results

We used matrix population models to compare the population growth rates among lines of Anopheles gambiae that had been selected for resistance or high susceptibility to the rodent malaria parasite, Plasmodium yoelii nigeriensis. The population growth rate of the resistant line was significantly lower than that of the highly susceptible and the unselected control lines, regardless of whether mosquitoes were infected with Plasmodium or not. The lower population growth of malaria-resistant mosquitoes was caused by reduced post blood-feeding survival of females and poor egg hatching.

Conclusion

With respect to eradicating malaria, the strategy of releasing Plasmodium-resistant Anopheles mosquitoes is unlikely to be successful if the costs of Plasmodium-resistance in the field are as great as the ones measured in this study. High densities of malaria-resistant mosquitoes would have to be maintained by continuous release from captive breeding facilities.

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.

Clarification of anomalies in the application of a 2La molecular karyotyping method for the malaria vector Anopheles gambiae

Background

Chromosomal inversions have been considered to be potentially important barriers to gene flow in many groups of animals through their effect on recombination suppression in heterokaryotypic individuals. Inversions can also enhance local adaptation in different groups of organisms and may often represent species-specific differences among closely related taxa. We conducted a study to characterize the 2La inversion karyotypes of An. gambiae sensu stricto mosquitoes sampled from the Kilombero Valley (Tanzania) using a newly designed PCR assay.

Results

We frequently encountered a (687 bp) fragment which was only present in the Kilombero Valley populations. Laboratory crossing between An. gambiae s.s. from Njage (Tanzania) and Kisumu (Western Kenya) populations resulted in F₁ offspring carrying the observed fragment. Karyotype analysis did not indicate differences in 2La region chromosome morphology between individuals carrying the PCR fragments, the 207 bp fragment, or the 687 bp fragement.

Conclusion

The observed insertion/deletion polymorphism within the region amplified by the 2La PCR diagnostic test may confound the interpretation of this assay and should be well considered in order to maintain an acceptable level of reliability in studies using this assay to describe the distribution and frequency of the 2La inversion among natural populations of An. gambiae s.s.

Anopheles mosquitoes: not just flying malaria vectors... especially in the field

The polymorphism of genes involved in the immunity of malaria vectors has been the subject of several recent studies with mosquitoes from natural populations. Most of the genes examined are known for their role against Plasmodium berghei and not necessarily for their role against Plasmodium falciparum. It seems, therefore, to be highly important not only to be cautious when linking natural selection with malaria epidemiology but also to consider the importance of other parasites and the environment on the mosquito genome.

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.

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.

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.

Adaptive evolution of a stress response protein

BACKGROUND

Some cancers are mediated by an interplay between tissue damage, pathogens and localised innate immune responses, but the mechanisms that underlie these linkages are only beginning to be unravelled.

METHODS AND PRINCIPAL FINDINGS

Here we identify a strong signature of adaptive evolution on the DNA sequence of the mammalian stress response gene SEP53, a member of the epidermal differentiation complex fused-gene family known for its role in suppressing cancers. The SEP53 gene appears to have been subject to adaptive evolution of a type that is commonly (though not exclusively) associated with coevolutionary arms races. A similar pattern of molecular evolution was not evident in the p53 cancer-suppressing gene.

CONCLUSIONS

Our data thus raises the possibility that SEP53 is a component of the mucosal/epithelial innate immune response engaged in an ongoing interaction with a pathogen. Although the pathogenic stress mediating adaptive evolution of SEP53 is not known, there are a number of well-known candidates, in particular viruses with established links to carcinoma.