Organizational analysis of elav gene and functional analysis of ELAV protein of Drosophila melanogaster and Drosophila virilis.

Drosophila virilis genomic DNA corresponding to the D. melanogaster embryonic lethal abnormal visual system (elav) locus was cloned. DNA sequence analysis of a 3.8-kb genomic piece allowed identification of (i) an open reading frame (ORF) with striking homology to the previously identified D. melanogaster ORF and (ii) conserved sequence elements of possible regulatory relevance within and flanking the second intron. Conceptual translation of the D. virilis ORF predicts a 519-amino-acid-long ribonucleoprotein consensus sequence-type protein. Similar to D. melanogaster ELAV protein, it contains three tandem RNA-binding domains and an alanine/glutamine-rich amino-terminal region. The sequence throughout the RNA-binding domains, comprising the carboxy-terminal 346 amino acids, shows an extraordinary 100% identity at the amino acid level, indicating a strong structural constraint for this functional domain. The amino-terminal region is 36 amino acids longer in D. virilis, and the conservation is 66%. In in vivo functional tests, the D. virilis ORF was indistinguishable from the D. melanogaster ORF. Furthermore, a D. melanogaster ORF encoding an ELAV protein with a 40-amino-acid deletion within the alanine/glutamine-rich region was also able to supply elav function in vivo. Thus, the divergence of the amino-terminal region of the ELAV protein reflects lowered functional constraint rather than species-specific functional specification.     

Invasion of Drosophila virilis by the Penelope transposable element

The Penelope family of transposable elements (TEs) is broadly distributed in most species of the virilis species group of Drosophila. This element plays a pivotal role in hybrid dysgenesis in Drosophila virilis, in which at least four additional TE families are also activated. Here we present evidence that the Penelope family of elements has recently invaded D. virilis. This evidence includes: (1) a patchy geographical distribution, (2) genomic locations mainly restricted to euchromatic chromosome arms in various geographical strains, and (3) a high level of nucleotide similarity among members of the family. Two samples from a Tashkent (Middle Asia) population of D. virilis provide further support for the invasion hypothesis. The 1968 Tashkent strain is free of Penelope sequences, but all individuals collected from a 1997 population carry at least five Penelope copies. Furthermore, a second TE, Ulysses, has amplified and spread in this population. These results provide evidence for the Penelope invasion of a D. virilis natural population and the mobilization of unrelated resident transposons following the invasion.     

The intron boundaries and flanking rRNA coding sequences of Calliphora erythrocephala rDNA

We have sequenced the available cloned examples of the intron-coding sequence junctions for the rDNA of the higher Dipteran, Calliphora erythrocephala. The introns interrupt the rDNA at the same position as the type 1 intron family detected in Drosophila melanogaster and D. virilis (10,11). A duplication of 14 base pairs of the 28S rRNA coding sequence surrounds a short version of the major genomic length class of introns. This same duplication is associated with boundaries of the type 1 introns in D. virilis and D. melanogaster (10, 13,14). We have detected considerable homology between the 3' intron sequences of C. erythrocephala and D. virilis. The rRNA coding sequences flanking the introns are extremely homologous in C. erythrocephala, D. melanogaster and D. virilis, with only one small region of significant divergence. This corresponds to a variable stem region previously identified in eukaryotic 28S rRNA at a site analogous to the L1 ribosomal protein binding site of prokaryotic 23S rRNA (27).     

Drosophila virilis has long and highly polymorphic microsatellites

Comparative genomics is a powerful approach to inference of the dynamics of genome evolution. Most information about the evolution of microsatellites in the genus Drosophila has been obtained from Drosophila melanogaster. For comparison, we collected microsatellite data for the distantly related species Drosophila virilis. Screening about 0.5 Mb of nonredundant genomic sequence from GenBank, we identified 239 dinucleotide microsatellites. On average, D. virilis dinucleotides were significantly longer than D. melanogaster microsatellites (7.69 repeats vs. 6.75 repeats). Similarly, direct cloning of microsatellites resulted in a higher mean repeat number in D. virilis than in D. melanogaster (12.7 repeats vs. 12.2 repeats). Characterization of 11 microsatellite loci mapping to division 40-49 on the fourth chromosome of D. virilis indicated that D. virilis microsatellites are more variable than those of D. melanogaster.     

