Evolution of the olfactory code in the Drosophila melanogaster subgroup

The Drosophila melanogaster subgroup has been the focus of numerous studies about evolution. We address the question of how the olfactory code has evolved among the nine sister species. By using in vivo electrophysiological measurements, so called single-cell recordings, we have established the ligand affinity of a defined subset of olfactory receptor neurons (ORNs) across all nine species. We show that the olfactory code as relayed by the investigated subset of ORNs is conserved to a striking degree. Distinct shifts in the code have occurred only within the simulans clade. However, these shifts are restricted to an altered tuning profile of the same single ORN type in all three of the simulans siblings and a more drastic change unique to D. sechellia, involving a complete loss of one sensillum type in favour of another. The alterations observed in D. sechellia may represent a novel host-specific adaptation to its sole host, morinda fruit (Morinda citrifolia). The overall high degree of similarity of the code within the subgroup is intriguing when considering the great variety in distributions as well as in habitat and host choice of the siblings, factors that could greatly affect the olfactory system.     

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.     

Behavioral analyses of mutants for two odorant-binding protein genes, Obp57d and Obp57e, in Drosophila melanogaster

The odorant-binding protein (OBP) functions in chemosensation in insects. Two OBP genes, Obp57d and Obp57e, are involved in the evolution of the unique host-plant preference in Drosophila sechellia. Comparative analysis of the Obp57d/e genomic sequences in the Drosophila melanogaster species group has revealed that the rapid evolution of Obp57d and Obp57e has resulted in functional divergence between the two genes. Here, using D. melanogaster knockout strains generated by gene targeting, we examined the roles of Obp57d and Obp57e in behavioral response to a series of fatty acids. In the taste-based oviposition-site preference assay, the knockout flies showed stronger preference for acids than wild-type flies, indicating that the normal functions of Obp57d and Obp57e are associated with the suppression of positive preference for C6-C9 acids. Heterozygotes for each knockout haplotype also exhibited a significantly different behavioral response compared with wild-type flies, suggesting that Obp57d and Obp57e have a gene dosage effect on behavior. In contrast, the wild-type and knockout flies exhibited similar responses in the feeding assay and the odor-based free-walking assay, suggesting that the two OBPs' contribution to feeding and olfactory behaviors is small. Taken together, our results demonstrated that each of Obp57d and Obp57e in D. melanogaster contributes to the determination of reproductive sites, suggesting that the two OBP genes play an important role in ecological adaptation of Drosophila.     

Patterns of Microsatellite Variability in the Drosophila Melanogaster Complex

Forty-seven microsatellite loci were amplified in Drosophila melanogaster, Drosophila simulans, Drosophila mauritiana and Drosophila sechellia. The two cosmopolitan species D. melanogaster and D. simulans were found to be the most variable ones, followed by D. mauritiana and D. sechellia. A model based clustering algorithm was applied to the population samples of D. melanogaster, D. simulans and D. sechellia. No evidence for population substructure was detected within species--most likely due to insufficient power. A Markov chain Monte Carlo method developed for demographic inference based on microsatellites provided unambiguous evidence for population contraction in D. melanogaster, D. simulans and D. sechellia, despite that the D. melanogaster and D. simulans population samples were of non-African origin and represented recently expanded populations.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

Mate choice in Drosophila melanogaster and D. sechellia: criteria and their variation depending on courtship song

Courtship song is one of the most important male signals in Drosophila mating. A female D. melanogaster copulates more readily when given a conspecific type of courtship song. Female D. melanogaster accepted winged D. sechellia more than wingless ones in a no choice test. Copulation frequencies varied between strains and were significantly correlated with male mating propensity of D. sechellia. Females from three of five strains of D. sechellia accepted winged D. melanogaster less than wingless ones in no choice tests, suggesting that D. melanogaster songs reduce mating in D. sechellia females. Multiple choice tests showed that males prefer conspecifics. In female choice tests, D. melanogaster females copulated with winged heterospecific D. sechellia males more than with wingless conspecific D. melanogaster males in the confined condition, suggesting that song differences between D. melanogaster and D. sechellia affect D. melanogaster females less than the absence of courtship song. Copyright 2000 The Association for the Study of Animal Behaviour.     

Genetic variability and robustness of host odor preference in Drosophila melanogaster

Chemosensory stimuli play a crucial role for host selection in insects, including the fruit fly Drosophila melanogaster. Drosophila has been instrumental in unraveling the neurological basis of olfactory processing in insects. Basic knowledge regarding chemical ecology and thorough studies of olfactory preferences are still lacking to a great extent in D. melanogaster, however. We have characterized repeatable variation in olfactory preference between five classical D. melanogaster wild-type strains toward a large array of natural host odors and synthetic compounds. By recording the rate of attraction over up to 24 hr, we could compare stimuli varying in attractiveness and characterize phenotypic parameters on the basis of individual stimuli and the whole stimulus array. Behavioral differences between strains were predominantly due to variation in a single phenotypic parameter: their overall responsiveness toward optimal and suboptimal olfactory stimuli. These differences were not explained by variation in olfactory sensitivity, locomotory activity, or general vigor monitored by survival. Comparisons with three recently established wild-type strains indicated that a high behavioral threshold against accepting suboptimal olfactory stimuli is the characteristic phenotype of wild D. melanogaster.     

