Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.     

DNA variability and divergence at the notch locus in Drosophila melanogaster and D. simulans: a case of accelerated synonymous site divergence

DNA diversity in two segments of the Notch locus was surveyed in four populations of Drosophila melanogaster and two of D. simulans. In both species we observed evidence of non-steady-state evolution. In D. simulans we observed a significant excess of intermediate frequency variants in a non-African population. In D. melanogaster we observed a disparity between levels of sequence polymorphism and divergence between one of the Notch regions sequenced and other neutral X chromosome loci. The striking feature of the data is the high level of synonymous site divergence at Notch, which is the highest reported to date. To more thoroughly investigate the pattern of synonymous site evolution between these species, we developed a method for calibrating preferred, unpreferred, and equal synonymous substitutions by the effective (potential) number of such changes. In D. simulans, we find that preferred changes per "site" are evolving significantly faster than unpreferred changes at Notch. In contrast we observe a significantly faster per site substitution rate of unpreferred changes in D. melanogaster at this locus. These results suggest that positive selection, and not simply relaxation of constraint on codon bias, has contributed to the higher levels of unpreferred divergence along the D. melanogaster lineage at Notch.     

Rates of DNA sequence evolution are not sex-biased in Drosophila melanogaster and D. simulans

To determine whether male- or female-biased mutation rates have affected the molecular evolution of Drosophila melanogaster and D. simulans, we calculated the male-to-female ratio of germline cell divisions ([symbol: see text]) from germline generation data and the male-to-female ratio of mutation rate ([symbol: see text]) by comparing chromosomal levels of nucleotide divergence. We found that the ratio of germline cell divisions changes from indicating a weak female bias to indicating a weak male bias as the age of reproduction increases. The range of [symbol: see text] values that we observed, however, does not lead us to expect much, if any, difference in mutation rate between the sexes. Silent and intron nucleotide divergence were compared between nine loci on the X chromosome and nine loci on the second and third chromosomes. The average levels of nucleotide divergence were not significantly different across the chromosomes, although both silent and intron sites show a trend toward slightly more divergence on the X. These results indicate a lack of sex- or chromosome-biased molecular evolution in D. melanogaster and D. simulans.      

Comparative studies of allozyme loci in Drosophila simulans and Drosophila melanogaster. I. Three dipeptidase loci

Genetic variation at three dipeptidase loci (Dip-A, Dip-B, and Dip-C) in Drosophila simulans was analyzed by starch gel electrophoresis. Dip-A was found to be polymorphic in four populations, while Dip-B and Dip-C were found to be polymorphic in one. The numbers of different alleles found at each respective locus were: Dip-A, two; Dip-B, two; and Dip-C, three. Dip-A was genetically mapped at 57.9 on the second chromosome, and Dip-B and Dip-C at 80.9 and 87.9 on the third chromosome, respectively. Neither Dip-B nor Dip-C has been mapped in D. melanogaster because both loci are apparently monomorphic. Their map positions in D. simulans with respect to flanking markers whose homologous genes have been cytogenetically localized in D. melanogaster suggested that they might be mapped cytogenetically by using available deficiencies in D. melanogaster. Accordingly, by the construction of interspecific hybrids which carried deficiencies of melanogaster and an allele of simulans with a mobility different from that of the fixed melanogaster allele, Dip-B and Dip-C were localized between 87F12-14 and 88C1-3 and between 87B5-6 and 87B8-10, respectively, in the salivary gland chromosomes of D. melanogaster. The similarity between these two species is discussed on the basis of these findings.     

hobo transposable elements in Drosophila melanogaster and D. simulans.

