The molecular evolution of the alcohol dehydrogenase and alcohol dehydrogenase-related genes in the Drosophila melanogaster species subgroup

The DNA sequences of the Adh genes of three members of the Drosophila melanogaster species subgroup have been determined. This completes the Adh sequences of the eight species of this subgroup. Two species, D. yakuba and D. teissieri, possess processed Adh pseudogenes. In all of the species of the subgroup, a gene of unknown function, Adhr, is located about 300 bp 3' to Adh. Although this gene is experiencing a higher rate of synonymous substitution than Adh, it is more constrained at the amino acid level. Phylogenetic relationships between all eight members of the melanogaster subgroup have been analyzed using a variety of methods. All analyses suggested that the D. yakuba and D. teissieri pseudogenes have a single common ancestor, rather than evolving independently in each species, and that D. melanogaster is the sister species to D. simulans, D. sechellia, and D. mauritiana. The evolutionary relationships of the latter three species remain equivocal.      

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.     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.     

Nucleotide Variation and Divergence in the Histone Multigene Family in Drosophila Melanogaster

Nucleotide differences in the histone H3 gene family in Drosophila melanogaster were studied on three levels: (1) within a chromosome, (2) within a population and (3) between species (D. melanogaster and Drosophila simulans). The average difference within the H3 gene within a chromosome was 0.0040 per nucleotide site, about 52% of that within a population (0.0077). The proportion of divergent sites between the two species was 0.0575, which is about 8.5 times the difference within a species. The distribution of divergence between species was similar to that of variation within a species. Divergence and variation were noted to be greatest in the 3' noncoding region and least in the coding region. Values intermediate between these were found for the 5' noncoding region. Divergence and variation in silent sites exceeded those in the total coding region, thus indicating possible purifying selection for amino-acid-altering change. Phylogenetic relations among H3 genes and genetic differences on these three levels are evidence for the concerted evolution of the histone gene family. The molecular mechanism by which variation is produced and maintained is discussed.     

Light wavelength dependency of mating activity in the Drosophila melanogaster species subgroup

The action spectra of mating activity among the six species of the Drosophila melanogaster species subgroup were compared to understand how light wavelength affects mating activity. The species fell into three groups with respect to the action spectrum of mating activity. We chose one representative species from each of the three types for detailed study: D. melanogaster, D. sechellia and D. yakuba. The mating activities were investigated under three different light intensities of three monochromatic lights stimulus. Each species showed a unique spectral and intensity response. To know the evolutionary meaning of the light wavelength dependency of mating activity, we superimposed the type of action spectrum of mating activity in these six species on a cladogram. Mating inhibition under UV was conserved in evolution among these species. Furthermore we clarified that D. melanogaster showed low mating activity under UV because males courted less under UV.     

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.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

Microevolutionary divergence pattern of the segmentation gene hunchback in Drosophila

To study the microevolutionary processes shaping the evolution of the segmentation gene hunchback (hb) from Drosophila melanogaster, we cloned and sequenced the gene from 12 isofemale lines representing wild-type populations of D. melanogaster, as well as from the closely related species Drosophila sechellia, Drosophila orena, and Drosophila yakuba. We find a relatively low degree of sequence variation in D. melanogaster (theta = 0.0017), which is, however, consistent with its chromosomal location in a region of low recombination. Tests of neutrality do not reject a neutral-evolution model for the whole region. However, pairwise tests with different subregions indicate that there is a relative excess of polymorphic sites in the leader and the intron. Codon usage pattern analysis shows a particularly biased codon usage in the highly conserved regions, which is in line with the hypothesis that selection on translational accuracy is the driving force behind such a bias. A comparison of the expression pattern of hb in different sibling species of D. melanogaster reveals some regulatory changes in D. yakuba, which could be interpreted as changes in the timing of secondary expression domains.     

Expansion and evolution of the X-linked testis specific multigene families in the melanogaster species subgroup

The testis specific X-linked genes whose evolution is traced here in the melanogaster species subgroup are thought to undergo fast rate of diversification. The CK2?tes and NACβtes gene families encode the diverged regulatory β-subunits of protein kinase CK2 and the homologs of β-subunit of nascent peptide associated complex, respectively. We annotated the CK2βtes-like genes related to CK2?tes family in the D. simulans and D. sechellia genomes. The ancestor CK2βtes-like genes preserved in D. simulans and D. sechellia are considered to be intermediates in the emergence of the D. melanogaster specific Stellate genes related to the CK2?tes family. The CK2?tes-like genes are more similar to the unique autosomal CK2?tes gene than to Stellates, taking into account their peculiarities of polymorphism. The formation of a variant the CK2?tes gene Stellate in D. melanogaster as a result of illegitimate recombination between a NAC?tes promoter and a distinct polymorphic variant of CK2?tes-like ancestor copy was traced. We found a close nonrandom proximity between the dispersed defective copies of DINE-1 transposons, the members of Helitron family, and the CK2βtes and NACβtes genes, suggesting an involvement of DINE-1 elements in duplication and amplification of these genes.     

