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

Transcriptional variation of sensory-related genes in natural populations of Aedes albopictus

BACKGROUND

The Asian tiger mosquito, Aedes albopictus, is a highly dangerous invasive vector of numerous medically important arboviruses including dengue, chikungunya and Zika. In four decades it has spread from tropical Southeast Asia to many parts of the world in both tropical and temperate climes. The rapid invasion process of this mosquito is supported by its high ecological and genetic plasticity across different life history traits. Our aim was to investigate whether wild populations, both native and adventive, also display transcriptional genetic variability for functions that may impact their biology, behaviour and ability to transmit arboviruses, such as sensory perception.

RESULTS

Antennal transcriptome data were derived from mosquitoes from a native population from Ban Rai, Thailand and from three adventive Mediterranean populations: Athens, Greece and Arco and Trento from Italy. Clear inter-population differential transcriptional activity was observed in different gene categories related to sound perception, olfaction and viral infection. The greatest differences were detected between the native Thai and the Mediterranean populations. The two Italian populations were the most similar. Nearly one million quality filtered SNP loci were identified.

CONCLUSION

The ability to express this great inter-population transcriptional variability highlights, at the functional level, the remarkable genetic flexibility of this mosquito species. We can hypothesize that the differential expression of genes, including those involved in sensory perception, in different populations may enable Ae. albopictus to exploit different environments and hosts, thus contributing to its status as a global vector of arboviruses of public health importance. The large number of SNP loci present in these transcripts represents a useful addition to the arsenal of high-resolution molecular markers and a resource that can be used to detect selective pressure and adaptive changes that may have occurred during the colonization process.

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole-genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination-reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non-Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour-joining trees based on whole-genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.

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.

Genome-Wide Divergence in the West-African Malaria Vector Anopheles melas

Anopheles melas is a member of the recently diverged An. gambiae species complex, a model for speciation studies, and is a locally important malaria vector along the West-African coast where it breeds in brackish water. A recent population genetic study of An. melas revealed species-level genetic differentiation between three population clusters. An. melas West extends from The Gambia to the village of Tiko, Cameroon. The other mainland cluster, An. melas South, extends from the southern Cameroonian village of Ipono to Angola. Bioko Island, Equatorial Guinea An. melas populations are genetically isolated from mainland populations. To examine how genetic differentiation between these An. melas forms is distributed across their genomes, we conducted a genome-wide analysis of genetic differentiation and selection using whole genome sequencing data of pooled individuals (Pool-seq) from a representative population of each cluster. The An. melas forms exhibit high levels of genetic differentiation throughout their genomes, including the presence of numerous fixed differences between clusters. Although the level of divergence between the clusters is on a par with that of other species within the An. gambiae complex, patterns of genome-wide divergence and diversity do not provide evidence for the presence of pre- and/or postmating isolating mechanisms in the form of speciation islands. These results are consistent with an allopatric divergence process with little or no introgression.

Patterns of Genome-Wide Variation in Glossina fuscipes fuscipes Tsetse Flies from Uganda

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r(2) max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r(2) max/2 at 708 bp, an order of magnitude slower than in Drosophila Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Different evolutionary patterns of SNPs between domains and unassigned regions in human protein-coding sequences

Protein evolution plays an important role in the evolution of each genome. Because of their functional nature, in general, most of their parts or sites are differently constrained selectively, particularly by purifying selection. Most previous studies on protein evolution considered individual proteins in their entirety or compared protein-coding sequences with non-coding sequences. Less attention has been paid to the evolution of different parts within each protein of a given genome. To this end, based on PfamA annotation of all human proteins, each protein sequence can be split into two parts: domains or unassigned regions. Using this rationale, single nucleotide polymorphisms (SNPs) in protein-coding sequences from the 1000 Genomes Project were mapped according to two classifications: SNPs occurring within protein domains and those within unassigned regions. With these classifications, we found: the density of synonymous SNPs within domains is significantly greater than that of synonymous SNPs within unassigned regions; however, the density of non-synonymous SNPs shows the opposite pattern. We also found there are signatures of purifying selection on both the domain and unassigned regions. Furthermore, the selective strength on domains is significantly greater than that on unassigned regions. In addition, among all of the human protein sequences, there are 117 PfamA domains in which no SNPs are found. Our results highlight an important aspect of protein domains and may contribute to our understanding of protein evolution.

Genetic Structure of a Local Population of the Anopheles gambiae Complex in Burkina Faso

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).

An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Anopheles gambiae and Anopheles coluzzii are the primary mosquito vectors of human malaria in West Africa. Both of these closely related species efficiently transmit the disease, although they display ecological differences. Previous work showed that A. coluzzii displays distinct genetic patterns in genes important for mosquito immunity. Here, we use this genetic pattern as a filter to examine a panel of potential immune genes, and show that the genetic pattern is strongly predictive for genes that play a functional role in immunity when tested with malaria parasites.

Ancient DNA reveals that the genetic structure of the northern Han Chinese was shaped prior to 3,000 years ago

The Han Chinese are the largest ethnic group in the world, and their origins, development, and expansion are complex. Many genetic studies have shown that Han Chinese can be divided into two distinct groups: northern Han Chinese and southern Han Chinese. The genetic history of the southern Han Chinese has been well studied. However, the genetic history of the northern Han Chinese is still obscure. In order to gain insight into the genetic history of the northern Han Chinese, 89 human remains were sampled from the Hengbei site which is located in the Central Plain and dates back to a key transitional period during the rise of the Han Chinese (approximately 3,000 years ago). We used 64 authentic mtDNA data obtained in this study, 27 Y chromosome SNP data profiles from previously studied Hengbei samples, and genetic datasets of the current Chinese populations and two ancient northern Chinese populations to analyze the relationship between the ancient people of Hengbei and present-day northern Han Chinese. We used a wide range of population genetic analyses, including principal component analyses, shared mtDNA haplotype analyses, and geographic mapping of maternal genetic distances. The results show that the ancient people of Hengbei bore a strong genetic resemblance to present-day northern Han Chinese and were genetically distinct from other present-day Chinese populations and two ancient populations. These findings suggest that the genetic structure of northern Han Chinese was already shaped 3,000 years ago in the Central Plain area.

Detecting SNP combinations discriminating human populations from HapMap data

The genomes of different human beings are similar. There are only a relatively small number of genetic differences between people. The genetic differences between people are very worthy of study. Researchers have proposed the fixation index FST measurement to find the single nucleotide polymorphisms (SNPs) which can reflect human population differences. However, most SNPs have interactions and they work together, which leads to the differences among human populations. The number of all possible m-locus combinations chosen from n SNPs grows exponentially. Most methods concern on 2-locus interactions. In this paper, we propose a novel method to find a new coordinate system under which the energy distributions of different populations are quite different. We select out candidate SNPs from n SNPs by using the information of the axes in the coordinate system. The number of candidate SNPs is small, thus SNP-SNP interactions can be searched efficiently. The method can also find interactions of more than two loci. These interactions should be able to reflect the evolution of human populations from another way. The numbers of SNP-SNP interactions are regarded as the differences between pairwise populations and a hierarchical clustering algorithm is used to construct the evolutionary tree. In the experiments, we apply the method to SNP data of four chromosomes separately and the trees constructed on these four chromosomes are highly consistent. Furthermore, the trees are also consistent with previous studies, which indicates that evolutionary information is well mined. The method provides a new insight to analyze the human population differences.

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.

No more non-model species: the promise of next generation sequencing for comparative immunology

Next generation sequencing (NGS) allows for the rapid, comprehensive and cost effective analysis of entire genomes and transcriptomes. NGS provides approaches for immune response gene discovery, profiling gene expression over the course of parasitosis, studying mechanisms of diversification of immune receptors and investigating the role of epigenetic mechanisms in regulating immune gene expression and/or diversification. NGS will allow meaningful comparisons to be made between organisms from different taxa in an effort to understand the selection of diverse strategies for host defence under different environmental pathogen pressures. At the same time, it will reveal the shared and unique components of the immunological toolkit and basic functional aspects that are essential for immune defence throughout the living world. In this review, we argue that NGS will revolutionize our understanding of immune responses throughout the animal kingdom because the depth of information it provides will circumvent the need to concentrate on a few "model" species.

