Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Transcriptomic and functional analysis of the Anopheles gambiae salivary gland in relation to blood feeding

BMC Genomics volume 11, Article number: 566 (2010) Cite this article

Abstract

Background

The Anopheles gambiae salivary glands play a major role in malaria transmission and express a variety of bioactive components that facilitate blood-feeding by preventing platelet aggregation, blood clotting, vasodilatation, and inflammatory and other reactions at the probing site on the vertebrate host.

Results

We have performed a global transcriptome analysis of the A. gambiae salivary gland response to blood-feeding, to identify candidate genes that are involved in hematophagy. A total of 4,978 genes were found to be transcribed in this tissue. A comparison of salivary gland transcriptomes prior to and after blood-feeding identified 52 and 41 transcripts that were significantly up-regulated and down-regulated, respectively. Ten genes were further selected to assess their role in the blood-feeding process using RNAi-mediated gene silencing methodology. Depletion of the salivary gland genes encoding D7L2, anophelin, peroxidase, the SG2 precursor, and a 5'nucleotidase gene significantly increased probing time of A. gambiae mosquitoes and thereby their capacity to blood-feed.

Conclusions

The salivary gland transcriptome comprises approximately 38% of the total mosquito transcriptome and a small proportion of it is dynamically changing already at two hours in response to blood feeding. A better understanding of the salivary gland transcriptome and its function can contribute to the development of pathogen transmission control strategies and the identification of medically relevant bioactive compounds.

Background

Adult mosquitoes feed on sugar to obtain energy for flight and other activities, while anautogenous females need a blood meal to develop eggs. Salivary glands and the saliva of insect disease vectors have attracted considerable attention because of their role in pathogen transmission and their production of pharmacologically active factors [1–4].

It is during the blood-feeding process that the Plasmodium parasite is taken up from an infected A. gambiae host. Once inside the mosquito, Plasmodium undergoes several developmental transitions and eventually becomes a sporozoite, which invades the salivary glands. This invasion represents a critical step in the transmission of the parasite to the vertebrate host. Completion of the infection depends on the injection of sporozoites, through the saliva, into the host's skin and leaving the inoculation site rapidly to enter and invade the liver for further development [5, 6].

Unlike male salivary glands, female mosquito salivary glands possess anti-hemostatic, vasodilatory and immune-modulatory factors to facilitate the acquisition of blood, while salivary glands of both sexes have activity related to the digestion of the sugar meal as well as antimicrobials to prevent microbial growth [7, 8]. Mosquitoes have been shown to require longer probing times during blood-feeding when an apyrase gene (an enzyme that counteracts hemostasis) has been silenced, or if they are deprived of salivation by removal of the salivary duct through which the saliva is transported to the probing site [9, 10]. Silencing of another A. gambiae salivary gland gene, SG6 (a small protein with unknown function), results in increased probing time and reduced blood-feeding ability [11]. Previous studies have shown that several Anopheles salivary gland proteins are reduced after blood-feeding, suggesting that these major polypeptides may have been introduced into the vertebrate hosts during the blood meal [12, 13]. There is evidence that the pharmacological activity of arthropod saliva affects pathogen transmission and the local inflammatory response of the host. For example, the salivary gland lysate from the sand fly Lutzomyia longipalpis facilitates the infection of mice by the protozoan parasite Leishmania major[14, 15]. It has also been shown that the L. longipalpis salivary gland lysate inhibits neutrophil migration and the Th1 immune inflammatory response. These findings suggest that the compounds responsible for such activities could be used for the development of novel anti-inflammatory drugs [4].

While earlier sialo-transcriptomic studies have identified a variety of salivary gland genes, [1, 16–21], we present the first global microarray transcriptome analysis of the A. gambiae salivary gland under conditions related to feeding. Earlier studies have identified some 3,000 adult female A. gambiae salivary gland-transcribed sequence tags and 4,719 genes were found to be transcribed in the larval gland [22], of which 747 were specific for this tissue. Here we report 4,978 adult female A. gambiae salivary gland transcripts, as defined by oligonucleotide microarray gene transcription analysis. We show that 52 and 41 salivary gland-expressed transcripts were up-regulated and down-regulated, respectively, at 2 hours after blood-feeding when compared to salivary glands of unfed mosquitoes. We have also used an RNAi-mediated gene silencing approach to assess the potential involvement of 10 selected salivary gland genes in regulating mosquito blood-feeding capacity. Silencing of several salivary gland transcripts; D7L2, anophelin, peroxidase, 5'nucleotidase and SG2 precursor, produced a significantly lowered blood-feeding phenotype and increased probing time, confirming that these genes may be playing an important role in blood-feeding. The updated list of the A. gambiae salivary gland transcriptome, together with the results of our comprehensive functional analyses, can facilitate the discovery of novel pharmacologically active compounds and provide tools for the development of malaria control strategies.

Results and Discussion

The A. gambiae salivary gland transcriptome

In order to characterize the A. gambiae female salivary gland transcriptome, we employed a microarray-based genome transcription approach to compare the transcript abundance in the salivary glands at 2 hours after blood feeding to salivary glands of unfed female mosquitoes. The relative transcript abundance of salivary gland expressed genes was analyzed by sorting the blood-fed versus unfed salivary gland transcripts into three different categories based on their corresponding spot fluorescence intensities (Figure 1A,B,C; see Additional file 1). Of the total predicted transcriptome (13,133 genes of A. gambiae, according to version 48, ENSEMBL); 4,978 genes (38%) genes were found to be transcribed in the salivary gland of 4-day-old A. gambiae female mosquitoes. Of these, 267 transcripts (2.1%) were considered as highly abundant (with fluorescence values ranging from 5,000 to the maximum threshold of 65,000; Figure 1A and Sheet 1 of Additional file 1); 482 transcripts (3.6%) were considered as moderately present (fluorescence values ranging from 1,000 to 4,999; Figure 1B and Sheet 2 of Additional file 1), and 4,229 transcripts (32.2%) were considered as being present at low levels based on spot intensity (fluorescence values ranging from 100 to 999; Figure 1C and Sheet 3 of Additional file 1).

Figure 1
figure 1

The adult female A. gambiae salivary gland transcriptome and its regulation upon blood-feeding. Genes that correspond to replica spots with a mean salivary gland RNA probe fluorescence hybridization intensity value >100 were grouped into the following three categories: (A) highly abundant transcripts (n = 267), with fluorescence values from 5,000 to 65,000; (B) moderately present transcripts (n = 482), with fluorescence values from 1,000 to 4,999; and (C) weakly present transcripts (n = 4,229), with fluorescence values from 100 to 999. The three categories were further sub-grouped into 11 functional groups according to their predicted functions, as shown in the pie charts. (D) The bar graph shows the number of genes that were regulated upon blood-feeding in different functional groups; those that were up-regulated (52 genes) are above zero (indicated by a red arrow), and those that were down-regulated (41 genes) are below zero (indicated by a green arrow). The color code for each group is shown. Those genes belonging to the unknown functional group are not shown in any of the diagrams.

Full size image

We analyzed the predicted functional class composition of each category (discussed below) to gain insight into the major processes that take place in the female mosquito's salivary glands upon blood-feeding.

Cytoskeletal and structural genes

Several cytoskeletal and structural gene transcripts (9 highly, 25 moderately, and 133 weakly transcribed) were identified in the salivary glands (Figure 1A, B, C and see Additional file 1). One of the most abundant transcripts encoded a mucin, a secretory protein that has been shown to be a major protein in the salivary gland lysate of several mosquito species [1, 16, 19]. This protein may function as a lubricant of the mosquito mouthparts and has been proposed to bind to the chitinous linings of the salivary ducts and mouthparts [1, 23]. Three mucin encoding transcripts have been identified in the A. gambiae larval salivary gland [22], suggesting the importance of mucins at multiple developmental stages. We have further investigated the function of one of the mucin gene (AGAP001192-RA) by assaying changes in the feeding and probing behavior of mosquitoes after its depletion through RNAi (discussed below). In addition, Saglin (AGAP000610-RA) was shown to be an abundant transcript in the salivary gland. Saglin encodes for a secretory protein (100 kDa) that has been shown to bind to the Plasmodium berghei sporozoite protein thrombospondin related anonymous protein (TRAP), and play a crucial role in salivary gland invasion by sporozoites, [24, 25]. The presence of transcripts such as actin, chitin binding, cuticlar protein, myosin, laminin, tropomysin, tubulin, annexin, dynein light chain, connectin, dynactin, cadherin, and myosin muscle suggests that these genes play a significant housekeeping role in the A. gambiae salivary gland.

Genes involved in protein and sugar digestion

As many as 175 genes putatively involved in protein and sugar digestion were found to be transcribed in the A. gambiae salivary glands; 26 of these transcripts were highly abundant, 24 were moderately present and 125 were present at low levels (Additional file1). Many of these transcripts coding for proteases (peptidase C1, peptidase S1, peptidase S6, TRY2, TRY3, TRY4, TRY5, TRY6, TRY7, several serine proteases, carboxypeptidases, and aminopeptidases) and glycosidases (four maltase transcripts AGAP011939-RA, AGAP012401-RA, AGAP002102-RA, and AGAP008965-RA and a chitinase transcript AGAP008061-RA) have a secretory signal peptide sequence.

Salivary gland glycosidases have been shown to be involved in the digestion of sugar, whereas the proteolytic enzymes could be involved in host-specific proteolytic events (clot prevention or the complement and coagulation cascades, or in the digestion of extracellular matrix components) that occur during the blood-feeding process. It is also possible that these proteolytic enzymes, which are secreted in the saliva, are ingested along with the blood, and might be involved in its digestion. Proteins encoded by transcripts that lack signal peptide sequences are most likely involved in housekeeping functions in the salivary gland or may be secreted by an alternative non-classical method [26]. The majority of these genes are also transcribed in the A. gambiae midgut and other tissues as well as during other developmental stages as has been demonstrated in earlier studies [27–29]. The proteolytic digestive enzymes in the midgut are mainly involved in the digestion of blood.

Circadian and chemosensory genes

Two circadian system-related genes, takeout and timeless, were weakly transcribed in the salivary gland. Takeout has been shown to act as a molecular link between circadian rhythm and feeding behavior in Drosophila melanogaster and A. gambiae[30, 31], and timeless is a transcription factor that plays a major role in circadian regulation [32]. The timeless transcript has been identified previously in the A. gambiae female salivary gland [1]; however, we report for the first time the expression of takeout in this tissue.

Several genes with putative chemosensory functions, such as odorant binding proteins (OBPs), odorant receptors, antennal carrier protein, arrestins, gustatory receptors, and pheromone binding proteins, were transcribed in the A. gambiae salivary gland (3 were highly expressed, 4 were moderately expressed, and 19 were weakly expressed; see Additional file 1). OBPs are mainly localized in the mosquito head and within the antennae, where they function as carriers of odorants to the olfactory receptors [33]. The D7 family of proteins are distantly related to the OBP super-family, because of the conserved 6-cysteine and 10-cysteine residues and a characteristic fold structure [34]. However, the OBPs function as odorant carriers [33], whereas D7 proteins have been proposed to inhibit hemostasis by trapping agonists of hemostasis [35], and thereby facilitate blood-feeding by hematophagous vectors (discussed in detail in the next section).

D7 protein family genes

The four highly transcribed members of the A. gambiae D7 family were the D7r-1, 2, 3 precursors (short form) and the D7 long form L2 (AGAP008279-RA) (fluorescence values above 5,000), whereas the D7r-4, 5 precursors (short form) and another D7 long form L1 (AGAP008278-RA) were found to be only moderately/weakly transcribed in the salivary glands (fluorescence values between 100 and 4,999) (see Additional file 1). The D7 gene (coding for salivary biogenic-amine binding protein) was first identified in A. aegypti[36] as a major salivary gland secretory protein. The D7 protein exists in two forms: a long form (~30-35 kDa), which is found exclusively in mosquitoes and sand-flies, and the short forms (~15 kDa), which are found in other insects [37]. Five different D7-related (D7r1, 2, 3, 4, and 5) short forms and three D7 long forms have been identified in the A. gambiae adult salivary gland [1, 38]. In 2006, Calvo et al. proposed that D7-related proteins from blood-sucking insects act as anti-hemostatic factors by trapping agonists of hemostasis. One short D7 protein from A. stephensi, hamadarin (D7r1), has been shown to inhibit the plasma contact system by preventing the activation of kallikrein by Factor XIIa [39]. The D7r1, 2, 3, 4 and D7 long forms have been shown to bind to the biogenic amines serotonin, histamine, and norepinephrine [1, 35]. The crystal structures of one member of the D7 family, D7r4 (short form), and those of the ligand complexes were determined in order to elucidate the mechanisms governing its anti-hemostatic and anti-inflammatory properties [34]. It has been suggested that the ancestral D7 has originated from the proteins of the OBP family, which were adapted to bind small hydrophobic ligands [34, 40]. We have further investigated the function of this family by looking at the changes in feeding and probing behavior of mosquitoes after silencing of the two D7 long forms, D7 L1 and D7 L2 (discussed below).

Immunity- and stress-related genes

Transcripts representing 338 immune genes (20 highly transcribed, 51 moderately transcribed, and 267 weakly transcribed) were found to be expressed in the A. gambiae female salivary glands (see Additional file 1). These genes included: lysozyme, gram negative bacteria binding protein (GNBP), the anti-microbial peptides cecropin and defensin, the NFκB transcription factors Rel1 and Rel2, a c-type lectin (CTL), galectins, several fibrinogen binding proteins (FBN), serpins, thioester containing proteins (TEP), a caspase, SOCS, TOLL, myd-88, a leucine rich receptor (LRR), a scavenger receptor, a prophenoloxidase (PPO), and others. Immune genes such as c-type lectins, prophenoloxidase activating enzymes, galectins, V5 allergens, and lysozyme were identified in the A. gambiae salivary gland [1, 16, 19] and in another study, it was shown that GNBP, defensin, putative gal-lectin, and putative serine protease transcripts were up-regulated in the A. gambiae salivary gland after Plasmodium berghei infection [41]. In a more recent SAGE analysis of the A. gambiae salivary gland transcriptome upon P. berghei infection; 37 immune related genes were identified and four of them (defensin1, GNBP, serpin, cecropin2) were up-regulated during salivary gland invasion by the sporozoites [42].

Several of these immune genes encode secretory proteins, including GNBPs, several FBNs, the short antimicrobial peptides, PGRPs, chitinase, and lectins. Till date, little is known about the role of these immune genes in the A. gambiae female salivary gland, or whether they are directly involved in anti-Plasmodium or anti-microbial defenses. The secretory immune proteins may also be involved in protective immunity against microbes or other pathogens with which the mosquitoes come into direct contact during feeding. In the larval salivary gland, three defensins, three fibrinogen binding proteins, one CLIPB8, and four serine proteases have been detected [22]. Three different chitinase transcripts have been detected in A. gambiae larval salivary glands [22] and have been proposed to protect against pathogenic bacteria and fungi derived from the mosquitoes' diet [43]. Of the three different lysozyme transcripts detected in this study, one (AGAP007386-RA) was also detected in the salivary gland of larvae [22]. Lysozyme might be involved in preventing the growth of microbes in the sugar meal that is ingested with saliva and stored in the mosquito crop [44]. Most of these immune genes are also transcribed in other tissues, where they are involved in limiting bacterial and parasite growth [29, 45–49]. Warr and coworkers [29] have shown that the cardia and anterior part of the midgut transcribe most of the antimicrobial peptides and other anti-Plasmodium factors.

Among the genes involved in oxidoreductive, antioxidant, and stress-related processes in A. gambiae, 270 transcripts were found to be present in the salivary glands (14 highly transcribed, 48 moderately transcribed, and 208 weakly transcribed) (see Additional file 1). This group was composed of several cytochrome P450 enzymes, glutathione s-transferase, cytochrome c oxidase, NADH ubiquinone, DNAJ homolog, oxidoreductase, animal heme peroxidase, thiorededoxin 1, mitochondrial inner membrane protein, and many others. Some of these transcripts have been identified earlier in salivary gland EST libraries or in other studies of different mosquito species [1, 16, 19, 50] and are involved in housekeeping functions of the salivary gland. Several transcripts coding for heat-shock proteins (HSPs) were found to be transcribed in the female salivary gland. HSPs belong to a family of proteins that function as chaperones or are involved in cellular defense against external stress from various sources [51].

Metabolism and other enzyme genes

As many as 392 transcripts (11 highly transcribed, 48 moderately transcribed, and 333 weakly transcribed) of genes putatively involved in metabolism and other enzymatic activities were found to be transcribed in the salivary gland of A. gambiae (see Additional file 1). Among the transcripts encoding putative secretory enzymes, three members of the 5' nucleotidase gene family; an apyrase precursor (AGAP011971-RA), and two 5' nucleotidase precursors (AGAP011026-RA and AGAP003629-RA) were found to be transcribed in the adult salivary gland (see Additional file 1). The apyrase and one of the 5' nucleotidase precursors (AGAP011026-RA) have been shown to be components of saliva [23, 52]. The second 5' nucleotidase precursor gene contains a secretory signal sequence, suggesting potential involvement in feeding. Of these three 5' nucleotidase gene family transcripts, the AGAP011026-RA was highly expressed, and the other two were moderately expressed. Enzymes of this family facilitate the acquisition of blood meals by removing and degrading pharmacologically active ATP/ADP nucleotides from the site of the injury [53].

Silencing of an apyrase gene in the salivary gland of A. gambiae resulted in an increase in the mosquito probing time as well as increase in vitro platelet aggregation [9]. We have further investigated the function of this gene family by assaying changes in the feeding and probing behavior of mosquitoes after silencing the 5' nucleotidase precursor transcript (AGAP011026-RA) (discussed below). Some of the metabolic enzymes with a signal peptide sequence include alpha amylase (AGAP002317-RA) and alpha glucosidase (AGAP000862-RA), which are weakly transcribed in the female salivary glands. The transcript glycohydrolase (AGAP008124-RA), which is weakly transcribed in the A. gambiae salivary gland, also encodes a secretory protein and belongs to a large family with diverse functions reflecting their amino acid sequences: chitosanase, and NAD+ glycohydrolase, among many others [54].

Genes of other functional groups

Several genes involved in housekeeping functions and belonging to the replication-transcription-translation (389 genes in total) and transport functional groups (326 genes in total) were found to be transcribed in the A. gambiae salivary gland (see Additional file 1). In the replication-transcription-translation group, genes including several ribosomal proteins, histones, elongation factors, repressors, activators, proteasome subunits, and enzymes such as helicase, ligase, polymerase, DNA repair, ribonuclease, and many other enzymes were found to be transcribed. In the transport group, transcripts such as vesicle transport protein, solute carrier protein, ion transporter, ion channels, apolipophorin precursor, aquaporin water channel, and many others were found to be transcribed. Several transcripts (including cyclin, ferritin, zeelin 1, anophelin, trio, zinc finger domain, SG2, SG3, SG5, SG6 and SG10, many hypothetical precursors/peptides, and others) were found to be transcribed and were placed in the diverse group (1,277 transcripts in total) because of their assorted roles and various functions. The presence of vitellogenin precursor transcripts is likely to indicate contamination with fat body which is attached to the salivary glands.

A study recently showed that silencing of SG6 transcript results in an increased probing time and reduced feeding ability in A. gambiae female mosquitoes [11]. In this study we have investigated the function of the SG2 transcript (the function of the SG2 protein is as yet unknown), by assaying changes in the feeding and probing behavior of mosquitoes after its depletion (discussed below). The SG1 protein precursor, a secretory protein, was found to be highly transcribed in the salivary gland; SG1 protein family members have molecular masses around 44·kDa and do not exhibit significant similarities [55] to other proteins in the NCBI database, with the exception of the Anopheline distantly related TRIO protein (an abundant salivary gland gene) [3, 19], when analyzed by BLAST analysis. Transcripts of unknown function or with no homology to known protein were assigned to the group of unknown (1,609 transcripts in total). Further studies are required to link these transcripts with specific roles of salivation, feeding, blood digestion, and storage in A. gambiae.

Regulation of A. gambiae salivary gland transcripts upon blood-feeding

We have performed a microarray-based gene transcription analysis of the blood-feeding-regulated genes in the A. gambiae salivary gland. We compared the transcriptome of the A. gambiae salivary glands at 2 h after blood-feeding to salivary glands of unfed mosquitoes to identify salivary gland transcripts that are likely to play a role, either directly or indirectly, in the feeding process. Previous studies have shown that salivary gland proteins are rapidly depleted during feeding and that this depletion is proportional to the duration of feeding. Since blood fed mosquitoes are capable of taking a second blood-meal already at 24 hours after the first feeding, we hypothesized that the replenishment of the salivary gland proteins must initiate soon after feeding, and therefore selected a 2 hour time point after feeding to assay blood feeding-related changes in the salivary gland transcriptome [23, 56]. Genome-wide microarray analysis of gene expression responses to blood feeding has been performed previously but for whole female A. gambiae mosquitoes [57], whereas our study is exclusively for the female salivary glands.

The resulting catalog of transcripts that exhibited differences in abundance prior and after blood feeding in the A. gambiae salivary gland is presented in Additional file 2. In total, 93 transcripts were differentially regulated; 52 transcripts were more abundant and 41 transcripts were less abundant after blood-feeding (Figure 1D and see Additional file 2). Marinotti et al. showed in 2005 that some 4,924 transcripts, derived from whole mosquitoes displayed changes in abundance within 24 hour after blood-meal; 2,388 genes were up-regulated while 2,536 genes were down-regulated [57]. The salivary gland expressed transcripts that were significantly regulated by blood-feeding are presented in Table 1 (excluding the diverse and unknown groups).

Table 1 The female A. gambiae salivary gland expressed genes that are significantly differentially regulated upon blood-feeding.
Full size table

The down-regulation of several female A. gambiae salivary gland transcripts can be attributed to the depletion of transcripts during protein translation or to the degradation of transcripts following blood-feeding. It is possible that their transcription may have been shut down because these transcripts were no longer required until the next blood-meal. Down-regulation of several A. aegypti salivary gland transcripts upon blood-feeding has been shown in a recent publication; the list of transcripts involved coded for secretory proteins, such as odorant binding protein, protease inhibitors, immune genes and other [58].

The secretory protein encoding salivary gland transcripts that displayed a lower abundance after blood-feeding included two OBPs (OBP 10 and OBP 7), two D7 long-form precursors (L1 and L2), two aminopeptidases, a trypsin 6 precursor, a salivary lipase, a 5' nucleotidase precursor, an apyrase, E1 protein, cecropin 3, defensin 1, and a hypothetical 6.2 precursor. In earlier studies, several OBP transcripts were shown to be down-regulated after blood-feeding in the head and antennae of A. gambiae[30] and in the salivary gland of A. aegypti[58]. In Marinotti et al., 2005, transcripts of several serine protease, aminopeptidase, odorant binding protein and transporter genes were found to be down-regulated in whole female mosquitoes 24 hrs after blood feeding. The OBP and odorant receptors are required by the female mosquitoes to sense host olfactory cues. The fact that the OBPs are not required after feeding until the next blood-meal may account for their lower transcript abundance. Both the secretory proteins aminopeptidases and trypsin 6 may be required for the digestion of host proteins or extracellular matrix components [59]. The apyrase and 5' nucleotidase proteins are known to facilitate the acquisition of a blood meal by removing pharmacologically active nucleotides that are important for platelet aggregation at the site of the injury [53]. The D7 proteins inhibit hemostasis, and one short A. stephensi D7 protein, hamadarin, has been shown to inhibit the plasma contact system [39]. The secretory immune proteins like cecropin and defensin may provide protection against microbes or other pathogens that the mosquitoes encounter during feeding. In 2007, Rosinski-Chupin et al. showed that the immune peptide defensin1 and cecropin 2 transcripts were up-regulated in P. berghei infected A. gambiae salivary glands, thereby demonstrating that the salivary epithelium responds to the presence of pathogens [42].

Several transcripts were up-regulated after blood-feeding, and those that encode putative secretory proteins included a trypsin 2 precursor, lysozyme, an OBP, an ATP binding protein, a vitellogenin precursor, and a DNAJ precursor. In contrast, two lysozyme transcripts were found to be down-regulated in whole female A. gambiae mosquitoes at 24 hours after blood-meal [57]. Lysozyme has been proposed to prevent bacterial growth in the sugar and blood meals [1, 16]. The DNAJ family of proteins is involved in protein folding, protein transport, and the cellular response to stress [60]. The exact role of these genes in the salivary gland and their function in the blood-feeding process are still not clear. We have further investigated the function of some of these differentially regulated genes by looking at how they influence the mosquito feeding and probing behavior upon depletion.

Functional analysis of salivary gland gene implication in blood-feeding and probing behavior

To identify putative components of the A. gambiae salivary gland that play a role in blood-feeding, we performed transient RNAi gene-silencing assays on 10 salivary gland transcribed genes and studied the effect of this treatment on mosquito feeding behavior on a vertebrate host (Figure 2A). In our first assay, instead of assaying probing time of individual mosquitoes we determined the blood feeding capacity (4 days after dsRNA injections) as a measure of their feeding status after 20 minutes exposure to the vertebrate host since probing time may not be the only determining factor of feeding capacity, nor an exclusive indicator for a salivary gland function that facilitates feeding. It is conceivable to hypothesize that the amount of ingested blood is dependent on probing time, feeding time and the flow rate of blood which may depend on its viscosity. We specifically compared blood-feeding capacity between gene silenced and control group mosquitoes that had been injected with GFP dsRNA according to previously established methodology [30]. Four of these genes, two D7 long forms (D7L1 and D7L2), one salivary gland peroxidase 5B (AGAP010735-RA), and one 5' nucleotidase (AGAP011026-RA), were selected because of their predicted anti-inflammatory, vasodilatory, and anti-platelet activities, respectively [39, 53, 61] as well as their down-regulation after blood-feeding (see Additional file 2). A salivary gland lipase enzyme (AGAP005822-RA) whose function is not clear and down-regulated after blood-feeding (see Additional file 2) was selected for functional analysis. Five other genes were selected based on their predicted involvement in anti-hemostatic processes and for representing major A. gambiae salivary gland secretory proteins: anophelin/cEF (AGAP008004-RA), which acts as an α-thrombin inhibitor to prevent blood coagulation [62]; a salivary 30 kDa protein (AGAP009974-RA), a member of the aegyptin family, and whose members have been shown to bind to collagen and block platelet adhesion to collagen [63, 64]; an SG2 precursor (AGAP006506-RA) with unknown function; and salivary mucin (AGAP001192-RA) and TRIO (AGAP001374-RA), which are two highly abundant salivary gland transcripts in A. gambiae[1, 19, 62] (see Additional file 1). In Drosophila melanogaster, the TRIO protein has been shown to be involved in the cytoskeletal rearrangements necessary for cell migration [65]; however, its exact function in A. gambiae is not known.

