The Evolutionary History of the Transposable Element Penelope in the Drosophila virilis Group of Species

We have used phylogenetic techniques to study the evolutionary history of the Penelope transposable element in the Drosophila virilis species group. Two divergent types of Penelope have been detected, one previously described, clade I, and a new one which we have termed clade III. The phylogeny of some copies of the Penelope clade I element was partially consistent with the species phylogeny of the D. montana subphylad, suggesting cospeciation and allowing the estimation of the evolutionary rate of Penelope. Divergence times of elements found in different species are younger than the age of the species, suggesting horizontal transfer events. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov]     

Distribution and evolution of mobile elements in the virilis species group of Drosophila

The distributions of Penelope and Ulysses, two transposable elements that can induce hybrid dysgenesis, were studied in several species groups of Drosophila. No significant hybridization to Penelope and Ulysses probes was detected by Southern blot analyses of species outside the virilis group. In contrast, both element families have had a long residence in all species of the virilis species group, as indicated by their strong presence in the heterochromatic chromocenter. Except for D. kanekoi, D. lummei, and some strains of D. virilis, species of the group carry full-sized, and at least potentially functional, copies of both element families. Consistent with the occurrence of recent transposition, Penelope and Ulysses elements are located at different chromosomal sites in different geographical strains of the same species. A total of 79 Penelope and 47 Ulysses euchromatic insertion sites were localized to chromosomal subsections in species of the virilis group. Highly significant deviations from independence of the distributions of Penelope and Ulysses and previously established inversion breakpoints were documented, suggesting that these transposable elements may have played an important role in genomic reorganization and evolution of the virilis species group, which is especially rich in karyotypic variation. Received: 13 April 1999; in revised form: 20 July 1999 / Accepted: 27 July 1999     

Identification of Two Penelope-Like Elements with Different Structures and Chromosome Localization in Kuruma Shrimp Genome

Penelope, originally found as a key element responsible for the hybrid dysgenesis in Drosophila virilis, has been widely conserved throughout eukaryotic genomes. In other organisms, they are often referred to as Penelope-like elements or PLEs. In this study, we found two types of PLEs, designated MjPLE01 and MjPLE02, from kuruma shrimp, Marsupenaeus japonicus. There was no observed nucleotide similarity between MjPLE01 and 02, and both elements differed from each other in terms of their structure; MjPLE02 has a distinctive endonuclease (EN) domain at the C-terminus while MjPLE01 do not. A phylogenetic tree that includes publicly available PLEs and TERTs showed that MjPLE01 and 02 were closely related to Coprina elements, which have been reported as an EN-deficient PLE, and to Penelope–Poseidon group, which possess an EN domain, respectively. Genomic Southern blot analysis using MjPLE01 as a probe showed several multiple bands that differ among individual shrimps. On the other hand, two major identical bands were observed when MjPLE02 was used. Colony hybridization showed co-localization of MjPLE01 and GGTTA repeats, suggesting that MjPLE01 might be prevalent in subtelomeric regions of kuruma shrimp genome. These results suggest that the kuruma shrimp genome has at least two types of PLEs with different domain compositions, phylogenetic positions, and probably chromosomeal localization. Such distinctive types of PLEs in an organism have never been described and hence could be a potential source to understand how multiple PLE types evolved.     

Nucleotide and repeat length variation at the nonA gene of the Drosophila virilis group species and its effects on male courtship song

Genes found to affect male courtship song characters in Drosophila melanogaster are good candidates when tracing genes responsible for species-specific songs in other Drosophila species. It has previously been shown that Thr–Gly repeat length variation at the period gene affects song traits in D. melanogaster, which gives the repetitive regions a special interest. In this work, we have characterised the patterns of nucleotide variation for gene regions containing two Gly and one Gln–Ala repeat in another D. melanogaster song gene, no-on-transient A, in D. virilis group species. The levels of nucleotide variability in D. virilis nonA were similar to those found for other genes of the species, and the gene sequences showed no signs of deviation from neutrality. The Gly 2 repeat preceding the central domain of the gene exhibited length variation, which did not, however, correlate with song variation either within D. virilis or between the species of D. virilis group. The Gly 3 repeat located on the other side of the central domain showed amino acid divergence parallel to the consensus phylogeny of the D. virilis group species. The species of the virilis subgroup having Asn after the first three glycines in this repeat have simple songs with no species-specificity, while the species of the montana subgroup having two Gly or Asn–Ser in this site have unique courtship songs. Amino acid differences between the species in this repeat may, however, reflect species phylogeny rather than have an effect on song divergence per se.     

