Polymorphism and Divergence at a Drosophila Pseudogene Locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcpψ) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcpψ fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcpψ revealed lower levels of polymorphism ( & = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcpψ within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcpψ lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

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.     

The structure of the Adh locus of Drosophila mettleri: an intermediate in the evolution of the Adh locus in the repleta group of Drosophila.

Members of species of the mulleri and hydei subgroups of the repleta group of Drosophila have duplicate Adh genes. The Adh regions of D. mojavensis, D. mulleri, and D. hydei contain three genes--a pseudogene, Adh-2, and Adh-1--arranged 5' to 3'. To understand the evolution of the triplicate Adh structure, we have cloned and sequenced the Adh locus of D. mettleri. This region consists of a 5' pseudogene and a 3' functional Adh gene. On the basis of the structure and nucleotide sequence comparisons of Adh genes of D. mettleri and other species, we propose that an initial duplication of the ancestral Adh gene generated two Adh genes arranged in tandem. The more 5' Adh gene became a pseudogene, while the more 3' gene remained functional through all the developmental stages. A second duplication of this 3' gene resulted in Adh regions with three genes--a pseudogene, Adh-2, and Adh-1.     

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.     

Drosophila Alcohol Dehydrogenase Polymorphism and Carbon-13 Fluxes: Opportunities for Epistasis and Natural Selection

The influence of genetic variations in Drosophila alcohol dehydrogenase (ADH) on steady-state metabolic fluxes was studied by means of (13)C NMR spectroscopy. Four pathways were found to be operative during 8 hr of ethanol degradation in third instar larvae of Drosophila. Seven strains differed by 18-25% in the ratio between two major pathway fluxes, i.e., into glutamate-glutamine-proline vs. lactate-alanine-trehalose. In general, Adh genotypes with higher ADH activity exhibit a twofold difference in relative carbon flux from malate into lactate and alanine vs. α,α-trehalose compared to low ADH activity genotypes. Trehalose was degraded by the pentose-phosphate shunt. The pentose-phosphate shunt and malic enzyme could supply NADPH necessary for lipid synthesis from ethanol. Lactate and/or proline synthesis may maintain the NADH/NAD(+) balance during ethanol degradation. After 24 hr the flux into trehalose is increased, while the flux into lipids declines in Adh(F) larvae. In Adh(S) larvae the flux into lipids remains high. This co-ordinated nature of metabolism and the genotype-dependent differences in metabolic fluxes may form the basis for various epistatic interactions and ultimately for variations in organismal fitness.     

Correct developmental expression of a cloned alcohol dehydrogenase gene transduced into the Drosophila germ line

We have used P-element-mediated transformation to introduce a cloned Drosophila alcohol dehydrogenase (Adh) gene into the germ line of ADH null flies. Six independent transformants expressing ADH were identified by their acquired resistance to ethanol. Each transformant carries a single copy of the cloned Adh gene in a different chromosomal location. Four of the six transformant lines exhibit normal Adh expression by the following criteria: quantitative levels of ADH enzyme activity in larvae and adults; qualitative tissue specificity; the size of stable Adh mRNA; and the characteristic developmental switch in utilization of two different Adh promoters. The remaining two transformants express ADH enzyme activity with the correct tissue specificity, but at a lower level than wild type. These results demonstrate that an 11.8 kb chromosomal fragment containing the Adh gene includes the cis-acting sequences necessary for its correct developmental expression, and that a variety of chromosomal sites permit proper Adh gene function.     

Drosophila fat body protein P6 and alcohol dehydrogenase are derived from a common ancestral protein

Summary Drosophila melanogaster alcohol dehydrogenase is an example of convergent evolution: it is not related to the ADHs of other organisms, but to short-chain dehydrogenases, which until now have been found only in bacteria and in mammalian steroid hormone metabolism. We present evidence that theDrosophila ADH is phylogenetically more closely related to P6, another highly expressed protein from the fat body ofDrosophila, than it is to the short-chain dehydrogenases. The polypeptide sequence of P6 was inferred from DNA sequence analysis. Both ADH and P6 polypeptides have retained a high structural similarity with respect to the Chou-Fasman prediction of secondary structure and hydropathy. P6 is also homologous to the 25-kd protein from the fat body ofSarcophaga peregrina, whose sequence we have reexamined. The evolution of the P6-ADH family of proteins is characterized by a dramatic increase in the methionine content of P6. Methionine accounts for 20% of P6 amino acids. This is in contrast with the absence of this amino acid in mature ADH. There is evidence that P6 and the 25-kd protein have undergone a parallel and independent enrichment in methionine. When corrected for this, the rate of amino acid replacement shows that the P6-25-kd lineage diverged from insect ADH shortly before the divergence of the ADH gene (Adh) from its 3-duplication (Adhdup).     

The Role of Gene Duplication and Unconstrained Selective Pressures in the Melanopsin Gene Family Evolution and Vertebrate Circadian Rhythm Regulation

Melanopsin is a photosensitive cell protein involved in regulating circadian rhythms and other non-visual responses to light. The melanopsin gene family is represented by two paralogs, OPN4x and OPN4m, which originated through gene duplication early in the emergence of vertebrates. Here we studied the melanopsin gene family using an integrated gene/protein evolutionary approach, which revealed that the rhabdomeric urbilaterian ancestor had the same amino acid patterns (DRY motif and the Y and E conterions) as extant vertebrate species, suggesting that the mechanism for light detection and regulation is similar to rhabdomeric rhodopsins. Both OPN4m and OPN4x paralogs are found in vertebrate genomic paralogons, suggesting that they diverged following this duplication event about 600 million years ago, when the complex eye emerged in the vertebrate ancestor. Melanopsins generally evolved under negative selection (ω = 0.171) with some minor episodes of positive selection (proportion of sites = 25%) and functional divergence (θ_I = 0.349 and θ_II = 0.126). The OPN4m and OPN4x melanopsin paralogs show evidence of spectral divergence at sites likely involved in melanopsin light absorbance (200F, 273S and 276A). Also, following the teleost lineage-specific whole genome duplication (3R) that prompted the teleost fish radiation, type I divergence (θ_I = 0.181) and positive selection (affecting 11% of sites) contributed to amino acid variability that we related with the photo-activation stability of melanopsin. The melanopsin intracellular regions had unexpectedly high variability in their coupling specificity of G-proteins and we propose that Gq/11 and Gi/o are the two G-proteins most-likely to mediate the melanopsin phototransduction pathway. The selection signatures were mainly observed on retinal-related sites and the third and second intracellular loops, demonstrating the physiological plasticity of the melanopsin protein group. Our results provide new insights on the phototransduction process and additional tools for disentangling and understanding the links between melanopsin gene evolution and the specializations observed in vertebrates, especially in teleost fish.     

Genomic Organization and Evolution of Alternative Exons in a Drosophila Calcium Channel Gene

The genomic organization of a gene coding for an α1 subunit of a voltage-gated calcium channel of Drosophila melanogaster (Dmca1A) was determined. Thirty-four exons, distributed over 45 kb of genomic sequence, have been identified and mapped, including exons in three regions involved in alternative splicing and new sites potentially involved in RNA editing. The comparison of the intron/exon boundaries of this channel with a mammalian counterpart shows that the genomic structure of these two genes has remained fairly similar during evolution, with more than half of the Drosophila intron positions being perfectly conserved compared to the human channel. Phylogenetic analysis of the mutually exclusive alternative exons revealed that they have diverged considerably. It is suggested that this divergence, rather than reflecting evolutionary age, is the likely result of accelerated rates of evolution following duplication.     

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.     

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.     

The evolution of an alpha-esterase pseudogene inactivated in the Drosophila melanogaster lineage

Previous analyses of the alpha-esterase cluster of Drosophila melanogaster revealed 10 active genes and the DmalphaE4a-Psi pseudogene. Here, we reconstruct the evolution of the pseudogene from the sequences of 12 alleles from widely scattered D. melanogaster populations and single alleles from Drosophila simulans and Drosophila yakuba. All of the DmalphaE4a-Psi alleles contain numerous inactivating mutations, suggesting that pseudogene alleles are fixed in natural populations. Several lines of evidence also suggest that DmalphaE4a is now evolving without selective constraint in the D. melanogaster lineage. There are three polymorphic indels which result in frameshifts; a key nucleotide of the intron splice acceptor is polymorphic; the neutral mutation parameter is the same for replacement and silent sites; one of the nonsilent polymorphisms results in a stop codon; only 1 of the 13 replacement polymorphisms is biochemically conservative; residues that are conserved among active esterases have different states in DmalphaE4a-Psi; and there are about half as many transitional polymorphisms as transversional ones. In contrast, the D. simulans and D. yakuba orthologs DsalphaE4a and DyalphaE4a do not have the inactivating mutations of DmalphaE4a-Psi and appear to be evolving under the purifying selection typical of protein- encoding genes. For instance, there have been more substitutions in the introns than in the exons, and more in silent sites than in replacement sites. Furthermore, most of the amino acid substitutions that have occurred between DyalphaE4a and DsalphaE4a are located in sites that typically vary among active alpha-esterases rather than those that are usually conserved. We argue that the original alphaE4a gene had a function which it has lost since the divergence of the D. melanogaster and D. simulans lineages.