Figure 2
figure 2

Implication of salivary gland genes in blood-feeding capacity. A. The percentage of gene-silenced mosquitoes that fed on a mouse after gene-specific silencing through RNAi compared to the GFP dsRNA injected control mosquitoes. Data were pooled from six replicates of 25 mosquitoes each. Genes displaying significant capacity to alter blood feeding propensity upon silencing are indicated with an asterisk, and the p-values are shown for each gene (from the Mann-Whitney test). The error bars indicate the standard error. B. Determination of RNAi-mediated gene silencing efficiency, by qRT-PCR at 4 days after dsRNA injection. The cDNA templates were normalized using the AgS7 gene specific primers. The graph shows the relative transcript abundance of each gene after knockdown, as compared to its abundance in the GFP dsRNA-treated control mosquitoes (set to 1.0). The error bars indicate the standard error.

Full size image

RNAi-mediated gene silencing assays in adult mosquitoes depend on the direct injection of dsRNAs into the hemolymph and we used the dsRNA injection procedure established by Boisson et al., 2006 (see Materials and Methods) and monitored gene silencing efficiency in the salivary glands 4 days after dsRNA injections (Figure 2B). Silencing efficiency for the selected genes ranged from 19% to 46% (Figure 2B and see Additional file 3). The effect of gene silencing on the mosquito's blood-feeding capacity is shown in Figure 2A, and the statistical analysis of the data is presented in Additional file 4.

Only mosquitoes in which D7L2, anophelin, peroxidase 5B, or SG2 precursor was silenced displayed a significantly altered feeding capacity (by ~15-22%) when compared to the GFP dsRNA-injected control mosquitoes (Figure 2A). The D7 proteins have been shown to bind biogenic amines such as serotonin and histamine, which may antagonize the vasoconstriction, platelet aggregation, and pain induction that occur during feeding [34]. The lower feeding capacity observed after the silencing of D7L2 strongly supports the involvement of its encoded protein (and of other members of this family) in the blood-feeding process. This protein may act as an anti-hemostatic factor by trapping agonists of hemostasis, as suggested by Calvo et al. (2006). The salivary peroxidase 5B is known to have vasodilatory activity. Therefore, we hypothesize that the lower blood-feeding capacity of peroxidase-silenced mosquitoes is due to the resulting increased vasoconstriction and diminished blood-flow at the probing site. The anophelin gene encodes an α-thrombin inhibitor that prevents blood coagulation during the blood-feeding process. Depletion of the anophelin protein is expected to decrease the blood-flow at the site of probing; an effect that may account for the lower feeding propensity of the anophelin-silenced mosquitoes. Silencing of a small secretory peptide SG2 precursor also resulted in a decreased blood-feeding capacity, although its role during blood-feeding is not yet clear. In earlier studies, silencing of a salivary apyrase gene (involved in the inhibition of platelet aggregation in A. gambiae) and a SG6 precursor led to an increased probing time and decreased blood-feeding ability in A. gambiae mosquitoes [9, 11]. However, further studies are needed to characterize the biological properties of this unknown peptide SG2.

The D7L1 and 5' nucleotidase genes are known to have anti-inflammatory and anti-platelet activity, respectively; however, the silencing of these two genes resulted in only an 8-10% decrease in their blood-feeding capacity (when compared to the GFP control). The lack of any significant effects on feeding after silencing of the remaining six genes (D7L1, TRIO, 5' nucleotidase, salivary mucin, salivary lipase, and 30 kDa) may be related to the relatively low silencing efficiency obtained in the salivary gland as well as to the presence of unknown factors that perform redundant functions in the mosquito salivary gland.

In a second set of assays we determined the effect of salivary gland transcript depletion on the mosquitoes probing time (Figure 3), as a measure of time, or interval, between the insertion of mouthparts into the host skin and the initiation of blood ingestion, according to methodology established in previous studies [9, 11]. Of the six genes that were selected for probing time assays, depletion of D7L2, anophelin, peroxidase and SG2 precursor (Figure 2A) resulted in decreased blood feeding capacity and as well as increased probing time (Figure 3). These results strongly suggest a vital role in the blood-feeding process for these genes. Depletion of the 5' nucleotidase resulted in a significantly longer probing time (Figure 3) but not a significantly altered blood-feeding capacity (Figure 2A). This indicates that blood-feeding capacity is not exclusively dependent on probing time.

Figure 3
figure 3

