| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Genetics, Vol. 177, 1655-1665, November 2007, Copyright © 2007
doi:10.1534/genetics.107.079731
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



,
,1
* Department of Genetics,
W. M. Keck Center for Behavioral Biology and
Department of Zoology, North Carolina State University, Raleigh, North Carolina 27695 and
National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland 20894
1 Corresponding author: W. M. Keck Center for Behavioral Biology, Campus Box 7617, North Carolina State University, Raleigh, NC 27695-7617.
E-mail: anholt{at}ncsu.edu
| ABSTRACT |
|---|
|
|
|---|
43 glomeruli in the antennal lobe (VOSSHALL et al. 2000). Individual neurons express one, or rarely two, odorant receptors from a repertoire of 62 Or genes (CLYNE et al. 1999; VOSSHALL et al. 1999). Uniquely expressed odorant receptors dimerize with the common Or83b receptor, which is essential for transport and insertion of odorant receptors in the chemosensory dendritic membrane (LARSSON et al. 2004; BENTON et al. 2006). Olfactory sensory neurons that express the same odorant receptor converge on the same antennal lobe glomerulus (GAO et al. 2000; VOSSHALL et al. 2000; BHALERAO et al. 2003). Odorants must dissolve in the aqueous perilymph to reach their cognate membrane-associated odorant receptors. Their solubilization and transport is thought to be mediated by odorant-binding proteins (OBPs) that are secreted by support cells. The Drosophila genome encodes 51 OBPs with different spatial patterns of expression (MCKENNA et al. 1994; PIKIELNY et al. 1994; GALINDO and SMITH 2001), which contain a characteristic structural signature of conserved cysteines (GRAHAM and DAVIES 2002; HEKMAT-SCAFE et al. 2002).
Altered regulation of expression of different subsets of OBPs has been observed following mating (MCGRAW et al. 2004), exposure to starvation stress (HARBISON et al. 2005), during the development of alcohol tolerance after exposure to alcohol (MOROZOVA et al. 2006), as a correlated response to artificial selection for divergent levels of copulation latency (MACKAY et al. 2005) and aggression (EDWARDS et al. 2006), and as a consequence of pleiotropic effects arising from single P-element-induced mutations that affect olfactory behavior (ANHOLT et al. 2003).
Whereas the role of OBPs in pheromone recognition has been clearly defined for several insect systems, the precise functions of these abundantly expressed proteins in olfaction remain obscure. Ligand specificities (whether broadly or narrowly tuned), interactions with odorant receptors (for which there remains scant evidence to date), interrelationships among OBPs with overlapping molecular-response profiles and their functional correspondence (if any) with odorant receptors, the significance of altered expression of some OBPs in aggression (EDWARDS et al. 2006), mating behavior (MCGRAW et al. 2004; MACKAY et al. 2005), and alcohol sensitivity (MOROZOVA et al. 2006), all pose as yet unresolved questions. To date, only one OBP in Drosophila melanogaster, encoded by Lush, has been characterized functionally (KIM et al. 1998; XU et al. 2005). Flies homozygous for a deletion of the Lush gene do not avoid repellant concentrations of short-chain alcohols (KIM et al. 1998) and do not respond behaviorally or electrophysiologically to the aggregation pheromone 11-cis-vaccenyl acetate (XU et al. 2005). The Or67d receptor, expressed in a subset of trichoid sensilla, has been identified as the receptor for 11-cis-vaccenyl acetate (HA and SMITH 2006; KURTOVIC et al. 2007). Lush appears to be essential for delivering this pheromone to its receptor. Other insights into the functions of OBPs come from a recent study reporting that a polymorphism in Obp57e in D. sechellia determines preference for its host plant, Morinda citrifolia, and that D. melanogaster knock-out flies for Obp57e and Obp57d showed altered behaviors to hexanoic and octanoic acid produced by this plant (MATSUO et al. 2007).
Functional studies on OBPs have been hampered by the lack of Obp mutants, with the exception of Lush. Furthermore, if odorant recognition by OBPs is combinatorial, as is the case for odorant recognition by mammalian odorant receptors (MALNIC et al. 1999), functional redundancy may render a laborious "one-gene-at-a-time" approach less than satisfactory, as it would provide only partial insights into the role of any one member of this multigene family in mediating olfactory behavior.
We devised a strategy to overcome these challenges by taking advantage of naturally occurring mutations that have arisen during evolution and that segregate as polymorphic variants in nature. We established isofemale lines from a natural population and inbred them for 20 generations, thus minimizing genetic variation within lines while retaining naturally occurring variation among the lines. We sequenced 13 Obp genes, located in two chromosomal clusters on the second and third chromosome, and used statistical tests for deviations from neutrality to assess patterns of selection. We then assessed whether polymorphisms in these 13 Obp genes were associated with naturally occurring variation in olfactory response to a standard test odorant, benzaldehyde. Finally, we show that single nucleotide polymorphisms (SNPs) in regulatory and coding regions that are associated with variation in chemosensory behavior can impact the predicted structure of pre-mRNA.
| MATERIALS AND METHODS |
|---|
|
|
|---|
OBP sequences:
Genomic DNA was extracted with the Puregene DNA extraction kit (Gentra, Minneapolis) and PCR primers were designed to amplify overlapping coding regions and 5'- and 3'-untranslated regions of the Obp genes (Obp56a–i and Obp99a–d) for Drosophila melanogaster and for a D. simulans line originally collected in Florida City (FL) by Jerry Coyne. PCR products were purified using Qiaquick columns (QIAGEN, Valencia, CA) and amplified samples were sequenced. Sequences were aligned with the Vector NTI Suite 9.0 program (Informax, Frederick, MD) to identify polymorphic sites. Singletons were excluded from the association and linkage disequilibrium (LD) analyses.
Molecular population genetics:
Neutrality tests were performed using the DnaSP 4.10.3 program (ROZAS et al. 2003) (http://www.ub.es/dnasp). D. simulans sequences were compared to the sequences from the D. melanogaster population for the HKA test (HUDSON et al. 1987) and the McDonald–Kreitman test (MCDONALD and KREITMAN 1991). Estimates of Tajima's D (TAJIMA 1993), Fu and Li's D* and F* (FU and LI 1993), and Fay and Wu's H (FAY and WU 2000) take into account the calculated population recombination rate (HUDSON 1987). Coalescent simulation was used to estimate P-values (two-tailed tests) with 103 coalescent simulations of an infinite site locus conditioned on the sample size; these simulations are implemented for a fixed number of segregating sites. LD between SNPs was analyzed using TASSEL 2.0 software (http://sourceforge.net/projects/tassel). Fisher's exact test was used to determine whether the pairs of sites were in significant LD.
Behavioral assays:
Olfactory behavior was quantified by measuring responses to the standard odorant, benzaldehyde, in a well established "dipstick" assay that we (ANHOLT et al. 1996; MACKAY et al. 1996; FANARA et al. 2002) and others (DEVAUD 2003; STOCKINGER et al. 2005) have used previously. Pilot experiments on 5 of the lines over a range of benzaldehyde concentrations established that a concentration of 3.5% (v/v) provided optimal resolution for evaluating variation in olfactory behavior in these lines. We measured olfactory behavior of 4–10 day-old non-virgin flies from 193 wild-derived inbred lines in single-sex groups of five flies/replicate and 10 replicates/sex. All assays were conducted between 2:00 and 4:00 PM in a behavioral room at 25° and 70% humidity under white light. The experimental design was randomized such that measurements on individual lines were collected over several days to average environmental variation. Theoretically a score of 2.5 reflects indifference to the odorant. Note, however, that the precise determination of the boundary between indifference and attraction or avoidance is determined statistically when for example the distribution of scores from mutants is compared to that of a control.
Locomotor reactivity was assessed by JORDAN et al. (2007) by subjecting single flies to a mechanical disturbance by tapping the vial twice against a table and recording the amount of time the fly is active in the 45 sec immediately following the disturbance.
Quantitative genetic analysis of olfactory behavior:
We used ANOVA to partition sources of variation in olfactory behavior according to the model Y = µ + L + S + L x S + E, where µ is the overall mean, L is the random effect of line, S is the fixed effect of sex, L x S is the random effect of line by sex, and E is environmental error. The total genotypic variance among lines was estimated as
, where
is the among-line variance component and
is the variance for the line-by-sex interaction. The total phenotypic variance was estimated as
, where
is the environmental variance component. Broad-sense heritability was estimated as
(CARBONE et al. 2006). Narrow-sense heritability was estimated as
, where
(FALCONER and MACKAY 1996). The genetic correlation between males and females was calculated as rMF = covMF/
M
F, where covMF is the covariance of line means for the two sexes and
M and
F are the square roots of the genotypic variances for each sex. Analyses of variance and tests of significance were calculated using the Proc GLM procedure, and variance components were estimated using the Proc VARCOMP procedure in SAS (SAS Institute, Cary, NC).
Genotype–phenotype associations:
Association between polymorphisms and line means for olfactory behavior were analyzed using two way factorial ANOVA with the model Y = µ + S + M + M x S + E, in which µ is the overall mean, marker (M) and sex (S) are fixed effects, and E indicates error variance. We used permutation tests to determine random distributions under the null hypotheses of no association between Obp genotypes and olfactory behavior in response to benzaldehyde (CHURCHILL and DOERGE 1994). We performed two different permutation tests. To assess whether we observed more significant associations with olfactory behavior for each gene than expected by chance, we permuted the phenotypes among the markers 1000 times and recorded the number of significant associations at P < 0.05 for each permuted data set. To identify particular polymorphic sites with significant associations with behavior, we similarly permuted the data 1000 times and recorded the lowest P-value of each permuted data set. In both cases we used the 5% significance threshold of the permuted data sets to give an empirical type I error rate that accounts for multiple tests.
In cases where more than one polymorphism in a gene was associated with olfactory behavior, we tested for associations between haplotypes of these variants and line means for olfactory behavior by two-way ANOVA with the model Y = µ + S + H + H x S + E, in which haplotype (H) and sex (S) are fixed factors and E indicates error. We conducted post hoc analyses using least-square means to assess the effect of haplotypes and Tukey tests to control the experimentwise error rate.
The additive variance attributable to a marker polymorphism (
) was estimated as
, where a is one-half the difference in mean olfactory behavior between homozygous genotypes for the marker, q is the frequency of the rare marker allele, and p = 1 – q (FALCONER and MACKAY 1996).
Prediction of RNA secondary structure:
Secondary structures of the full length mRNA and pre-mRNA molecules transcribed from Obp99a, Obp99c, and Obp99d genes were predicted using Mfold and Afold, as described previously (NACKLEY et al. 2006; SHABALINA et al. 2006). Energy minimization was performed by a dynamic programming method that finds the secondary structure with the minimum free energy with sums composed of factors that include stacking and loop length (NACKLEY et al. 2006; SHABALINA et al. 2006). The RNA folding parameters were published by the Turner group (MATHEWS et al. 1999). Suboptimal stem-loop structures were analyzed by the Hybrid program (NAZIPOVA et al. 1995).
| RESULTS |
|---|
|
|
|---|
To address this question, we sequenced alleles of 13 Obp genes organized in two representative chromosomal clusters, Obp56a–i on the second chromosome, and Obp99a–d on the third chromosome, in wild-derived inbred lines of D. melanogaster. We initially sequenced 50 alleles of each gene. Preliminary analyses suggested associations with members of the Obp99 gene family and olfactory behavior (see below); therefore, we obtained additional sequences of 143 alleles for the Obp99 genes. We observed 299 SNPs and 18 insertion/deletion (indel) polymorphisms in this sample of Obp genes, with 154 and 163 polymorphic sites in the Obp56 and Obp99 gene clusters, respectively (Table 1). SNP numbers were highly variable, ranging from only a single SNP in Obp56f to 76 SNPs in Obp99c (Table 1).
|
) and the number of segregating sites (
w)] vary over an order of magnitude among the Obp genes (Table 1), suggesting that OBPs have experienced different evolutionary histories. We applied tests for deviation from selective neutrality to members of the Obp56 and Obp99 gene clusters (HUDSON et al. 1987; MCDONALD and KREITMAN 1991; FU and LI 1993; TAJIMA 1993; FAY and WU 2000) [corrected for the estimated recombination rate, R, (HUDSON 1987)] that are based on detecting reduction in genetic diversity over different evolutionary time scales. Null alleles of Obp56c that contain a premature stop codon at predicted amino acid position 17 segregate in this population with an allele frequency of 0.06. Since we expected these null alleles not to be under selection, we did not include them in tests for deviation from neutrality (our results, however, do not change when the null alleles are included in the analysis, because the software performing the tests treats the missing sequence as missing data in any case). Seven of the 13 Obp genes exhibited signatures of departure from neutrality (Obp56a, Obp56c, Obp56g, Obp56h, Obp99b, Obp99c, and Obp99d; Tables 1 and 2), although none of the P-values of these tests survives a strict Bonferroni correction for multiple tests. Tajima's D was significant for Obp99c, Fu and Li's D* for Obp56g and Obp56h, F* for Obp99c and Obp99d, and Fay and Wu's H test statistic indicated a recent selective sweep for Obp56c (Table 1). No significant deviation from neutral expectations was detected by the HKA test (Table 1), but the McDonald–Kreitman test showed deviations from neutrality for Obp56a, Obp99b, and Obp99d (Table 2). For 5 Obp genes (Obp56g, Obp56h, Obp99b, Obp99c, and Obp99d) the departure from neutrality was such that there were more polymorphisms segregating at intermediate frequency than expected under neutral mutation-drift balance, which could be attributable to balancing selection. In contrast, Obp56a and Obp56c show evidence of rapid evolution. Since there is no LD between Obp56a and Obp56c, the signatures of positive selection experienced by these 2 genes are likely not due to hitchhiking (supplemental Table S1 at http://www.genetics.org/supplemental/).
|
|
Quantitative genetics of olfactory behavior:
We quantified naturally occurring variation in olfactory behavior among the 193 Raleigh inbred lines and observed broad variation in responsiveness to benzaldehyde, which appeared normally distributed and ranged from attractant (<<2.5) to repellant responses (>>2.5) (Figure 2a). To test whether the low responses were specific to olfactory behavior or due to a general deficit in locomotion, we assessed the correlation between locomotor behavior in response to a mechanical stimulus, which exhibits considerable inbreeding depression (JORDAN et al. 2007) and olfactory behavior in these lines. The correlation between locomotion scores and olfactory response scores was not significantly different from zero (Figure 2b). Thus, the variation in olfactory behavior cannot be explained by inbreeding depression for general locomotor impairment.
|
) is high (0.621), demonstrating the sensitivity of olfactory behavior to uncontrollable environmental variance. The line variance component (
) is also high (0.435), reflecting substantial genetic variation among lines. The estimate of broad-sense heritability (H2) is H2 = 0.441. Assuming strict additivity, the estimate of narrow-sense heritability (h2) is h2 = 0.283. This value is larger than previously estimated in a different population (MACKAY et al. 1996; h2 = 0.084 and 0.134 for chromosome 1 and chromosome 3 substitution lines, respectively). The substantial genetic component to variation in olfactory behavior provides a favorable scenario for association analysis. The significant sex-by-line interaction term indicates that there is variation in sexual dimorphism in the response to benzaldehyde, although the cross-sex genetic correlation (rMF = 0.893; Table 3) is higher than observed in a previous study of chromosome substitution lines (MACKAY et al. 1996).
|
Association analyses with 193 alleles for each of the four Obp99 genes revealed individual genotype–phenotype associations with olfactory responsiveness that exceeded the permutation threshold for multiple testing in Obp99a, Obp99c, and Obp99d (Table 4, Figure 3). Since we found a significant line-by-sex interaction when we analyzed phenotypic variation in olfactory behavior in this population, we included a marker-by-sex interaction term in the association model. However, none of the significant SNPs had significant marker-by-sex interactions. The effects of the SNPs on olfactory behavior are quite large, ranging from 0.40 to 1.33 genetic standard deviation units and 0.26–0.88 phenotypic standard deviation units in the population of inbred lines. However, if we assume strict additivity, the individual polymorphisms explain only 3–6% of the total additive variance in olfactory behavior (Table 4).
|
|
Our results implicate Obp99a, Obp99c, and Obp99d in the recognition of benzaldehyde and suggest that these OBPs are either redundant for the recognition of this odorant, or more likely, that OBPs recognize odorants in a combinatorial manner similar to odorant recognition by odorant receptors (MALNIC et al. 1999). Thus, a given odorant would interact with multiple OBPs with different affinities in a concentration-dependent manner and a given OBP would recognize multiple odorants on the basis of its molecular response profile.
Effects of polymorphisms associated with olfactory behavior on predicted secondary pre-mRNA structure:
To gain insights into the mechanism by which synonymous or noncoding SNPs could affect phenotypic variation in olfactory behavior, we predicted secondary structures for pre-mRNA molecules transcribed from alternative SNP genotypes. The C75G polymorphism in Obp99a is located in an intron. Frequent local secondary structures predicted for sequences with C and G in position 75 are very similar except for the pairing in position 75, where C is paired and G is not paired (Figure 4). Position 75 in Obp99a pre-mRNA is close to the splice site (35 bp upstream from the exon/intron boundary), and the different predicted RNA structures could potentially affect splicing (SOLNICK 1985).
|
The three common haplotypes formed by SNPs at positions 67 and 78 of Obp99d produce different optimal local secondary structures (Figure 4), although the free energy for the local secondary structures in this region does not differ dramatically (from –133.9 to –126.5 kcal/mol, Figure 4). Structure predictions for Obp99d mRNAs showed that nucleotides in positions 67 (G67A, a nonsynonymous polymorphism) and 78 (T78G, a synonymous polymorphism) are frequently involved in the neighboring stem-loop structures and modulate base pairing of neighboring nucleotides (Figure 4). Nucleotides G and A in position 67 have different potentials for base pairing for the different haplotypes (P = 0.7 for G67/T78, P = 0.61 for G67/G78, and P = 0.05 for A67/78G, respectively, supplemental Table S4 at http://www.genetics.org/supplemental/). The rarest haplotype (A67/T78) has a dramatically different potential for base pairing in position 78 (P = 0.2 for T78, supplemental Table S4) compared to the three more common haplotypes. Since these polymorphisms are in strong LD, we cannot infer whether the effect on olfactory behavior is due to a structural change in the protein or the mRNA.
| DISCUSSION |
|---|
|
|
|---|
Sequence analyses showed that OBPs have different evolutionary histories and statistical tests for deviations from neutrality identify different signatures of selection for eight OBPs. The neutrality tests used detect signatures of selection over different evolutionary time scales (SABETI et al. 2006); this is most likely the reason that inferences regarding selection were not consistent for the different tests. Tajima's D (TAJIMA 1993) and Fu and Li's D* (FU and LI 1993) tests are based on different sensitivities of summary measures of nucleotide diversity within a species to a selective sweep and have different powers to detect departure from neutrality. Fay and Wu's H test (FAY and WU 2000) analyzes high-frequency-derived alleles as a signature of a recent selective sweep and, thus, extends over a recent evolutionary period as high-frequency-derived alleles rapidly reach fixation (SABETI et al. 2006). The results of these tests may be confounded by demographic history. Changes of population size can affect Tajima's D (TAJIMA 1993) and Fu and Li's D* (FU and LI 1993) tests. Populations that have undergone a recent bottleneck in population size have a similar compressed genealogy to populations that are under positive selection; Tajima's D (TAJIMA 1993) and Fu and Li's D* (FU and LI 1993) statistics are expected to be positive in both cases. Negative Tajima's D (TAJIMA 1993) and Fu and Li's D* (FU and LI 1993) statistics are expected for populations undergoing recent expansion of population size, which increases in
w. Population subdivision can confound Fay and Wu's test by generating a nonsignificant H statistic. In contrast, the McDonald–Kreitman test (MCDONALD and KREITMAN 1991) compares the ratios of nonsynonymous and synonymous substitutions within and among species and detects deviations from neutrality that persist over a long evolutionary time. Perhaps the best indicators of the diverse patterns of evolution of Obp genes are the different levels of nucleotide diversity in closely linked genes that presumably experience similar mutation and recombination rates (Figure 1).
Behavioral responses to benzaldehyde showed a greater range of phenotypic variation than observed previously with standard inbred laboratory stocks, and were elicited at a higher concentration of benzaldehyde. Olfactory responses were not correlated with locomotor reactivity scores, indicating that low olfactory-response scores did not result from locomotor impairments due to inbreeding depression. ANOVA showed significant variation in sexual dimorphism but a high genetic correlation between the sexes (rMF = 0.893; Table 3). Variation in sexual dimorphism in avoidance responses to benzaldehyde has also been observed previously in chromosome substitution lines (MACKAY et al. 1996) and among co-isogenic P-element insertion lines that affect olfactory behavior (ANHOLT et al. 1996; SAMBANDAN et al. 2006).
The extensive phenotypic variation in our wild-derived inbred population provided a favorable scenario for association analyses. We detected associations with SNPs in Obp99a, Obp99c, and Obp99d which implicate these OBPs in recognition of benzaldehyde. However, we examined only
25% of the Obp gene family and it is likely that additional OBPs would contribute to phenotypic variation in the response to this odorant. It should be noted that OBPs that interact with odorants, such as benzaldehyde, but for which there is no segregating variation in the population under study, would go undetected by our approach. Furthermore, the detection power for associations depends on the sample size, and additional associations that make a smaller contribution to the observed phenotypic variation might be detected if the population size were expanded. Similarly, our analysis to date has focused on only the commonly used test odorant benzaldehyde. Expanding this analysis to include the entire family of Obp genes with a battery of odorants would enable a comprehensive characterization of ligand specificities of the OBP family. However, our experiments with a single odorant and a limited number of Obp genes already show that the recognition of benzaldehyde by OBPs is redundant and likely combinatorial, reminiscent of odorant recognition by odorant receptors (MALNIC et al. 1999). Functional redundancy may allow the persistence of segregating null alleles in the population, observed by us and others (TAKAHASHI and TAKANO-SHIMIZU 2005).
Previous association studies in Drosophila have implicated SNPs in noncoding regions of Catsup in phenotypic variation in sternopleural bristle number, environmental plasticity of abdominal bristle number, and starvation resistance (CARBONE et al. 2006). SNPs associated with variation in longevity, locomotor behavior, starvation resistance, and bristle number have been identified also in functional regions of the protein (CARBONE et al. 2006). Tests for association of SNPs in Obp genes and responsiveness to benzaldehyde revealed four polymorphisms implicating three OBPs, all within the Obp99 cluster, in the recognition of this odorant. Whereas nonsynonymous SNPs in coding regions of Obp genes can affect ligand binding by introducing variation in protein structure, the most parsimonious explanation for the phenotypic effects of synonymous SNPs and SNPs in regulatory regions would be alterations in mRNA structure. Indeed, structure predictions of mRNAs encoded by alternative haplotypes of Obp99d show that a single base substitution can have a profound effect on secondary mRNA structure (Figure 4), which could affect its transport, splicing, ribosome binding, or translation efficiency (KIMCHI-SARFATY et al. 2007). Our theoretical predictions of causal effects of altered mRNA structures on the behavioral phenotype, however, will need to be supported experimentally in the future.
Elegant electrophysiological studies have generated molecular response profiles of odorant receptors in the Drosophila antennae and maxillary palps and shown that a single odorant can activate multiple odorant receptors (DE BRUYNE et al. 1999; HALLEM et al. 2004; CARLSON and HALLEM 2006). Thus far, the function of OBPs in odorant recognition has remained enigmatic, as there is no clear correlation between expression patterns of OBPs and odorant receptors. The population genetics approach described here represents a first step toward defining molecular recognition profiles of the OBP family. Such information together with the expression patterns of odorant receptors of known response profiles will, ultimately, clarify how these olfactory gene families interact in enabling the fly to sense its chemical environment.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
ANHOLT, R. R. H., R. F. LYMAN and T. F. C. MACKAY, 1996 Effects of single P-element insertions on olfactory behavior in Drosophila melanogaster. Genetics 143: 293–301.[Abstract]
ANHOLT, R. R. H., C. L. DILDA, S. CHANG, J. J. FANARA, N. H. KULKARNI et al., 2003 The genetic architecture of odor-guided behavior in Drosophila: epistasis and the transcriptome. Nat. Genet. 35: 180–184.[CrossRef][Medline]
BENTON, R., S. SACHSE, S. W. MICHNICK and L. B. VOSSHALL, 2006 Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 4: e20.[CrossRef][Medline]
BHALERAO, S., A. SEN, R. STOCKER and V. RODRIGUES, 2003 Olfactory neurons expressing identified receptor genes project to subsets of glomeruli within the antennal lobe of Drosophila melanogaster. J. Neurobiol. 54: 577–592.[CrossRef][Medline]
CARBONE, M. A., K. W. JORDAN, R. F. LYMAN, S. T. HARBISON, J. LEIPS et al., 2006 Phenotypic variation and natural selection at Catsup, a pleiotropic quantitative trait gene in Drosophila. Curr. Biol. 16: 912–919.[CrossRef][Medline]
CARLSON, J. R., and E. A. HALLEM, 2006 Coding of odors by a receptor repertoire. Cell 125: 143–160.[CrossRef][Medline]
CHURCHILL, G. A., and R. W. DOERGE, 1994 Empirical threshold values for quantitative trait mapping. Genetics 138: 963–971.[Abstract]
CLYNE, P. J., C. G. WARR, M. R. FREEMAN, D. LESSING, J. KIM et al., 1999 A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron 22: 327–338.[CrossRef][Medline]
DE BRUYNE, M., P. J. CLYNE and J. R. CARLSON, 1999 Odor coding in a model olfactory organ: the Drosophila maxillary palp. J. Neurosci. 19: 4520–4532.
DEVAUD, J. M., 2003 Experimental studies of adult Drosophila chemosensory behaviour. Behav. Processes 64: 177–196.[CrossRef][Medline]
EDWARDS, A. C., S. M. ROLLMANN, T. J. MORGAN and T. F. C. MACKAY, 2006 Quantitative genomics of aggressive behavior in Drosophila melanogaster. PLoS Genet. 2: e154.[CrossRef][Medline]
FALCONER, D. S., and T. F. C. MACKAY, 1996 Introduction to Quantitative Genetics, Ed. 4. Prentice Hall, Englewood Cliffs, NJ.
FANARA, J. J., K.O. ROBINSON, S.M. ROLLMANN, R. R. H. ANHOLT and T. F. C. MACKAY, 2002 Vanaso is a candidate quantitative trait gene for olfactory behavior in Drosophila melanogaster. Genetics 162: 1321–1328.
FAY, J. C., and C. I. WU, 2000 Hitchhiking under positive Darwinian selection. Genetics 155: 1405–1413.
FU, Y. X., and W. H. LI, 1993 Statistical tests of neutrality of mutations. Genetics 133: 693–709.[Abstract]
GALINDO, K., and D. P. SMITH, 2001 A large family of divergent Drosophila odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics 159: 1059–1072.
GAO, Q., B. YUAN and A. CHESS, 2000 Convergent projections of Drosophila olfactory neurons to specific glomeruli in the antennal lobe. Nat. Neurosci. 3: 780–785.[CrossRef][Medline]
GRAHAM, L. A., and P. L. DAVIES, 2002 The odorant-binding proteins of Drosophila melanogaster: annotation and characterization of a divergent gene family. Gene 292: 43–55.[CrossRef][Medline]
HA, T. S., and D. P. SMITH, 2006 A pheromone receptor mediates 11-cis-vaccenyl acetate-induced responses in Drosophila. J. Neurosci. 26: 8727–8733.
HALLEM, E. A., M. G. HO and J. R. CARLSON, 2004 The molecular basis of odor coding in the Drosophila antenna. Cell 117: 965–979.[CrossRef][Medline]
HEKMAT-SCAFE, D. S., C. R. SCAFE, A. J. MCKINNEYAND and M. A. TANOUYE, 2002 Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 12: 1357–1369.
HARBISON, S.T., S. CHANG, K. P. KAMDAR and T. F. C. MACKAY, 2005 Quantitative genomics of starvation stress resistance in Drosophila. Genome Biol. 6: R36.[CrossRef][Medline]
HUDSON, R. R., 1987 Estimating the recombination parameter of a finite population model without selection. Genet. Res. 50: 245–250.[Medline]
HUDSON, R. R., M. KREITMAN and M. AGUADÉ, 1987 A test of neutral molecular evolution based on nucleotide data. Genetics 116: 153–159.
JORDAN, K. W., M. A. CARBONE, A. YAMAMOTO, T. J. MORGAN and T. F. C. MACKAY, 2007 Quantitative genomics of locomotor behavior in Drosophila melanogaster. Genome Biol. 8: R172.[CrossRef][Medline]
KIM, M. S., A. REPP and D. P. SMITH, 1998 LUSH odorant-binding protein mediates chemosensory responses to alcohols in Drosophila melanogaster. Genetics 150: 711–721.
KIMCHI-SARFATY, C., J. M. OH, I. W. KIM, Z. E. SAUNA, A. M. CALCAGNO et al., 2007 A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science 315: 525–528.
KURTOVIC, A., A. WIDMER and B. J. DICKSON, 2007 A single class of olfactory neurons mediates behavioural responses to a Drosophila sex pheromone. Nature 446: 542–546.[CrossRef][Medline]
LARSSON, M. C., A. I. DOMINGOS, W. D. JONES, M. E. CHIAPPE, H. AMREIN et al., 2004 Or83b encodes a broadly expressed odorant receptor essential for Drosophila olfaction. Neuron 43: 703–714.[CrossRef][Medline]
MACKAY, T. F. C., J. B. HACKETT, R. F. LYMAN, M. L. WAYNE and R. R. H. ANHOLT, 1996 Quantitative genetic variation of odor-guided behavior in a natural population of Drosophila melanogaster. Genetics 144: 727–735.[Abstract]
MACKAY, T. F. C., S. L. HEINSOHN, R. F. LYMAN, A. J. MOEHRING, T. J. MORGAN et al., 2005 Genetics and genomics of Drosophila mating behavior. Proc. Natl. Acad. Sci. USA 102: 6622–6629.
MALNIC, B., J. HIRONO, T. SATO and L. B. BUCK, 1999 Combinatorial receptor codes for odors. Cell 96: 713–723.[CrossRef][Medline]
MATHEWS, D. H., J. SABINA, M. ZUKER and D. H. TURNER, 1999 Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J. Mol. Biol. 288: 911–940.[CrossRef][Medline]
MATSUO, T., S. SUGAYA, J. YASUKAWA, T. AIGAKI and Y. FUYAMA, 2007 Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol. 5: e118.[CrossRef][Medline]
MCDONALD, J. H., and M. KREITMAN, 1991 Adaptive protein evolution at the Adh locus in Drosophila. Nature 351: 652–654.[CrossRef][Medline]
MCGRAW, L. A., G. GIBSON, A. G. CLARK and M. F. WOLFNER, 2004 Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster. Curr. Biol. 14: 1509–1514.[CrossRef][Medline]
MCKENNA, M.P., D. S. HEKMAT-SCAFE, P. GAINES and J. R. CARLSON, 1994 Putative Drosophila pheromone-binding proteins expressed in a subregion of the olfactory system. J. Biol. Chem. 269: 16340–16347.
MOROZOVA, T. V., R. R. H. ANHOLT and T. F. C. MACKAY, 2006 Transcriptional response to alcohol exposure in Drosophila melanogaster. Genome Biol. 7: R95.[CrossRef][Medline]
NACKLEY, A. G., S. A. SHABALINA, I. E. TCHIVILEVA, K. SATTERFIELD, O. KORCHYNSKYI et al., 2006 Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science 314: 1930–1933.
NAZIPOVA, N. N., S. A. SHABALINA, A. Y. OGURTSOV, A. S. KONDRASHOV, M. A. ROYTBERG et al., 1995 SAMSON: a software package for the biopolymer primary structure analysis. Comput. Appl. Biosci. 11: 423–426.
PIKIELNY, C. W., G. HASAN, F. ROUYER and M. ROSBASH, 1994 Members of a family of Drosophila putative odorant-binding proteins are expressed in different subsets of olfactory hairs. Neuron 12: 35–49.[CrossRef][Medline]
ROBERTSON, H. M., C. G. WARR and J. R. CARLSON, 2003 Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 100: 14537–14542.
ROZAS, J., J. C. SáNCHEZ-DELBARRIO, X. MESSEGUERAND and R. ROZAS, 2003 DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496–2497.
SABETI, P. C, S. F. SCHAFFNER, B. FRY, J. LOHMUELLER, P. VARILLY et al., 2006 Positive natural selection in the human lineage. Science 312: 1614–1620.
SAMBANDAN, D., A. YAMAMOTO, J. J. FANARA, T. F. C. MACKAY and R. R. H. ANHOLT, 2006 Dynamic genetic interactions determine odor-guided behavior in Drosophila melanogaster. Genetics 174: 1349–1363.
SHABALINA, S. A., A. Y. OGURTSOV and N. A. SPIRIDONOV, 2006 A periodic pattern of mRNA secondary structure created by the genetic code. Nucleic Acids Res. 34: 2428–2437.
SOLNICK, D., 1985 Alternative splicing caused by RNA secondary structure. Cell 43: 667–676.[CrossRef][Medline]
STOCKINGER, P., D. KVITSIANI, S. ROTKOPF, L. TIRIAN and B. J. DICKSON, 2005 Neural circuitry that governs Drosophila male courtship behavior. Cell 121: 795–807.[CrossRef][Medline]
TAJIMA, F., 1993 Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595.
TAKAHASHI, A., and T. TAKANO-SHIMIZU, 2005 A high-frequency null mutant of an odorant-binding protein gene, Obp57e, in Drosophila melanogaster. Genetics 170: 709–718.
VOSSHALL, L., A. WONG and R. AXEL, 2000 An olfactory sensory map in the fly brain. Cell 102: 147–159.[CrossRef][Medline]
VOSSHALL, L., H. AMREIN, P. MOROZOV, A. RZHETSKY and R. AXEL, 1999 A spatial map of olfactory receptor expression in the Drosophila antenna. Cell 96: 725–736.[CrossRef][Medline]
XU, P., R. ATKINSON, D. N. JONESAND and D. P. SMITH, 2005 Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron 45: 193–200.[CrossRef][Medline]
YOUNG, J. M., C. FRIEDMAN, E. M. WILLIAMS, J. A. ROSS, L. TONNES-PRIDDY et al., 2002 Different evolutionary processes shaped the mouse and human olfactory receptor gene families. Hum. Mol. Genet. 11: 535–546.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |