Drosophila melanogaster Cu, Zn superoxide dismutase gene sequence

Patterns of DNA sequence polymorphism at Sod vicinities in Drosophila melanogaster: unraveling the footprint of a recent selective sweep

We survey DNA sequence polymorphisms at the Sod locus and four neighboring regions of Drosophila melanogaster, spanning 55,513 base pairs (bp), in 15 strains from a natural population, plus one reference laboratory strain and one strain of Drosophila simulans. Our objective is to characterize a proposed selective sweep that occurred at a locus close to Sod in D. melanogaster and to characterize the strength of the selection event, its time, and the size of the hitchhiked region. Two regions, 1819 and 6kbr3r, show a pattern of polymorphism very similar to the one of Sod, implying that they have been affected by the same evolutionary process that impacted Sod. A third fragment, 2021 seems unaffected by the event. A fourth one, 4039, on the opposite flank of Sod in relation to 2021, is only partially affected. We estimate that the length of the chromosomal segment impacted by the selective sweep is 41-54 kb, the age of the selective sweep is 2,600-22,000 years, and the selective advantage is 0.020 < s < 0.103.

Reevaluation of phylogeny in the Drosophila obscura species group based on combined analysis of nucleotide sequences

The Drosophila obscura species group has served as an important model system in many evolutionary and population genetic studies. Despite the amount of study this group has received, some phylogenetic relationships remain unclear. While individual analysis of different nuclear, mitochondrial, allozyme, restriction fragment, and morphological data partitions are able to discern relationships among closely related species, they are unable to resolve relationships among the five obscura species subgroups. A combined analysis of several nucleotide data sets is able to provide resolution and support for some nodes not seen or well supported in analyses of individual loci. A phylogeny of the obscura species group based on combined analysis of nucleotide sequences from six mitochondrial and five nuclear loci is presented here. The results of several different combined analyses indicate that the Old World obscura and subobscura subgroups form a monophyletic clade, although they are unable to resolve the relationships among the major lineages within the obscura species group.

Characterization of a Cu/Zn superoxide dismutase-encoding gene region in Drosophila willistoni

A Cu/Zn superoxide dismutase-encoding gene (Sod) from Drosophila willistoni was cloned and sequenced. The gene shows a typical structure for a fruit-fly Sod gene, with a coding region of 462 bp in two exons separated by a 417-bp intron. Comparison of the Sod sequences from D. willistoni and D. melanogaster suggests that these species are only remotely related. Downstream from the Sod gene, there is an ORF on the opposite strand that putatively encodes the last exon of an unidentified gene. The polyadenylation signals of the two genes are separated by only 61 bp in D. willistoni, conforming to the common picture of compact dipteran genomes.

Phylogeny of Drosophila and related genera inferred from the nucleotide sequence of the Cu,Zn Sod gene

The phylogeny and taxonomy of the drosophilids have been the subject of extensive investigations. Recently, Grimaldi (1990) has challenged some common conceptions, and several sets of molecular data have provided information not always compatible with other taxonomic knowledge or consistent with each other. We present the coding nucleotide sequence of the Cu,Zn superoxide dismutase gene (Sod) for 15 species, which include the medfly Ceratitis capitata (family Tephritidae), the genera Chymomyza and Zaprionus, and representatives of the subgenera Dorsilopha, Drosophila, Hirtodrosophila, Scaptodrosophila, and Sophophora. Phylogenetic analysis of the Sod sequences indicates that Scaptodrosophila and Chymomyza branched off the main lineage before the major Drosophila radiations. The presence of a second intron in Chymomyza and Scaptodrosophila (as well as in the medfly) confirms the early divergence of these two taxa. This second intron became deleted from the main lineage before the major Drosophila radiations. According to the Sod sequences, Sophophora (including the melanogaster, obscura, saltans, and willistoni species groups) is older than the subgenus Drosophila; a deep branch splits the willistoni and saltans groups from the melanogaster and obscura groups. The genus Zaprionus and the subgenera Dorsilopha and Hirtodrosophila appear as branches of a prolific "bush" that also embraces the numerous species of the subgenus Drosophila. The Sod results corroborate in many, but not all, respects Throckmorton's (King, R.C. (ed) Handbook of Genetics. Plenum Press, New York, pp. 421-469, 1975) phylogeny; are inconsistent in some important ways with Grimaldi's (Bull. Am. Museum Nat. Hist. 197: 1-139, 1990) cladistic analysis; and also are inconsistent with some inferences based on mitochondrial DNA data. The Sod results manifest how, in addition to the information derived from nucleotide sequences, structural features (i.e., the deletion of an intron) can help resolve phylogenetic issues.