Functional implication of salivary gland genes in the probing behavior of mosquitoes. Probing time spent by gene silenced mosquitoes (depleted of D7-L2, anophelin, peroxidase, trio, SG2 precursor and 5' nucleotidase transcripts), and GFP dsRNA injected control mosquitoes, from the insertion of mouthparts into the mouse skin and the first sign of blood within the midgut. Four days after dsRNA injection, individual mosquitoes were exposed to an anesthetized mouse for 500 seconds. Probing times were only measured for those mosquitoes that initiated ingestion of blood within the 500 seconds interval. GFP (n = 37), D7-L2 (n = 20; p = 0.002), anophelin (n = 19; p = 0.012), peroxidase (n = 19; p = 0.007), trio (n = 23; p = 0.7), SG2 precursor (n = 20; p = 0.02), 5' nucleotidase (n = 23; p = 0.05), where n is the total number of mosquitoes tested for each gene and p value is determined by the Mann-Whitney test. Results from three independent sets of experiments are shown in the graph. The rank order (the cumulative number of mosquitoes from different replicas that completed their probing within this time-interval) was plotted against the time.

Full size image

The depletion of TRIO (p = 0.7) did not result in any significant changes either in blood-feeding behavior or probing time.

Conclusions

The mosquito salivary glands play essential roles in blood-feeding and pathogen transmission. We have used a genome-wide microarray approach to characterize the A. gambiae salivary gland transcriptome, which corresponds to 4,978 genes (Figure 1), thereby extending earlier work based on cDNA/EST library clones [1, 21, 57]. Through the use of a full-genome microarray, we have obtained information on salivary gland transcript abundance and assessed the impact of blood feeding on salivary gland gene transcript abundance. The 93 transcripts showing significant differential abundance between the salivary glands of unfed and blood-fed female mosquitoes were mainly involved in metabolic, immunity, and oxidative stress-related processes (Figure 1D). It will be interesting to know whether the salivary gland expressed immune genes are involved in the interaction with, or killing of, the Plasmodium sporozoites. Interestingly, the majority of the secretory salivary gland-expressed transcripts displayed a lower abundance after blood-feeding, suggesting that their products may have been employed for protein synthesis during feeding and eventually degraded, or their transcription had been shut down. The lower blood-feeding capacity and increased probing time of mosquitoes upon silencing of the D7L2, anophelin, peroxidase 5B, SG2 precursor, and 5' nucleotidase genes is consistent with their role as anti-hemostatic factors (Figure 2A and Figure 3) [1, 62].

A better understanding of salivary gland gene expression and function can contribute to the development of malaria control strategies based on blocking the Plasmodium parasite in genetically modified mosquitoes or transmission blocking vaccines that could inhibit infection of the glands [66]. Furthermore, our studies indicate that the salivary glands of Anopheles mosquitoes express a variety of pharmacologically active substances that may be of value for designing novel anti-thrombotic or anti-platelet drugs for the treatment of cardiovascular disorders.

Methods

Rearing of A. gambiae mosquitoes

A. gambiae Keele strain mosquitoes were reared at 27°C in 70% humidity, and adults were maintained on a 10% sucrose solution on a 12-h light/dark cycle, according to standard rearing conditions [67].

RNA isolation and quantitative real-time PCR (qRT-PCR)

RNA was extracted and quantified (in triplicate samples) from dissected salivary glands (~80-100) using an RNeasy kit (Qiagen, Valencia, California, USA) according to standard methodology [68]. The qRT-PCR assays were performed as previously described [68], and the ribosomal protein S7 gene was used for normalization of the cDNA templates. The gene silencing efficiency (see below) was calculated according to the standard EΔΔCt method [69].

Microarray analysis

A. gambiae mosquitoes were allowed to feed for 20 min on Swiss Webster strain mice that had been anesthetized with a ketamine solution. For the microarray transcription analysis, approximately 100 salivary glands in triplicate were dissected 2 h after blood feeding (from both blood-fed and unfed mosquitoes). The blood-fed mosquitoes were not provided sugar solution for 2 hrs after blood-feeding, after which the salivary glands were dissected. The control unfed mosquitoes were treated in the same way. It should be noted that the triplicate assays were performed on different cohorts and generations. For the transcriptomic studies involving blood-fed regulated salivary gland gene expression, probes prepared from the RNA of salivary glands from blood-fed mosquitoes (labeled with Cy-5-dUTP fluorescent nucleotides) were co-hybridized with probes prepared from the RNA (labeled with Cy-3-dUTP fluorescent nucleotides) obtained from salivary glands of unfed mosquitoes. An Agilent custom-made 60-mer oligonucleotide microarray representing the full A. gambiae transcriptome was used for the assays. The Cy-5- and Cy-3-labeled probes were synthesized using the Low RNA Input Fluorescent Linear Amplification Kit from Agilent Technologies (Santa Clara, California, USA), and the slides were scanned using a GenePix 7000 autoloader scanner (Axon Instruments), as previously described in Dong et al., 2006. Scanning and data analyses were performed as previously described, with a signal cut-off intensity of 100 to remove low-intensity/poorly hybridized spots from the analysis. Loc-Fit normalization (LOWESS) was performed independently for each data set. The TMEV software, available online (http://www.tm4.org/mev/), was then used for t-test analysis at a significance level of p < 0.05, to evaluate the variability between biological replicates. Only transcripts that had signal values above or below a log2 cutoff value of ± 0.75 were considered significantly regulated with blood-feeding and used for further analysis (Figure 1D). If a particular gene was not regulated in the same direction in all the three biological replicates, it was considered to be non-significantly regulated. The GEO-NCBI accession numbers are: GSM573605, GSM573606 and GSM573607.

For the determination of the global A. gambiae female salivary gland transcriptome, the microarray spot hybridization fluorescence intensities were used as an indicator of transcript abundance/expression level, since each spot contained a similar amount of nucleic acid. First, the spot fluorescence intensity values from the three replicate assays were averaged for each spot and then for each gene (some genes were represented by multiple spots), generating a list of transcripts that were expressed in the A. gambiae female salivary gland. Transcripts were then categorized in groups on the basis of their relative transcript abundance, as measured by the mean fluorescence value (at 635 nm) of each transcript. This analysis produced the following three groups: i) highly abundant transcripts (fluorescence intensity of 5,000 to the maximum threshold of 65,000), ii) moderately expressed transcripts (1,000 to 4,999) and iii) poorly expressed transcripts (100 to 999) (listed as three different Excel spreadsheets in Additional file 1). The average fluorescence values at both 635 and 532 nm from the three GPR files (Genepix) are listed in Additional file 1, for direct comparison of each transcript before and after blood-feeding. The lists of genes were then further sub-grouped into different functional groups according to their known functions (Figure 1). The microarray validation by independent quantitative real-time PCR (qRT-PCR) on some of the genes was performed (data not shown) as described in our earlier publications [68, 70]. The validation data showed a high degree of correlation between the microarray and qRT-PCR values (Pearson correlation coefficient, p = 0.91; best-fit linear regression, R2 = 0.84; slope of the regression line = 0.85).

RNAi gene-silencing assays in adult mosquitoes

RNA interference (RNAi) assays in adult female mosquitoes were performed according to established RNAi methodology [71], using GFP dsRNA as the control. At least three independent RNAi-mediated gene-silencing assays were performed for each gene on different days, with two replicates each time; ~ 35-40 mosquitoes were injected for each replicate, and 25 randomly selected surviving mosquitoes from each replicate were used for further blood-feeding assays. It should be noted that the triplicate assays were conducted on different days, and the mosquitoes were from different cohorts/generations.

To silence the transcripts in the salivary glands, at least 8-9 times more dsRNA (~1.6-1.8 μg per mosquito) than normal (~0.2 μg) was injected [9] in order to obtain a better silencing efficiency in this body part, which is less accessible to diffusion of dsRNA from the hemolymph. The primer sequences of the selected genes and the GFP used to generate dsRNA are listed in Additional file 5. For gene silencing verification by qRT-PCR, ~40-50 mosquitoes were injected with dsRNA, and the salivary glands were dissected 4 days later. The ribosomal protein S7 gene was used for normalization of the cDNA templates (primer sequences are provided in Additional file 5). The fold difference in transcript abundance levels (the silencing efficiency) after silencing was calculated according to the standard EΔΔCt method [69]. The primers used for silencing verification are presented in Additional file 6.

Blood-feeding assays

To study the blood-feeding propensity upon gene silencing, mosquitoes were exposed to anesthetized (with ketamine solution) Swiss Webster mice 4 days after the dsRNA injections. The mosquitoes were allowed to feed for 20 min, after which we scored the number of mosquitoes that had fed. Two mice were required for each feeding experiment (GFP control and experimental), and were used simultaneously by placing each mouse with its head over one cage and its tail over the other. Every 5 min, the positions of the mice were flipped (with minimum interruption) to ensure that the two mosquito sets drew blood from the identical body parts of the two mice. For statistical analysis of the gene silencing phenotypes (blood-feeding propensity) with respect to GFP control, a Mann Whitney test was used, as shown in our earlier publications [68, 72]; the resulting p-values are given in Figure 2A and Additional file 4.

Probing assays

Three days-old adult female mosquitoes were injected with dsRNA as described above and their probing time was assayed four days later. The mosquitoes were starved from sugar for 4-5 hours before conducting the probing time measurements as described in Lombardo et al., 2009. A single mosquito was transferred in a small cage and allowed to rest for 30 min before being offered an anesthetized (with ketamine solution) Swiss Webster mouse. Probing time is defined as the time taken from the initial insertion of the mouthparts in the skin until the initial engorgement of blood [9, 11, 73]. If the mosquito terminates the probing unsuccessfully and tries again elsewhere, the second probing time is progressively added to the first one until blood is ingested. The time between two probings (interprobing time), is not included. The assay was conducted for 500 seconds for each mosquito. Three independent sets of dsRNA injections for probing time were performed for each of the six genes and were compared to the probing time of GFP dsRNA injected control mosquitoes. The results from all replica assays were pooled and the rank order of the mosquitoes (the cumulative number of mosquito pooled from different replicas that completes their probing, at any particular time within the 500 seconds period) was plotted against time (seconds), as performed in earlier studies [9, 11]. The data was then analyzed by the Mann-Whitney test for statistical significance as described in Lombardo et al., 2009. The total number of mosquitoes tested (n) and the resulting p-values for each silenced transcript and GFP dsRNA injected control are shown in Figure 3 legend.

Abbreviations

GFP:

green fluorescent protein

qRT-PCR:

quantative real-time polymerase chain reaction

RNAi:

RNA interference

OBP:

odorant binding proteins

dsRNA:

double-strand RNA.

