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

RETRACTED ARTICLE: Species-specific chemosensory gene expression in the olfactory organs of the malaria vector Anopheles gambiae

This article was retracted on 03 August 2015

Abstract

Background

The malaria mosquito Anopheles gambiae has a high preference for human hosts, a characteristic that contributes greatly to its capacity for transmitting human malaria. A sibling species, An. quadriannulatus, has a quite different host preference and feeds mostly on bovids. For this reason it does not contribute to human malaria transmission. Host seeking in mosquitoes is modulated by the olfactory system, which is primarily housed in the antennae and maxillary palps. Therefore, the detection of differing host odors by sibling species may be reflected in the expression level of the olfactory genes involved. Accordingly, we compared the transcriptomes of the antennae and maxillary palps of An. gambiae and An. quadriannulatus.

Results

We identified seven relatively abundant olfactory receptors, nine ionotropic receptors and three odorant binding proteins that are substantially up-regulated in An. gambiae antennae. Interestingly, we find that the maxillary palps of An. gambiae contain a species-specific olfactory receptor, Or52, and five An. gambiae-specific gustatory receptors (AgGr48-52) that are relatively abundant. These five gustatory receptors are also expressed in An. gambiae antennae, although at lower level, indicating a likely role in olfaction, rather than gustation. We also document an approximately three-fold higher overall expression of olfaction genes in the maxillary palps of An. quadriannulatus, indicating an important role of this organ in the olfaction system of this species. Finally, the expression of the CO2 receptor genes is five to six-fold higher in the zoophilic An. quadriannulatus, implying a much higher sensitivity for detecting CO2.

Conclusions

These results identify potential human host preference genes in the malaria vector An. gambiae. Interestingly, species-specific expression of several gustatory receptors in the olfactory organs indicate a role in olfaction rather than gustation. Additionally, a more expansive role for maxillary palps in olfaction is implicated than previously thought, albeit more so in the zoophilic An. quadriannulatus.

Background

The malaria mosquitoes within the Anopheles gambiae complex vary considerably in their host preference. Africa’s main malaria vector An. gambiae s.s. is highly anthropophilic, whereas the zoophilic An. quadriannulatus rarely if ever attacks humans [1]. This preference of An. gambiae for human hosts is a major factor in its high vectorial capacity for human malaria parasites. Conversely, although the zoophilic An. quadriannulatus is a competent malaria vector [2], this species does not contribute to malaria transmission because it rarely feeds on human hosts in the field, although it does so readily in the lab [3, 4].

Mosquitoes’ host attraction is primarily modulated by the olfaction system and An. gambiae females are strongly attracted to emanations from human sweat. Volatiles produced by microflora on the surface of human skin are believed to be responsible for the uniqueness of human odor [5, 6]. Over 350 volatiles are found in human sweat [7], and while not all of these play a role in allowing An. gambiae to differentiate human hosts from others, it is likely that a blend of human volatiles is involved. For example, An. gambiae females are attracted to a mixture of ammonia, lactic acid, as well as a synergistic blend of ammonia, lactic acid and carboxylic acids [810]. Anopheles gambiae and An. quadriannulatus also show different sensitivities to various compounds found in human and animal sweat and/or breath. Therefore, the relative quantities of the constituents of host odor blends, rather than the presence or absence of specific volatiles, could be important in determining attractiveness to various species [11, 12].

The antennae and the maxillary palps, the two main olfactory appendages of An. gambiae[13, 14], are lined with sensilla that house the olfactory sensory neurons that express olfactory receptors (ORs) [15] or ionotropic receptors (IRs) [16, 17]. The binding of odorants to the ORs and IRs triggers the transduction cascade that sends a signal to the olfactory lobes in the cerebral ganglion of the insects [18]. Because of this direct interaction between the receptors and the odorants, differences in host preference between species may be reflected in differences in the expression or molecular structure of the receptors.

Currently, 76 Ors ([19, 20] and 44 Irs have been identified [21]). ORs are heteromeric ligand-gated ion channels encoded by the highly conserved co-receptor Orco and a specific Or. ORs differ in their tuning breadth and some ORs respond to either a single or small number of odorants, while others respond to a variety of volatiles [2226]. IRs are also heteromeric ligand-gated ion channels, but these can contain up to three different subunits that include one or two of the broadly expressed co-receptors Ir25a and Ir8a[16, 26].

In addition to the ORs and IRs, odorant binding proteins (OBPs) play a role in odorant recognition and interact directly with odorants. OBPs are small, water-soluble transport molecules abundant in the lymph of the sensilla. They transport hydrophobic odorants through the haemolymph to the receptors (reviewed in [27]). Currently, 57 putative Obps have been identified [2830], but only 34 Obps are expressed in female antennae [31]. Some OBPs almost certainly play a role in the transport of molecules outside the olfaction system, as two Obps are known to be expressed only in female heads [17].

Differences in the expression level of olfaction genes have been observed between closely related species feeding on different hosts. Expression levels of as many as 53% of the ORs and 55% of the OBPs in the antennae differed between the generalists D. melanogaster, D. simulans, and their specialist sister-species D. sechellia which feeds on the toxic Morinda citrifolia pairs. This is a significantly higher number than observed in other genes [32]. Although some of these changes may be due to neutral evolution, several genes have undergone a major change in expression level along the D. sechellia branch, and are thought to be associated with host shifts [32]. For example, Or22a is strongly up-regulated in D. sechellia. This receptor is sensitive to a compound emitted by the fruit of D. sechellia’s host plant Morinda citrifolia[33]. Additionally, D. sechellia lost six Or genes since its split from its generalist sister-species D. simulans, which lost none [34]. Furthermore, an increase in olfactory receptor loss was also associated with host specialization in D. erecta[35]. Recently, a comparison between the day-time transcriptome of An. gambiae and An. quadriannulatus antennae identified differences in olfaction gene expression that may be related to the difference in host preference between these sibling species [36]. It has been shown however that olfactory gene expression fluctuates across the circadian cycle [37]. Here, we compare the transcriptome of both the female antennae and palps of the anthropophilic An. gambiae and the zoophilic An. quadriannulatus, during the early dark cycle, when both species are actively seeking hosts [38, 39]. These comparisons further show the divergence of the olfactory organs of these two species, and allows us to identify species-specific chemosensory genes in An. gambiae that may be responsible for human host preference.

Results

Host choice assay

The attraction of An. gambiae and An. quadriannulatus laboratory strains to human odor vs cow odor was examined in a dual choice olfactometer. Consistent with the host preference of these species in the field and with recent work on laboratory colonies [40], An. gambiae was significantly attracted to human odor (77%, N = 770, p < 0.0001), whereas An. quadriannulatus significantly prefers cow odor (67%, N = 330, p = 0.0029). Therefore, the natural host preference of these species is largely preserved in strains kept in laboratory conditions for many generations.

Gene expression analyses

Three replicate female antennae RNAseq data sets and two replicate maxillary palps RNAseq data sets were obtained for both An. gambiae and An. quadriannulatus. After quality control screening, 91.0% of antennal reads from An. gambiae and 87.0% of antennal reads from An. quadriannulatus mapped to the An. gambiae reference genome. For palps, 86.8% of the An. gambiae reads and 84.3% of the An. quadriannulatus reads mapped back to the genome. A higher percentage of total reads obtained for the antennae mapped to a single location in An. gambiae vs An. quadriannulatus (83.7% vs 78.7%), whereas fewer reads from the palps mapped to only one location for this species (65.6% vs 76.4%). Additionally, the mapping software reported that no An. quadriannulatus reads remained unmapped due to mismatches with the reference genome, hence the difference in read mapping is not due to a divergence between the genomes of the two species.

We obtained 58.7 to 79.8 million mapped reads for each of the six antennal samples, for a total of 429.5 million mapped reads. Between 52.3 and 75.0 million mapped reads were obtained for each of the four maxillary palp samples, for a total of 261.9 million. Clustering of the variance-stabilized transformed counts shows that for both antennae and maxillary palp samples there was relatively little variation among biological replicates relative to differences among tissues and species (Additional file 1: Figure S1).

A total of 9,258 and 9,385 annotated genes were detected in the antennae of An. gambiae and An. quadriannulatus respectively. Of these, 2,593 (28.0%) are significantly higher expressed (q value of < 0.05) in antennae of An. gambiae and 2,778 (29.6%) in the antennae of An. quadriannulatus (Figure 1A). In the maxillary palps, 9,824 and 9,994 genes are expressed in An. gambiae and An. quadriannulatus, respectively. Of these, 1,243 (12.6%) genes are significantly up-regulated in An. gambiae and 1,517 (15.2%) in An. quadriannulatus respectively (Figure 1B).

Figure 1
figure 1

Differential gene expression between An. gambiae and An. quadriannulatus in antennae (A) and maxillary palps (B). The expression of genes indicated in red is statistically significant (q < 0.05).

A gene ontology analysis (GO) was conducted to recover descriptions of molecular and biological function. For this analysis only significantly enhanced genes that were more than 2-fold expressed were considered. This resulted in 564 antennal genes in An. gambiae and 608 antennal genes in An. quadriannulatus (Figure 2, Additional file 2: Figure S2). For the maxillary palps, 870 and 787 genes met these criteria in the two respective species (Additional file 3: Figure S3). Of these, 217 genes are shared between the antenna and palps of An. gambiae, and 292 of these genes are shared between tissues in An. quadriannulatus (Additional file 4: Figure S4). Not surprisingly, some of the gene ontology (GO) terms recovered in the significantly enhanced genes are connected to olfaction (e.g., “odorant binding”) and signal transduction (e.g., “response to stimulus”, “signal transducer activity”). Additionally, we found strong representation of terms connected to enzymatic activity. For example, “transferase activity” represents 7% of the up-regulated genes in the antennae of An. gambiae and 5% of those in the maxillary palps of this species (Figure 2, Additional file 3: Figure S3).

Figure 2
figure 2

GO analysis of 564 An. gambiae genes with > 2-fold antennal expression compared to An. quadriannulatus , predicting their involvement in molecular functions (A) and biological processes (B). Data are presented as level 3 GO categorization. Categories with less than 1% of representation were grouped in “others”.

Olfactory receptors

Out of the 76 annotated olfactory receptors (Ors), 65 were detected above the threshold in the antennae of at least one species (Additional file 5: Table S1). As expected, the co-receptor Orco is highly expressed in both An. gambiae and An. quadriannulatus female antennae (1,429 and 1,756 RPKM, respectively), but is significantly higher in An. quadriannulatus (q = 0.025) (Additional file 5: Table S1). Consistent with this observation, the total expression of the specific Ors is higher in this species as well (1,738 vs 2,183 RPKM), which is also reflected in the regression slope (1.20) for Or expression between species.

A total of 17 Ors are expressed at a significantly higher level in An. gambiae female antennae (Table 1, Figure 3]. The expression level of Ors with significantly enhanced expression in An. gambiae ranged from 2.5 to 42.8 RPKM, but Or36, 45, 66, 69, 70, 73 and 75 are noteworthy for being both relatively abundant (>12.3 RPKM) and substantially up-regulated (>1.9-fold) in An. gambiae. Two expressed Ors (Or8 and 51) were not expressed in the antennae of An. quadriannulatus, but these were among the least abundant Ors in An. gambiae as well (2.5 and 2.7 RPKM, respectively).

Table 1 Olfactory and gustatory genes that are significantly enhanced in the female antennae of An. gambiae vs An. quadriannulatus
Figure 3
figure 3

Regression plot of gene expression between the antennae of An. gambiae and An. quadriannulatus for odorant receptors (A), ionotropic receptors (B), odorant binding proteins (C) and gustatory receptors (D). Axis represent Ln(RPKM) values. Inset box shows regression line based on non-transformed data. Genes whose expression was significantly different are indicated in black.

In contrast to An. gambiae, twenty-eight specific Ors are expressed at a significantly higher level in the female antennae of An. quadriannulatus. In this species Or1, 9, 23, 33, 46, 61, and 63 stand out by being both highly expressed (RPKM >17.3) and substantially enhanced (>2.0-fold). Four Ors (Or18, 20, 30 and 74) are expressed only in the antennae of An. quadriannulatus, although at low levels (1.1 < RPKM < 6.42).

Although no abundant Ors are uniquely expressed in the antennae of either species, our analysis identifies a set of Ors that show clear species-specific enhancement of their antennal expression. Despite these specific differences, a linear regression analysis shows that overall antennal Or expression is highly correlated between the two species with (R2 = 0.937, slope = 1.20, Figure 3A).

Strikingly, the overall expression of Ors is much higher in the palps of An. quadriannulatus compared to An. gambiae. For example, Orco is expressed at 187.3 and 700.1 RPKM in An. gambiae and An. quadriannulatus respectively, and the regression slope for Or expression between the two species is 2.84 (Figure 4A). This 2.8 to 3.7-fold enhancement of olfactory receptors expression implies a relatively larger importance of the maxillary palps in the olfaction system of An. quadriannulatus.

Figure 4
figure 4

Regression plot of gene expression between the maxillary palps of An. gambiae and An. quadriannulatus for odorant receptors (A), ionotropic receptors (B), odorant binding proteins (C) and gustatory receptors (D). Axis represent Ln(RPKM) values. Inset box shows regression line based on non-transformed data. Genes whose expression was significantly different are indicated in black.

That being said, Ors are expressed at much lower level in the palps than in the antennae for both An. gambiae (slope = 0.105, Figure 5A) and An. quadriannulatus (slope = 0.236, Figure 6A). The number of detected Ors is also substantially less in the maxillary palps. Only 45 specific Ors were detected in the palps of An. gambiae, whereas 53 are present in the An. quadriannulatus female palps (Additional file 6: Table S2). Interestingly, one olfactory receptor, Or52, is unique to the maxillary palps of An. gambiae. This gene is among the seven most abundant Ors in this species (9.55 RPKM), but did not reach our detection threshold in An. quadriannulatus (0.83 RPKM). Furthermore, this gene is all but undetectable in the antennae of either species (0.13 and 0.18 RPKM, Additional file 5: Table S1).

Figure 5
figure 5

Regression plot of gene expression between the antennae and maxillary palps of An. gambiae and for odorant receptors (A), ionotropic receptors (B), odorant binding proteins (C) and gustatory receptors (D). Axis represent Ln(RPKM) values. Inset box shows regression line based on non-transformed data. Genes whose expression was significantly different are indicated in black.

Figure 6
figure 6

Regression plot of gene expression between the antennae and maxillary palps of An. quadriannulatus and for odorant receptors (A), ionotropic receptors (B), odorant binding proteins (C) and gustatory receptors (D). Axis represent Ln(RPKM) values. Inset box shows regression line based on non-transformed data. Genes whose expression was significantly different are indicated in black.

Or expression in the antennae and palps is highly correlated in An. gambiae (R2 = 0.80, Figure 5A), but considerably less so in An. quadriannulatus (R2 = 0.57, Figure 6A). Of the 17 Ors significantly enhanced in An. gambiae antennae, most are enhanced in the palps of this species compared to An. quadriannulatus as well. However, there are two notable exceptions; Or8, which shows 3.8-fold up-regulation in the An. gambiae antennae, is expressed 4.9-fold higher in the palps of An. quadriannulatus. It is also one of most abundant Ors in the palps of this latter species. Similarly, Or28 is significantly 2.0-fold enhanced in An. gambiae antennae, but 6.1-fold enhanced in the palps of An. quadriannulatus where it is the second most highly expressed specific Or.

Not surprisingly, given the 2.8-fold higher overall level of Or expression, 20 specific Ors are significantly enhanced in the palps of An. quadriannulatus vs. An. gambiae. With the exception of Or8 and Or28, these are also enhanced in the antennae of this species, or are not expressed in either (Or3 and 5). Overall, the correlation between Or expression in the palps (R2 = 0.80, Figure 4A) of the two species is less than for the antennae (R2 = 0.94, Figure 3A).

Ionotropic receptors

The total antennal expression of Irs is similar between species (1,283 vs 1,232 RPKM in An. gambiae and An. quadriannulatus respectively), with a regression slope of 1.07 (Figure 3B). Of the 44 annotated Irs, 29 are expressed in the female antennae of An. gambiae and 32 are expressed in An. quadriannulatus (Additional file 5: Table S1). Ir25a, one of the co-receptors, is by far the most highly expressed Ir in both species. However, both species also express Ir76b at very high levels. This gene has been considered a putative co-receptor [26], but which was more recently proposed to encode a Na+ leak channel which in Drosophila also plays a role in salt detection [41]. A total of 15 Irs are significantly up-regulated in An. gambiae antennae, with nine (Ir75h.2, 7t, 7w, 41c, 41n, 75g, 75k, 93a, and 100a) standing out by being both considerably enhanced (>1.8-fold), as well as among the more abundant Irs (Figure 3A, Additional file 5: Table S1). Twelve Irs are significantly enhanced in the female antennae of An. quadriannulatus. Only two specific ionotropic receptors, Ir75d and 75h.1, are considerably up-regulated (>1.9-fold) and abundant (>19.5 RPKM) in this species. Similarly to Or expression, Ir expression is highly correlated between species (R2 = 0.79, Figure 3B).

The total Ir expression is much lower in the palps than the antennae (slope = 0.115 for An. gambiae, Figure 5B, and 0.384 for An. quadriannulatus, Figure 6B), but like the Ors is much higher in An. quadriannulatus (regression slope = 3.7, Figure 4B). Twenty-four Irs are expressed in the palps of An. gambiae, of which only Ir101 is significantly 3.1-fold enhanced in this species, but it is expressed at very low levels in the palps (and antennae) of both species (Table 2, Additional file 6: Table S2). None of the 15 significantly enhanced Irs in An. gambiae antennae are significantly enhanced in its palps. The expression of eight Irs is significantly enhanced in the palps of An. quadriannulatus. Six of these are also significantly enhanced in the antennae of this species. Ir expression in the palps is highly correlated between the two species (R2 = 0.89, Figure 4B), as well as between the antennae and palps of the two species, with R2 = 0.92 for An. gambiae and R2 = 0.98 An. quadriannulatus (Figures 5B and 6B) (See Table 3).

Table 2 Olfactory and gustatory genes that are significantly enhanced in the female
Table 3 Summary of olfaction/gustation genes whose expression is significantly different between An. gambiae and An. quadriannulatus , and those who are also > 2-fold expressed

Odorant binding proteins

As expected based on OBP function, Obp expression in the antennae is considerably higher than that of the Ors and Irs (as much as 51,541.1 RPKM in An. gambiae and 61,872.7 RPKM in An. quadriannulatus, Figures 3C and 4C). In fact, Obp48 is by far the most highly expressed gene in the antennae of both species, and nine of the top 15 most highly expressed genes are Obps (Additional file 5: Table S1). That being said, only 27 and 29 of the 57 putative Obps were detected in An. gambiae and An. quadriannulatus female antennae.

The overall level of antennal Obp expression is similar in An. gambiae and An. quadriannulatus (183,218 vs 172,629 RPKM, slope = 1.08, Figure 3C). Expression of eleven Obps is significantly enhanced in An. gambiae female antennae, and three abundant odorant binding proteins, Obp1, 3 and 10 are more than 1.8-fold enhanced in this species (Table 1). The expression of five Obps was significantly higher in An. quadriannulatus, but of the three that were considerably up-regulated (>1.9-fold), only Obp19 was expressed at any appreciable level (464.8 RPKM). Similarly to Or and Ir expression, antennal Obp expression was highly correlated between species (R2 = 0.95, Figure 3C).

Also consistent with Or and Ir expression, the Obp abundance in the palps is considerably lower than in the antennae (slope = 0.133 for An. gambiae and 0.411 for An. quadriannulatus, Figures 5C and 6C). Again similar to Or and Ir expression, the overall Obp expression is several folds higher in the palps of An. quadriannulatus as compared to An. gambiae (slope = 3.39, Figure 4C). Furthermore, Obp48, the most abundant gene in the antennae of both species, ranks 8th in An. gambiae maxillary palps in abundance, whereas it is also the most highly expressed of all genes in An. quadriannulatus.

One relatively abundant odorant binding protein, Obp26 (337.2 RPKM), is expressed at significantly higher levels in An. gambiae palps (2.1-fold), and is also significantly up-regulated in the antennae of An. gambiae. Not surprisingly given the 3.4-fold higher level of overall Obp expression in An. quadriannulatus palps, 14 Obps are significantly enhanced in this species, including Obp48 (Figure 4C, Additional file 6: Table S2). However, the correlation between Obp expression in the palps of the two species is high (R2 = 0.95, Figure 4C), and there is also a strong correlation between Obp expression in the antennae and palps for both species (R2 = 0.90 for An. gambiae and R2 = 0.88 for An. quadrianulatus, Figures 5C and 6C).

Gustatory receptors

The gustatory receptors (AgGrs) are expressed at very low levels in the antennae of both species, with total RPKM values of 57 and 71 in An. gambiae and An. quadriannulatus respectively (Figures 3D and 4D, Additional file 5: Table S1). None-the-less, the expression of seven AgGrs is significantly enhanced in An. gambiae, and the same number is significantly up-regulated in An. quadriannulatus. Interestingly, five AgGrs that are significantly enhanced in An. gambiae (AgGr48-52) are not expressed in An. quadriannulatus, although the expression of these genes in An. gambiae is low as well, ranging from 1.8 to 7.8 RPKM. To compare, ranking these with the Ors in level of abundance would place AgGr52 in 44th position. Two AgGrs, AgGr1 and AgGr26, stand out in the An. quadriannulatus antennal dataset. They are relatively abundant in this species (8.7 and 21.3 RPKM) and significantly up-regulated (>2.6-fold). In contrast to the olfaction gene families, little correlation exists between AgGr expression in the antennae of the two species (R2 = 0.01, Figure 3D).

Only 18 out of 60 annotated AgGrs are expressed in the palps of An. gambiae, and 15 are expressed in An. quadriannulatus palps (Figure 4D, Additional file 6: Table S2). The three AgGrs responsible for CO2 detection in mosquito palps (AgGr 22, 23 and 24) are by far the most highly expressed AgGrs in An. quadriannulatus. Interestingly, the expression level of these CO2 receptor genes is between 5.1 and 6.1-fold higher in this species, but is on par with the expression of other AgGrs, such as AgGr48 and 49, in An. gambiae. Overall the regression slope between the two species for the palps is 3.03 (Figure 4D), indicating a much higher expression in An. quadriannulatus, but this is entirely due to the CO2 receptors.