The hobo transposable element has transposase-dependent and-independent excision activity in drosophilid species

Mobility of the hobo transposable element was determined for several strains of Drosophila melanogaster and several Drosophila species. Mobility was assessed by use of an in vivo transient assay in the soma of developing embryos, which monitored hobo excision from injected indicator plasmids. Excision was detected in a D. melanogaster strain (cn; ry ⁴²) devoid of endogenous hobo elements only after co-injection of a helper plasmid containing functional hobo transposase under either heat shock or normal promoter regulation. Excision was also detected in D. melanogaster without helper in strains known to contain genomic copies of hobo. In Drosophila species confirmed not to contain hobo, hobo excision occurred at significant rates both in the presence and absence of co-injected helper plasmid. In four of the seven species tested, excision frequencies were two- to fivefold lower in the presence of plasmid-borne hobo. hobo excision donor sites were sequenced in indicator plasmids extracted from D. melanogaster cn; ry ⁴² and D. virilis embryos. In the presence of hobo transposase, the predominant excision sites were identical in both species, having breakpoints at the hobo termini with an inverted duplication of proximal insertion site DNA. However, in the absence of hobo transposase in D. virilis, excision breakpoints were apparently random and occurred distal to the hobo termini. The data indicate that hobo is capable of functioning in the soma during embryogenesis, and that its mobility is unrestricted in drosophilids. Furthermore, drosophilids not containing hobo are able to mobilize hobo, presumably by a hobo-related cross-mobilizing system. The cross-mobilizing system in D. virilis is not functionally identical to hobo with respect to excision sequence specificity.     

Histone H1 Genes and Histone Gene Clusters in the Genus Drosophila

Whereas the genomes of many organisms contain several nonallelic types of linker histone genes, one single histone H1 type is known in Drosophila melanogaster that occurs in about 100 copies per genome. Amplification of H1 gene sequences from genomic DNA of wild type strains of D. melanogaster from Oregon, Australia, and central Africa yielded numerous clones that all exhibited restriction patterns identical to each other and to those of the known H1 gene sequence. Nucleotide sequences encoding the evolutionarily variable domains of H1 were determined in two gene copies of strain Niamey from central Africa and were found to be identical to the known H1 sequence. Most likely therefore, the translated sequences of D. melanogaster H1 genes do not exhibit intragenomic or intergenomic variations. In contrast, three different histone H1 genes were isolated from D. virilis and found to encode proteins that differ remarkably from each other and from the H1 of D. melanogaster and D. hydei. About 40 copies of H1 genes are organized in the D. virilis genome with copies of core histone genes in gene quintets that were found to be located in band 25F of chromosome 2. Another type of histone gene cluster is present in about 15 copies per genome and contains a variable intergenic sequence instead of an H1 gene. The H1 heterogeneity in D. virilis may have arisen from higher recombination rates than occur near the H1 locus in D. melanogaster and might provide a basis for formation of different chromatin subtypes. Received: 2 March 2000 / Accepted: 1 June 2000     

Detection and location of three simple sequence DNAs in polytene chromosomes from virilis group species of Drosophila

In vitro synthesized RNAs complementary to the three satellite DNAs of Drosophila virilis have been used in a series of in situ hybridization experiments with polytene chromosomes from virilis group species. Gall and Atherton (1974) demonstrated that each of the satellites of D. virilis is comprised of many repeats of a distinct, seven base pair long, simple sequence. With few exceptions, copies of each of these simple sequences are detected in the chromocenters of all virilis group species. This is true even in species which do not possess satellite DNAs at buoyant densities corresponding to those of the satellite DNAs of D. virilis. Small quantities of the three simple sequences are also detected in euchromatic arms of several different species. The same euchromatic location may contain detectable copies of one, two, or all three simple sequence DNAs. The amounts of simple sequences at each location in the euchromatin may vary between species, between different stocks of the same species, and even between individuals of the same stock. The simple sequences located in the euchromatin appear to undergo DNA replication during formation of polytene chromosomes unlike those in heterochromatin. The locations of the euchromatic sequences are not the results of single chromosomal inversion events involving heterochromatic and euchromatic breakpoints.     