References

  1. Arca B, Lombardo F, Valenzuela JG, Francischetti IM, Marinotti O, Coluzzi M, Ribeiro JM: An updated catalogue of salivary gland transcripts in the adult female mosquito, Anopheles gambiae. J Exp Biol. 2005, 208: 3971-3986. 10.1242/jeb.01849.

    Article  CAS  PubMed  Google Scholar 

  2. Bolognesi ML, Minarini A, Budriesi R, Cacciaguerra S, Chiarini A, Spampinato S, Tumiatti V, Melchiorre C: Universal template approach to drug design: polyamines as selective muscarinic receptor antagonists. J Med Chem. 1998, 41: 4150-4160. 10.1021/jm981038d.

    Article  CAS  PubMed  Google Scholar 

  3. Lanfrancotti A, Lombardo F, Santolamazza F, Veneri M, Castrignano T, Coluzzi M, Arca B: Novel cDNAs encoding salivary proteins from the malaria vector Anopheles gambiae. FEBS Lett. 2002, 517: 67-71. 10.1016/S0014-5793(02)02578-4.

    Article  CAS  PubMed  Google Scholar 

  4. Monteiro RQ: Targeting exosites on blood coagulation proteases. An Acad Bras Cienc. 2005, 77: 275-280.

    Article  CAS  PubMed  Google Scholar 

  5. Baldacci P, Menard R: The elusive malaria sporozoite in the mammalian host. Mol Microbiol. 2004, 54: 298-306. 10.1111/j.1365-2958.2004.04275.x.

    Article  CAS  PubMed  Google Scholar 

  6. Frevert U, Engelmann S, Zougbede S, Stange J, Ng B, Matuschewski K, Liebes L, Yee H: Intravital observation of Plasmodium berghei sporozoite infection of the liver. PLoS Biol. 2005, 3: e192-10.1371/journal.pbio.0030192.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Ribeiro JM: Vector salivation and parasite transmission. Mem Inst Oswaldo Cruz. 1987, 82 (Suppl 3): 1-3.

    Article  PubMed  Google Scholar 

  8. Stark KR, James AA: Anticoagulants in vector arthropods. Parasitol Today. 1996, 12: 430-437. 10.1016/0169-4758(96)10064-8.

    Article  CAS  PubMed  Google Scholar 

  9. Boisson B, Jacques JC, Choumet V, Martin E, Xu J, Vernick K, Bourgouin C: Gene silencing in mosquito salivary glands by RNAi. FEBS Lett. 2006, 580: 1988-1992. 10.1016/j.febslet.2006.02.069.

    Article  CAS  PubMed  Google Scholar 

  10. Ribeiro JM, Rossignol PA, Spielman A: Role of mosquito saliva in blood vessel location. J Exp Biol. 1984, 108: 1-7.

    Article  CAS  PubMed  Google Scholar 

  11. Lombardo F, Ronca R, Rizzo C, Mestres-Simon M, Lanfrancotti A, Curra C, Fiorentino G, Bourgouin C, Ribeiro JM, Petrarca V: The Anopheles gambiae salivary protein gSG6: an anopheline-specific protein with a blood-feeding role. Insect Biochem Mol Biol. 2009, 39: 457-466. 10.1016/j.ibmb.2009.04.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Golenda CF, Klein T, Coleman R, Burge R, Ward RA, Seeley DC: Depletion of total salivary gland protein in blood-fed Anopheles mosquitoes. J Med Entomol. 1995, 32: 300-305.

    Article  CAS  PubMed  Google Scholar 

  13. Siriyasatien P, Tangthongchaiwiriya K, Kraivichian K, Nuchprayoon S, Tawatsin A, Thavara U: Decrease of mosquito salivary gland proteins after a blood meal: an implication for pathogenesis of mosquito bite allergy. J Med Assoc Thai. 2005, 88 (Suppl 4): S255-S259.

    PubMed  Google Scholar 

  14. Belkaid Y, Valenzuela JG, Kamhawi S, Rowton E, Sacks DL, Ribeiro JM: Delayed-type hypersensitivity to Phlebotomus papatasi sand fly bite: An adaptive response induced by the fly?. Proc Natl Acad Sci USA. 2000, 97: 6704-6709. 10.1073/pnas.97.12.6704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Valenzuela JG, Belkaid Y, Garfield MK, Mendez S, Kamhawi S, Rowton ED, Sacks DL, Ribeiro JM: Toward a defined anti-Leishmania vaccine targeting vector antigens: characterization of a protective salivary protein. J Exp Med. 2001, 194: 331-342. 10.1084/jem.194.3.331.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Arca B, Lombardo F, de Lara CM, della TA, Dimopoulos G, James AA, Coluzzi M: Trapping cDNAs encoding secreted proteins from the salivary glands of the malaria vector Anopheles gambiae. Proc Natl Acad Sci USA. 1999, 96: 1516-1521. 10.1073/pnas.96.4.1516.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Arca B, Lombardo F, Francischetti IM, Pham VM, Mestres-Simon M, Andersen JF, Ribeiro JM: An insight into the sialome of the adult female mosquito Aedes albopictus. Insect Biochem Mol Biol. 2007, 37: 107-127. 10.1016/j.ibmb.2006.10.007.

    Article  CAS  PubMed  Google Scholar 

  18. Calvo E, Andersen J, Francischetti IM, deL CM, deBianchi AG, James AA, Ribeiro JM, Marinotti O: The transcriptome of adult female Anopheles darlingi salivary glands. Insect Mol Biol. 2004, 13: 73-88. 10.1111/j.1365-2583.2004.00463.x.

    Article  CAS  PubMed  Google Scholar 

  19. Francischetti IM, Valenzuela JG, Pham VM, Garfield MK, Ribeiro JM: Toward a catalog for the transcripts and proteins (sialome) from the salivary gland of the malaria vector Anopheles gambiae. J Exp Biol. 2002, 205: 2429-2451.

    Article  CAS  PubMed  Google Scholar 

  20. Ribeiro JM, Charlab R, Pham VM, Garfield M, Valenzuela JG: An insight into the salivary transcriptome and proteome of the adult female mosquito Culex pipiens quinquefasciatus. Insect Biochem Mol Biol. 2004, 34: 543-563. 10.1016/j.ibmb.2004.02.008.

    Article  CAS  PubMed  Google Scholar 

  21. Ribeiro JM, Arca B, Lombardo F, Calvo E, Phan VM, Chandra PK, Wikel SK: An annotated catalogue of salivary gland transcripts in the adult female mosquito, Aedes aegypti. BMC Genomics. 2007, 8: 6-10.1186/1471-2164-8-6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Neira OM, Ribeiro JM, Heyland A, Vanekeris L, Moroz T, Linser PJ: The salivary transcriptome of Anopheles gambiae (Diptera: Culicidae) larvae: A microarray-based analysis. Insect Biochem Mol Biol. 2009, 39 (5-6): 382-94. 10.1016/j.ibmb.2009.03.001.

    Article  CAS  Google Scholar 

  23. Choumet V, Carmi-Leroy A, Laurent C, Lenormand P, Rousselle JC, Namane A, Roth C, Brey PT: The salivary glands and saliva of Anopheles gambiae as an essential step in the Plasmodium life cycle: a global proteomic study. Proteomics. 2007, 7: 3384-3394. 10.1002/pmic.200700334.

    Article  CAS  PubMed  Google Scholar 

  24. Ghosh AK, Devenport M, Jethwaney D, Kalume DE, Pandey A, Anderson VE, Sultan AA, Kumar N, Jacobs-Lorena M: Malaria parasite invasion of the mosquito salivary gland requires interaction between the Plasmodium TRAP and the Anopheles saglin proteins. PLoS Pathog. 2009, 5: e1000265-10.1371/journal.ppat.1000265.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Okulate MA, Kalume DE, Reddy R, Kristiansen T, Bhattacharyya M, Chaerkady R, Pandey A, Kumar N: Identification and molecular characterization of a novel protein Saglin as a target of monoclonal antibodies affecting salivary gland infectivity of Plasmodium sporozoites. Insect Mol Biol. 2007, 16: 711-722. 10.1111/j.1365-2583.2007.00765.x.

    Article  CAS  PubMed  Google Scholar 

  26. Nickel W: The mystery of nonclassical protein secretion. A current view on cargo proteins and potential export routes. Eur J Biochem. 2003, 270: 2109-2119. 10.1046/j.1432-1033.2003.03577.x.

    Article  CAS  PubMed  Google Scholar 

  27. Koutsos AC, Blass C, Meister S, Schmidt S, Maccallum RM, Soares MB, Collins FH, Benes V, Zdobnov E, Kafatos FC: Life cycle transcriptome of the malaria mosquito Anopheles gambiae and comparison with the fruitfly Drosophila melanogaster. Proc Natl Acad Sci USA. 2007, 104: 11304-11309. 10.1073/pnas.0703988104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Neira OM, Vanekeris L, Corena-McLeod MD, Linser PJ: A microarray-based analysis of transcriptional compartmentalization in the alimentary canal of Anopheles gambiae (Diptera: Culicidae) larvae. Insect Mol Biol. 2008, 17: 61-72. 10.1111/j.1365-2583.2008.00779.x.

    Article  Google Scholar 

  29. Warr E, Aguilar R, Dong Y, Mahairaki V, Dimopoulos G: Spatial and sex-specific dissection of the Anopheles gambiae midgut transcriptome. BMC Genomics. 2007, 8: 37-10.1186/1471-2164-8-37.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Das S, Dimopoulos G: Molecular analysis of photic inhibition of blood-feeding in Anopheles gambiae. BMC Physiol. 2008, 8: 23-10.1186/1472-6793-8-23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. So WV, Sarov-Blat L, Kotarski CK, McDonald MJ, Allada R, Rosbash M: Takeout, a novel Drosophila gene under circadian clock transcriptional regulation. Mol Cell Biol. 2000, 20: 6935-6944. 10.1128/MCB.20.18.6935-6944.2000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Hardin PE: The circadian timekeeping system of Drosophila. Curr Biol. 2005, 15: R714-R722. 10.1016/j.cub.2005.08.019.

    Article  CAS  PubMed  Google Scholar 

  33. Xu PX, Zwiebel LJ, Smith DP: Identification of a distinct family of genes encoding atypical odorant-binding proteins in the malaria vector mosquito, Anopheles gambiae. Insect Mol Biol. 2003, 12: 549-560. 10.1046/j.1365-2583.2003.00440.x.

    Article  CAS  PubMed  Google Scholar 

  34. 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: 36626-36633. 10.1074/jbc.M706410200.

    Article  CAS  PubMed  Google Scholar 

  35. Calvo E, Mans BJ, Andersen JF, Ribeiro JM: Function and evolution of a mosquito salivary protein family. J Biol Chem. 2006, 281: 1935-1942. 10.1074/jbc.M510359200.

    Article  CAS  PubMed  Google Scholar 

  36. James AA, Blackmer K, Marinotti O, Ghosn CR, Racioppi JV: Isolation and characterization of the gene expressing the major salivary gland protein of the female mosquito, Aedes aegypti. Mol Biochem Parasitol. 1991, 44: 245-253. 10.1016/0166-6851(91)90010-4.

    Article  CAS  PubMed  Google Scholar 

  37. Valenzuela JG, Charlab R, Gonzalez EC, de Miranda-Santos IK, Marinotti O, Francischetti IM, Ribeiro JM: The D7 family of salivary proteins in blood sucking diptera. Insect Mol Biol. 2002, 11: 149-155. 10.1046/j.1365-2583.2002.00319.x.

    Article  CAS  PubMed  Google Scholar 

  38. Arca B, Lombardo F, Lanfrancotti A, Spanos L, Veneri M, Louis C, Coluzzi M: A cluster of four D7-related genes is expressed in the salivary glands of the African malaria vector Anopheles gambiae. Insect Mol Biol. 2002, 11: 47-55. 10.1046/j.0962-1075.2001.00309.x.

    Article  CAS  PubMed  Google Scholar 

  39. Isawa H, Yuda M, Orito Y, Chinzei Y: A mosquito salivary protein inhibits activation of the plasma contact system by binding to factor XII and high molecular weight kininogen. J Biol Chem. 2002, 277: 27651-27658. 10.1074/jbc.M203505200.

    Article  CAS  PubMed  Google Scholar 

  40. Calvo E, Mans BJ, Ribeiro JM, Andersen JF: Multifunctionality and mechanism of ligand binding in a mosquito antiinflammatory protein. Proc Natl Acad Sci USA. 2009, 106: 3728-3733. 10.1073/pnas.0813190106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Dimopoulos G, Seeley D, Wolf A, Kafatos FC: Malaria infection of the mosquito Anopheles gambiae activates immune-responsive genes during critical transition stages of the parasite life cycle. EMBO J. 1998, 17: 6115-6123. 10.1093/emboj/17.21.6115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Rosinski-Chupin I, Briolay J, Brouilly P, Perrot S, Gomez SM, Chertemps T, Roth CW, Keime C, Gandrillon O, Couble P: SAGE analysis of mosquito salivary gland transcriptomes during Plasmodium invasion. Cell Microbiol. 2007, 9: 708-724. 10.1111/j.1462-5822.2006.00822.x.

    Article  CAS  PubMed  Google Scholar 

  43. Van Steijn GJ, Amerongen AV, Veerman EC, Kasanmoentalib S, Overdijk B: Chitinase in whole and glandular human salivas and in whole saliva of patients with periodontal inflammation. Eur J Oral Sci. 1999, 107: 328-337. 10.1046/j.0909-8836.1999.eos107503.x.

    Article  CAS  PubMed  Google Scholar 

  44. Rossignol PA, Lueders AM: Bacteriolytic factor in the salivary glands of Aedes aegypti. Comp Biochem Physiol B. 1986, 83: 819-822. 10.1016/0305-0491(86)90153-7.

    Article  CAS  PubMed  Google Scholar 

  45. Christophides GK, Zdobnov E, Barillas-Mury C, Birney E, Blandin S, Blass C, Brey PT, Collins FH, Danielli A, Dimopoulos G: Immunity-related genes and gene families in Anopheles gambiae. Science. 2002, 298: 159-165. 10.1126/science.1077136.

    Article  CAS  PubMed  Google Scholar 

  46. Dimopoulos G, Christophides GK, Meister S, Schultz J, White KP, Barillas-Mury C, Kafatos FC: Genome expression analysis of Anopheles gambiae: responses to injury, bacterial challenge, and malaria infection. Proc Natl Acad Sci USA. 2002, 99: 8814-8819. 10.1073/pnas.092274999.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Osta MA, Christophides GK, Vlachou D, Kafatos FC: Innate immunity in the malaria vector Anopheles gambiae: comparative and functional genomics. J Exp Biol. 2004, 207: 2551-2563. 10.1242/jeb.01066.

    Article  CAS  PubMed  Google Scholar 

  48. Vlachou D, Schlegelmilch T, Christophides GK, Kafatos FC: Functional genomic analysis of midgut epithelial responses in Anopheles during Plasmodium invasion. Curr Biol. 2005, 15: 1185-1195. 10.1016/j.cub.2005.06.044.

    Article  CAS  PubMed  Google Scholar 

  49. Waterhouse RM, Kriventseva EV, Meister S, Xi Z, Alvarez KS, Bartholomay LC, Barillas-Mury C, Bian G, Blandin S, Christensen BM: Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science. 2007, 316: 1738-1743. 10.1126/science.1139862.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. David JP, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, Louis C, Hemingway J, Ranson H: The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proc Natl Acad Sci USA. 2005, 102: 4080-4084. 10.1073/pnas.0409348102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Lund AA, Rhoads DM, Lund AL, Cerny RL, Elthon TE: In vivo modifications of the maize mitochondrial small heat stress protein, HSP22. J Biol Chem. 2001, 276: 29924-29929. 10.1074/jbc.M103373200.

    Article  CAS  PubMed  Google Scholar 

  52. Champagne DE, Smartt CT, Ribeiro JM, James AA: The salivary gland-specific apyrase of the mosquito Aedes aegypti is a member of the 5'-nucleotidase family. Proc Natl Acad Sci USA. 1995, 92: 694-698. 10.1073/pnas.92.3.694.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Ribeiro JM, Francischetti IM: Role of arthropod saliva in blood feeding: sialome and post-sialome perspectives. Annu Rev Entomol. 2003, 48: 73-88. 10.1146/annurev.ento.48.060402.102812.

    Article  CAS  PubMed  Google Scholar 

  54. Robertus JD, Monzingo AF, Marcotte EM, Hart PJ: Structural analysis shows five glycohydrolase families diverged from a common ancestor. J Exp Zool. 1998, 282: 127-132. 10.1002/(SICI)1097-010X(199809/10)282:1/2<127::AID-JEZ14>3.0.CO;2-R.

    Article  CAS  PubMed  Google Scholar 

  55. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Calvo E, Dao A, Pham VM, Ribeiro JM: An insight into the sialome of Anopheles funestus reveals an emerging pattern in anopheline salivary protein families. Insect Biochem Mol Biol. 2007, 37: 164-175. 10.1016/j.ibmb.2006.11.005.

    Article  CAS  PubMed  Google Scholar 

  57. Marinotti O, Nguyen QK, Calvo E, James AA, Ribeiro JM: Microarray analysis of genes showing variable expression following a blood meal in Anopheles gambiae. Insect Mol Biol. 2005, 14: 365-373. 10.1111/j.1365-2583.2005.00567.x.

    Article  CAS  PubMed  Google Scholar 

  58. Thangamani S, Wikel SK: Differential expression of Aedes aegypti salivary transcriptome upon blood feeding. Parasit Vectors. 2009, 2: 34-10.1186/1756-3305-2-34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Kerlin RL, Hughes S: Enzymes in saliva from four parasitic arthropods. Med Vet Entomol. 1992, 6: 121-126. 10.1111/j.1365-2915.1992.tb00587.x.

    Article  CAS  PubMed  Google Scholar 

  60. Qiu XB, Shao YM, Miao S, Wang L: The diversity of the DnaJ/Hsp40 family, the crucial partners for Hsp70 chaperones. Cell Mol Life Sci. 2006, 63: 2560-2570. 10.1007/s00018-006-6192-6.

    Article  CAS  PubMed  Google Scholar 

  61. Ribeiro JM, Valenzuela JG: Purification and cloning of the salivary peroxidase/catechol oxidase of the mosquito Anopheles albimanus. J Exp Biol. 1999, 202: 809-816.

    Article  CAS  PubMed  Google Scholar 

  62. Francischetti IM, Valenzuela JG, Ribeiro JM: Anophelin: kinetics and mechanism of thrombin inhibition. Biochemistry. 1999, 38: 16678-16685. 10.1021/bi991231p.

    Article  CAS  PubMed  Google Scholar 

  63. Calvo E, Tokumasu F, Marinotti O, Villeval JL, Ribeiro JM, Francischetti IM: Aegyptin, a novel mosquito salivary gland protein, specifically binds to collagen and prevents its interaction with platelet glycoprotein VI, integrin alpha2beta1, and von Willebrand factor. J Biol Chem. 2007, 282: 26928-26938. 10.1074/jbc.M705669200.

    Article  CAS  PubMed  Google Scholar 

  64. Yoshida S, Sudo T, Niimi M, Tao L, Sun B, Kambayashi J, Watanabe H, Luo E, Matsuoka H: Inhibition of collagen-induced platelet aggregation by anopheline antiplatelet protein, a saliva protein from a malaria vector mosquito. Blood. 2008, 111: 2007-2014. 10.1182/blood-2007-06-097824.

    Article  CAS  PubMed  Google Scholar 

  65. Lin MZ, Greenberg ME: Orchestral maneuvers in the axon: trio and the control of axon guidance. Cell. 2000, 101: 239-242. 10.1016/S0092-8674(00)80833-8.

    Article  CAS  PubMed  Google Scholar 

  66. Terenius O, Marinotti O, Sieglaff D, James AA: Molecular genetic manipulation of vector mosquitoes. Cell Host Microbe. 2008, 4: 417-423. 10.1016/j.chom.2008.09.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Benedict MQ: Care and maintenance of Anopheline mosquitoes. The molecular biology of disease vectors: A methods manual. Edited by: Crampton JM. 1997, London: Champman & Hall, 3-12.

    Chapter  Google Scholar 

  68. Dong Y, Aguilar R, Xi Z, Warr E, Mongin E, Dimopoulos G: Anopheles gambiae immune responses to human and rodent Plasmodium parasite species. PLoS Pathog. 2006, 2: e52-10.1371/journal.ppat.0020052.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Aguilar R, Das S, Dong Y, Dimopoulos G: Continuous exposure to Plasmodium results in decreased susceptibility and transcriptomic divergence of the Anopheles gambiae immune system. BMC Genomics. 2007, 8: 451-10.1186/1471-2164-8-451.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Blandin S, Moita LF, Kocher T, Wilm M, Kafatos FC, Levashina EA: Reverse genetics in the mosquito Anopheles gambiae: targeted disruption of the Defensin gene. EMBO Rep. 2002, 3: 852-856. 10.1093/embo-reports/kvf180.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Warr E, Das S, Dong Y, Dimopoulos G: The Gram-negative bacteria-binding protein gene family: its role in the innate immune system of Anopheles gambiae and in anti-Plasmodium defence. Insect Mol Biol. 2008, 17: 39-51. 10.1111/j.1365-2583.2008.00778.x.

    Article  CAS  PubMed  Google Scholar 

  73. Ribeiro JM: Blood-feeding in mosquitoes: probing time and salivary gland anti-haemostatic activities in representatives of three genera (Aedes, Anopheles, Culex). Med Vet Entomol. 2000, 14: 142-148. 10.1046/j.1365-2915.2000.00227.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the microarray core facility and the insectary personnel at the Johns Hopkins Malaria Research Institute. We also thank Dr. Fotis C. Kafatos and Dr. Jose M. Ribeiro for their discussions and fruitful suggestions and Dr. Deborah McClellan for editorial assistance. This work was supported by the Ellison Medical Foundation, the National Institutes of Health/National Institute for Allergy and Infectious Disease 1R01AI061576, the United Nations Development Program/World Bank/World Health Organization Special Program for Research and Training in Tropical Diseases, and the Johns Hopkins School of Public Health.

