Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

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.     

The high turnover Drosophila multidrug resistance-associated protein shares the biochemical features of its human orthologues

DMRP, an ABC transporter encoded by the dMRP/CG6214 gene, is the Drosophila melanogaster orthologue of the “long” human multidrug resistance-associated proteins (MRP1/ABCC1, MRP2/ABCC2, MRP3/ABCC3, MRP6/ABCC6, and MRP7/ABCC10). In order to provide a detailed biochemical characterisation we expressed DMRP in Sf9 insect cell membranes. We demonstrated DMRP as a functional orthologue of its human counterparts capable of transporting several human MRP substrates like β-estradiol 17-β-d-glucuronide, leukotriene C₄, calcein, fluo3 and carboxydichlorofluorescein. Unexpectedly, we found DMRP to exhibit an extremely high turnover rate for the substrate transport as compared to its human orthologues. Furthermore, DMRP showed remarkably high basal ATPase activity (68-75 nmol P_i/mg membrane protein/min), which could be further stimulated by probenecid and the glutathione conjugate of N-ethylmaleimide. Surprisingly, this high level basal ATPase activity was inhibited by the transported substrates. We discussed this phenomenon in the light of a potential endogenous substrate (or activator) present in the Sf9 membrane.     

A 9.6 kb intervening sequence in D. virilis rDNA, and sequence homology in rDNA interruptions of diverse species of Drosophila and other diptera.

A large proportion of the 28S ribosomal RNA genes in Drosophila virilis are interrupted by a DNA sequence 9.6 kilobase pairs long. As regards both its presence and its position in the 28S gene (about two thirds of the way in), the D. virilis rDNA intervening sequence is similar to that found in D. melanogaster rDNA, but lengths differ markedly between the two species. Degrees of nucleotide sequence homology have been detected bewteen rDNA interruptions of the two species. This homology extends to putative rDNA intervening sequences in diverse higher diptera (other Drosophila species, the house fly and the flesh fly), but hybridization of cloned D. melanogaster and D. virilis rDNA interruption segments to DNA of several lower diptera has been negative. As is the case with melanogaster rDNA interruptions, segments of the virilis rDNA intervening sequence hybridize with non-rDNA components of the virilis genome, and interspecific homology may involve these non-rDNA sequences as well as rDNA interruptions. There is, however, evidence from buoyant density fractionation of DNA that the distributions of interruption-related sequences are distinct in D. melanogaster and D. virilis genomes. Moreover, thermal denaturation studies have indicated differing extents of homology between hybridizable sequences in D. virilis DNA and different segments of the D. melanogaster rDNA intervening sequence. We infer from our studies that rDNA intervening sequences are prevalent among higher diptera; that in the course of the evolution of these organisms, elements of the intervening sequences have been moderately to highly conserved; and that this conservation extends in at least two distantly related species of Drosophila to similar sequences found elsewhere in the genomes.     

Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

Evolution of gypsy endogenous retrovirus in the Drosophila obscura species group

The Ty3/gypsy family of retroelements is closely related to retroviruses, and some of their members have an open reading frame resembling the retroviral gene env. Sequences homologous to the gypsy element from Drosophila melanogaster are widely distributed among Drosophila species. In this work, we report a phylogenetic study based mainly on the analysis of the 5' region of the env gene from several species of the obscura group, and also from sequences already reported of D. melanogaster, Drosophila virilis, and Drosophila hydei. Our results indicate that the gypsy elements from species of the obscura group constitute a monophyletic group which has strongly diverged from the prototypic D. melanogaster gypsy element. Phylogenetic relationships between gypsy sequences from the obscura group are consistent with those of their hosts, indicating vertical transmission. However, D. hydei and D. virilis gypsy sequences are closely related to those of the affinis subgroup, which could be indicative of horizontal transmission.     

Phylogenetic and genomewide analyses suggest a functional relationship between kayak, the Drosophila fos homolog, and fig, a predicted protein phosphatase 2c nested within a kayak intron

A gene located within the intron of a larger gene is an uncommon arrangement in any species. Few of these nested gene arrangements have been explored from an evolutionary perspective. Here we report a phylogenetic analysis of kayak (kay) and fos intron gene (fig), a divergently transcribed gene located in a kay intron, utilizing 12 Drosophila species. The evolutionary relationship between these genes is of interest because kay is the homolog of the proto-oncogene c-fos whose function is modulated by serine/threonine phosphorylation and fig is a predicted PP2C phosphatase specific for serine/threonine residues. We found that, despite an extraordinary level of diversification in the intron-exon structure of kay (11 inversions and six independent exon losses), the nested arrangement of kay and fig is conserved in all species. A genomewide analysis of protein-coding nested gene pairs revealed that approximately 20% of nested pairs in D. melanogaster are also nested in D. pseudoobscura and D. virilis. A phylogenetic examination of fig revealed that there are three subfamilies of PP2C phosphatases in all 12 species of Drosophila. Overall, our phylogenetic and genomewide analyses suggest that the nested arrangement of kay and fig may be due to a functional relationship between them.     

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).     

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.     

Circular DNA Molecules in the Genus Drosophila

The satellite DNA's from the embryos of five species of Drosophila (D. melanogaster, D. simulans, D. nasuta, D. virilis and D. hydei) have been analyzed for the presence of closed circular duplex DNA molecules, as determined by CsCl-EBr gradients. Circular DNA molecules were found in every species but D. melanogaster. Analyses of cell fractions from adult Drosophila and organ fractions from Drosophila larvae show that fractions containing mitochondria are highly enriched in these molecules.     