Structural and functional analysis of a new retrotransposon class in Drosophila species

Several copies of the Penelope transposable element, previously described in Drosophila virilis, have been studied in different D. virilis strains and D. melanogaster strains transformed with P-based constructs bearing a full-size Penelope copy. Most Penelope copies in both species have large terminal inverted repeats (TIRs) and deletions of various sizes at the 5′ ends of their ORFs. Junctions between TIRs and ORFs usually have microhomologies of various lengths, which allowed a hypothesis explaining the emergence of these complex structures at the molecular level to be put forward. Most Penelope copies have a 34 bp long direct repeat at the ORF ends. Full-size and truncated Penelope copies are usually surrounded by target site duplications of various lengths.     

Aurora, a non-mobile retrotransposon in Drosophila melanogaster heterochromatin

A novel retrotransposon, aurora, containing 324 by long terminal repeats (LTRs) was detected in Drosophila melanogaster as a 5 kb insertion in the heterochromatic Stellate gene. This insertion causes a 5 bp duplication of the integration site. Southern analysis and in situ hybridization data show that all detectable copies of aurora are immobilized in the D. melanogaster heterochromatin. However, mobile copies of aurora were revealed in the cuchromatin of D. simulans. The element was also found in various species of the melanogaster subgroup and in the D. virilis genome.The nucleotide sequence data reported in this paper will appear in the EMBL, GenBank and DDBJ Nucleotide Sequence Databases under the accession numbers X70361 and X70362     

Evolution and arrangement of the <Emphasis Type="Italic">hsp70</Emphasis> gene cluster in two closely related species of the <Emphasis Type="Italic">virilis</Emphasis> group of <Emphasis Type="Italic">Drosophila</Emphasis>

To investigate the genetic basis of differing thermotolerance in the closely related species Drosophila virilis and Drosophila lummei, which replace one another along a latitudinal cline, we characterized the hsp70 gene cluster in multiple strains of both species. In both species, all hsp70 copies cluster in a single chromosomal locus, 29C1, and each cluster includes two hsp70 genes arranged as an inverted pair, the ancestral condition. The total number of hsp70 copies is maximally seven in the more thermotolerant D. virilis and five in the less tolerant D. lummei, with some strains of each species exhibiting lower copy numbers. Thus, maximum hsp70 copy number corresponds to hsp70 mRNA and Hsp70 protein levels reported previously and the size of heat-induced puffs at 29C1. The nucleotide sequence and spacing of the hsp70 copies are consistent with tandem duplication of the hsp70 genes in a common ancestor of D. virilis and D. lummei followed by loss of hsp70 genes in D. lummei. These and other data for hsp70 in Drosophila suggest that evolutionary adaptation has repeatedly modified hsp70 copy number by several different genetic mechanisms.     

Diversity in the Integrase Coding Domain of a gypsy-like Retrotransposon Among Wild Relatives of Rice in the Oryza officinalis Complex

The Oryza officinalis complex is a genetically diverse, tertiary genepool of rice. We analyzed part of the primary structure of the integrase coding domain (ICD) of a gypsy-like retrotransposon from species of the O. officinalis species complex. PCR was performed with degenerate primers that hybridized to conserved sequences in the integrase genes of gypsy-type retrotransposons, using total DNA from different species of the O. officinalis complex as templates. Cloning and sequencing of the PCR products showed that the amplified fragments are highly homologous to each other (75–90%) and belong to one family of retrotransposons that is related to the previously studied RIRE-2 element from rice. Two main subfamilies of 292 and 351 bp were distinguished. Analysis of primary sequence data supports previous reports that sequence divergence during vertical transmission has been the major influence on the evolution of gypsy-type retrotransposons in Oryza species. Based on sequence data phylogenetic relationships among species of the O. officinalis complex were estimated. The data suggests that O. eichingeri is more closely related to the ancestral species of the complex. This revised version was published online in July 2006 with corrections to the Cover Date.     

Gypsy homologous sequences inDrosophila subobscura (gypsyDS)