DNA Sequence Variation at the Period Locus within and among Species of the Drosophila Melanogaster Complex

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced ``sister' taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.     

Chromosome localization of the lambda20p1.4 clone of the Drosophila melanogaster nuclear lamina DNA in the melanogaster species subgroup of the genus Drosophila (Sophophora)

Chromosome localization of sequences homologous to lambda 20p1.4 of the Drosophila melanogaster nuclear lamina DNA (nlDNA) was established by in situ hybridization in species of the melanogaster subgroup. DNA of the lambda 20p1.4 clone was shown to be located in the chromocenter in all the species examined. Laboratory strains of D. simulans, D. mauritiana, and D. sechellia exhibited interspecific differences in localization of lambda 20p1.4 nlDNA on chromosome arms. In eight natural populations, intraspecific polymorphism of lambda 20p1.4 nlDNA chromosome localization was shown to be present in D. simulans but absent in D. melanogaster. The possible participation of transposable elements in nlDNA relocation is discussed.     

Genetics of egg production in Drosophila sechellia

Drosophila sechellia, an island endemic that specializes on a single host plant, has a lower rate of egg production than its generalist sister species D. melanogaster, D. simulans, and D. mauritiana. Earlier work showed that part of this difference in egg production was due to a reduction in the number of ovarioles in D. sechellia relative to its sister species. Here, I extend this earlier work by genetically analyzing the difference in egg production between D. sechellia and D. simulans. In all, 10 genetic markers were used in several interspecific backcrosses to identify chromosome regions that affected the rate of egg production. While previously mapped factors affecting ovariole number appear to impact the rate of egg production, new, non-ovariole factors were also identified. Overall, the difference in egg production between D. sechellia and D. simulans appears to be a polygenic trait. The relationship between these factors and genes involved the adaptation of D. sechellia to its host plant is not yet clear. The data are consistent with the hypothesis that decline in egg production is, in part, a negative pleiotropic effect of genetic changes required for host specialization in D. sechellia, although finer-scale genetic analysis of both traits is needed to truly test this hypothesis.     

Patterns of mutation and selection at synonymous sites in Drosophila

That natural selection affects molecular evolution at synonymous sites in protein-coding sequences is well established and is thought to predominantly reflect selection for translational efficiency/accuracy mediated through codon bias. However, a recently developed maximum likelihood framework, when applied to 18 coding sequences in 3 species of Drosophila, confirmed an earlier report that the Notch gene in Drosophila melanogaster was evolving under selection in favor of those codons defined as unpreferred in this species. This finding opened the possibility that synonymous sites may be subject to a variety of selective pressures beyond weak selection for increased frequencies of the codons currently defined as "preferred" in D. melanogaster. To further explore patterns of synonymous site evolution in Drosophila in a lineage-specific manner, we expanded the application of the maximum likelihood framework to 8,452 protein coding sequences with well-defined orthology in D. melanogaster, Drosophila sechellia, and Drosophila yakuba. Our analyses reveal intragenomic and interspecific variation in mutational patterns as well as in patterns and intensity of selection on synonymous sites. In D. melanogaster, our results provide little statistical evidence for recent selection on synonymous sites, and Notch remains an outlier. In contrast, in D. sechellia our findings provide evidence in support of selection predominantly in favor of preferred codons. However, there is a small subset of genes in this species that appear to be evolving under selection in favor of unpreferred codons, which indicates that selection on synonymous sites is not limited to the preferential fixation of mutations that enhance the speed or accuracy of translation in this species.     

Highly tissue specific expression of Sphinx supports its male courtship related role in Drosophila melanogaster

Sphinx is a lineage-specific non-coding RNA gene involved in regulating courtship behavior in Drosophila melanogaster. The 5' flanking region of the gene is conserved across Drosophila species, with the proximal 300 bp being conserved out to D. virilis and a further 600 bp region being conserved amongst the melanogaster subgroup (D. melanogaster, D. simulans, D. sechellia, D. yakuba, and D. erecta). Using a green fluorescence protein transformation system, we demonstrated that a 253 bp region of the highly conserved segment was sufficient to drive sphinx expression in male accessory gland. GFP signals were also observed in brain, wing hairs and leg bristles. An additional ～800 bp upstream region was able to enhance expression specifically in proboscis, suggesting the existence of enhancer elements. Using anti-GFP staining, we identified putative sphinx expression signal in the brain antennal lobe and inner antennocerebral tract, suggesting that sphinx might be involved in olfactory neuron mediated regulation of male courtship behavior. Whole genome expression profiling of the sphinx knockout mutation identified significant up-regulated gene categories related to accessory gland protein function and odor perception, suggesting sphinx might be a negative regulator of its target genes.     

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.     