Coding region deletions associated with the major form of rDNA interruption in Drosophila

The nucleotide sequences at and around the termini of 5 kb type 1 interruptions in three separate clones of D. melanogaster rDNA repeats have been determined, and have been compared with the sequence of the corresponding region of an insertion-free rDNA repeat. All three interrupted rDNA repeats contain a small deletion of 28S rRNA coding material at the left coding/insertion sequence junction. A second deletion was found in one of the three clones, ad other aberrations were suggested by the results of restriction enzyme digestions of unfractionated rDNA. The termini of 5 kb type 1 rDNA insertions in D. melanogaster were also compared with the corresponding regions of 28S rDNA interruptions in D. virilis: the insertion site is identical in the two species, but the termini of the two species' interruptions show no homology. I sequenced a 1.1 kb region of the 5 kb type 1 D. melanogaster rDNA interruption that covers the sequences of the 1 kb and 0.5 kb insertions. There is 98% homology between the rightmost 1 kb of the 5 kb interruption and the sequences of the shorter insertions. Data suggest that Drosophila rDNA interruptions arose as a transposable element, and that divergence had included length alterations generated by unequal crossing over.     

Evolution of the Autosomal Chorion Locus in Drosophila. I. General Organization of the Locus and Sequence Comparisons of Genes S15 and S19 in Evolutionarily Distant Species

We have isolated clones corresponding to the autosomal chorion locus of Drosophila melanogaster, from two distantly (D. virilis and D. grimshawi) and one closely (D. subobscura) related species. In all the species the locus is unique within the genome and encompasses the same four chorion genes and an adjacent nonchorion gene, in the same order. In all species the locus specifically amplifies in the ovary, as in D. melanogaster. We present the nucleotide sequences of DNA segments that total 8.3 kb in length and include gene s15-1 from D. subobscura, D. virilis, and D. grimshawi as well as gene s19-1 from D. subobscura and D. grimshawi. They show clearly nonuniform rates of divergence, both within and outside the limits of the genes. Highlighted by a background of extensive sequence divergence elsewhere in the extragenic region, highly conserved elements are observed in the 5' flanking DNA and might represent regulatory elements.     

Promoter analysis of the Drosophila melanogaster gene encoding the pterin 4alpha-carbinolamine dehydratase.

Pterin-4alpha-carbinolamine dehydratase (PCD) is a key enzyme in the regeneration pathway of tetrahydrobiopterin. Previously, we isolated and reported the Drosophila melanogaster gene encoding PCD. In the present study, we isolated and characterized the Drosophila virilis gene encoding PCD. The Drosophila virilis PCD gene has two introns and an open reading frame to encode a protein of 101 amino acids. The amino acid sequence of Drosophila virilis PCD shows a 83% homology to that of the Drosophila melanogaster PCD protein. From the alignment of the nucleotide sequence in the 5'-flanking region of the Drosophila melanogaster and Drosophila virilis PCD genes, we found four conserved sequences. Using a transient transfection assay, we showed that one of the conserved sequences (-127 to approximately -115) is critical for expression, also the minimal promoter region between -127 and +51 is necessary for the efficient expression of Drosophila melanogaster PCD.     

Colonization of heterochromatic genes by transposable elements in Drosophila

