Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

Molecular evolution of the Sex-Ratio inversion complex in Drosophila pseudoobscura: analysis of the Esterase-5 gene region

The Sex-Ratio chromosome in Drosophila pseudoobscura is subject to meiotic drive. It is associated with a series of three nonoverlapping paracentric inversions on the right arm of the X chromosome. The esterase-5 gene region has been localized to section 23 within the subbasal inversion of the Sex-Ratio inversion complex, making esterase- 5 a convenient locus for molecular evolutionary analyses of the Sex- Ratio inversion complex and the associated drive system. A 504-bp fragment of noncoding, intergenic DNA from the esterase-5 gene region was amplified and sequenced from 14 Sex-Ratio and 14 Standard X chromosomes of D. pseudoobscura, and from 9 X chromosomes of its two sibling species, Drosophila persimilis and Drosophila miranda. There is extensive sequence differentiation between the Sex-Ratio and Standard chromosomal types. The common Standard chromosome is highly polymorphic, while, as expected from either the neutral mutation theory or the selective sweep hypothesis, the rarer Sex-Ratio chromosome has much less within-chromosome nucleotide polymorphism. We estimate that the Standard and Sex-Ratio chromosomes in D. pseudoobscura diverged between 700,000 and 1.3 Mya, or at least 2 million generations ago. The clustering of D. pseudoobscura Sex-Ratio chromosomes in a neighbor- joining phylogeny indicates a fairly old, monophyletic origin in this species. It appears from these data that Sex-Ratio genes were present prior to the divergence of D. pseudoobscura and D. persimilis and that both the Standard and Sex-Ratio chromosomes of D. persimilis were derived from the Standard chromosome of D. pseudoobscura after the inversion events that isolated the D. pseudoobscura Sex-Ratio chromosome.      

Nucleotide polymorphism at the xanthine dehydrogenase locus in Drosophila pseudoobscura.

Sequential polyacrylamide electrophoresis has revealed 20 allozymes of xanthine dehydrogenase (XDH) in Drosophila pseudoobscura. DNA sequence determination of seven isolates of the Xdh locus that represent six allozyme classes are presented here. Of the 5,456 sites examined, 180 are polymorphic, with 27 polymorphisms occurring at nonsynonymous, or replacement, sites. An average of nine amino acids differ between XDH allozyme classes, with 85% of the polymorphic amino acids singly represented. The level and pattern of variation observed at Xdh argue that the effective population size of the species is quite large--i.e., on the order of 2 x 10(6)--and that the populations sampled are quite ancient. In addition, as judged by two statistical tests, the levels of nucleotide polymorphism observed at Xdh are compatible with predictions from the neutral theory of molecular evolution.     

Molecular population genetics of Xdh and the evolution of base composition in Drosophila

Few loci have been measured for DNA polymorphism and divergence in several species. Here we report such data from the protein-coding region of xanthine dehydrogenase (Xdh) in 22 species of Drosophila. Many of our samples were from closely related species, allowing us to confidently assign substitutions to individual lineages. Surprisingly, Xdh appears to be fixing more A/T mutations than G/C mutations in most lineages, leading to evolution of higher A/T content in the recent past. We found no compelling evidence for selection on protein variation, though some aspects of the data support the notion that a significant fraction of amino acid polymorphisms are slightly deleterious. Finally, we found no convincing evidence that levels of silent heterozygosity are associated with rates of protein evolution.     

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.      

Molecular evolution of X-linked accessory gland proteins in Drosophila pseudoobscura

In Drosophila melanogaster and Drosophila simulans, positive Darwinian selection drives high rates of evolution of male reproductive genes, and accessory gland proteins (Acps) in particular. Here, we tested whether 13 X-linked male-specific genes, 4 Acps and 9 non-Acps, are under selective forces in the Drosophila pseudoobscura species group, much as those in the D. melanogaster group. We observed a statistically significant correlation in relative rates of nonsynonymous evolution between the two species groups tested. One Acp examined had a higher rate of nonsynonymous substitution than predicted by a neutral model in both species groups, suggesting its divergence was driven by positive Darwinian selection. To further test for the signature of selection, we examined polymorphism of three Acps within D. pseudoobscura. From this test, no Acp individually bore the signature of positive selection, but the 3 Acps together possessed an excess of nonsynonymous differences between species, relative to polymorphism within species. We conclude that faster evolution of Acps in the D. pseudoobscura group appears to be driven by positive selection, as previously suggested in the D. melanogaster group.     

Sexual selection and experimental evolution of chemical signals in Drosophila pseudoobscura

