Selective pressures on Drosophila chemosensory receptor genes

Genome-wide identification and characterization of the chemosensory relative protein genes in Rhus gall aphid Schlechtendalia chinensis

BACKGROUND

The Rhus gall aphid Schlechtendalia chinensis specially uses the only species Rhus chinensis and certain moss species (Mniaceae) as its primary host plant and secondary host plants, respectively. Rhus galls are formed on the primary host by the sucking of aphids, and used in traditional medicine as well as other various areas due to their high tannin contents. Chemoreception is critical for insect behaviors such as host searching, location and identification of mates and reproductive behavior. The process of chemoreception is mediated by a series of protein gene families, including odorant-binding proteins (OBPs), chemosensory proteins (CSPs), olfactory receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), and sensory neuron membrane proteins (SNMPs). However, there have been no reports on the analysis of molecular components related to the chemoreception system of S. chinensis at the genome level.

RESULTS

We examined the genes of eight OBPs, nine CSPs, 24 ORs, 16 GRs, 22 IRs, and five SNMPs in the S. chinensis genome using homological searches, and these chemosensory genes appeared mostly on chromosome 1. Phylogenetic and gene number analysis revealed that the gene families, e.g., ORs, GRs, CSPs and SNMPs in S. chinensis, have experienced major contractions by comparing to Myzus persicae, while the two gene families OBPs and IRs had slight expansion. The current results might be related to the broader host range of M. persicae versus the specialization of S. chinensis on only a host plant. There were 28 gene pairs between genomes of S. chinensis and Acyrthosiphon pisum in the chemoreceptor gene families by collinear comparison. Ka/Ks ratios (< 1) indicated that the genes of S. chinensis were mainly affected by purification selection during evolution. We also found the lower number and expression level of chemoreception genes in S. chinensis than in other 11 aphid species, such as ORs, GRs and IRs, which play an important role in host search.

CONCLUSION

Our study firstly identified the genes of the different chemosensory protein gene families in the S. chinensis genome, and analyzed their general features and expression profile, demonstrating the importance of chemoreception in the aphid and providing new information for further functional research.

Genome-wide identification and expression profiling of odorant receptor genes in the malaria vector Anopheles sinensis

Background

The olfactory system plays a crucial role in regulating insect behaviors. The detection of odorants is mainly mediated by various odorant receptors (ORs) that are expressed in the dendrites of olfactory neurons of chemosensilla. Anophelessinensis is a major malaria vector in Eastern Asia and its genome has recently been successfully sequenced and annotated. In this study, we present genome-wide identification and expression profiling of OR genes in different chemosensory tissues of An.sinensis.

Methods

The OR genes were identified using the available genome sequences of An.sinensis. A series of bioinformatics analyses were conducted to investigate the structure, genome distribution, selective pressure and phylogenetic relationships of OR genes, the conserved domains and specific functional sites in the OR amino acid sequences. The expression levels of OR genes were analyzed from transcriptomic data from An.sinensis antennae, proboscis and maxillary palps of both sexes.

Results

A total of 59 putative OR genes have been identified and characterized in An.sinensis. This number is significantly less than that in An.gambiae. Whether this difference is caused by the contraction or expansion of OR genes after divergence of the two species remains unknown. The RNA-seq analysis showed that AsORs have obvious tissue- and sex-specific expression patterns. Most AsORs are highly expressed in the antennae and the expression pattern and number of AsORs expressed in antennae are similar in males and females. However, the relative levels of AsOR transcripts are much higher in female antennae than in male antennae, which indicates that the odor sensitivity is likely to be increased in female mosquitoes. Based on the expression patterns and previous studies, we have speculated on the functions of some OR genes but this needs to be validated by further behavioral, molecular and electrophysiological studies. Further studies are necessary to compare the olfactory-driven behaviors and identify receptors that respond strongly to components of human odors that may act in the process of human recognition.

Conclusions

