BACKGROUND
most adult female mosquitoes are hematophagous, in addition to taking sugar meals.
saliva helps blood feeding by interfering with host reactions that could disrupt the blood flow, assists sugar meals with glycosidases, and contain antimicrobials that may control microbial growth in their meals  <cit> .
with the advent of tissue transcriptomics, we can postulate that these functions are mediated by large numbers  of polypeptides, many of which are expressed solely in the adult female salivary glands  <cit> .
unique protein families have been found in anopheles, aedes or culex mosquitoes, as well as a group of common proteins or enzymes  <cit> .
functional characterization of these proteins uncovers scavengers of biogenic amines  <cit>  or leukotrienes  <cit> , inhibitors of blood clotting  <cit> , bradykinin formation  <cit> , platelet aggregation  <cit>  and vasodilators  <cit> .
other molecularly uncharacterized activities include inhibitors of mast cell tnf production  <cit>  and inhibition of t cell activation  <cit> .
it is apparent that the complexity of the salivary components affecting host hemostasis and inflammation mirrors the complexity of host hemostasis and inflammation itself, which must be disarmed for successful blood feeding.
indeed mosquitoes lacking salivation by salivary duct ablation feed less and take more dangerous time of exposure to their hosts  <cit> .
perhaps due to the potent pharmacological activities of saliva, or the immune reactions to it, mosquito saliva also plays a role in pathogen transmission, including in arboviral transmission  <cit> .
accordingly, determination of the salivary composition of vector mosquitoes not only discovers new potentially pharmacologically active molecules, but also can also help generating vaccine targets for disruption of arboviral transmission.
these proteins may also be of epidemiological significance as selective human or animal markers of vector exposure  <cit> .
we have previously described the sialotranscriptome of c. quinquefasciatus  <cit> , where several unique protein families were discovered, many of which are products of gene duplications.
indications of horizontal gene transfer from bacteria to mosquitoes were also pointed out as participating in the generation of mosquito sialomes  <cit> .
to this date, no other sialotranscriptome from a member of the culex genus has been described.
we currently portray the sialotranscriptome of culex tarsalis, a north american bird feeding mosquito  <cit> , and, like c. pipiens, a good vector of west nile virus  <cit> .
RESULTS
characteristics of the assembled salivary est set
a total of  <dig>  cdna clones were used to assemble a database  that yielded  <dig> clusters of related sequences,  <dig> of which contained only one expressed sequence tag .
the  <dig> clusters were compared, using the programs blastx, blastn, or rpsblast  <cit> , to the nonredundant protein database of the national center of biological information , to a gene ontology database  <cit>  to the conserved domains database of the ncbi  <cit> , and to a custom prepared subset of the ncbi nucleotide database containing either mitochondrial or rrna sequences.
three categories of expressed genes derived from the manual annotation of the contigs .
the putatively secreted  category contained 44% of the sequences, the housekeeping  category had  <dig> % and  <dig> % of the ests could not be classified and belong to the unknown  class.
the transcripts of the unknown class could represent novel proteins or derive from the less conserved 3' or 5' untranslated regions of genes, as was indicated for the sialotranscriptome of an.
gambiae  <cit> .
this distribution is typical of previous mosquito sialotranscriptomes  <cit> .
housekeeping  genes
the  <dig> ests attributed to h genes expressed in the salivary glands of c. tarsalis were further characterized into  <dig> subgroups according to function .
transcripts associated with the protein synthesis machinery represented 44% of all transcripts associated with a housekeeping function, an expected result for the secretory nature of the organ.
nineteen percent of the transcripts were classified as either "unknown conserved" or "conserved secreted proteins".
these represent highly conserved proteins of unknown function, presumably associated with cellular function, but still uncharacterized.
these sets may help functional identification of the "conserved hypothetical" proteins as previously reviewed in  <cit> .
transporters accounted for 7% and included several members of the vacuolar atpase machinery, important for mosquito salivary secretion  <cit> .
energy metabolism and protein modification related enzymes accounted for nearly 5% each of the transcripts of the h class.
possibly secreted  class of expressed genes
a total of  <dig> ests represent putative c. tarsalis salivary components .
these include ubiquitous known gene families, e.g, families seen in other organs and organisms, as well as protein families unique to hematophagous diptera, or to culicines, or to culex genus alone.
table  <dig> also indicates our degree of knowledge, or ignorance, regarding these protein families,  <dig> of which we have no hint for their function, more than half of the  <dig> classes of proteins represented in table  <dig>  many of these putatively secreted protein families of unknown function are multigenic, such as the cwrp family found previously only in c. quinquefasciatus  <cit> , which contributes to 58% of the ests associated with the s class in c. tarsalis.
the d7/obp-like families contribute to 9% of all transcripts associated with secreted products, which is less than previously seen in anopheles or aedes sialotranscriptomes, probably due to the dilution effect following abundant expression of the cwrp family.
* y for yes and n for no - yes is probably an overestimation as proteins can be multi functional
the salivary secretome of culex tarsalis
from the sequenced cdnas, a total of  <dig> novel c. tarsalis protein sequences were derived,  <dig> of which code for putative secreted products .
because cdna sequences coding for many of these proteins were found as singletons, this secretome is to be considered preliminary and incomplete, but many parallels with previous mosquito sialotranscriptomes can be drawn, as follows:
proteins with presumed or experimentally validated function
the d7/odorant binding protein-like family
the d <dig> proteins constitute a unique multi gene family found expressed in the salivary glands of blood sucking diptera  <cit> , belonging to the superfamily of odorant binding proteins   <cit> .
long and short versions exist, the long versions containing two and the short versions containing one obp domain.
three genes codes for long d <dig> proteins and  <dig> genes code for short d <dig> proteins in an.
gambiae  <cit> , while in ae.
aegypti  <dig> genes code for long and  <dig> genes code for short d <dig> proteins  <cit> .
short versions of anopheline d <dig> mosquitoes were shown to bind biogenic amines such as serotonin and histamine  <cit> , thus counteracting hemostasis and inflammation.
more recently, the amino terminal obp domain of a d <dig> long form of ae.
aegypti was shown to bind peptidic leukotrienes with high affinity.
the crystal structures of a short d <dig> protein from an.
gambiae and a long d <dig> protein from ae.
aegypti revealed that the d <dig> obp domains have  <dig> alpha helices,  <dig> more than the canonical obp family  <cit> .
in addition to these inflammatory agonist binding functions, a short d <dig> protein from an.
stephensi, named hamadarin, was shown to inhibit bradykinin formation by inhibiting the fxii/kallikrein pathway  <cit> .
eight members of the d <dig> protein family were found in the c. tarsalis sialotranscriptome, in addition to two canonical obp proteins.
alignment of these proteins with their close blast matches found in the nr protein database, and construction of a bootstrapped phylogram  shows several clades with strong bootstrap support, and allows the following inferences: the eight c. tarsalis proteins are the probable results of  <dig> genes.
ctar- <dig>  ctar- <dig> and ctar- <dig> are possible alleles of a single gene related to other long c. quinquefasciatus d <dig> proteins shown in clade i, while ctar- <dig>  ctar- <dig> and ctar- <dig> are probably alleles of a second c. tarsalis gene coding for yet another long d <dig> protein, and are associated within clade ii.
ctar- <dig>  results from a gene coding for a third long d <dig> protein grouping in clade v with other mosquito long d <dig> proteins, including, anopheles and aedes proteins.
ctar- <dig> is the only c. tarsalis sequence for a short d <dig> protein, grouping in clade iv with other culex pipiens d <dig> proteins.
notice also clades iii and clade vi, which exclusively contain long d <dig> anopheline or phlebotomine proteins, respectively.
ctar- <dig>  ctar- <dig>  ctar- <dig> and ctar- <dig> possess the signature - -x-c-x-w-x-w-x- -x-c- -x- -c-x-q-x-c-x-  found in lipid binding d <dig> domain of ae.
aegypti  <cit> .
the serotonin binding motif found on short anopheline d <dig> proteins as well as the carboxydomain of the long d <dig> protein of a. aegypti  <cit>  was not found on any c. tarsalis d <dig> protein , indicating this motif may have evolved beyond recognition, or that other proteins in culex may undertake this task .
ctar- <dig> and ctar- <dig> represent  <dig> canonical obp protein sequences, transcripts of which were found expressed in the salivary glands of c. tarsalis.
homologs of these canonical obp proteins were also previously found expressed in c. quinquefasciatus and an.
gambiae sialotranscriptomes.
these proteins tend to be much more conserved in sequence to mosquito and even drosophila proteins and it is possible that they may play an endogenous or housekeeping function.
the  <dig> kda antigen/ge rich/aegyptin family
this protein family is typical of the salivary glands of adult female mosquitoes, and was first identified as a salivary antigen in ae.
aegypti  <cit> , and later found in salivary transcriptomes and proteomes of both culicine and anopheline mosquitoes  <cit> , where it was named ge rich protein.
distantly related proteins also exist in black flies and sand flies.
more recently, proteins of this family from aedes and anopheles were shown to prevent platelet aggregation by collagen  <cit> , indicating conservation of function after the split of the culicidae into the culicines and anophelines, approximately  <dig> mya  <cit> .
the sialotranscriptome of c. tarsalis allowed the identification of  <dig> protein sequences from the  <dig> kda antigen/aegyptin family , all represented by 1- <dig> ests found in the library.
blast searches against the nr database retrieved closely related sequences.
those that were more than 95% identical and of the same species were excluded to avoid inclusion of alleles or very closely related genes.
clustal alignment of this subset was used to produce a bootstrapped phylogram  that allows the following inferences: the  <dig> proteins from c. tarsalis most probably derive from  <dig> genes, ctar- <dig>  ctar- <dig> and ctar- <dig> being probable alleles of a gene closely related to a c. quinquefasciatus salivary gene shown in clade i. clade i includes aedes aegypti salivary proteins that are short versions of the canonical  <dig> kda protein, as are all culex proteins in this clade.
a second c. tarsalis gene possibly codes for ctar- <dig>  ctar- <dig> and ctar- <dig>  and a third gene for ctar- <dig> and ctar- <dig>  these two gene products groups with strong bootstrap support to their c. quinquefasciatus homologs, and these in turn group relatively weakly  in clade ii as shown in fig  <dig>  fig  <dig> also shows the uniquely aedes containing proteins in clade iii, represented by pairs of homologs between ae.
aegypti and ae.
albopictus, and the uniquely anopheline clade iv, with one sequence each of  <dig> different mosquito species.
this phylogram suggests the single gene status of this family in the anopheles genus, and the multi gene character in culicines.
at least  <dig> genes for canonical  <dig> kda ag/aegyptin exists in the aedes genus, and at least one more gene for the shorter protein shown in clade i. culex also have the shorter version gene, plus at least  <dig> genes of the canonical type.
mucins and peritrofins
serine and threonine rich proteins are a common finding in sialotranscriptomes.
these proteins are generally modified post-translationally and their mature forms have n-acetyl galactosamine residues, typical of mucins  <cit> .
they probably have a function to lubricate the food canals and may also have antimicrobial function.
ctar- <dig> and ctar- <dig> encode related protein sequences that might reflect splice variants of the same gene.
these  <dig> truncated protein sequences have over  <dig> predicted glycosylation sites, and are similar to a previously described salivary mucin of c. quinquefasciatus.
ctar- <dig> is related to aedes and culex mucins that, in the case of the aedes protein, were associated with induction by viral infection suggesting an immune function for this protein with  <dig> putative glycosylation sites.
ctar- <dig> encodes a 5' truncated protein sequence similar to salivary mucins of c. quinquefasciatus and aedes albopictus.
enzymes
additional file  <dig>  table s <dig> indicate the presence of transcripts coding for enzymes possibly associated with the sugar meal and blood meal.
related to the sugar meal are several transcripts coding for glycosidases similar to proteins annotated as maltase and amylase.
additional file  <dig>  table s <dig> provides one full length sequence of a salivary alpha-glucosidase similar to other culicine salivary maltases, plus  <dig> other truncated sequences coding for different sugar hydrolases.
regarding blood meal-related enzymes,  <dig> est's coding for fragments of adenosine deaminase were found, as well as for an endonuclease.
in c. quinquefasciatus a salivary endonuclease has been previously characterized and associated with helping forming the feeding hematoma  <cit> .
however, the c. tarsalis enzyme is most closely related to sand flies salivary endonucleases and only distantly related to the salivary endonuclease of c. quinquefasciatus, although it is closely related to a non-salivary c. quinquefasciatus enzyme, indicating this c. tarsalis enzyme may be playing a housekeeping function.
remarkably absent from the c. tarsalis sialotranscriptome are ests coding for members of the 5' nucleotidase, which functions as the salivary apyrase of mosquitoes.
apyrase hydrolyzes atp and adp to amp and orthophosphate, destroying these important agonists of inflammation and platelet aggregation  <cit> .
it has been previously noticed that c. quinquefasciatus has little salivary apyrase activity when compared to aedes and anopheles mosquitoes, and this observation was postulated to be a consequence of the lack of platelets in birds, the most common host of culex  <cit> , while aedes and anopheles are mostly mammalian feeders  <cit> .
this may explain the absence of 5'-nucleotidase/apyrase coding transcripts in c. tarsalis, although an increased sequencing effort could produce 5'-nucleotidases that might be secreted.
as an example of the abundance of this type of transcript in mosquito sialotranscriptomes, only one transcript out of  <dig> clones coded for apyrase in the c. quinquefasciatus sialotranscriptome  <cit> , while  <dig> out of  <dig>  were found in a similarly made library from an.
gambiae, and  <dig> out of  <dig>  ests were found in the ae.
aegypti sialotranscriptome  <cit> .
from these observed frequencies, the expected number of est's for c. tarsalis and c. quinquefasciatus are  <dig> and  <dig>  respectively, producing a chi square of  <dig>  and a p <  <dig> , indicating the low expression of apyrase transcripts in culex.
three ests coding for a phosphatase produce the predicted truncated sequence encoded by ctar- <dig>  which is similar to a salivary alkaline phosphatase previously described in ae.
aegypti.
the function of this enzyme in feeding, if any, is unknown.
transcripts coding for at least  <dig> different serine proteases were found in c. tarsalis sialotranscriptome .
these enzymes may function in immunity as prophenoloxidase activators, or in digesting skin matrix components, such as in an elastase function, or hydrolysing host blood clotting enzymes such as fibrinogen/fibrin, or activating plasminogen.
immunity related products
antimicrobial peptides, lysozyme, and pathogen pattern recognition polypeptides are commonly found in the sialotranscriptome of blood sucking arthropods.
additional file  <dig>  table s <dig> shows the full length sequence of c. tarsalis salivary lysozyme, which is 91% identical to the c. quinquefasciatus homolog, and 75% identical to the salivary homolog of ae.
albopictus.
truncated orf's of a c-type lectin and a gram-negative binding protein were also found.
they both match previously described salivary proteins of aedes and culex.
secreted proteins with unknown function
promiscuous antigen  <dig>  family
this is a ubiquitous protein family found in animals and plants  <cit> , and in all sialotranscriptomes of blood sucking diptera analyzed so far.
the function of these proteins in mosquito saliva is not known, but in blood sucking brachycera two proteins of this family have been functionally characterized.
remarkably, in a tabanid fly, a member of the ag- <dig> family acquired a typical rgd domain surrounded by cys residues and acts as a main platelet aggregation  <cit> , and in the stable fly a salivary ag- <dig> protein binds immunoglobulins and may function as an inhibitor of the classical complement pathway  <cit> .
we present evidence, in the form of truncated transcripts, for the expression of at least two members of the family in c. tarsalis salivary glands; ctar- <dig>  assembled from  <dig> ests, matches with 94% identity the salivary secreted antigen- <dig> precursor ag5- <dig> from culex quinquefasciatus while ctar- <dig> matches with 91% identity another c. quinquefasciatus protein of the same family.
insect specific cys-rich polypeptide family
two transcripts were found in the c. tarsalis sialotranscriptome coding for basic  peptides of mature mw of  <dig>  kda containing  <dig> cys residues.
this peptide family was previously found in the salivary transcriptome of c. quinquefasciatus, but close relatives of the same size exist in drosophila, bombyx, tribolium and apis.
the ubiquity of this protein family in insects, its size and pi suggests an antimicrobial role.
hematophagous diptera specific  <dig>  kda family
the first  <dig>  kda family member was characterized in the sialotranscriptome of ae.
aegypti, and later found in c. quinquefasciatus and in ae.
albopictus  <cit> .
although not present in the sialotranscriptomes of members of the anopheline cellia subgenus an.
funestus, an.
stephensi and an.
gambiae , it was recently found in an.
darlingi, a member of the nyssorhynchus subgenus  <cit> , characterizing this family as uniquely culicidae.
additional file  <dig>  table s <dig> provides evidence of a member of this protein family in c. tarsalis, encoded by ctar- <dig>  producing a predicted protein of mature mw of  <dig> kda, being 69% identical to its c. quinquefasciatus homolog.
psiblast of ctar- <dig> against the nr protein database retrieves on its first blastp cycle only mosquito salivary proteins, as expected.
on the second cycle, it retrieves with lower significance  mostly bacterial proteins, but also salivary proteins from simulium vittatum and culicoides sonorensis.
further iterations of psiblast retrieves with high significance bacterial proteins of the methyl-accepting chemotaxis receptor family , which may suggest this bacterial family to be originated from horizontal gene transfer to an ancestral blood feeding nematocera.
the gene structure of this protein in c. quinquefasciatus and ae.
aegypti is similar, containing  <dig> exons with a short intron of ca.
<dig> nt.
inclusion of simulium and culicoides sequences in this family indicates an ancient origin before the nematocera split  <cit> .
mosquito specific  <dig>  kda protein
ctar- <dig> codes for a protein related to mosquito proteins so far only found expressed in mosquito adult salivary glands.
psiblast of ctar- <dig> against the nr database does not retrieve any additional protein with an e value better than  <dig> .
mosquito specific hyp <dig> family
this member of this mosquito salivary protein family was identified in an.
stephensi  <cit> , and later found also an.
gambiae  <cit> , but not previously in other mosquito transcriptomes, including c. quinquefasciatus and ae.
aegypti.
ctar- <dig>  represented by a single est in our database, codes for a member of this family as indicated by blastp comparisons to the nr protein database, where it retrieves only culex and anopheles proteins.
psiblast of ctar- <dig> against the nr database does not increase finding matches beyond additional mosquito proteins, including proteins deducted from the uncovering of the c. quinquefasciatus genome with and without signal peptide, and additional proteins from an.
gambiae lacking a signal peptide indicative of secretion.
interestingly, the anopheles gambiae gene found expressed in the salivary glands was shown to reside as a single exon in chromosome arm 2r  <cit> .
the an.
gambiae gene was also found to be selectively expressed in the adult female salivary glands suggesting a role in blood feeding  <cit> .
culicine specific  <dig>  kda family
members of the  <dig>  kda protein were previously discovered in ae.
aegypti, ae.
albopictus and c. quinquefasciatus sialotranscriptomes.
additional related proteins were also deducted from the genomes of ae.
aegypti and c. quinquefasciatus, but not an.
gambiae.
the sialotranscriptome of c. tarsalis provides additional evidence for this multi copy family exclusive of culicines.
clustal alignment of these protein sequences  shows many conserved amino acids and a conserved cys framework on the second half of the protein.
the bootstrapped phylogram shows  <dig> robust clades.
clade i contains a sub clade of  <dig> ae.
aegypti and  <dig> ae.
albopictus proteins, and a second sub clade of a c. quinquefasciatus and ctar- <dig>  clade ii has  <dig> c. tarsalis sequences, possibly alleles of a single gene or closely related genes, and one c. quinquefasciatus sequence.
clade iii shows a possible gene expansion in c. quinquefasciatus containing  <dig> gene products, clustering with one ae.
aegypti sequence.
overall the phylogram indicates that aedes aegypti has at least  <dig> genes coding for this protein family, while c. quinquefasciatus has at least  <dig>   <dig> of which possibly arrived by further gene duplications within clade iii.
culex tarsalis and ae.
albopictus, have at least two genes expressing these proteins in their adult female salivary glands.
psiblast of members of this protein family against the nr protein database converges after  <dig> iterations, retrieving solely aedes and culex proteins .
the function of this protein family is unknown, but transcripts for this family were found enriched in the salivary glands of adult female ae.
aegypti  <cit>  suggesting a blood feeding role.
culicine specific  <dig>  kda family
ctar- <dig> codes for a mature peptide of predicted mw =  <dig>  kda and pi =  <dig> .
it is similar to previously described salivary peptides found in sialotranscriptomes of ae.
aegypti and ae.
albopictus, but not culex.
psiblast using an inclusion e value of  <dig>   against the nr protein database additionally identifies a salivary peptide previously found in c. quinquefasciatus, annotated as  <dig>  kda salivary peptide.
inclusion of this peptide with the previous to generate the reverse position search model  <cit>  additionally identifies other salivary peptides from c. quinquefasciatus, with convergence of the search and identification of the  <dig> non redundant peptides with an e value of  <dig> e- <dig> or smaller.
alignment of these  <dig> peptides with ctar- <dig> shows that the aedes sequences are more compact than those from culex, containing a conserved framework of  <dig> cysteines, while those of culex have  <dig> additional cys residues .
only  <dig> other residues are identical in all sequences, indicating a high evolutionary rate of this peptide family, if they indeed share a common ancestor.
the bootstrapped phylogram  shows  <dig> robust clades, clade i containing only aedes sequences in two subclades, each containing a pair of ae.
aegypti or ae.
albopictus sequences.
clade ii has  <dig> sequences from c. quinquefasciatus.
ctar- <dig> does not cluster with either of the clades, indicating again the possible fast evolutionary rate of this peptide family.
culicine specific  <dig>  kda peptide family
ctar- <dig>  ctar- <dig> and ctar- <dig> encode peptides with similarities to c. quinquefasciatus  <dig>  kda salivary peptide, which has similarity solely to another salivary peptide previously found in ae.
albopictus.
culex genus specific cysteine and tryptophan rich protein family
this protein family was previously discovered in the c. quinquefasciatus sialotranscriptome  <cit> , where  <dig> proteins were deduced from the est sequences.
these have a signal peptide indicative of secretion,  <dig> conserved cysteins,  <dig> conserved tryptophans and one conserved asparagine and leucine residues.
inspection of the c. quinquefasciatus genome indicates a total of  <dig> genes coding for this protein family, most  members being single exonic.
expression of these genes in c. quinquefasciatus salivary glands was confirmed by edman degradation of electrophoretically separated salivary homogenate  <cit> .
the c. tarsalis sialotranscriptome led to identification of  <dig> full length members of this protein family, plus  <dig> other fragments.
alignment of the members of this family  shows a conserved cysteine framework in most members, two conserved tryptophan residues , as well at the conserved asparagine and leucine residues.
this protein family is by far the most transcribed in c. tarsalis adult female salivary glands, accounting for 25% of the transcripts that were sequenced.
interestingly, the most transcribed member in c. tarsalis, ctar- <dig>  is a truncated member of the family by relatively recent acquisition of a stop codon just after the second conserve cysteine.
this abundant expression suggests that this family may be responsible for host serotonin or histamine antagonism, as is performed by the abundantly expressed d <dig> proteins of aedes or culex.
this possibility is further supported by the relatively low expression of the d <dig> protein family in c. tarsalis and their absent serotonin binding motifs which are found both in aedes and anopheles d <dig> proteins .
psiblast of members of the cwrp protein family retrieves, interestingly, sugar binding proteins including bacterial glycosidases, lectins such as ricin and hemolysins, the crystal structure of which have a trefoil structure  <cit> .
the single exon structure of this protein family in c. quinquefasciatus suggests either acquition of this gene family by horizontal transfer, and/or gene duplication occurring by retroposition of an rna template.
culex genus specific  <dig>  peptide family
ctar- <dig> and ctar- <dig> encode truncated peptides matching the previously described putative  <dig>  kda salivary peptide of c. quinquefasciatus, which has no other matches to known proteins.
culex genus specific gqp repeat protein family
ctar- <dig>  ctar- <dig> and ctar- <dig> produce matches to a low complexity protein previously described in the c. quinquefasciatus sialome characterized by a poly histidine repeat in the mature aminoterminal region, and a series of gqp/gqg repeats.
the polyhistidine domain may confer a bacteriostatic role for this protein if it chelates zn ions, a bacterial growth factor, and a domain existing in various antimicrobial peptides  <cit> .
comparison of protein sequence identities between c. tarsalis and c. quinquefasciatus gene products
over  <dig> deduced protein sequences coding for putative housekeeping  products are presented in additional file  <dig>  table s <dig>  these proteins allow comparison of the evolutionary rate of the s proteins compared to the h proteins, using the c. quinquefasciatus proteome as a reference set, as done before for comparing an.
stephensi salivary proteins to those of an.
gambiae  <cit> .
for this comparison, we used only protein sequences from c. tarsalis that had at least  <dig> aa of alignment by the blastp tool to a c. quinquefasciatus protein, and excluded from this set possible alleles or closely related gene duplications by removing the smaller sequence that had 80% or more similarity to another one within the set.
we thus compared  <dig> putative secreted c. tarsalis proteins with their c. quinquefasciatus proteome, obtaining an average of  <dig> % protein identity, while  <dig> putative housekeeping proteins from c. tarsalis were  <dig> % identical to c. quinquefasciatus predicted proteins .
this significant difference further supports the concept that the evolution of mosquito salivary secreted proteins occur at a faster pace than housekeeping proteins, as indicated before for anopheline proteins  <cit> .
it has been suggested that codon volatility  could be a measure of selection for fast evolution of proteins, as could occur in pathogens in a constant avoidance of antibody recognition  <cit> .
although this intuitive idea has created strong opposition  <cit> , it is supported by published models  <cit> .
we accordingly decided to measure the average codon volatility for the  <dig> sequences coding for housekeeping and  <dig> sequences coding for putative secreted proteins shown in additional file  <dig>  table s <dig>  the average codon volatility for the h class genes was  <dig>  +  <dig>  while the s class had an average volatility of  <dig>  +  <dig>  , a highly significant result .
whatever the discussion regarding the value of this index, it indicates that a single point mutation on an s class gene has a significantly higher chance of producing a non-synonymous amino acid substitution than in an h class gene.
CONCLUSIONS
from a conservative perspective, the sialotranscriptome of c. tarsalis confirms the presence of ubiquitous salivary mosquito protein families, such as the d <dig>   <dig> kda, hyp <dig>   <dig> kda antigen/aegyptin,  <dig>  kda, mucins, cysteine rich peptide, antigen  <dig>  amylase/maltase, and the immunity related proteins lysozyme, pathogen recognition molecules, and serine proteases possibly associated with immunity or with host matrix/fibrin degradation.
adenosine deaminase has also been found, but not the 5'nucleotidase/apyrase that has been abundantly found in mammalian feeding mosquitoes, but possibly absent or reduced in the bird feeding c. tarsalis.
from another stand point, the c. tarsalis sialotranscriptome confirmed the presence of proteins so far known exclusively in culicine mosquitoes, such as endonucleases and the  <dig> kda,  <dig>  kda,  <dig>  kda and  <dig>  kda families, all without a known function.
additionally, further confirmation of unique proteins of the culex genus were found for members of the highly expanded cwrp,  <dig>  kda and the gpq repeat families.
it has been previously indicated that many unique mosquito protein families appear to be related to bacterial proteins,  suggesting their acquisition by horizontal transfer.
interestingly, the genes coding for such proteins, now available for an.
gambiae, c. quinquefasciatus and ae.
aegypti are often single exonic, as is the case of bacterial genes.
in the case of c. tarsalis, this appears to be the case with the  <dig> kda  and the cwrp  families.
however, other possibilities for single exonic genes are their acquisition as gene duplication resulting from retrotransposition of an mrna precursor deriving from an endogenous, multi exonic gene.
at any rate, the frequency of the relatively unusual processes of horizontal transfer and/or retrotransposition in the acquisition of new genes associated with blood feeding appear to occur at a high rate in the formation of mosquito sialomes.
since our transcriptome was based on ~ <dig>  ests from a non-normalized library, the question arises as to what extent novel putative secreted proteins can be discovered with a more extensive sequencing, or the use of a normalized library.
as indicated in a recent review  <cit> , it has been our experience that, perhaps due to the relatively low complexity of the salivary gland proteome , 1000- <dig> sequenced clones are enough to display the majority of the sialome.
indeed, ~ <dig>  clones sequenced from an ae.
aegypti salivary gland cdna library  <cit>  discovered virtually all those found in another effort that obtained ~ <dig>  sequences  <cit> .
a similar situation was encountered with the an.
gambiae sialome, where ~ <dig>  sequences identified basically all those in a large sequencing effort, also ~ <dig>  sequences done by the pasteur institute  <cit> .
sequencing of normalized libraries, or more intensive sequencing of existing libraries, perhaps with newer generation of pyrosequencing methodology, may increase the number of salivary secreted peptides, but possibly to no more than 10% of the number found with random sequencing of ~ <dig>  clones.
it should be indicated, however, that these additional low abundance transcripts may account for pharmacologically potent peptides.
finally, the fast divergence of salivary proteins allows the possibility of using c. tarsalis proteins as specific markers of vector exposure, as is being attempted for ticks  <cit> , mosquitoes  <cit>  and sand flies  <cit> .

 most adult female mosquitoes are hematophagous, in addition to taking sugar meals. saliva helps blood feeding by interfering with host reactions that could disrupt the blood flow, assists sugar meals with glycosidases, and contain antimicrobials that may control microbial growth in their meals. with the advent of tissue transcriptomics, we can postulate that these functions are mediated by large numbers of polypeptides, many of which are expressed solely in the adult female salivary glands. unique protein families have been found in anopheles, aedes or culex mosquitoes, as well as a group of common proteins or enzymes. functional characterization of these proteins uncovers scavengers of biogenic amines or leukotrienes, inhibitors of blood clotting, bradykinin formation, platelet aggregation and vasodilators. other molecularly uncharacterized activities include inhibitors of mast cell tnf production and inhibition of t cell activation. it is apparent that the complexity of the salivary components affecting host hemostasis and inflammation mirrors the complexity of host hemostasis and inflammation itself, which must be disarmed for successful blood feeding. indeed mosquitoes lacking salivation by salivary duct ablation feed less and take more dangerous time of exposure to their hosts. perhaps due to the potent pharmacological activities of saliva, or the immune reactions to it, mosquito saliva also plays a role in pathogen transmission, including in arboviral transmission. accordingly, determination of the salivary composition of vector mosquitoes not only discovers new potentially pharmacologically active molecules, but also can also help generating vaccine targets for disruption of arboviral transmission. these proteins may also be of epidemiological significance as selective human or animal markers of vector exposure. we have previously described the sialotranscriptome of c. quinquefasciatus, where several unique protein families were discovered, many of which are products of gene duplications. indications of horizontal gene transfer from bacteria to mosquitoes were also pointed out as participating in the generation of mosquito sialomes. to this date, no other sialotranscriptome from a member of the culex genus has been described. we currently portray the sialotranscriptome of culex tarsalis, a north american bird feeding mosquito, and, like c. pipiens, a good vector of west nile virus.  characteristics of the assembled salivary est set a total of cdna clones were used to assemble a database that yielded clusters of related sequences, of which contained only one expressed sequence tag. the clusters were compared, using the programs blastx, blastn, or rpsblast, to the nonredundant protein database of the national center of biological information, to a gene ontology database to the conserved domains database of the ncbi, and to a custom prepared subset of the ncbi nucleotide database containing either mitochondrial or rrna sequences. three categories of expressed genes derived from the manual annotation of the contigs. the putatively secreted category contained 44% of the sequences, the housekeeping category had % and % of the ests could not be classified and belong to the unknown class. the transcripts of the unknown class could represent novel proteins or derive from the less conserved 3' or 5' untranslated regions of genes, as was indicated for the sialotranscriptome of an. gambiae. this distribution is typical of previous mosquito sialotranscriptomes. housekeeping genes the ests attributed to h genes expressed in the salivary glands of c. tarsalis were further characterized into subgroups according to function. transcripts associated with the protein synthesis machinery represented 44% of all transcripts associated with a housekeeping function, an expected result for the secretory nature of the organ. nineteen percent of the transcripts were classified as either "unknown conserved" or "conserved secreted proteins". these represent highly conserved proteins of unknown function, presumably associated with cellular function, but still uncharacterized. these sets may help functional identification of the "conserved hypothetical" proteins as previously reviewed in. transporters accounted for 7% and included several members of the vacuolar atpase machinery, important for mosquito salivary secretion. energy metabolism and protein modification related enzymes accounted for nearly 5% each of the transcripts of the h class. possibly secreted class of expressed genes a total of ests represent putative c. tarsalis salivary components. these include ubiquitous known gene families, e.g, families seen in other organs and organisms, as well as protein families unique to hematophagous diptera, or to culicines, or to culex genus alone. table also indicates our degree of knowledge, or ignorance, regarding these protein families, of which we have no hint for their function, more than half of the classes of proteins represented in table many of these putatively secreted protein families of unknown function are multigenic, such as the cwrp family found previously only in c. quinquefasciatus, which contributes to 58% of the ests associated with the s class in c. tarsalis. the d7/obp-like families contribute to 9% of all transcripts associated with secreted products, which is less than previously seen in anopheles or aedes sialotranscriptomes, probably due to the dilution effect following abundant expression of the cwrp family. * y for yes and n for no - yes is probably an overestimation as proteins can be multi functional the salivary secretome of culex tarsalis from the sequenced cdnas, a total of novel c. tarsalis protein sequences were derived, of which code for putative secreted products. because cdna sequences coding for many of these proteins were found as singletons, this secretome is to be considered preliminary and incomplete, but many parallels with previous mosquito sialotranscriptomes can be drawn, as follows: proteins with presumed or experimentally validated function the d7/odorant binding protein-like family the d proteins constitute a unique multi gene family found expressed in the salivary glands of blood sucking diptera, belonging to the superfamily of odorant binding proteins. long and short versions exist, the long versions containing two and the short versions containing one obp domain. three genes codes for long d proteins and genes code for short d proteins in an. gambiae, while in ae. aegypti genes code for long and genes code for short d proteins. short versions of anopheline d mosquitoes were shown to bind biogenic amines such as serotonin and histamine, thus counteracting hemostasis and inflammation. more recently, the amino terminal obp domain of a d long form of ae. aegypti was shown to bind peptidic leukotrienes with high affinity. the crystal structures of a short d protein from an. gambiae and a long d protein from ae. aegypti revealed that the d obp domains have alpha helices, more than the canonical obp family. in addition to these inflammatory agonist binding functions, a short d protein from an. stephensi, named hamadarin, was shown to inhibit bradykinin formation by inhibiting the fxii/kallikrein pathway. eight members of the d protein family were found in the c. tarsalis sialotranscriptome, in addition to two canonical obp proteins. alignment of these proteins with their close blast matches found in the nr protein database, and construction of a bootstrapped phylogram shows several clades with strong bootstrap support, and allows the following inferences: the eight c. tarsalis proteins are the probable results of genes. ctar- ctar- and ctar- are possible alleles of a single gene related to other long c. quinquefasciatus d proteins shown in clade i, while ctar- ctar- and ctar- are probably alleles of a second c. tarsalis gene coding for yet another long d protein, and are associated within clade ii. ctar- results from a gene coding for a third long d protein grouping in clade v with other mosquito long d proteins, including, anopheles and aedes proteins. ctar- is the only c. tarsalis sequence for a short d protein, grouping in clade iv with other culex pipiens d proteins. notice also clades iii and clade vi, which exclusively contain long d anopheline or phlebotomine proteins, respectively. ctar- ctar- ctar- and ctar- possess the signature - -x-c-x-w-x-w-x- -x-c- -x- -c-x-q-x-c-x- found in lipid binding d domain of ae. aegypti. the serotonin binding motif found on short anopheline d proteins as well as the carboxydomain of the long d protein of a. aegypti was not found on any c. tarsalis d protein, indicating this motif may have evolved beyond recognition, or that other proteins in culex may undertake this task. ctar- and ctar- represent canonical obp protein sequences, transcripts of which were found expressed in the salivary glands of c. tarsalis. homologs of these canonical obp proteins were also previously found expressed in c. quinquefasciatus and an. gambiae sialotranscriptomes. these proteins tend to be much more conserved in sequence to mosquito and even drosophila proteins and it is possible that they may play an endogenous or housekeeping function. the kda antigen/ge rich/aegyptin family this protein family is typical of the salivary glands of adult female mosquitoes, and was first identified as a salivary antigen in ae. aegypti, and later found in salivary transcriptomes and proteomes of both culicine and anopheline mosquitoes, where it was named ge rich protein. distantly related proteins also exist in black flies and sand flies. more recently, proteins of this family from aedes and anopheles were shown to prevent platelet aggregation by collagen, indicating conservation of function after the split of the culicidae into the culicines and anophelines, approximately mya. the sialotranscriptome of c. tarsalis allowed the identification of protein sequences from the kda antigen/aegyptin family, all represented by 1- ests found in the library. blast searches against the nr database retrieved closely related sequences. those that were more than 95% identical and of the same species were excluded to avoid inclusion of alleles or very closely related genes. clustal alignment of this subset was used to produce a bootstrapped phylogram that allows the following inferences: the proteins from c. tarsalis most probably derive from genes, ctar- ctar- and ctar- being probable alleles of a gene closely related to a c. quinquefasciatus salivary gene shown in clade i. clade i includes aedes aegypti salivary proteins that are short versions of the canonical kda protein, as are all culex proteins in this clade. a second c. tarsalis gene possibly codes for ctar- ctar- and ctar- and a third gene for ctar- and ctar- these two gene products groups with strong bootstrap support to their c. quinquefasciatus homologs, and these in turn group relatively weakly in clade ii as shown in fig fig also shows the uniquely aedes containing proteins in clade iii, represented by pairs of homologs between ae. aegypti and ae. albopictus, and the uniquely anopheline clade iv, with one sequence each of different mosquito species. this phylogram suggests the single gene status of this family in the anopheles genus, and the multi gene character in culicines. at least genes for canonical kda ag/aegyptin exists in the aedes genus, and at least one more gene for the shorter protein shown in clade i. culex also have the shorter version gene, plus at least genes of the canonical type. mucins and peritrofins serine and threonine rich proteins are a common finding in sialotranscriptomes. these proteins are generally modified post-translationally and their mature forms have n-acetyl galactosamine residues, typical of mucins. they probably have a function to lubricate the food canals and may also have antimicrobial function. ctar- and ctar- encode related protein sequences that might reflect splice variants of the same gene. these truncated protein sequences have over predicted glycosylation sites, and are similar to a previously described salivary mucin of c. quinquefasciatus. ctar- is related to aedes and culex mucins that, in the case of the aedes protein, were associated with induction by viral infection suggesting an immune function for this protein with putative glycosylation sites. ctar- encodes a 5' truncated protein sequence similar to salivary mucins of c. quinquefasciatus and aedes albopictus. enzymes additional file table s indicate the presence of transcripts coding for enzymes possibly associated with the sugar meal and blood meal. related to the sugar meal are several transcripts coding for glycosidases similar to proteins annotated as maltase and amylase. additional file table s provides one full length sequence of a salivary alpha-glucosidase similar to other culicine salivary maltases, plus other truncated sequences coding for different sugar hydrolases. regarding blood meal-related enzymes, est's coding for fragments of adenosine deaminase were found, as well as for an endonuclease. in c. quinquefasciatus a salivary endonuclease has been previously characterized and associated with helping forming the feeding hematoma. however, the c. tarsalis enzyme is most closely related to sand flies salivary endonucleases and only distantly related to the salivary endonuclease of c. quinquefasciatus, although it is closely related to a non-salivary c. quinquefasciatus enzyme, indicating this c. tarsalis enzyme may be playing a housekeeping function. remarkably absent from the c. tarsalis sialotranscriptome are ests coding for members of the 5' nucleotidase, which functions as the salivary apyrase of mosquitoes. apyrase hydrolyzes atp and adp to amp and orthophosphate, destroying these important agonists of inflammation and platelet aggregation. it has been previously noticed that c. quinquefasciatus has little salivary apyrase activity when compared to aedes and anopheles mosquitoes, and this observation was postulated to be a consequence of the lack of platelets in birds, the most common host of culex, while aedes and anopheles are mostly mammalian feeders. this may explain the absence of 5'-nucleotidase/apyrase coding transcripts in c. tarsalis, although an increased sequencing effort could produce 5'-nucleotidases that might be secreted. as an example of the abundance of this type of transcript in mosquito sialotranscriptomes, only one transcript out of clones coded for apyrase in the c. quinquefasciatus sialotranscriptome, while out of were found in a similarly made library from an. gambiae, and out of ests were found in the ae. aegypti sialotranscriptome. from these observed frequencies, the expected number of est's for c. tarsalis and c. quinquefasciatus are and respectively, producing a chi square of and a p <, indicating the low expression of apyrase transcripts in culex. three ests coding for a phosphatase produce the predicted truncated sequence encoded by ctar- which is similar to a salivary alkaline phosphatase previously described in ae. aegypti. the function of this enzyme in feeding, if any, is unknown. transcripts coding for at least different serine proteases were found in c. tarsalis sialotranscriptome. these enzymes may function in immunity as prophenoloxidase activators, or in digesting skin matrix components, such as in an elastase function, or hydrolysing host blood clotting enzymes such as fibrinogen/fibrin, or activating plasminogen. immunity related products antimicrobial peptides, lysozyme, and pathogen pattern recognition polypeptides are commonly found in the sialotranscriptome of blood sucking arthropods. additional file table s shows the full length sequence of c. tarsalis salivary lysozyme, which is 91% identical to the c. quinquefasciatus homolog, and 75% identical to the salivary homolog of ae. albopictus. truncated orf's of a c-type lectin and a gram-negative binding protein were also found. they both match previously described salivary proteins of aedes and culex. secreted proteins with unknown function promiscuous antigen family this is a ubiquitous protein family found in animals and plants, and in all sialotranscriptomes of blood sucking diptera analyzed so far. the function of these proteins in mosquito saliva is not known, but in blood sucking brachycera two proteins of this family have been functionally characterized. remarkably, in a tabanid fly, a member of the ag- family acquired a typical rgd domain surrounded by cys residues and acts as a main platelet aggregation, and in the stable fly a salivary ag- protein binds immunoglobulins and may function as an inhibitor of the classical complement pathway. we present evidence, in the form of truncated transcripts, for the expression of at least two members of the family in c. tarsalis salivary glands; ctar- assembled from ests, matches with 94% identity the salivary secreted antigen- precursor ag5- from culex quinquefasciatus while ctar- matches with 91% identity another c. quinquefasciatus protein of the same family. insect specific cys-rich polypeptide family two transcripts were found in the c. tarsalis sialotranscriptome coding for basic peptides of mature mw of kda containing cys residues. this peptide family was previously found in the salivary transcriptome of c. quinquefasciatus, but close relatives of the same size exist in drosophila, bombyx, tribolium and apis. the ubiquity of this protein family in insects, its size and pi suggests an antimicrobial role. hematophagous diptera specific kda family the first kda family member was characterized in the sialotranscriptome of ae. aegypti, and later found in c. quinquefasciatus and in ae. albopictus. although not present in the sialotranscriptomes of members of the anopheline cellia subgenus an. funestus, an. stephensi and an. gambiae, it was recently found in an. darlingi, a member of the nyssorhynchus subgenus, characterizing this family as uniquely culicidae. additional file table s provides evidence of a member of this protein family in c. tarsalis, encoded by ctar- producing a predicted protein of mature mw of kda, being 69% identical to its c. quinquefasciatus homolog. psiblast of ctar- against the nr protein database retrieves on its first blastp cycle only mosquito salivary proteins, as expected. on the second cycle, it retrieves with lower significance mostly bacterial proteins, but also salivary proteins from simulium vittatum and culicoides sonorensis. further iterations of psiblast retrieves with high significance bacterial proteins of the methyl-accepting chemotaxis receptor family, which may suggest this bacterial family to be originated from horizontal gene transfer to an ancestral blood feeding nematocera. the gene structure of this protein in c. quinquefasciatus and ae. aegypti is similar, containing exons with a short intron of ca. nt. inclusion of simulium and culicoides sequences in this family indicates an ancient origin before the nematocera split. mosquito specific kda protein ctar- codes for a protein related to mosquito proteins so far only found expressed in mosquito adult salivary glands. psiblast of ctar- against the nr database does not retrieve any additional protein with an e value better than. mosquito specific hyp family this member of this mosquito salivary protein family was identified in an. stephensi, and later found also an. gambiae, but not previously in other mosquito transcriptomes, including c. quinquefasciatus and ae. aegypti. ctar- represented by a single est in our database, codes for a member of this family as indicated by blastp comparisons to the nr protein database, where it retrieves only culex and anopheles proteins. psiblast of ctar- against the nr database does not increase finding matches beyond additional mosquito proteins, including proteins deducted from the uncovering of the c. quinquefasciatus genome with and without signal peptide, and additional proteins from an. gambiae lacking a signal peptide indicative of secretion. interestingly, the anopheles gambiae gene found expressed in the salivary glands was shown to reside as a single exon in chromosome arm 2r. the an. gambiae gene was also found to be selectively expressed in the adult female salivary glands suggesting a role in blood feeding. culicine specific kda family members of the kda protein were previously discovered in ae. aegypti, ae. albopictus and c. quinquefasciatus sialotranscriptomes. additional related proteins were also deducted from the genomes of ae. aegypti and c. quinquefasciatus, but not an. gambiae. the sialotranscriptome of c. tarsalis provides additional evidence for this multi copy family exclusive of culicines. clustal alignment of these protein sequences shows many conserved amino acids and a conserved cys framework on the second half of the protein. the bootstrapped phylogram shows robust clades. clade i contains a sub clade of ae. aegypti and ae. albopictus proteins, and a second sub clade of a c. quinquefasciatus and ctar- clade ii has c. tarsalis sequences, possibly alleles of a single gene or closely related genes, and one c. quinquefasciatus sequence. clade iii shows a possible gene expansion in c. quinquefasciatus containing gene products, clustering with one ae. aegypti sequence. overall the phylogram indicates that aedes aegypti has at least genes coding for this protein family, while c. quinquefasciatus has at least of which possibly arrived by further gene duplications within clade iii. culex tarsalis and ae. albopictus, have at least two genes expressing these proteins in their adult female salivary glands. psiblast of members of this protein family against the nr protein database converges after iterations, retrieving solely aedes and culex proteins. the function of this protein family is unknown, but transcripts for this family were found enriched in the salivary glands of adult female ae. aegypti suggesting a blood feeding role. culicine specific kda family ctar- codes for a mature peptide of predicted mw = kda and pi =. it is similar to previously described salivary peptides found in sialotranscriptomes of ae. aegypti and ae. albopictus, but not culex. psiblast using an inclusion e value of against the nr protein database additionally identifies a salivary peptide previously found in c. quinquefasciatus, annotated as kda salivary peptide. inclusion of this peptide with the previous to generate the reverse position search model additionally identifies other salivary peptides from c. quinquefasciatus, with convergence of the search and identification of the non redundant peptides with an e value of e- or smaller. alignment of these peptides with ctar- shows that the aedes sequences are more compact than those from culex, containing a conserved framework of cysteines, while those of culex have additional cys residues. only other residues are identical in all sequences, indicating a high evolutionary rate of this peptide family, if they indeed share a common ancestor. the bootstrapped phylogram shows robust clades, clade i containing only aedes sequences in two subclades, each containing a pair of ae. aegypti or ae. albopictus sequences. clade ii has sequences from c. quinquefasciatus. ctar- does not cluster with either of the clades, indicating again the possible fast evolutionary rate of this peptide family. culicine specific kda peptide family ctar- ctar- and ctar- encode peptides with similarities to c. quinquefasciatus kda salivary peptide, which has similarity solely to another salivary peptide previously found in ae. albopictus. culex genus specific cysteine and tryptophan rich protein family this protein family was previously discovered in the c. quinquefasciatus sialotranscriptome, where proteins were deduced from the est sequences. these have a signal peptide indicative of secretion, conserved cysteins, conserved tryptophans and one conserved asparagine and leucine residues. inspection of the c. quinquefasciatus genome indicates a total of genes coding for this protein family, most members being single exonic. expression of these genes in c. quinquefasciatus salivary glands was confirmed by edman degradation of electrophoretically separated salivary homogenate. the c. tarsalis sialotranscriptome led to identification of full length members of this protein family, plus other fragments. alignment of the members of this family shows a conserved cysteine framework in most members, two conserved tryptophan residues, as well at the conserved asparagine and leucine residues. this protein family is by far the most transcribed in c. tarsalis adult female salivary glands, accounting for 25% of the transcripts that were sequenced. interestingly, the most transcribed member in c. tarsalis, ctar- is a truncated member of the family by relatively recent acquisition of a stop codon just after the second conserve cysteine. this abundant expression suggests that this family may be responsible for host serotonin or histamine antagonism, as is performed by the abundantly expressed d proteins of aedes or culex. this possibility is further supported by the relatively low expression of the d protein family in c. tarsalis and their absent serotonin binding motifs which are found both in aedes and anopheles d proteins. psiblast of members of the cwrp protein family retrieves, interestingly, sugar binding proteins including bacterial glycosidases, lectins such as ricin and hemolysins, the crystal structure of which have a trefoil structure. the single exon structure of this protein family in c. quinquefasciatus suggests either acquition of this gene family by horizontal transfer, and/or gene duplication occurring by retroposition of an rna template. culex genus specific peptide family ctar- and ctar- encode truncated peptides matching the previously described putative kda salivary peptide of c. quinquefasciatus, which has no other matches to known proteins. culex genus specific gqp repeat protein family ctar- ctar- and ctar- produce matches to a low complexity protein previously described in the c. quinquefasciatus sialome characterized by a poly histidine repeat in the mature aminoterminal region, and a series of gqp/gqg repeats. the polyhistidine domain may confer a bacteriostatic role for this protein if it chelates zn ions, a bacterial growth factor, and a domain existing in various antimicrobial peptides. comparison of protein sequence identities between c. tarsalis and c. quinquefasciatus gene products over deduced protein sequences coding for putative housekeeping products are presented in additional file table s these proteins allow comparison of the evolutionary rate of the s proteins compared to the h proteins, using the c. quinquefasciatus proteome as a reference set, as done before for comparing an. stephensi salivary proteins to those of an. gambiae. for this comparison, we used only protein sequences from c. tarsalis that had at least aa of alignment by the blastp tool to a c. quinquefasciatus protein, and excluded from this set possible alleles or closely related gene duplications by removing the smaller sequence that had 80% or more similarity to another one within the set. we thus compared putative secreted c. tarsalis proteins with their c. quinquefasciatus proteome, obtaining an average of % protein identity, while putative housekeeping proteins from c. tarsalis were % identical to c. quinquefasciatus predicted proteins. this significant difference further supports the concept that the evolution of mosquito salivary secreted proteins occur at a faster pace than housekeeping proteins, as indicated before for anopheline proteins. it has been suggested that codon volatility could be a measure of selection for fast evolution of proteins, as could occur in pathogens in a constant avoidance of antibody recognition. although this intuitive idea has created strong opposition, it is supported by published models. we accordingly decided to measure the average codon volatility for the sequences coding for housekeeping and sequences coding for putative secreted proteins shown in additional file table s the average codon volatility for the h class genes was + while the s class had an average volatility of +, a highly significant result. whatever the discussion regarding the value of this index, it indicates that a single point mutation on an s class gene has a significantly higher chance of producing a non-synonymous amino acid substitution than in an h class gene.  from a conservative perspective, the sialotranscriptome of c. tarsalis confirms the presence of ubiquitous salivary mosquito protein families, such as the d kda, hyp kda antigen/aegyptin, kda, mucins, cysteine rich peptide, antigen amylase/maltase, and the immunity related proteins lysozyme, pathogen recognition molecules, and serine proteases possibly associated with immunity or with host matrix/fibrin degradation. adenosine deaminase has also been found, but not the 5'nucleotidase/apyrase that has been abundantly found in mammalian feeding mosquitoes, but possibly absent or reduced in the bird feeding c. tarsalis. from another stand point, the c. tarsalis sialotranscriptome confirmed the presence of proteins so far known exclusively in culicine mosquitoes, such as endonucleases and the kda, kda, kda and kda families, all without a known function. additionally, further confirmation of unique proteins of the culex genus were found for members of the highly expanded cwrp, kda and the gpq repeat families. it has been previously indicated that many unique mosquito protein families appear to be related to bacterial proteins, suggesting their acquisition by horizontal transfer. interestingly, the genes coding for such proteins, now available for an. gambiae, c. quinquefasciatus and ae. aegypti are often single exonic, as is the case of bacterial genes. in the case of c. tarsalis, this appears to be the case with the kda and the cwrp families. however, other possibilities for single exonic genes are their acquisition as gene duplication resulting from retrotransposition of an mrna precursor deriving from an endogenous, multi exonic gene. at any rate, the frequency of the relatively unusual processes of horizontal transfer and/or retrotransposition in the acquisition of new genes associated with blood feeding appear to occur at a high rate in the formation of mosquito sialomes. since our transcriptome was based on ~ ests from a non-normalized library, the question arises as to what extent novel putative secreted proteins can be discovered with a more extensive sequencing, or the use of a normalized library. as indicated in a recent review, it has been our experience that, perhaps due to the relatively low complexity of the salivary gland proteome, 1000- sequenced clones are enough to display the majority of the sialome. indeed, ~ clones sequenced from an ae. aegypti salivary gland cdna library discovered virtually all those found in another effort that obtained ~ sequences. a similar situation was encountered with the an. gambiae sialome, where ~ sequences identified basically all those in a large sequencing effort, also ~ sequences done by the pasteur institute. sequencing of normalized libraries, or more intensive sequencing of existing libraries, perhaps with newer generation of pyrosequencing methodology, may increase the number of salivary secreted peptides, but possibly to no more than 10% of the number found with random sequencing of ~ clones. it should be indicated, however, that these additional low abundance transcripts may account for pharmacologically potent peptides. finally, the fast divergence of salivary proteins allows the possibility of using c. tarsalis proteins as specific markers of vector exposure, as is being attempted for ticks, mosquitoes and sand flies. 
