Analysis of salivary transcripts and antigens of the sand fly Phlebotomus arabicus
© Hostomská et al; licensee BioMed Central Ltd. 2009
Received: 17 March 2009
Accepted: 25 June 2009
Published: 25 June 2009
Sand fly saliva plays an important role in blood feeding and Leishmania transmission as it was shown to increase parasite virulence. On the other hand, immunity to salivary components impedes the establishment of infection. Therefore, it is most desirable to gain a deeper insight into the composition of saliva in sand fly species which serve as vectors of various forms of leishmaniases. In the present work, we focused on Phlebotomus (Adlerius) arabicus, which was recently shown to transmit Leishmania tropica, the causative agent of cutaneous leishmaniasis in Israel.
A cDNA library from salivary glands of P. arabicus females was constructed and transcripts were sequenced and analyzed. The most abundant protein families identified were SP15-like proteins, ParSP25-like proteins, D7-related proteins, yellow-related proteins, PpSP32-like proteins, antigen 5-related proteins, and 34 kDa-like proteins. Sequences coding for apyrases, hyaluronidase and other putative secreted enzymes were also represented, including endonuclease, phospholipase, pyrophosphatase, amylase and trehalase. Mass spectrometry analysis confirmed the presence of 20 proteins predicted to be secreted in the salivary proteome. Humoral response of mice bitten by P. arabicus to salivary antigens was assessed and many salivary proteins were determined to be antigenic.
This transcriptomic analysis of P. arabicus salivary glands is the first description of salivary proteins of a sand fly in the subgenus Adlerius. Proteomic analysis of P. arabicus salivary glands produced the most comprehensive account in a single sand fly species to date. Detailed information and phylogenetic relationships of the salivary proteins are provided, expanding the knowledge base of molecules that are likely important factors of sand fly-host and sand fly-Leishmania interactions. Enzymatic and immunological investigations further demonstrate the value of functional transcriptomics in advancing biological and epidemiological research that can impact leishmaniasis.
Phlebotomine sand flies are the arthropod vectors of Leishmania parasites, the causative agents of leishmaniasis. During the feeding process sand flies inject saliva into the site of the bite to facilitate successful acquisition of a blood meal . An infected sand fly regurgitates infective metacyclic promastigote stage Leishmania while feeding; thus, the parasite is always introduced to the host as a mixture with sand fly saliva. Sand fly saliva facilitates the transmission of Leishmania parasites to mammalian hosts; at the same time, immune response to salivary components was shown to partially protect the host from Leishmania infection . Therefore, salivary components essential for parasite transmission and/or eliciting protective immune response are sought-after. Salivary proteins from Phlebotomus papatasi, the vector of Leishmania major, and Lutzomyia longipalpis, the vector of L. infantum, have been extensively studied. In addition, cDNA libraries from several other sand fly species were characterized and include other sand flies that vector L. major (P. duboscqi), L. infantum (P. ariasi and P. perniciosus) and L. donovani (P. argentipes).
Phlebotomus (Adlerius) arabicus is distributed in certain parts of East Africa and the Middle East. In Ethiopia, P. arabicus females infected with uncharacterized Leishmania sp. were reported , and in northern Israel P. arabicus is the proven vector of cutaneous leishmaniasis caused by L. tropica [4, 5]. Cutaneous leishmaniasis caused by L. tropica is found in a vast discontinuous area reaching from the south-western Mediterranean to Turkey, north-western India and Sub-Saharan Africa . Long supposed to circulate in anthroponotic foci exclusively, L. tropica was recently shown to occur as an anthropozoonosis as well [4, 5]. Laboratory experiments demonstrated that P. arabicus is a permissive vector, meaning it is susceptible to development of more than one species of Leishmania, including L. major and L. infantum .
In the present study, salivary gland transcripts and proteins of P. (Adlerius) arabicus were studied by cDNA sequencing, electrophoretic and proteomic methods. This is the first study of the repertoire of salivary molecules of a vector of L. tropica and it is the first report of the composition of salivary proteins in the subgenus Adlerius.
Results and Discussion
Sequencing of salivary gland cDNA library
A cDNA library was constructed from salivary glands of Phlebotomus arabicus females dissected one day after emergence. From this cDNA library, 1152 random transcripts were selected and sequenced, resulting in 985 high quality sequences. Sequences were clustered together based on sequence homology and produced 107 clusters and 288 sequences were assessed as singletons (only one sequence). The term cluster may refer to either singletons or multiple homologous sequences. Similar to other sand flies studied so far, the most abundant transcripts were those coding for putative secretory proteins and resulted in 74 clusters with an average number of 7.65 sequences per cluster. Predicted proteins containing retention signals for endoplasmic reticulum and/or transmembrane domains were not treated as putatively secreted. An example of such proteins is the translocon-associated protein complex, δ subunit (PabSP91; GenBank accession number FJ427208), which has homologs previously designated as 16 kDa or 16.1 kDa salivary protein in P. ariasi or L. longipalpis, respectively.
Members of 21 different families were found among putative secretory proteins. BLAST comparison of translated nucleotide sequences with the NR protein database showed that overall, high similarity was observed namely with salivary proteins of L. infantum vectors P. (Larroussius) ariasi and P. (L.) perniciosus. The expected values of these matches were highly significant at values lower than E-60. To a lesser extent, similarity to sequences of salivary proteins of P. (Euphlebotomus) argentipes, the vector of L. donovani in India, was also observed. These findings are fully in concert with the close evolutionary relationship of Larroussius and Adlerius subgenera reported by Aransay et al. .
Some of the protein families contained multiple members. The observed variability among individual protein family members might be explained by intraspecific polymorphism, as the tested sample was heterogeneous (a pool of salivary glands from 30 female sand flies). Nevertheless, analysis of genetic variation of SP15 salivary protein in P. papatasi brought strong evidence that SP-15 is a multicopy gene . While individual intraspecific variability of sand fly salivary proteins awaits broader analysis, we propose that the multiple homologous transcripts within protein families observed in this P. arabicus salivary gland cDNA library may reflect gene duplication events or allelic variation.
Most abundant salivary protein transcripts from Phlebotomus arabicus
Best match to NR protein database
GenBank accession #
Seq per cluster
14 kDa salivary protein A
D7-related salivary protein A
27 kDa salivary protein similar to ParSP25
yellow-related salivary protein
antigen 5-related protein
25 kDa salivary protein A similar to ParSP25
25 kDa salivary protein similar to ParSP02
putative salivary phospholipase A2
putative salivary endonuclease
46 kDa salivary protein
34 kDa-like salivary protein B
14 kDa salivary protein B
26 kDa salivary protein similar to ParSP02
salivary apyrase A
2.7 kDa salivary peptide
34 kDa-like salivary protein A
28 kDa salivary protein similar to ParSP02
D7-related salivary protein C
D7-related salivary protein B
25 kDa salivary protein B similar to ParSP25
25 kDa salivary protein C similar to ParSP25
salivary apyrase C
10 kDa salivary protein similar to ParSP13
16 kDa salivary protein A
14 kDa salivary protein C
16 kDa salivary protein B
10 kDa salivary protein similar to ParSP21
25 kDa salivary protein D similar to ParSP25
27 kDa salivary protein B similar to ParSP25
salivary apyrase B
D7-related salivary protein D
22 kDa salivary protein
17 kDa salivary protein
27 kDa and 25 kDa proteins
Six clusters coding for proteins related to P. ariasi 27 kDa salivary protein (ParSP25; GenBank accession number AAX55664) and P. perniciosus 29 kDa salivary protein (PpeSP08, GenBank accession number ABA43056) were found in the P. arabicus salivary gland cDNA library. There are no other homologs of these proteins in accessible databases, no conserved domains were found in the translated sequences, and no function has been assigned to these proteins. However, in P. arabicus cDNA library they represent the second most abundant protein family. Transcripts coding for ParSP25-like proteins occurred in long (PabSP15 and PabSP11; GenBank accession numbers FJ538100 and FJ538101) and short forms (PabSP14, PabSP16, PabSP13 and PabSP12; GenBank accession numbers FJ538102, FJ538103, FJ538104 and FJ538105, respectively), with very little variability among individual clusters. The mature proteins coded by these transcripts have a predicted Mw of 27 kDa and 25 kDa, respectively, and are composed of 90.5% identical amino acids and 0.9% similar amino acids. The predicted pI of the proteins is acidic (average pI = 5.03), unlike most sand fly salivary proteins described thus far.
Yellow-related proteins are common in insects; in bloodsucking Diptera, yellow-related proteins were described from mosquitoes and sand flies. In Ae. aegypti, a dopachrome-converting enzyme shares homology with Drosophila melanogaster yellow proteins  and, according to Li et al. , it might play a role in melanotic encapsulation of parasites in the hemocoel. In sand fly salivary gland samples, however, dopachrome-converting enzyme activity could not be detected (Hostomska, unpublished observations), while yellow protein of P. duboscqi was detected in midgut and salivary glands and shown to have lectin properties . Sand fly yellow proteins were previously proposed as potential antigens recognized by sera of experimentally bitten mice and dogs, and naturally exposed humans [24–26]. In L. longipalpis this was also suggested by mass spectrometry . In the P. arabicus salivary gland cDNA library a single homolog of yellow-related proteins was found (PabSP26; GenBank accession number FJ410293). The predicted Mw of the protein is 42.9 kDa with a pI of 8.4. No N-glycosylation sites were predicted in the protein sequence by amino acid submission to the NetNGlyc server.
Antigen 5-related proteins
Apyrases are widespread in saliva of bloodsucking insects. The antihemostatic effects of saliva are, for a great part, due to apyrase anti-platelet activity . Sand fly apyrases belong to the Cimex apyrase family . Three very similar apyrase clusters coding for apyrases were found in P. arabicus cDNA library (98.4% identity) (PabSP41, PabSP40 and PabSP39; GenBank accession numbers EZ000631, EZ000632 and EZ000633, respectively). The predicted average pI of P. arabicus apyrases is 8.85 and the predicted Mw is 35.3 kDa.
A cluster encoding a putative endonuclease was identified in the P. arabicus cDNA library (PabSP49; GenBank accession number FJ439531). Similar sequences were reported from P. ariasi, P. perniciosus, P. argentipes, and L. longipalpis salivary glands. Cluster PabSP49 encodes an endonuclease domain, which is typical for DNA/RNA non-specific endonucleases. Since all residues composing the active site, the substrate binding site and the Mg2+ binding site are conserved in this cluster; we suggest that PabSP49 might possess endonuclease activity. The predicted pI of the mature protein is 9.45 and the predicted Mw of the mature protein is 40.5 kDa. Possible roles for a salivary endonuclease include reducing the viscosity of the blood pool during feeding and liberating nucleosides. Exogenous nucleosides, primarily adenosine, can exhibit regulatory effects on blood clotting, immune and inflammatory responses, and Leishmania pathogenesis .
Additional putative enzymes
In the amino acid translation of sequence cluster 52 (PabSP52, GenBank accession number EZ000627), a phospholipase A2 (PLA2) domain is present, containing all conserved residues of both catalytic and metal-binding sites of PLA2. In hymenopteran venoms, PLA2 represents a major allergen. In the salivary glands of blood-feeding insects, sequences coding for PLA2-like proteins were reported only from sand flies of the subgenus Larroussius [10, 30]. We tested salivary gland samples of P. arabicus specifically for PLA2 activity and did not detect any positive reaction. Cluster 52 contains an exceptionally long 5' untranslated region (5' UTR) compared to other clusters from this cDNA library coding for secreted proteins. The 5' UTR in this cluster is more than 500 nucleotides long. Thus, the regulation of expression of this transcript might be different from other transcripts reported herein.
Other sequences coding for other putative enzymes could not be obtained as full-length clones. These included a pyrophosphatase-like protein (PabSP288, GenBank accession number EZ000634), amylase (PabSP47, GenBank accession number EZ000626), an enzyme involved in digestion of dietary starch , and trehalase (PabSP315, GenBank accession number EZ000625). Previously, sequences coding for pyrophosphatase-like proteins were reported in P. argentipes and P. duboscqi sand flies [10, 37]. These proteins, as well as their P. arabicus homolog reported herein, contain a conserved phosphodiesterase domain, typical for enzymes cleaving phosphodiester and phosphosulphate bonds in NAD, deoxynucleotides and nucleotide sugars . Transcripts coding for α-amylase were found in L. longipalpis salivary glands and midguts as well as P. papatasi midguts [28, 39, 40]. Amylase activity was shown in L. longipalpis and P. papatasi salivary gland samples [28, 36] and it is likely that the enzymatic activity is present also in P. arabicus salivary glands. The putative trehalase enzyme from P. arabicus salivary glands might either be an intrinsic component of insect metabolism, or might be related to sugar feeding and digestion. Trehalose is the main energy source in insect hemolymph in general. Trehalases are involved in its hydrolysis, yielding glucose molecules which are then readily available to various cells of the insect body . So far, trehalase enzyme or sequence has not been reported from salivary glands of any blood-feeding insects, but sequences coding for sand fly trehalase have been found in midgut cDNA libraries of P. papatasi .
Putative secreted proteins of unknown function
There were a number of transcripts with no homology to known enzymes or structural proteins; however, eight of these transcripts encode potentially secreted proteins with high homology to other sand fly salivary molecules. P. arabicus salivary transcripts code for 34 kDa proteins homologous to ParSP09. Polymorphisms resulting in different translations of the transcripts were observed (PabSP32 and PabSP34; GenBank accession numbers FJ489241 and FJ489242, respectively). These proteins are seemingly sand fly-specific; apart from 5 sand fly species (P. ariasi, P. perniciosus, P. argentipes, P. duboscqi and L. longipalpis); no other related proteins from any organism were reported. Another family of putative sand fly-specific proteins contain homologs of P. arabicus 46 kDa salivary protein (PabSP53; GenBank accession number FJ538106). Transcripts coding for such proteins were reported from P. perniciosus, P. ariasi, P. duboscqi and L. longipalpis. Homologs of P. arabicus 2.7 kDa peptide (PabSP56; GenBank accession number FJ538099) were only found in P. ariasi and P. perniciosus [10, 30]. Our finding contradicts the suggestion that 2.7 kDa peptides are specific for the subgenus Larroussius . Likewise, two unrelated 10 kDa proteins (PabSP75 and PabSP97; GenBank accession numbers FJ474087 and FJ474088, respectively) were found in P. arabicus; homologous molecules were identified in P. ariasi . Additionally, two transcripts putatively encoding 16 kDa salivary proteins (two polymorphic variants PabSP63 and PabSP64; GenBank accession numbers FJ474085 and FJ474086, respectively) are homologous to molecules identified in P. argentipes .
Two sequence clusters coding for putatively secreted proteins in the P. arabicus cDNA library show no similarity with known sand fly sequences. Cluster 107 (PabSP107; GenBank accession number EZ000635) is homologous to Ae. aegypti putative salivary secreted mucin 3, as well as the IgE binding protein icarapin from honeybee venom . The predicted molecular weight of the mature protein is 22.2 kDa and the putative protein would have an acidic pI of 4.4. There are 2 potential N-glycosylation sites and 9 potential O-glycosylation sites in the amino acid sequence of cluster 107, as predicted by submission to the NetNGlyc and NetOGlyc servers. Similarly, putative extracellular proteins of Anopheles gambiae (XP_001230739) and Aedes aegypti (XP_001650286) were also predicted to contain multiple O-glycosylation sites. These proteins might be involved in hypersensitivity to bites of blood-feeding insects. The second cluster 126 (PabSP126; GenBank accession number EZ000636), encodes a homolog of conserved hypothetical proteins of culicine as well as anopheline mosquitoes. The predicted molecular weight of the mature protein from P. arabicus is 17.2 kDa and the predicted pI is 5.34. No N- or O-glycosylation sites were predicted in cluster 126 protein and nothing is known about these hypothetical proteins.
Proteome analysis of P. arabicus salivary glands
For the proteome analysis, P. arabicus salivary gland samples separated by SDS-PAGE were subjected to Edman degradation and mass spectrometry. Edman degradation resulted in the identification of 7 different N-terminal sequences. These were representative of two 14 kDa proteins (PabSP2 and PabSP45; GenBank accession numbers FJ538111 and FJ538112, respectively), yellow-related protein (PabSP26; FJ410293), and endonuclease (PabSP49; FJ439531). An N-terminal sequence common to all six variants of salivary proteins similar to ParSP25 was also detected by Edman degradation (PabSP11-16; GenBank accession numbers FJ538100–FJ538105), as well as N-terminal sequences common to apyrases (PabSP39-41; EZ000631-EZ000633) and to D7-related proteins A and C (PabSP20 and PabSP54; FJ538107 and FJ538109, respectively). From the data obtained by Edman degradation analysis it could not be concluded which variants of polymorphic salivary proteins were present in the proteome.
Humoral response to P. arabicus saliva
In this study we generated a transcriptome of female sand fly Phlebotomus arabicus salivary glands using a PCR-based cDNA library. This is the first reported salivary gland transcriptome of a sand fly from the subgenus Adlerius. The most abundant transcripts were represented in the 985 high quality sequences. Many of the transcripts encoded full- or partial-length proteins; most of which are homologous to other sand fly species saliva molecules. Phylogenetic analysis consistently shows a strong relationship between P. arabicus with sand flies from the Larroussius subgenus; specifically, P. ariasi and P. perniciousus. The phylogenetic analyses of sand fly salivary proteins reaffirm the taxonomy of phlebotomines .
Some of the most abundant molecules identified in the transcriptome that have a predicted signal secretion peptide include a 14 kDa protein (PabSP2), a D7-related protein (PabSP20), a yellow-related protein (PabSP26), an Antigen 5-related protein (PabSP4) and 25 kDa and 27 kDa proteins similar to P. ariasi ParSP25 (PabSP14 and PabSP15, respectively). A number of paralogous transcripts were identified, such as those in the SP15 and D7 families. The presence of duplicate gene copies has been observed in other blood feeding arthropods [10, 44] and can serve several potential functions including increased transcript abundance and rapid evolution of blood feeding strategies while retaining intrinsic proteins. Proteomic analysis by N-terminal sequencing or tryptic digestion followed by mass spectrometry identified 20 proteins in the salivary gland homogenate of P. arabicus that were characterized in the transcriptome. In addition, one protein was identified by mass spectrometry that did not match any of the characterized transcripts. This is the most comprehensive description of sand fly salivary proteome to date and also demonstrates that the transcriptome represents >95% of the most abundant proteins present in the salivary gland.
In the analysis of the P. arabicus salivary gland transcriptome four sequences were identified as encoding a putative hyaluronidase. Hyaluronidase is an enzyme that has been identified in a number of phlebotomine salivary glands including Lutzomyia longipalpis, P. Phlebotomus papatasi, P. Phlebotomus duboscqi, P. Paraphlebotomus sergenti and P. Adlerius halepensis . The zymographic analysis of salivary gland extract confirms the presence of an active hyaluronidase enzyme and demonstrates the effectiveness of a transcriptomic approach to identifying disease vector salivary components.
Having described the repertoire of saliva molecules opens more doors in the research of vector-host and vector-parasite interactions, pharmacology and insect biochemistry. The antigenicity of sand fly saliva is one important aspect of the vector-host interaction. Evaluating the humoral response of mice bitten by P. arabicus demonstrates the abundance and diversity of antigenic molecules in the saliva. Future work may focus on the use of functional transcriptomics (expression of recombinant protein and biological and biochemical assays) to use these P. arabicus salivary proteins to evaluate the role of these molecules in the epidemiology of leishmaniasis.
Sand flies and salivary gland dissection
The colony of P. arabicus (Israel) was reared in the insectary of Charles University in Prague in standard conditions as described by Benkova and Volf . For mRNA extraction salivary glands of 1-day-old females were dissected in saline and stored in RNA later (Ambion). For proteome analysis and Western blot analysis, salivary glands from 5- to 7-day-old P. arabicus females were dissected and stored in Tris buffer (20 mM Tris, 150 mM NaCl, pH 7.5).
Construction of salivary gland cDNA library
Salivary gland mRNA was isolated from 30 pairs of glands using Micro-FastTrack mRNA isolation kit (Invitrogen). PCR-based cDNA library was made following the manufacturer's instructions for SMART™ cDNA library Construction Kit (BD Clontech) with some modifications described by Chmelar et al. . The cDNA library was fractionated into three sets of cDNAs containing large, medium and small fragments. Gigapack® III Gold Packaging Extract (Stratagene) was used for packaging the phage particles. The libraries were plated by infecting log-phase XL-1 blue Escherichia coli (Clontech). Several plaques from each plate were selected and a PCR with vector primers flanking the inserted cDNA was performed. The presence of recombinants was checked by visualisation the PCR products on 1.1% agarose gel with ethidium bromide.
Sequencing of Selected cDNA Clones
Plaques were randomly selected from the plated libraries and transferred to 96-well polypropylene plate containing 75 μl of water per well. The PCR reaction amplifying randomly selected cDNAs was performed using FastStart PCR Master mix (Roche), 3 μl of the phage sample as a template and primers described elsewhere . Amplification conditions were as follows: 1 hold of 75°C for 3 min, 1 hold of 94°C for 2 min and 34 cycles of 94°C for 1 min, 49°C for 1 min and 72°C for 2 min. Final elongation step lasted for 10 min at 72°C. Reaction products were cleaned using ExcelaPure 96-Well UF PCR Purification Plates (EdgeBio) and used as templates for cycle-sequencing reaction using BigDye Terminator v3.1 cycle sequencing kit (Applied Biosystems) and a forward primer described elsewhere . Cycle-sequencing reaction products were cleaned using sephadex and MultiScreen HV Plates (Millipore), dried and stored at -20°C. Sequencing was performed on an ABI 3730Xl DNA sequencer (Applied Biosystems).
Detailed description of the bioinformatic treatment of the data can be found elsewhere [29, 46]. Briefly, EST trace files were analyzed using a customized program based on the Phred algorithm [47, 48]. Sequences with Phred quality scores lower than 20 were removed, as well as primer and vector sequences. Resulting sequences were grouped into clusters using a customized program based on identity (95% identity, 64 word size) and aligned into contiguous sequences (contigs) using the CAP3 sequence assembly program . BLASTX, BLASTN or RPS BLAST programs  were used to compare contigs and singletons (contigs with a single sequence) to the non-redundant (NR) protein database of the NCBI, the gene ontology (GO) fasta subset , to the conserved domains database (CDD) of NCBI  which contains KOG (), Pfam  and Smart databases, and to mitochondrial and rRNA nucleotide sequences available from NCBI. The three frame translations of each dataset were submitted to the SignalP server  to detect signal peptides. The grouped and assembled sequences, BLAST results and SignalP results were combined in an Excel spreadsheet and manually verified and annotated. N- and O-glycosylation site prediction was performed for selected sequences using NetNGlyc 1.0 and NetOGlyc 3.1 software (www.cbs.dtu.dk/services/NetNGlyc, www.cbs.dtu.dk/services/NetOGlyc) .
Protein sequences of the members of identified protein families were compared with related sequences of other sand fly species obtained from GenBank. Sequences were aligned using ClustalW version 1.4  running under BioEdit sequence-editing software, version 7, and manually refined in BioEdit. For each alignment, best substitution matrix was determined by ProtTest software, version 1.4 . This matrix was then used by TREE-PUZZLE 5.2  to reconstruct phylogenetic trees from the protein alignments by maximum likelihood. TREE-PUZZLE implements quartet puzzling (QP) tree search; at the same time, the algorithm estimates support values for each internal branch. The number of puzzling steps was 1000 in each phylogenetic analysis. Resulting trees were viewed in MEGA 4 .
Salivary glands from 5-day-old P. arabicus females were homogenized by 5 freeze-thaw cycles. Samples were reduced using sample buffer with 2-mercaptoethanol, and electrophoretically separated in 12% polyacrylamide SDS gel. Gels were stained for total proteins with Coomassie G-250 (SimplyBlue SafeStain, Invitrogen) or for glycoproteins with Pro-Q Emerald 300 glycoprotein stain (Invitrogen). Mass spectrometric analysis was performed with individual bands cut from the Coomassie-stained gel. The individual bands were placed in microtubes and covered with 100 μl 50 mM ammonium bicarbonate (ABC) buffer in 50% acetonitrile (ACN) with 50 mM dithiothreitol (DTT). The samples were subjected to sonication in an ultrasonic bath for 5 minutes. After 15 minutes the supernatant was discarded and the gel was covered with 100 μl of 50 mM ABC/50% ACN with 50 mM iodoacetamide and sonicated for 5 minutes. After 25 minutes, the supernatant was discarded and exchanged for 100 μl 50 mM ABC/50% ACN with 50 mM DTT and sonicated for 5 minutes to remove any excess iodoacetamide. The supernatant was discarded and samples were sonicated for 5 minutes in 100 μl of HPLC water. The water was discarded and samples were sonicated for another 5 minutes in 100 μl of ACN. ACN was discarded and microtubes with samples were left open for a couple of minutes to allow the rest of ACN to evaporate. Five ng of trypsin (Promega) in 10 μl of 50 mM ABC were added to the gel. Samples were incubated at 37°C overnight. Trifluoroacetic acid (TFA) and ACN were added to reach final concentration 1% TFA, 30% ACN. Samples were sonicated for 10 minutes and 0.5 μl drop was transferred onto MALDI target and let to dry. Dried droplets were covered with 0.5 μl drop of alpha-cyano-hydroxycinnamic acid solution (2 mg/ml in 80% ACN) and let to dry. Samples were measured using a 4800 Plus MALDI TOF/TOF analyzer (Applied Biosystems/MDS Sciex) equipped with a Nd:YAG laser (355 nm, firing rate 200 Hz).
Peak lists from the MS spectra were generated by 4000 Series Explorer V 3.5.3 (Applied Biosystems/MDS Sciex) without smoothing, peaks with local signal to noise ratio greater than 5 were picked and searched by local Mascot v. 2.1 (Matrix Science) against a database of proteins sequences derived from cDNA library. Database search criteria were as follows – enzyme: trypsin, taxonomy: none, fixed modification: carbamidomethylation, variable modification: methionine oxidation, peptide mass tolerance: 120 ppm, one missed cleavage allowed. Only hits that were scored as significant (p < 0.0001) are included.
For the Edman degradation analysis, Phlebotomus arabicus salivary glands were electrophoretically separated on 1 mm thick 4–20% NuPAGE Novex Bis-Tris gels using MES SDS buffer (Invitrogen). A sample containing 30 glands was reduced with NuPAGE Sample Reducing Agent (Invitrogen) and run in parallel with non-reduced samples (50 glands) on the same gel. Wet blotting on a PVDF membrane was performed using XCell II™ Blot Module (Invitrogen). SeeBlue® Pre-Stained Standards (Invitrogen) were used to estimate molecular weight (Mw) of separated proteins and assess transfer efficiency. The membrane was stained with 0.025% Coomassie blue without acetic acid. Stained bands were cut and subjected to Edman degradation using a Procise 494cLC sequencer (Applied Biosystems). cDNA sequences corresponding to obtained N-terminal amino acid sequences of salivary proteins were identified using an in-house search program . This program compared three possible translations of each cDNA sequence obtained in the P. arabicus cDNA sequencing project with the amino acid sequences.
Salivary gland samples from 5-day-old P. arabicus females were tested for the activities of hyaluronidase and phospholipase A2. Salivary glands were dissected in Tris buffer (20 mM Tris, 150 mM NaCl, pH 7.8) and stored at -20°C. Before use, the glands were mechanically disrupted, samples were centrifuged at 12000 g for 5 minutes, and the supernatant was used in the assays.
The zymographic hyaluronidase assay was performed as described by Volfova et al. . Briefly, salivary gland samples were separated by SDS electrophoresis in 8% polyacrylamide slab gels with incorporated hyaluronan (0.002%). Prior to electrophoresis, aliquots of the sample were subjected to treatment with reducing agents (dithiothreitol 50 mM, 45 minutes at 25°C, or 2-mercaptoethanol 7 mM, 45 minutes at 40°C). An equivalent of 1 salivary gland per lane was used for electrophoresis under reducing conditions using dithiothreitol, an equivalent of 2 salivary glands per lane for electrophoresis under reducing conditions using 2-mercaptoethanol, and for electrophoresis under non-reducing conditions, an equivalent of 1/12 of a gland per lane was used. After the electrophoresis, SDS was washed out of the gels; the gels were equilibrated with 0.1 M acetate buffer and incubated at 37°C for 120 minutes. Gels stained with Stains-all (Sigma) were rinsed in distilled water and documented.
The phospholipase A2 activity was tested in the same salivary gland samples using the EnzChek® Phospholipase A2 Assay Kit (Invitrogen). The manufacturer's instructions were followed. Samples containing 0.5, 1, 2, 5, 10 and 20 salivary glands in 50 μl were assayed; as a positive control, PLA2 from honey bee venom supplied with the kit was used. The assay was performed at 26°C, and fluorescence emission at 515 nm was detected after an incubation step of 12 minutes (Tecan Infinite 200).
Immunization of mice
Experiments on mice were done in accordance with Czech Act No. 246/1992 and approved by IACUC of the Faculty of Science, Charles University in Prague. Female BALB/c mice, 8 weeks old (Charles River Deutschland, Sulzfeld, Germany) were used for the exposure to P. arabicus sand flies. In biweekly intervals, the mice were anaesthetized with ketamine (1.5 mg/10 g body weight) and xylazin (0.15 mg/10 g body weight) and exposed to sand fly females. Each time, an average of 60 females fed on each mouse (SE = 14.9). Ten days after the sixth exposure the mice were bled and obtained sera were kept at -20°C.
Phlebotomus arabicus salivary gland proteins were disrupted by three freeze-thaw cycles in liquid nitrogen and separated by SDS-PAGE on 12% polyacrylamide gel, 0.75 mm thick using Mini-Protean III apparatus (BioRad). As a reducing agent, 2.5% 2-mercaptoethanol was used in sample buffer. Biotinylated low range standards (BioRad) were run on the same gel. Separated proteins were electro-transferred onto a nitrocellulose membrane by iBlot™ Gel Transfer Device (Invitrogen), and the membrane was cut into strips corresponding to sample load of 5 salivary glands per strip. These strips carrying salivary gland samples were either stained with amidoblack (Merck; 0.1% solution in 25% isopropanol and 10% acetic acid) or blocked with 5% BSA in phosphate-buffered saline, pH 7.4, with 0.05% Tween-20 (PBS-Tw) overnight. The sample-carrying strips were incubated for 1 hour with pre-immune or immune sera of BALB/c mice exposed to P. arabicus (diluted 1:250 in PBS-Tw), and then with horseradish peroxidase-conjugated swine anti-mouse antibody (Sevapharma, 1:1000 in PBS-Tw) for 45 minutes. Streptavidin-conjugated peroxidase (Sigma, 1 μg/ml in PBS-Tw) was used for the biotinylated standards. The colour reaction was developed using H2O2 and diaminobenzidine in PBS. Washes between individual steps were done with PBS-Tw.
We thank Petr Jedelský for help with the mass spectrometry analysis, Dr. Helena Kulíková and Dušan Eremiáš for helpful technical and administrative assistance. Participation of Iva Rohoušová was supported by the Fulbright Scholar Program. The study was partially supported by projects MSM 0021620828, LC 06009, and The Division of Intramural Research, The National Institute of Allergy and Infectious Diseases.
- Ribeiro JM: Role of saliva in blood-feeding by arthropods. Annu Rev Entomol. 1987, 32: 463-478. 10.1146/annurev.en.32.010187.002335.View ArticlePubMedGoogle Scholar
- Rohousova I, Volf P: Sand fly saliva: effects on host immune response and Leishmania transmission. Folia Parasitol (Praha). 2006, 53 (3): 161-171.View ArticleGoogle Scholar
- Gebre-Michael T, Balkew M, Ali A, Ludovisi A, Gramiccia M: The isolation of Leishmania tropica and L. aethiopica from Phlebotomus (Paraphlebotomus) species (Diptera: Psychodidae) in the Awash Valley, northeastern Ethiopia. Trans R Soc Trop Med Hyg. 2004, 98 (1): 64-70. 10.1016/S0035-9203(03)00008-7.View ArticlePubMedGoogle Scholar
- Jacobson RL, Eisenberger CL, Svobodova M, Baneth G, Sztern J, Carvalho J, Nasereddin A, El Fari M, Shalom U, Volf P, et al: Outbreak of cutaneous leishmaniasis in northern Israel. J Infect Dis. 2003, 188 (7): 1065-1073. 10.1086/378204.View ArticlePubMedGoogle Scholar
- Svobodova M, Votypka J, Peckova J, Dvorak V, Nasereddin A, Baneth G, Sztern J, Kravchenko V, Orr A, Meir D, et al: Distinct transmission cycles of Leishmania tropica in 2 adjacent foci, Northern Israel. Emerg Infect Dis. 2006, 12 (12): 1860-1868.PubMed CentralView ArticlePubMedGoogle Scholar
- Jacobson RL: Leishmania tropica (Kinetoplastida: Trypanosomatidae) – a perplexing parasite. Folia Parasitol (Praha). 2003, 50 (4): 241-250.View ArticleGoogle Scholar
- Myskova J, Svobodova M, Beverley SM, Volf P: A lipophosphoglycan-independent development of Leishmania in permissive sand flies. Microbes Infect. 2007, 9 (3): 317-324. 10.1016/j.micinf.2006.12.010.PubMed CentralView ArticlePubMedGoogle 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 (2): 157-168. 10.1046/j.1365-2583.2000.00168.x.View ArticlePubMedGoogle 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 (2): 145-150. 10.1111/j.1365-2583.2004.00539.x.View ArticlePubMedGoogle 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 CentralView ArticlePubMedGoogle 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 (3): 331-342. 10.1084/jem.194.3.331.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersen JF, Pham VM, Meng Z, Champagne DE, Ribeiro JM: Insight into the Sialome of the Black Fly, Simulium vittatum. J Proteome Res. 2009, 8 (3): 1474-1488. 10.1021/pr8008429.PubMed CentralView ArticlePubMedGoogle Scholar
- Campbell CL, Vandyke KA, Letchworth GJ, Drolet BS, Hanekamp T, Wilson WC: Midgut and salivary gland transcriptomes of the arbovirus vector Culicoides sonorensis (Diptera: Ceratopogonidae). Insect Mol Biol. 2005, 14 (2): 121-136. 10.1111/j.1365-2583.2004.00537.x.View ArticlePubMedGoogle 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 (2): 245-253. 10.1016/0166-6851(91)90010-4.View ArticlePubMedGoogle 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 (2): 149-155. 10.1046/j.1365-2583.2002.00319.x.View ArticlePubMedGoogle Scholar
- Calvo E, Mans BJ, Andersen JF, Ribeiro JM: Function and evolution of a mosquito salivary protein family. J Biol Chem. 2006, 281 (4): 1935-1942. 10.1074/jbc.M510359200.View ArticlePubMedGoogle Scholar
- Isawa H, Orito Y, Iwanaga S, Jingushi N, Morita A, Chinzei Y, Yuda M: Identification and characterization of a new kallikrein-kinin system inhibitor from the salivary glands of the malaria vector mosquito Anopheles stephensi. Insect Biochem Mol Biol. 2007, 37 (5): 466-477. 10.1016/j.ibmb.2007.02.002.View ArticlePubMedGoogle Scholar
- Mans BJ, Calvo E, Ribeiro JM, Andersen JF: The crystal structure of D7r4, a salivary biogenic amine-binding protein from the malaria mosquito Anopheles gambiae. J Biol Chem. 2007, 282 (50): 36626-36633. 10.1074/jbc.M706410200.View ArticlePubMedGoogle 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. 10.1016/S0091-6749(98)70357-4.View ArticlePubMedGoogle Scholar
- Bahia D, Gontijo NF, Leon IR, Perales J, Pereira MH, Oliveira G, Correa-Oliveira R, Reis AB: Antibodies from dogs with canine visceral leishmaniasis recognise two proteins from the saliva of Lutzomyia longipalpis. Parasitol Res. 2007, 100 (3): 449-454. 10.1007/s00436-006-0307-8.View ArticlePubMedGoogle 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 (11): 1125-1135. 10.1016/S0965-1748(01)00072-8.View ArticlePubMedGoogle Scholar
- Li J, Zhao X, Christensen BM: Dopachrome conversion activity in Aedes aegypti: significance during melanotic encapsulation of parasites and cuticular tanning. Insect Biochem Mol Biol. 1994, 24 (10): 1043-1049. 10.1016/0965-1748(94)90142-2.View ArticlePubMedGoogle Scholar
- Volf P, Skarupova S, Man P: Characterization of the lectin from females of Phlebotomus duboscqi sand flies. Eur J Biochem. 2002, 269 (24): 6294-6301. 10.1046/j.1432-1033.2002.03349.x.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Hostomska J, Rohousova I, Volfova V, Stanneck D, Mencke N, Volf P: Kinetics of canine antibody response to saliva of the sand fly Lutzomyia longipalpis. Vector Borne Zoonotic Dis. 2008, 8 (4): 443-450. 10.1089/vbz.2007.0214.View ArticlePubMedGoogle 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 (Pt 5): 493-499. 10.1017/S003118200400681X.View ArticlePubMedGoogle 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 (3): 895-899. 10.1073/pnas.85.3.895.PubMed CentralView ArticlePubMedGoogle 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 (26): 15155-15160. 10.1073/pnas.96.26.15155.PubMed CentralView ArticlePubMedGoogle 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 (9): 1101-1122. 10.1016/S0965-1748(02)00047-4.View ArticlePubMedGoogle 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. 2006, 24 (3): 374-390. 10.1016/j.vaccine.2005.07.085.View ArticlePubMedGoogle Scholar
- Champagne DE: Antihemostatic molecules from saliva of blood-feeding arthropods. Pathophysiol Haemost Thromb. 2005, 34 (4–5): 221-227. 10.1159/000092428.View ArticlePubMedGoogle 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
- Sansom FM, Robson SC, Hartland EL: Possible effects of microbial ecto-nucleoside triphosphate diphosphohydrolases on host-pathogen interactions. Microbiol Mol Biol Rev. 2008, 72 (4): 765-781. 10.1128/MMBR.00013-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Cerna P, Mikes L, Volf P: Salivary gland hyaluronidase in various species of phlebotomine sand flies (Diptera: psychodidae). Insect Biochem Mol Biol. 2002, 32 (12): 1691-1697. 10.1016/S0965-1748(02)00109-1.View ArticlePubMedGoogle Scholar
- Volfova V, Hostomska J, Cerny M, Votypka J, Volf P: Hyaluronidase of bloodsucking insects and its enhancing effect on leishmania infection in mice. PLoS Negl Trop Dis. 2008, 2 (9): e294-10.1371/journal.pntd.0000294.PubMed CentralView ArticlePubMedGoogle Scholar
- Jacobson RL, Schlein Y: Phlebotomus papatasi and Leishmania major parasites express alpha-amylase and alpha-glucosidase. Acta Trop. 2001, 78 (1): 41-49. 10.1016/S0001-706X(00)00164-9.View ArticlePubMedGoogle Scholar
- Kato H, Anderson JM, Kamhawi S, Oliveira F, Lawyer PG, Pham VM, Sangare CS, Samake S, Sissoko I, Garfield M, et al: High degree of conservancy among secreted salivary gland proteins from two geographically distant Phlebotomus duboscqi sandflies populations (Mali and Kenya). BMC Genomics. 2006, 7: 226-10.1186/1471-2164-7-226.PubMed CentralView ArticlePubMedGoogle 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 (17): 9849-9855. 10.1074/jbc.270.17.9849.View ArticlePubMedGoogle Scholar
- Jochim RC, Teixeira CR, Laughinghouse A, Mu J, Oliveira F, Gomes RB, Elnaiem DE, Valenzuela JG: The midgut transcriptome of Lutzomyia longipalpis: comparative analysis of cDNA libraries from sugar-fed, blood-fed, post-digested and Leishmania infantum chagasi-infected sand flies. BMC Genomics. 2008, 9: 15-10.1186/1471-2164-9-15.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramalho-Ortigao M, Jochim RC, Anderson JM, Lawyer PG, Pham VM, Kamhawi S, Valenzuela JG: Exploring the midgut transcriptome of Phlebotomus papatasi: comparative analysis of expression profiles of sugar-fed, blood-fed and Leishmania-major-infected sandflies. BMC Genomics. 2007, 8: 300-10.1186/1471-2164-8-300.PubMed CentralView ArticlePubMedGoogle Scholar
- Clements AN: The biology of mosquitoes. development, nutrition and reproduction. 1992, Chapman & Hall; London, 1: 1Google Scholar
- Peiren N, de Graaf DC, Brunain M, Bridts CH, Ebo DG, Stevens WJ, Jacobs FJ: Molecular cloning and expression of icarapin, a novel IgE-binding bee venom protein. FEBS Lett. 2006, 580 (20): 4895-4899. 10.1016/j.febslet.2006.08.005.View ArticlePubMedGoogle Scholar
- Volf P, Tesarova P, Nohynkova EN: Salivary proteins and glycoproteins in phlebotomine sandflies of various species, sex and age. Med Vet Entomol. 2000, 14 (3): 251-256. 10.1046/j.1365-2915.2000.00240.x.View ArticlePubMedGoogle Scholar
- Ribeiro JM, Alarcon-Chaidez F, Francischetti IM, Mans BJ, Mather TN, Valenzuela JG, Wikel SK: An annotated catalog of salivary gland transcripts from Ixodes scapularis ticks. Insect Biochem Mol Biol. 2006, 36 (2): 111-129. 10.1016/j.ibmb.2005.11.005.View ArticlePubMedGoogle Scholar
- Benkova I, Volf P: Effect of temperature on metabolism of Phlebotomus papatasi (Diptera: Psychodidae). J Med Entomol. 2007, 44 (1): 150-154. 10.1603/0022-2585(2007)44[150:EOTOMO]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Chmelar J, Anderson JM, Mu J, Jochim RC, Valenzuela JG, Kopecky J: Insight into the sialome of the castor bean tick, Ixodes ricinus. BMC Genomics. 2008, 9: 233-10.1186/1471-2164-9-233.PubMed CentralView ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.View ArticlePubMedGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.View ArticlePubMedGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle 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 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Panchenko AR, Shoemaker BA, Thiessen PA, Geer LY, Bryant SH: CDD: a database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res. 2002, 30 (1): 281-283. 10.1093/nar/30.1.281.PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, et al: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41-10.1186/1471-2105-4-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2000, 28 (1): 263-266. 10.1093/nar/28.1.263.PubMed CentralView ArticlePubMedGoogle 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 (1): 231-234. 10.1093/nar/28.1.231.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Julenius K, Molgaard A, Gupta R, Brunak S: Prediction, conservation analysis, and structural characterization of mammalian mucin-type O-glycosylation sites. Glycobiology. 2005, 15 (2): 153-164. 10.1093/glycob/cwh151.View ArticlePubMedGoogle Scholar
- Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene. 1988, 73 (1): 237-244. 10.1016/0378-1119(88)90330-7.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18 (3): 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle 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.