Stress dependent infection cost of the human malaria agent Plasmodium falciparum on its natural vector Anopheles coluzzii

Unraveling selective forces that shape vector-parasite interactions has critical implications for malaria control. However, it remains unclear whether Plasmodium infection induces a fitness cost to their natural mosquito vectors. Moreover, environmental conditions are known to affect infection outcome and may impact the effect of infection on mosquito fitness. We investigated in the laboratory the effects of exposition to and infection by field isolates of Plasmodium falciparum on fecundity and survival of a major vector in the field, Anopheles coluzzii under different conditions of access to sugar resources after blood feeding. The results evidenced fitness costs induced by exposition and infection. When sugar was available after blood meal, infected and exposed mosquitoes had either reduced or equal to survival to unexposed mosquitoes while fecundity was either increased or decreased depending on the blood donor. Under strong nutritional stress, survival was reduced for exposed and infected mosquitoes in all assays. We therefore provide here evidence of an environmental-dependant reduced survival in mosquitoes exposed to infection in a natural and one of the most important parasite-mosquito species associations for human malaria transmission.

Positive selection drives accelerated evolution of mosquito salivary genes associated with blood-feeding

The saliva of bloodsucking animals contains dozens to hundreds of proteins that counteract their hosts' haemostasis, inflammation and immunity. It was previously observed that salivary proteins involved in haematophagy are much more divergent in their primary sequence than those of housekeeping function, when comparisons were made between closely related organisms. While this pattern of evolution could result from relaxed selection or drift, it could alternatively be the result of positive selection driven by the intense pressure of the host immune system. We investigated the polymorphism of five different genes associated with blood-feeding in the mosquito Anopheles gambiae and obtained evidence in four genes for sites with signatures of positive selection. These results add salivary gland genes from bloodsucking arthropods to the small list of genes driven by positive selection.

Diversity, differentiation, and linkage disequilibrium: prospects for association mapping in the malaria vector Anopheles arabiensis

Association mapping is a widely applied method for elucidating the genetic basis of phenotypic traits. However, factors such as linkage disequilibrium and levels of genetic diversity influence the power and resolution of this approach. Moreover, the presence of population subdivision among samples can result in spurious associations if not accounted for. As such, it is useful to have a detailed understanding of these factors before conducting association mapping experiments. Here we conducted whole-genome sequencing on 24 specimens of the malaria mosquito vector, Anopheles arabiensis, to further understanding of patterns of genetic diversity, population subdivision and linkage disequilibrium in this species. We found high levels of genetic diversity within the An. arabiensis genome, with ~800,000 high-confidence, single- nucleotide polymorphisms detected. However, levels of nucleotide diversity varied significantly both within and between chromosomes. We observed lower diversity on the X chromosome, within some inversions, and near centromeres. Population structure was absent at the local scale (Kilombero Valley, Tanzania) but detected between distant populations (Cameroon vs. Tanzania) where differentiation was largely restricted to certain autosomal chromosomal inversions such as 2Rb. Overall, linkage disequilibrium within An. arabiensis decayed very rapidly (within 200 bp) across all chromosomes. However, elevated linkage disequilibrium was observed within some inversions, suggesting that recombination is reduced in those regions. The overall low levels of linkage disequilibrium suggests that association studies in this taxon will be very challenging for all but variants of large effect, and will require large sample sizes.

An analysis of two island groups as potential sites for trials of transgenic mosquitoes for malaria control

Considerable technological advances have been made towards the generation of genetically modified mosquitoes for vector control. In contrast, less progress has been made towards field evaluations of transformed mosquitoes which are critical for evaluating the success of, and hazards associated with, genetic modification. Oceanic islands have been highlighted as potentially the best locations for such trials. However, population genetic studies are necessary to verify isolation. Here, we used a panel of genetic markers to assess for evidence of genetic isolation of two oceanic island populations of the African malaria vector, Anopheles gambiae s.s. We found no evidence of isolation between the Bijagós archipelago and mainland Guinea-Bissau, despite separation by distances beyond the known dispersal capabilities of this taxon. Conversely, the Comoros Islands appear to be genetically isolated from the East African mainland, and thus represent a location worthy of further investigation for field trials. Based on assessments of gene flow within and between the Comoros islands, the island of Grande Comore was found to be genetically isolated from adjacent islands and also exhibited local population structure, indicating that it may be the most suitable site for trials with existing genetic modification technologies.

