Comparative salivary gland transcriptomics of sandfly vectors of visceral leishmaniasis
- Jennifer M Anderson†1,
- Fabiano Oliveira†1, 2,
- Shaden Kamhawi3,
- Ben J Mans4,
- David Reynoso1,
- Amy E Seitz1,
- Phillip Lawyer3,
- Mark Garfield5,
- MyVan Pham1 and
- Jesus G Valenzuela1Email author
© Anderson et al; licensee BioMed Central Ltd. 2006
Received: 06 January 2006
Accepted: 15 March 2006
Published: 15 March 2006
Immune responses to sandfly saliva have been shown to protect animals against Leishmania infection. Yet very little is known about the molecular characteristics of salivary proteins from different sandflies, particularly from vectors transmitting visceral leishmaniasis, the fatal form of the disease. Further knowledge of the repertoire of these salivary proteins will give us insights into the molecular evolution of these proteins and will help us select relevant antigens for the development of a vector based anti-Leishmania vaccine.
Two salivary gland cDNA libraries from female sandflies Phlebotomus argentipes and P. perniciosus were constructed, sequenced and proteomic analysis of the salivary proteins was performed. The majority of the sequenced transcripts from the two cDNA libraries coded for secreted proteins. In this analysis we identified transcripts coding for protein families not previously described in sandflies. A comparative sandfly salivary transcriptome analysis was performed by using these two cDNA libraries and two other sandfly salivary gland cDNA libraries from P. ariasi and Lutzomyia longipalpis, also vectors of visceral leishmaniasis. Full-length secreted proteins from each sandfly library were compared using a stand-alone version of BLAST, creating formatted protein databases of each sandfly library. Related groups of proteins from each sandfly species were combined into defined families of proteins. With this comparison, we identified families of salivary proteins common among all of the sandflies studied, proteins to be genus specific and proteins that appear to be species specific. The common proteins included apyrase, yellow-related protein, antigen-5, PpSP15 and PpSP32-related protein, a 33-kDa protein, D7-related protein, a 39- and a 16.1- kDa protein and an endonuclease-like protein. Some of these families contained multiple members, including PPSP15-like, yellow proteins and D7-related proteins suggesting gene expansion in these proteins.
This comprehensive analysis allows us the identification of genus- specific proteins, species-specific proteins and, more importantly, proteins common among these different sandflies. These results give us insights into the repertoire of salivary proteins that are potential candidates for a vector-based vaccine.
Phlebotomine sandflies are vectors of Leishmania parasites, causal agents of leishmaniasis in at least 88 countries. The manifestations of this disease range from the self-healing cutaneous and mucocutaneous forms to the potentially fatal visceral form. The incidence of leishmaniasis is 2 million cases annually, of which 500,000 cases are visceral and potentially lethal . Visceral leishmaniasis is caused by parasites of the Leishmania donovani complex: L. donovani, L. infantum and L. chagasi (L. infantum chagasi). There are a limited number of competent sandfly vectors that can transmit parasites within this complex. For example, Phlebotomus argentipes transmits only L. donovani in the India sub-continent, P. ariasi and P. perniciosus transmit L. infantum within southern Europe, and Lutzomyia longipalpis exclusively transmits L. chagasi (L. infantum chagasi) in Central and South America.
Infected sandflies deliver the Leishmania parasite to a mammalian host during blood meal acquisition. Together with the parasite, sandflies inject saliva into the host. This saliva contains potent pharmacologically active components that facilitate blood feeding . Additionally, the saliva affects the establishment of the parasite within the vertebrate host; small amount of L. longipalpis saliva exacerbates L. major infection in mice . On the other hand, immune responses to sandfly saliva have been shown to protect against Leishmania infection [4, 5]. Antibodies to maxadilan, a salivary protein from the sandfly L. longipalpis protected mice against L. major infection , while cellular immune response to PpSP15, a protein from the sandfly P. papatasi was sufficient to control L. major infection in mice . Therefore, immune responses to salivary protein have promise as an effective vector-based vaccine to control Leishmania infection.
The repertoire of sandfly salivary proteins that have been studied is largely limited to three sandflies: P. papatasi, P. ariasi and L. longipalpis, vectors of L major, L. infantum and L. chagasi (L. infantum chagasi), respectively. Only two salivary proteins have been extensively studied: maxadilan from the sandfly L. longipalpis, and PpSP15 from the sandfly P. papatasi. Maxadilan, a potent vasodilatory, immunomodulatory and protective molecule was shown to be very polymorphic . On the other hand, PpSP15, a protective molecule with unknown biological function was shown to be highly conserved among colonized or field-collected P. papatasi sandflies .
Because of the potential of sandfly salivary proteins as anti-Leishmania vaccines, it is important to understand the diversity and degree of similarity between salivary proteins from various sandflies. More importantly, understanding the evolutionary relatedness of salivary proteins will help us to identify proteins that can be used as a global or general anti-Leishmania vaccine within a complex of vector species.
Here we explored the proteins and transcripts encoded in the salivary glands of the sandflies P. argentipes and P. perniciosus and studied the repertoire of proteins on these sandflies and compared them with the P. ariasi and L. longipalpis salivary proteins, also vectors of visceral leishmaniasis. We studied their molecular characteristics in the relation to molecular evolution of sandfly salivary proteins, and in the context of possible scenarios for global vector-based vaccines strategies.
Results and discussion
Sequencing of P. argentipes and P. perniciosus salivary gland cDNA libraries
From the P. argentipes salivary gland cDNA library, we sequenced 603 randomly selected clones from which 135 unique clusters of related sequences were obtained. Out of the 135 clusters, we found 45 clusters (1.11 sequences per cluster) of transcripts coding for housekeeping genes. We found 111 sequences, arranged in 55 clusters (1.8 sequences per cluster) that were not similar to other genes in the NCBI databank and lacked a secretory signal peptide. The most abundant transcripts in this cDNA library contained putative secretory proteins. We found 438 cDNA with potentially secreted proteins signals arranged in 30 clusters (an average of 14.36 sequences per cluster). The number of cDNA coding for secretory proteins is 9 times greater than the cDNA coding for housekeeping genes and 4 times greater than the cDNA coding for non-secreted proteins with unknown function. The transcripts coding for secretory proteins represent 73% of the total transcripts sequenced in the P. argentipes salivary gland library.
Putative secreted proteins from the salivary glands of Phlebotomus argentipes.
NCBI accession number
Signal P site
Best match to NR protein database
Present in proteome
SL1 protein L. longipalpis
Similar to PpSP15
SL1 protein L. longipalpis
Similar to PpsP15
44 kDa salivary protein
Antigen 5 L. longipalpis
Antigen 5 protein
32 kDa salivary protein
Similar to PpSP32
SL1 protein L. longipalpis
Similar to PpSP15
32 kDa protein L. longipalpis
28 kDa salivary protein
D7 related protein
Endonuclease L. longipalpis
SL1 protein L. longipalpis
Similar to PpSP15
14 kDa salivary protein
Similar to PpSP15
ebiP3881 An gambiae
30 kDa salivary protein
putative histone promoter
A. gambiae unknown
agCP4255 An. gambiae
Putative secreted proteins from the salivary glands of Phlebotomus perniciosus.
NCBI accession number
Signal P site
Best match to NR protein database
Present in proteome
Salivary apyrase P. papatasi
Salivary apyrase P. papatasi
SL1 protein L. longipalpis
SP15 like protein
42 kDa salivary prot. P. papatasi
44 kDa salivary prot. P. papatasi
28 kDa salivary prot. P. papatasi
28 kDa salivary prot. P. papatasi
29 kDa salivary prot. L. longipalpis
32 kDa salivary prot. L. longipalpis
Antigen 5 prot. L. longipalpis
Antigen 5 protein
Salivary prot. C. sonorensis
14 kDa salivary prot. P. papatasi
30 kDa salivary prot. P. papatasi
SL1 prot. L. longipalpis
SP15 like protein
Phospholipase A2, Drosophila
37 kDa prot.. L. longipalpis
L. longipalpis endonuclease
In addition to the identification of proteins previously reported from other sandflies, we found a number of transcripts coding for proteins not previously shown to be present in the salivary glands of sandflies. A protein homologous to lipases from Anopheles gambiae, Drosophila melanogaster and other organisms was found in the P. argentipes cDNA library. We also found in this library and in the P. perniciosus cDNA library, a transcript coding for a protein homologous to a pyrophosphatase. The predicted 47-kDa protein named PagSP132 contains a phosphodiesterase type I, phosphodiesterase/nucleotide pyrophosphatase motif. This type of enzymes cleaves the phosphodiester and phosphosulfate bonds in NAD, deoxynucleotides and nucleotide sugars . BLAST search of this protein identified protein orthologs found in mammals as well as in A. gambiae. To our knowledge, pyrophosphatases have not been described in the saliva of other sandflies.
We found one cluster in the P. perniciosus cDNA library coding for a phospholipase A2 (PLA2) protein (PpeSP18). This type of protein has never been reported from the saliva of a blood-feeding insect. PLA2 (Phosphatidylcholine-2-acylhydrolase, E.C. 188.8.131.52) are well known for their ability to cleave the arachidonic acid and lysophosphatidylcholine from the sn-2 position of membrane glycerol-3-phospholipids. Also PLA2 are known to work as toxins by blocking the release of neurotransmitters .
We identified transcripts coding for secreted proteins that did not match any reported proteins in accessible databases. P. argentipes contained six unknown proteins that ranged from 15 to 30 kDa, while only three were found in the P. perniciosus library and all were relatively small ranging from 10 to 27 kDa (Tables 1 and 2).
Proteome analysis of P. argentipes and P. perniciosus salivary proteins
From P. perniciosus salivary gland protein analysis we found 13 N-terminus sequences (Figure 1B). The identified proteins included: three PpSP15-like proteins (PpeSP11, PpeSP02 and PpeSP09), three D7-related proteins (PpeSP04, PpeSP04B, PpeSP10), a 37-kDa-like protein described previously in the saliva of L. longipalpis (PpeSP19), an antigen-5 related protein (PpeSP07), two apyrase-like proteins (PpeSP01, PpeSP01B), a 32-kDa-like salivary protein described on L. longipalpis (PpeSP06) and two yellow-related proteins (PpeSP03, PpeSP03B). Not all attempted Edman degradation experiments resulted in a sequence, either because of insufficient protein amount or because N-terminal ends were blocked.
Comparative analysis of salivary transcripts from vectors of visceral leishmaniasis: P. argentipes, P. perniciosus, P. ariasi and L. longipalpis
In an attempt to understand the relationship of salivary proteins among different sandflies and to gain insights into the evolution of sandfly salivary proteins, we compared cDNA libraries from four different sandflies. We selected the sandflies based on their availability and their significance for this study. First, we selected sandflies from the two different genera, Phlebotomus and Lutzomyia. Secondly, from the Phlebotomus genus we selected two different subgenus, Euphlebotomus (P. argentipes) and Larroussius (P. ariasi and P. perniciosus). The phylogenetic relationship among these sandflies was previously studied using the small subunit nuclear ribosomal DNA .
Salivary transcripts shared by Phlebotomus and Lutzomyia sandflies.
Family of proteins
PagSP01, 02, 07, 12, 13
PpeSP02, 09, 11
LJM17, LJM11, LJM111
Antigen 5-related protein
33 kDa, unknown function
ParSP07, 12, 16
PpeSP04, 04B, 10
39 kDa, unknown function
16.1 kDa, unknown function.
It is interesting to note the amount of variation that exists in the number of members of the different protein families found in the four sandflies (Table 3). The PpSP15-like family has five members identified in P. argentipes (PagSP01, PagSP02, PagSP07, PagSP12 and PagSP13), two in P. ariasi (ParSP03, ParSP08), three in P. perniciosus (PpeSP02, PpeSP09 and PpeSP11), yet only one (LloSP05) in L. longipalpis (Table 3). On the other hand, only one member of the apyrase family of proteins has been found in each sandfly except for P. perniciosus, which has two members (Table 3). Other families of salivary proteins, such as antigen 5, PpSP32, 32 kDa and endonuclease-like protein were represented by only one member from each of the different sandflies. The yellow-related protein has one member found in P. argentipes, two members in the P. perniciosus and P. ariasi and three members in the Lutzomyia longipalpis sandfly. The D7-related protein was represented by three members in Phlebotomus sandflies, while only one member is present in the L. longipalpis sandfly.
Salivary transcripts shared by Phlebotomus sandflies.
32 kDa, unknown function
2 kDa, unknown function
5 kDa, unknown function
Molecular characteristics of salivary proteins shared among the analysed sandflies
In order to understand the relationship among salivary proteins from different sandflies, we performed multiple sequence alignment followed by phylogenetic analysis of the salivary transcripts shared by the vectors of visceral leishmaniasis studied (Table 3). Following is a description of the shared proteins:
SL1/PpSP15 related proteins
The SL1/PpSP15 group of proteins is similar to the SL1 salivary protein from L. longipalpis, which has no known function , and to PpSP15, a 15- kDa salivary protein from Phlebotomus papatasi that was previously shown to confer protection against L. major infection . The predicted molecular weight of these transcripts is approximately 14 kDa and is in agreement with the observed MW found through the proteome analysis (Figure 1). This group represents the most abundant transcripts in the salivary gland cDNA library of P. argentipes (Table 1). The fact that only one PpSP15 member was found in L. longipalpis, suggests that a number of lineage-specific gene expansions (gene duplication events) occurred in the Phlebotomus lineage at various periods in the evolution of these sandflies.
The PpSP15 family of proteins has only been found in species of sandflies suggesting that this family was a specific invention that occurred during sandfly evolution, most likely during their adaptation to a blood-feeding environment. Although PSI-BLAST analysis using PpSP15 proteins retrieved only members of the PpSP15 family, the PHYRE prediction servers indicated that members of the PpSP15 family possess an EF-hand fold most closely related to members of the odorant-binding protein (OBP) family to which the D7-proteins belong. It is thus likely that PpSP15 members were derived from an OBP ancestral protein. Characteristically, the OBP family in Drosophila has a low degree of sequence similarity among its members with only six conserved cysteines among the thirty-four members of this family .
Phylogenetic analysis of PpSP15 transcripts from the different sandflies including P. papatasi resulted in the formation of 5 distinct clades (Figure 2B). As such, three major clusters of orthologous groups of proteins (COGs), and hence gene duplication events, can be identified that possibly occurred in the ancestor to the Phlebotomus lineage. COG1 includes members from P. perniciosus, P. argentipes and P. papatasi with a second related clade that includes a lineage-specific expansion (LSE) in P. papatasi. COG2 includes members from P. perniciosus and P. argentipes. COG3 includes members from P. perniciosus, P. argentipes and P. ariasi. Another clade composed solely of members of P. argentipes indicates another LSE. Evolutionarily, the P. papatasi group was basal to the other Phlebotomus members analysed in this study, followed by P. argentipes, with P. perniciosus and P. ariasi forming the terminal clade. Given this, the topology of the cladogram obtained for the PpSP15 family suggests that COG2 resulted from gene duplication event that occurred in the ancestor to P. perniciosus and P. ariasi. COG1 follows the expected phylogenetic grouping but suggests that this gene was lost in P. ariasi or we failed to detect the ortholog in our library. COG3 again suggests that this specific gene duplication event occurred after divergence from the shared ancestor with P. papatasi. The PpSP15 family found in sandflies is thus characterised by both gene duplication and possibly gene loss events, both restricting an accurate reconstruction of its phylogeny. It is currently impossible to say which proteins share a conserved function with the PpSP15 from L. longipalpis as the major clades create a polytomy. As such, this family might bind related or similar pharmacologic components so that they all have, in fact, a similar function. This might explain the seemingly haphazard acquisition and loss of genes.
D7 family of proteins
D7-related proteins are found in the saliva of different diptera, including Anopheles , Aedes  and Culex  mosquitoes as well as in the sandflies P. papatasi , L. longipalpis  and P. ariasi  and does not appear to occur in non-dipteran species. Two forms of the protein have been described; a long and short form ; and appear to be distantly related to an odorant binding super family of proteins. Interestingly, the OBP family seems to be the ancestral molecule of the PpSP15 family (see above). Therefore, it may be possible that both D7 and PpSP15 related proteins have a common ancestor.
The D7 protein named hamadarin from Anopheles stephensi acts as an anticoagulant affecting the plasma contact system by inhibiting the activation of Factor XII and kallikrein . Recently, a biological function of four short members of the D7 family from A. gambiae and a long D7 from Aedes aegypti was described . These salivary proteins were shown to bind biogenic amines such as serotonin, histamine and norepinephrine. This function is relevant for blood-feeding because of the inhibition of the vasoconstrictor, platelet aggregating, and pain inducing properties of these biogenic amines . The exact function of D7 proteins in sandflies is largely unknown, but it may be related to the function observed in mosquito D7 proteins, either as an anticoagulant or binding biogenic amines.
The D7 family is represented in the P. argentipes cDNA library by two members, PagSP10 and PagSP25, and in the P. perniciosus cDNA library by three members, PpeSP04, PpeSP04B and PpeSP10. Only one member of this family is present in L. longipalpis sandfly, suggesting a case of gene duplication of this protein that probably occurred more recently in the Phlebotomus genus.
Phylogenetic analysis of D7 proteins from different organisms shows 2 distinct clades with the sandfly proteins branching from the long form (Figure 3B). All of the sandfly D7 members are clustered within one clade, distinct from the long form members. The clade containing the mosquito D7 short-form proteins contains two clusters, one containing the Anopheles mosquitoes and the other containing the Culex and Aedes species. The sandfly clade subdivides into two distinct clades. The lower clade shows 5 distinct groups, two of these groups represent COGs of P. perniciosus and P. ariasi, and the group of P. papatasi seems to be a case of lineage expansion. The lower clade shows a COG that includes P. argentipes, P. ariasi and P. perniciosus proteins, PagSP25, ParSP07 and PpeSP10.
Apyrase family of proteins
Both the P. perniciosus and P. argentipes libraries contained transcripts homologous to the Cimex family of apyrases , a protein also present in the saliva of P. papatasi  and L. longipalpis . Apyrases are enzymes that function as potent anti-platelet factors by destroying or hydrolysing the platelet activator ADP. An orthologue was found in humans and the recombinant protein was shown to hydrolyse a variety of nucleoside di- and triphosphates, preferentially UDP, followed by GDP, UTP, GTP, ADP, and ATP [25, 26].
Phylogenetic analysis of apyrases from different organisms indicates three main clades with the sandfly apyrases in a distinct clade, apart from vertebrates, yet closely related to other insects (Figure 4B). Interestingly sandfly apyrases share a common ancestor with Cimex lectularius apyrase. The two insects appear to have evolved to the blood feeding mode independently . Within the clade containing the sandflies, one of the P. perniciosus apyrases, PpeSP01, is more closely related to the P. ariasi apyrase (ParSP01) than the second apyrase from P. perniciosus (PpSP01B). This may be the result of a gene duplication event in P. perniciosus and subsequent loss in P. ariasi. When searching databases we found a transcript from A. gambiae similar to sandfly apyrases. This is interesting because the known mosquito apyrases belong to the 5'-nucleotidase family of proteins. The known mosquito apyrase is very distinct, in size and sequence, from the Cimex family of apyrases also present in sandflies [28, 24, 7], thus it is possible that mosquitoes in addition to a functional 5'-nucleotidase type apyrase, may have a non-functional sandfly/bedbug-like apyrase gene or it may have a house keeping function such as hydrolysing UDP formed after transglycosylation reactions in the Golgi . Alternatively, the mosquito may have a similar apyrase but with different substrate specificity or this protein is not present in their salivary gland.
Interestingly, the mosquito apyrase gene seems to be ancestral to the sandfly apyrase based on phylogenetic association (Figure 4B). Then it may be possible that mosquitoes have lost the function of this gene and kept the active form of the 5'-nucleotidase gene, which is the active apyrase in these insects.
The gene coding for the yellow protein was first described in Drosophila melanogaster . The proteins in this family appear to be derived from a common ancestor of the major royal jelly proteins (MRJPs) from honeybees and the yellow protein from Drosophila spp.
Drosophila yellow protein is related to pigmentation and male sexual behavior. In the family Culicidae a yellow protein was identified in Ae. aegypti whole-larvae extract and was associated to a dopachrome converting enzyme activity found in this insect . The function of this protein in the saliva of sandflies and its importance for blood feeding remains to be elucidated. The yellow protein family is one of the most abundant proteins found in the sandfly saliva.
We identified transcripts coding for secreted proteins of approximately 45 kDa, previously described in the saliva of L. longipalpis, P. papatasi and P. ariasi as yellow-related proteins [13, 14, 21]. In the P. argentipes cDNA library we found only one cluster (PagSP04) coding for this protein, yet in P. papatasi and L. longipalpis salivary glands there are multiple members of this family of proteins . Proteomic analysis (Figure 1A) revealed that the PagSP04 transcript (YHVEREYAWRNVTFEGVN) was one of the most abundant proteins found in the salivary glands of P. argentipes. Interestingly, three proteins with different mobilities coded for the same N-terminus sequence (Figure 1A) suggesting they may represent the same protein with different post-translational modifications.
Based on their MW, the yellow proteins from L. longipalpis appear to be the most recognised proteins from the sera of individuals living in endemic areas of visceral leishmaniasis and from individuals that have anti-Leishmania immunity . The antibody response against these salivary proteins appears to correlate with protection against leishmaniasis.
Antigen-5 family of proteins
This cluster codes for a secreted protein of 29 kDa similar to antigen 5-related protein found in wasp venom . Similar proteins have been isolated from the salivary glands of Aedes aegypti , An. gambiae  and from the salivary glands of L. longipalpis . We found only one cluster coding for this protein in the cDNA library of P. argentipes (PagSP05) and in the cDNA library of P. perniciosus (PpeSP07). The N-terminal sequence corresponding to these transcripts was identified in the salivary glands of P. argentipes (Figure 1A) and P. perniciosus (Figure 1B).
This family of proteins belong to the CAP family (CRISP, Ag5,PR-1) of proteins [37, 35]. A remarkable feature of this family is the large number of cysteine residues, particularly at the carboxy-terminal region. X-ray structure of Na-AS-2, a member of this family from the human hookworm Necator americanus, was recently reported  and showed structural similarities to chemokines. Thus, it is possible that this type of protein in sandflies or other insects may bind cytokines with potential effects on the host immune response.
33-kDa protein family
Endonuclease family of proteins
Endonuclease proteins were found in all sandflies studied. The multiple alignment of sandfly endonuclease showed various regions of identity among amino acids, and many regions of conserved amino acids, even when comparing endonucleases from different sandfly genera (Figure 10A). Phylogenetic analysis indicated that the sandfly endonucleases clustered (92% bootstrap support) with other arthropods, including two non-insect arthropods (Paralithode camtschaticus and Marsupenaeus japonicus) (Figure 10B), however, sandfly endonucleases formed a distinct clade within the arthropod cluster (90% bootstrap support). Additionally, these endonucleases were clearly distant from other endonucleases. This suggests that this endonuclease may represent a common antigen for different sandflies and that the distant relationship to other endonucleases may avoid potential cross reactivity with non-insect organisms. Since this type of enzyme can cleave double- and single-stranded DNA, the role of this protein in the saliva of sandflies should be investigated further.
The role of a collagen-like protein in sandfly salivary glands is unknown, yet intriguing. Because this protein may bind matrix protein, it is possible that this type of protein may form strong associations with basal matrix proteins. Within the Phlebotomus genera the relationship to collagen appears to be limited to PpSP32 kDa-like proteins from P. pernicious. P. papatasi and P. ariasi share only weak (non-significant) homology with Type VII collagen (E = 2.8 and 4.5, respectively) and P. argentipes did not match any collagen proteins upon multiple BLAST searches.
One explanation for the apparent homology between the four sandflies studied, yet lack of identity between collagen-related proteins in all sandflies, was revealed in the comparative analysis of the four sandflies. The majority of the homology between the four sandflies was found at the N- and carboxy terminus whereas the middle section of the protein appeared to be less conserved. The N-terminal and carboxy terminus has 43% and 23% identity between the four flies, respectively, whereas the midsection of the protein contains only 4% identity (Figure 11A). The region of homology between P. perniciosus and Type VII collagen was found in the non-homologous region of the protein (Figure 11B).
Additionally, BLAST analysis of L. longipalpis LJL04 (AAS16906) revealed a significant homology to collagen adhesion proteins from B. thuringiensis (ZP_00739782) and B. cereus (NP_830673) (Figure 11C). Again, the region of homology was found in the divergent section of the PpSP32-kDa protein. The most perplexing aspect of this family of proteins is that, although highly conserved among the four sandflies studied here, the N-terminal and carboxy regions of the protein have no homology to proteins of known function.
Perspectives on the evolution of hematophagy in sandflies
The presence of a specific gene in distantly related species could indicate a true orthologous relationship (i.e., the presence of the gene in the ancestral species) or acquisition of the gene by one of the species by horizontal gene transfer. While horizontal gene transfer is a common event in prokaryotes its provenance in eukaryotes is not well established and the general consensus is that its occurrence is rare, if present at all. In the case of sandflies this would imply that proteins present in evolutionary distant species (insects or other metazoans) were present in the ancestral sandfly. As such, two general trends can be observed for the sandfly protein families found in their transcriptomes. Firstly, those proteins that occur throughout the metazoans or insecta tend to be found as single members in all sandfly genera with no extensive gene duplication events occurring in sandflies. These proteins are generally well conserved and possibly share the same or a very similar function to those found in the main family. These proteins include the apyrase, antigen 5 and endonuclease families and probably consist of the core repertoire of the ancestral sandfly proteins that develop during adaptation to hematophagous behavior. Although yellow-related proteins and D7 proteins are generally well conserved, we observed some gene expansion of these proteins in sandflies (Figures 3 and 7).
Alternatively, a number of protein families are limited to sandflies as a group or specific species and show low levels of similarity to family members found in other insects. Members of this group show high levels of divergence, more gene duplication events and were probably evolved specifically during adaptation to a blood-feeding lifestyle and specific host species. Proteins in this category include the PSP15 like, PSP32 like, 32 kDa, 39 kDa and 16.1 kDa protein families (Table 3) Within this group can also be placed the singletons, which are proteins limited exclusively to single species (Table 4). In the case of the 39 kDa and 16.1 kDa families (Table 3), it is possible that gene losses occurred among selected sandfly species or the sequences or transcripts were missed on this analysis due to low representation of these transcripts in the sandflies studied. In the case of the 2- and 5-kDa salivary proteins only found in P. ariasi and P. perniciosus (Table 4), these proteins may be specific for the subgenus Larroussius. We found many transcripts in the P. argentipes and L. longipalpis cDNA libraries coding for proteins of similar MW, however with no significant homology to these two proteins.
As expected, we found that salivary proteins from sand flies belonging to the same subgenus are more closely related to each other than proteins from different subgenus. We observed that proteins belonging to P. ariasi (subgenus Larroussius) in all phylogenetic tree analysis based on protein sequence were more closely related to P. perniciosus (subgenus Larroussius) than to P. argentipes salivary proteins (subgenus Euphlebotomus). These results are in agreement with previously studies using the small subunit nuclear ribosomal DNA .
Can we use the salivary proteins common to these four sandflies as global antigens for a vector-based vaccine?
The large degree of divergence found in the majority of the most abundant sandfly salivary proteins suggests that a specific salivary protein may not be used as common vaccine target or as a common marker for sandfly exposure to different sandfly genera. This is supported by the recent findings by Rohousova et al , who compared antibody responses to salivary proteins from P. papatasi, P. sergenti and L. longipalpis and also demonstrated a lack of cross-reactivity between animals bitten by two different sandfly genera. This lack of cross-reactivity is possibly due to the low degree of similarities found in the Lutzomyia and Phlebotomus salivary proteins, as reported in the present work, and to the lack of recognition to specific molecules present exclusively in the different genera. We cannot exclude at this point potential cross-reactivity, or lack thereof, with components of the cellular immune response to sandfly salivary proteins. This area has been under-studied and should be evaluated experimentally. Although there is an overall low degree of identity between most salivary proteins across the genera, there are small regions of identities that may contain common T cell epitopes conserved between Lutzomyia and Phlebotomus salivary proteins. These small regions are observed in proteins such as D7, apyrases, yellow-related protein, antigen 5, a 33-kDa protein and endonucleases from these sandflies studied thus far.
In contrast, we found salivary proteins that have the potential to be a common vaccine target within the genus Phlebotomus. We identified three proteins that are highly conserved in different Phlebotomus species, the yellow-related protein, the apyrases and the antigen 5-related proteins. Additionally, these proteins have the potential to be markers of exposure for Phlebotomus sandflies in general. This is supported by observations by Rohousova et al  and Volf and Rohousova , which showed some cross-reactivity in animals to salivary proteins between different Phlebotomus species. These observations need to be expanded and evaluated experimentally for the potential cross-reactivity of these salivary proteins on specific cellular immune responses that may protect against Leishmania infection.
Overall, this study led us to the identification of novel salivary proteins from two sandfly vectors of visceral leishmaniasis, P. argentipes and P. perniciosus, and the identification of the overall repertoire of secreted proteins present in their salivary glands. Additionally, this study allowed the discovery of the salivary proteins common among four different sandflies, from two different genera (Lutzomyia and Phlebotomus) and from two different subgenera (Euphlebotomus and Larroussius). This comparative study is giving us insight into the evolution of sandfly salivary proteins, their relationship and their molecular characteristics.
Moreover, this study is providing a better understanding of the overall sequence identity of sandfly salivary proteins across genus and species, while suggesting that a global vector-based vaccine may not be possible across different genera and that possibly genus- or species-specific salivary proteins may need to be used for this type of vaccine.
Adult Phlebotomus argentipes (NIH colony) and P. perniciosus (kindly obtained by Dr. Michele Maroli, Italy) sandflies were kept with free access to a 20% solution of sucrose. Salivary glands from recently emerged and 1- to 2-day-old adult female flies were dissected and transferred to 10 or 20 ul 10 mM HEPES pH 7.0, 0.15 M NaCl in 1.5 ml polypropylene vials, usually in groups of 10 pairs of glands in 20 ul of HEPES saline. Salivary glands were kept at -75°C until needed.
Salivary Gland cDNA Library
Salivary gland mRNA from both species was isolated from 40–50 salivary gland pairs, respectively, using the Micro-FastTrack mRNA isolation kit (Invitrogen, San Diego, CA). The PCR-based cDNA library was made using the SMART cDNA library construction kit (BD-Clontech, Palo Alto, CA), following the manufacturer's recommendation with some modifications . The obtained cDNA libraries (large, medium and small sizes) were plated by infecting log phase XL1- blue cells (Clontech); insert size was determined with PCR using vector primers flanking the inserted cDNA and visualised on a 1.1 % agarose gel with ethidium bromide (1.5 ug/ml). Inserts were sequenced as previously described using a CEQ 2000XL DNA sequencing instrument (Beckman Coulter) .
Detailed description of the bioinformatic treatment of the data appear in . Briefly, primer and vector sequences were removed from raw sequences and sequences shorter than 50 nucleotides or containing more than 15% N were removed from further analysis. Sequences were compared to the GenBank non-redundant (nr) protein database using the standalone BlastX program  using a cut-off E-value of 1 × 10-5. Related sequences containing less than 5% N were clustered into rated groups based on 90% homology over a continuous stretch of 90 nucleotides using the CAP3 assembler program . Sequences were then grouped into contigs and aligned. Contigs and singletons (contig containing only one sequence) were compared using the program BlastX, BlastN, or rpsBlast  to the non-redundant (nr) protein database of the National Center of Biological Information (NCBI), to the gene ontology database (GO) , to the Conserved Domains Database (CDD) that includes all Pfam , Smart [47, 48] and COG protein domains in the NCBI .
Additionally, contigs were compared to a customised subset of the NCBI nucleotide database containing either mitochondrial (mit-pla) or rRNA (rrna) sequences. Identification of putative secreted proteins was conducted using the SignalP server . The three-frame translation of each dataset was used to determine open reading frames (ORF). Only ORFs that started with a methionine and were longer than 40 amino acid (AA) residues were submitted to the SignalP server. The grouped and assembled sequences, BLAST results and signal peptide results were combined in an Excel spreadsheet and manually verified and annotated.
Sequence contigs containing signal peptides were selected from both species for further analysis and compared with secreted proteins from L. longipalpis  and P. ariasi . Full-length secreted proteins from each sandfly library were compared using a stand alone version of BLAST  creating formatted protein databases of each sandfly library. Related groups of proteins from each sandfly species were combined into defined families of proteins.
To ensure the fidelity of the sandfly library comparative analysis, the original nucleotide sequence data files from each library (including secreted and non-secreted sequences) were combined and compared to known families of salivary gland proteins using BLAST analysis. Results of this analysis were compared to the original individual analysis.
Predictions of protein secondary structures were performed using the PHYRE prediction server .
The evolutionary relatedness of the protein families identified through the bioinformatics analysis was evaluated using phylogenetics. Consensus protein sequences of the identified protein families from each of the sandflies used in this analysis were compared with related sequences from non-visceral Leishmania sandfly vectors as well as non-sandfly species obtained from GenBank. Sequences were aligned using ClustalX  and manually refined using BioEdit sequence-editing software . Alignments were analysed using ProtTest version 1.2.6  to determine the best fit model of protein evolution for each particular alignment. Phylogenetic analysis was conducted on protein alignments using Tree Puzzle version 5.2  incorporating the appropriate model of evolution defined by ProtTest. Tree Puzzle constructs phylogenetic trees by maximum likelihood using quartet puzzling, automatically estimating internal branch node support (1000 replications). Derived trees were visualised using TreeView .
Full-length Sequencing of Selected cDNA Clones
An aliquot (4 μl) of the λ-phage containing the cDNA of interestwas amplified using the PT2F1 and PT2R1 primers, as described previously . The PCR samples were cleaned using the multiscreen PCR 96-well filtration system (Millipore). Cleaned samples were sequenced first with PT2F3 primer (5'-TCT CGG GAA GCG CGC CAT TGT-3') and subsequently with custom primers.
For P. argentipes salivary glands, Tris-glycine gels (4–20%), 1 mm thick (Invitrogen), were used. Gels were run with Tris-glycine SDS buffer according to the manufacturer's instructions. To estimate the molecular weight of the samples, SeeBlue™ MW markers from Invitrogen were used. SGH were treated with equal parts of 2× SDS sample buffer (8% SDS in Tris-HCl buffer, 0.5 M, pH 6.8, 10% glycerol and 1% bromophenol blue dye). Each lane contained 20 pairs of homogenised P. argentipes salivary glands (20 ug protein). Protein were visualised with Coomassie blue stain. For aminoterminal sequencing of proteins, 20 pairs of homogenised salivary glands were electrophoresed and transferred to polyvinylidene difluoride (PVDF) membrane using 10 mM CAPS, pH 11.0, 10% methanol as the transfer buffer on a Blot-Module for the XCell II Mini-Cell (Invitrogen). The membrane was stained with Coomassie blue without acetic acid. Stained bands were cut from the PVDF membrane and subjected to Edman degradation using a Procise sequencer (Perkin-Elmer Corp.). Similar procedure was used to perform the amino-terminal sequence of P. perniciosus with the exception that 10% polyacylamide NuPAGE Bis-Tris (Invitrogen) was used for protein separation.
To locate the cDNA sequence cluster that corresponds to the amino acid sequence obtained by Edman degradation, we used a search program that compared the amino acid sequences against the three possible protein translations of each cDNA sequence obtained in the DNA sequencing project .
List of abbreviations
cluster of orthologous groups
phospholipase A2, MRJP, major royal jelly proteins
We want to thank Dr. Jose' M.C. Ribeiro and Dr. William C. Black, IV for critical review of this manuscript. We want to thank Dr. Robert Gwadz and Dr. Thomas Wellems for their continuous support and Nancy Shulman for editorial assistance. This project was funded by Grant 1Z01AI000932-02 to JGV from the Division of Intramural Research, National Institutes of Allergy and Infectious Diseases, National Institutes of Health.
- Guerin PJ, Olliaro P, Sundar S, Boelaert M, Croft SL, Desjeux P, Wasunna WK, Bryceson AD: Visceral leishmaniasis: current status of control, diagnosis, and treatment, and a proposed research and development agenda. The Lancet. 2002, 2: 494-501.PubMedView ArticleGoogle Scholar
- Ribeiro JMC, 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.PubMedView ArticleGoogle Scholar
- Belkaid Y, Kamhawi S, Modi G, Valenzuela JG, Noben-Trauth N, Rowton E, Ribeiro JMC, Sacks DL: Development of a natural model of cutaneous leishmaniasis: powerful effects of vector saliva and saliva preëxposure on the long-term outcome of Leishmania major infection in the mouse ear dermis. J Exp Med. 1998, 188: 941-953. 10.1084/jem.188.10.1941.View ArticleGoogle Scholar
- Kamhawi S, Belkaid Y, Modi G, Rowton E, Sacks D: Protection against cutaneous leishmaniasis resulting from bites of uninfected sand flies. Science. 2000, 290: 1351-1354. 10.1126/science.290.5495.1351.PubMedView ArticleGoogle Scholar
- Titus RG, Ribeiro JM: Salivary gland lysates from the sand fly Lutzomyia longipalpis enhance Leishmania infectivity. Science. 1988, 239: 1306-1308.PubMedView ArticleGoogle Scholar
- Morris RV, Shoemaker CB, David JR, Lanzaro GC, Titus RG: Sandfly maxadilan exacerbates infection with Leishmania major and vaccinating against it protects against L. major infection. J Immunol. 2001, 167: 5226-5230.PubMedView ArticleGoogle Scholar
- Valenzuela JG, Belkaid Y, Rowton E, Ribeiro JM: The salivary apyrase of the blood-sucking sand fly Phlebotomus papatasi belongs to the novel Cimex family of apyrases. J Exp Biol. 2001, 204: 229-237.PubMedGoogle Scholar
- Milleron RS, Mutebi JP, Valle S, Montoya A, Yin H, Soong L, Lanzaro GC: Antigenic diversity in maxadilan, a salivary protein from the sand fly vector of American visceral leishmaniasis. Am J Trop Med Hyg. 2004, 70: 286-293.PubMedGoogle Scholar
- Elnaiem DE, Meneses C, Slotman M, Lanzaro GC: Genetic variation in the sand fly salivary protein, SP-15, a potential vaccine candidate against Leishmania major. Insect Mol Biol. 2005, 14: 145-150. 10.1111/j.1365-2583.2004.00539.x.PubMedView ArticleGoogle Scholar
- Deissler H, Lottspeich F, Rajewsky MF: Affinity purification and cDNA cloning of rat neural differentiation and tumor cell surface antigen gp130RB13-6 reveals relationship to human and murine PC-1. J Biol Chem. 1995, 270: 9849-9855. 10.1074/jbc.270.17.9849.PubMedView ArticleGoogle Scholar
- Rossetto O, Rigoni M, Montecucco C: Different mechanism of blockade of neuroexocytosis by presynaptic neurotoxins. Toxicol Lett. 2004, 149: 91-101. 10.1016/j.toxlet.2003.12.023.PubMedView ArticleGoogle Scholar
- Aransay AM, Scoulica E, Tselentis Y, Ready PD: Phylogenetic relationships of phlebotomine sandflies inferred from small subunit nuclear ribosomal DNA. Insect Mol Biol. 2000, 9: 157-168. 10.1046/j.1365-2583.2000.00168.x.PubMedView ArticleGoogle Scholar
- Charlab R, Valenzuela JG, Rowton ED, Ribeiro JM: Toward an understanding of the biochemical and pharmacological complexity of the saliva of a hematophagous sand fly Lutzomyia longipalpis. Proc Natl Acad Sci USA. 1999, 96: 15155-15160. 10.1073/pnas.96.26.15155.PubMedPubMed CentralView ArticleGoogle Scholar
- Valenzuela JG, Belkaid Y, Garfield MK, Mendez S, Kamhawi S, Rowton ED, Sacks DL, Ribeiro JM: Toward a defined anti-Leishmania vaccine targeting vector antigens: characterization of a protective salivary protein. J Exp Med. 2001, 194: 331-342. 10.1084/jem.194.3.331.PubMedPubMed CentralView ArticleGoogle Scholar
- Galindo K, Smith DP: A large family of divergent Drosophila odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics. 2001, 159: 1059-1072.PubMedPubMed CentralGoogle Scholar
- Calvo E, Mans BJ, Andersen JF, Ribeiro JMC: Function and evolution of a mosquito salivary protein family. J Biol Chem. 2005, PMID: 16301315.,Google 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 USA. 1999, 96: 1516-1521. 10.1073/pnas.96.4.1516.PubMedPubMed CentralView ArticleGoogle Scholar
- James AA, Blackmer K, Marinotti O, Ghosn CR, Racioppi JV: Isolation and characterization of the gene expressing the major salivary gland protein of the female mosquito, Aedes aegypti. Mol Biochem Parasitol. 1991, 44: 245-253. 10.1016/0166-6851(91)90010-4.PubMedView ArticleGoogle Scholar
- Valenzuela JG, Francischetti IM, Pham VM, Garfield MK, Mather TN, Ribeiro JM: Exploring the sialome of the tick Ixodes scapularis. J Exp Biol. 2002, 205: 2843-2864.PubMedGoogle Scholar
- Valenzuela JG, Pham VM, Garfield MK, Francischetti IM, Ribeiro JM: Toward a description of the sialome of the adult female mosquito Aedes aegypti. Insect Biochem Mol Biol. 2002, 32: 1101-1122. 10.1016/S0965-1748(02)00047-4.PubMedView ArticleGoogle Scholar
- Oliveira F, Kamhawi S, Seitz AE, Pham VM, Guigal PM, Fischer L, Ward J, Valenzuela JG: From transcriptome to immunome: Identification of DTH inducing proteins from a Phlebotomus ariasi salivary gland cDNA library. Vaccine. 2005,Google Scholar
- Valenzuela JG, Charlab R, Gonzalez EC, de Miranda-Santos IK, Marinotti O, Francischetti IM, Ribeiro JM: The D7 family of salivary proteins in blood sucking diptera. Insect Mol Biol. 2002, 11: 149-155. 10.1046/j.1365-2583.2002.00319.x.PubMedView ArticleGoogle Scholar
- Isawa H, Yuda M, Orito Y, Chinzei Y: Mosquito salivary protein inhibits activation of the plasma contact system by binding to factor XII and high molecular weight kininogen. J Biol Chem. 2002, 277: 27651-27658. 10.1074/jbc.M203505200.PubMedView ArticleGoogle Scholar
- Valenzuela JG, Ribeiro JM: Purification and cloning of the salivary nitrophorin from the hemipteran Cimex lectularius. J Exp Biol. 1998, 201: 2659-2664.PubMedGoogle Scholar
- Smith TM, Hicks-Berger CA, Kim S, Kirley TL: Cloning, expression, and characterization of a soluble calcium-activated nucleotidase, a human enzyme belonging to a new family of extracellular nucleotidases. Arch Biochem Biophys. 2002, 406: 105-115. 10.1016/S0003-9861(02)00420-4.PubMedView ArticleGoogle Scholar
- Murphy DM, Ivanenkov VV, Kirley TL: Bacterial expression and characterization of a novel, soluble, calcium-binding, and calcium-activated human nucleotidase. Biochemistry. 2003, 42: 2412-2421. 10.1021/bi026763b.PubMedView ArticleGoogle Scholar
- Grimaldi D, Engel MS: Evolution of the insects. Evolution of Ectoparasites and Blood Feeders of Vertebrates. 2005, Cambridge University Press. New York, 489-496.Google Scholar
- Champagne DE, Smartt CT, Ribeiro JMC, James AA: The salivary gland-specific apyrase of the mosquito Aedes aegypti is a member of the 5'-nucleotidase family. Proc Natl Acad Sci USA. 1995, 92: 694-698.PubMedPubMed CentralView ArticleGoogle Scholar
- Sarkis JJ, Guimaraes JA, Ribeiro JMC: Salivary apyrase of Rhodnius prolixus. Kinetics and purification. Biochem J. 1986, 233: 885-891.PubMedPubMed CentralView ArticleGoogle Scholar
- Dai J, Liu J, Deng Y, Smith TM, Lu M: Structure and protein design of a human platelet function inhibitor. Cell. 2004, 116: 649-659. 10.1016/S0092-8674(04)00172-2.PubMedView ArticleGoogle Scholar
- Geyer PK, Spana C, Corces VG: On the molecular mechanism of gypsy-induced mutations at the yellow locus of Drosophila melanogaster. EMBO J. 1986, 5: 2657-2662.PubMedPubMed CentralGoogle Scholar
- Johnson JK, Li J, Christensen BM: Cloning and characterization of a dopachrome conversion enzyme from the yellow fever mosquito, Aedes aegypti. Insect Biochem Mol Biol. 2001, 31: 1125-1135. 10.1016/S0965-1748(01)00072-8.PubMedView ArticleGoogle Scholar
- Valenzuela JG, Garfield M, Rowton ED, Pham VM: Identification of the most abundant secreted proteins from the salivary glands of the sand fly Lutzomyia longipalpis, vector of Leishmania chagasi. J Exp Biol. 2004, 207: 3717-3729. 10.1242/jeb.01185.PubMedView ArticleGoogle Scholar
- Gomes RB, Brodskyn C, de Oliveira CI, Costa J, Miranda JC, Caldas A, Valenzuela JG, Barral-Netto M, Barral A: Seroconversion against Lutzomyia longipalpis saliva concurrent with the development of anti-Leishmania chagasi delayed-type hypersensitivity. J Infect Dis. 2002, 186 (10): 1530-1534. 10.1086/344733.PubMedView ArticleGoogle Scholar
- Lu G, Villalba M, Coscia MR, Hoffman DR, King TP: Sequence analysis and antigenic cross-reactivity of a venom allergen, antigen 5, from hornets, wasps, and yellow jackets. J Immunol. 1993, 150: 2823-2830.PubMedGoogle Scholar
- Francischetti IM, Valenzuela JG, Pham VM, Garfield MK, Ribeiro JM: Toward a catalog for the transcripts and proteins (sialome) from the salivary gland of the malaria vector Anopheles gambiae. J Exp Biol. 2002, 205: 2429-2451.PubMedGoogle Scholar
- Fang KS, Vitale M, Fehlner P, King TP: cDNA cloning and primary structure of a white-face hornet venom allergen, antigen 5. Proc Natl Acad Sci USA. 1988, 85: 895-899.PubMedPubMed CentralView ArticleGoogle Scholar
- Asojo OA, Goud G, Dhar K, Loukas A, Zhan B, Deumic V, Liu S, Borgstahl GE, Hotez PJ: X-ray structure of Na-ASP-2, a pathogenesis-related-1 protein from the nematode parasite, Necator americanus, and a vaccine antigen for human hookworm infection. J Mol Biol. 2005, 346: 801-814. 10.1016/j.jmb.2004.12.023.PubMedView ArticleGoogle Scholar
- Li S, Kwon J, Aksoy S: Characterization of genes expressed in the salivary glands of the tsetse fly, Glossina morsitans morsitans. Insect Mol Biol. 2001, 10: 69-76. 10.1046/j.1365-2583.2001.00240.x.PubMedView ArticleGoogle 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: 543-563. 10.1016/j.ibmb.2004.02.008.PubMedView ArticleGoogle Scholar
- Rohousova I, Ozensoy S, Ozbel Y, Volf P: Detection of species-specific antibody response of humans and mice bitten by sand flies. Parasitology. 2005, 130: 493-499. 10.1017/S003118200400681X.PubMedView ArticleGoogle Scholar
- Volf P, Rohousova I: Species-specific antigens in salivary glands of phlebotomine sandflies. Parasitology. 2001, 122: 37-41. 10.1017/S0031182000007046.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang M, Madan A: CAP3: A DNA sequence assembly program. Genome Research. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
- Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2000, 28: 263-266. 10.1093/nar/28.1.263.PubMedPubMed CentralView ArticleGoogle Scholar
- Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci USA. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.PubMedPubMed CentralView ArticleGoogle Scholar
- Schultz J, Copley RR, Doerks T, Ponting CP, Bork P: SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000, 28: 231-234. 10.1093/nar/28.1.231.PubMedPubMed CentralView ArticleGoogle Scholar
- Marchler-Bauer A, Panchenko AR, Ariel N, Bryant SH: Comparison of sequence and structure alignments for protein domains. Proteins. 2002, 48: 439-446. 10.1002/prot.10163.PubMedView ArticleGoogle 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: 581-599. 10.1142/S0129065797000537.PubMedView ArticleGoogle Scholar
- BLAST. [ftp://ftp.ncbi.nih.gov/blast/executables/]
- PHYRE. [http://www.sbg.bio.ic.ac.uk/~phyre]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL-X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralView ArticleGoogle Scholar
- Bioedit. [http://www.mbio.ncsu.edu/BioEdit/page2.html]
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.PubMedView ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar