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.     

Distribution of transposable elements in neotropical species of Drosophila

The phylogenetic distribution of transposable families, P, gypsy, hobo, I, and mariner has been analyzed in 33 species of 11 groups of neotropical Drosophila and a Drosophilidae species Zygotrica vittimaculosa, using squash blot and dot blot. Genomic DNA of almost all neotropical species tested hybridized with gypsy probe and some species showed a particularly strong hybridization signal, as D. gaucha, D. virilis, and species of flavopilosa group. The hobo element was restricted to melanogaster group and some strains of D. willistoni. Only D. simulans DNA showed hybridization to mariner probe in all species tested and D. simulans and D. melanogaster showed hybridization with I element probe. P element homologous sequence was present in D. melanogaster and all species and strains of the willistoni and saltans groups tested. The presence of at least one P-homologous sequence was detected in Drosophila mediopunctata. This one was the only P-bearing species of all six tested from the tripunctata group. Four different pairs of primers homologous to segments of the canonical sequence of D. melanogaster's P were used to amplify specific sequences from D. mediopunctata DNA, showing the occurrence of seemingly well-conserved P-homologous sequences. This revised version was published online in July 2006 with corrections to the Cover Date.     

Molecular structure of a gypsy element of Drosophila subobscura (gypsyDs) constituting a degenerate form of insect retroviruses.

We have determined the nucleotide sequence of a 7.5 kb full-size gypsy element from Drosophila subobscura strain H-271. Comparative analyses were carried out on the sequence and molecular structure of gypsy elements of D.subobscura (gypsyDs), D.melanogaster (gypsyDm) and D.virilis (gypsyDv). The three elements show a structure that maintains a common mechanism of expression. ORF1 and ORF2 show typical motifs of gag and pol genes respectively in the three gypsy elements and could encode functional proteins necessary for intracellular expansion. In the three ORF1 proteins an arginine-rich region was found which could constitute a RNA binding motif. The main differences among the gypsy elements are found in ORF3 (env-like gene); gypsyDm encodes functional env proteins, whereas gypsyDs and gypsyDv ORF3s lack some motifs essential for functionality of this protein. On the basis of these results, while gypsyDm is the first insect retrovirus described, gypsyDs and gypsyDv could constitute degenerate forms of these retroviruses. In this context, we have found some evidence that gypsyDm could have recently infected some D.subobscura strains. Comparative analyses of divergence and phylogenetic relationships of gypsy elements indicate that the gypsy elements belonging to species of different subgenera (gypsyDs and gypsyDv) are closer than gypsy elements of species belonging to the same subgenus (gypsyDs and gypsyDm). These data are congruent with horizontal transfer of gypsy elements among different Drosophila spp.     

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.     

Comparative and functional studies of Drosophila species invasion by the gypsy endogenous retrovirus.

Gypsy is an endogenous retrovirus of Drosophila melanogaster. Phylogenetic studies suggest that occasional horizontal transfer events of gypsy occur between Drosophila species. gypsy possesses infective properties associated with the products of the envelope gene that might be at the origin of these interspecies transfers. We report here the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster, suggesting that potentially infective gypsy copies able to spread between sexually isolated species can occur. The ability of gypsy to invade the genome of a new species is conditioned by its capacity to be expressed in the naive genome. The genetic basis for the regulation of gypsy activity in D. melanogaster is now well known, and it has been assigned to an X-linked gene called flamenco. We established an experimental simulation of the invasion of the D. melanogaster genome by gypsy elements derived from other Drosophila species, which demonstrates that these non- D. melanogaster gypsy elements escape the repression exerted by the D. melanogaster flamenco gene.     

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.     

Determination of some biochemical and structural features of alcohol dehydrogenases from Drosophila simulans and Drosophila virilis. Comparison of their properties with the Drosophila melanogaster Adhs enzyme.