As a further step toward understanding transposable element-host genome interactions, we investigated the molecular anatomy of introns from five heterochromatic and 22 euchromatic protein-coding genes of Drosophila melanogaster. A total of 79 kb of intronic sequences from heterochromatic genes and 355 kb of intronic sequences from euchromatic genes have been used in Blast searches against Drosophila transposable elements (TEs). The results show that TE-homologous sequences belonging to 19 different families represent about 50% of intronic DNA from heterochromatic genes. In contrast, only 0.1% of the euchromatic intron DNA exhibits homology to known TEs. Intraspecific and interspecific size polymorphisms of introns were found, which are likely to be associated with changes in TE-related sequences. Together, the enrichment in TEs and the apparent dynamic state of heterochromatic introns suggest that TEs contribute significantly to the evolution of genes located in heterochromatin.     

The chorion genes of the medfly, Ceratitis capitata, I: Structural and regulatory conservation of the s36 gene relative to two Drosophila species.

We have used low stringency screening with the Drosophila melanogaster s36 chorion gene to recover its homologue from genomic and cDNA libraries of the medfly, Ceratitis capitata. The same gene has also been recovered from a genomic library of D. virilis. The medfly s36 gene shows similar developmental specificity as in Drosophila (early choriogenesis). It is also specifically amplified in ovarian follicles; this is the first report of chorion gene amplification outside the genus Drosophila. Alignments of s36 sequences from three species show that, in addition to its regulatory conservation, the s36 gene is extensively conserved in sequence, in a region corresponding to a central protein domain, and in short regions of 5' flanking DNA that might correspond to cis-regulatory elements.     

Conserved Alternative Splicing Patterns and Splicing Signals in the Drosophila Sodium Channel Gene Para

We cloned genomic DNA corresponding to the Drosophila virilis homologue of para, a gene encoding a sodium channel α-subunit, and obtained many partial cDNA clones from embryos and adults. Para protein has been well conserved, and the optional elements at six different sites of alternative splicing in D. melanogaster are present in D. virilis, in addition to one new optional exon. Among 31 different splice-types observed in D. virilis, the stage-specific pattern of alternative splicing seen in D. melanogaster is also conserved. Comparison of genomic DNA sequence revealed three aspects that vary between alternatively and constitutively used exon sequences. Sixteen short blocks (10-75 bp), the only recognizably conserved intron sequence, were disproportionately associated with alternatively used splice sites. Silent site substitutions were found much less frequently in alternative than constitutive exon elements, and the degree of match to the Drosophila splice site consensus tended to be lower at less frequently selected alternative splice junctions. This study shows that the developmentally regulated variability of para products is highly conserved and therefore likely to be of functional significance and suggests that a variety of different sequence-dependent mechanisms may regulate this pattern of alternative splicing.     

Sequence and expression of the Drosophila phenylalanine hydroxylase mRNA

We report the cloning, nucleotide (nt) sequence and expression of the cDNA (pah) encoding phenylalanine hydroxylase (PAH) of Drosophila melanogaster. The strong hybridization signals observed in genomic blots when D. melanogaster DNA was probed with 32P-labeled human pah cDNA, indicated the existence of a high degree of sequence similarity between the pah genes of both species. The length of the pah genomic fragment is about 30 to 40 kb. The cDNA contains 84 bp of the 5'-untranslated region, 1359 bp of the protein-coding region and 87 bp of the 3' region, with only one polyadenylation signal. The isolated cDNA is probably full-length, since the size of the D. melanogaster PAH mRNA is 1.5 kb. At the nt level, the similarity of the D. melanogaster cDNA with human and rat pah cDNAs is 57.9% and 58.1%, respectively. The highest similarities are restricted to the nt sequence coding for the presumed hydroxylation domain. There is no nt sequence similarity between the first three exons of the human pah gene and an equivalent fraction of the D. melanogaster pah gene. At the amino acid (aa) level, the similarity in the presumed hydroxylation domain is 88.5%, in which two motifs of the structure AGLLSSXXXL are found, where X represents any aa. It was interesting to notice the conservation of aa 408, 311 and 280, where mutations are associated with phenylketonuria in humans. We observed, moreover, that, as it occurs in humans and rats, the expression of the D. melanogaster pah gene is tissue-specific and temporally regulated.