Evidence for positive selection in the superoxide dismutase (Sod) region of Drosophila melanogaster

DNA sequence variation in a 1410-bp region including the Cu,Zn Sod locus was examined in 41 homozygous lines of Drosophila melanogaster. Fourteen lines were from Barcelona, Spain, 25 were from California populations and the other two were from laboratory stocks. Two common electromorphs, SODS and SODF, are segregating in the populations. Our sample of 41 lines included 19 SodS and 22 SodF alleles (henceforward referred to as Slow and Fast alleles). All 19 Slow alleles were identical in sequence. Of the 22 Fast alleles sequenced, nine were identical in sequence and are referred to as the Fast A haplotypes. The Slow allele sequence differed from the Fast A haplotype at a single nucleotide site, the site that accounts for the amino acid difference between SODS and SODF. There were nine other haplotypes among the remaining 13 Fast alleles sequenced. The overall level of nucleotide diversity (pi) in this sample is not greatly different than that found at other loci in D. melanogaster. It is concluded that the Slow/Fast polymorphism is a recently arisen polymorphism, not an old balanced polymorphism. The large group of nearly identical haplotypes suggests that a recent mutation, at the Sod locus or tightly linked to it, has increased rapidly in frequency to around 50%, both in California and Spain. The application of a new statistical test demonstrates that the occurrence of such large numbers of haplotypes with so little variation among them is very unlikely under the usual equilibrium neutral model. We suggest that the high frequency of some haplotypes is due to natural selection at the Sod locus or at a tightly linked locus.

Compositional heterogeneity and patterns of molecular evolution in the Drosophila genome

The rates and patterns of molecular evolution in many eukaryotic organisms have been shown to be influenced by the compartmentalization of their genomes into fractions of distinct base composition and mutational properties. We have examined the Drosophila genome to explore relationships between the nucleotide content of large chromosomal segments and the base composition and rate of evolution of genes within those segments. Direct determination of the G + C contents of yeast artificial chromosome clones containing inserts of Drosophila melanogaster DNA ranging from 140-340 kb revealed significant heterogeneity in base composition. The G + C content of the large segments studied ranged from 36.9% G + C for a clone containing the hunchback locus in polytene region 85, to 50.9% G + C for a clone that includes the rosy region in polytene region 87. Unlike other organisms, however, there was no significant correlation between the base composition of large chromosomal regions and the base composition at fourfold degenerate nucleotide sites of genes encompassed within those regions. Despite the situation seen in mammals, there was also no significant association between base composition and rate of nucleotide substitution. These results suggest that nucleotide sequence evolution in Drosophila differs from that of many vertebrates and does not reflect distinct mutational biases, as a function of base composition, in different genomic regions. Significant negative correlations between codon-usage bias and rates of synonymous site divergence, however, provide strong support for an argument that selection among alternative codons may be a major contributor to variability in evolutionary rates within Drosophila genomes.

Sequence analysis of the indirect flight muscle actin-encoding gene of Drosophila simulans