A New Hybrid Rescue Allele in Drosophila melanogaster

Crosses of Drosophila melanogaster females to males of its sibling species Drosophila simulans, Drosophila mauritiana and Drosophila sechellia produce no sons and daughters that are viable only at low temperatures. We describe here a novel rescue allele Df(1)EP307-1-2 isolated on the basis of its suppression of high temperature hybrid female lethality. Df(1)EP307-1-2 also rescues hybrid males to the pharate adult stage, the same stage at which it is lethal to D. melanogaster pure species males. Molecular analysis indicates that Df(1)EP307-1-2 is associated with a deletion of about 61 kb in the 9D region of the X chromosome. The structure of Df(1)EP307-1-2 suggests that it was formed by a process similar to P-element induced male recombination.     

The genetic basis of larval resistance to a host plant toxin in Drosophila sechellia

The larvae of Drosophila sechellia are highly resistant to octanoic acid, a toxin found in D. sechellia's host plant, Morinda citrifolia. In contrast, close relatives of D. sechellia, D. simulans and D. melanogaster, are not resistant. In a series of interspecific backcrosses, 11 genetic markers were used to map factors affecting egg-to-adult ('larval') resistance in D. sechellia. The third chromosome harbours at least one partially dominant resistance factor. The second chromosome carries at least two mostly dominant resistance factors but no recessive factors. However, neither the X chromosome--which contains 20% of D. sechellia's genome--nor the fourth chromosome appear to affect resistance. These data suggest that larval resistance to Morinda toxin may involve only a handful of genes. These results, when compared with a previous analysis of adult resistance to Morinda toxin in D. sechellia, suggest that larval resistance may involve a subset of the genes underlying adult resistance.     

Variation in fiber number of a male-specific muscle between Drosophila species: a genetic and developmental analysis

We characterize a newly discovered morphological difference between species of the Drosophila melanogaster subgroup. The muscle of Lawrence (MOL) contains about four to five fibers in D. melanogaster and Drosophila simulans and six to seven fibers in Drosophila mauritiana and Drosophila sechellia. The same number of nuclei per fiber is present in these species but their total number of MOL nuclei differs. This suggests that the number of muscle precursor cells has changed during evolution. Our comparison of MOL development indicates that the species difference appears during metamorphosis. We mapped the quantitative trait loci responsible for the change in muscle fiber number between D. sechellia and D. simulans to two genomic regions on chromosome 2. Our data eliminate the possibility of evolving mutations in the fruitless gene and suggest that a change in the twist might be partly responsible for this evolutionary change.     

Wolbachia Infection and Cytoplasmic Incompatibility in Drosophila Species

Forty-one stocks from 30 Drosophila species were surveyed for Wolbachia infection using PCR technology. D. sechellia and two strains of D. auraria were found to be infected and were tested for the expression of cytoplasmic incompatibility, along with D. ananassae and D. melanogaster strains, which are already known to be infected. D. ananassae and D. melanogaster show levels of incompatibility up to 25%, while D. auraria and D. sechellia exhibit levels of egg mortality ~60%. A dot-blot assay using the dnaA sequence as probe was developed to assess the infection levels in individual males that were used in incompatibility crosses. A positive correlation between bacterial density and cytoplasmic incompatibility was observed. The stocks examined can be clustered into at least two groups, depending on the levels of infection relative to the degree of cytoplasmic incompatibility exhibited. One group, containing D. simulans Hawaii, D. sechellia, and D. auraria, exhibits high levels of cytoplasmic incompatibility relative to levels of infection; all the other species and D. simulans Riverside exhibit significantly lower levels of cytoplasmic incompatibility relative to levels of infection. These data show that, in addition to bacterial density, bacterial and/or host factors also affect the expression of cytoplasmic incompatibility.     

Genetics of Differences in Pheromonal Hydrocarbons between Drosophila Melanogaster and D. Simulans

Females of Drosophila melanogaster and its sibling species D. simulans have very different cuticular hydrocarbons, with the former bearing predominantly 7,11-heptacosadiene and the latter 7-tricosene. This difference contributes to reproductive isolation between the species. Genetic analysis shows that this difference maps to only the third chromosome, with the other three chromosomes having no apparent effect. The D. simulans alleles on the left arm of chromosome 3 are largely recessive, allowing us to search for the relevant regions using D. melanogaster deficiencies. At least four nonoverlapping regions of this arm have large effects on the hydrocarbon profile, implying that several genes on this arm are responsible for the species difference. Because the right arm of chromosome 3 also affects the hydrocarbon profile, a minimum of five genes appear to be involved. The large effect of the thrid chromosome on hydrocarbons has also been reported in the hybridization between D. simulans and its closer relative D. sechellia, implying either an evolutionary convergence or the retention in D. sechellia of an ancestral sexual dimorphism.     

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1- kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58- 0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.