Evolution of nucleotide substitutions and gene regulation in the amylase multigenes in Drosophila kikkawai and its sibling species

In order to determine evolutionary changes in gene regulation and the nucleotide substitution pattern in a multigene family, the amylase multigenes were characterized in Drosophila kikkawai and its sibling species. The nucleotide substitution pattern was investigated. Drosophila kikkawai has four amylase genes. The Amy1 and Amy2 genes are a head-to-head duplication in the middle of the B arm of the second chromosome, while the Amy3 and Amy4 genes are a tail-to-tail duplication near the centromere of the same chromosome. In the sibling species of D. kikkawai (Drosophila bocki, Drosophila leontia, and Drosophila lini), sequencing of the Amy1, Amy2, Amy3, and Amy4 genes revealed that the Amy1 and Amy2 gene group diverged from Amy3 and Amy4 after duplication. In the Amy1 and Amy2 genes, the divergent evolution occurred in the flanking regions; in contrast, the coding regions have evolved in concerted fashion. The electrophoretic pattern of AMY isozymes was also examined. In D. kikkawai and its siblings, two or three electrophoretically different isozymes are encoded by the Amy1 and Amy2 genes (S isozyme) and by the Amy3 and Amy4 genes (F (M) isozymes). The S and F (M) isozymes show different patterns of band intensity when larvae and flies were fed in different media. Amy1 and Amy2, which encode the S isozyme, are more strikingly regulated than Amy3 and Amy4, which encode the F (M) isozyme. The GC content and codon usage bias were higher for the Amy1 and Amy2 genes than for the Amy3 and Amy4 genes. Although the ratio of synonymous and replacement substitutions within the Amy1 and Amy2 gene group was not significantly different from that within the Amy3 and Amy4 gene group, the synonymous substitution rate in the lineage of Amy1 and Amy2 was lower than that of Amy3 and Amy4. In conclusion, after the first duplication but before speciation of four species, the synonymous substitution rate between the two lineages and the electrophoretic pattern of the isozymes encoded by them changed, although we do not know whether there was any evolutionary relationship between the two.     

Molecular evolution of duplicated amylase gene regions in Drosophila melanogaster: evidence of positive selection in the coding regions and selective constraints in the cis-regulatory regions

In this study, we randomly sampled Drosophila melanogaster from Japanese and Kenyan natural populations. We sequenced duplicated (proximal and distal) Amy gene regions to test whether the patterns of polymorphism were consistent with neutral molecular evolution. F(st) between the two geographically distant populations, estimated from Amy gene regions, was 0.084, smaller than reported values for other loci, comparing African and Asian populations. Furthermore, little genetic differentiation was found at a microsatellite locus (DROYANETSB) in these samples (G'st = -0.018). The results of several tests (Tajima's, Fu and Li's, and Wall's tests) were not significantly different from neutrality. However, a significantly higher level of fixed replacement substitutions was detected by a modified McDonald and Kreitman test for both populations. This indicates that positive selection occurred during or immediately after the speciation of D. melanogaster. Sliding-window analysis showed that the proximal region 1, a part of the proximal 5' flanking region, was conserved between D. melanogaster and its sibling species, D. simulans. An HKA test was significant when the proximal region 1 was compared with the 5' flanking region of Alcohol dehydrogenase (Adh), indicating a severe selective constraint on the Amy proximal region 1. These results suggest that natural selection has played an important role in the molecular evolution of Amy gene regions in D. melanogaster.     

Molecular Evolution of the Duplicated Amy Locus in the Drosophila Melanogaster Species Subgroup: Concerted Evolution Only in the Coding Region and an Excess of Nonsynonymous Substitutions in Speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a ``specialist' species for host plants and probably also occurs in the process of speciation in general.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.      

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.      

Evidence for gene conversion in the amylase multigene family of Drosophila pseudoobscura