This is the first genome-wide analysis of the entire repertoire of OR genes in An.sinensis. Characterized features and profiled expression patterns of ORs suggest their involvement in the odorous reception of this species. Our findings provide a basis for further research on the functions of OR genes and additional genetic and behavioral targets for more sustainable management of An.sinensis in the future.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05259-x.

The Evolutionary Dynamics of the Odorant Receptor Gene Family in Corbiculate Bees

Insects rely on chemical information to locate food, choose mates, and detect potential predators. It has been hypothesized that adaptive changes in the olfactory system facilitated the diversification of numerous insect lineages. For instance, evolutionary changes of Odorant Receptor (OR) genes often occur in parallel with modifications in life history strategies. Corbiculate bees display a diverse array of behaviors that are controlled through olfaction, including varying degrees of social organization, and manifold associations with floral resources. Here we investigated the molecular mechanisms driving the evolution of the OR gene family in corbiculate bees in comparison to other chemosensory gene families. Our results indicate that the genomic organization of the OR gene family has remained highly conserved for ∼80 Myr, despite exhibiting major changes in repertoire size among bee lineages. Moreover, the evolution of OR genes appears to be driven mostly by lineage-specific gene duplications in few genomic regions that harbor large numbers of OR genes. A selection analysis revealed that OR genes evolve under positive selection, with the strongest signals detected in recently duplicated copies. Our results indicate that chromosomal translocations had a minimal impact on OR evolution, and instead local molecular mechanisms appear to be main drivers of OR repertoire size. Our results provide empirical support to the longstanding hypothesis that positive selection shaped the diversification of the OR gene family. Together, our results shed new light on the molecular mechanisms underlying the evolution of olfaction in insects.

Odorant Receptors of the New Zealand Endemic Leafroller Moth Species Planotortrix octo and P. excessana

Moths use their sense of smell to find food sources, mating partners and oviposition sites. For this they possess a family of odorant receptors (ORs). Some ORs are used by both sexes whereas others have sex-specific roles. For example, male moths possess ORs specifically tuned to sex pheromones produced by conspecific females. Here we identify sets of ORs from the antennae of New Zealand endemic leafroller moths Planotortrix octo (48 ORs) and P. excessana (47 ORs) using an RNA-Seq approach. Two orthologous ORs show male-biased expression in the adult antennae of both species (OR7 and OR30) and one other OR in each species was female-biased in its expression (PoctOR25, PexcOR14) by qPCR. PAML analysis conducted on male-biased ORs indicated positive selection acting on the male-biased OR7. The fact that OR7 is likely under positive selection, that it is male-biased in its expression and that its orthologue in C. obliquana, CoblOR7, responds to sex pheromone components also utilised by Planotortrix species, suggests that this receptor may also be important in sex pheromone reception in Planotortrix species.

Rapid evolution of chemosensory receptor genes in a pair of sibling species of orchid bees (Apidae: Euglossini)

Background

Insects rely more on chemical signals (semiochemicals) than on any other sensory modality to find, identify, and choose mates. In most insects, pheromone production is typically regulated through biosynthetic pathways, whereas pheromone sensory detection is controlled by the olfactory system. Orchid bees are exceptional in that their semiochemicals are not produced metabolically, but instead male bees collect odoriferous compounds (perfumes) from the environment and store them in specialized hind-leg pockets to subsequently expose during courtship display. Thus, the olfactory sensory system of orchid bees simultaneously controls male perfume traits (sender components) and female preferences (receiver components). This functional linkage increases the opportunities for parallel evolution of male traits and female preferences, particularly in response to genetic changes of chemosensory detection (e.g. Odorant Receptor genes). To identify whether shifts in pheromone composition among related lineages of orchid bees are associated with divergence in chemosensory genes of the olfactory periphery, we searched for patterns of divergent selection across the antennal transcriptomes of two recently diverged sibling species Euglossa dilemma and E. viridissima.

Results

