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.     

Synonymous Substitutions in the Xdh Gene of Drosophila: Heterogeneous Distribution along the Coding Region

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Evolutionary conservation of 5' upstream sequence of nine genes between human and great apes

Nucleotide sequences of nine 5' upstream gene regions for human, chimpanzee, gorilla, and orangutan were determined. We estimated nucleotide differences (d) for each region between human and great apes. The overall d was 0.027 (ranged from 0.004 to 0.052). Rates of nucleotide substitution were estimated by using d and divergence times of human, chimpanzee, gorilla, and orangutan. The overall rate of nucleotide substitution between human and other hominoids was estimated to be 0.52-0.85 x 10(-9). This rate in 5' upstream regions was lower than that of synonymous sites, suggesting that 5' upstream regions have evolved under some functional constraints. Because lower rates have been reported for coding sequences in primates compared to rodents, we also estimated the rate (1.17-1.76 x 10(-9)) of nucleotide substitutions for the corresponding 5' upstream regions in rodents (mouse/rat comparison). Thus the primate rate was lower than rodent rate also for the 5' upstream regions.     

Molecular drift of the bride of sevenless (boss) gene in Drosophila

DNA sequences were determined for three to five alleles of the bride-of- sevenless (boss) gene in each of four species of Drosophila. The product of boss is a transmembrane receptor for a ligand coded by the sevenless gene that triggers differentiation of the R7 photoreceptor cell in the compound eye. Population parameters affecting the rate and pattern of molecular evolution of boss were estimated from the multinomial configurations of nucleotide polymorphisms of synonymous codons. The time of divergence between D. melanogaster and D. simulans was estimated as approximately 1 Myr, that between D. teissieri and D. yakuba as approximately 0.75 Myr, and that between the two pairs of sibling species as approximately 2 Myr. (The boss genes themselves have estimated divergence times approximately 50% greater than the species divergence times.) The effective size of the species was estimated as approximately 5 x 10(6), and the average mutation rate was estimated as 1-2 x 10(-9)/nucleotide/generation. The ratio of amino acid polymorphisms within species to fixed differences between species suggests that approximately 25% of all possible single-step amino acid replacements in the boss gene product may be selectively neutral or nearly neutral. The data also imply that random genetic drift has been responsible for virtually all of the observed differences in the portion of the boss gene analyzed among the four species.      

Molecular Analysis of the Hot Spot Region Related to Length Mutations in Wheat Chloroplast Dnas. I. Nucleotide Divergence of Genes and Intergenic Spacer Regions Located in the Hot Spot Region

The nucleotide divergence of chloroplast DNAs around the hot spot region related to length mutation in Triticum (wheat) and Aegilops was analyzed. DNA sequences (ca. 4.5 kbp) of three chloroplast genome types of wheat complex were compared with one another and with the corresponding region of other grasses. The sequences region contained rbcL and psaI, two open reading frames, and a pseudogene, rpl23' (pseudogene for ribosomal protein L23) disrupted by AT-rich intergic spacer regions. The evolution of these genes in the closely related wheat complex is characterized by nonbiased nucleotide substitutions in terms of being synonymous/nonsynonymous, having A-T pressure transitions over transversions, and frequent changes at the third codon position, in contrast with the gene evolution among more distant plant groups where biased nucleotide substitutions have frequently occurred. The sequences of these genes had diverged almost in proportion to taxonomic distance. The sequence of the pseudogene rpl23' changed approximately two times faster than that of the coding region. Sequence comparison between the pseudogene and its protein-coding counterpart revealed different degrees of nucleotide homology in wheat, rice and maize, suggesting that the transposition timing of the pseudogene differed and/or that different rates of gene conversion operated on the pseudogene in the cpDNA of the three plant groups in Gramineae. The intergenic spacer regions diverged approximately ten times faster than the genes. The divergence of wheat from barley, and that from rice are estimated based on the nucleotide similarity to be 1.5, 10 and 40 million years, respectively.     

Isolation and characterization of the genomic region from Drosophila kuntzei containing the Adh and Adhr genes