Evolutionary conservation of regulatory strategies for the sex determination factor transformer-2.

Sex determination in Drosophila melanogaster is regulated by a cascade of splicing factors which direct the sex-specific expression of gene products needed for male and female differentiation. The splicing factor TRA-2 affects sex-specific splicing of multiple pre-mRNAs involved in sexual differentiation. The tra-2 gene itself expresses a complex set of mRNAs generated through alternative processing that collectively encode three distinct protein isoforms. The expression of these isoforms differs in the soma and germ line. In the male germ line the ratio of two isoforms present is governed by autoregulation of splicing. However, the functional significance of multiple TRA-2 isoforms has remained uncertain. Here we have examined whether the structure, function, and regulation of tra-2 are conserved in Drosophila virilis, a species diverged from D. melanogaster by over 60 million years. We find that the D. virilis homolog of tra-2 produces alternatively spliced RNAs encoding a set of protein isoforms analogous to those found in D. melanogaster. When introduced into the genome of D. melanogaster, this homolog can functionally replace the endogenous tra-2 gene for both normal female sexual differentiation and spermatogenesis. Examination of alternative mRNAs produced in D. virilis testes suggests that germ line-specific autoregulation of tra-2 function is accomplished by a strategy similar to that used in D. melanogaster. The similarity in structure and function of the tra-2 genes in these divergent Drosophila species supports the idea that sexual differentiation in D. melanogaster and D. virilis is accomplished under the control of similar regulatory pathways.     

Components Acting in Localization of Bicoid mRNA Are Conserved among Drosophila Species

Substantial insights into basic strategies for embryonic body patterning have been obtained from genetic analyses of Drosophila melanogaster. This knowledge has been used in evolutionary comparisons to ask if genes and functions are conserved. To begin to ask how highly conserved are the mechanisms of mRNA localization, a process crucial to Drosophila body patterning, we have focused on the localization of bcd mRNA to the anterior pole of the embryo. Here we consider two components involved in that process: the exuperantia (exu) gene, required for an early step in localization; and the cis-acting signal that directs bcd mRNA localization. First, we use the cloned D. melanogaster exu gene to identify the exu genes from Drosophila virilis and Drosophila pseudoobscura and to isolate them for comparisons at the structural and functional levels. Surprisingly, D. pseudoobscura has two closely related exu genes, while D. melanogaster and D. virilis have only one each. When expressed in D. melanogaster ovaries, the D. virilis exu gene and one of the D. pseudoobscura exu genes can substitute for the endogenous exu gene in supporting localization of bcd mRNA, demonstrating that function is conserved. Second, we reevaluate the ability of the D. pseudoobscura bcd mRNA localization signal to function in D. melanogaster. In contrast to a previous report, we find that function is retained. Thus, among these Drosophila species there is substantial conservation of components acting in mRNA localization, and presumably the mechanisms underlying this process.     

Transgenic tools for members of the genus Drosophila with sequenced genomes

The sequencing of the genomes of 12 Drosophila species has created an opportunity for much in the way of comparative molecular analyses amongst these species. To aid that endeavor, we have made several transformation vectors based on the piggyBac transposon with 3xP3-EGFP and -ECFP transgenic markers that should be useful for mutagenesis and establishing the GAL4/UAS system in these species. We have tested the ability of mini-white to be used as a marker for insertional mutagenesis, and have observed mini-white derived pigmentation of the testes sheath in a subset of lines from D. pseudoobscura and D. virilis. We have incorporated a source of piggyBac transposase into nine Drosophila species, and have demonstrated the functionality of these transposase lines for mobilization of marked inserts in vivo. Additionally, we tested the ability of a D. melanogaster nanos enhancer element to drive expression of GAL4 in D. melanogaster, D. simulans, D. erecta, D. yakuba, D. pseudoobscura, and D. virilis. The efficacy of the nos-Gal4 transgene was determined by measuring the response of UAS-EGFPtub in all six species. Our results show that D. melanogaster nos-Gal4 drives expression in other species, to varying degrees, in similar spatiotemporal domains in the ovaries, testes, and embryos as seen in D. melanogaster. However, expression levels are variable, demonstrating the possible need to use species-specific promoters in some cases. In summary, we hope to provide a set of guidelines and basic tools, based upon this work, for both insertional mutagenesis and GAL4/UAS system-based experiments in multiple species of Drosophila.     

A telomeric satellite in Drosophila virilis and its sibling species

Telomere elongation by telomerase is the most widespread mechanism among eukaryotes. However, alternative mechanisms such as homologous recombination between terminal satellite DNAs are probably used in lower dipteran insects and in some plants. Drosophila melanogaster uses the very unusual telomere elongation pathway of transposition of telomere-specific retrotransposable elements. The uniqueness of this telomere elongation mechanism raises the question of its origin. We, therefore, analyzed sequences located at telomeres of fairly distantly related Drosophila species, and in this paper we describe the characterization of complex satellite DNA sequences located at the telomeres of D. virilis and other species in the virilis group. We suggest an involvement of these DNA satellites in telomere elongation by homologous recombination similar to that found in lower dipterans. Our findings raise the possibility that telomere elongation by specific retrotransposons as found in D. melanogaster and its sibling species is a recent event in the evolution of dipteran insects.     

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.