The biochemical properties of the enzyme alcohol dehydrogenase of two different Drosophila species, Drosophila simulans and Drosophila virilis, were studied and compared with those of Drosophila melanogaster Adhs enzyme. All of them consist of two identical subunits of molecular weight 27800 and share significant similarities in function. The substrate specificities of these enzymes were characterized and Km(app.) and Vmax.(app.) values were calculated. All these alcohol dehydrogenases show greater affinity for secondary rather than for primary alcohols. The amino acid compositions of the three enzymes were determined, and there is a close similarity between the D. simulans and the D. melanogaster enzymes, but there are significant differences from the alcohol dehydrogenase of D. virilis. The N-terminal amino acid is blocked and the C-terminal amino acid is the same for all three alcohol dehydrogenases. The enzymes from the three species were carboxymethylated and digested with trypsin. The peptide 'maps' reveal, as expected, more homologies between the enzymes of D. simulans and D. melanogaster than with the enzyme of D. virilis.     

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.     

Regulative interactions between cells of wing discs from different dipteran species

The mechanism by which patterns are produced appears to be repeated in each segment of an animal, and it has been proposed that it may even have been conserved in evolution so that different species would have the same system of positional information. This idea has been tested by mixing cells of a defined fragment of the wing disc of Drosophila melanogaster with wing disc fragments of five other dipteran species to assay the ability of these disc fragments to stimulate intercalary regeneration of the D. melanogaster cells. The genetically marked (y; mwh) D. melanogaster fragment was mechanically mixed with wing discs or wing disc fragments of four drosophilids (D. melanogaster as a control, D. virilis, D. hydei, Zaprionus vittiger), of Musca domestica, and of Piophila casei. The mixed aggregates were cultured in vivo for 7 days, then metamorphosed in D. melanogaster larval hosts. The D. melanogaster fragments were only stimulated to regenerate when combined with complementary fragments from D. melanogaster or D. virilis wing discs. In the combination between D. melanogaster and D. hydei, the tissue formed integrated mosaic patterns, but no regeneration ensued. The one positive result (D. melanogaster mixed with D. virilis) shows that positional cues can be exchanged and correctly interpreted between cells of different species. The negative results do not prove that the mechanism for establishing patterns is different in the tested species, but may be due to incompatibilities that are not related to pattern formation.     

Satellite Ic: a possible link between the satellite DNAs of D. virilis and D. melanogaster.

In this study, we isolated and characterized a previously undetected cryptic satellite DNA comprising 0.1% of the total nuclear genome of D. virilis. This satellite is hidden from detection in neutral CsCl by satellite I and is therefore designated cryptic satellite I or Ic. Sequence analysis reveals that Ic is the repeating heptanucleotide [poly d(AATATAG): d(CTATATT)]. It is more closely related to the three simple sequence satellite DNAs of D. melanogaster, a distantly related species, than it is to any of the major D. virilis satellite DNA sequences. Ic may therefore be a link between the simple sequence satellites of D. virilis and D. melanogaster. As an extension of this theory, we have constructed a "family tree" linking the satellites of D. virilis and D. melanogaster by a series of "simple" operations. Only one intermediate required by this evolutionary scheme has not yet been identified.     

Retrotransposon mdg3 is transferred between somatic cells of unrelated species in coculture

Since retrovirus-like particles of gypsy (mdg4) are capable of interspecific transfer, other Drosophila melanogaster gypsy-related retrotransposons were tested for this property. As a donor and a recipient, D. melanogaster and D. virilis cultured cells were used. Recipient cell DNA was analyzed with probes directed to mdg1, mdg3, 17.6, 297, 412, or B104/roo. Transfer was demonstrated for mdg3, which lacks env. The possible mechanism of transfer is discussed.     

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.     

Distribution and conservation of mobile elements in the genus Drosophila

Essentially nothing is known of the origin, mode of transmission, and evolution of mobile elements within the genus Drosophila. To better understand the evolutionary history of these mobile elements, we examined the distribution and conservation of homologues to the P, I, gypsy, copia, and F elements in 34 Drosophila species from three subgenera. Probes specific for each element were prepared from D. melanogaster and hybridized to genomic DNA. Filters were washed under conditions of increasing stringency to estimate the similarity between D. melanogaster sequences and their homologues in other species. The I element homologues show the most limited distribution of all elements tested, being restricted to the melanogaster species group. The P elements are found in many members of the subgenus Sophophora but, with the notable exception of D. nasuta, are not found in the other two subgenera. Copia-, gypsy-, and F-element homologues are widespread in the genus, but their similarity to the D. melanogaster probe differs markedly between species. The distribution of copia and P elements and the conservation of the gypsy and P elements is inconsistent with a model that postulates a single ancient origin for each type of element followed by mating-dependent transmission. The data can be explained by horizontal transmission of mobile elements between reproductively isolated species.      