The nucleotide sequences of the Adh and Adhr genes of Drosophila kuntzei were derived from combined overlapping sequences of clones isolated from a genomic library and from cloned PCR and inverse-PCR fragments. Only a proximal promoter was detected upstream of the Adh gene, indicating that D. kuntzei Adh is regulated by a one-promoter system. Further upstream of the Adh structural gene, an adult enhancer region (AAE) was found that contains most of the regulatory sequences described for AAEs of other Drosophila species. Analysis of the ADH protein showed an amino acid change from valine to threonine in the active site at position 189 which is also found in D. funebris but is otherwise unique among DROSOPHILA: This difference alone may be responsible for the very low ADH activity found in this species and may cause a difference in substrate usage pattern. Codon bias in Adh and Adhr was comparable and found to be very low compared with other species. Phylogenetic analysis showed that D. kuntzei is closest related to D. funebris and D. immigrans. The time of divergence between D. kuntzei and D. funebris was estimated to be 14.2-20.2 Myr and that between D. kuntzei-D. funebris and D. immigrans to be 30.8-44.0 Myr. An analysis of the genetic variation in the Adh gene and upstream sequences of four European strains showed that this gene was highly variable. Overall nucleotide diversity (pi) was 0.0139, which is two times higher than that in D. melanogaster.     

Molecular Population Genetics of Sex Determination Genes: The Transformer Gene of Drosophila Melanogaster

The transformer locus (tra) produces an RNA processing protein that alternatively splices the doublesex pre-mRNA in the sex determination hierarchy of Drosophila melanogaster. Comparisons of the tra coding region among Drosophila species have revealed an unusually high degree of divergence in synonymous and nonsynonymous sites. In this study, we tested the hypothesis that the tra gene will be polymorphic in synonymous and nonsynonymous sites within species by investigating nucleotide sequence variation in eleven tra alleles within D. melanogaster. Of the 1063 nucleotides examined, two synonymous sites were polymorphic and no amino acid variation was detected. Three statistical tests were used to detect departures from an equilibrium neutral model. Two tests failed to reject a neutral model of molecular evolution because of low statisitical power associated with low levels of genetic variation (Tajima/Fu and Li). The Hudson, Kreitman, and Aguade test rejected a neutral model when the tra region was compared to the 5'-flanking region of alcohol dehydrogenase (Adh). The lack of variability in the tra gene is consistent with a recent selective sweep of a beneficial allele in or near the tra locus.     

Evolution of hominoid mitochondrial DNA with special reference to the silent substitution rate over the genome

Summary Focusing on the synonymous substitution rate, we carried out detailed sequence analyses of hominoid mitochondrial (mt) DNAs of ca. 5-kb length. Owing to the outnumbered transitions and strong biases in the base compositions, synonymous substitutions in mtDNA reach rapidly a rather low saturation level. The extent of the compositional biases differs from gene to gene. Such changes in base compositions, even if small, can bring about considerable variation in observed synonymous differences and may result in the region-dependent estimate of the synonymous substitution rate. We demonstrate that such a region dependency is due to a failure to take proper account of heterogeneous compositional biases from gene to gene but that the actual synonymous substitution rate is rather uniform. The synonymous substitution rate thus estimated is 2.37 ± 0.11 × 10^–8 per site per year and comparable to the overall rate for the noncoding region. On the other hand, the rate of nonsynonymous substitutions differs considerably from gene to gene, as expected under the neutral theory of molecular evolution. The lowest rate is 0.8 × 10^–9 per site per year forCOI and the highest rate is 4.5 × 10^–9 forATPase 8, the degree of functional constraints (measured by the ratio of the nonsynonymous to the synonymous substitution rate) being 0.03 and 0.19, respectively. Transfer RNA (tRNA) genes also show variability in the base contents and thus in the nucleotide differences. The average rate for 11 tRNAs contained in the 5-kb region is 3.9 × 10^–9 per site per year. The nucleotide substitutions in the genome suggest that the transition rate is about 17 times faster than the transversion rate.     

The origin and evolution of Ebola and Marburg viruses

