Polymorphism and divergence at a Drosophila pseudogene locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcp psi) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcp psi fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcp psi revealed lower levels of polymorphism (theta = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcp psi within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcp psi lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

Evolution of the mitochondrial ATPase 6 gene in Drosophila: unusually high level of polymorphism in D. melanogaster

We have determined 1990 bp mitochondrial DNA sequence which extends from 3' end of the cytochrome oxidase subunit I (COI) gene to 5' end of the COIII gene from two sibling species of Drosophila, D. simulans and D. mauritiana. Analyses of the sequences and part of the NADH dehydrogenase subunit 2 gene and the COI gene together with those from D. melanogaster and D. yakuba revealed that amino-acid substitution rate of the ATPase 6 gene seems to be higher in some strains of D. melanogaster than in the other species. High level of amino-acid polymorphism in this gene was observed in D. melanogaster. Synonymous substitution rate is relatively constant in all the genes examined, suggesting that mutation rate is not higher in the ATPase 6 gene of D. melanogaster. The amino-acid substitutions found specifically in D. melanogaster are at the sites which are not conserved among mammals, yeast and E. coli. These sites of the ATPase 6 gene might lose the selective constraint in D. melanogaster, and the amino-acid substitutions can be explained by neutral mutations and random genetic drift.     

Evolutionary inferences from DNA variation at the 6-phosphogluconate dehydrogenase locus in natural populations of drosophila: selection and geographic differentiation

Several allozyme-coding genes in Drosophila melanogaster show patterns suggesting that polymorphisms at these loci are targets of balancing selection. An important question is whether these genes have similar distributions of underlying DNA sequence variation which would indicate similar evolutionary processes occurring in this class of loci. One such locus, 6-phosphogluconate dehydrogenase (Pgd), has previously been shown to exhibit clinal variation for Fast/Slow electromorph variation in the United States and Australia, unusually large electromorph frequency differences between the United States and Africa, and other patterns indicative of selection. We measured four-cutter DNA restriction site and allozyme variation at Pgd among 142 D. melanogaster X chromosomes collected from several geographic regions including North Carolina, California, and Zimbabwe (Africa). We also sequenced a representative sample of 13 D. melanogaster Pgd genes collected in North Carolina and a single copy of Pgd from the sibling species, Drosophila simulans. While some population genetic models predict excess DNA polymorphism in genes which are targets of balancing selection, the D. melanogaster samples from the United States had significantly reduced levels of DNA polymorphism and extraordinarily high levels of linkage disequilibrium, providing evidence of hitchhiking effects of advantageous mutants at Pgd or at linked sites. Therefore, while selection has probably influenced the distribution of DNA variation at Pgd, the precise nature of these selective events remains obscure. Since the Pgd region appears to have low rates of crossing over, the reduced level of variation at this locus supports the idea that recombination rates are important determinants of levels of DNA polymorphism in natural populations. Furthermore, while patterns of allozyme variation are very similar at Pgd and Adh, the DNA data show that the evolutionary histories of these genes are dramatically different. We observed extensive differences in the amount and distribution of variation in D. melanogaster Pgd samples from the United States and Zimbabwe which cannot be explained by differential selection on the Fast/Slow polymorphism in these two geographic regions. Thus, genetic drift among partially isolated populations has also been an important factor in determining the distribution of variation at Pgd in D. melanogaster. Finally, we assayed four-cutter variation at Pgd in a sample of 19 D. simulans X chromosomes and observed reduced levels of DNA variability and high levels of linkage disequilibrium. These patterns are consistent with predictions of some hitchhiking models.     

Contrasting patterns of nucleotide sequence variation at the glucose dehydrogenase (Gld) locus in different populations of Drosophila melanogaster

We have analyzed nucleotide sequence variation at the Glucose dehydrogenase (Gld) locus from four populations of Drosophila melanogaster from four continents. All four population samples show a significant reduction in silent variation compared to the neutral expectation. The levels of silent variation across all four populations are consistent with the predictions of the background selection model; however, Zimbabwe has a remarkably low level of variation. In the face of dramatically reduced silent polymorphism, an amino acid variant, leading to the common allozyme polymorphism at Gld, remains in low to intermediate frequency in all non-African samples. In the Chinese population sample, the ratio of replacement to silent variation is significantly elevated compared to the neutral expectation. The difference in patterns of variation across these population samples suggests that selection on Gld (or the Gld region) has been different in the Chinese population than in the other three.     

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.     

Loss of notum macrochaetae as an interspecific hybrid anomaly between Drosophila melanogaster and D. simulans

With the aim of revealing genetic variation accumulated among closely related species during the course of evolution, this study focuses on loss of macrochaetae on the notum as one of the developmental anomalies seen in interspecific hybrids between Drosophila melanogaster and its closely related species. Interspecific hybrids between a line of D. melanogaster and D. simulans isofemale lines exhibited a wide range in the number of missing bristles. By contrast, D. mauritiana and D. sechellia lines showed almost no reduction in bristle number in hybrids with D. melanogaster. Genetic analysis showed that the D. simulans X chromosome confers a large effect on hybrid bristle loss, although X-autosome interaction may be involved. This suggests that at least one genetic factor contributing to hybrid anomalies arose recently on a D. simulans X chromosome. Moreover, the results indicate sex dependency: the male hybrids were more susceptible to bristle loss than the female hybrids were. Use of cell type markers suggests that the defect does not lie in cell fate decisions during bristle development, but in the maintenance of neural fate and/or differentiation of the descendants of sensory mother cells.     

Conserved regions of the timeless (tim) clock gene in Drosophila analyzed through phylogenetic and functional studies

Circadian (approximately 24-hr) rhythms in Drosophila melanogaster depend upon cyclic expression of the period (per) and timeless (tim) genes, which encode interacting components of the endogenous clock. The per gene has been isolated from other insects and, more recently, a per ortholog was found in mammals where its expression oscillates in a circadian fashion. We report here the complete sequence of a tim gene from another species, Drosophila virilis. TIM is better conserved than the PER protein is between these two species (76 vs. 54% overall amino acid identity), and putative functional domains, such as the PER interaction domains and the nuclear localization signal, are highly conserved. The acidic domain and the cytoplasmic localization domain, however, are within the least conserved regions. In addition, the initiating methionine in the D. virilis gene lies downstream of the proposed translation start for the original D. melanogaster tim cDNA and corresponds to the one used by D. simulans and D. yakuba. Among the most conserved parts of TIM is a region of unknown function near the N terminus. We show here that deletion of a 32 amino acid segment within this region affects rescue of rhythms in arrhythmic tim01 flies. Flies carrying a full-length tim transgene displayed rhythms with approximately 24-hr periods, indicating that a fully functional clock can be restored in tim01 flies through expression of a tim transgene. Deletion of the segment mentioned above resulted in very long activity rhythms with periods ranging from 30.5 to 48 hr.     

Isolation of protease-free alcohol dehydrogenase (ADH) from Drosophila simulans and several homozygous and heterozygous Drosophila melanogaster variants

The enzyme alcohol dehydrogenase (ADH) from several naturally occurring ADH variants of Drosophila melanogaster and Drosophila simulans was isolated. Affinity chromatography with the ligand Cibacron Blue and elution with NAD+ showed similar behavior for D. melanogaster ADH-FF, ADH-71k, and D. simulans ADH. Introduction of a second Cibacron Blue affinity chromatography step, with gradient elution with NAD+, resulted in pure and stable enzymes. D. melanogaster ADH-SS cannot be eluted from the affinity chromatography column at a high concentration of NAD+ and required a pH gradient for its purification, preceded by a wash step with a high concentration of NAD+. Hybrid Drosophila melanogaster alcohol dehydrogenase FS has been isolated from heterozygous flies, using affinity chromatography with first elution at a high concentration NAD+, directly followed by affinity chromatography elution with a pH gradient. Incubation of equal amounts of pure homodimers of Drosophila melanogaster ADH-FF and ADH-SS, in the presence of 3 M urea at pH 8.6, for 30 min at room temperature, followed by reassociation yielded active Drosophila melanogaster ADH-FS heterodimers. No proteolytic degradation was found after incubation of purified enzyme preparations in the absence or presence of SDS, except for some degradation of ADH-SS after very long incubation times. The thermostabilities of D. melanogaster ADH-71k and ADH-SS were almost identical and were higher than those of D. melanogaster ADH-FF and D. simulans ADH. The thermostability of D. melanogaster ADH-FS was lower than those of D. melanogaster ADH-FF and ADH-SS. D. melanogaster ADH-FF and ADH-71k have identical inhibition constants with the ligand Cibacron Blue at pH 8.6, which are two times higher at pH 9.5. The Ki values for D. simulans ADH are three times lower at both pH values. D. melanogaster ADH-SS and ADH-FS have similar Ki values, which are lower than those for D. melanogaster ADH-FF at pH 8.6. But at pH 9.5 the Ki value for ADH-FS is the same as at pH 8.6, while that of ADH-SS is seven times higher. Kinetic parameters of Drosophila melanogaster ADH-FF, ADH-SS, and ADH-71k and Drosophila simulans ADH, at pH 8.6 and 9.5, showed little or no variation in K(m)eth values. The K(m)NAD values measured at pH 9.5 for Drosophila alcohol dehydrogenases are all lower than those measured at pH 8.6. The rate constants (kcat) determined for all four Drosophila alcohol dehydrogenases are higher at pH 9.5 than at pH 8.6. D. melanogaster ADH-FS showed nonlinear kinetics.     

Genetics of food preference in Drosophila sechellia. I. Responses to food attractants

To reveal the genetic mechanism of host selection in a monophagous fruit fly Drosophila sechellia, olfactory responses and oviposition preferences of this species were compared with those of closely related polyphagous species, D. simulans and D. melanogaster. Adult flies of D. sechellia were strongly attracted to the ripe fruit of Morinda citrifolia which is known to be the sole breeding site of this species. They were also attracted to the odor of n-caproic acid which is contained in the ripe fruit of M. citrifolia and is presumably responsible for the characteristic odor of the fruit. In contrast, D. simulans and D. melanogaster showed a strong repulsion to n-caproic acid. In parallel with the olfactory responses, D. sechellia females laid eggs preferentially on a medium containing n-caproic acid, to which the other two species showed an aversion. Genetic analyses using the hybrid progeny between D. sechellia and D. simulans suggested that the species differences in these behaviors are controlled by gene(s) located on the second chromosome.     

Reproductive isolation and morphogenetic evolution in Drosophila analyzed by breakage of ethological barriers

This article reports the breaking of ethological barriers through the constitution of soma-germ line chimeras between species of the melanogaster subgroup of Drosophila, which are ethologically isolated. Female Drosophila yakuba and D. teissieri germ cells in a D. melanogaster ovary produced functional oocytes that, when fertilized by D. melanogaster sperm, gave rise to sterile yakuba-melanogaster and teissieri-melanogaster male and female hybrids. However, the erecta-melanogaster and orena-melanogaster hybrids were lethal, since female D. erecta and D. orena germ cells in a D. melanogaster ovary failed to form oocytes with the capacity to develop normally. This failure appears to be caused by an altered interaction between the melanogaster soma and the erecta and orena germ lines. Germ cells of D. teissieri and D. orena in a D. melanogaster testis produced motile sperm that was not stored in D. melanogaster females. This might be due to incompatibility between the teissieri and orena sperm and the melanogaster seminal fluid. A morphological analysis of the terminalia of yakuba-melanogaster and teissieri-melanogaster hybrids was performed. The effect on the terminalia of teissieri-melanogaster hybrids of a mutation in doublesex, a regulatory gene that controls the development of the terminalia, was also investigated.     

Survey of malathion resistance and avermectin susceptibility in field populations of Drosophila melanogaster (Diptera: Drosophilidae) and D. simulans

Populations of Drosophila melanogaster (Meigen) and Drosophila simulans (Sturtevant) from various geographic locations in North America and elsewhere were sampled to assess the distribution of malathion resistance and avermectin tolerance. Comparisons are made to previous reports of a differential response to malathion selection in these species. Results indicate that for malathion, resistance is widespread in D. simulans but varies considerably in D. melanogaster. For avermectins, although the pattern of tolerance is similar in both species, there exists a significant amount of variation in these levels between geographic regions. We consider how these different geographic patterns of resistance distribution may shed insight on the relative roles of population structure and exposure history in affecting the spread of resistance. In addition, we consider further the use of D. melanogaster and D. simulans as a model for examining the effects of insecticide exposure on sympatric populations.     

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.     

Hybrid lethal systems in the Drosophila melanogaster species complex

Lethal phases of the hybrids between Drosophila melanogaster and its sibling species, D. simulans are classified into three types: (1) embryonic lethality in hybrids carrying D. simulans cytoplasm and D. melanogaster X chromosome, (2) larval lethality in hybrids not carrying D. simulans X, and (3) temperature-sensitive pupal lethality in hybrids carrying D. simulans X. The same lethal phases are also observed when either of the two other sibling species, D. mauritiana or D. sechellia, is employed for hybridization with D. melanogaster. Here, we describe genetic analyses of each hybrid lethality, and demonstrate that these three types of lethality are independent phenomena. We then propose two models to interpret the mechanisms of each hybrid lethality. The first model is a modification of the conventional X/autosome imbalance hypothesis assuming a lethal gene and a suppressor gene are involved in the larval lethality, while the second model is for embryonic lethality assuming an interaction between a maternal-effect lethal gene and a suppressor gene.     

Functional analysis of eve stripe 2 enhancer evolution in Drosophila: rules governing conservation and change

Experimental investigations of eukaryotic enhancers suggest that multiple binding sites and trans-acting regulatory factors are often required for wild-type enhancer function. Genetic analysis of the stripe 2 enhancer of even-skipped (eve), an important developmental gene in Drosophila, provides support for this view. Given the importance of even-skipped expression in early Drosophila development, it might be predicted that many structural features of the stripe 2 enhancer will be evolutionarily conserved, including the DNA sequences of protein binding sites and the spacing between them. To test this hypothesis, we compared sequences of the stripe 2 enhancer between four species of Drosophila: D. melanogaster, D. yakuba, D. erecta and D. pseudoobscura. Our analysis revealed a large number of nucleotide substitutions in regulatory protein binding sites for bicoid, hunchback, Kruppel and giant, as well as a systematic change in the size of the enhancer. Some of the binding sites in D. melanogaster are either absent or modified in other species. One functionally important bicoid-binding site in D. melanogaster appears to be recently evolved. We, therefore, investigated possible functional consequences of sequence differences among these stripe 2 enhancers by P-element-mediated transformation. This analysis revealed that the eve stripe 2 enhancer from each of the four species drove reporter gene expression at the identical time and location in D. melanogaster embryos. Double staining of native eve protein and transgene mRNA in early embryos showed that the reporter gene mimicked native eve expression and, in every case, produced sharply defined stripes at the blastoderm stage that were coincident with eve stripe 2 protein. We argue that stripe 2 eve expression in Drosophila evolution can be viewed as being under constant stabilizing selection with respect to the location of the anterior and posterior borders of the stripe. We further hypothesize that the stripe 2 enhancer is functionally robust, so that its evolution may be governed by the fixation of both slightly deleterious and adaptive mutations in regulatory protein binding sites as well as in the spacing between binding sites. This view allows for a slow but continual turnover of functionally important changes in the stripe 2 enhancer.     

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.     

Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex

The Acp26Aa and Acp26Ab genes that code for male accessory gland proteins are tandemly arranged in the species of the Drosophila melanogaster complex. An approximately 1.6-kb region encompassing both genes has been sequenced in 10, 24, and 18 lines from Spain, Ivory Coast, and Malawi, respectively; the previously studied 10 lines from North Carolina have also been included in the analyses. A total of 110 nucleotide and 4 length polymorphisms were detected. Silent variation for the whole Acp26A region was slightly higher in African than in non-African populations, while for both genes nonsynonymous variation was similar in all populations studied. Based on Fst estimates no major genetic differentiation was detected between East and West Africa, while in general non-African populations were strongly differentiated from both African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites revealed that directional selection acting on amino acid replacement changes has driven the evolution of the Acp26Aa protein in the last 2.5 myr.     

Mutation and evolution of microsatellites in Drosophila melanogaster

Levels of nucleotide polymorphism in the Drosophila melanogaster genome are correlated with rates of recombination. This relationship may be due to hitchhiking of advantageous mutations (selective sweeps) or to continual removal of deleterious mutations from the genome (background selection). One test of the relative contributions of selective sweeps and background selection to the observed levels of variation in the genome of D. melanogaster is to compare levels of nucleotide variability (with a mutation rate on the order of 10(-9) per nucleotide per generation) with more rapidly evolving DNA loci such as microsatellites. This test depends critically on details of the mutational process of microsatellites. In this paper, we summarize our studies of microsatellite characteristics and mutation rates in D. melanogaster. We find that D. melanogaster microsatellites are short and have a mutation rate (6.5 x 10(-6) per locus per generation) several orders of magnitude lower than mammals studied to date. We further show that genetic variation at 18 dinucleotide repeat microsatellites in a population of D. melanogaster from Maryland is correlated with regional rates of recombination. These and other microsatellite data suggest that both background selection and selective sweeps may contribute to the correlation between DNA sequence variation and recombination in Drosophila.     

Distribution of transposable elements in neotropical species of Drosophila

The phylogenetic distribution of transposable families, P, gypsy, hobo, I, and mariner has been analyzed in 33 species of 11 groups of neotropical Drosophila and a Drosophilidae species Zygotrica vittimaculosa, using squash blot and dot blot. Genomic DNA of almost all neotropical species tested hybridized with gypsy probe and some species showed a particularly strong hybridization signal, as D. gaucha, D. virilis, and species of flavopilosa group. The hobo element was restricted to melanogaster group and some strains of D. willistoni. Only D. simulans DNA showed hybridization to mariner probe in all species tested and D. simulans and D. melanogaster showed hybridization with I element probe. P element homologous sequence was present in D. melanogaster and all species and strains of the willistoni and saltans groups tested. The presence of at least one P-homologous sequence was detected in Drosophila mediopunctata. This one was the only P-bearing species of all six tested from the tripunctata group. Four different pairs of primers homologous to segments of the canonical sequence of D. melanogaster's P were used to amplify specific sequences from D. mediopunctata DNA, showing the occurrence of seemingly well-conserved P-homologous sequences.     

Intra- and interspecies variation among Bari-1 elements of the melanogaster species group

We have investigated the distribution of sequences homologous to Bari-1, a Tc1-like transposable element first identified in Drosophila melanogaster, in 87 species of the Drosophila genus. We have also isolated and sequenced Bari-1 homologues from D. simulans, D. mauritiana, and D. sechellia, the species constituting with D. melanogaster the melanogaster complex, and from D. diplacantha and D. erecta, two phylogenetically more distant species of the melanogaster group. Within the melanogaster complex the Bari-1 elements are extremely similar to each other, showing nucleotide identity values of at least 99.3%. In contrast, Bari-1-like elements from D. diplacantha and D. erecta are on average only 70% similar to D. melanogaster Bari-1 and are usually defective due to nucleotide deletions and/or insertions in the ORFs encoding their transposases. In D. erecta the defective copies are all located in the chromocenter and on chromosome 4. Surprisingly, while D. melanogaster Bari-1 elements possess 26-bp inverted terminal repeats, their D. diplacantha and D. erecta homologues possess long inverted terminal repeats similar to the terminal structures observed in the S elements of D. melanogaster and in several other Tc1-like elements of different organisms. This finding, together with the nucleotide and amino acid identity level between D. diplacantha and D. erecta elements and Bari-1 of D. melanogaster, suggests a common evolutionary origin and a rapid diversification of the termini of these Drosophila Tc1-like elements.