Retrotransposon gtwin: structural analysis and distribution in Drosophila melanogaster strains

A search for noncanonical variants of the gypsy retrotransposon (MDG4) in the genome of the Drosophila melanogaster strain G32 led to the cloning of four copies of the poorly studied 7411-bp gtwin element. Sequence analysis showed that gtwin belongs to a family of endogeneous retroviruses, which are widespread in the Drosophila genome and have recently been termed insect erantiviruses. The gtwin retrotransposon is evolutionarily closest to MDG4, as evident from a good alignment of their nucleotide sequences including ORF1 (the pol gene) and ORF3 (the env gene), as well as the amino acid sequences of their protein products. These regions showed more than 75% homology. The distribution of gtwin was studied in several strains of the genus Drosophila. While strain G32 contained more than 20 copies of the element, ten other D. melanogaster strains carried gtwin in two to six copies per genome. The gtwin element was not detected in D. hydei or D. virilis. Comparison of the cloned gtwin sequences with the gtwin sequence available from the D. melanogaster genome database showed that the two variants of the mobile element differ by the presence or absence of a stop codon in the central region of ORF3. Its absence from the gtwin copies cloned from the strain G32 may indicate an association between the functional state of ORF3 and amplification of the element.     

Evidence for recent horizontal transfer of gypsy-homologous LTR-retrotransposon gtwin into Drosophila erecta followed by its amplification with multiple aberrations

Long terminal repeat (LTR) retrotransposon gtwin was initially discovered in silico, and then it was isolated as gypsy-homologous sequence from Drosophila melanogaster strain, G32. The presence of ORF3 suggests, that gtwin, like gypsy, may be an endogenous retrovirus, which can leave the cell and infect another one. Therefore, in this study we decided to investigate the distribution of gtwin in different species of the melanogaster subgroup in order to find out whether gtwin can be transferred horizontally as well as vertically. Gtwin was found in all 9 species of this subgroup, hence it seems to have inhabited the host genomes for a long time. In addition, we have shown that in the Drosophila erecta genome two gtwin families are present. The first one has 93% of identity to D. melanogaster element and is likely to be a descendant of gtwin that existed in Drosophila before the divergence of the melanogaster subgroup species. The other one has >99% of identity to D. melanogaster gtwin. The most reasonable explanation is that this element has been recently horizontally transferred between D. melanogaster and D. erecta. The number and variety of gtwin copies from the "infectious" family suggest that after the horizontal transfer into D. erecta genome, gtwin underwent amplification and aberrations, leading to the rise of its diverse variants.     

Evidence for horizontal transfer of the LTR retrotransposon mdg3, which lacks an env gene

Horizontal (interspecific) transfer is regarded as a possible strategy for the propagation of transposable elements through evolutionary time. To date, however, conclusive evidence that transposable elements are capable of horizontal transfer from one species to another has been limited to class II or DNA-type elements. We tested the possibility of such transfer for several Drosophila melanogaster LTR retrotransposons of the gypsy group in an experiment in which D. melanogaster and D. virilis somatic cell lines were used as donor and recipient cells, respectively. This approach was chosen in light of the high levels of LTR retrotransposon amplification and expression observed in cultured D. melanogaster cells. In the course of the experiment, parallel analysis for mdg1, mdg3, 17.6, 297, 412 and B104/roo retrotransposons was performed to detect their presence in the genome of recipient cells. Only the mdg3 retrotransposon, which lacks an env gene, was found to be transmitted into recipient cells. This model, based on the use of cultured cells, is a promising system for further investigating the mechanisms of LTR retrotransposon transfer.