Probably the most striking result from the palp data set is a set of species-specific AgGrs in An. gambiae. The five up-regulated gustatory receptors from the An. gambiae antennae (AgGr48-52) are in fact highly abundant in the palps of this species (11.1 < RPKM < 24.2, Table 2). Interestingly, they are all but absent from the palps of An. quadriannulatus (<0.84 RPKM).

Similarly, four AgGrs (15, 17, 21, and 26) are abundant and up-regulated in An. quadriannulatus palps (9.5 < RPKM < 24.7, > 4.9-fold). These genes are all at very low levels in An. gambiae palps, with AgGr26 being unique to An. quadriannulatus. Interestingly, this gene is also expressed at high levels in the antennae of An. quadriannulatus (21.3 RPKM).

Discussion

Because of An. gambiae’s odor-mediated host seeking behavior [1], it is expected that its preference for human hosts has a strong genetic basis in its olfactory system. This system is primarily housed in the antennae, but the maxillary palps are also involved [14]. In this study a comparison of the olfactory organ transcriptomes of the anthropophilic An. gambiae and its zoophilic sibling species An. quadriannulatus identified species-specific patterns of olfaction gene expression. Even though the expression profiles of olfaction genes are highly correlated between species, clear differences were observed which identify olfaction genes that may play an important role in differential host preference.

The olfactory system of Anopheles mosquitoes plays a role in at least two other aspects of their biology; finding a sugar source, most often nectar, and identifying oviposition sites. No data is available on how often or from what source An. quadriannulatus females obtain sugar meals. However, An. gambiae starts ignoring honey volatiles five days after emergence and responds almost exclusively to human odor at that point [42]. Our experience in the laboratory indicates that An. quadriannulatus also switches to host seeking around this time. The larval ecology of both species appears to be similar, with both breeding in shallow, open, sunlit fresh water pools [43, 44], and in any case oviposition-site searching does not commence until 48 hours post-blood feeding. Therefore, although we cannot rule out that the differences in olfaction gene expression between the two species are due to biological differences other than host-seeking, there is no data to suggest that such differences are substantial. In addition, the use of 6-day old females in our study optimizes our ability to detect differences in olfaction gene expression that are related to host-seeking [45].

In the antennal transcriptome, seven Ors (Or36, 45, 66, 69, 70, 73, 75) and nine Irs (Ir75h.2, 7t, 7w, 41c, 41n, 75g, 75k, 93a, 100a) stand out by being among the more highly expressed receptor genes, while also being considerably up-regulated in An. gambiae (1.8 to 4.7-fold). We speculate that the enhanced expression of some of these genes in An. gambiae contributes to an increased sensitivity to human odor. Divergence in olfaction gene expression associated with host specialization has been observed between closely related Drosophila species that feed on different host plants. Antennal expression of Ors differed markedly between the generalists D. melanogaster, D. simulans, and their specialist sister-species D. sechellia, and as many as 53% of the Ors were differentially expressed between species pairs [32].

Previous studies have examined the response of 56 AgOrs to a wide range of odorants [24, 25]. These included 11 of the 17 Ors enhanced in An. gambiae in the present study. Only four receptors showed a positive response to any of the tested odorants. The exposure of Or75 to eight human volatiles led to a small to moderate increase in the firing rate of the neuron [24]. Or36 has a very narrow tuning curve and responded strongly to only two of the tested odorants, however neither of which are of human origin [25]. Or65 responded mildly to one human odorant (4 methylphenol), as well as to several other chemicals [24]. Finally, Or8 on the other hand responded strongly to two known human odorants, 1-octen-3-ol and 1-hepten-3-ol [24, 25]. However, for many of these human volatiles it is not known if they are unique to humans. For example 1-octen-3-ol is exhaled by bovids as well, and is a common compound produced by mushrooms [46]. Finally, it should be kept in mind that 346 volatiles have been identified in human sweat [7], and only a small subset of these volatiles were tested on these Ors.

Expression differences in Obps may also play a role in the human host preference of An. gambiae. Several highly expressed Obps (1, 3, 10) are enhanced (1.8 to 2.4-fold) in the antennae of An. gambiae, whereas no abundant Obps are substantially higher expressed in An. quadriannulatus. The presence of specific odorant-binding proteins is well-known to impact behavior in Drosophila. Flies carrying LUSH, a mutant OBP, are defective in detecting an aggregation pheromone [47], and Obp57d and Obp57e are involved in differences in oviposition behavior between Drosophila species [48].

Importantly, our data shed new light on the role of the maxillary palps in odor detection in these species. Lu et al. [14] concluded that a relatively small repertoire of Ors is responsible for olfaction coding in the maxillary palps, although it was found that in Culex quinquefasciatus the maxillary palps are broad spectrum odorant detectors [49]. A previous analysis of the transcriptome of An. gambiae maxillary palps indicated the expression of relatively small Or repertoire [20], with only four Ors expressed at >1 RPKM. Those data contrast with our results, in which 49 Ors were detected in the palps of this species. Although it is not clear at what level the expression of receptors is biologically relevant, 19 of these Ors are expressed at > 4 RPKM, suggesting the palps may be able to detect a suite of odors. Possible reasons for these contrasting results between the two studies may be that we conducted our dissections during the early dark cycle, and included replicates, which is recommended for obtaining reliable RNAseq data [50].

Another interesting feature of our data is the several fold higher overall expression of the three olfaction gene families in the maxillary palps of An. quadriannulatus compared to An. gambiae. Clearly, the maxillary palps of An. quadriannulatus are considerably more important component of this species’ olfactory olfaction system than is the case for An. gambiae.

While no An. gambiae specific olfaction genes were identified in the antennae, an analysis of the maxillary palp transcriptome revealed several An. gambiae specific chemosensory genes. Or52, the seventh most abundant Or in the An. gambiae palps, is not expressed in An. quadriannulatus. Interestingly, this gene is also absent from the antennae of An. gambiae. This indicates that this maxillary palp receptor could play a species-specific role in An. gambiae’s biology, and thus may possibly be involved in human host preference. Unfortunately, this Or was not including in the odorant affinity studies discussed above [24, 25].

Additionally, the expression pattern of several AgGrs indicates that they play a species-specific role in olfaction. AgGr48-52 are specific to An. gambiae, and these five AgGrs are expressed at relatively high levels in the maxillary palps, indicating a functional role of the receptors encoded by these genes. With the exception of AgGr22-24, which together encode the heteromeric CO2 receptor [14, 51], gustatory receptors are generally considered to be primarily involved in gustation. However, the fact that AgGr48-52 are expressed in both the antennae and maxillary palps of An. gambiae suggests a role in olfaction for these genes. These genes are located in tandem on the chromosome 2R and each pair is separated by only 46 to 326 bp. Therefore, the expression of these genes is likely controlled by the same regulatory elements. Additionally, AgGr26 is specific to An. quadriannulatus. It is expressed at high levels only in the maxillary palps and the antennae of this species, suggesting a species-specific role. Therefore, these gustatory receptors may play a significant role in the behavioral differentiation between An. gambiae and An. quadriannulatus.

The several fold higher expression of olfaction genes in the palps of An. quadriannulatus dominates the comparison between the palps of the two species. Nevertheless, Obp26 and Obp56 are more than 2-fold enhanced in the maxillary palps of An. gambiae. Interestingly, both are enhanced in the antennae of this species as well, although Obp56 is expressed at very low levels in this organ.