A genomic clone of the Drosophila simulans indirect flight muscle actin-encoding gene (Act88F) has been isolated, and the sequence of a 3358-nucleotide segment corresponding to the Drosophila melanogaster Act88F transcription unit is presented. The ACt88F homologs in these two sibling species encode identical proteins and the general genomic organization of the Act88F locus is highly conserved, including the location of the transcription start point, and the size and position of intron/exon junctions. Substitutions within the 5' flanking region, however, are clearly nonuniform and the regions of lowest divergence coincide with regions that have been implicated in transcript accumulation and the regulation of tissue-specific expression. Silent substitutions within the coding regions have been compared to five other gene homologs in these sibling species. The rate of silent substitution at these loci varies more than threefold, suggesting selection at the codon level.

Rates of synonymous substitution and base composition of nuclear genes in Drosophila

We compared the rates of synonymous (silent) substitution among various genes in a number of species of Drosophila. First, we found that even for a particular gene, the rate of synonymous substitution varied considerably with Drosophila lineages. Second, we showed a large variation in synonymous substitution rates among nuclear genes in Drosophila. These rates of synonymous substitution were correlated negatively with C content and positively with A content at the third codon positions. Nucleotide sequences were also compared between pseudogenes and their functional homologs. The C content of the pseudogenes was lower than that of the functional genes and the A content of the former was higher than that of the latter. Because the synonymous substitution for functional genes and the nucleotide substitution for pseudogenes are exempted from any selective constraint at the protein level, these observations could be explained by a biased pattern of mutation in the Drosophila nuclear genome. Such a bias in the mutation pattern may affect the molecular clock (local clock) of each nuclear gene of each species. Finally, we obtained the average rates of synonymous substitution for three gene groups in Drosophila; 11.0 x 10(-9), 17.5 x 10(-9) and 27.1 x 10(-9)/site/year.

Structure and sequence of the Cu,Zn Sod gene in the Mediterranean fruit fly, Ceratitis capitata: intron insertion/deletion and evolution of the gene

We have cloned a 4-kb region encompassing the Cu,Zn superoxide dismutase (Sod) gene from a genomic library of the Mediterranean fruit fly, Ceratitis capitata, using a cDNA probe from Drosophila melanogaster. The coding sequence of 462 bases is equally as long as that in Drosophila species. The rate of amino acid replacement over the past 100 million years is approximately the same in the Diptera and in mammals, thus excluding the hypothesis (proposed to account for an apparent acceleration in rate of evolution of Sod over geological time) that the evolution of the SOD protein is much higher in the mammals than in other organisms. The coding region is interrupted by two introns in Ceratitis, whereas only one occurs in Drosophila. Phylogenetic comparisons indicate that the second intron was present in the common dipteran ancestor, but was lost shortly after the divergence of the Drosophila lineage from other Diptera. Analysis of the exon/intron structure of Sod in various animal phyla, plants, and fungi indicates that intron insertions as well as deletions have occurred in the evolution of the Sod gene.

A truncated P element is inserted in the transcribed region of the Cu,Zn SOD gene of an SOD "null" strain of Drosophila melanogaster

The decreased Cu,Zn SOD activity (less than 5%) in a "null" SODCA1 Drosophila melanogaster strain isolated in our laboratory is due to the insertion of a truncated P element into the transcribed region of the Cu,Zn SOD gene. Using a cDNA Cu,Zn SOD probe from a wild type D. melanogaster (F allele) we isolated an EcoRI Cu,Zn SOD clone from an EMBL3 genomic library of the SODCA1 strain, subcloned it, restriction-mapped and partially sequenced it. The 2.5 kb clone consists of a wild-type 1.84 kb EcoRI fragment containing the Cu,Zn SOD gene previously isolated in our laboratory, with an insertion of 0.68 kb derived (by an internal deletion) from an autonomous, 2.9 kb P element. The insertion starts 21 bp upstream from the coding sequence and causes an 8 bp target site duplication characteristic of P elements. A point mutation in the second exon results in a replacement of Asn by Lys at position 96, confirming that the mature protein encoded by the SOCCA1 is the same one encoded by the S allele, commonly found in natural populations. The diminished expression of SODCA1 allele is most possibly due to a reduction of the rate of transcription attributable to the insertion of the P element.