Our expectations for the evolution of chemical signals in response to sexual selection are uncertain. How are chemical signals elaborated? Does sexual selection result in complexity of the composition or in altered quantities of expression? We addressed this in Drosophila pseudoobscura by examining male and female cuticular hydrocarbons (CHs) after 82 generations of elevated (E) sexual selection or relaxed sexual selection through monogamy (M). The CH profile consisted of 18 different components. We extracted three eigenvectors using principal component analysis that explained 72% of the variation. principal component (PC)1 described the amount of CHs produced, PC2 the trade‐off between short‐ and long‐chain CHs and PC3 the trade‐off between apparently arbitrary CHs. In both sexes, the amount of CHs produced was greater in flies from the E treatment. PC3 was also higher, indicating that sexual selection also influenced the evolution of CH composition. The sexes differed in all three PCs, indicating substantial sexual dimorphism in this species, although the magnitude of this dimorphism was not increased as a result of our experimental evolution. Collectively, our work provides direct evidence that sexual selection plays an important role in the evolution of CHs in D. pseudoobscura and that both increased quantity and overall composition are targeted.     

Nucleotide sequence analysis of the delta beta-globin gene region in humans

The continuous DNA sequence of a 16.5-kilobase pair region encompassing the linked delta beta-globin gene cluster in humans is presented with a detailed restriction endonuclease map. There are 38 differences (0.5%) in comparison with published sequence data, corrected for errors in sequencing, resulting in polymorphic rates of 0.2% in exons and 0.76% in 5'-gene flanking regions. Fifteen changes result in the generation or elimination of restriction sites which may be useful in linkage disequilibrium studies. Two pairs of inverted Alu repeats, a pyrimidine-rich region 5' to delta, and (TG)n, (Pu/Py)n, and (ATTTT)n tracts 5' to beta are described. Dinucleotide frequencies and deviation from expected values approximated those found in total human genomic DNA. Regions of less than 50% A + T content were found associated with Alu sequences, a 150-base pair region immediately 5' to the beta gene, exon regions from both genes, and an area 3' to the beta gene. These regions also contained significantly lower than expected CpG levels compared to other regions, suggesting a possible relationship between DNA organizational patterns and functionally important regions. In addition, strand asymmetries in base composition in this region differ from those associated with the fetal globin genes.     

Nucleotide sequence and analysis of the Vibrio alginolyticus sucrose uptake-encoding region

The nucleotide sequence of the Vibrio alginolyticus sucrose uptake-encoding region was determined, and contained two genes, scrA and scrK. The scrA gene encodes an enzyme IISucrose (EIIScr) protein of the phosphoenolpyruvate dependent phosphotransferase system and the scrK gene encodes a fructokinase. The deduced amino acid (aa) sequence for the V. alginolyticus EIIScr protein was homologous with the EIIScr proteins from Streptococcus mutans, Salmonella typhimurium (pUR400 system) and Bacillus subtilis. The deduced aa sequence for the V. alginolyticus fructokinase was homologous with the Escherichia coli enzymes, 6-phosphofructokinase (isoenzyme 2) and ribokinase. Transposon phoA mutagenesis experiments indicated that the EIIScr protein was a membrane-bound protein with a region that extended into the periplasm.     

Sequence and evolution of the Drosophila pseudoobscura glycerol-3-phosphate dehydrogenase locus

The Gpdh genomic region has been cloned and sequenced in Drosophila pseudoobscura. A total of 6.8 kb of sequence was obtained, encompassing all eight exons of the gene. The exons have been aligned with the sequence from D. melanogaster, and the rates of synonymous and nonsynonymous substitution have been compared to those of other genes sequenced in these two species. Gpdh has the lowest rate of nonsynonymous substitution yet seen in genes sequenced in both D. pseudoobscura and D. melanogaster. No insertion/deletion events were observed, and the overall architecture of the gene (i.e., intron sites, etc.) is conserved. An interesting amino acid reversal was noted between the D. melanogaster Fast allele and the D. pseudoobscura gene.     

No accelerated rate of protein evolution in male-biased Drosophila pseudoobscura genes

Sexually dimorphic traits are often subject to diversifying selection. Genes with a male-biased gene expression also are probably affected by sexual selection and have a high rate of protein evolution. We used SAGE to measure sex-biased gene expression in Drosophila pseudoobscura. Consistent with previous results from D. melanogaster, a larger number of genes were male biased (402 genes) than female biased (138 genes). About 34% of the genes changed the sex-related expression pattern between D. melanogaster and D. pseudoobscura. Combining gene expression with protein divergence between both species, we observed a striking difference in the rate of evolution for genes with a male-biased gene expression in one species only. Contrary to expectations, D. pseudoobscura genes in this category showed no accelerated rate of protein evolution, while D. melanogaster genes did. If sexual selection is driving molecular evolution of male-biased genes, our data imply a radically different selection regime in D. pseudoobscura.     

Positive Selection on Nucleotide Substitutions and Indels in Accessory Gland Proteins of the Drosophila pseudoobscura Subgroup

