Evolution of single and double Wolbachia symbioses during speciation in the Drosophila simulans complex.

Maternally inherited bacteria of the genus Wolbachia are responsible for the early death of embryos in crosses between uninfected females and infected males in several insect species. This phenomenon, known as cytoplasmic incompatibility, also occurs between strains infected by different symbionts in some species, including Drosophila simulans. Wolbachia was found in two species closely related to D. simulans, Drosophila mauritiana, and Drosophila sechellia, and shown to cause incompatibility in the latter species but not in D. mauritiana. Comparison of bacterial and mtDNA history clarifies the origins of bacterial and incompatibility polymorphisms in D. simulans. Infection in D. mauritiana is probably the result of introgression of an infected D. simulans cytoplasm. Some D. simulans and D. sechellia cytoplasmic lineages harbor two bacteria as a consequence of a double infection which probably occurred in a common ancestor. The descendant symbionts in each species are associated with similar incompatibility relationships, which suggests that little variation of incompatibility types has occurred within maternal lineages beyond that related to the density of symbionts in their hosts.

The phylogeny of closely related species as revealed by the genealogy of a speciation gene, Odysseus

Molecular differentiation between races or closely related species is often incongruent with the reproductive divergence of the taxa of interest. Shared ancient polymorphism and/or introgression during secondary contact may be responsible for the incongruence. At loci contributing to speciation, these two complications should be minimized (1, 2); hence, their variation may more faithfully reflect the history of the species' reproductive differentiation. In this study, we analyzed DNA polymorphism at the Odysseus (OdsH) locus of hybrid sterility between Drosophila mauritiana and Drosophila simulans and were able to verify such a prediction. Interestingly, DNA variation only a short distance away (1.8 kb) appears not to be influenced by the forces that shape the recent evolution of the OdsH coding region. This locus thus may represent a test case of inferring phylogeny of very closely related species.

Molecular heterochrony in the early development of Drosophila

Heterochrony, the relative change of developmental timing, is one of the major modes of macroevolutionary change; it identifies temporally disassociated units of developmental evolution. Here, we report the results of a fine-scale temporal study for the expression of the developmental gene hairy and morphological development in three species of Drosophila, D. melanogaster, D. simulans, and D. pseudoobscura. The results suggest that between and among closely related species, temporal displacement of ontogenetic trajectory is detected even at the earliest stage of development. Overall, D. simulans shows the earliest expression, followed by D. melanogaster, and then by D. pseudoobscura. Setting D. melanogaster as the standard, we find the approximate time to full expression is accelerated by 13 min, 48 s in D. simulans and retarded by 24 min in D. pseudoobscura. Morphologically, again with D. melanogaster setting the standard, initiation of cellularization is faster in D. simulans by 15 min, 42 s; and initiation of morphogenesis is faster in D. simulans by 18 min, 7 s. These results seem to be consistent with the finding that the approximate time to full expression of hairy is accelerated by 13 min, 48 s in D. simulans. On the other hand, the same morphological events are delayed by 5 min, 32 s, and by 11 min, 32 s, respectively, in D. pseudoobscura. These delays are small, compared with the 24-min delay in full expression. The timing changes, in total, seem consistent with continuous phyletic evolution of temporal trajectories. Finally, we speculate that epigenetic interactions of hairy expression timing and cell-cycle timing may have led to morphological differences in the terminal system of the larvae.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Evolutionary EST analysis identifies rapidly evolving male reproductive proteins in Drosophila

Sequence comparisons of genomes or expressed sequence tags (ESTs) from related organisms provide insight into functional conservation and diversification. We compare the sequences of ESTs from the male accessory gland of Drosophila simulans to their orthologs in its close relative Drosophila melanogaster, and demonstrate rapid divergence of many of these reproductive genes. Nineteen (∼11%) of 176 independent genes identified in the EST screen contain protein-coding regions with an excess of nonsynonymous over synonymous changes, suggesting that their divergence has been accelerated by positive Darwinian selection. Genes that encode putative accessory gland-specific seminal fluid proteins had a significantly elevated level of nonsynonymous substitution relative to nonaccessory gland-specific genes. With the 57 new accessory gland genes reported here, we predict that ∼90% of the male accessory gland genes have been identified. The evolutionary EST approach applied here to identify putative targets of adaptive evolution is readily applicable to other tissues and organisms.

Naturally occurring variation in copia expression is due to both element (cis) and host (trans) regulatory variation.

Significant differences in levels of copia [Drosophila long terminal repeat (LTR) retrotransposon] expression exist among six species representing the Drosophila melanogaster species complex (D. melanogaster, Drosophila mauritiana, Drosophila simulans, Drosophila sechellia, Drosophila yakuba, and Drosophila erecta) and a more distantly related species (Drosophila willistoni). These differences in expression are correlated with major size variation mapping to putative regulatory regions of the copia 5' LTR and adjacent untranslated leader region (ULR). Sequence analysis indicates that these size variants were derived from a series of regional duplication events. The ability of the copia LTR-ULR size variants to drive expression of a bacterial chloramphenicol acetyltransferase reporter gene was tested in each of the seven species. The results indicate that both element-encoded (cis) and host-genome-encoded (trans) genetic differences are responsible for the variability in copia expression within and between Drosophila species. This finding indicates that models purporting to explain the dynamics and distribution of retrotransposons in natural populations must consider the potential impact of both element-encoded and host-genome-encoded regulatory variation to be valid. We propose that interelement selection among retrotransposons may provide a molecular drive mechanism for the evolution of eukaryotic enhancers which can be subsequently distributed throughout the genome by retrotransposition.

Evolutionary changes in the expression pattern of a developmentally essential gene in three Drosophila species.

The hypothesis that morphological evolution may largely result from changes in gene regulation rather than gene structure has been difficult to test. Morphological differences among insects are often apparent in the cuticle structures produced. The dopa decarboxylase (Ddc) and alpha-methyldopa hypersensitive (amd) genes arose from an ancient gene duplication. In Drosophila, they have evolved nonoverlapping functions, including the production of distinct types of cuticle, and for Ddc, the production of the neurotransmitters, dopamine and serotonin. The amd gene is particularly active in the production of specialized flexible cuticles in the developing embryo. We have compared the pattern of amd expression in three Drosophila species. Several regions of expression conserved in all three species but, surprisingly, a unique domain of expression is found in Drosophila simulans that does occur in the closely related (2-5 million years) Drosophila melanogaster or in the more remote (40-50 million years) Drosophila virilis. The "sudden" appearance of a completely new and robust domain of expression provides a glimpse of evolutionary variation resulting from changes in regulation of structural gene expression.

Localization of pheromonal sexual dimorphism in Drosophila melanogaster and its effect on sexual isolation.

Drosophila melanogaster is sexually dimorphic for cuticular hydrocarbons, with males and females having strikingly different profiles of the long-chain compounds that act as contact pheromones. Gas-chromatographic analysis of sexual mosaics reveals that the sex specificity of hydrocarbons is located in the abdomen. This explains previous observations that D. melanogaster males display the strongest courtship toward mosaics with female abdomens. We also show that males of the sibling species Drosophila simulans preferentially court D. melanogaster mosaics with male abdomens. Because the primary male hydrocarbon in D. melanogaster is also the primary female hydrocarbon in D. simulans, this supports the idea that interspecific differences in cuticular hydrocarbons contribute to sexual isolation.