Several other noteworthy observations result from these data. The expression of Orco and the Ors is 18% to 26% higher in the antennae of An. quadriannulatus. Although the antennal sensilla of An. quadriannulatus were found to outnumber those of An. gambiae, their density is actually similar in both species [52], indicating that this difference is not explained by differences in the antennal morphology. Furthermore, the overall level of Ir and Obp expression is actually slightly higher in An. gambiae antennae. This implies that the overall sensitivity to the odorants detected by Ors is higher in An. quadriannulatus antennae.

The expression of the CO2 receptor genes AgGr22-24 in the maxillary palps of both species differs markedly. These genes are expressed between 5.1-6.1 fold higher in An. quadriannulatus, which greatly exceeds the overall higher level of expression of olfaction genes in the palps of this species. By itself, CO2 is a poor attractant to An. gambiae. It does activate and guide it towards a human odor source, but at this point other semiochemicals become important [53]. In contrast, CO2 is highly attractive to An. quadriannulatus, which has a more catholic host preference [38]. Our data suggest that the lesser attraction of An. gambiae to CO2 is accompanied by a lower sensitivity to CO2. It has been suggested that anthropophilic mosquitoes primarily use CO2 to detect hosts as a long distance cue [54]. Given the smaller amount of CO2 produced by a human vs. the preferred host of An. quadriannulatus; a bovid, and the relatively low expression level of the CO2 receptor genes in An. gambiae, this species either relies little on CO2 for its long range attraction, or is incapable of detecting hosts from the same distances as An. quadriannulatus.

Although the An. quadriannulatus strain examined in this study showed a preference for bovine hosts, this species did not distinguish between human and cow sweat in an olfactometer in a previous study [4]. Furthermore, when offered a choice of a human or equal sized calf, it blood fed equally on both [3]. This suggests that An. quadriannulatus has a wider host preference and is more of a generalist than An. gambiae. This is consistent with its much higher level of expression of the CO2 receptor genes. However, the large number of olfaction genes with enhanced expression in An. quadriannulatus indicates that it, like An. gambiae, likely responds to a complex blend during host seeking, similar to An. gambiae[12, 55].

Rinker et al. [36] recently compared the daytime transcriptomes of the antennae of An. gambiae and An. quadriannulatus. Their results correspond roughly with ours. A linear regression analyses of the day-time olfaction gene expression with the early dark phase reported here, resulted in R2 values of 0.74-0.78 for the three olfaction gene families in An. gambiae, and 0.76-0.91 in An. quadriannulatus. That being said, a few notable exceptions were observed. For example, the expression differences for Or66, 73 and Obp26 was much less pronounced during the daytime. For a few genes, e.g. Ir75k, 7t and 75g the expression pattern was even reversed. Similarly, several olfaction genes showed differential expression during the day-time, but not during dark cycle (e.g. Obp2 and 13). The expression of Orco and a variety female antennae Obps in female Anopheles gambiae fluctuates throughout the circadian cycle. Expression of these olfaction genes was generally found to be highest during the early stage of the dark phase, and therefore seems to be correlated with the female’s host seeking activity [37]. However, a comparison between light and dark cycle transcriptomes, suggests that the expression pattern of olfaction genes do not all follow the same diel expression cycle.

Ionotropic receptors have been divided into “antennal” and “divergent” IRs depending on whether they are expressed in the antennae of Drosophila[16]. It was suggested that this distinction held across a wide range of insects and that divergent IRs play a predominant role in gustation rather than olfaction. Additionally, it was found that antennal IRs tend to be more conserved than the divergent IRs. However, this classification has limited relevance to the expression pattern observed in the antennae of An. gambiae. Of the expressed Irs, 17 were classified as antennal, and 14 as divergent IRs, although of the 14 Irs not expressed in the antennae, 12 are divergent IRs.

Conclusion

Our data identifies potential human host preference genes in the malaria vector An. gambiae, but also provides new insight into the importance of the maxillary palps in the olfactory system. The palps are where the most dramatic difference in chemosensory gene expression is observed between the anthropophilic An. gambiae and the zoophilic An. quadriannulatus, with several highly expressed receptor genes that are specific to either species. Finally, the expression patterns of several AgGrs strongly suggest a species-specific role for them in the olfaction system of An. gambiae.

Ethics statement

Colonies of Anopheles mosquitoes were kept following the Arthropod Containment Guidelines established by The American Committee of Medical Entomology of the American Society of Tropical Medicine and Hygiene. The behavioral experiments were conducted in the Laboratory of Entomology at Wageningen University in the Netherlands. Approval to obtain an odor sample from a cow was obtained from the Animal Use Committee of Wageningen University.

Methods

Mosquito rearing

Laboratory strains of An. gambiae M form (GASUA), recently proposed to be named An. coluzzii[56], and originally collected in Suakoko, Liberia, as well as An. quadriannulatus (SANQUA) established from female mosquitoes collected in Sangwe, Zimbabwe were reared in the insectaries at Wageningen University, The Netherlands (host choice experiment) and Texas A&M University, College Station, TX, USA (RNAseq analyses). Rearing conditions were 25°C, 75-85% relative humidity and a light:dark photoperiod of 12 hours. Female mosquitoes were blood fed on defibrinated rabbit blood using a membrane feeding system. Larvae were maintained at densities of approx. 150 per 2 L container and fed finely ground fish food (Tetramin, Melle, Germany). Pupae were collected and placed into cages at densities of two cups of 150 pupae per cage.

For mosquitoes used in RNAseq analyses, cages were checked daily for newly emerged mosquitoes. To ensure mosquitoes were the same age, pupae that did not eclose were transferred to new cages. Male and females mosquitoes were kept together in a cage and fed a 5-10% sucrose solution for six days until tissue dissections. Hence, females were given an opportunity to mate, but not to blood-feed. We checked the insemination rate in 50 An. gambiae females at day 6 and found it to be high (82%).

Dual odor-choice assay

A total of 750 female An. gambiae and 330 An. quadriannulatus females were tested to determine the odor preference of the two species in laboratory conditions. Mosquitoes were put at a density of 75–80 in release cages for use in a dual-choice olfactometer [57] the night before experiments, and provided with a wet cotton ball for hydration. Human and cow odor traps were prepared on the morning of the experiments. Human odor was derived from the socks worn by volunteers for 24 hours, and cow odor was derived from a panty hose tied around the leg of a cow for 24 hours. Odor sources were switched between the left and right port of the olfactometer between runs. A single, centrally placed CO2 plume was used as activator. Conditions during the experiments were as follows: temperature = 26-28°C, humidity 55-75% inside olfactometer, 80% in front of port holes, air-speed 018–0.22 ms−1, and released [CO2] = 4.5%. Mosquitoes were released into the olfactometer during the dark-cycle for 15 min. under semi-dark conditions. Mosquitoes remaining in the wind tunnel after the experiment were disposed of.

Molecular methods

Female mosquitoes were killed shortly after the start of the dark cycle by placing them at −20°C. This is when anophelines begin their host searching activity and when Orco expression in An. gambiae peaks during the circadian rhythm [37]. The antennae and maxillary palps were removed from frozen mosquitoes placed on dry ice and were stored in RNAlater® (Ambion). Between 600 to 800 6-day old females were dissected for each replicate and three replicates per species were included for a total of six samples per species. Samples were stored at 4°C for 24 hours, before RNAlater was removed and stored at −80°C until RNA extraction.

Total RNA was isolated from each sample using miRNeasy (Qiagen) columns according to the protocol supplied by Qiagen. RNA quantity was initially verified using a Qubit fluorometer (Life Technologies). Next, RNA was further quantified using a NanoDrop spectrophotometer (Thermo Scientific) and the quality assessed using RNA Pico LabChip analysis on an Agilent BioAnalyzer 2100 (Agilent Technologies) by the Agrilife Genomics Center at Texas A&M University.

