Nonfixed duplication containing the Adh gene and a truncated form of the Adhr gene in the Drosophila funebris species group: different modes of evolution of Adh relative to Adhr in Drosophila

Binding site number variation and high-affinity binding consensus of Myb-SANT-like transcription factor Adf-1 in Drosophilidae

There is a growing interest in the evolution of transcription factor binding sites and corresponding functional change of transcriptional regulation. In this context, we have examined the structural changes of the ADF-1 binding sites at the Adh promoters of Drosophila funebris and D. virilis. We detected an expanded footprinted region in D. funebris that contains various adjacent binding sites with different binding affinities. ADF-1 was described to direct sequence-specific DNA binding to sites consisting of the multiple trinucleotide repeat . The ADF-1 recognition sites with high binding affinity differ from this trinucleotide repeat consensus sequence and a new consensus sequence is proposed for the high-affinity ADF-1 binding sites. In vitro transcription experiments with the D. funebris and D. virilis ADF-1 binding regions revealed that stronger ADF-1 binding to the expanded D. funebris ADF-1 binding region only moderately lead to increased transcriptional activity of the Adh gene. The potential of this regional expansion is discussed in the context of different ADF-1 cellular concentrations and maintenance of the ADF-1 stimulus. Altogether, evolutionary change of ADF-1 binding regions involves both, rearrangements of complex binding site cluster and also nucleotide substitutions within sites that lead to different binding affinities.

Potential gene flow in natural populations of the Drosophila ananassae species cluster inferred from a nuclear mitochondrial pseudogene