We identified 3185 orthologous genes including 94 chemosensory loci from five different gene families (Odorant Receptors, Ionotropic Receptors, Gustatory Receptors, Odorant Binding Proteins, and Chemosensory Proteins). Our results revealed that orthologs with signatures of divergent selection between E. dilemma and E. viridissima were significantly enriched for chemosensory genes. Notably, elevated signals of divergent selection were almost exclusively observed among chemosensory receptors (i.e. Odorant Receptors).

Conclusions

Our results suggest that rapid changes in the chemosensory gene family occurred among closely related species of orchid bees. These findings are consistent with the hypothesis that strong divergent selection acting on chemosensory receptor genes plays an important role in the evolution and diversification of insect pheromone systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0451-9) contains supplementary material, which is available to authorized users.

Pheromone receptor evolution in the cryptic leafroller species, Ctenopseustis obliquana and C. herana

How new mate recognition systems evolve when changes are required in both the male and female components remains a conundrum. Here, we investigated the molecular basis of pheromone reception in two species of tortricid (leafroller) moth, Ctenopseustis obliquana and C. herana. Male C. obliquana are attracted to a 90:10 blend of (Z)-8-tetradecenyl acetate (Z8-14:OAc) and (Z)-5-tetradecenyl acetate (Z5-14:OAc), whereas C. herana males are attracted to Z5-14:OAc alone. We used a transcriptome sequencing approach from adult male and female antennae to identify 47 olfactory receptors (ORs) from each species and assessed their expression levels in male and female antennae using RNA-Seq counting and quantitative RT-PCR. Three male-biased and one female-biased OR were identified in C. obliquana by quantitative RT-PCR, and four male-biased and one female-biased receptor in C. herana. The male-biased receptors, CoblOR7, CoblOR30, CherOR7, CherOR30, CherOR1a and CherOR1b were tested for their ability to respond to sex pheromone components in a HEK293 cell calcium assay. CoblOR7 and CherOR7 responded to Z8-14:OAc, however, no receptor for Z5-14:OAc was identified. In addition to Z8-14:OAc, CherOR7 also responded to Z7-14:OAc, indicating that this receptor may be under relaxed constraint. Of the 29 amino acid differences between CoblOR7 and CherOR7, significantly more are located in the third and the sixth transmembrane domain regions. Overall, these findings are consistent with studies revealing the presence of neurons tuned to both Z8-14:OAc and Z5-14:OAc in both species, but that for C. herana males, the ability to detect Z8-14:OAc is currently not required.

A single gene affects both ecological divergence and mate choice in Drosophila

Evolutionary changes in traits involved in both ecological divergence and mate choice may produce reproductive isolation and speciation. However, there are few examples of such dual traits, and the genetic and molecular bases of their evolution have not been identified. We show that methyl-branched cuticular hydrocarbons (mbCHCs) are a dual trait that affects both desiccation resistance and mate choice in Drosophila serrata. We identify a fatty acid synthase mFAS (CG3524) responsible for mbCHC production in Drosophila and find that expression of mFAS is undetectable in oenocytes (cells that produce CHCs) of a closely related, desiccation-sensitive species, D. birchii, due in part to multiple changes in cis-regulatory sequences of mFAS. We suggest that ecologically influenced changes in the production of mbCHCs have contributed to reproductive isolation between the two species.

Functional evolution of a multigene family: orthologous and paralogous pheromone receptor genes in the turnip moth, Agrotis segetum

Lepidopteran pheromone receptors (PRs), for which orthologies are evident among closely related species, provide an intriguing example of gene family evolution in terms of how new functions may arise. However, only a limited number of PRs have been functionally characterized so far and thus evolutionary scenarios suffer from elements of speculation. In this study we investigated the turnip moth Agrotis segetum, in which female moths produce a mixture of chemically related pheromone components that elicit specific responses from receptor cells on male antennae. We cloned nine A. segetum PR genes and the Orco gene by degenerate primer based RT-PCR. The nine PR genes, named as AsegOR1 and AsegOR3-10, fall into four distinct orthologous clusters of known lepidopteran PRs, of which one contains six paralogues. The paralogues are under relaxed selective pressure, contrasting with the purifying selection on other clusters. We identified the receptors AsegOR9, AsegOR4 and AsegOR5, specific for the respective homologous pheromone components (Z)-5-decenyl, (Z)-7-dodecenyl and (Z)-9-tetradecenyl acetates, by two-electrode voltage clamp recording from Xenopus laevis oocytes co-expressing Orco and each PR candidate. These receptors occur in three different orthologous clusters. We also found that the six paralogues with high sequence similarity vary dramatically in ligand selectivity and sensitivity. Different from AsegOR9, AsegOR6 showed a relatively large response to the behavioural antagonist (Z)-5-decenol, and a small response to (Z)-5-decenyl acetate. AsegOR1 was broadly tuned, but most responsive to (Z)-5-decenyl acetate, (Z)-7-dodecenyl acetate and the behavioural antagonist (Z)-8-dodecenyl acetate. AsegOR8 and AsegOR7, which differ from AsegOR6 and AsegOR1 by 7 and 10 aa respectively, showed much lower sensitivities. AsegOR10 showed only small responses to all the tested compounds. These results suggest that new receptors arise through gene duplication, and relaxed evolutionary constraints or positive selection among paralogues allow functional divergence to occur in spite of purifying selection being the norm.

Genome-wide patterns of natural variation reveal strong selective sweeps and ongoing genomic conflict in Drosophila mauritiana

Although it is well understood that selection shapes the polymorphism pattern in Drosophila, signatures of classic selective sweeps are scarce. Here, we focus on Drosophila mauritiana, an island endemic, which is closely related to Drosophila melanogaster. Based on a new, annotated genome sequence, we characterized the genome-wide polymorphism by sequencing pooled individuals (Pool-seq). We show that the interplay between selection and recombination results in a genome-wide polymorphism pattern characteristic for D. mauritiana. Two large genomic regions (>500 kb) showed the signature of almost complete selective sweeps. We propose that the absence of population structure and limited geographic distribution could explain why such pronounced sweep patterns are restricted to D. mauritiana. Further evidence for strong adaptive evolution was detected for several nucleoporin genes, some of which were not previously identified as genes involved in genomic conflict. Since this adaptive evolution is continuing after the split of D. mauritiana and Drosophila simulans, we conclude that genomic conflict is not restricted to short episodes, but rather an ongoing process in Drosophila.

Sequence comparisons of odorant receptors among tortricid moths reveal different rates of molecular evolution among family members

In insects, odorant receptors detect volatile cues involved in behaviours such as mate recognition, food location and oviposition. We have investigated the evolution of three odorant receptors from five species within the moth genera Ctenopseustis and Planotrotrix, family Tortricidae, which fall into distinct clades within the odorant receptor multigene family. One receptor is the orthologue of the co-receptor Or83b, now known as Orco (OR2), and encodes the obligate ion channel subunit of the receptor complex. In comparison, the other two receptors, OR1 and OR3, are ligand-binding receptor subunits, activated by volatile compounds produced by plants - methyl salicylate and citral, respectively. Rates of sequence evolution at non-synonymous sites were significantly higher in OR1 compared with OR2 and OR3. Within the dataset OR1 contains 109 variable amino acid positions that are distributed evenly across the entire protein including transmembrane helices, loop regions and termini, while OR2 and OR3 contain 18 and 16 variable sites, respectively. OR2 shows a high level of amino acid conservation as expected due to its essential role in odour detection; however we found unexpected differences in the rate of evolution between two ligand-binding odorant receptors, OR1 and OR3. OR3 shows high sequence conservation suggestive of a conserved role in odour reception, whereas the higher rate of evolution observed in OR1, particularly at non-synonymous sites, may be suggestive of relaxed constraint, perhaps associated with the loss of an ancestral role in sex pheromone reception.

Evolutionary Genomics of Genes Involved in Olfactory Behavior in the Drosophila melanogaster Species Group

Previous comparative genomic studies of genes involved in olfactory behavior in Drosophila focused only on particular gene families such as odorant receptor and/or odorant binding proteins. However, olfactory behavior has a complex genetic architecture that is orchestrated by many interacting genes. In this paper, we present a comparative genomic study of olfactory behavior in Drosophila including an extended set of genes known to affect olfactory behavior. We took advantage of the recent burst of whole genome sequences and the development of powerful statistical tools to analyze genomic data and test evolutionary and functional hypotheses of olfactory genes in the six species of the Drosophila melanogaster species group for which whole genome sequences are available. Our study reveals widespread purifying selection and limited incidence of positive selection on olfactory genes. We show that the pace of evolution of olfactory genes is mostly independent of the life cycle stage, and of the number of life cycle stages, in which they participate in olfaction. However, we detected a relationship between evolutionary rates and the position that the gene products occupy in the olfactory system, genes occupying central positions tend to be more constrained than peripheral genes. Finally, we demonstrate that specialization to one host does not seem to be associated with bursts of adaptive evolution in olfactory genes in D. sechellia and D. erecta, the two specialists species analyzed, but rather different lineages have idiosyncratic evolutionary histories in which both historical and ecological factors have been involved.

Antennal transcriptome of Manduca sexta

In recent years, considerable progress has been made in understanding the molecular mechanisms underlying olfaction in insects. Because of the diverse nature of the gene families involved, this process has largely relied on genomic data. As a consequence, studies have focused on a small subset of species with extensive genomic information. For Lepidoptera, a large order historically crucial to olfactory research, this circumstance has mostly limited advances to the domesticated species Bombyx mori, with some progress in the noctuid Heliothis virescens based on a nonpublic partial genome database. Because of the limited behavioral repertoire and nonexistent ecological importance of Bombyx, molecular data on the tobacco hornworm Manduca sexta are of utmost importance, especially with regards to its position as a classical olfactory model and its complex natural behavior. Here we present the use of transcriptomic and microarray data to identify members of the main olfactory gene families of Manduca. To assess the quality of our data, we correlate information on expressed receptor genes with detailed morphological data on the antennal lobe. Finally, we compare the expression of the near-complete transcript sets in male and female antennae.

The insect chemoreceptor superfamily of the parasitoid jewel wasp Nasonia vitripennis

Chemoreception is important for locating food, mates and other resources in many insects, including the parasitoid jewel wasp Nasonia vitripennis. In the insect chemoreceptor superfamily, Nasonia has 58 gustatory receptor (Gr) genes, of which 11 are pseudogenes, leaving 47 apparently intact proteins encoded. No carbon dioxide receptors, two candidate sugar receptors, a DmGr43a orthologue, and several additional Gr lineages were identified, including significant gene subfamily expansions related to the 10 Grs found in the honey bee Apis mellifera. Nasonia has a total of 301 odorant receptor (Or) genes, of which 76 are pseudogenes, leaving 225 apparently intact Ors. Phylogenetic comparison with the 174 honey bee Ors reveals differential gene subfamily expansion in each hymenopteran lineage, along with a few losses from each species. The only simple orthologous relationship is the expected single DmOr83b orthologue. The large number of Nasonia Ors is the result of several major subfamily expansions, including one of 55 genes. Nasonia does not have the elaborate social chemical communication of honey bees, nor the diversity of floral odours honey bees detect, however, Nasonia wasps might need to detect a diversity of odours to find potential mates and hosts or avoid harmful substances in its environment.

Functional and molecular evolution of olfactory neurons and receptors for aliphatic esters across the Drosophila genus

Insect olfactory receptor (Or) genes are large, rapidly evolving gene families of considerable interest for evolutionary studies. They determine the responses of sensory neurons which mediate critical behaviours and ecological adaptations. We investigated the evolution across the genus Drosophila of a subfamily of Or genes largely responsible for the perception of ecologically relevant aliphatic esters; products of yeast fermentation and fruits. Odour responses were recorded from eight classes of olfactory receptor neurons known to express this Or subfamily in D. melanogaster and from homologous sensilla in seven other species. Despite the fact that these species have diverged over an estimated 40 million years, we find that odour specificity is largely maintained in seven of the eight species. In contrast, we observe extensive changes in most neurons of the outgroup species D. virilis, and in two neurons across the entire genus. Some neurons show small shifts in specificity, whilst some dramatic changes correlate with gene duplication or loss. An olfactory receptor neuron response similarity tree did not match an Or sequence similarity tree, but by aligning Or proteins of likely functional equivalence we identify residues that may be relevant for odour specificity. This will inform future structure-function studies of Drosophila Ors.

Drosophila sechellia as a model in chemosensory neuroecology

The olfactory system directly interfaces with the environment, and thus changes in the environment or an animal's habits would presumably also lead to changes in the olfactory system. Comparative studies on specialized animals with known generalist ancestors could hence be a way of revealing general processes shaping olfactory systems as well as highlighting the importance and function of specific chemosensory genes. Drosophila sechellia, a close relative of D. melanogaster, and highly specialized towards the use of the toxic Morinda fruit, has been the subject of several studies aimed at pinpointing the consequences of host choice to the chemosensory system, from the molecular level to behavior. This paper highlights some recent (and not so recent) advances in this area and concludes with a brief discussion on the antiquity of the D. sechellia-Morinda connection.

The insect SNMP gene family

SNMPs are membrane proteins observed to associate with chemosensory neurons in insects; in Drosophila melanogaster, SNMP1 has been shown to be essential for the detection of the pheromone cis-vaccenyl acetate (CVA). SNMPs are one of three insect gene clades related to the human fatty acid transporter CD36. We previously characterized the CD36 gene family in 4 insect Orders that effectively cover the Holometabola, or some 80% of known insect species and the 300 million years of evolution since this lineage emerged: Lepidoptera (e.g. Bombyx mori, Antheraea polyphemus, Manduca sexta, Heliothis virescens, Helicoverpa assulta, Helicoverpa armigera, Mamestra brassicae); Diptera (D. melanogaster, Drosophila pseudoobscura, Aedes aegypti, Anopheles gambiae, Culex pipiens quinquefasciatus); Hymenoptera (Apis mellifera); and Coleoptera (Tribolium castaneum). This previous study suggested a complex topography within the SNMP clade including a strongly supported SNMP1 sub-clade plus additional SNMP genes. To further resolve the SNMP clade here, we used cDNA sequences of SNMP1 and SNMP2 from various Lepidoptera species, D. melanogaster and Ae. aegypti, as well as BAC derived genomic sequences from Ae. aegypti as models for proposing corrected sequences of orthologues in the D. pseudoobscura and An. gambiae genomes, and for identifying orthologues in the B. mori and C. pipiens q. genomes. We then used these sequences to analyze the SNMP clade of the insect CD36 gene family, supporting the existence of two well supported sub-clades, SNMP1 and SNMP2, throughout the dipteran and lepidopteran lineages, and plausibly throughout the Holometabola and across a broad evolutionary time scale. We present indirect evidence based on evolutionary selection (dN/dS) that the dipteran SNMPs are expressed as functional proteins. We observed expansions of the SNMP1 sub-clade in C. pipiens q. and T. castaneum suggesting that the SNMP1s may have an expanded functional role in these species.

The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessity

Chemosensory receptors are essential for the survival of organisms that range from bacteria to mammals. Recent studies have shown that the numbers of functional chemosensory receptor genes and pseudogenes vary enormously among the genomes of different animal species. Although much of the variation can be explained by the adaptation of organisms to different environments, it has become clear that a substantial portion is generated by genomic drift, a random process of gene duplication and deletion. Genomic drift also generates a substantial amount of copy-number variation in chemosensory receptor genes within species. It seems that mutation by gene duplication and inactivation has important roles in both the adaptive and non-adaptive evolution of chemosensation.

On the scent of speciation: the chemosensory system and its role in premating isolation

Chemosensory speciation is characterized by the evolution of barriers to genetic exchange that involve chemosensory systems and chemical signals. Here, we review some representative studies documenting chemosensory speciation in an attempt to evaluate the importance and the different aspects of the process in nature and to gain insights into the genetic basis and the evolutionary mechanisms of chemosensory trait divergence. Although most studies of chemosensory speciation concern sexual isolation mediated by pheromone divergence, especially in Drosophila and moth species, other chemically based behaviours (habitat choice, pollinator attraction) can also play an important role in speciation and are likely to do so in a wide range of invertebrate and vertebrate species. Adaptive divergence of chemosensory traits in response to factors such as pollinators, hosts and conspecifics commonly drives the evolution of chemical prezygotic barriers. Although the genetic basis of chemosensory speciation remains largely unknown, genomic approaches to chemosensory gene families and to enzymes involved in biosynthetic pathways of signal compounds now provide new opportunities to dissect the genetic basis of these complex traits and of their divergence among taxa.

Evolution of gene expression in the Drosophila olfactory system

Host plant shifts by phytophagous insects play a key role in insect evolution and plant ecology. Such shifts often involve major behavioral changes as the insects must acquire an attraction and/or lose the repulsion to the new host plant's odor and taste. The evolution of chemotactic behavior may be due, in part, to gene expression changes in the peripheral sensory system. To test this hypothesis, we compared gene expression in the olfactory organs of Drosophila sechellia, a narrow ecological specialist that feeds on the fruit of Morinda citrifolia, with its close relatives Drosophila simulans and Drosophila melanogaster, which feed on a wide variety of decaying plant matter. Using whole-genome microarrays and quantitative polymerase chain reaction, we surveyed the entire repertoire of Drosophila odorant receptors (ORs) and odorant-binding proteins (OBPs) expressed in the antennae. We found that the evolution of OR and OBP expression was accelerated in D. sechellia compared both with the genome average in that species and with the rate of OR and OBP evolution in the other species. However, some of the gene expression changes that correlate with D. sechellia's increased sensitivity to Morinda odorants may predate its divergence from D. simulans. Interspecific divergence of olfactory gene expression cannot be fully explained by changes in the relative abundance of different sensilla as some ORs and OBPs have evolved independently of other genes expressed in the same sensilla. A number of OR and OBP genes are upregulated in D. sechellia compared with its generalist relatives. These genes include Or22a, which likely responds to a key odorant of M. citrifolia, and several genes that are yet to be characterized in detail. Increased expression of these genes in D. sechellia may have contributed to the evolution of its unique chemotactic behavior.

Five Drosophila genomes reveal nonneutral evolution and the signature of host specialization in the chemoreceptor superfamily

The insect chemoreceptor superfamily comprises the olfactory receptor (Or) and gustatory receptor (Gr) multigene families. These families give insects the ability to smell and taste chemicals in the environment and are thus rich resources for linking molecular evolutionary and ecological processes. Although dramatic differences in family size among distant species and high divergence among paralogs have led to the belief that the two families evolve rapidly, a lack of evolutionary data over short time scales has frustrated efforts to identify the major forces shaping this evolution. Here, we investigate patterns of gene loss/gain, divergence, and polymorphism in the entire repertoire of approximately 130 chemoreceptor genes from five closely related species of Drosophila that share a common ancestor within the past 12 million years. We demonstrate that the overall evolution of the Or and Gr families is nonneutral. We also show that selection regimes differ both between the two families as wholes and within each family among groups of genes with varying functions, patterns of expression, and phylogenetic histories. Finally, we find that the independent evolution of host specialization in Drosophila sechellia and D. erecta is associated with a fivefold acceleration of gene loss and increased rates of amino acid evolution at receptors that remain intact. Gene loss appears to primarily affect Grs that respond to bitter compounds while elevated Ka/Ks is most pronounced in the subset of Ors that are expressed in larvae. Our results provide strong evidence that the observed phenomena result from the invasion of a novel ecological niche and present a unique synthesis of molecular evolutionary analyses with ecological data.