Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

Divergence and heterogeneity of the histone gene repeating units in the Drosophila melanogaster species subgroup

The repeating units of the histone gene cluster containing the H1, H2A, H2B and H4 genes were amplified by PCR from the Drosophila melanogaster species subgroup, i.e., D. yakuba, D. erecta, D. sechellia, D. mauritiana, D. teissieri and D. orena. The PCR products were cloned and their nucleotide sequences of about 4.6-4.8kbp were determined to elucidate the mechanism of molecular evolution of the histone gene family. The heterogeneity among the histone gene repeating units was 0.6% and 0.7% for D. yakuba and D. sechellia, respectively, indicating the same level of heterogeneity as in the H3 gene region of D. melanogaster. Divergence of the genes among species even in the most closely related ones was much greater than the heterogeneity among family members, indicating a concerted mode of evolution for the histone gene repeating units. Among the species in the D. melanogaster species subgroup, the histone gene regions as well as 3rd codon position of the coding region showed nearly the same GC contents. These results suggested that the previous conclusion on analysis of the H3 gene regions, the gene family evolution in a concerted fashion, holds true for the whole histone gene repeating unit.     

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.     

Divergence of the Yellow Gene between Drosophila Melanogaster and D. Subobscura: Recombination Rate, Codon Bias and Synonymous Substitutions

The yellow (y) gene maps near the telomere of the X chromosome in Drosophila melanogaster but not in D. subobscura. Thus the strong reduction in the recombination rate associated with telomeric regions is not expected in D. subobscura. To study the divergence of a gene whose recombination rate differs between two species, the y gene of D. subobscura was sequenced. Sequence comparison between D. melanogaster and D. subobscura revealed several elements conserved in noncoding regions that may correspond to putative cis-acting regulatory sequences. Divergence in the y gene coding region between D. subobscura and D. melanogaster was compared with that found in other genes sequenced in both species. Both, yellow and scute exhibit an unusually high number of synonymous substitutions per site (p(s)). Also for these genes, the extent of codon bias differs between both species, being much higher in D. subobscura than in D. melanogaster. This pattern of divergence is consistent with the hitchhiking and background selection models that predict an increase in the fixation rate of slightly deleterious mutations and a decrease in the rate of fixation of slightly advantageous mutations in regions with low recombination rates such as in the y-sc gene region of D. melanogaster.     

Intraspecific nuclear DNA variation in Drosophila

We have summarized and analyzed all available nuclear DNA sequence polymorphism studies for three species of Drosophila, D. melanogaster (24 loci), D. simulans (12 loci), and D. pseudoobscura (5 loci). Our major findings are: (1) The average nucleotide heterozygosity ranges from about 0.4% to 2% depending upon species and function of the region, i.e., coding or noncoding. (2) Compared to D. simulans and D. pseudoobscura (which are about equally variable), D. melanogaster displays a low degree of DNA polymorphism. (3) Noncoding introns and 3' and 5' flanking DNA shows less polymorphism than silent sites within coding DNA. (4) X-linked genes are less variable than autosomal genes. (5) Transition (Ts) and transversion (Tv) polymorphisms are about equally frequent in non-coding DNA and at fourfold degenerate sites in coding DNA while Ts polymorphisms outnumber Tv polymorphisms by about 2:1 in total coding DNA. The increased Ts polymorphism in coding regions is likely due to the structure of the genetic code: silent changes are more often Ts's than are replacement substitutions. (6) The proportion of replacement polymorphisms is significantly higher in D. melanogaster than in D. simulans. (7) The level of variation in coding DNA and the adjacent noncoding DNA is significantly correlated indicating regional effects, most notably recombination. (8) Surprisingly, the level of polymorphism at silent coding sites in D. melanogaster is positively correlated with degree of codon usage bias. (9) Three proposed tests of the neutral theory of DNA polymorphisms have been performed on the data: Tajima's test, the HKA test, and the McDonald-Kreitman test. About half of the loci fail to conform to the expectations of neutral theory by one of the tests. We conclude that many variables are affecting levels of DNA polymorphism in Drosophila, from properties of nucleotides to population history and, perhaps, mating structure. No simple, all encompassing explanation satisfactorily accounts for the data.      