A pseudogene with 94% similarity to mitochondrial cytochrome c oxidase subunit I (COI) was identified and localized to chromosome 4 of Drosophila ananassae. Because this chromosome is believed to have reduced recombination, its history can be traced using the pseudo-COI sequence. Pseudo-COI sequences were obtained from 27 iso-female lines of six taxa belonging to the D. ananassae species cluster in which reproductive isolation is incomplete. The phylogenetic network constructed from seven recognized haplotypes (#0-#6) indicated that different taxa inhabiting the same geographic area share the haplotypes: #1 from Papua New Guinean populations of D. ananassae and pallidosa-like-Wau; #2 from Papua New Guinean populations of D. ananassae, pallidosa-like, and papuensis-like; and #4 from South Pacific populations of D. ananassae and D. pallidosa. Taxon-K has a unique haplotype (#6), and 18 mutation steps separate it from the closest haplotype, #2. We discuss the possibility of chromosome 4 introgression beyond taxon boundaries.

Partial characterization and evolution of Adh-Adhr in Drosophila dunni

We sequenced 2123 bp of the Adh-Adhr genomic region of Drosophila dunni of the cardini group from two cloned DNA PCR fragments and from two cDNA clones of an Adh transcript. This comprises the Adh coding region and introns, 3' UTR, intergenic sequence, and most of Adhr, which is 260 bp downstream of Adh. Both genes have the typical Drosophila melanogaster Adh structure of three exons and two introns, except for changes in the putative 8 bp sequence involved in downregulation within the 3' UTR of Adh. Two amino acid substitutions could explain the low activity previously reported for this enzyme in D. dunni: Thr --> Lys at position 191 and Val --> Thr at position 189. D. dunni's Adh has the lowest codon bias reported so far for Drosophila species, and based on analysis of the nucleotide substitution rate, it is less conserved than Adhr.

Origin and evolution of a chimeric fusion gene in Drosophila subobscura, D. madeirensis and D. guanche

An understanding of the mutational and evolutionary mechanisms underlying the emergence of novel genes is critical to studies of phenotypic and genomic evolution. Here we describe a new example of a recently formed chimeric fusion gene that occurs in Drosophila guanche, D. madeirensis, and D. subobscura. This new gene, which we name Adh-Twain, resulted from an Adh mRNA that retrotransposed into the Gapdh-like gene, CG9010. Adh-Twain is transcribed; its 5' promoters and transcription patterns appear similar to those of CG9010. Population genetic and phylogenetic analyses suggest that the amino acid sequence of Adh-Twain evolved rapidly via directional selection shortly after it arose. Its more recent history, however, is characterized by slower evolution consistent with increasing functional constraints. We present a model for the origin of this new gene and discuss genetic and evolutionary factors affecting the evolution of new genes and functions.

Sequences upstream of the homologous cis-elements of the Adh adult enhancer of Drosophila are required for maximal levels of Adh gene transcription in adults of Scaptodrosophila lebanonensis

The evolution of cis-regulatory elements is of particular interest for our understanding of the evolution of gene regulation. The Adh gene of Drosophilidae shows interspecific differences in tissue-specific expression and transcript levels during development. In Scaptodrosophila lebanonensis adults, the level of distal transcripts is maximal between the fourth and eighth day after eclosion and is around five times higher than that in D. melanogaster Adh(S). To examine whether these quantitative differences are regulated by sequences lying upstream of the distal promoter, we performed in vitro deletion mutagenesis of the Adh gene of S. lebanonensis, followed by P-element-mediated germ-line transformation. All constructs included, as a cotransgene, a modified Adh gene of D. melanogaster (dAdh) in a fixed position and orientation that acted as a chromosomal position control. Using this approach, we have identified a fragment of 1.5 kb in the 5' region, 830 bp upstream of the distal start site, which is required to achieve maximal levels of distal transcript in S. lebanonensis. The presence of this fragment produces a 3.5-fold higher level of distal mRNA (as determined by real time quantitative PCR) compared with the D. melanogaster dAdh cotransgene. This region contains the degenerated end of a minisatellite sequence expanding farther upstream and does not correspond to the Adh adult enhancer (AAE) of D. melanogaster. Indeed, the cis-regulatory elements of the AAE have been identified by phylogenetic footprinting within the region 830 bp upstream of the distal start site of S. lebanonensis. Furthermore, the deletions Delta-830 and Delta-2358 yield the same pattern of tissue-specific expression, indicating that all tissue-specific elements are contained within the region 830 bp upstream of the distal start site.

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.

The probability of preservation of a newly arisen gene duplicate

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.

Evolutionary change in the structure of the regulatory region that drives tissue and temporally regulated expression of alcohol dehydrogenase gene in Drosophila funebris

The Adh locus of Drosophilidae is organized as a single gene transcribed from two spatially and temporally regulated promoters except in species of the repleta group, which have two single promoter genes. Here we show that in Drosophila funebris the Adh gene is transcribed from a single promoter, in both larva and adult, with qualitative and quantitative species specific-differences in tissue distribution. The gene is expressed in larval fat body but in other tissues such as gastric caeca, midgut and Malpighian tubules its expression is reduced compared to most Drosophilidae species, and in adults it is almost limited to the fat body. The comparative analysis of gene expression of two strains, which differ by a duplication, indicates that the cis elements necessary for this pattern of expression in larvae are included in the region of 1.55 kb upstream of the transcription initiation site. This new organization reveals the evolution of a different regulatory strategy to express the Adh gene in the subgenus Drosophila.

Mobile elements and the genesis of microsatellites in dipterans

Factors that influence the genesis and genomic distribution of microsatellite DNA are poorly understood. We have identified a novel class of Dipteran mobile elements, mini-me elements, which help elucidate both of these issues. These retroposons contain two internal proto-microsatellite regions that commonly expand into lengthy microsatellite repeats. These elements are highly abundant, accounting for approximately 1.2% of the Drosophila melanogaster genome, giving them the potential to be a prolific source of microsatellite DNA variation. They also give us the opportunity to observe the outcomes of multiple microsatellite genesis events (initiating from the same proto-microsatellite) at separate mini-me loci. Based on these observations, we determined that the genesis of microsatellites within mini-me elements occurs through two separate mutational processes: the expansion of preexisting tandem repeats and the conversion of sequence with high cryptic simplicity into tandemly repetitive DNA. These modes of microsatellite genesis can be generalized beyond the case of mini-me elements and help to explain the genesis of microsatellites in any sequence region that is not constrained by selection.

Molecular organization of the Drosophila melanogaster Adh chromosomal region in D. repleta and D. buzzatii, two distantly related species of the Drosophila subgenus

The molecular organization of a 1.944-Mb chromosomal region of Drosophila melanogaster around the Adh locus has been analyzed in two repleta group species: D. repleta and D. buzzatii. The extensive genetic and molecular information about this region in D. melanogaster makes it a prime choice for comparative studies of genomic organization among distantly related species. A set of 26 P1 phages from D. melanogaster were successfully hybridized using fluorescence in-situ hybridization (FISH) to the salivary gland chromosomes of both repleta group species. The results show that the Adh region is distributed in D. repleta and D. buzatii over six distant sites of chromosome 3, homologous to chromosomal arm 2L of D. melanogaster (Muller's element B). This observation implies a density of 2.57 fixed breakpoints per Mb in the Adh region and suggests a considerable reorganization of this chromosomal element via the fixation of paracentric inversions. Nevertheless, breakpoint density in the Adh region is three times lower than that estimated for D. repleta chromosome 2, homologous to D. melanogaster 3R (Muller's element E). Differences in the rate of evolution among chromosomal elements are seemingly persistent in the Drosophila genus over long phylogenetic distances.

Duplication, dicistronic transcription, and subsequent evolution of the Alcohol dehydrogenase and Alcohol dehydrogenase-related genes in Drosophila

It has recently been discovered that the Alcohol dehydrogenase and Alcohol dehydrogenase-related genes of Drosophila melanogaster and closely related species constitute a single transcription unit and that the Alcohol dehydrogenase-related gene is exclusively expressed from a dicistronic mRNA. Here, we show that in Drosophila lebanonensis, subgenus Scaptodrosophila, Adhr: is also transcribed as a dicistronic transcript with Adh Using degenerate primers designed on the sequence of the known Adhr proteins, we have been able to amplify and sequence a partial sequence of Adhr: in species representative of the whole subgenus Drosophila. This has allowed the study of the organization and expression of Adhr: in Drosophila buzzatii. We find that in D. buzzatii Adhr is transcribed as a monocistronic transcript. Adh and Adhr are believed to originate by duplication, and our data suggest that the cotranscription of these two genes was the primitive state, and that their independent transcription in the subgenus Drosophila is derived. We can rationalize the D. buzzatii condition as being correlated with the two genes evolving independent transcriptional control. However, why these two genes with clear divergence in the functions of their proteins should remain cotranscribed in groups as divergent as the subgenus Sophophora and the subgenus Scaptodrosophila remains a mystery.