Evolutionary conservation of the chromosomal configuration and regulation of amylase genes among eight species of the Drosophila melanogaster species subgroup

Nuclear DNA was extracted from each of the eight species comprising the Drosophila melanogaster species subgroup. Southern hybridization of this DNA by using a molecular probe specific for the alpha-amylase coding region showed that the duplicated structure of the amylase locus, first found in D. melanogaster, is conserved among all species of the melanogaster subgroup. Evidence is also presented for the concerted evolution of the duplicated genes within each species. In addition, it is shown that the glucose repression of amylase gene expression, which has been extensively studied in D. melanogaster, is not confined to this species but occurs in all eight members of the species subgroup. Thus, both the duplicated gene structure and the glucose repression of Drosophila amylase gene activity are stable over extended periods of evolutionary time.      

tRNA derived insertion element in histone gene repeating unit of Drosophila melanogaster.

Analysis of 41 histone homologous clones from an isogenic gene library of Drosophila melanogaster showed that non-histone fragments interrupt the histone repetitive clusters at several sites. Long (L) and short (S) forms of the repeating units are distinguished by the insertion of 240 bp into the spacer between H1 and H3 of the L units; Each form appears to be clustered with its own kind. The complete DNA sequence of the histone 5.0 kb repeating unit was determined. Five histone genes (H1, H2A, H2B, H3, H4) were identified in a repeating unit and several sequence blocks common to the five histone genes were found in the 5'- and 3'-regions. The insertion sequence of 240 bp was found to be similar to the Alu family, an element derived from tRNA.     

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

Patterns of mutation and selection at synonymous sites in Drosophila

That natural selection affects molecular evolution at synonymous sites in protein-coding sequences is well established and is thought to predominantly reflect selection for translational efficiency/accuracy mediated through codon bias. However, a recently developed maximum likelihood framework, when applied to 18 coding sequences in 3 species of Drosophila, confirmed an earlier report that the Notch gene in Drosophila melanogaster was evolving under selection in favor of those codons defined as unpreferred in this species. This finding opened the possibility that synonymous sites may be subject to a variety of selective pressures beyond weak selection for increased frequencies of the codons currently defined as "preferred" in D. melanogaster. To further explore patterns of synonymous site evolution in Drosophila in a lineage-specific manner, we expanded the application of the maximum likelihood framework to 8,452 protein coding sequences with well-defined orthology in D. melanogaster, Drosophila sechellia, and Drosophila yakuba. Our analyses reveal intragenomic and interspecific variation in mutational patterns as well as in patterns and intensity of selection on synonymous sites. In D. melanogaster, our results provide little statistical evidence for recent selection on synonymous sites, and Notch remains an outlier. In contrast, in D. sechellia our findings provide evidence in support of selection predominantly in favor of preferred codons. However, there is a small subset of genes in this species that appear to be evolving under selection in favor of unpreferred codons, which indicates that selection on synonymous sites is not limited to the preferential fixation of mutations that enhance the speed or accuracy of translation in this species.     

Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements

Horizontal transmission has been well documented as a major mechanism for the dissemination of mariner-like elements (MLEs) among species. Less well understood are mechanisms that limit vertical transmission of MLEs resulting in the "spotty" or discontinuous distribution observed in closely related species. In this article we present evidence that the genome of the common ancestor of the melanogaster species subgroup of Drosophila contained an MLE related to the mellifera (honey bee) subfamily. Horizontal transmission, approximately 3-10 MYA, is strongly suggested by the observation that the sequence of the MLE in Drosophila erecta is 97% identical in nucleotide sequence with that of an MLE in the cat flea, Ctenocephalides felis. The D. erecta MLE has a spotty distribution among species in the melanogaster subgroup. The element has a high copy number in D. erecta and D. orena, a moderate copy number in D. teissieri and D. yakuba, and was apparently lost ("stochastic loss") in the lineage leading to D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In D. erecta, most copies are concentrated in the heterochromatin. Two copies from D. erecta, denoted De12 and De19, were cloned and sequenced, and they appear to be nonfunctional ("vertical inactivation"). It therefore appears that the predominant mode of MLE evolution is vertical inactivation and stochastic loss balanced against occasional reinvasion of lineages by horizontal transmission.