Author information

Authors and Affiliations

  1. W. Harry Feinstone Department of Molecular Microbiology and Immunology, Bloomberg School of Public Health, Johns Hopkins University, 615 N Wolfe Street, Baltimore, MD, 21205-2179, USA

    Suchismita Das, Andrea Radtke, Antonio M Mendes & George Dimopoulos

  2. Department of Pathobiological Sciences, University of Wisconsin-Madison, 1656 Linden Drive, Madison, WI, 53706, USA

    Young-Jun Choi

  3. Imperial College London, Division of Cell and Molecular Biology, Faculty of Natural Sciences, South Kensington Campus, London, UK

    Antonio M Mendes

  4. Laboratory of Malaria and Vector Research, NIAID, National Institutes of Health, Rockville, Maryland, 20852, USA

    Jesus G Valenzuela

Authors
  1. Suchismita Das
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Andrea Radtke
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Young-Jun Choi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Antonio M Mendes
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jesus G Valenzuela
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. George Dimopoulos
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to George Dimopoulos.

Additional information

Authors' contributions

SD and AR conducted the microarray and gene silencing assays and SD contributed to the writing of the manuscript. YJC, AM and JGV performed preliminary and confirmatory analyses and assays and edited the manuscript. GD designed the experiments and the oligonucleotide microarrays and contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2010_9957_MOESM1_ESM.XLS

