Research | Open | Published:
Identification of three single nucleotide polymorphisms in Anopheles gambiae immune signaling genes that are associated with natural Plasmodium falciparum infection
Malaria Journalvolume 9, Article number: 160 (2010)
Laboratory studies have demonstrated that a variety of immune signaling pathways regulate malaria parasite infection in Anopheles gambiae, the primary vector species in Africa.
To begin to understand the importance of these associations under natural conditions, an association mapping approach was adopted to determine whether single nucleotide polymorphisms (SNPs) in selected immune signaling genes in A. gambiae collected in Mali were associated with the phenotype of Plasmodium falciparum infection.
Three SNPs were identified in field-collected mosquitoes that were associated with parasite infection in molecular form-dependent patterns: two were detected in the Toll5B gene and one was detected in the gene encoding insulin-like peptide 3 precursor. In addition, one infection-associated Toll5B SNP was in linkage disequilibrium with a SNP in sequence encoding a mitogen-activated protein kinase that has been associated with Toll signaling in mammalian cells. Both Toll5B SNPs showed divergence from Hardy-Weinberg equilibrium, suggesting that selection pressure(s) are acting on these loci.
Seven of these eight infection-associated and linked SNPs alter codon frequency or introduce non-synonymous changes that would be predicted to alter protein structure and, hence, function, suggesting that these SNPs could alter immune signaling and responsiveness to parasite infection.
The causative agents of malaria are protozoan parasites of the genus Plasmodium, which are transmitted to humans by anopheline mosquitoes. The highest number of cases occurs in sub-Saharan Africa, where the most deadly parasite Plasmodium falciparum is transmitted primarily by Anopheles gambiae sensu stricto. With over 500 million new cases a year  and over half of the world's population at risk for malaria, increased understanding of the complex factors governing transmission is a pressing need.
Extensive genetic structuring among natural populations of A. gambiae is likely to have an impact on parasite transmission. In particular, frequencies of paracentric inversions on the right arm of the second chromosome (2R) have revealed that many A. gambiae populations deviate strongly from Hardy-Weinberg Equilibrium. On the basis of these findings, Mopti, Savanna, and Bamako chromosomal forms were described . Genetic differentiation among these chromosomal forms revealed the presence of molecular forms characterized by fixed nucleotide differences in the intergenic spacer of the X-linked ribosomal DNA . In Mali, the M molecular form corresponds to the Mopti chromosomal form and the S molecular form corresponds to the Savanna and Bamako forms. The distribution of the M molecular form appears mostly limited to West and Central Africa, whereas the S molecular form is found throughout the range of A. gambiae . Gene flow between the S and Mopti-M forms is severely restricted [5, 6].
In addition to chromosomal complexity, the genome of A. gambiae is characterized by high level nucleotide polymorphism. The published genome of A. gambiae PEST strain was reported to have a frequency of single nucleotide polymorphisms (SNPs) of 1.6 × 10-3 or approximately 1 SNP per 625 bp . Morlais et al  found similarly high levels of polymorphism in other laboratory strains of A. gambiae. The most recent genome assembly of A. gambiae in the Ensembl database  reports a SNP density of 1 per 247 bp.
In 35 genes putatively involved in mosquito-pathogen interactions and behaviour, Morlais et al  identified 460 SNPs, 140 of which encode nonsynonymous substitutions. The authors also examined SNP frequencies in genes encoding fibronectin, thioester-containing protein 3 (TEP3), and peptidoglycan recognition protein (PGRP) in field-caught A. gambiae from sites in Senegal, Burkina Faso and Cameroon where the M and S molecular forms are sympatric . Although the genes selected by Morlais et al  have been shown to function in anti-pathogen defense in A. gambiae , no data on SNP associations with parasite infection and molecular form - a key variable based on aforementioned predictions of gene flow [5, 6]- were presented.
In addition to studies that have reported SNPs in genes implicated in anti-parasite immunity, a number of studies have taken a genome-wide approach to identify factors that regulate parasite transmission. In particular, Riehle et al  determined that colocalized quantitative trait loci (QTL) within a region of chromosome 2L in A. gambiae collected from a single site in Mali accounted for a significant amount of the variation in P. falciparum infection. The authors designated these QTL as a Plasmodium-resistance island or PRI, but also acknowledged that there were numerous other loci throughout the genome that likely contributed to variation in infection . Riehle et al  also hypothesized that resistance is the ancestral phenotype, whereas susceptibility to infection with P. falciparum results from mutations in the A. gambiae genome that result in failure of anti-parasite immunity.
Based on these general concepts - that A. gambiae populations exhibit significant genetic structuring, that SNPs are present in anti-pathogen genes in natural populations of A. gambiae, that susceptibility to infection may result from a failure of anti-parasite immunity- and the use of association mapping to link SNPs in functional genes with disease states , it was hypothesized that SNPs predicted to alter immune signaling could be linked to natural P. falciparum infection in A. gambiae. Furthermore, the strength of these associations may be a reflection of the different genetic backgrounds of A. gambiae, a phenomenon that will dictate the success of genetic interventions (e.g., creation of refractory mosquito strains) to block parasite transmission in natural populations .
To begin testing these hypotheses, SNPs in a subset of immune signaling genes were identified using direct sequencing of conserved domains, and the association of these SNPs with P. falciparum infection in A. gambiae collected in Mali were analysed. This association study revealed significant molecular form-dependent associations of P. falciparum infection with three SNPs in genes encoding Toll5B and an insulin-like peptide 3 precursor. In addition, these analyses revealed an epistatically associated SNP in a gene encoding MKK4 - a mitogen-activated protein kinase that has been predicted to function in the Toll signaling pathway [14–16]. Hence, these data support the proposed hypotheses and provide new insights on regulatory factors that are associated with control of P. falciparum infection in A. gambiae.
Blood fed female mosquitoes were collected from five villages in Mali in 2006 and 2007 (Table 1). These villages were selected based on previous data, which identified them as populations in which M and S forms are sympatric. This facilitated acquisition of adequate numbers of each form from a single environment. Mosquitoes were collected between July 28 and September 2 in 2006 and between August 8 and September 8 in 2007 to reduce genotype frequency changes related to season. Collected mosquitoes were dissected so that the head and thorax were separated from the abdomen. Head and thorax samples were subjected to enzyme-linked immunosorbent assay (ELISA) to determine P. falciparum infection status and the paired abdomens were used for species identification and SNP genotyping.
Species identification and ELISA for P. falciparum infection
Mosquitoes were identified to species using the PCR assay described by Scott et al  and to molecular form using the restriction fragment length polymorphism protocol of Fanello et al . Briefly, following HhaI restriction of an amplified fragment of the X-linked ribosomal DNA, M molecular form individuals were identified by the presence of a 367 bp fragment, while S molecular form individuals were identified by fragments of 110 bp and 257 bp . The presence of P. falciparum in heads and thoraces was determined using a "sandwich" circumsporozoite protein (CSP)-based ELISA [18, 19]. All ELISA-positive specimens were rescreened against a standard curve of recombinant CSP to estimate relative infection levels.
Genes encoding Toll5B, MKK4 and insulin-like peptide precursor 3 were selected for analysis based on established roles in anti-pathogen signaling in A. gambiae or in other mosquito species [20–25]. Primers were designed to amplify fragments of the encoding sequences of the genes of interest in the size range of 120-260 bp, which was optimal for Luminex genotyping (see below).
For PCR, genomic DNA was isolated from A. gambiae abdominal tissue using the Qiagen BioSprint 96 (Valencia, CA) following the manufacturer's protocol. Each amplification reaction contained 1× buffer with MgCl2 (Roche), 67 μM of each dNTP (Applied Biosystems, Foster City, CA), 0.22 μM primers (Sigma-Aldrich, St. Louis, MO), 0.7 U Taq Polymerase (Roche), and 40 ng of genomic DNA. Cycling conditions were as follows: 94°C for 2 min; 2 cycles of 94°C for 30 sec, 62°C for 30 sec, 72°C for 30 sec; 6 cycles of 94°C for 30 sec, 55-61°C (the range of target-specific annealing temperatures) for 30 sec, 72°C for 30 sec; 27 cycles of 94°C for 30 sec, 55°C for 30 sec, 72°C for 30 sec, followed by a 7 minute extension at 72°C. Amplimer sizes were confirmed by electrophoresis prior to purification with Exo SapIT (; USB, Cleveland, OH) for DNA sequencing.
Initial amplifications were performed using six sets of pooled genomic DNA samples from each of four P. falciparum-infected A. gambiae (n = 24) and six sets of pooled genomic DNA samples from each of four uninfected control A. gambiae (n = 24) collected in the same village on the same day. This pooling strategy has been shown to be very accurate (>98%) for detecting SNPs in individuals . Pooled sequences were aligned using MegAlign (DNAStar Lasergene 8.0, Madison, WI) software and the ClustalW algorithm . Each potential SNP was confirmed manually with the chromatogram of the pool. When a polymorphism was detected in the pooled A. gambiae amplimers, each of the four genomic DNAs in the pool was subjected to separate re-amplification and sequencing to determine if the locus was indeed polymorphic, and if so, which individuals were polymorphic and whether the individuals were homozygous or heterozygous.
Individual sequences from the 48 genomic DNA samples used for SNP discovery were aligned using MegAlign (DNAStar Lasergene 8.0, Madison, WI) software and the ClustalW algorithm . Loci demonstrating divergence from the consensus sequences were designated as potential SNPs. The existence of each SNP was confirmed by manual examination of individual chromatograms. All SNPs were analysed by SIFT  and pMUT  to determine the likelihood that a particular SNP would affect protein function.
SNPs for Luminex analysis  were chosen based on the presence of a minimum of 20 bp of flanking sequence that was free of SNPs. Based on this requirement, we selected 11 SNP loci out of 96 SNPs for analysis: insulin-like peptide 3 precursor loci 2-5 (AGAP010602; Ins32, Ins33, Ins34, Ins35); MKK4 loci 1 and 3 (AGAP003365; MKK41, MKK43) and Toll5B loci 1-4 and 6 (AGAP010669; Toll5B1, Toll5B2, Toll5B3, Toll5B4, Toll5B6). Allele-specific primers were designed for each of the 11 SNPs (Additional file 1). A total of 134 P. falciparum-infected A. gambiae and 134 uninfected A. gambiae collected in the same village on the same day were analysed with Luminex (Additional file 1). The 48 individuals used for SNP discovery were included among these samples, providing internal quality control checks for Luminex performance in detecting heterozygotes and homozygotes at each of the dimorphic loci.
Allele-specific primer extension (ASPE)
For each Luminex assay, SNP-containing amplimers were combined in a final volume of 10 μ l. This mixture was treated with shrimp alkaline phosphatase (SAP) and exonuclease (Exo I) to remove unused PCR primers and dNTPs. For each reaction, 1 U of SAP (Amersham Biosciences, Inc., Piscataway, New Jersey) and 5 U Exo I (USB Corp., Cleveland, Ohio) were mixed with 10 μ l of amplified product and incubated at 37°C for 30 minutes, followed by 80°C for 15 minutes to inactivate the enzymes.
Ten μ l aliquots of ASPE master mix containing 500 nM of each ASPE primer, 50 mM MgCl2, 10× buffer (supplied with enzyme), 100 μ M each dATP, dCTP and dGTP, 400 μ M biotin-dCTP (Invitrogen, Carlsbad, California) and 0.8 U Platinum Tsp DNA Polymerase (Invitrogen) were dispensed to the pooled, SAP/Exo I-treated PCR products using a Biomek® 2000 Laboratory Automation Workstation (Beckman Coulter, Fullerton, California). The ASPE reaction was performed in a PTC-225 Peltier Thermal Cycler (MJ Research, Watertown, Massachusetts) and cycling parameters consisted of an initial denaturation at 96°C for 2 minutes, followed by 94°C for 30 seconds and 56°C for 2 minutes, repeated 49 times, ending with 72°C for 5 minutes.
Hybridization to FlexMAP beads
Individual FlexMAP™ bead types (MiraiBio, Alameda, California) were obtained at a concentration of 2.5 × 105 microspheres/ml. After resuspension by vortexing for 5 minutes, 2.5 μ l of each bead type per sample was used to make up the appropriate bead mix for each of the three multiplex assays. Each bead mix was concentrated by centrifugation at 10,000 × g for three minutes followed by careful removal of the supernatant. Beads were then resuspended in 2× TM hybridization buffer (0.4 M NaCl, 0.2 M Tris HCl pH 8.0, 0.16% Triton X-100) and dH2O, and added to each sample for a final volume of 50 μ l (resulting in 625 beads/allele in 1× TM). Samples were denatured at 96°C for 90 seconds, followed by hybridization at 52°C, 47°C and 37°C for 30 minutes at each temperature. The hybridized beads were washed twice by centrifugation at 3,000 × g for 3 minutes, removal of the supernatant and resuspension in 70 μ l of 1× TM buffer. After centrifugation and removal of the supernatant for a third time, beads were resuspended in 70 μ l 1× TM buffer containing 8 μ g/ml streptavidin-phycoerythrin (ProZyme, San Leandro, California), a fluorescent reporter molecule used to detect the ASPE-incorporated biotin.
Sample analysis on the Luminex 100 System
Samples were analysed on the Luminex 100 System using Data Collection Software Version 1.7 with settings specified by the manufacturer; the median fluorescence intensity (signal) was measured over 100 independent events (beads). The genotype of each SNP locus was determined by the ratio of fluorescence intensity of allele A, I A , and that of allele B, I B . If I A /I B was greater than 3.5, genotype was set as AA; if I A /I B was less than 0.5, genotype was set as BB; for other ratio values, the genotype was set as AB. These thresholds were confirmed by comparison to the direct sequencing data available for each of the polymorphic loci.
The pwr package for R http://cran.r-project.org/web/packages/pwr/index.html was used for power analysis. This package calculated effect size, required sample size, and power for each SNP using the methods of Cohen . Significant differences in allele frequency (1) between infected and uninfected samples within a population and (2) between any two populations were determined using Fisher's exact test (for 2 by 2 contingency tables) or Chi-square test implemented in R statistics package http://www.r-project.org. Significance thresholds were adjusted for multiple comparisons .
Unaccounted population structure can lead to the discovery of spurious associations or dilute true associations [34, 35]. To eliminate associations due to population structure, we tested whether our SNP genotypes showed significant divergence based on geographic location and molecular form. Samples were first divided into eight groups based on collection site and molecular form. Next, pairwise FST values were calculated using Arlequin v 3.11 . The two groups with the highest p-values were then combined and pair-wise FST values were recalculated. If genetic divergence, as described by FST, between a pair of groups is minimal, the p-value is close to 1. This was repeated until all between-group FST values were significant. The eight initial groups (based on collection site and molecular form) were ultimately reduced to 3 distinct groups: M, S1, and S2, as described in the Results. Within each of the three groups, linkage disequilibrium was calculated using Arlequin. Phylip v 3.68 [37, 38] was used for phylogenetic tree construction.
Among the 11 SNPs analysed, the frequency of exact matches between direct sequencing and Luminex for individual genotypes (n = 48) was greater than 98%, indicating that Luminex was highly accurate for genotyping.
Based on SNP genotypes at the 11 loci selected, at least three genetically distinct populations in Mali (Figure 1) were identified. In agreement with population genetic studies based on microsatellites [39, 40], significant genetic differentiation between M and S molecular forms was observed. Further, Pimperena S (S2) form mosquitoes were differentiated from S forms collected in other villages (S1). This outcome may be due to temporal variation in the relative abundance of the Bamako and Savanna chromosomal forms, both of which are the S molecular form and therefore indistinguishable using the Fanello et al  molecular form diagnostic.
The association between SNP genotypes at each locus and infection status within each of the three populations indicated in Figure 1 was evaluated using Chi-square and Fisher's exact tests. Infection rates of the genotyped A. gambiae ranged from 3.27% in Selinkenyi to a high of 15.09% in Doneguebougou (Table 1). Genotype frequencies at the Ins35 locus were significantly different between the Pimperena (S2) and S1 groups, while genotype frequencies at the Toll5B1 locus were significantly different between S2 and all other groups (M and S1 in Table 2). All SNPs identified to be associated with infection status had an effect size (h) greater than 1 (Table 3) and power greater than 0.85. Those SNPs that were determined to be associated with population groups had power greater than 0.83 (Table 3). SNPs with smaller h, such as Ins32 and Ins35 (0.75 < h < 1), could be associated with infection if a larger sample size was scored. For example, if 30 infected specimens and 62 uninfected specimens for Ins32 were genotyped and a similar genotype frequency distribution was observed, Ins32 could be identified as significantly associated with infection status. The sample sizes of this study were not sufficient to confirm involvement of SNPs with smaller h - such as Ins32 or Ins35 - in the regulation of P. falciparum development in A. gambiae.
Three SNPs - Ins34 (3L), Toll5B1 (3L), and Toll5B6 (3L) (Table 2) - were found to be associated with P. falciparum infection status. Ins34 is a synonymous SNP resulting in a change from GGC to GGT at nucleotide position 462 in the Insulin-like peptide 3 precursor gene. Toll5B1 introduces a synonymous SNP at nucleotide position 129, changing the codon from ATC, a common codon (28.3%), to ATT, a rarer codon (12.9%). Toll5B6 results in a non-synonymous codon change from AGC (Ser, S) to AAC (Asn, N) at amino acid position 454 in the predicted translation of Toll5B. These SNPs were not identified in either the NCBI or Anobase SNP databases generated by reads from the PEST strain (10× coverage) or from the MOPTI strain (1.2× coverage) aligned to Celera A. gambiae contigs, although Ins34 and Toll5B6 were both present in the M and S form scans available from VectorBase . However, these loci were not identified to be polymorphic in either scan . The CC genotype at the Ins34 locus in M form mosquitoes was more common in samples that were not infected with P. falciparum. In both M and S1 populations, Toll5B1 C alleles were significantly associated with infection status. While the Toll5B1 TT genotype was rare, TC heterozygotes were more common in uninfected A. gambiae than in P. falciparum-infected A. gambiae. Within the S1 population, individuals with the AG genotype at the Toll5B6 locus were less likely to be infected with P. falciparum.
Significant differences in median P. falciparum sporozoite infection intensities were observed in the collected samples (Kruskal-Wallis rank sum test, p = 1.88x10-5, Figure 2, Table 4). However, no significant association between P. falciparum infection intensity and any SNP locus genotype was observed. The sporozoite infection intensities were highly variable in M and S1 populations as indicated by the standard deviations being more than double their associated means shown in Table 4. Among infected mosquitoes, samples from Pimperena (S2) had higher median sporozoite intensities than those from other villages. VectorBase A. gambiae population data indicated that samples from Pimperena are mostly Savanna chromosomal form while other sites like Doneguebougou and Selinkenyi are composed of a mixture of Bamako and Savanna forms. Thus, the possibility of chromosomal form affecting the density of P. falciparum in the S1 group cannot be ruled out.
Many factors contribute to mosquito infection and successful transmission of malaria parasites, including innate immunity. As such, considerable efforts have been focused on understanding the mosquito immune system and natural variability in these defenses. In humans, numerous reports have demonstrated that natural variability to infection can be associated with single nucleotide polymorphisms (SNPs) in genes that regulate host immunity [42–44]. Although SNP associations have been described for a variety of human infections and diseases, these have not been identified for P. falciparum infection in genetically defined natural populations of A. gambiae until this study.
The Toll and Imd signaling pathways are important regulators of innate immunity in A. gambiae. In particular, Garver et al  reported that, under laboratory conditions, the Toll pathway controlled P. falciparum infection intensity in A. gambiae, while the Imd pathway appeared to regulate resistance to infection. Two Toll5B SNPs - Toll5B1 and Toll5B6 - were significantly associated with P. falciparum infection status among field-collected A. gambiae. Anopheles gambiae Toll5B is orthologous to Drosophila melanogaster Toll-5, also known as Tehao , and Toll5 in Aedes aegypti . In Ae. aegypti, Toll5B was inducibly expressed in the mosquito fat body following fungal infection . In A. gambiae, Pinto et al  reported that Toll5B expression is significantly upregulated in hemocytes of adult female mosquitoes infected with Plasmodium berghei at 24-28 hours post-infection, a period associated with active parasite invasion of the midgut epithelium. Based on these observations and these data, Toll5B is likely to be responsive to P. falciparum infection in A. gambiae under natural conditions.
Both Toll5B SNPs showed divergence from Hardy-Weinberg equilibrium. In particular, Toll5B1 showed significant divergence from Hardy-Weinberg expectation in both the M (p = 0.00016) and S1 (p = 0.00077) forms, while Toll5B6 showed significant divergence from Hardy-Weinberg expectation in the M form (p = 0.00035; Table 5), suggesting that selection pressure(s) acting on these loci may skew the genotype frequencies of this gene. For Toll5B1, the observed heterozygosity was greater than expected, both in M and S1 groups (Table 5). For this SNP, CT heterozygotes were more common (96%) in uninfected samples while 72% of infected samples were CC homozygotes. TT homozygotes were very rare and only one sample of this genotype was identified in the entire collection. These circumstances indicate the possibility of heterosis or balancing selection on this locus. For Toll5B6, the GG homozygote was the most common form in all groups and the proportion of GG was highest in the M uninfected group. The effect size of this SNP was 0.514, so a larger sample size (83 infected and 183 uninfected) would be needed (Table 3) to confirm involvement of this SNP with respect to infection status in the M form population.
Toll5B1 occurs in a highly polymorphic region, containing 13 polymorphic loci within a stretch of 166 bp. Six of these SNPs occur in clusters of rare codons or result in the change from a common codon to a rare one. Clusters of rare codons have been shown to alter protein production in a synergistic manner [49–51]. Genotype frequencies of Toll5B1 and infection associations of this SNP were significantly different between Pimperena (S2) and A. gambiae M and S1 groups (Table 2). Taken together with observations on codon frequency and protein function, these data suggest that within the M and S1 molecular forms, this SNP may have a functional effect on P. falciparum infection, but this association may be population-specific, as indicated by the lack of any Toll5B1-P. falciparum infection association in the Pimperena (S2) population. Toll5B6 is a non-synonymous SNP that introduces the amino acid change S454N. The replacement of serine by the more bulky asparagine could introduce subtle changes in the 3-dimensional structure of the Toll leucine-rich repeat (LRR) in which this mutation occurs.
Additional analyses revealed that MKK43, a SNP in the MAPK kinase MKK4 gene, located on chromosome 2, was in linkage disequilibrium with infection-associated Toll5B1, a SNP on chromosome 3, in the S1 population (p = 0.0001). In mammals, MKK4 can be activated by TLR3  and TLR2 signaling [14, 15], indicating that MKK4 activation is functionally linked to Toll signaling. Downstream of this activation, MKK4 functions with MKK7 as the primary activator of c-Jun N-terminal kinase or JNK, one of three immunity-associated MAPKs. In D. melanogaster, Toll activation of MKK4/7 and JNK during septic injury regulates cytoskeletal genes typically associated with a wound healing response to infection . A potential Toll/MKK4/JNK signaling module - if biologically functional in A. gambiae - could be linked to the profound cytoskeletal changes in the midgut epithelium that have been described during P. berghei and P. falciparum infection of laboratory and field specimens of A. gambiae, respectively [53–55]. This possibility is currently being investigated.
In mammals, the insulin/insulin-like growth factor signaling cascade (IIS) has been shown to regulate innate immunity through Toll- and NF-κB-dependent pathways [56, 57]. In particular, IIS activation can induce or inhibit NF-κB-dependent signaling and is, thereby, capable of exerting both pro- and anti-inflammatory effects on the host immune response. In Anopheles stephensi, control of malaria parasite development is regulated by signaling proteins associated with the IIS cascade [, Corby-Harris et al. unpublished]. The IIS is highly conserved  and critical components of the cascade, as well as a variety of insulin-like peptides, including insulin-like peptide 3 precursor, are expressed in the midgut of A. gambiae [; Luckhart, unpublished].
Four SNPs in the insulin-like peptide 3 precursor gene - Ins32, Ins33, Ins34 and Ins35 - were analysed and significant infection (Ins34) and molecular form (Ins35) associations were found for two of these (Table 2). Ins34 introduces a synonymous mutation, with no remarkable change in codon frequency. Although Ins34 is predicted to have little to no effect on protein function, this infection-associated SNP was in linkage disequilibrium in the S1 population with Toll5B2 (p = 0.03307, Additional file 2), a SNP that introduces the non-synonymous mutation D56A into Toll5B. This SNP encodes a mutation close to the N-terminus and outside of the predicted LRRs, in an ectodomain region identified in other Toll proteins as the cysteine-rich capping structure . The N-terminal capping structure may participate in protein-protein interactions that are critical for Toll receptor function .
Within the S1 population, Ins32, Ins33, and Ins35 showed levels of linkage disequilibrium with Toll5B2 similar to that observed for Ins34 (p = 0.01238, 0.02931, 0.00851, respectively, Additional file 2), which likely reflects the fact that these loci are closely linked physically. Although Ins35 was not significantly associated with infection status (Table 2), this SNP was also in linkage disequilibrium with Toll5B2 in the S2 population (p-value = 0.00406, Additional file 2). However, in the S2 population, the neighboring SNPs - Ins32, Ins33, Ins34 - did not show similar levels of linkage disequilibrium with Toll5B2 (p = 0.09584, 0.14396, 0.13238, respectively, Additional file 2), suggesting that the associations of Ins34 and Ins35 with Toll5B2 may, in fact, represent novel biological functionality that is population-specific in A. gambiae.
The positional effects of SNPs in the architecture of signaling cascades have been investigated in a series of relevant studies. Riley et al  examined SNPs in the D. melanogaster Ras-mediated signal transduction pathway. The least polymorphic signaling protein genes (e.g., Ras, Dr, and Polehole) were those that were proximal to the origin of signaling at the cell surface, while the most polymorphic genes were located farther downstream in the signaling cascade (e.g., Dsor1, Csw and Ksr; ). Computational simulations of MAPK signaling predicted similar constraints in that if signal amplification was crucial, the upstream signaling proteins were more constrained than were the downstream components . Together these findings suggest that a SNP in an upstream component of a signaling cascade, such as Toll or an ILP, could have pleiotropic effects on the downstream components of the signaling cascade and an increased potential for a deleterious outcome . In studies of the Tor signaling pathway, however, selection constraints were found to be greater in the downstream components rather than the upstream components of the cascade . Regardless of the pattern of selection constraints, the polarity of a signal transduction pathway is an integral part of determining the effect of a SNP. As such, the molecular cell biology of the signaling pathway protein networks highlighted herein can complement studies of phenotype-associated SNPs in natural A. gambiae populations.
In summary, the infection-associated SNPs that have been identified here could be used as genetic markers for the susceptibility of an anopheline population to malaria parasite infection. More importantly, however, these SNPs confirm previous laboratory associations of the target genes with the regulation of parasite infection in A. gambiae, extending possible functional linkages among these target genes and confirming the influence of population structuring on biological associations of the SNPs under study. In particular, these data together with associated laboratory data that implicate these gene products in anti-parasite immunity suggest that Toll5B, MKK4 and insulin-like peptide 3 alone or perhaps in some combination may be involved in the regulation of P. falciparum development in A. gambiae under natural conditions. The alternative hypothesis - that the true functional genes are linked to these SNPs - will be evaluated as haplotype maps are developed. In addition, these data have revealed that SNPs are not equally distributed among A. gambiae molecular forms and that selection pressures may be driving deviations from Hardy-Weinberg equilibrium. These insights suggest that gene flow may impede the spread of transgenes under field conditions and that selection may alter the distribution of transgenes, factors that must be accommodated in the development of any field strategies.
The World Health Organization: The World Malaria Report. 2005
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maiga HM, Sankaré O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40: 477-511.
Fanello C, Santolamazza F, della Torre A: Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001, 10: 9-18. 10.1046/j.1365-2583.2001.00235.x.
della Torre A, Tu ZJ, Petrarca V: On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms. Insect Biochem Mol Biol. 2005, 35: 755-769. 10.1016/j.ibmb.2005.02.006.
Taylor C, Toure YT, Carnahan J, Norris DE, Dolo G, Traore SF, Edillo FE, Lanzaro GC: Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali, west Africa. Genetics. 2001, 157: 743-750.
Tripet F, Toure YT, Taylor CE, Norris DE, Dolo G, Lanzaro GC: DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae. Mol Ecol. 2001, 10: 1725-1732. 10.1046/j.0962-1083.2001.01301.x.
Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, Cai S, Center A, Chaturverdi K, Christophides GK, Chrystal MA, Clamp M, Cravchik A, Curwen V, Dana A, Delcher A, Dew I, Evans CA, Flanigan M, Grundschober-Freimoser A, Friedli L, Gu Z, Guan P, Guigo R, Hillenmeyer ME, Hladun SL, Hogan JR, Hong YS, Hoover J, Jaillon O, Ke Z, Kodira C, Kokoza E, Koutsos A, Letunic I, Levitsky A, Liang Y, Lin JJ, Lobo NF, Lopez JR, Malek JA, McIntosh TC, Meister S, Miller J, Mobarry C, Mongin E, Murphy SD, O'Brochta DA, Pfannkoch C, Qi R, Regier MA, Remington K, Shao H, Sharakhova MV, Sitter CD, Shetty J, Smith TJ, Strong R, Sun J, Thomasova D, Ton LQ, Topalis P, Tu Z, Unger MF, Walenz B, Wang A, Wang J, Wang M, Wang X, Woodford KJ, Wortman JR, Wu M, Yao A, Zdobnov EM, Zhang H, Zhao Q, Zhao S, Zhu SC, Zhimulev I, Coluzzi M, della Torre A, Roth CW, Louis C, Kalush F, Mural RJ, Myers EW, Adams MD, Smith HO, Broder S, Gardner MJ, Fraser CM, Birney E, Bork P, Brey PT, Venter JC, Weissenbach J, Kafatos FC, Collins FH, Hoffman SL: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298: 129-149. 10.1126/science.1076181.
Morlais I, Poncon N, Simard F, Cohuet A, Fontenille D: Intraspecific nucleotide variation in Anopheles gambiae: New insights into the biology of malaria vectors. Am J Trop Med Hyg. 2004, 71: 795-802.
Christophides GK, Zdobnov E, Barillas-Mury C, Birney E, Blandin S, Blass C, Brey PT, Collins FH, Danielli A, Dimopoulos G, Hetru C, Hoa NT, Hoffmann JA, Kanzok SM, Letunic I, Levashina EA, Loukeris TG, Lycett G, Meister S, Michel K, Moita LF, Müller HM, Osta MA, Paskewitz SM, Reichhart JM, Rzhetsky A, Troxler L, Vernick KD, Vlachou D, Volz J, von Mering C, Xu J, Zheng L, Bork P, Kafatos FC: Immunity-related genes and gene families in Anopheles gambiae. Science. 2002, 298: 159-165. 10.1126/science.1077136.
Riehle MM, Markianos K, Niare O, Xu JN, Li J, Toure AM, Podiougou B, Oduol F, Diawara S, Diallo M, Coulibaly B, Ouatara A, Kruglyak L, Traoré SF, Vernick KD: Natural malaria infection in Anopheles gambiae is regulated by a single genomic control region. Science. 2006, 312: 577-579. 10.1126/science.1124153.
Wang WYS, Barratt BJ, Clayton DG, Todd JA: Genome-wide association studies: Theoretical and practical concerns. Nat Rev Genet. 2005, 6: 109-118. 10.1038/nrg1522.
Tripet F, Aboagye-Antwi F, Hurd H: Ecological immunology of mosquito-malaria interactions. Trends Parasitol. 2008, 24: 219-227. 10.1016/j.pt.2008.02.008.
Vasselon T, Hanlon WA, Wright SD, Detmers PA: Toll-like receptor 2 (TLR2) mediates activation of stress-activated MAP kinase p38. J Leukoc Biol. 2002, 71: 503-510.
Sakai A, Han JH, Cato ACB, Akira S, Li JD: Glucocorticoids synergize with IL-1 beta to induce TLR2 expression via MAP kinase phosphatase-1-dependent dual inhibition of MAPK JNK and p38 in epithelial cells. Bmc Mol Bio. 2004, 5: 2-10.1186/1471-2199-5-2.
Yoshizawa T, Hanunaker D, Sweeney SE, Boyle DL, Firestein GS: Synoviocyte innate immune responses: I. Differential regulation of interferon responses and the JNK pathway by MAPK kinases. J Immunol. 2008, 181: 3252-3258.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain-reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Burkot TR, Williams JL, Schneider I: Identification of Plasmodium-falciparum-infected mosquitoes by a double antibody enzyme-linked immunosorbent-assay. Am J Trop Med Hyg. 1984, 33: 783-788.
Wirtz RA, Zavala F, Charoenvit Y, Campbell GH, Burkot TR, Schneider I, Esser KM, Beaudoin RL, Andre RG: Comparative testing of monoclonal antibodies against Plasmodium falciparum sporozoites for ELISA development. Bull World Health Organ. 1987, 65: 39-45.
Mizutani T, Kobayashi M, Eshita Y, Shirato K, Kimura T, Ako Y, Miyoshi H, Takasaki T, Kurane I, Kariwa H, Umemura T, Takashima I: Involvement of the JNK-like protein of the Aedes albopictus mosquito cell line, C6/36, in phagocytosis, endocytosis and infection of West Nile virus. Insect Mol Biol. 2003, 12: 491-499. 10.1046/j.1365-2583.2003.00435.x.
Mizutani T, Kobayashi M, Eshita Y, Inanami O, Yamamori T, Goto A, Ako Y, Miyoshi H, Miyamoto H, Kariwa H, Kuwabara M, Takashima I: Characterization of JNK-like protein derived from a mosquito cell line, C6/36. Insect Mol Biol. 2003, 12: 61-66. 10.1046/j.1365-2583.2003.00387.x.
Lim JH, Gowda DC, Krishnegowda G, Luckhart S: Induction of nitric oxide synthase in Anopheles stephensi by Plasmodium falciparum: Mechanism of signaling and the role of parasite glycosylphosphatidylinositols. Infect Immun. 2005, 73: 2778-2789. 10.1128/IAI.73.5.2778-2789.2005.
Luckhart S, Riehle MA: The insulin signaling cascade from nematodes to mammals: Insights into innate immunity of Anopheles mosquitoes to malaria parasite infection. Dev Comp Immunol. 2007, 31: 647-656. 10.1016/j.dci.2006.10.005.
Pinto SB, Koutsos AC, Waterhouse RM, McKay K, An C, Ramakrishnan C, Kafatos FC, Michel K: Discovery of Plasmodium modulators by genome-wide analysis of circulating hemocytes in Anopheles gambiae. Proc Natl Acad Sci USA. 2009, 106: 21270-21275. 10.1073/pnas.0909463106.
Surachetpong W, Singh N, Cheung KW, Luckhart S: MAPK ERK signaling regulates the TGF-beta 1-dependent mosquito response to Plasmodium falciparum. PloS Pathog. 2009, 5: e1000366-10.1371/journal.ppat.1000366.
Werle E, Schneider C, Renner M, Volker M, Fiehn W: Convenient single-step, one tube purification of PCR products for direct sequencing. Nucleic Acids Res. 1994, 22: 4354-4355. 10.1093/nar/22.20.4354.
Van Deynze A, Stoffel K, Buell CR, Kozik A, Liu J, van der Knaap E, Francis D: Diversity in conserved genes in tomato. BMC Genomics. 2007, 8: 465-10.1186/1471-2164-8-465.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W - improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Ng PC, Henikoff S: Predicting the effects of amino acid substitutions on protein function. Annu Rev Genomics Hum Genet. 2006, 7: 61-80. 10.1146/annurev.genom.7.080505.115630.
Ferrer-Costa C, Orozco M, de la Cruz X: Sequence-based prediction of pathological mutations. Proteins. 2004, 57: 811-819. 10.1002/prot.20252.
Dunbar SA: Applications of Luminex (R) xMAP (TM) technology for rapid, high-throughput multiplexed nucleic acid detection. Clin Chim Acta. 2006, 363: 71-82. 10.1016/j.cccn.2005.06.023.
Cohen J: Statistical Power Analysis for the Behavioral Sciences. 1988, Lawrence Erlbaum Associates, 2
Miller RG: Simultaneous Statistical Inference. 1981, New York: Springer Verlag, 2
Pritchard JK, Rosenberg NA: Use of unlinked genetic markers to detect population stratification in association studies. Am J Hum Genet. 1999, 65: 220-228. 10.1086/302449.
Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, Eskin E: Efficient control of population structure in model organism association mapping. Genetics. 2008, 178: 1709-1723. 10.1534/genetics.107.080101.
Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform. 2005, 47-50.
Felsenstein J: Mathematical Evolutionary-Theory. Science. 1989, 246: 941-942. 10.1126/science.246.4932.941.
Felsenstein J: Using the quantitative genetic threshold model for inferences between and within species. Philos Trans R Soc B Biol Sci. 2005, 360: 1427-1434. 10.1098/rstb.2005.1669.
Slotman MA, Tripet F, Cornel AJ, Meneses CR, Lee Y, Reimer LJ, Thiemann TC, Fondjo E, Fofana A, Traore SF, Lanzaro GC: Evidence for subdivision within the M molecular form of Anopheles gambiae. Mol Ecol. 2007, 16: 639-649. 10.1111/j.1365-294X.2006.03172.x.
Lee Y, Cornel AJ, Meneses CR, Fofana A, Andrianarivo AG, McAbee RD, Fondjo E, Traore SF, Lanzaro GC: Ecological and genetic relationships of the Forest-M form among chromosomal and molecular forms of the malaria vector Anopheles gambiae sensu stricto. Malar J. 2009, 8: 75-10.1186/1475-2875-8-75.
Badger SA, Soong CV, O'Donnell ME, Sharif MA, Makar RR, Hughes AE: Common polymorphisms of Fibulin-5 and the risk of abdominal aortic aneurysm development. Vasc Med. 2009, 15: 113-117. 10.1177/1358863X09355667.
Lamsyah H, Rueda B, Baassi L, Elaouad R, Bottini N, Sadki K, Martin J: Association of PTPN22 gene functional variants with development of pulmonary tuberculosis in Moroccan population. Tissue Antigens. 2009, 74: 228-232. 10.1111/j.1399-0039.2009.01304.x.
Wang W, Yuasa T, Tsuchiya N, Ma ZY, Maita S, Narita S, Kumazawa T, Inoue T, Tsuruta H, Horikawa Y, Saito M, Hu W, Ogawa O, Habuchi T: The novel tumor-suppressor Mel-18 in prostate cancer: Its functional polymorphism, expression and clinical significance. Int J Cancer. 2009, 125: 2836-2843. 10.1002/ijc.24721.
Garver LS, Dong YM, Dimopoulos G: Caspar controls resistance to Plasmodium falciparum in diverse anopheline species. PloS Pathog. 2009, 5: e1000335-10.1371/journal.ppat.1000335.
Luna C, Hoa NT, Zhang J, Kanzok SM, Brown SE, Imler JL, Knudson DL, Zheng LB: Characterization of three Toll-like genes from mosquito Aedes aegypti. Insect Mol Biol. 2003, 12: 67-74. 10.1046/j.1365-2583.2003.00388.x.
Shin SW, Bian GW, Raikhel AS: A toll receptor and a cytokine, Toll5A and Spz1C, are involved in toll antifungal immune signaling in the mosquito Aedes aegypti. J Biol Chem. 2006, 281: 39388-39395. 10.1074/jbc.M608912200.
Pinto SB, Koutsos AC, Waterhouse RM, McKay K, An C, Ramakrishnan C, Kafatos FC, Michel K: Discovery of Plasmodium modulators by genome-wide analysis of circulating hemocytes in Anopheles gambiae. Proc Natl Acad Sci USA. 2009, 106: 21270-21275. 10.1073/pnas.0909463106.
Varenne S, Lazdunski C: Effect of distribution of unfavorable codons on the maximum rate of gene-expression by an heterologous organism. J Theor Bio. 1986, 120: 99-110. 10.1016/S0022-5193(86)80020-0.
Varenne S, Baty D, Verheij H, Shire D, Lazdunski C: The maximum rate of gene-expression is dependent on the downstream context of unfavorable codons. Biochimie. 1989, 71: 1221-1229. 10.1016/0300-9084(89)90027-8.
Clarke T, Clark PL: Rare codons cluster. PLoS One. 2008, 3: e3412-10.1371/journal.pone.0003412.
Boutros M, Agaisse H, Perrimon N: Sequential activation of signaling pathways during innate immune responses in Drosophila. Dev Cell. 2002, 3: 711-722. 10.1016/S1534-5807(02)00325-8.
Mendes AM, Schlegelmilch T, Cohuet A, Awono-Ambene P, De Iorio M, Fontenille D, Morlais I, Christophides GK, Kafatos FC, Vlachou D: Conserved mosquito/parasite interactions affect development of Plasmodium falciparum in Africa. PloS Pathog. 2008, 4: e1000069-10.1371/journal.ppat.1000069.
Shiao SH, Whitten MMA, Zachary D, Hoffmann JA, Levashina EA: Fz2 and Cdc42 mediate melanization and actin polymerization but are dispensable for Plasmodium killing in the mosquito midgut. PloS Pathog. 2006, 2: 1152-1164. 10.1371/journal.ppat.0020133.
Han YS, Thompson J, Kafatos FC, Barillas-Mury C: Molecular interactions between Anopheles stephensi midgut cells and Plasmodium berghei: the time bomb theory of ookinete invasion of mosquitoes. EMBO J. 2001, 20: 1483-1483. 10.1093/emboj/20.23.6909.
Dandona P, Aljada A, Mohanty P, Ghanim H, Hamouda W, Assian E, Ahmad S: Insulin inhibits intranuclear nuclear factor kappa B and stimulates I kappa B in mononuclear cells in obese subjects: Evidence for an anti-inflammatory effect?. J Clin Endocrinol Metab. 2001, 86: 3257-3265. 10.1210/jc.86.7.3257.
Vallabhapurapu S, Karin M: Regulation and function of NF-kappa B transcription factors in the immune system. Annu Rev Immunol. 2009, 27: 693-733. 10.1146/annurev.immunol.021908.132641.
Krieger MJB, Jahan N, Riehle MA, Cao C, Brown MR: Molecular characterization of insulin-like peptide genes and their expression in the African malaria mosquito, Anopheles gambiae. Insect Mol Biol. 2004, 13: 305-315. 10.1111/j.0962-1075.2004.00489.x.
Zimmerman JM, Eliezer N, Simha R: Characterization of amino acid sequences in proteins by statistical methods. J Theor Biol. 1968, 21: 170-201. 10.1016/0022-5193(68)90069-6.
Riley RM, Jin W, Gibson G: Contrasting selection pressures on components of the Ras-mediated signal transduction pathway in Drosophila. Mol Ecol. 2003, 12: 1315-1323. 10.1046/j.1365-294X.2003.01741.x.
Nijhout HF, Berg AM, Gibson WT: A mechanistic study of evolvability using the mitogen-activated protein kinase cascade. Evol Dev. 2003, 5: 281-294. 10.1046/j.1525-142X.2003.03035.x.
Alvarez-Ponce D, Aguade M, Rozas J: Network-level molecular evolutionary analysis of the insulin/TOR signal transduction pathway across 12 Drosophila genomes. Genome Res. 2009, 19: 234-242. 10.1101/gr.084038.108.
We thank Sekou F. Traore and the Malaria Research and Training Center, Faculty of Medicine, University of Mali, Bamako, Mali for support for mosquito collections; Allen van Deynze of the Seed Biotechnology Center, UC Davis, for assistance with design of SNP discovery; Charles M. Nicolet, Director, DNA Technologies Core Facility, UC Davis, for support of genotyping analyses; Kong Wai Cheung and Susan House for assistance with primer design and SNP discovery. Funding for this work was provided by NIH 1R01AI078183-01A2 to SL and GCL; NIH T32 AI074550 Biology of Disease Vectors Training Fellowship to AAH; NIH Fogarty International Center D43 TW007390-01; Jastro-Shields Fellowship and Hazeltine Fellowship to AAH; and a UC Davis Genome Center Core Facility Pilot Project grant. The work at UC Davis was conducted in a facility constructed with support from Research Facilities Improvement Program Grant Number C06 RR-12088-01 from the NIH National Center for Research Resources.
The authors declare that they have no competing interests.
AAH performed the SNP discovery, assisted with the Luminex genotyping, and prepared the manuscript; YL performed the statistical analyses and edited the manuscript; CAC performed the P. falciparum CSP ELISA and assisted with the SNP discovery; VKR performed the Luminex genotyping and edited the manuscript; AJC planned and directed the mosquito collections and edited the manuscript; GCL assisted with the design of the studies, planned the collections and edited the manuscript; SL assisted with the design of the studies and experimental plans, and helped draft the manuscript. All authors read and approved the final manuscript.