Accelerated evolution of innate immunity proteins in social insects: adaptive evolution or relaxed constraint?

The genomes of eusocial insects have a reduced complement of immune genes-an unusual finding considering that sociality provides ideal conditions for disease transmission. The following three hypotheses have been invoked to explain this finding: 1) social insects are attacked by fewer pathogens, 2) social insects have effective behavioral or 3) novel molecular mechanisms for combating pathogens. At the molecular level, these hypotheses predict that canonical innate immune pathways experience a relaxation of selective constraint. A recent study of several innate immune genes in ants and bees showed a pattern of accelerated amino acid evolution, which is consistent with either positive selection or a relaxation of constraint. We studied the population genetics of innate immune genes in the honey bee Apis mellifera by partially sequencing 13 genes from the bee's Toll pathway (∼10.5 kb) and 20 randomly chosen genes (∼16.5 kb) sequenced in 43 diploid workers. Relative to the random gene set, Toll pathway genes had significantly higher levels of amino acid replacement mutations segregating within A. mellifera and fixed between A. mellifera and A. cerana. However, levels of diversity and divergence at synonymous sites did not differ between the two gene sets. Although we detect strong signs of balancing selection on the pathogen recognition gene pgrp-sa, many of the genes in the Toll pathway show signatures of relaxed selective constraint. These results are consistent with the reduced complement of innate immune genes found in social insects and support the hypothesis that some aspect of eusociality renders canonical innate immunity superfluous.

Snpdat: easy and rapid annotation of results from de novo snp discovery projects for model and non-model organisms

Background

Single nucleotide polymorphisms (SNPs) are the most abundant genetic variant found in vertebrates and invertebrates. SNP discovery has become a highly automated, robust and relatively inexpensive process allowing the identification of many thousands of mutations for model and non-model organisms. Annotating large numbers of SNPs can be a difficult and complex process. Many tools available are optimised for use with organisms densely sampled for SNPs, such as humans. There are currently few tools available that are species non-specific or support non-model organism data.

Results

Here we present SNPdat, a high throughput analysis tool that can provide a comprehensive annotation of both novel and known SNPs for any organism with a draft sequence and annotation. Using a dataset of 4,566 SNPs identified in cattle using high-throughput DNA sequencing we demonstrate the annotations performed and the statistics that can be generated by SNPdat.

Conclusions

SNPdat provides users with a simple tool for annotation of genomes that are either not supported by other tools or have a small number of annotated SNPs available. SNPdat can also be used to analyse datasets from organisms which are densely sampled for SNPs. As a command line tool it can easily be incorporated into existing SNP discovery pipelines and fills a niche for analyses involving non-model organisms that are not supported by many available SNP annotation tools. SNPdat will be of great interest to scientists involved in SNP discovery and analysis projects, particularly those with limited bioinformatics experience.

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.

Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth

The fight against diseases spread by mosquitoes and other insects has enormous environmental, economic and social consequences. Chemical insecticides remain the first line of defence but the control of diseases, especially malaria and dengue fever, is being increasingly undermined by insecticide resistance. Mosquitoes have a large repertoire of P450s (over 100 genes). By pinpointing the key enzymes associated with insecticide resistance we can begin to develop new tools to aid the implementation of control interventions and reduce their environmental impact on Earth. Recent technological advances are helping us to build a functional profile of the P450 determinants of insecticide metabolic resistance in mosquitoes. Alongside, the cross-responses of mosquito P450s to insecticides and pollutants are also being investigated. Such research will provide the means to produce diagnostic tools for early detection of P450s linked to resistance. It will also enable the design of new insecticides with optimized efficacy in different environments.

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.

The Antimicrobial Defense of the Pacific Oyster, Crassostrea gigas. How Diversity may Compensate for Scarcity in the Regulation of Resident/Pathogenic Microflora

Healthy oysters are inhabited by abundant microbial communities that vary with environmental conditions and coexist with immunocompetent cells in the circulatory system. In Crassostrea gigas oysters, the antimicrobial response, which is believed to control pathogens and commensals, relies on potent oxygen-dependent reactions and on antimicrobial peptides/proteins (AMPs) produced at low concentrations by epithelial cells and/or circulating hemocytes. In non-diseased oysters, hemocytes express basal levels of defensins (Cg-Defs) and proline-rich peptides (Cg-Prps). When the bacterial load dramatically increases in oyster tissues, both AMP families are driven to sites of infection by major hemocyte movements, together with bactericidal permeability/increasing proteins (Cg-BPIs) and given forms of big defensins (Cg-BigDef), whose expression in hemocytes is induced by infection. Co-localization of AMPs at sites of infection could be determinant in limiting invasion as synergies take place between peptide families, a phenomenon which is potentiated by the considerable diversity of AMP sequences. Besides, diversity occurs at the level of oyster AMP mechanisms of action, which range from membrane lysis for Cg-BPI to inhibition of metabolic pathways for Cg-Defs. The combination of such different mechanisms of action may account for the synergistic activities observed and compensate for the low peptide concentrations in C. gigas cells and tissues. To overcome the oyster antimicrobial response, oyster pathogens have developed subtle mechanisms of resistance and evasion. Thus, some Vibrio strains pathogenic for oysters are equipped with AMP-sensing systems that trigger resistance. More generally, the known oyster pathogenic vibrios have evolved strategies to evade intracellular killing through phagocytosis and the associated oxidative burst.

Cloning and characterization of prophenoloxidase A3 (proPOA3) from Culex pipiens pallens

The prophenoloxidase subunit A3 (proPOA3) gene was cloned from Culex pipiens pallens, which had an open reading frame of 2061 bp encoding a putative 686 amino acid protein. The deduced amino acid sequence shares 98% with proPOA3 from Culex quinquefasciatus. ProPOA3 is expressed at all developmental stages of C. pipiens pallens. Significant negative correlation was observed between proPOA3 expression and deltamethrin resistance in resistant C. pipiens pallens. Furthermore, proPOA3 expression levels were significantly lower in deltamethrin-resistant mosquitoes than in susceptible mosquitoes collected at four locations in Eastern China. However, we did not find any substantial change in proPOA3 expression in field-collected resistant Anopheles mosquitoes. Moreover, overexpressing proPOA3 in C6/36 cells led to more sensitivity to deltamethrin treatment. In laboratory and field-collected resistant C. pipiens pallens, a valine to isoleucine mutation (769G>A) and two synonymous mutations (1116G>C and 1116G>A) were identified in proPOA3. In addition, the mutation frequency of 769G>A and 1116G>C increased gradually, which corresponded with raised deltamethrin resistance levels. Taken together, our study provides the first evidence that proPOA3 may play a role in the regulation of deltamethrin-resistance in C. pipiens pallens.

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.

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.

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.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiaemolecular form identification

BACKGROUND

Anopheles gambiae M and S molecular forms, the major malaria vectors in the Afro-tropical region, are ongoing a process of ecological diversification and adaptive lineage splitting, which is affecting malaria transmission and vector control strategies in West Africa. These two incipient species are defined on the basis of single nucleotide differences in the IGS and ITS regions of multicopy rDNA located on the X-chromosome. A number of PCR and PCR-RFLP approaches based on form-specific SNPs in the IGS region are used for M and S identification. Moreover, a PCR-method to detect the M-specific insertion of a short interspersed transposable element (SINE200) has recently been introduced as an alternative identification approach. However, a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification was never carried out to evaluate whether they could be used interchangeably, as commonly assumed.

RESULTS

The genotyping of more than 400 A. gambiae specimens from nine African countries, and the sequencing of the IGS-amplicon of 115 of them, highlighted discrepancies among results obtained by the different approaches due to different kinds of biases, which may result in an overestimation of MS putative hybrids, as follows: i) incorrect match of M and S specific primers used in the allele specific-PCR approach; ii) presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches; iii) incomplete cleavage during the restriction reactions; iv) presence of different copy numbers of M and S-specific IGS-arrays in single individuals in areas of secondary contact between the two forms.

CONCLUSIONS

The results reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae M and S forms are not fully interchangeable as usually assumed, and highlight limits of the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. These limits are discussed and operational suggestions on the choice of the most convenient method for large-scale M- and S-form identification are provided, also taking into consideration technical aspects related to the epidemiological characteristics of different study areas.

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.

Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature--where larvae may be the more vulnerable life stage--the process of adaptive divergence between M and S has potential consequences for malaria transmission.

De novo transcriptome sequencing in Anopheles funestus using Illumina RNA-seq technology

BACKGROUND

Anopheles funestus is one of the primary vectors of human malaria, which causes a million deaths each year in sub-Saharan Africa. Few scientific resources are available to facilitate studies of this mosquito species and relatively little is known about its basic biology and evolution, making development and implementation of novel disease control efforts more difficult. The An. funestus genome has not been sequenced, so in order to facilitate genome-scale experimental biology, we have sequenced the adult female transcriptome of An. funestus from a newly founded colony in Burkina Faso, West Africa, using the Illumina GAIIx next generation sequencing platform.

METHODOLOGY/PRINCIPAL FINDINGS

We assembled short Illumina reads de novo using a novel approach involving iterative de novo assemblies and "target-based" contig clustering. We then selected a conservative set of 15,527 contigs through comparisons to four Dipteran transcriptomes as well as multiple functional and conserved protein domain databases. Comparison to the Anopheles gambiae immune system identified 339 contigs as putative immune genes, thus identifying a large portion of the immune system that can form the basis for subsequent studies of this important malaria vector. We identified 5,434 1:1 orthologues between An. funestus and An. gambiae and found that among these 1:1 orthologues, the protein sequence of those with putative immune function were significantly more diverged than the transcriptome as a whole. Short read alignments to the contig set revealed almost 367,000 genetic polymorphisms segregating in the An. funestus colony and demonstrated the utility of the assembled transcriptome for use in RNA-seq based measurements of gene expression.

CONCLUSIONS/SIGNIFICANCE

We developed a pipeline that makes de novo transcriptome sequencing possible in virtually any organism at a very reasonable cost ($6,300 in sequencing costs in our case). We anticipate that our approach could be used to develop genomic resources in a diversity of systems for which full genome sequence is currently unavailable. Our An. funestus contig set and analytical results provide a valuable resource for future studies in this non-model, but epidemiologically critical, vector insect.

Association mapping of insecticide resistance in wild Anopheles gambiae populations: major variants identified in a low-linkage disequilbrium genome

Background

Association studies are a promising way to uncover the genetic basis of complex traits in wild populations. Data on population stratification, linkage disequilibrium and distribution of variant effect-sizes for different trait-types are required to predict study success but are lacking for most taxa. We quantified and investigated the impacts of these key variables in a large-scale association study of a strongly selected trait of medical importance: pyrethroid resistance in the African malaria vector Anopheles gambiae.

Methodology/Principal Findings

We genotyped ≈1500 resistance-phenotyped wild mosquitoes from Ghana and Cameroon using a 1536-SNP array enriched for candidate insecticide resistance gene SNPs. Three factors greatly impacted study power. (1) Population stratification, which was attributable to co-occurrence of molecular forms (M and S), and cryptic within-form stratification necessitating both a partitioned analysis and genomic control. (2) All SNPs of substantial effect (odds ratio, OR>2) were rare (minor allele frequency, MAF<0.05). (3) Linkage disequilibrium (LD) was very low throughout most of the genome. Nevertheless, locally high LD, consistent with a recent selective sweep, and uniformly high ORs in each subsample facilitated significant direct and indirect detection of the known insecticide target site mutation kdr L1014F (OR≈6; P<10⁻⁶), but with resistance level modified by local haplotypic background.

Conclusion