Summary Characterization of sequences homologous to theDrosophila melanogaster gypsy transposable element was carried out inDrosophila subobscura (gypsyDS). They were found to be widely distributed among natural populations of this species. From Southern blot and in situ analyses, these sequences appear to be mobile in this species.GypsyDS sequences are located in both euchromatic and heterochromatic regions. A completegypsyDS sequence was isolated from aD. subobscura genomic library, and a 1.3-kb fragment which aligns with the ORF2 of theD. melanogaster gypsy element was sequenced. Comparisons of this sequence in three species (D. subobscura, D. melanogaster, and D. virilis) indicate that there is greater similarity between theD. subobscura-D. virilis sequences than betweenD. subobscura andD. melanogaster. Molecular divergence ofgypsy sequences betweenD. virilis andD. subobscura is estimated at 16 MY, whereas the most likely divergence time of these two species is more than 60 MY. These data strongly suggest thatgypsy sequences have been horizontally transferred between these species.Offprint requests to: T.M. Alberola     

Cloning and molecular characterization of a novel x-type glutenin gene Glu-D(t)1 from Aegilops tauschii

A novel gene encoding an x-type high molecular weight glutenin subunit (HMW-GS), designated 1Dx1.1 ^t , was isolated from Aegilops tauschii. It is the largest HMW-GS gene reported so far in this species and its product has a slower mobility than that of subunit 1Ax1 in SDS-PAGE. The open reading frame (ORF) of the gene was 2,628 bp, encoding a protein of 874 amino acid residues. Comparisons of amino acid sequences showed that subunit 1Dx1.1^t had high similarity with other 1Dx subunits but also had two unique characteristics. Firstly, a tripeptide of consensus LQE present in the N-terminal domains of other 1Dx subunits was absent from subunit Dx1.1^t. Secondly, three copies of tandem duplications of the tripeptide motif GQQ and a novel tripeptide sequence (GQL) were present in its central repetitive domain. Phylogenetic analysis showed that subunit 1Dx1.1^t clustered with other known 1Dx subunits.     

One- and two-dimensional polyacrylamide gel analysis of the heat shock proteins of the virilis group of Drosophila

The heat shock proteins of the virilis group of Drosophila are analyzed by one- and two-dimensional polyacrylamide gel analysis. This group consists of the two closely related but distinct virilis and montana phylads. The analysis reveals that some of the heat shock proteins are highly conserved among the two phylads while others are not. The 83-, 72-, and 69-kdalton proteins comigrate in all species examined. There is, however, a noticable trend toward greater molecular weight variability in the smaller heat shock proteins. In general, the heat shock protein patterns within each phylad follow the proposed phylogenetic relationships with some exceptions. D. ezoana and D. littoralis, both members of the montana phylad, exhibit heat shock protein patterns more similar to those of the virilis phylad. The data also demonstrate that the montana phylad has almost two times the heat shock allele members that the virilis phylad has. It is also shown that F₁ and F₂ hybrid flies of crosses between Drosophila species having different patterns of heat shock proteins show Mendelian segregation of alleles. After several generations of inbred growth, however, the pattern of heat shock protein synthesis in reciprocal hybrids each resembles that of the paternal parent. The implications of these findings are discussed.This research was supported in part by Damon Runyon-Walter Winchell Grant DRG-233F to R.M.S. and NIH Grant GM 27611 to R.V.S. R.V.S. is the recipient of an NIH Research Career Development Award.     

Evolutionary divergence of the cytochrome b5 gene of Drosophila

Cytochrome proteins perform a broad spectrum of biological functions ranging from oxidative metabolism to electron transport and are thus essential to all organisms. The b-type cytochrome proteins bind heme noncovalently, are expressed in many different forms and are localized to various cellular compartments. We report the characterization of the cytochrome b₅ (Cyt-b) gene of Drosophila virilis and compare its structure to the Cyt-b gene of Drosophila melanogaster. As in D. melanogaster, the D. virilis gene is nuclear encoded and single copy. Although the intron/exon structures of these homologues differ, the Cyt-b proteins of D. melanogaster and D. virilis are approximately 75% identical and share the same size coding regions (1,242 nucleotides) and protein products (414 amino acids). The Drosophila Cyt-b proteins show sequence similarity to other b-type cytochromes, especially in the N-terminal heme-binding domain, and may be targeted to the mitochondrial membrane. The greatest levels of similarity are observed in areas of potential importance for protein structure and function. The exon sequences of the D. virilis Cyt-b gene differ by a total of 292 base changes. However, 62% of these changes are silent. The high degree of conservation between species separated by 60 million years of evolution in both the DNA and amino acid sequences suggests this nuclear cytochrome b₅ locus encodes an essential product of the Drosophila system.Correspondence to: C.E. Rozek