Genes encoding reproductive proteins often diverge rapidly due to positive selection on nucleotide substitutions. While this general pattern is well established, the extent to which specific reproductive genes experience similar selection in different clades has been little explored, nor have possible targets of positive selection other than nucleotide substitutions, such as indels, received much attention. Here, we inspect for the signature of positive selection in the genes encoding five accessory gland proteins (Acps) (Acp26Aa, Acp32CD, Acp53Ea, Acp62F, and Acp70A) originally described from Drosophila melanogaster but with recognizable orthologues in the D. pseudoobscura subgroup. We compare patterns of selection within the D. psuedoobscura subgroup to those in the D. melanogaster subgroup. Similar patterns of positive selection were found in Acp26Aa and Acp62F in the two subgroups, while Acp53Ea and Acp70A experienced purifying selection in both subgroups. These proteins have thus remained targets for similar types of selection over long (>21-MY) periods of time. We also found several indel substitutions and polymorphisms in Acp26Aa and Acp32CD. These indels occur in the same regions as positively selected nucleotide substitutions for Acp26Aa in the D. pseudoobscura subgroup but not in the D. melanogaster subgroup. Rates of indel substitution within Acp26Aa in the D. pseudoobscura subgroup were up to several times those in noncoding regions of the Drosophila genome. This suggests that indel substitutions may be under positive selection and may play a key role in the divergence of some Acps. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Willis Swanson]     

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.     

Nucleotide sequence of the rat gamma-crystallin gene region and comparison with an orthologous human region

The sequences of a 51-kb region containing the cluster of five rat gamma-crystallin-coding genes (CRYG) and of a 7-kb region surrounding the sixth rat CRYG gene were determined. Approximately 78% of the total sequence represents intergenic DNA. We also sequenced 22 kb of DNA from the human CRYG gene cluster. All CRYG genes are associated with CpG-rich regions. The sequence similarity between the human and rat gene regions drops sharply (to 65%) in intronic and 3'-flanking regions but decreases only gradually in the 5'-flanking region. Highly conserved regions (greater than 80%) are found as far upstream as 1.5 kb. Overall intergenic distances are conserved. The human region contains much more repetitive DNA (24% vs. 10%) but less simple-sequence (sps) DNA (0.7% vs. 4%) than the rat region. Almost all repeats and spsDNA elements are located in the intergenic region. The location of repetitive and spsDNA differs between the orthologous regions and these elements were probably inserted after the evolutionary separation of rat and man. The Alu repeats in man and the B3 repeats in the rat are close copies of their respective consensus sequences and bordered by virtually perfect repeats. In contrast, the B1 and B2 repeats in the rat have diverged considerably from the consensus sequence and the surrounding direct repeats are usually imperfect. Thus the dispersion of the B1 and B2 repeats in the rat probably preceded that of the B3 repeats. Within the rat genomic region the spacing of Z-DNA elements is surprisingly regular, they are located about 12 kb apart. A search for putative matrix-associated regions suggests that the rat CRYG gene cluster is organized into two chromosomal domains.     

Nucleotide sequence comparison of the rp49 gene region between Drosophila subobscura and D. melanogaster

A 1.6-kb fragment encompassing the rp49 gene, which codes for a ribosomal protein, has been cloned and sequenced in Drosophila subobscura. The rp49 coding region has accumulated 46 nucleotide differences out of 402 bp since D. subobscura diverged from D. melanogaster. Forty-three percent of the effectively silent sites have changed since both species diverged. Both silent and replacement differences are distributed at random between the two exons of the gene. The frequency of silent differences in exons does not differ from that observed in the 5' leader sequence and in the intron. The frequency of silent differences in exon and intron sites is much greater than the number of amino acid replacement differences. This observation indicates strong purifying selection against amino acid replacements.      

Nucleotide polymorphism in the RpII215 gene region of the insular species Drosophila guanche: reduced efficacy of weak selection on synonymous variation

An approximately 6.9-kb region encompassing the RpII215 gene was sequenced for 24 individuals of the island endemic species Drosophila guanche. The comparative analysis of synonymous polymorphism and divergence in D. guanche and D. subobscura, two species with pronounced differences in population size, allows contrasting the nearly neutral character of synonymous mutations. In D. guanche, unlike in D. subobscura, (1) the ratio of preferred to unpreferred synonymous changes was similar for polymorphic and fixed changes, (2) the numbers of preferred and unpreferred changes, both polymorphic and fixed, could be explained by the mutational process, and (3) the estimated scaled selection coefficient for unpreferred mutations did not differ significantly from zero. Additionally, the comparative analysis revealed that both the ratio of preferred to unpreferred synonymous changes and the frequency spectrum of unpreferred polymorphic mutations differed significantly between species. All these results indicate that a large fraction of synonymous mutations in the RpII215 gene behave as effectively neutral in D. guanche, whereas they are weakly selected in D. subobscura. The reduced efficacy of selection in the insular species constitutes strong evidence of the nearly neutral character of synonymous mutations and, therefore, of the role of weak selection in maintaining codon bias.