Molecular evolutionary analysis of a histone gene repeating unit from Drosophila simulans

A repeating unit of the histone gene cluster from Drosophila simulans containing the H1, H2A, H2B and H4 genes (the H3 gene region has already been analyzed) was cloned and analyzed. A nucleotide sequence of about 4.6 kbp was determined to study the nucleotide divergence and molecular evolution of the histone gene cluster. Comparison of the structure and nucleotide sequence with those of Drosophila melanogaster showed that the four histone genes were located at identical positions and in the same directions. The proportion of different nucleotide sites was 6.3% in total. The amino acid sequence of H1 was divergent, with a 5.1% difference. However, no amino acid change has been observed for the other three histone proteins. Analysis of the GC contents and the base substitution patterns in the two lineages, D. melanogaster and D. simulans, with a common ancestor showed the following. 1) A strong negative correlation was found between the GC content and the nucleotide divergence in the whole repeating unit. 2) The mode of molecular evolution previously found for the H3 gene was also observed for the whole repeating unit of histone genes; the nucleotide substitutions were stationary in the 3' and spacer regions, and there was a directional change of the codon usage to the AT-rich codons. 3) No distinct difference in the mode or pattern of molecular evolution was detected for the histone gene repeating unit in the D. melanogaster and D. simulans lineages. These results suggest that selectional pressure for the coding regions of histones, which eliminate A and T, is less effective in the D. melanogaster and D. simulans lineages than in the other GC-rich species.     

Excess Polymorphism at the Adh Locus in DROSOPHILA MELANOGASTER

The evolutionary history of a region of DNA encompassing the Adh locus is studied by comparing patterns of variation in Drosophila melanogaster and its sibling species, D. simulans. An unexpectedly high level of silent polymorphism in the Adh coding region relative to the 5' and 3' flanking regions in D. melanogaster is revealed by a populational survey of restriction polymorphism using a four-cutter filter hybridization technique as well as by direct sequence comparisons. In both of these studies, a region of the Adh gene encompassing the three coding exons exhibits a frequency of polymorphism equal to that of a 4-kb 5' flanking region. In contrast, an interspecific sequence comparison shows a two-fold higher level of divergence in the 5' flanking sequence compared to the structural locus. Analysis of the patterns of variation suggest an excess of polymorphism within the D. melanogaster Adh locus, rather than lack of polymorphism in the 5' flanking region. An approach is outlined for testing neutral theory predictions about patterns of variation within and between species. This approach indicates that the observed patterns of variation are incompatible with an infinite site neutral model.     

Patterns of natural selection at the Alcohol dehydrogenase gene of Drosophila americana

Similar outcomes are often observed in species exposed to similar selective regimes, but it is unclear how often the same mechanism of adaptive evolution is followed. Here we present an analysis of selection affecting sequence variation in the Alcohol dehydrogenase (Adh) gene of Drosophila americana, a species endemic to a large climate range that has been colonized by D. melanogaster. Unlike D. melanogaster, there is no evidence of selection on allozymes of ADH across the sampled range. This indicates that if there has been a similar adaptive response to climate in D. americana, it is not within the coding region of Adh. Instead, analyses of a combined dataset containing 86 alleles of Adh reveal purifying selection on the Adh gene, especially within its intron sequences. Frequency spectra of derived unpreferred variants at synonymous sites indicate that these sites are affected by weak purifying selection, but the deviation from neutrality is less drastic than observed for derived variants in noncoding introns. This contrast further supports the notion that noncoding sites in Drosophila are often subject to stronger selection pressures than synonymous sites.     

Regulation of histone H4 Lys16 acetylation by predicted alternative secondary structures in roX noncoding RNAs