the average amino acid identity among salivary proteins is %, while that for housekeeping proteins is %, and the codon volatility of secreted proteins is significantly higher than those of housekeeping proteins. comparison of the c. tarsalis sialotranscriptome with that of c. quinquefasciatus reveals accelerated evolution of salivary proteins as compared to housekeeping proteins. mosquito saliva also potentiates the transmission of vector borne pathogens, including arboviruses. culex tarsalis is a bird feeding mosquito vector of west nile virus closely related to c. quinquefasciatus, a mosquito relatively recently adapted to feed on humans, and the only mosquito of the genus culex to have its sialotranscriptome so far described. a total of clones randomly selected from an adult female c. tarsalis salivary glands cdna library were sequenced and used to assemble a database that yielded clusters of related sequences, of which were singletons. several protein families previously found exclusive of mosquitoes, including only in the aedes genus have been identified in c. tarsalis. interestingly, a protein family so far unique to c. quinquefasciatus, with genes, is also found in c. tarsalis, indicating it was not a specific c. quinquefasciatus acquisition in its evolution to optimize mammal blood feeding. saliva of adult female mosquitoes help sugar and blood feeding by providing enzymes and polypeptides that help sugar digestion, control microbial growth and counteract their vertebrate host hemostasis and inflammation.
 saliva of adult female mosquitoes help sugar and blood feeding by providing enzymes and polypeptides that help sugar digestion, control microbial growth and counteract their vertebrate host hemostasis and inflammation. mosquito saliva also potentiates the transmission of vector borne pathogens, including arboviruses. culex tarsalis is a bird feeding mosquito vector of west nile virus closely related to c. quinquefasciatus, a mosquito relatively recently adapted to feed on humans, and the only mosquito of the genus culex to have its sialotranscriptome so far described.  a total of clones randomly selected from an adult female c. tarsalis salivary glands cdna library were sequenced and used to assemble a database that yielded clusters of related sequences, of which were singletons. primer extension experiments were performed in selected clones to further extend sequence coverage, allowing for the identification of protein sequences, of which code for putative secreted proteins.  comparison of the c. tarsalis sialotranscriptome with that of c. quinquefasciatus reveals accelerated evolution of salivary proteins as compared to housekeeping proteins. the average amino acid identity among salivary proteins is %, while that for housekeeping proteins is %, and the codon volatility of secreted proteins is significantly higher than those of housekeeping proteins. several protein families previously found exclusive of mosquitoes, including only in the aedes genus have been identified in c. tarsalis. interestingly, a protein family so far unique to c. quinquefasciatus, with genes, is also found in c. tarsalis, indicating it was not a specific c. quinquefasciatus acquisition in its evolution to optimize mammal blood feeding. 