mRNA was isolated from 1 μg of total RNA and cDNA libraries were prepared using an Illumina TruSeq RNA Library kit (Illumina). Each single-end library contained two/three replicates that had been given a unique tag using barcode sequences supplied by the library kit. Each library was sequenced on a single lane of an Illumina flow cell and using 50 cycles on an Illumina HiSeq 2000. Preparation and sequencing of libraries were both performed by the University at Buffalo Next-Generation Sequencing and Expression Analysis Core Facility. Approximately 50–70 million reads with an average read of 51 base pairs were generated for each replicate sample and used for further analysis.

RNA sequencing analyses

Read quality was assessed using FastQC (ver 0.10.0) and processed using NGS QC toolkit [58] with at least 80% of the reads had Phred > 30 (raw reads Phred quality score 0–40). Reads were trimmed and then filtered by length, discarding reads < 40 bp. Sequencing reads were mapped to the reference An. gambiae genome (AgamP3; December 2013) using the software package STAR [59]. Alignments were discarded if they had more than two mismatches. Read counts were conducted with HTSeq-count (ver 0.5.4) (http://www.huber.embl.de/users/anders/HTSeq/doc/count.html). Only reads that aligned to a unique location in the genome were used to calculate the expression levels. Sequence data was obtained for three replicates of the antennae for each species, and for three replicates for the palps. One palp replicate for each species provided poor quality data, and these were therefore discarded from further analyses. Tests for differential expression in the female antennae or palps from An. gambiae versus An. quadriannulatus were performed in the R package DESeq2 [60]. Size factors for each dataset were calculated to normalize library sizes across replicates, and overall means and variances were determined using a negative binomial distribution model. Genes were considered to be differentially expressed if q < 0.05 after correcting for multiple testing. Genes were considered not expressed if RPKM <1.

To compare the tissue and species effect on the overall gene expression, we computed the correlation coefficient (R2) and slope from a linear regression between An. gambiae versus An. quadriannulatus data sets, as well as between maxillary palps and antennae. Scripts used to run RNAseq analyses are presented in Additional file 7. Reads for Obp6 and Obp29 mapped to multiple locations in the genome, therefore expression data for these two genes are unreliable, and not considered in this analysis.

Gene ontology analysis

Genes that met the following criteria: q < 0.05, FoldChange > 2, RPKM > 1 at least in one sample, between antenna and maxillary palps of An. gambiae versus An. quadriannulatus were used for gene ontology (GO) analyses. GO Annotation was performed using Blast2GO [61]. The gene sequences were retrieved from Ensembl Genomes release 22 via Biomart (http://metazoa.ensembl.org/index.html). GO annotation was used for assessment of the genes differentially expressed in each sample. GO annotation associates analyzed transcripts with terms from hierarchical vocabularies describing, e.g., molecular function or biological process.

Availability of supporting data

All fastq files containing the raw data were deposited at the NCBI Sequence Read Archive. [http://www.ncbi.nlm.nih.gov/sra?term=SRP050131]. The full gene expression data are available in “Additional file 5” (antennae) and “Additional file 6” (maxillary palps).

References

  1. Takken W, Knols BG: Odor-mediated behavior of Afrotropical malaria mosquitoes. Annu Rev Entomol. 1999, 44: 131-157. 10.1146/annurev.ento.44.1.131.

    Article  CAS  PubMed  Google Scholar 

  2. Takken W, Eling W, Hooghof J, Dekker T, Hunt R, Coetzee M: Susceptibility of Anopheles quadriannulatus Theobald (Diptera: Culicidae) to Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1999, 93: 578-580. 10.1016/S0035-9203(99)90054-8.

    Article  CAS  PubMed  Google Scholar 

  3. Pates HV, Takken W, Curtis CF, Huisman PW, Akinpelu O, Gill GS: Unexpected anthropophagic behaviour in Anopheles quadriannulatus. Med Vet Entomol. 2001, 15 (3): 293-298. 10.1046/j.0269-283x.2001.00310.x.

    Article  CAS  PubMed  Google Scholar 

  4. Pates HV, Takken W, Curtis CF: Laboratory studies on the olfactory behavior of An. quadriannulatus. Entomol Exp Appl. 2005, 114: 153-159. 10.1111/j.1570-7458.2005.00249.x.

    Article  Google Scholar 

  5. Braks MAH, Takken W: Incubated human sweat but not fresh sweat attracts the malaria mosquito Anopheles gambiae sensu stricto. J Chem Ecol. 1999, 25 (3): 663-672. 10.1023/A:1020970307748.

    Article  CAS  Google Scholar 

  6. Verhulst NO, Andriessen R, Groenhagen U, Bukovinszkiné Kiss G, Schulz S, Takken W, van Loon JJA, Schraa G, Smallegange RC: Differential attraction of malaria mosquitoes to volatile blends produced by human skin bacteria. PLoS One. 2010, 5 (12): e15829-10.1371/journal.pone.0015829.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Bernier UR, Kline DL, Barnard DR, Schreck CE, Yost RA: Analysis of human skin emanations by gas chromatography/ mass spectrometry. 2. Identification of volatile compounds that are candidate attractants for the yellow fever mosquito (Aedes aegypti). Anal Chem. 2000, 72: 747-756. 10.1021/ac990963k.

    Article  CAS  PubMed  Google Scholar 

  8. Braks MAH, Meijerink J, Takken W: The response of the malaria mosquito, Anopheles gambiae to two components of human sweat, ammonia and L-lactic acid. Physiol Entomol. 2001, 26: 142-148. 10.1046/j.1365-3032.2001.00227.x.

    Article  CAS  Google Scholar 

  9. Dekker T, Steib B, Carde RT, Geier M: L-lactic acid: a human-signifying host cue for the anthropophilic mosquito Anopheles gambiae. Med Vet Entomol. 2002, 16: 91-98. 10.1046/j.0269-283x.2002.00345.x.

    Article  CAS  PubMed  Google Scholar 

  10. Smallegange RC, Qiu YT, van Loon JJ, Takken W: Synergism between ammonia, lactic acid and carboxylic acids as kairomones in the host-seeking behaviour of the malaria mosquito Anopheles gambiae sensu stricto (Diptera: Culicidae). Chem Senses. 2005, 30: 145-152. 10.1093/chemse/bji010.

    Article  CAS  PubMed  Google Scholar 

  11. Van den Broek IVF, den Otter CJ: Olfactory sensitivities of mosquitoes with different host preferences (Anopheles gambiae s.s., An. arabiensis, An. quadriannulatus, An. m. atroparvus) to synthetic odours. J Insect Physiol. 1999, 45: 1001-1010. 10.1016/S0022-1910(99)00081-5.

    Article  CAS  PubMed  Google Scholar 

  12. Mukabana WR, Mweresa CK, Otieno B, Omusula P, Smallegange RC, Loon JJA, Takken W: A novel synthetic odorant blend for trapping of malaria and other African mosquito species. J Chem Ecol. 2012, 38 (3): 235-244. 10.1007/s10886-012-0088-8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Pitts RJ, Fox AN, Zwiebel LJ: A highly conserved candidate chemoreceptor expressed in both olfactory and gustatory tissues in the malaria vector, Anopheles gambiae. Proc Natl Acad Sci U S A. 2004, 101: 5058-5063. 10.1073/pnas.0308146101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Lu T, Qiu YT, Wang G, Kwon JY, Rutzler M, Kwon HW, Pitts RJ, van Loon JJA, Takken W, Carlson JR, Zwiebel LJ: Odor coding in the maxillary palp of the malaria vector mosquito Anopheles gambiae. Curr Biol. 2007, 17: 1533-1544. 10.1016/j.cub.2007.07.062.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Zwiebel LJ, Takken W: Olfactory regulation of mosquito-host interactions. Insect Biochem Mol Biol. 2004, 34: 645-652. 10.1016/j.ibmb.2004.03.017.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, Toby J, Gibson TJ, Benton R: Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010, 6 (8): e1001064-10.1371/journal.pgen.1001064.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Li Z-X, Pickett JA, Field LM, Zhou JJ: Identification and expression of odorant binding proteins of the malaria-carrying mosquito Anopheles gambiae and Anopheles arabiensis. Archiv Insect Biochem Physiol. 2005, 58: 175-189. 10.1002/arch.20047.

    Article  CAS  Google Scholar 

  18. Siju KP, Reifenrath A, Scheiblich H, Neupert S, Predel R, Hansson BS, Schachtner J, Ignell R: Neuropeptides in the antennal lobe of the yellow fever mosquito Aedes aegypti. J Comp Neurol. 2014, 522 (3): 592-608. 10.1002/cne.23434.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Hill CA: G protein-coupled receptors in Anopheles gambiae. Science. 2002, 298 (5591): 176-178. 10.1126/science.1076196.

    Article  CAS  PubMed  Google Scholar 

  20. Pitts RJ, Rinker DC, Jones PL, Rokas A, Zwiebel LJ: Transcriptome profiling of chemosensory appendages in the malaria vector Anopheles gambiae reveals tissue- and sex-specific signatures of odor coding. BMC Genomics. 2011, 12 (1): 271-10.1186/1471-2164-12-271.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Liu C, Pitts RJ, Bohbot JD, Jones PL, Wang G, Zwiebel LJ: Distinct olfactory signaling mechanisms in the malaria vector mosquito Anopheles gambiae. PLoS Biol. 2010, 8 (8): e1000467-10.1371/journal.pbio.1000467.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Hallem E, Fox AN, Zwiebel LJ, Carlson JR: A mosquito odorant receptor tuned to a component of human sweat. Nature. 2004, 427: 212-213. 10.1038/427212a.

    Article  CAS  PubMed  Google Scholar 

  23. Xia Y, Wang G, Buscariollo D, Pitts RJ, Wenger H, Zwiebel LJ: The molecular and cellular basis of olfactory-driven behavior in Anopheles gambiae larvae. Proc Natl Acad Sci U S A. 2008, 105: 6433-6438. 10.1073/pnas.0801007105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Carey AF, Wang G, Su C-Y, Zwiebel LJ, Carlson JR: Odorant reception in the malaria mosquito Anopheles gambiae. Nature. 2010, 464 (7285): 66-71. 10.1038/nature08834.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Wang G, Carey AF, Carlson JR, Zwiebel LJ: Molecular basis of odor coding in the malaria vector mosquito Anopheles gambiae. Proc Natl Acad Sci U S A. 2010, 107 (9): 4418-4423. 10.1073/pnas.0913392107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Abuin L, Bargeton B, Ulbrich MH, Isacoff EY, Kellenberger S, Benton R: Functional architecture of olfactory ionotropic glutamate receptors. Neuron. 2011, 69 (1): 44-60. 10.1016/j.neuron.2010.11.042.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Leal WS: Odorant reception in insects: Roles of receptors, binding proteins, and degrading enzymes. Annu Rev Entomol. 2013, 58 (1): 373-391. 10.1146/annurev-ento-120811-153635.

    Article  CAS  PubMed  Google Scholar 

  28. Vogt RG: Odorant binding protein homologues of the malaria mosquito Anopheles gambiae; Possible orthologues of the OS-E and OS-F OBPs of Drosophila melanogaster. J Chem Ecol. 2002, 28 (11): 2371-2376. 10.1023/A:1021009311977.

    Article  CAS  PubMed  Google Scholar 

  29. 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 (6): 549-560. 10.1046/j.1365-2583.2003.00440.x.

    Article  CAS  PubMed  Google Scholar 

  30. Zhou J-J, Huang W, Zhang G-A, Pickett JA, Field LM: “Plus-C” odorant-binding protein genes in two Drosophila species and the malaria mosquito Anopheles gambiae. Gene. 2004, 327 (1): 117-129. 10.1016/j.gene.2003.11.007.

    Article  CAS  PubMed  Google Scholar 

  31. Biessmann H, Nguyen OK, Le D, Walter MF: Microarray-based survey of a subset of putative olfactory genes in the mosquito Anopheles gambiae. Insect Mol Biol. 2005, 14: 575-589. 10.1111/j.1365-2583.2005.00590.x.

    Article  CAS  PubMed  Google Scholar 

  32. Kopp A, Barmina O, Hamilton AM, Higgins L, McIntyre LM, Jones CD: Evolution of gene expression in the Drosophila olfactory system. Mol Biol Evol. 2008, 25: 1081-1092. 10.1093/molbev/msn055.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Dekker T, Ibba I, Siju KP, Stensmyr MC, Hansson BS: Olfactory shifts parallel superspecialism for toxic fruit in Drosophila melanogaster sibling D sechellia. Curr Biol. 2006, 16: 101-109. 10.1016/j.cub.2005.11.075.

    Article  CAS  PubMed  Google Scholar 

  34. McBride CS: Rapid evolution of smell and taste receptor genes during host specialization in Drosophila sechellia. Proc Natl Acad Sci U S A. 2007, 104: 4996-5001. 10.1073/pnas.0608424104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. McBride CS, Arguello JR: Five Drosophila genomes reveal nonneutral evolution and the signature of host specialization in the chemoreceptor super familiy. Genetics. 2007, 177: 1395-1416. 10.1534/genetics.107.078683.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Rinker DC, Zhou X, Pitts RJ, Rokas A, Zwiebel LJ: Antennal transcriptome profiles of anopheline mosquitoes reveal human host olfactory specialization in Anopheles gambiae. BMC Genomics. 2013, 14 (1): 1-10.1186/1471-2164-14-1.

    Article  Google Scholar 

  37. Rund SSC, Hou TY, Ward SM, Collins FH, Duffield GE: Genome-wide profiling of diel and circadian gene expression in the malaria vector Anopheles gambiae. Proc Natl Acad Sci U S A. 2011, 108 (32): E421-E430. 10.1073/pnas.1100584108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Dekker T, Takken W: Differential responses of mosquito sibling species Anopheles arabiensis and An. quadriannulatus to carbon dioxide, a man or a calf. Med Vet Entomol. 1998, 12 (2): 136-140. 10.1046/j.1365-2915.1998.00073.x.

    Article  CAS  PubMed  Google Scholar 

  39. Sharp BL, Quicke FC, Jansen EJ: Aspects of the behaviour of five anopheline species. J Entomol Soc Southern Africa. 1984, 47: 251-258.

    Google Scholar 

  40. Pates HV, Curtis CF, Takken W: Hybridization studies to modify the host preference of Anopheles gambiae. Med Vet Entomol. 2014, 28 (S1): 68-74. 10.1111/mve.12070.

    Article  PubMed  Google Scholar 

  41. Zhang YV, Ni J, Montell C: The molecular basis for attractive salt-taste coding in Drosophila. Science. 2013, 340: 1334-1338. 10.1126/science.1234133.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Foster WA, Takken W: Nectar-related vs. human-related volatiles: behavioral response and choice by female and male Anopheles gambiae (Diptera: Culicidae) between emergence and first feeding. Bull Entomol Res. 2004, 94: 145-157. 10.1079/BER2003288.

    Article  CAS  PubMed  Google Scholar 

  43. Coluzzi M: Heterogeneities of the malaria vectorial system in tropical Africa and their significance in malaria epidemiology and control. Bull World Health Organ. 1984, 62 (suppl): 107-113.

    PubMed Central  PubMed  Google Scholar 

  44. Gillies MT, Coetzee M:Publication of the South African Institute for Medical Research, Johannesburg A supplement to the Anophilinae of Africa south of the Sahara 1987, 55-

  45. Takken W, Klowden MJ, Chambers G: Effect of body size on host seeking and blood meal utilization in Anopheles gambiae sensu Stricto (Diptera: Culicidae): the disadvantage of being small. J Med Entomol. 1998, 35: 639-645.

    Article  CAS  PubMed  Google Scholar 

  46. Ramoni R, Vincent F, Grolli S, Conti V, Malosse C, Boyer F-D, Nagnan-Le Meillour P, Spinelle S, Cambillau C, Tegoni M: The insect attractant 1-octen-3-ol is the natural ligand of the bovine odorant binding protein. J Biol Chem. 2001, 276 (10): 7150-7155. 10.1074/jbc.M010368200.

    Article  CAS  PubMed  Google Scholar 

  47. Xu P, Atkinson R, Jones DNM, Smith DP: Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron. 2005, 45 (2): 193-200. 10.1016/j.neuron.2004.12.031.

    Article  CAS  PubMed  Google Scholar 

  48. Matsuo T, Sugaya S, Yasukawa J, Aigaki T, Fuyama Y: Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol. 2007, 5: e118-10.1371/journal.pbio.0050118.

    Article  PubMed Central  PubMed  Google Scholar 

  49. Syed Z, Leal WS: Maxillary palps are broad spectrum odorant detectors in Culex quinquefasciatus. Chem Senses. 2007, 32 (8): 727-738. 10.1093/chemse/bjm040.

    Article  CAS  PubMed  Google Scholar 

  50. Auer PL, Doerge RW: Statistical Design and Analysis of RNA Sequencing Data. Genetics. 2010, 185 (2): 405-416. 10.1534/genetics.110.114983.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Jones WD, Cayirlioglu P, Grunwald Kadow I, Vosshall LB: Two chemosensory receptors together mediate carbon dioxide detection in Drosophila. Nature. 2006, 445 (7143): 86-90.

    PubMed  Google Scholar 

  52. Pitts RJ, Zwiebel LJ: Antennal sensilla of two female Anopheline sibling species with differing host ranges. Malar J. 2006, 5: 26-10.1186/1475-2875-5-26.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Spitzen J, Smallegange RC, Takken W: Effect of human odours and positioning of CO2 release point on trap catches of the malaria mosquito Anopheles gambiae sensu stricto in an olfactometer. Physiol Entomol. 2008, 33 (2): 116-122. 10.1111/j.1365-3032.2008.00612.x.

    Article  CAS  Google Scholar 

  54. Takken W: The role of olfaction in host-seeking of mosquitoes: a review. Insect Sci Applic. 1991, 12: 287-290.

    Google Scholar 

  55. Verhulst NO, Mbadi PA, Kiss GB, Mukabana WR, van Loon JJ, Takken W, Smallegange RC: Improvement of a synthetic lure for Anopheles gambiae using compounds produced by humanskin microbiota. Malaria J. 2011, 10 (1): 28-10.1186/1475-2875-10-28.

    Article  CAS  Google Scholar 

  56. Coetzee M, Hunt RH, Wilkerson R, della Torre A, Coulibaly MB, Besansky NJ: Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013, 3619 (3): 246-274.

    Article  PubMed  Google Scholar 

  57. Pates HV, Takken W, Stuke K, Curtis CF: Differential behaviour of Anopheles gambiae sensu stricto (Diptera : Culicidae) to human and cow odours in the laboratory. B Entomol Res. 2001, 91 (4): 289-296. 10.1079/BER200198.

    Article  CAS  Google Scholar 

  58. Patel RK, Jain M: NGS QC toolkit: A toolkit for quality control of next generation sequencing data. PLoS One. 2012, 7 (2): e30619-10.1371/journal.pone.0030619.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Sonali J, Batut P, Chaisson M, Gingeras TR: STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013, 29: 15-21. 10.1093/bioinformatics/bts635.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Love MI, Huber W, Anders S: Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. BioRxiv. 2014, doi:10.1101/002832

    Google Scholar 

  61. Conesa A, Gotz S, Garcia-Gomez J, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to Leon Westerd for assistance with colony maintenance, to Jeroen Spitzen for providing odor samples, and to Craig Coates for providing rearing space. We also thank two anonymous reviewers for their help in improving an earlier version of this manuscript. This work was supported by NIH/NIAID grant 1R01 AI085079 to MAS.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Michel A Slotman.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TKH: data collection, analysis and manuscript preparation; LVC: data collection, data analysis and manuscript preparation, GA; data collection, data analysis and manuscript preparation; Sharmila Pathikonda; data analysis and manuscript preparation; WT: supervision and planning of experiments, and manuscript preparation; MAS: conceived, planned and supervised experiments, manuscript preparation. All authors read and approved the final manuscript.

A retraction note to this article can be found online at http://dx.doi.org/10.1186/s12864-015-1802-z.

Electronic supplementary material

12864_2014_6806_MOESM1_ESM.pdf

Additional file 1: Figure S1: PCA plot showing the antennal and maxillary palps data sets of An. gambiae and An. quadriannulatus in the 2D plane spanned by their first two principal components. (PDF 41 KB)

12864_2014_6806_MOESM2_ESM.pdf

Additional file 2: Figure S2: GO analysis of 608 An. quadriannulatus genes with >2-fold antennal expression compared to An. gambiae, predicting their involvement in molecular functions (A) and biological processes (B). Data are presented as level 3 GO categorization. Categories with less than 1% of representation were grouped in “others”. (PDF 128 KB)

12864_2014_6806_MOESM3_ESM.pdf

Additional file 3: Figure S3: GO analysis of 870 An. gambiae genes with >2-fold maxillary palp expression compared to An. quadriannulatus, predicting their involvement in molecular functions (A) and biological processes (B), and of 787 An. quadriannulatus genes with >2-fold maxillary palp expression compared to An. gambiae, predicting their involvement in molecular functions (C) and biological processes (D). (PDF 247 KB)

12864_2014_6806_MOESM4_ESM.pdf

Additional file 4: Figure S4: Venn diagram showing the overlap in the number of genes significantly and more than 2-fold higher expressed in the tissues of either species. (PDF 79 KB)

12864_2014_6806_MOESM5_ESM.xlsx

Additional file 5: Table S1: Gene expression data for all genes in the antennae of An. gambiae and An. quadriannulatus.(XLSX 2 MB)

12864_2014_6806_MOESM6_ESM.xlsx

Additional file 6: Table S2: Gene expression data for all genes in the maxillary palps of An. gambiae and An. quadriannulatus.(XLSX 2 MB)

Additional file 7: Scripts used to run RNAseq analyses.(DOCX 110 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  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 https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://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.

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hodges, T.K., Cosme, L.V., Athrey, G. et al. RETRACTED ARTICLE: Species-specific chemosensory gene expression in the olfactory organs of the malaria vector Anopheles gambiae. BMC Genomics 15, 1089 (2014). https://doi.org/10.1186/1471-2164-15-1089

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-1089

Keywords