Despite differences in size and sequence, the two noncoding roX1 and roX2 RNAs are functionally redundant for dosage compensation of the Drosophila melanogaster male X chromosome. Consistent with functional conservation, we found that roX RNAs of distant Drosophila species could complement D. melanogaster roX mutants despite low homology. Deletion of a conserved predicted stem-loop structure in roX2, containing a short GUb (GUUNUACG box) in its 3' stem, resulted in a defect in histone H4K16 acetylation on the X chromosome in spite of apparently normal localization of the MSL complex. Two copies of the GUb sequence, newly termed the "roX box," were functionally redundant in roX2, as mutants in a single roX box had no phenotype, but double mutants showed reduced H4K16 acetylation. Interestingly, mutation of two of three roX boxes in the 3' end of roX1 RNA also reduced H4K16 acetylation. Finally, fusion of roX1 sequences containing a roX box restored function to a roX2 deletion RNA lacking its cognate roX box. These results support a model in which the functional redundancy between roX1 and roX2 RNAs is based, at least in part, on short GUUNUACG sequences that regulate the activity of the MSL complex.     

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.     

Sites of P element insertion and structures of P element deletions in the 5'' region of Drosophila melanogaster RpII215.

Several P element insertion and deletion mutations near the 5' end of Drosophila melanogaster RpII215 have been examined by nucleotide sequencing. Two different sites of P element insertion, approximately 90 nucleotides apart, have been detected in this region of the gene. Therefore, including an additional site of P element insertion within the coding region, there are at least three distinct sites of P element insertion at RpII215. Both 5' sites are within a noncoding portion of transcribed sequences. The sequences of four revertants of one P element insertion mutation (D50) indicate that the P element is either precisely deleted or internally deleted to restore RpII215 activity. Partial internal deletions of the P element result in different RpII215 activity levels, which appear to depend on the specific sequences that remain after excision.     

Molecular Population Genetics of the Distal Portion of the X Chromosome in Drosophila: Evidence for Genetic Hitchhiking of the Yellow-Achaete Region

We have estimated DNA sequence variation and differentiation within and between Drosophila melanogaster and its sibling species, Drosophila simulans, using six-cutter restriction site variation at yellow-achaete (y-ac), phosphogluconate dehydrogenase (Pgd), and period (per). These three gene regions are of varying distance from the telomere of the X chromosome and range from very low to moderate rates of recombination in D. melanogaster. According to Tajima's test of neutrality, the Pgd region has been influenced by balancing selection in D. melanogaster. This is consistent with previous data suggesting the allozyme polymorphism at this locus is visible to selection. The Hudson, Kreitman, Aguadé test of neutrality reveals a significant departure from neutrality for the y-ac region compared to the per or rosy regions in D. simulans. There is also a significant departure for the y-ac region compared to the Adh 5' flanking region in D. melanogaster. In both species the departure appears to be due to reduced variation at y-ac compared to that expected from divergence between D. simulans and D. melanogaster. We conclude that recent hitchhiking associated with the selective fixation of one or more advantageous mutants in the y-ac region is the best explanation for reduced variation at y-ac.     

Nucleotide Variation and Conservation at the Dpp Locus, a Gene Controlling Early Development in Drosophila

A study of polymorphism and species divergence of the dpp gene of Drosophila has been made. Eighteen lines from a population of D. melanogaster were sequenced for 5200 bp of the Hin region of the gene, coding for the dpp polypeptide. A comparison was made with sequence from D. simulans. Ninety-six silent polymorphisms and three amino acid replacement polymorphisms were found. The overall silent polymorphism (0.0247) is low, but haplotype diversity (0.0066 for effectively silent sites and 0.0054 for all sites) is in the range found for enzyme loci. Amino acid variation is absent in the N-terminal signal peptide, the C-terminal TGF-β peptide and in the N-terminal half of the pro-protein region. At the nucleotide level there is strong conservation in the middle half of the large intron and in the 3' untranslated sequence of the last exon. The 3' untranslated conservation, which is perfect for 110 bp among all the divergent species, is unexplained. There is strong positive linkage disequilibrium among polymorphic sites, with stretches of apparent gene conversion among originally divergent sequences. The population apparently is a migration mixture of divergent clades.