Primarily as a result of very low LD in wild A. Gambiae, LD-based association mapping is challenging, but is feasible at least for major effect variants, especially where LD is enhanced by selective sweeps. Such variants will be of greatest importance for predictive diagnostic screening.

Polymorphisms in Anopheles gambiae immune genes associated with natural resistance to Plasmodium falciparum

Many genes involved in the immune response of Anopheles gambiae, the main malaria vector in Africa, have been identified, but whether naturally occurring polymorphisms in these genes underlie variation in resistance to the human malaria parasite, Plasmodium falciparum, is currently unknown. Here we carried out a candidate gene association study to identify single nucleotide polymorphisms (SNPs) associated with natural resistance to P. falciparum. A. gambiae M form mosquitoes from Cameroon were experimentally challenged with three local wild P. falciparum isolates. Statistical associations were assessed between 157 SNPs selected from a set of 67 A. gambiae immune-related genes and the level of infection. Isolate-specific associations were accounted for by including the effect of the isolate in the analysis. Five SNPs were significantly associated to the infection phenotype, located within or upstream of AgMDL1, CEC1, Sp PPO activate, Sp SNAKElike, and TOLL6. Low overall and local linkage disequilibrium indicated high specificity in the loci found. Association between infection phenotype and two SNPs was isolate-specific, providing the first evidence of vector genotype by parasite isolate interactions at the molecular level. Four SNPs were associated to either oocyst presence or load, indicating that the genetic basis of infection prevalence and intensity may differ. The validity of the approach was verified by confirming the functional role of Sp SNAKElike in gene silencing assays. These results strongly support the role of genetic variation within or near these five A. gambiae immune genes, in concert with other genes, in natural resistance to P. falciparum. They emphasize the need to distinguish between infection prevalence and intensity and to account for the genetic specificity of vector-parasite interactions in dissecting the genetic basis of Anopheles resistance to human malaria.

Anopheles gambiae is the main malaria vector in Africa, transmitting the parasite when it blood feeds on human hosts. The parasite undergoes several developmental stages in the mosquito to complete its life cycle, during which time it is confronted by the mosquito's immune system. The resistance of mosquitoes to malaria infection is highly variable in wild populations and is known to be under strong genetic control, but to date the specific genes responsible for this variation remain to be identified. The present study uncovers variations in A. gambiae immune genes that are associated with natural resistance to Plasmodium falciparum, the deadliest human malaria parasite. The association of some mosquito genetic loci with the level of infection depended on the P. falciparum isolate, suggesting that resistance is determined by interactions between the genome of the mosquito and that of the parasite. This finding highlights the need to account for the natural genetic diversity of malaria parasites in future research on vector-parasite interactions. The loci uncovered in this study are potential targets for developing novel malaria control strategies based on natural mosquito resistance mechanisms.

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.

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.

Molecular diversity of antimicrobial effectors in the oyster Crassostrea gigas

Background

To gain insight into the molecular diversity of antimicrobial peptides and proteins in the oyster Crassostrea gigas, we characterized and compared the sequence polymorphism of the antimicrobial peptides (AMPs), Cg-Defensins (Cg-Defs) and Cg-Proline Rich peptide (Cg-Prp), and of the bactericidal permeability increasing protein, Cg-BPI. For that, we analyzed genomic and transcript sequences obtained by specific PCR amplification and in silico searches.

Results

High diversification among the three antimicrobial effectors was evidenced by this polymorphism survey. On the basis of sequence phylogenies, each AMP aggregates into clearly defined groups of variants and is the product of a multigenic family displaying a variety of gene structures. In contrast, Cg-bpi forms a single group and is encoded by a single gene copy. Moreover, we identified for both AMPs several genetic mechanisms of diversification such as recombination, parallel mutations leading to phylogenetic homoplasy and indel events. In addition, the non synonymous to synonymous substitutions ratio by codon (dN/dS) revealed several negatively and positively selected sites for both AMPs, suggesting that directional selection pressures have shaped their sequence variations.

Conclusions

This study shows for the first time in a mollusc that antimicrobial peptides and proteins have been subject to distinct patterns of diversification and we evidence the existence of different evolutionary routes leading to such sequence variability.

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.

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.

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.