Genomic patterns of occurrence of the transposable element hobo are polymorphic in the sibling species Drosophila melanogaster and D. simulans. Most tested strains of both species have apparently complete (3.0 kb) and smaller hobo elements (H lines), but in both species some strains completely lack such canonical hobo elements (E lines). The occurrence of H and E lines in D. simulans as well as in D. melanogaster implies that an hypothesis of recent introduction in the latter species is inadequate to explain the phylogenetic occurrence of hobo. Particular internally deleted elements, the approximately 1.5 kb Th1 and Th2 elements, are abundant in many lines of D. melanogaster, and an analogous 1.1 kb internally deleted element, h del sim, is abundant in most lines of D. simulans. Besides the canonical hobo sequences, both species (and their sibling species D. sechellia and D. mauritiana) have many hobo-hybridizing sequences per genome that do not appear to be closely related to the canonical hobo sequence.     

Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans.

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.     

Nucleotide sequence divergence in the A+T-rich region of mitochondrial DNA in Drosophila simulans and Drosophila mauritiana

We have determined the nucleotide sequences of two regions within the A+T-rich region of mitochondrial DNA (mtDNA) in the siIII type of Drosophila simulans and the maI type of D. mauritiana. The sequences of the two regions in siIII and maI are almost identical. The sequences include elements corresponding to the type I and type II repeats elements and the T-stretches as reported in D. melanogaster; an approximately 340-bp region (A region) adjacent to the tRNA(Ile) gene includes a part of the type II repeat element, and an approximately 440- bp region (B region) includes a central portion of the A+T-rich region between the type I and type II repeat arrays. Each sequence of the two species was compared with those of D. melanogaster and D. yakuba. The sequences of the A region are relatively well conserved among the four species. The alignment of the two sequences of the B region with those of D. melanogaster and D. yakuba requires numerous insertions/deletions. For both regions, nucleotide differences between D. simulans or D. mauritiana and D. melanogaster are similar to those between the two and D. yakuba. The tendency is obvious in a subregion within the type II repeat element in the A region. These findings suggest that the rate of nucleotide substitution in the subregion is accelerated in the lineage leading to D. melanogaster. Loss of functional constraint in the stem-loop-forming sequence is proposed for this acceleration.      

High level of divergence of male-reproductive-tract proteins, between Drosophila melanogaster and its sibling species, D. simulans

We compared male-reproductive-tract polypeptides of Drosophila melanogaster and D. simulans by using two-dimensional gel electrophoresis. Approximately 64% of male-reproductive-tract polypeptides were identical between two randomly chosen isofemale lines from these two species, compared with 83% identity for third-instar imaginal wing-disc polypeptides. Qualitatively similar differences were found between reproductive tracts and imaginal discs when D. sechellia was compared with D. melanogaster and with D. simulans. When genic polymorphism was taken into account, approximately 10% of male- reproductive-tract polypeptides were apparently fixed for different alleles between D. melanogaster and D. simulans; this proportion is the same as that found for soluble enzymes by one-dimensional gel electrophoresis. Strikingly, approximately 20% of male-reproductive- tract polypeptides of either D. melanogaster or D. simulans had no detectable homologue in the other species. We propose that proteins of the Drosophila male reproductive tract may have diverged more extensively between species than have other types of proteins and that much of this divergence may involve large changes in levels of polypeptide expression.      

Comparative Life Histories and Ecophysiology of Drosophila melanogaster and D. simulans

Numerous laboratory investigations have compared Drosophila melanogaster and D. simulans for various life history traits and fitness related ecophysiological parameters. From presently available information, it is however difficult to get a general comparative pattern describing the divergence of their ecological niches and understanding their demographic success. Two environmental factors seem however to have played a major role: temperature and alcoholic resources. From an ecophysiological approach, D. simulans may be described as generally more sensitive to stresses; other results point to this species as more cold adapted than its sibling; in some cases, however, D. simulans may appear as better adapted to a warm environment. When investigated, ecophysiological traits show a lesser geographic variability in D. simulans than in D. melanogaster. Presently available information does not explain the ecological prevalence of D. simulans in many places with a mild temperate or subtropical climate. This is presumably due to the fact that most comparisons have been done at a single, standard temperature of 25 degrees C. Comparative studies should be undertaken, spanning the thermal ranges of the two species, and the phenotypic plasticity of ecophysiological traits should now be considered.     

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.      