The alpha-amylase (Amy) multigene family in Drosophila pseudoobscura is located on the third chromosome, which is polymorphic for more than 40 inverted gene arrangements. The number of copies in this family ranges from one to three, depending on the arrangement in question. A previous study of the three Amy genes from the Standard (ST) arrangement suggested either that duplicated copies (Amy2 and Amy3) are functionally constrained or that they are undergoing gene conversion with Amy1. In order to elucidate further the pattern of molecular evolution in this family, we cloned and sequenced four additional Amy genes, two from the Santa Cruz (SC) and two from the Chiricahua (CH) gene arrangement. Of the two alternatives, only the hypothesis of gene conversion is supported by the sequence analysis. The homogenization effect of gene conversion has been strongest in SC, whose copies differ by only two nucleotides, less noticeable in ST, and negligible in the CH. Furthermore, the action of gene conversion is apparently localized, occurring only in the coding region. Interestingly, these results concur with the findings of other workers for the duplicated Amy genes in the Drosophila melanogaster group. Thus, the occurrence of gene conversion in the Amy multigene family seems to be a common feature in the Drosophila species studied so far.      

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.      

Nucleotide variation of seven genes in Drosophila kikkawai

We examined levels and patterns of nucleotide variation in 21 strains of Drosophila kikkawai from Miyako island, Japan for the partial regions of the following seven nuclear genes: Adh, Ddc, esc, ksr, Pgi, su(f), and Tpi. The nucleotide variation at total sites (pi(t)) ranged from 0.0013 in the ksr, to 0.0173 in the Adh. The nucleotide divergence at total sites (K(t)) between D. kikkawai and D. lini ranged from 0.0286 in the Tpi to 0.0687 in the su(f). The levels of nucleotide polymorphism and divergence were heterogeneous among the investigated gene regions. The HKA test, which tests imbalance between the intra and interspecific nucleotide variation, showed that the intraspecific nucleotide variation in the Pgi region was much lower than the interspecific variation, while intraspecific variation in the Tpi region was only slightly lower than interspecific variation. The MK test showed an excess of low frequency replacement polymorphic changes in the Adh region, suggesting that most replacement mutations are deleterious. Fay and Wu's test detected an excess of newly arisen variants in the Ddc region. In total, four of the seven gene regions showed significant deviation from the neutrality.     

Evolutionary Relationships and Sequence Variation of α-Amylase Variants Encoded by Duplicated Genes in the Amy Locus of Drosophila Melanogaster

To infer the genealogical relationships of α-amylase electromorphs of Drosophila melanogaster, we determined the nucleotide sequences of a collection of electromorphs sampled throughout the world. On average there were 1.0 amino acid substitutions between identical electromorphs and 3.9 between different electromorphs, respectively. We found that the evolution of AMY(1) through AMY(6) electromorphs occurred by sequential accumulation of single amino acid substitutions each causing one charge difference. The nucleotide diversities at synonymous sites within Amy(1),Amy(2),Amy(3),Amy(4) and Amy(6) were 0.0321, 0.0000, 0.0355, 0.0059 and 0.0030, respectively. We also obtained evidence of genetic exchanges, such as intrachromosomal recombination, interchromosomal recombination or gene conversion, between the two duplicated Amy genes as well as among the alleles.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

Inference of positive and negative selection on the 5' regulatory regions of Drosophila genes

Both positive selection and negative selection have been shown to drive the evolution of coding regions. It is of interest to know if the corresponding 5' regions of genes may be subjected to selection of comparable intensities. For such a comparison, we chose the Accessory gland protein (Acp) genes as our test case. About 700 bp and 600 bp for the 5' and coding regions, respectively, of eight previously unstudied genes were sequenced from 21 isogenic lines of D. melanogaster and one line from D. simulans. The ratio of divergence at the amino-acid replacement sites (A) over that at the synonymous sites (S) was twice the ratio for common polymorphism. Interestingly, the 5' region shows the same trend, with the 5'/S divergence ratio being 1.8 times higher than the 5'/S ratio for common polymorphism. There are several possible explanations for the 5'/S ratios, including demography, negative selection, and positive selection. Under normal conditions, positive selection is the most likely explanation. If that is true, about 45 to 50 percent of all fixed differences at both the replacement and 5' sites were adaptive, even though the substitution rate in the former is only half that of the latter (K(A)/K(S) approximately 0.3 vs. K(5')/K(S) approximately 0.6). As previous analyses have indicated, the inclusion of slightly deleterious polymorphism confounds the inference of positive selection. The analysis of published polymorphism data covering 97 verified 5' regions of Drosophila suggests more pronounced selective constraint on the 5' untranslated region and the core promoter (together corresponding to approximately 200 bp in this data set) when compared to the more distal portion of the 5' region of genes.     

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.     

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.