Molecular evolutionary analyses for Ebola and Marburg viruses were conducted with the aim of elucidating evolutionary features of these viruses. In particular, the rate of nonsynonymous substitutions for the glycoprotein gene of Ebola virus was estimated to be, on the average, 3.6 x 10(-5) per site per year. Marburg virus was also suggested to be evolving at a similar rate. Those rates were a hundred times slower than those of retroviruses and human influenza A virus, but were of the same order of magnitude as that of the hepatitis B virus. When these rates were applied to the degree of sequence divergence, the divergence time between Ebola and Marburg viruses was estimated to be more than several thousand years ago. Moreover, most of the nucleotide substitutions were transitions and synonymous for Marburg virus. This suggests that purifying selection has operated on Marburg virus during evolution.      

Molecular Population Genetics of an Electrophoretically Monomorphic Protein in the Alcohol Dehydrogenase Region of Drosophila Pseudoobscura

Nucleotide sequence data from the alcohol dehydrogenase (Adh) region of 18 isochromosomal strains of Drosophila pseudoobscura were used to determine whether the lack of amino acid polymorphism in ADH results from a low neutral mutation rate or a recent directional selection event. We estimated the neutral mutation parameter, 4Nmu, in synonymous sites for 17 subregions of Adh. The nucleotide diversity data were tested for departures from an equilibrium neutral model with two statistical tests. The Tajima test and the Hudson, Kreitman and Aguade test each failed to reject a neutral model. These results suggest that the ADH enzyme of D. pseudoobscura lacks amino acid polymorphisms because the neutral mutation rate of nonsynonymous sites is low. The neutral mutation parameter for synonymous sites is heterogeneous between domains of the Adh region. These data indicate that selective constrains on synonymous sites can vary between functional domains.     

Synonymous rates at the RpII215 gene of Drosophila: variation among species and across the coding region

The region encompassing the RpII215 gene that encodes the largest component of the RNA polymerase II complex (1889 amino acids) has been sequenced in Drosophila subobscura, D. madeirensis, D. guanche, and D. pseudoobscura. Nonsynonymous divergence estimates (Ka) indicate that this gene has a very low rate of amino acid replacements. Given its low Ka and constitutive expression, synonymous substitution rates are, however, unexpectedly high. Sequence comparisons have allowed the molecular clock hypothesis to be tested. D. guanche is an insular species and it is therefore expected to have a reduced effective size relative to D. subobscura. The significantly higher rate of synonymous substitutions detected in the D. guanche lineage could be explained if synonymous mutations behave as nearly neutral. Significant departure from the molecular clock hypothesis for synonymous and nonsynonymous substitutions was detected when comparing the D. subobscura, D. pseudoobscura, and D. melanogaster lineages. Codon bias and synonymous divergence between D. subobscura and D. melanogaster were negatively correlated across the RpII215 coding region, which indicates that selection coefficients for synonymous mutations vary across the gene. The C-terminal domain (CTD) of the RpII215 protein is structurally and functionally differentiated from the rest of the protein. Synonymous substitution rates were significantly different in both regions, which strongly indicates that synonymous mutations in the CTD and in the non-CTD regions are under detectably different selection coefficients.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

Structure and Evolution of the Adh Genes of Drosophila Mojavensis

The nucleotide sequence of the Adh region of Drosophila mojavensis has been completed and the region found to contain a pseudogene, Adh-2 and Adh-1 arranged in that order. Comparison of the sequence divergence of these genes to one another and to the Adh region of Drosophila mulleri and other species has allowed the development of a model for the evolution of the duplication of the Adh genes. There have been two major events. An initial duplication of an Adh gene whose dual promoter structure was similar to Drosophila melanogaster, resulted in a species with two Adh genes, one of which may have had only a proximal promoter. A second duplication of this gene generated an Adh region containing three genes. It is proposed that one of these is the ancestral gene having dual promoters, while the other two possess only proximal promoters. Subsequent events have resulted in both a change in the regulation of Adh-2 such that it is expressed as if it had a "distal" type promoter and the mutational inactivation of the most upstream gene resulting in the creation of a pseudogene. The sequence of the D. mojavensis Adh region has also revealed the presence of an element which is composed of juxtaposed inverted imperfectly repeated elements. There is a surprising and not fully explainable strong similarity of the nucleotide sequence of the 5' flanking region of the pseudogene in D. mojavensis and D. mulleri.