Effects of female behavior on wing display and courtship of male Drosophila simulans and Drosophila melanogaster

Corrolations between female rejection behaviors and male wing display were calculated for both Drosophila simulans and Drosophila melanogaster intraspicific pair-matings. No significant correlations were found for D. melanogaster, but in D. simulans flicking by the female appeared to be associated with a shift in male wing display pattern resulting in higher levels of vibration. Flicking did not appear to discourage courtship by males in either species.     

Comparison of alcohol dehydrogenase expression in Drosophila melanogaster and D. simulans

Alcohol dehydrogenase (ADH) gene expression was analyzed in Drosophila melanogaster and its sibling species D. simulans. The levels of ADH activity, ADH-cross-reacting material (CRM), and ADH-mRNA were analyzed for several strains of each species, which derive from diverse geographic locations around the world. There is considerable quantitative variation in ADH activity, CRM level, and RNA level among strains within species at all developmental stages. However, the only consistent differences between the two species are in pupal RNA level and in late-adult activity and CRM level. Late-adult melanogaster flies that are homozygous for the Slow allozyme have approximately twice the level of ADH activity and CRM as do simulans flies. The regression of activity on CRM over strains is highly significant and essentially the same for each species, which means that most, if not all, of the activity difference between the species is due to a difference in concentration of the ADH protein. In contrast, there is no significant regression of CRM level on mRNA level in adults of either species; nor is there a significant difference in RNA level between species. Therefore, the difference in ADH protein concentration is not due to RNA template availability. Thus, the interspecific difference in ADH level in adults must be due either to a difference in the rate of translation of the two RNAs or to a difference in protein stability.     

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

Lack of nucleotide polymorphism in the Y-linked sperm flagellar dynein gene Dhc-Yh3 of Drosophila melanogaster and D. simulans

We studied levels of intra- and interspecific nucleotide variation associated with a Y-linked gene in five members of the Drosophila melanogaster subgroup. Using published sequence for 348 bp of the Dhc-Yh3 gene, and degenerate PCR primers designed from comparisons of the sea urchin and Chlamydomonas flagellar dynein genes, we recovered a 1738-bp region in D. melanogaster. Analyses of sequence variation in a worldwide collection of 11 lines of D. melanogaster and 10 lines of D. simulans found only a single silent polymorphism in the latter species. The synonymous site divergence per site for Dhc-Yh3 is comparable to values for X and autosomal genes. Assuming a Wright-Fisher population model, the lack of variation is statistically less than expected using appropriately reduced estimates of theta from the X and autosomes. Because the Y chromosome encodes only six known genes, genetic hitchhiking associated with background selection is unlikely to explain this low variation. Conversely, adaptive hitchhiking, as associated with sex-ratio chromosomes, or a large variance in male fertility may reduce the polymorphism on the Y chromosome. Codon bias is very low, as seen for other genes in regions of low recombination.     

Low genic variation in male-reproductive-tract proteins of Drosophila melanogaster and D. simulans

We report results, using two-dimensional gel electrophoresis (2DE), of natural population surveys of allelic variation in approximately 300 male-reproductive-tract polypeptides in both Drosophila melanogaster and its sibling species, D. simulans. Despite our efforts to maximize operational sensitivity of our 2DE gels to polymorphism, variation estimates in both species were low (proportion of polymorphic loci [P] = 9%, and average heterozygosity [H] = 1%-3%), compared with those by one-dimensional gel electrophoresis (1DE) (P = 29%-55%; H = 8%-19%) in the same populations. However, H of polymorphic loci was very similar for 2DE and 1DE proteins; and for 17 of a total of 54 polymorphic proteins, 2DE detected three or four distinct alleles. The results suggest that the differing levels of variability widely seen with 1DE and 2DE are real and reflect differing intensities of functional constraint between different classes of structural loci. However, the alternative possibility remains that 2DE has a greater between-locus unevenness of variant detection sensitivity than does 1DE.