Anopheles salivary gland proteomes from major malaria vectors
© Fontaine et al.; licensee BioMed Central Ltd. 2012
Received: 3 March 2012
Accepted: 29 October 2012
Published: 13 November 2012
Antibody responses against Anopheles salivary proteins can indicate individual exposure to bites of malaria vectors. The extent to which these salivary proteins are species-specific is not entirely resolved. Thus, a better knowledge of the diversity among salivary protein repertoires from various malaria vector species is necessary to select relevant genus-, subgenus- and/or species-specific salivary antigens. Such antigens could be used for quantitative (mosquito density) and qualitative (mosquito species) immunological evaluation of malaria vectors/host contact. In this study, salivary gland protein repertoires (sialomes) from several Anopheles species were compared using in silico analysis and proteomics. The antigenic diversity of salivary gland proteins among different Anopheles species was also examined.
In silico analysis of secreted salivary gland protein sequences retrieved from an NCBInr database of six Anopheles species belonging to the Cellia subgenus (An. gambiae, An. arabiensis, An. stephensi and An. funestus) and Nyssorhynchus subgenus (An. albimanus and An. darlingi) displayed a higher degree of similarity compared to salivary proteins from closely related Anopheles species. Additionally, computational hierarchical clustering allowed identification of genus-, subgenus- and species-specific salivary proteins. Proteomic and immunoblot analyses performed on salivary gland extracts from four Anopheles species (An. gambiae, An. arabiensis, An. stephensi and An. albimanus) indicated that heterogeneity of the salivary proteome and antigenic proteins was lower among closely related anopheline species and increased with phylogenetic distance.
This is the first report on the diversity of the salivary protein repertoire among species from the Anopheles genus at the protein level. This work demonstrates that a molecular diversity is exhibited among salivary proteins from closely related species despite their common pharmacological activities. The involvement of these proteins as antigenic candidates for genus-, subgenus- or species-specific immunological evaluation of individual exposure to Anopheles bites is discussed.
Several mosquito species of the Anopheles genus are vectors of Plasmodium parasites, the causal agents of malaria. This major vector-borne disease affects around 216 million individuals annually and leads to more than 600,000 deaths, mainly in tropical and sub-tropical countries . Among approximately 470 Anopheles species indexed worldwide [2, 3], 34 species found in different regions around the world are considered to be the main vectors of the four Plasmodium parasite species (P. falciparum, P. vivax, P. ovale and P. malariae) responsible for human malaria . An. funestus and two sister-species of the An. gambiae sensu lato (s.l.) species complex (i.e., An. gambiae and An. arabiensis) are primary vectors of P. falciparum malaria in sub-Saharan Africa , where 80% of malaria mortality and morbidity occur . Among other anopheline vectors of medical importance, An. stephensi plays a prominent role in urban malaria transmission in the Indo-Pakistan subcontinent [7, 8], and both An. albimanus and An. darlingi are primary vectors of malaria in Central America and various areas of South America [9–11].
In the absence of a licensed malaria vaccine [12, 13] and while Plasmodium drug resistance spreads across the world , vector control is still the most effective method to protect people from arthropod-borne diseases . Prevention of arthropod infective bites can be achieved by personal protective measures and vector control strategies [15, 16]. The evaluation of the effectiveness of these anti-vectorial measures is essentially based on entomological methods such as measuring the density of a mosquito species relative to human density . Catching human landing mosquitoes is currently the most reliable method to estimate host/vector contact [18, 19]. Entomological parameters are also a component of numerous indices used to monitor malaria transmission [17, 20, 21]. However, entomological methods for evaluating the risk of malaria transmission are mainly applicable to the population level and are poorly efficient at evaluating heterogeneity in exposure to vector bites among individuals due to considerable variation of exposure within small geographic areas [22, 23] or heterogeneity in socioeconomic and demographic factors (e.g. age of humans). Furthermore, the human landing catch method is labour-intensive, has budgetary and logistical constraints and is hampered by ethical limitations with the deliberate exposure of individuals to mosquito-borne pathogens. Thus, alternative cost-effective and convenient methods need to be developed to assess human exposure to malaria vectors. During their blood meal, mosquitoes inject saliva into the host’s skin. This saliva contains a cocktail of active components that counteract host haemostasis and modulate immune responses to ensure blood meal success [24, 25]. Secreted salivary proteins of mosquitoes have been reported to elicit antibody responses in people living in endemic areas [26–29] and among travellers transiently exposed to vector bites in tropical areas . These antibody responses were described as being short lived and linked to the level of exposure [28, 30, 31], highlighting the potential use of these responses to arthropod saliva antigens as immunological markers to evaluate individual exposure to arthropod bites  or assess the impact of vector control interventions . Several studies demonstrated the presence of cross-reactive antibody responses against salivary proteins from different hematophagous arthropod species [34–38]. This cross-reactivity was attributed to the existence of antigens shared among different vector species . Species-specific antibody responses against salivary proteins from arthropods have repeatedly been reported [40–42]. Thus, variable levels of homology between salivary protein sequences from different hematophagous arthropods could determine their specificity or cross-recognition . Recently, availability of the genome sequence of several arthropods of major health importance [44, 45], combined with transcriptomic and proteomic analyses of their salivary gland extracts (SGEs) [46–54] have provided new insight into the diversity of salivary molecules among various hematophagous arthropods . These studies revealed a number of secreted protein families, potentially involved in haematophagy or sugar digestion, that are ubiquitous in the Nematocera suborder. Completion of Culex quinquefasciatus, Aedes aegypti and Anopheles gambiae genome sequences also led to the discovery of genus-specific salivary proteins . For Plasmodium spp. vectors, salivary gland transcriptomes and proteomes of An. gambiae (Cellia sb.) [47, 56], An. stephensi (Cellia sb.) , An. funestus (Cellia sb.)  and An. darlingi (Nyssorhynchus sb.) [48, 49] have been examined to date, providing a thorough description of the salivary protein repertoire from these main malaria vectors throughout the world. Notably, secreted salivary proteins were found to be more divergent than housekeeping proteins, indicating a rapid evolution of these proteins within the Anopheles genus [48, 49, 52]. However, sialome diversity among these different anopheline species is poorly documented at the molecular and antigenic levels.
Six Anopheles species were selected according to their major role in human malaria parasite transmission (i.e., An. gambiae, An. arabiensis, An. stephensi, An. funestus, An. albimanus and An. darlingi). Selection of these species was also motivated by the various degrees of phylogenetic relationships among them (i.e., species from the same genus belonging to different subgenera and species complex) and access to their sialomes via salivary gland dissection or through protein sequences obtained by conceptual translation of mRNA sequences previously identified in sialotranscriptomic studies. Importantly, few protein sequences are available for An. arabiensis and An. albimanus due to the lack of sialotranscriptomic studies conducted on these species. Recently, assembly of transcriptional sequences derived from several body tissues including salivary glands of adult female An. albimanus was performed . However, merging the sequence data from the different tissues into a single assembly did not allow clustering secreted salivary protein sequences from the others mosquito body parts, thereby restricting the number of available salivary protein sequences for this species. Salivary gland protein (SGP) repertoires of these different Anopheles species were compared in the present study using in silico analysis and proteomics approaches to assess their diversity at the molecular and protein levels. Conceptual secreted salivary gland protein sequences retrieved from an NCBInr database of six Anopheles species were clustered according to their level of amino acid identity to identify both conserved protein families throughout the Anopheles genus and sub-genus- or species- specific salivary proteins. Proteins contained in salivary gland extracts from four Anopheles species were separated by 1-D SDS-PAGE and identified by tandem mass spectrometry (MS/MS). The antigenic diversity of SGPs was also examined by immunoblot analysis. Collectively, these data represent the first report of genus-, subgenus- and species-specific Anopheles secreted salivary proteins. These proteins could be used for immunological evaluation of the exposure to Anopheles bites.
Phylogenetic relationships between selected anopheline species
Comparison of secreted salivary protein sequences retrieved from public databases
A total of 401 salivary gland protein sequences from the six Anopheles species were retrieved from public databases according to their annotation. Focusing on proteins potentially injected into the human host during mosquito blood feeding, protein sequences were sorted based on signal peptide predictions [64, 65]. A total of 272 out of these 401 salivary proteins, heterogeneously distributed among the six Anopheles species (i.e., 71, 5, 44, 5, 117 and 30 protein sequences for An. gambiae, An. arabiensis, An. stephensi, An. albimanus, An. darlingi and An. funestus, respectively), were predicted to harbor a secretory-signal peptide following submission to SignalP server 3.0 and were thus retained for further analysis. A pairwise protein-protein sequence comparison was first performed on these selected SGPs using “BLAST 2 Sequences”  to find regions of local similarity between sequences from different Anopheles taxa. Briefly, secreted salivary protein sequences from each Anopheles species were gathered into six different databases in accordance with their species affiliation. Each protein from a database was then searched against all proteins from other database in a pairwise fashion. BLAST E-values were used as a parsing criterion in order to select best matches between two different salivary protein repertoires. In order to perform comparative analysis between salivary proteins repertoires, Raw BLAST score obtained from the match of a query protein sequence with a targeted protein sequence were divided against raw self-BLAST score from the match of the query protein sequence to itself to obtain normalised BLAST scores (or BLAST Score Ratio). Normalised BLAST scores range from 0 (no BLAST match) to 1 (perfect BLAST match between two salivary proteins) [66, 67]. The average normalised BLAST scores and average percentage identity between local alignments that estimate the degree of homology between salivary protein sequences from each pair of Anopheles species are presented in Figure 1C and Additional file 1. Only protein sequences from each Anopheles species matching at least one other salivary protein in another species at 40% identity (q.v. Additional file 2) were considered in this analysis, representing 50, 5, 33, 5, 53 and 27 protein sequences for An. gambiae, An. arabiensis, An. stephensi, An. albimanus, An. darlingi and An. funestus, respectively. The average normalised BLAST score between salivary protein sequences from An. gambiae and An. arabiensis was 0.98 ± 0.01 (mean ± SD). Based on the same criteria, salivary protein sequences from both species from the Nyssorhynchus subgenus had normalised BLAST scores of 0.75 ± 0.17 and sequences from species belonging to the Cellia subgenus had normalised BLAST scores larger than 0.45. Lower normalised BLAST scores were observed when comparing salivary protein sequences from species belonging to Cellia with those of the Nyssorhynchus subgenera (all normalised BLAST scores were inferior to 0.40) (Figure 1C, Additional file 1). Thus, secreted salivary protein sequence similarities were the highest among closely related anopheline species and decreased with increasing phylogenetic distance (Figure 1A, B).
Hierarchical clustering of secreted protein sequences
Analysis of salivary gland protein repertoires from four Anopheles species
Antigenic heterogeneity of Anopheles salivary gland proteins
Reduction in exposure to malaria vectors either by controlling Anopheles density or avoiding their bites remains one of the most effective methods to protect human individuals from Plasmodium infections . However, to determine the effectiveness of individual or collective anti-vector measures, it would be useful to develop new indicators that can measure the kinetic variations of individual exposure to specific malaria vectors within a population. By eliciting an antibody response linked to the level of exposure, antibody responses to salivary antigenic proteins from malaria vectors were proposed as valuable immunological markers to estimate host/vector contact [26, 28–30, 72]. Thus, the evaluation of the molecular diversity of salivary proteins from different Anopheles species can be used as a proxy to select genus-, subgenus- and species-specific salivary candidates for subsequent evaluation as immunological markers for Anopheles exposure . The present study assessed the level of molecular and antigenic diversity and relatedness of secreted salivary proteins from major malaria vector species.
Diversity of salivary protein among different species of the Anopheles genus
Vector species of the Anopheles genus throughout the world have different levels of phylogenetic relatedness . Variation in the mitochondrial cytochrome oxidase subunit 2 (COII) sequence has been widely used to display the phylogenetic relationships and population structure of anopheline mosquitoes [73, 74]. Molecular phylogenetic analysis using COII protein sequences from six selected malaria vectors (An. gambiae, An. arabiensis, An. stephensi, An. funestus, An. albimanus and An. darlingi) indicated that genetic distances between species were in agreement with their taxonomic classification. In addition, taking into account all of the secreted salivary proteins available from public databases for these six mosquito species, a pairwise protein-protein sequence analysis demonstrated that the proportion of salivary protein sequence homology decreases according to increasing phylogenetic distance between Anopheles species. Similar observations were reported for salivary proteins from phlebotomine sandflies , underlining this diversity of salivary secreted proteins inside a hematophagous arthropod genus. Comparative analyses of sialotranscriptomes between Anopheles species highlighted that secreted salivary proteins are less conserved than housekeeping proteins [48, 49, 52, 57]. Heterogeneous secreted salivary protein repertoires among Anopheles species are consistent with the existence of secreted salivary proteins occurring throughout the Anopheles genus and others limited to the subgenus or species level.
Genus-specific anopheline secreted proteins
Despite the low number of sequences available for some Anopheles species, several salivary protein families were found in all Anopheles species tested using agglomerative hierarchical clustering. These protein families included the apyrase/5’nucleotidase, antigen 5/gvag, GE-rich/30 kDa, long and short form D7, mucin/13.5 kDa, SG3, SG7, SG10 or hypothetical 6.2 kDa protein families. At the protein level, the clustering of salivary secreted proteins identified by MS confirmed that some proteins, including proteins from the apyrase/5’nucleotidase and 30 kDa/GE-rich/anti-platelet families, are present in all the investigated Anopheles species. Several of these protein families were also found in the saliva of other hematophagous arthropods, which could result from convergent evolution or a common ancestry [32, 76]. For instance, members of the apyrase/5’nucleotidase and antigen 5/gvag protein families are found in the salivary glands of diverse hematophagous insect and tick species across the Arthropoda phylum. Interestingly, members from the GE-rich/30 kDa protein family are exclusively found in salivary glands of both culicine and anopheline female mosquitoes [49, 57, 77]. Here, other protein families appeared to be exclusively found in saliva from anopheline mosquitoes, such as SG3, SG7 or hypothetical 6.2 kDa proteins, offering the opportunity to use these proteins as genus-specific immunological markers. However, cross-reactivity is likely to occur when antigenic proteins share more than 70% amino acid identity . Most of these orthologous protein sequences belonging to the Anopheles genus mostly shared less than 70% identity, minimising the probability of characterising conserved epitopes inside the Anopheles genus [43, 79].
Although an antigenic band was commonly detected at 54 kDa in all Anopheles species, MS did not identify common or orthologous secreted protein in these antigenic bands that could explain the observed cross-reactivity. The low protein abundance in the corresponding gel bands and incomplete molecular sequencing data for some of the Anopheles species could explain this lack of identification. However, members of the GE-rich/30 kDa/anti-platelet family were identified in antigenic bands from the four Anopheles species, and are thus potential candidates to serve as pan-Anopheles genus markers of immunological exposure. Complementary experiments are required to evaluate the lack of cross-reactivity against salivary proteins from other mosquito species from the Culicidae family. Interestingly, in cases of human allergic reactions involving the production of IgE antibodies in response to mosquito bites, some salivary allergens from Aedes aegypti mosquito species have been characterized [34, 80]. Among Aedes aegypti salivary allergens tested, recombinant forms of the 68 kDa salivary apyrase, the 37-kDa protein belonging to the D7 family and the 30 kDa salivary gland allergen were demonstrated to elicit an IgE responses in mosquito-allergic individuals. These data underline the antigenicity of some salivary proteins including 30 kDa family members, and the opportunity to use such proteins for the diagnosis and the desensitisation of mosquito allergic individuals.
Subgenus-specific anopheline secreted salivary proteins
Combinations of subgenus-specific immunological markers may be an alternative for assessing exposure to Anopheles bites. Moreover, such markers could be helpful in identifying the anopheline subgenus involved in malaria transmission in places where several vector subgenera with different behaviour and vectorial capacity are encountered within the same area, notably in Asia [81, 82]. In silico analysis of secreted salivary protein sequences from the six Anopheles species tested in the present study indicated that secreted salivary proteins sharing between 70% and 90% amino acid sequence identity are mainly clustered according to their subgenus affiliation. Thus, several salivary secreted proteins, such as TRIO, gSG2, gSG6 proteins and numerous hypothetical proteins with low molecular weights (Hypothetical 4.2/13, 8.2, 10, 12, or 15/17 kDa proteins), clustered only in the Cellia subgenus. The same observation was made for the Nyssorhynchus subgenus (e.g., salivary peroxidase). Some of the secreted salivary proteins identified by MS were uniquely detected in the Anopheles species belonging to the Cellia subgenus, such as salivary apyrase, anophensin, D7 proteins or TRIO proteins. In silico analysis revealed that members of the gSG2 and gSG6 protein families are well conserved inside the Cellia subgenus with at least 67% and 77% identity, respectively, among An. gambiae, An. funestus and An. stephensi. Among these two protein families, An. gambiae gSG6, which is a reliable marker for exposure to An. gambiae[83–85] bites, was recently reported to be a good indicator for exposure to bites from three main African malaria vectors from Cellia sb. (i.e., An. gambiae, An. arabiensis and An. funestus). This cross-recognition could result from shared epitopes among these orthologous SG6 proteins [38, 86]. According to their low molecular weights (approximately 11.7 and 13 kDa, respectively), gSG2 and gSG6 proteins should be contained in the large 11 kDa antigenic band detected in the Cellia subgenus. However, MS did not identify these last two proteins. Description of the salivary gland protein repertoire from An. gambiae performed by Kalume and colleagues identified gSG6 proteins only using a gel-free approach . Thus, unsuccessful identification of these small proteins could be attributed to the combination of a high number of salivary proteins of this molecular weight and a low number of MS spectra generated by these small proteins, rendering this complex protein mixture unidentifiable by MS. The application of gel-free proteomic methods might increase salivary proteome coverage, especially concerning secreted protein with low molecular weights . As orthologous secreted salivary proteins belonging to the same subgenus clustered only at important sequence identity levels, shared-epitopes are likely to occur, increasing the likelihood of observing cross-reactivity among these subgenus-specific proteins. Thus, identification of apyrase, anophensin and TRIO orthologs in antigenic bands from An. gambiae s.l. and An. stephensi SGEs highlight the potential of these salivary proteins to be Cellia subgenus-specific immunological markers.
Species complex-specific anopheline secreted salivary proteins
The Anopheles subgenera encompass several groups of closely related species that are morphologically indistinguishable, such as the An. gambiae s.l. species complex, which includes at least 6 species [89–91]. In silico analysis revealed a high degree of homology (> 90%) among salivary protein sequences from each selected mosquito of this complex (i.e., An. gambiae and An. arabiensis). Moreover, protein profiles and protein repertoires were highly similar between these two closely related Anopheles species, although one was collected in the field (i.e., An. arabiensis) and the other came from continuous laboratory rearing (i.e., An. gambiae). These results point to the likelihood that salivary protein candidates from either of these species would be able to assess exposure to Anopheles mosquitoes pertaining to the An. gambiae s.l. species complex. Using pooled sera from individuals mainly exposed to An. gambiae, An. arabiensis and An. funestus, most intense antigenic bands (40, 35 and 11 kDa) were revealed in SGEs from both species, suggesting that strong antibody responses against An. gambiae s.l. SGEs are elicited following exposure to these mosquitoes. Indeed, in a recent study conducted in the South of France, a positive association between mosquito exposure and the level of antibody response was reported and this immunological response appeared to be species-specific . The major antigenic bands observed in An. stephensi and An. albimanus SGEs were about 4-fold less intense than the most intense antigenic bands detected in An. gambiae s.l. SGEs. Potential common antigenic salivary proteins were identified in major antigenic bands from An. gambiae and An. arabiensis salivary gland protein profiles, including anti-platelets, anophensin and proteins from the D7 and SG1 families. Even if members of these protein families are present throughout the Cellia subgenus, low amino acid sequence identities can occur between An. gambiae s.l. and other species, which could drastically affect antibody binding. Cross-reactivity usually implies a lower affinity for the cross-reactive antigen compared to the primary antigen . The development of better quantitative methods, such as ELISA or Luminex®, combined with the production of salivary antigenic protein candidates may provide more distinct antibody responses to specific mosquito bites from the cross-reactivity attributed to partial shared-antigens.
Species-specific anopheline secreted salivary proteins
It would be interesting to use immunological tools to assess individual exposure to a specific Anopheles species, even to closely related phylogenetically species, in areas harbouring several Anopheles species with distinct vector competences and behaviours. For instance, the An. gambiae s.l. complex encompasses species that often occur in sympatry but profoundly differ in their ability to transmit malaria parasites, in host feeding preferences, larval habitat requirements or responses to select vector control measures [90, 93]. Thus, characterisation of Anopheles species-specific immunological markers may help to implement adapted vector-control strategies and assess their efficiency to prevent host/vector contact at the individual level. Additionally, the combination of different species-specific immunological markers using multiplex techniques, such as Luminex® , may increase sensitivity [38, 95] and specificity of the test. Such multiplexing of validated antigenic salivary proteins will distinguish exposure to the bite of malaria vectors from that of non-vector Anopheles. Species-specific immunological markers may provide a more detailed view of the history of exposure at an individual level in retrospective studies. Thus, these markers may be useful to determine the implication of different vector species in malaria epidemics. In silico analysis indicated that 73 (27%) salivary proteins have no orthologous proteins in other Anopheles species at the lowest amino acid sequence identity threshold tested (40% identity). In addition, MS only identified 18 secreted in unique Anopheles species. Some of these species-specific salivary proteins, including mainly proteins with low molecular weights, are not found in SGEs from other hematophagous arthropod species (Additional file 2). Thus, these can serve as potential immunological markers for the assessment of individual exposure to specific Anopheles species.
Although in silico analysis, protein patterns and protein repertoires indicated a low diversity of secreted salivary proteins between An. gambiae and An. arabiensis, the detection of three antigenic bands exclusively in An. arabiensis SGEs suggest that closely related species could be distinct at the antigenic level. Unfortunately, MS analysis could not successfully identify the corresponding proteins from the gel bands. To identify these antigenic bands, a two-dimensional immunoproteomic approach using a fluorescence-based method could be an alternative . Although the challenge appears more important for closely related mosquito species, our data suggest that identification of species-specific immunological markers seems reasonably conceivable. Moreover, this hypothesis is supported by our recent work on the existence of species-specific serological responses against Ae. caspius SGEs .
The present study assessed, for the first time, the sialome diversity among different Anopheles species at the molecular and antigenic levels by combining in silico, proteomic and immuno-proteomic approaches. Our results demonstrate that salivary protein sequence identities among different species from one Culicidae genus are heterogeneous, with salivary proteins present throughout the Anopheles genus, or specific at the subgenus or species levels. This result demonstrates that salivary proteins from closely related species exhibit molecular diversity despite their common pharmacological activities (e.g. anti-haemostatic or immunomodulatory activities). This work supports the idea that genus-, subgenus- and species-specific salivary proteins can be used to develop immunological markers of individual exposure to malaria vectors. In complement to entomological methods, such immunological markers of exposure may be useful in the evaluation of anti-vector intervention strategies. In addition, development of species-specific immunological markers may help determine the implication of different vector species in malaria epidemics and provide further understanding of the vectorial system in a given area.
Sera from 5 individuals living in the Senegalese village of Dielmo (13°45’N, 16°25’W), sampled in March 1995 when malaria was holoendemic were used in this study . These individuals were exposed to high levels of malaria transmission (about 200 infective bites/person/per year) and mosquito bites (human biting rate about 23.2), with An. gambiae (11%), An. arabiensis (56%) and An. funestus (33%) as principal vectors . These individuals did not travel outside Senegal country in the twelve months prior to blood sampling. The protocol was approved by the Senegal National Ethics Committee (Dakar, Senegal). The informed consent of each participant was obtained at the beginning of the study, after a thorough explanation of its purpose.
The cytochrome oxidase II (COII) protein sequences from An. gambiae, An. arabiensis, An. stephensi, An. albimanus, An. funestus, An. darlingi and Aedes aegypti were retrieved from NCBInr database (May 10th, 2011) and multiple sequence alignment was performed with Clustal W 1.7 multiple sequence alignment program  which is included in Molecular Evolutionary genetic Analysis 5 (MEGA 5) programs package . The evolutionary history was inferred using the Neighbour-Joining method in MEGA 5 . The Ae. aegypti sequence was taken as an outgroup. The bootstrap consensus tree inferred from 10,000 replicates  is taken to represent the evolutionary history of the taxa analyzed. Branches corresponding to partitions reproduced in less than 50% bootstrap replicates are collapsed. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (10,000 replicates) are shown next to the branches . The evolutionary distances were computed using the Poisson correction method and are in the units of the number of amino acid substitutions per site. All positions containing gaps and missing data were eliminated. Evolutionary analyses were conducted in MEGA5 . Protein sequence alignments were generated in Geneious Pro 5.6.4 http://www.geneious.com/.
In silico analysis
Sequences retrieval and pairwise protein sequence comparison
All salivary protein sequences were retrieved in FASTA format from the online non-redundant National Center for Biotechnology Information protein database (NCBInr, NIH, Bethesda, June 15th, 2011) under the taxonomies An. gambiae , An. arabiensis , An. stephensi , An. funestus , An. darlingi  and An. albimanus , using the search term “salivary” in any fields of their description text. Signal peptides were predicted by submission of the protein sequences to the SignalP server 3.0 , allowing the determination of putative secreted proteins. Putative secreted proteins sequences from each Anopheles species were blasted again each other using “BLAST 2 Sequences”  with default parameters (NCBI, National Library of Medicine, http://blast.ncbi.nlm.nih.gov/Blast.cgi). For each protein from species A, the raw BLAST score and percentage identity of the best match with species B was manually selected according to the E-value. Duplicates of protein sequences from species B that matched several proteins in species A were then removed according to their lowest E-value in order to select unique and best homologous protein sequence between two species. The BLAST Score Ratio approach was adopted to represent similarities between salivary proteins pertaining to different Anopheles species [66, 67]. In this approach, the raw BLAST score resulting from a comparison between a query saliva protein sequence from species A and a protein sequence from species B is divided by the self-BLAST score obtained with the BLAST of the query protein sequence from species A with itself. The resulting normalized BLAST score vary from 0 (no match) to 1 (perfect match). The use of such normalized scores overcomes several problems associated with the use of E-values, such as biases entailed in comparisons among different databases, falsely high E-values assigned to low-complexity proteins and falsely low E-values based on small regions of high similarities [66, 67]. The average normalized BLAST score as well as average percentage identity between salivary proteins from two species were then calculated.
All putative secreted sequences were merged in a single FASTA file and submitted to the CD-HIT server  for hierarchical clustering (H-CD-HIT) as describe elsewhere . Briefly, the program performs clustering three times in succession with decreasing similarity thresholds. First, clustering start with the input dataset at a high identity threshold (≥ 90%). The longest sequence becomes the representative of the first cluster. Then, each remaining sequence is compared to the representatives of all existing clusters. If the predefined similarity threshold is met, the sequence is grouped into the most similar cluster. Otherwise, a new cluster is defined with that sequence as the representative. The last two steps of the hierarchical clustering (≥ 70% and ≥ 40% similarity threshold) start with representatives of the previous clustering runs and the whole process produces a hierarchical structure. The percentage identities are calculated by counting the numbers of identical amino acids between two protein sequences by using a short word filter (For details see the user’s guide at the following web link (http://www.bioinformatics.org/cd-hit/cd-hit-user-guide.pdf). Each salivary protein sequence was further blasted against a database containing SGP sequences identified so far in 26 hematophagous arthropod species (i.e., An. gambiae , An. arabiensis , An. stephensi , An. funestus , An. albimanus , An. darlingi , Ae. aegypti , Ae. albopictus , Ochlerotatus triseriatus , Culex tarsalis , Cx quinquefasciatus , Triatoma brasiliensis , T. infestans , Rhodnius prolixus , Cimex lectularius , Glossina morsitans morsitans , Xenopsylla cheopis , Stomoxys calcitrans , Simulium vittatum , S. nigrimanum , Ornithodoros parkeri , O. coriaceus , Argas monolakensis , Ixodes pacificus , I. scapularis  and Rhipicephalus sanguineus ) by sialotranscriptomic studies . Proteins matching with an E-value < 1×10-10 were considered putative homologs .
Mosquitoes and salivary gland extraction
Uninfected 5-day-old adult females of the An. gambiae s.s., An. arabiensis, An. stephensi and An. albimanus species were used in this study. An. stephensi and An. albimanus species were reared at the Institut Pasteur (CEntre de Production et Infection des Anopheles, CEPIA, Paris). An. gambiae was reared at the Institut de Recherche pour le Développement (laboratoire de Lutte contre les Insectes Nuisibles, Montpellier). An. arabiensis species was collected on the field at the larvae stage at Dakar (Senegal) at the end of the rainy season in September 2008  and identified by polymerase chain reaction (PCR) . After emergence, adult mosquitoes were reared at the Institut de Recherche pour le Développement (UR24, Dakar). All mosquitoes used in the experiments were maintained under identical standard conditions: 26°C and 60% humidity, took no blood meals and were maintained on a diet of 10% syrup solution. The salivary glands from adult mosquito females were dissected under a stereomicroscope at 4X magnification as previously described . The salivary glands from each experiment were pooled by strains into a microcentrifuge tube on ice in phosphate-buffered saline (PBS) and then stored frozen at −20°C until needed.
Salivary glands were disrupted by ultrasonication (Vibracell 72412, Bioblock Scientific, Illkirch, France) for 5 min on ice at maximum amplitude. Salivary gland homogenates were then centrifuged for 15 min at 16,100 ×g  and protein concentration of the supernatant was determined in duplicate by the Lowry method (DC Protein assay Kit, Bio-Rad) according to the manufacturer’s instructions. Salivary gland proteins were then concentrated by precipitation with acetone (Sigma, St Louis, MI), and were suspended in a UTC buffer containing 8M urea (Sigma), 2M thiourea (Sigma), 4% (w/v) CHAPS (Sigma) and 30 mM Tris (Sigma), adjusted to pH 8.5 in order to obtain a protein concentration adjusted to 2.5μg/μL.
One-dimensional electrophoresis (SDS-PAGE)
For each species, 20 μg of salivary gland proteins were separated onto a 12% SDS-PAGE in a PROTEAN II xi or Mini (BioRad, Hercules, USA). A broad range molecular weight marker (BioRad) was loaded on each gel. Gels were stained with Sypro Ruby (BioRad) according to the manufacturer’s protocol and digitalized using the TyphoonTM Trio Image scanner (GE Healthcare, UK). Salivary glands densitometry profiles were analyzed using the ImageQuantTM TL software (GE Healthcare, UK). Background subtraction was performed and the densitometry profiles were normalized to take into account global differences .
For each species, 15 μg of SGPs were separated onto a 12% SDS-PAGE as described above. SGPs were minimally labelled with CyDye 5 according to the manufacturer’s protocol (GE Healthcare) prior electrophoresis, as described elsewhere . Gels were then transferred to a nitrocellulose membrane (0.45-μm, Amersham Pharmacia, Saclay, France) by semidry blotting (0.8 mA per cm2) . Blots were saturated 1 hour at room temperature with 5% w/v non-fat dried milk, and were carried out with human sera diluted at 1/50 in phosphate buffer saline (PBS) containing 0.1% v/v tween-20 with 5% w/v non-fat dried milk. After an overnight incubation, blots were incubated with anti-human Fcγ/IgG FITC conjugated antibody 1/1000 (Sigma, St Louis, MI). Immunoblots were directly digitalized using a TyphoonTM Trio Image scanner (GE Healthcare) and densitometric analysis of IgG immune profiles was performed using ImageQuantTM TL software (GE Healthcare), as previously described .
Protein band excision and in-gel digestion
Each loading tracks were excised covering the totality of the each lane as previously described [54, 68], using Shimadzu Biotech Xcise System (Champs sur Marne, France). Salivary protein identification was made in duplicate on two distinct gels. Protein bands were digested overnight at 37°C with sequencing-grade trypsin (12.5 μg/mL; Promega Madison, WI, USA) in 50 mM NH4HCO3 (Sigma). The resulting peptides were extracted with 25 mM NH4HCO3 for 15 min, dehydrated with acetonitrile (ACN) (Sigma), incubated with 5% acid formic (Sigma) for 15 min under agitation, then dehydrated with ACN, and finally completely dried using a SpeedVac. Samples were then stored at −20°C before analysis by mass spectrometry (MS).
Mass spectrometry analysis
For MS analysis, a LCQ DecaXPplus (ThermoFinnigan, San Jose, CA) ion trap was used. Nano-liquid separation of peptides was carried out using an Ettan MDLC chromatographic system (GE Healthcare) in high throughput configuration. Ten microliters of the digest were first trapped on a zorbax 300SB-C18 5 × 0.3 mm column and eluted at a flow rate of approximately 200 nl/min on a zorbax 300SB-C18, 3.5 μm, 150 × 0.075 mm by a linear gradient of eluant B (0.1% Formic acid, 84% ACN) in eluant A (1% Formic acid). The chromatographic system was piloted by the Unicorn 5.01 software (GE Healthcare). MS measurements were done on a LCQTM Deca XP Plus ion trap mass spectrometer (ThermoFinnigan) equipped with a LCQTM nanospray ionization source. A spray voltage of 1.8 kV was applied to the liquid junction via an in-union high voltage contact coupled to a silicaTip emitter (New Objective). Operation of the mass spectrometer was fully automated during the entire procedure using the Excalibur 1.3 data system (ThermoFinnigan). Continuous cycles of one full scan (m/z 500 to 1700) followed by three data-dependent MS/MS measurements at 35% normalized collision energy were done. MS/MS measurements were allowed for the three most intense precursor ions with a maximum rejection time limit of 1 min.
MS data analysis
All MS/MS spectra from one species were gathered and sequence database searched using the Bioworks 3.1 (ThermoFinnigan) or Mascot Daemon v.2.2.2 (Matrix Science, London, UK). The MS/MS spectra were searched against the non-redundant protein database (NCBInr, NIH, Bethesda) of An. gambiae , An. arabiensis , An. stephensi  and An. albimanus  together (July 27th, 2009, 16677 sequences). The following search parameters were used: precursor-ion mass tolerance of 0.8 Da, fragment ion tolerance of 0.8 Da with methionine oxidation and cysteine carboxyamidomethylation specified as differential modifications, and a maximum of one missed cleavage site allowed. Scaffold (version Scaffold_3.6.2, Proteome Software Inc., Portland, OR) was used to validate MS/MS based peptide and protein identifications. Peptide identifications were accepted if they could be established at greater than 95.0% probability as specified by the Peptide Prophet algorithm . Protein identifications were accepted if they could be established at greater than 95.0% probability and contained at least 1 identified peptide. Protein probabilities were assigned by the Protein Prophet algorithm .
The authors gratefully acknowledge D. Fontenille and M.N. Lacroix from IRD UR016, Montpellier and C. Bourgouin from CEPIA, Institut Pasteur, Paris for access to mosquitoes bred in insectary. The authors also acknowledge J.F. Trape and C. Bouganali from IRD UR24, Dakar for access to wild mosquito in Senegal and S. Audebert, L. Camoin, M. Belghazi and P. Fourquet for their help in proteomic analysis. We thank F. Pagès for his thoughtful comments and inputs. Authors also acknowledge the teams of IRD and Institut Pasteur de Dakar for the collect of sera and epidemiological and entomological data in Dielmo, Senegal. We are grateful to L. Lambrechts for linguistic improvements. This study was supported by the French Armed Forces Medical Service and the Délégation Générale pour l’Armement (DGA/SSA, ArthroSer project, Grant 10CO401).
- WHO: World malaria report. 2009,http://www.who.int/malaria/world_malaria_report_2009/en/index.html,Google Scholar
- Systematic Catalog of Culicidae.http://www.mosquitocatalog.org/default.aspx,
- Harbach RE: The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships. Bull Entomol Res. 2004, 94 (6): 537-553.PubMedGoogle Scholar
- Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004, 70 (5): 486-498.PubMedGoogle Scholar
- Kelly-Hope LA, Hemingway J, McKenzie FE: Environmental factors associated with the malaria vectors Anopheles gambiae and Anopheles funestus in Kenya. Malar J. 2009, 8: 268-10.1186/1475-2875-8-268.PubMed CentralPubMedGoogle Scholar
- RollBackMalariaPartnership: World Malaria Report. 2005, Geneva: World Health OrganizationGoogle Scholar
- Rowland M, Mohammed N, Rehman H, Hewitt S, Mendis C, Ahmad M, Kamal M, Wirtz R: Anopheline vectors and malaria transmission in eastern Afghanistan. Trans R Soc Trop Med Hyg. 2002, 96 (6): 620-626. 10.1016/S0035-9203(02)90331-7.PubMedGoogle Scholar
- Hati AK: Urban malaria vector biology. Indian J Med Res. 1997, 106: 149-163.PubMedGoogle Scholar
- Gutierrez LA, Naranjo NJ, Cienfuegos AV, Muskus CE, Luckhart S, Conn JE, Correa MM: Population structure analyses and demographic history of the malaria vector Anopheles albimanus from the Caribbean and the Pacific regions of Colombia. Malar J. 2009, 8: 259-10.1186/1475-2875-8-259.PubMed CentralPubMedGoogle Scholar
- Gutierrez LA, Naranjo N, Jaramillo LM, Muskus C, Luckhart S, Conn JE, Correa MM: Natural infectivity of Anopheles species from the Pacific and Atlantic Regions of Colombia. Acta Trop. 2008, 107 (2): 99-105. 10.1016/j.actatropica.2008.04.019.PubMedGoogle Scholar
- Zimmerman RH: Ecology of malaria vectors in the Americas and future direction. Mem Inst Oswaldo Cruz. 1992, 87 (Suppl 3): 371-383.PubMedGoogle Scholar
- Crompton PD, Pierce SK, Miller LH: Advances and challenges in malaria vaccine development. J Clin Invest. 2010, 120 (12): 4168-4178. 10.1172/JCI44423.PubMed CentralPubMedGoogle Scholar
- Rogier C, Orlandi-Pradines E, Fusai T, Pradines B, Briolant S, Almeras L: Malaria vaccines: prospects and reality. Med Mal Infect. 2006, 36 (8): 414-422. 10.1016/j.medmal.2006.05.016.PubMedGoogle Scholar
- Pradines B, Vial H, Olliaro P: Malaria prophylaxis and treatment: problems, recent developments and perspectives. Med Trop (Mars). 2003, 63 (1): 79-98.Google Scholar
- Enayati A, Lines J, Maharaj R, Hemingway J: Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA. San Francisco: The Global Health Group 2009, Chapter 9: Suppressing the Vector, 140-154.Google Scholar
- WHO: Malaria vector control and personal protection. 2006, Geneva, Switzerland: WHO Technical Report, No 936Google Scholar
- Shaukat AM, Breman JG, McKenzie FE: Using the entomological inoculation rate to assess the impact of vector control on malaria parasite transmission and elimination. Malar J. 2010, 9: 122-10.1186/1475-2875-9-122.PubMed CentralPubMedGoogle Scholar
- Mathenge EM, Misiani GO, Oulo DO, Irungu LW, Ndegwa PN, Smith TA, Killeen GF, Knols BG: Comparative performance of the Mbita trap, CDC light trap and the human landing catch in the sampling of Anopheles arabiensis, An. funestus and culicine species in a rice irrigation in western Kenya. Malar J. 2005, 4: 7-10.1186/1475-2875-4-7.PubMed CentralPubMedGoogle Scholar
- Mboera LE: Sampling techniques for adult Afrotropical malaria vectors and their reliability in the estimation of entomological inoculation rate. Tanzan Health Res Bull. 2005, 7 (3): 117-124.PubMedGoogle Scholar
- Smith DL, Smith TA, Hay SI: Shrinking the Malaria map: A prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA. The Global Health Group 2009, Chapter 7: Measuring Malaria for Elimination, San FranciscoGoogle Scholar
- Garrett-Jones C, Shidrawi GR: Malaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomology. Bull World Health Organ. 1969, 40 (4): 531-545.PubMed CentralPubMedGoogle Scholar
- Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, Swalm C, Keating J, Regens JL, Shililu JI, et al: Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003, 68 (6): 734-742.PubMedGoogle Scholar
- Orlandi-Pradines E, Rogier C, Koffi B, Jarjaval F, Bell M, Machault V, Pons C, Girod R, Boutin JP, Pages F: Major variations in malaria exposure of travellers in rural areas: an entomological cohort study in western Cote d'Ivoire. Malar J. 2009, 8: 171-10.1186/1475-2875-8-171.PubMed CentralPubMedGoogle Scholar
- Ribeiro JM: Blood-feeding arthropods: live syringes or invertebrate pharmacologists?. Infect Agents Dis. 1995, 4 (3): 143-152.PubMedGoogle Scholar
- Ribeiro JM, Francischetti IM: Role of arthropod saliva in blood feeding: sialome and post-sialome perspectives. Annu Rev Entomol. 2003, 48: 73-88. 10.1146/annurev.ento.48.060402.102812.PubMedGoogle Scholar
- Cornelie S, Remoue F, Doucoure S, Ndiaye T, Sauvage FX, Boulanger D, Simondon F: An insight into immunogenic salivary proteins of Anopheles gambiae in African children. Malar J. 2007, 6: 75-10.1186/1475-2875-6-75.PubMed CentralPubMedGoogle Scholar
- Remoue F, Alix E, Cornelie S, Sokhna C, Cisse B, Doucoure S, Mouchet F, Boulanger D, Simondon F: IgE and IgG4 antibody responses to Aedes saliva in African children. Acta Trop. 2007, 104 (2–3): 108-115.PubMedGoogle Scholar
- Remoue F, Cisse B, Ba F, Sokhna C, Herve JP, Boulanger D, Simondon F: Evaluation of the antibody response to Anopheles salivary antigens as a potential marker of risk of malaria. Trans R Soc Trop Med Hyg. 2006, 100 (4): 363-370. 10.1016/j.trstmh.2005.06.032.PubMedGoogle Scholar
- Londono-Renteria BL, Eisele TP, Keating J, James MA, Wesson DM: Antibody response against Anopheles albimanus (Diptera: Culicidae) salivary protein as a measure of mosquito bite exposure in Haiti. J Med Entomol. 2010, 47 (6): 1156-1163. 10.1603/ME09240.PubMedGoogle Scholar
- Orlandi-Pradines E, Almeras L, Denis De Senneville L, Barbe S, Remoue F, Villard C, Cornelie S, Penhoat K, Pascual A, Bourgouin C, et al: Antibody response against saliva antigens of Anopheles gambiae and Aedes aegypti in travellers in tropical Africa. Microbes Infect. 2007, 9 (12–13): 1454-1462.PubMedGoogle Scholar
- Palosuo K, Brummer-Korvenkontio H, Mikkola J, Sahi T, Reunala T: Seasonal increase in human IgE and IgG4 antisaliva antibodies to Aedes mosquito bites. Int Arch Allergy Immunol. 1997, 114 (4): 367-372. 10.1159/000237696.PubMedGoogle Scholar
- Fontaine A, Diouf I, Bakkali N, Misse D, Pages F, Fusai T, Rogier C, Almeras L: Implication of haematophagous arthropod salivary proteins in host-vector interactions. Parasit Vectors. 2011, 4: 187-10.1186/1756-3305-4-187.PubMed CentralPubMedGoogle Scholar
- Drame PM, Poinsignon A, Besnard P, Le Mire J, Dos-Santos MA, Sow CS, Cornelie S, Foumane V, Toto JC, Sembene M, et al: Human antibody response to Anopheles gambiae saliva: an immuno-epidemiological biomarker to evaluate the efficacy of insecticide-treated nets in malaria vector control. Am J Trop Med Hyg. 2010, 83 (1): 115-121. 10.4269/ajtmh.2010.09-0684.PubMed CentralPubMedGoogle Scholar
- Peng Z, Xu W, Lam H, Cheng L, James AA, Simons FE: A new recombinant mosquito salivary allergen, rAed a 2: allergenicity, clinical relevance, and cross-reactivity. Allergy. 2006, 61 (4): 485-490. 10.1111/j.1398-9995.2006.00985.x.PubMedGoogle Scholar
- Wheeler CM, Coleman JL, Benach JL: Salivary gland antigens of Ixodes dammini are glycoproteins that have interspecies cross-reactivity. J Parasitol. 1991, 77 (6): 965-973. 10.2307/3282751.PubMedGoogle Scholar
- Brummer-Korvenkontio H, Palosuo T, Francois G, Reunala T: Characterization of Aedes communis, Aedes aegypti and Anopheles stephensi mosquito saliva antigens by immunoblotting. Int Arch Allergy Immunol. 1997, 112 (2): 169-174. 10.1159/000237450.PubMedGoogle Scholar
- Cross ML, Cupp MS, Cupp EW, Ramberg FB, Enriquez FJ: Antibody responses of BALB/c mice to salivary antigens of hematophagous black flies (Diptera: Simuliidae). J Med Entomol. 1993, 30 (4): 725-734.PubMedGoogle Scholar
- Rizzo C, Ronca R, Fiorentino G, Mangano VD, Sirima SB, Nebie I, Petrarca V, Modiano D, Arca B: Wide cross-reactivity between Anopheles gambiae and Anopheles funestus SG6 salivary proteins supports exploitation of gSG6 as a marker of human exposure to major malaria vectors in tropical Africa. Malar J. 2011, 10: 206-10.1186/1475-2875-10-206.PubMed CentralPubMedGoogle Scholar
- Peng Z, Li H, Simons FE: Immunoblot analysis of salivary allergens in 10 mosquito species with worldwide distribution and the human IgE responses to these allergens. J Allergy Clin Immunol. 1998, 101 (4 Pt 1): 498-505.PubMedGoogle Scholar
- Marshall NA, Chapman MD, Saxon A: Species-specific allergens from the salivary glands of Triatominae (Heteroptera:Reduviidae). J Allergy Clin Immunol. 1986, 78 (3 Pt 1): 430-435.PubMedGoogle Scholar
- Volf P, Rohousova I: Species-specific antigens in salivary glands of phlebotomine sandflies. Parasitology. 2001, 122 (Pt 1): 37-41.PubMedGoogle Scholar
- Jeon SH, Park JW, Lee BH: Characterization of human IgE and mouse IgG1 responses to allergens in three mosquito species by immunoblotting and ELISA. Int Arch Allergy Immunol. 2001, 126 (3): 206-212. 10.1159/000049515.PubMedGoogle Scholar
- Aalberse RC, Akkerdaas J, van Ree R: Cross-reactivity of IgE antibodies to allergens. Allergy. 2001, 56 (6): 478-490. 10.1034/j.1398-9995.2001.056006478.x.PubMedGoogle Scholar
- Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, et al: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298 (5591): 129-149. 10.1126/science.1076181.PubMedGoogle Scholar
- Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, et al: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316 (5832): 1718-1723. 10.1126/science.1138878.PubMedGoogle Scholar
- Arca B, Lombardo F, Francischetti IM, Pham VM, Mestres-Simon M, Andersen JF, Ribeiro JM: An insight into the sialome of the adult female mosquito Aedes albopictus. Insect Biochem Mol Biol. 2007, 37 (2): 107-127. 10.1016/j.ibmb.2006.10.007.PubMedGoogle Scholar
- Arca B, Lombardo F, Valenzuela JG, Francischetti IM, Marinotti O, Coluzzi M, Ribeiro JM: An updated catalogue of salivary gland transcripts in the adult female mosquito, Anopheles gambiae. J Exp Biol. 2005, 208 (Pt 20): 3971-3986.PubMedGoogle Scholar
- Calvo E, Andersen J, Francischetti IM, De LCM, DeBianchi AG, James AA, Ribeiro JM, Marinotti O: The transcriptome of adult female Anopheles darlingi salivary glands. Insect Mol Biol. 2004, 13 (1): 73-88. 10.1111/j.1365-2583.2004.00463.x.PubMedGoogle Scholar
- Calvo E, Pham VM, Marinotti O, Andersen JF, Ribeiro JM: The salivary gland transcriptome of the neotropical malaria vector Anopheles darlingi reveals accelerated evolution of genes relevant to hematophagy. BMC Genomics. 2009, 10: 57-10.1186/1471-2164-10-57.PubMed CentralPubMedGoogle Scholar
- Calvo E, Sanchez-Vargas I, Favreau AJ, Barbian KD, Pham VM, Olson KE, Ribeiro JM: An insight into the sialotranscriptome of the West Nile mosquito vector. Culex tarsalis. BMC genomics. 2010, 11: 51-10.1186/1471-2164-11-51.PubMed CentralPubMedGoogle Scholar
- Ribeiro JM, Charlab R, Pham VM, Garfield M, Valenzuela JG: An insight into the salivary transcriptome and proteome of the adult female mosquito Culex pipiens quinquefasciatus. Insect Biochem Mol Biol. 2004, 34 (6): 543-563. 10.1016/j.ibmb.2004.02.008.PubMedGoogle Scholar
- Valenzuela JG, Francischetti IM, Pham VM, Garfield MK, Ribeiro JM: Exploring the salivary gland transcriptome and proteome of the Anopheles stephensi mosquito. Insect Biochem Mol Biol. 2003, 33 (7): 717-732. 10.1016/S0965-1748(03)00067-5.PubMedGoogle Scholar
- Ribeiro JM, Arca B, Lombardo F, Calvo E, Phan VM, Chandra PK, Wikel SK: An annotated catalogue of salivary gland transcripts in the adult female mosquito. Aedes aegypti. BMC genomics. 2007, 8: 6-10.1186/1471-2164-8-6.PubMed CentralPubMedGoogle Scholar
- Almeras L, Fontaine A, Belghazi M, Bourdon S, Boucomont-Chapeaublanc E, Orlandi-Pradines E, Baragatti M, Corre-Catelin N, Reiter P, Pradines B, et al: Salivary gland protein repertoire from Aedes aegypti mosquitoes. Vector Borne Zoonotic Dis. 2010, 10 (4): 391-402. 10.1089/vbz.2009.0042.PubMedGoogle Scholar
- Ribeiro JMC, Mans BJ, Arcà B: An insight into the sialome of blood-feeding Nematocera. Insect Biochem Mol Biol. 2010, 40: 767-784. 10.1016/j.ibmb.2010.08.002.PubMed CentralPubMedGoogle Scholar
- Calvo E, Pham VM, Lombardo F, Arca B, Ribeiro JM: The sialotranscriptome of adult male Anopheles gambiae mosquitoes. Insect Biochem Mol Biol. 2006, 36 (7): 570-575. 10.1016/j.ibmb.2006.04.005.PubMedGoogle Scholar
- Calvo E, Dao A, Pham VM, Ribeiro JM: An insight into the sialome of Anopheles funestus reveals an emerging pattern in anopheline salivary protein families. Insect Biochem Mol Biol. 2007, 37 (2): 164-175. 10.1016/j.ibmb.2006.11.005.PubMed CentralPubMedGoogle Scholar
- Martinez-Barnetche J, Gomez-Barreto RE, Ovilla-Munoz M, Tellez-Sosa J, Garcia-Lopez DE, Dinglasan RR, Ubaida Mohien C, Maccallum RM, Redmond SN, Gibbons JG, et al: Transcriptome of the adult female malaria mosquito vector Anopheles albimanus. BMC Genomics. 2012, 13 (1): 207-10.1186/1471-2164-13-207.PubMed CentralPubMedGoogle Scholar
- Arca B, Lombardo F, De Lara Capurro M, Della Torre A, Dimopoulos G, James AA, Coluzzi M: Trapping cDNAs encoding secreted proteins from the salivary glands of the malaria vector Anopheles gambiae. Proc Natl Acad Sci U S A. 1999, 96 (4): 1516-1521. 10.1073/pnas.96.4.1516.PubMed CentralPubMedGoogle Scholar
- Mohanty A, Swain S, Kar SK, Hazra RK: Analysis of the phylogenetic relationship of Anopheles species, subgenus Cellia (Diptera: Culicidae) and using it to define the relationship of morphologically similar species. Infect Genet Evol. 2009, 9 (6): 1204-1224. 10.1016/j.meegid.2009.06.021.PubMedGoogle Scholar
- Foley DH, Bryan JH, Yeates D, Saul A: Evolution and systematics of Anopheles: insights from a molecular phylogeny of Australasian mosquitoes. Mol Phylogenet Evol. 1998, 9 (2): 262-275. 10.1006/mpev.1997.0457.PubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedGoogle Scholar
- Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174 (2): 247-250. 10.1111/j.1574-6968.1999.tb13575.x.PubMedGoogle Scholar
- Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.PubMedGoogle Scholar
- Nielsen H, Engelbrecht J, Brunak S, von Heijne G: A neural network method for identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Int J Neural Syst. 1997, 8 (5–6): 581-599.PubMedGoogle Scholar
- Serazin AC, Dana AN, Hillenmeyer ME, Lobo NF, Coulibaly MB, Willard MB, Harker BW, Sharakhov IV, Collins FH, Ribeiro JM, et al: Comparative analysis of the global transcriptome of Anopheles funestus from Mali, West Africa. PLoS One. 2009, 4 (11): e7976-10.1371/journal.pone.0007976.PubMed CentralPubMedGoogle Scholar
- Rasko DA, Myers GS, Ravel J: Visualization of comparative genomic analyses by BLAST score ratio. BMC Bioinformatics. 2005, 6: 2-10.1186/1471-2105-6-2.PubMed CentralPubMedGoogle Scholar
- Huang Y, Niu B, Gao Y, Fu L, Li W: CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010, 26 (5): 680-682. 10.1093/bioinformatics/btq003.PubMed CentralPubMedGoogle Scholar
- Li W, Wooley JC, Godzik A: Probing metagenomics by rapid cluster analysis of very large datasets. PLoS One. 2008, 3 (10): e3375-10.1371/journal.pone.0003375.PubMed CentralPubMedGoogle Scholar
- Almeras L, Orlandi-Pradines E, Fontaine A, Villard C, Boucomont E, de Senneville LD, Baragatti M, Pascual A, Pradines B, Corre-Catelin N, et al: Sialome individuality between Aedes aegypti colonies. Vector Borne Zoonotic Dis. 2009, 9 (5): 531-541. 10.1089/vbz.2008.0056.PubMedGoogle Scholar
- Fontenille D, Lochouarn L, Diagne N, Sokhna C, Lemasson JJ, Diatta M, Konate L, Faye F, Rogier C, Trape JF: High annual and seasonal variations in malaria transmission by anophelines and vector species composition in Dielmo, a holoendemic area in Senegal. Am J Trop Med Hyg. 1997, 56 (3): 247-253.PubMedGoogle Scholar
- Waitayakul A, Somsri S, Sattabongkot J, Looareesuwan S, Cui L, Udomsangpetch R: Natural human humoral response to salivary gland proteins of Anopheles mosquitoes in Thailand. Acta Trop. 2006, 98 (1): 66-73. 10.1016/j.actatropica.2006.02.004.PubMedGoogle Scholar
- Watson TM, Kay BH: Vector competence of Aedes notoscriptus (Diptera: Culicidae) for Barmah Forest virus and of this species and Aedes aegypti (Diptera: Culicidae) for dengue 1–4 viruses in Queensland, Australia. J Med Entomol. 1999, 36 (4): 508-514.PubMedGoogle Scholar
- Aitken T: An in vitro feeding technique for artificially demonstrating virus transmission by mosquitoes. Mosq News. 1977, 37: 130-133.Google Scholar
- Anderson JM, Oliveira F, Kamhawi S, Mans BJ, Reynoso D, Seitz AE, Lawyer P, Garfield M, Pham M, Valenzuela JG: Comparative salivary gland transcriptomics of sandfly vectors of visceral leishmaniasis. BMC Genomics. 2006, 7: 52-10.1186/1471-2164-7-52.PubMed CentralPubMedGoogle Scholar
- Ribeiro JM, Arca B: From sialomes to the sialoverse: an insight into the salivary potion of blood feeding insects. Adv Insect Physiol. 2009, 37: 59-118.Google Scholar
- Francischetti IM, Sa-Nunes A, Mans BJ, Santos IM, Ribeiro JM: The role of saliva in tick feeding. Front Biosci. 2009, 14: 2051-2088.Google Scholar
- Aalberse RC: Structural biology of allergens. J Allergy Clin Immunol. 2000, 106 (2): 228-238. 10.1067/mai.2000.108434.PubMedGoogle Scholar
- Jacquenet S, Moneret-Vautrin DA: Allergy to peanut and all kinds of nuts. Revue française d’allergologie et d’immunologie clinique. 2007, 47: 487-491.Google Scholar
- Peng Z, Xu W, James AA, Lam H, Sun D, Cheng L, Simons FE: Expression, purification, characterization and clinical relevance of rAed a 1–a 68-kDa recombinant mosquito Aedes aegypti salivary allergen. Int Immunol. 2001, 13 (12): 1445-1452. 10.1093/intimm/13.12.1445.PubMedGoogle Scholar
- Cooper RD, Edstein MD, Frances SP, Beebe NW: Malaria vectors of Timor-Leste. Malar J. 2010, 9: 40-10.1186/1475-2875-9-40.PubMed CentralPubMedGoogle Scholar
- Garros C, Van Nguyen C, Trung HD, Van Bortel W, Coosemans M, Manguin S: Distribution of Anopheles in Vietnam, with particular attention to malaria vectors of the Anopheles minimus complex. Malar J. 2008, 7: 11-10.1186/1475-2875-7-11.PubMed CentralPubMedGoogle Scholar
- Rizzo C, Ronca R, Fiorentino G, Verra F, Mangano V, Poinsignon A, Sirima SB, Nebie I, Lombardo F, Remoue F, et al: Humoral response to the Anopheles gambiae salivary protein gSG6: a serological indicator of exposure to Afrotropical malaria vectors. PLoS One. 2011, 6 (3): e17980-10.1371/journal.pone.0017980.PubMed CentralPubMedGoogle Scholar
- Poinsignon A, Cornelie S, Ba F, Boulanger D, Sow C, Rossignol M, Sokhna C, Cisse B, Simondon F, Remoue F: Human IgG response to a salivary peptide, gSG6-P1, as a new immuno-epidemiological tool for evaluating low-level exposure to Anopheles bites. Malar J. 2009, 8: 198-10.1186/1475-2875-8-198.PubMed CentralPubMedGoogle Scholar
- Poinsignon A, Cornelie S, Mestres-Simon M, Lanfrancotti A, Rossignol M, Boulanger D, Cisse B, Sokhna C, Arca B, Simondon F, et al: Novel peptide marker corresponding to salivary protein gSG6 potentially identifies exposure to Anopheles bites. PLoS One. 2008, 3 (6): e2472-10.1371/journal.pone.0002472.PubMed CentralPubMedGoogle Scholar
- Poinsignon A, Samb B, Doucoure S, Drame PM, Sarr JB, Sow C, Cornelie S, Maiga S, Thiam C, Rogerie F, et al: First attempt to validate the gSG6-P1 salivary peptide as an immuno-epidemiological tool for evaluating human exposure to Anopheles funestus bites. Trop Med Int Health. 2010, 15 (10): 1198-1203. 10.1111/j.1365-3156.2010.02611.x.PubMedGoogle Scholar
- Kalume DE, Okulate M, Zhong J, Reddy R, Suresh S, Deshpande N, Kumar N, Pandey A: A proteomic analysis of salivary glands of female Anopheles gambiae mosquito. Proteomics. 2005, 5 (14): 3765-3777. 10.1002/pmic.200401210.PubMedGoogle Scholar
- Wu L, Han DK: Overcoming the dynamic range problem in mass spectrometry-based shotgun proteomics. Expert Rev Proteomics. 2006, 3 (6): 611-619. 10.1586/147894126.96.36.1991.PubMedGoogle Scholar
- 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 (4): 520-529.PubMedGoogle Scholar
- Bonizzoni M, Afrane Y, Yan G: Loop-mediated isothermal amplification (LAMP) for rapid identification of Anopheles gambiae and Anopheles arabiensis mosquitoes. Am J Trop Med Hyg. 2009, 81 (6): 1030-1034. 10.4269/ajtmh.2009.09-0333.PubMed CentralPubMedGoogle Scholar
- Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16 (2): 74-77. 10.1016/S0169-4758(99)01563-X.PubMedGoogle Scholar
- Fontaine A, Pascual A, Orlandi-Pradines E, Diouf I, Remoue F, Pages F, Fusai T, Rogier C, Almeras L: Relationship between exposure to vector bites and antibody responses to mosquito salivary gland extracts. PLoS One. 2011, 6 (12): e29107-10.1371/journal.pone.0029107.PubMed CentralPubMedGoogle Scholar
- Fontenille D, Cohuet A, Awono-Ambene P, Kengne P, Antonio-Nkondjio C, Wondji C, Simard F: Malaria vectors: from the field to genetics. Research in Africa. Rev Epidemiol Sante Publique. 2005, 53 (3): 283-290. 10.1016/S0398-7620(05)84605-X.PubMedGoogle Scholar
- Ambrosino E, Dumoulin C, Orlandi-Pradines E, Remoue F, Toure-Balde A, Tall A, Sarr JB, Poinsignon A, Sokhna C, Puget K, et al: A multiplex assay for the simultaneous detection of antibodies against 15 Plasmodium falciparum and Anopheles gambiae saliva antigens. Malar J. 2010, 9: 317-10.1186/1475-2875-9-317.PubMed CentralPubMedGoogle Scholar
- Souza AP, Andrade BB, Aquino D, Entringer P, Miranda JC, Alcantara R, Ruiz D, Soto M, Teixeira CR, Valenzuela JG, et al: Using recombinant proteins from Lutzomyia longipalpis saliva to estimate human vector exposure in visceral Leishmaniasis endemic areas. PLoS Negl Trop Dis. 2010, 4 (3): e649-10.1371/journal.pntd.0000649.PubMed CentralPubMedGoogle Scholar
- Fontaine A, Pophillat M, Bourdon S, Villard C, Belghazi M, Fourquet P, Durand C, Lefranc D, Rogier C, Fusai T, et al: Specific antibody responses against membrane proteins of erythrocytes infected by Plasmodium falciparum of individuals briefly exposed to malaria. Malar J. 2010, 9 (1): 276-10.1186/1475-2875-9-276.PubMed CentralPubMedGoogle Scholar
- 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 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.Google Scholar
- Transcriptome Resources.http://exon.niaid.nih.gov/transcriptome.html,
- Machault V, Gadiaga L, Vignolles C, Jarjaval F, Bouzid S, Sokhna C, Lacaux JP, Trape JF, Rogier C, Pages F: Highly focused anopheline breeding sites and malaria transmission in Dakar. Malar J. 2009, 8: 138-10.1186/1475-2875-8-138.PubMed CentralPubMedGoogle Scholar
- Fontenille D, Lochouarn L, Diatta M, Sokhna C, Dia I, Diagne N, Lemasson JJ, Ba K, Tall A, Rogier C, et al: Four years' entomological study of the transmission of seasonal malaria in Senegal and the bionomics of Anopheles gambiae and A. arabiensis. Trans R Soc Trop Med Hyg. 1997, 91 (6): 647-652. 10.1016/S0035-9203(97)90506-X.PubMedGoogle Scholar
- Fontaine A, Pascual A, Diouf I, Bakkali N, Bourdon S, Fusai T, Rogier C, Almeras L: Mosquito salivary gland protein preservation in the field for immunological and biochemical analysis. Parasit Vectors. 2011, 4: 33-10.1186/1756-3305-4-33.PubMed CentralPubMedGoogle Scholar
- Moneret-Vautrin DA: Allergic risk and role of the allergy vigilance network. Bull Acad Natl Med. 2007, 191 (4–5): 807-814. discussion 814PubMedGoogle Scholar
- Davey RJ, Lenes BL, Casper AJ, Demets DL: Adequate survival of red cells from units "undercollected" in citrate-phosphate-dextrose-adenine-one. Transfusion. 1984, 24 (4): 319-322. 10.1046/j.1537-2995.1984.24484275572.x.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.