Additional file 1: List of genes expressed in the female salivary gland of A. gambiae, along with their previous (ENSANGT) and recent (AGAP-R) transcript IDs. Note that the recent AGAP-R IDs for some transcripts are missing in VectorBase and have therefore been left blank in column B. The fluorescent values at 635 nm (Cy-5 dye) and 532 nm (Cy-3 dye) for each gene are shown in column E and F respectively. The excel file has three sheets: Sheet 1 lists the genes that are highly expressed, with fluorescent values (at 635 nm) from 5,000 to the maximum threshold of 65,000, Sheet 2 lists the moderately expressed genes with fluorescent values (at 635 nm) between 1,000 to 4,999 and Sheet 3 lists the weakly expressed genes with fluorescent values (at 635 nm) between 100 to 999. The gene list is further sorted into different functional groups, according to the predicted function of the genes [CS: cytoskeletal and structural; CIR: circadian; CSR: chemosensory; PROT/DIG: proteolytic digestion; IMM: immunity; MET/ENZ: metabolism and enzymes; RED/STR: redox/stress; R/T/T: replication/transcription and translation; TRP: transport]. The list of genes in each group is further sorted into descending order of the fluorescent value at 635 nm (column E). (XLS 846 KB)

12864_2010_9957_MOESM2_ESM.XLS

Additional file 2: List of genes that are regulated upon blood-feeding in the female A. gambiae salivary gland, along with their ENSANGT and more recent AGAP-R transcript IDs. Transcripts lacking AGAP-R IDs are only indicated with the ENSANT transcript ID. Genes are only listed if they display significant differential regulation between fed and unfed mosquitoes, at a log2 fold ratio above or below the ± 0.75 threshold. Predicted functional groups classification for genes are also presented [CS: cytoskeletal and structural; CIR: circadian; CSR: chemosensory; PROT/DIG: proteolytic digestion; IMM: immunity; MET/ENZ: metabolism and enzymes; RED/STR: redox/stress; R/T/T: replication/transcription and translation; TRP: transport]. (XLS 31 KB)

12864_2010_9957_MOESM3_ESM.DOC

Additional file 3: Gene silencing efficiency of the 10 selected genes. The gene transcript abundance in GFP dsRNA treated mosquitoes was set to 1.0, and the corresponding percentage silencing was determined by qRT-PCR. The AgS7 gene was used for normalization of cDNA templates. The standard error values are shown. (DOC 36 KB)

12864_2010_9957_MOESM4_ESM.DOC

Additional file 4: The blood-feeding propensity of mosquitoes in gene silenced mosquitoes. Ten salivary gland expressed genes were individually silenced in mosquitoes, and the mosquitoes were allowed to feed on blood, 4 days after dsRNA injections. The percentages of mosquitoes that feed ("Fed") were scored and are shown along with their standard error values. The p-value from statistical analysis of Mann Whitney Test is also shown. (DOC 38 KB)

12864_2010_9957_MOESM5_ESM.DOC

Additional file 5: Primers used for synthesis of dsRNAs for RNAi gene silencing assays. The first 20 bases in bold correspond to the T7 polymerase promoter site. The transcript ID numbers (AGAP-RA from VectorBase) are also shown for each gene. The list also includes the sequences of GFP primers used for generating dsRNA and the A. gambiae S7 primer (used for standardization of cDNA templates). (DOC 36 KB)

12864_2010_9957_MOESM6_ESM.DOC

Additional file 6: Primers used for verification of RNAi silencing. The sense primers were newly designed (VERI Sense), and the antisense primers used were the same as those used for making the corresponding dsRNAs. The transcript ID numbers (AGAP-RA from VectorBase) are also shown for each gene. (DOC 34 KB)

Authors’ original submitted files for images

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Das, S., Radtke, A., Choi, YJ. et al. Transcriptomic and functional analysis of the Anopheles gambiae salivary gland in relation to blood feeding. BMC Genomics 11, 566 (2010). https://doi.org/10.1186/1471-2164-11-566

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-566

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us