The evolutionary genetics of thehobo transposable element in theDrosophila melanogaster complex

Hobo elements are a family of transposable elements found inDrosophila melanogaster and its three sibling species:D. simulans, D. mauritiana andD. sechellia. Studies inD. melanogaster have shown thathobo may be mobilized, and that the genetic effects of such mobilizations included the general features of hybrid dysgenesis: mutations, chromosomal rearrangements and gonadal dysgenis in F1 individuals. At the evolutionary level somehobo-hybridizing sequences have also been found in the other members of themelanogaster subgroup and in many members of the relatedmontium subgroup. Surveys of older collected strains ofD. melanogaster suggest that completehobo elements were absent prior to 50 years ago and that they have recently been introduced into this species by horizontal transfer. In this paper we review our findings and those of others, in order to precisely describe the geographical distribution and the evolutionary history ofhobo in theD. melanogaster complex. Studies of the DNA sequences reveal a different level of divergence between the groupD. melanogaster, D. simulans andD. mauritiana and the fourth speciesD. sechellia. The hypothesis of multiple transfers in the recent past into theD. melanogaster complex from a common outside source is discussed.     

Evolution of thedec-1 eggshell locus inDrosophila. I. Restriction site mapping and limited sequence comparison in themelanogaster species subgroup

Summary We have analyzed 18 kb of DNA in and upstream of thedefective chorion-1 (dec-1) locus of the eight known species of themelanogaster species subgroup ofDrosophila. The restriction maps ofD. simulans, D. mauritiana, D. sechellia, D. erecta, andD. orena are shown to have basically the restriction map ofD. melanogaster, whereas the maps ofD. teissieri andD. yakuba were more difficult to align. However, the basic amount of DNA and sequence arrangement appear to have been conserved in these species. A small deletion of varying length (65–200 bp) is found in a repeated sequence of the central transcribed region ofD. melanogaster, D. simulans, andD. erecta. Restriction site mapping indicated that thedec-1 gene is highly conserved in themelanogaster species subgroup. However, sequence comparison revealed that the amount of nucleotide and amino acid substitution in the repeated region is much larger than in the 5 translated region. The 5 flanking region showed noticeable restriction site polymorphisms between species. Based on calculations from the restriction maps a dendrogram was derived that supports earlier published phylogenetic relationships within themelanogaster species subgroup except that theerecta-orena pair is placed closer to themelanogaster complex than toD. teissieri andD. yakuba.     

The distribution of the transposable elementBari-1 in theDrosophila melanogaster andDrosophila simulans genomes

The distribution of the transposable elementBari-1 inD. melanogaster andD. simulans was examined by Southern blot analysis and byin situ hybridization in a large number of strains of different geographical origins and established at different times.Bari-1 copies mostly homogeneous in size and physical map are detected in all strains tested. Both inD. melanogaster and inD. simulans a relatively high level of intraspecific insertion site polymorphism is detectable, suggesting that in both speciesBari-1 is or has been actively transposing. The main difference between the two sibling species is the presence of a large tadem array of the element in a well-defined heterochromatic location of theD. melanogaster genome, whereas such a cluster is absent inD. simulans. The presence ofBari-1 elements with apparently identical physical maps in allD. melanogaster andD. simulans strains examined suggests thatBari-1 is not a recent introduction in the genome of themelanogaster complex. Structural analysis reveals unusual features that distinguish it from other inverted repeat transposons, whereas many aspects are similar to the widely distributedTc1 element ofC. elegans.     

Length polymorphism in the threonine-glycine-encoding repeat region of theperiod gene inDrosophila

Summary Single-fly polymerase chain reaction amplification and direct DNA sequencing revealed high levels of length polymorphism in the threonine-glycine encoding repeat region of theperiod (per) gene in natural populations ofDrosophila melanogaster. DNA comparison of two alleles of identical lengths gave a high number of synonymous substitutions suggesting an ancient time of separation. However detailed examination of the sequences of different Thr-Gly length variants indicated that this divergence could be understood in terms of four deletion/insertion events. InDrosophila pseudoobscura a length polymorphism is observed in a five-amino acid degenerate repeat, which corresponds tomelanogaster's Thr-Gly domain. In spite of the differences betweenD. melanogaster andD. pseudoobscura in the amino acid sequence of the repeats, the predicted secondary structures suggest evolutionary and mechanistic constraints on theper protein of these two species.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

Thehermit transposable element of the Australian sheep blowfly,Lucilia cuprina, belongs to thehAT family of transposable elements

We report the cloning ofhermit, a member of thehAT family of transposable elements from the genome of the Australian sheep blowfly,Lucilia cuprina. Hermit is 2716 bp long and is 49% homologous to the autonomoushobo element,HFL1, at the nucleic acid level.Hermit has 15 bp terminal inverted repeats that share 10 bp with the terminal inverted repeats ofHFL1. Conceptual translation reveals a 583 residue open reading frame (ORF) that is 64% similar and 42% identical to theHFL1 ORF. However, the sequence of thehermit element contains two frameshifts within the putative ORF, indication thathermit is an inactive element. Analysis ofL. cuprina strains from within and outside Australia suggested thathermit is present as a single copy in all the genomes analysed.     

A multicopy suppressor ofnin1-1 of the yeastSaccharomyces cerevisiae is a counterpart of theDrosophila melanogaster diphenol oxidase A2 gene,Dox-A2

NIN1 is an essential gene for growth of the yeastSaccharomyces cerevisiae and was recently found to encode a component of the regulatory subunit of the 26S proteasome. Thenin1-1 mutant is temperature sensitive and its main defect is in G₁/S progression and G₂/M progression at non-permissive temperatures. One of the two multicopy suppressors ofnin1-1, SUN2 (SUppressor of Nin1-1), was found to encode a protein of 523 amino acids whose sequence is similar to those ofDrosophila melanogaster diphenol oxidase A2 and the mouse mast-cell Tum^– transplantation antigen, P91A. The C-terminal half of Sun2p was found to be functional as Sun2p at 25° C, 30° C, and 34° C but not at 37° C. The open reading frame (ORF) of theDrosophila diphenol oxidase A2 gene (Dox-A2) was obtained from a lambda phage cDNA library using the polymerase chain reaction technique. TheDox-A2 ORF driven by theTDH3 promoter complemented the phenotype of a strain deleted forsun2. ThisDox-A2-dependent strain was temperature sensitive and accumulated dumb-bell-shaped cells, with an undivided nucleus at the isthmus, after temperature upshift. This morphology is similar to that ofnin1-1 cells kept at a restrictive temperature. These results suggest thatSUN2 is a functional counterpart ofDox-A2 and that these genes play a pivotal role in the cell cycle in each organism.     

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     

Identification of replication origins in the genome of the methanogenic archaeon,<Emphasis Type="Italic"> Methanocaldococcus jannaschii</Emphasis>

Methanocaldococcus jannaschii has been notorious as an archaeon in which the replication origins are difficult to identify. Although extensive efforts have been exerted on this issue, the locations of replication origins still remain elusive 7 years after the publication of its complete genome sequence in 1996. Ambiguous results were obtained in identifying the replication origins of M. jannaschii based on all theoretical and experimental approaches. In the genome of M. jannaschii, we found that an ORF (MJ0774), annotated as a hypothetical protein, is a homologue of the Cdc6 protein. The position of the gene is at a global minimum of the x component of the Z curve, i.e., RY disparity curve, which has been used to identify replication origins in other Archaea. In addition, an intergenic region (694,540–695,226 bp) that is between the cdc6 gene and an adjacent ORF shows almost all the characteristics of known replication origins, i.e., it is highly rich in AT composition (80%) and contains multiple copies of repeat elements and AT stretches. Therefore, these lines of evidence strongly suggest that the identified region is a replication origin, which is designated as oriC1. The analysis of the y component of the Z curve, i.e., MK disparity curve, suggests the presence of another replication origin corresponding to one of the peaks in the MK disparity curve at around 1,388 kb of the genome.Communicated by G. Antranikian     

Cloning and characterization of thewhite andtopaz eye color genes from the sheep blowflylucilia cuprina

Summary Clones carrying thewhite andtopaz eye color genes have been isolated from genomic DNA libraries of the blowflyLucilia cuprina using cloned DNA from the homologouswhite andscarlet genes. respectively, ofDrosophila melanogaster as probes. On the basis of hybridization studies using adjacent restriction fragments, homologous fragments were found to be colinear between the genes from the two species. The nucleotide sequence of a short region of thewhite gene ofL. cuprina has been determined, and the homology to the corresponding region ofD. melanogaster is 72%; at the derived amino acid level the homology is greater (84%) due to a marked difference in codon usage between the species. A major difference in genome organization between the two species is that whereas the DNA encompassing theD. melanogaster genes is free of repeated sequences. that encompassing theirL. cuprina counterparts contains substantial amounts of repeated sequences. This suggests that the genome ofL. cuprina is organized on the short period interspersion pattern. Repeated sequence DNA elements, which appear generally to be short (less than 1 kb) and which vary in repetitive frequency in the genome from greater than 10⁴ copies to less than 10² copies, are found in at least two different locations in the clones carrying these genes. One type of repeat structure, found by sequencing, consists of tandemly repeating short sequences. Restriction site and restriction fragment length polymorphisms involving both thewhite andtopaz gene regions are found within and between populations ofL. cuprina.     

The occurrence of the transposable elementpogo inDrosophila melanogaster

We examined the genomic occurrence of the transposable elementpogo in over 120 strains ofDrosophila melanogaster, from around the world and from different eras. All had multiple copies of a 2.1 kilobase (kb)pogo element, and multiple copies of several size classes between 1.0 and 1.8 kb. There were differences between strains in intensities or presences of deletion-derivative size classes, suggesting current or recent mobility in the species. We were unable to find anypogo-hybridization in eight other species in the genus, in three subgenera, or in the relatedScaptomyza pallida. Thepogo element may be a ‘middle-aged’ element in the genome ofD. melanogaster, having entered the species since its divergence from its sibling species, but long before theP andhobo elements.     

P-related sequences inDrosophila bifasciata: A molecular clue to the understanding of P-element evolution in the genusDrosophila

Summary Two P-elements (bif1 and bif2) were isolated from a genomic library ofDrosophila bifasciata. Both elements are internally deleted and have lost the coding capacity for a functional transposase. One of the elements (bif2) contains an insert consisting of a repetitive sequence. The terminal inverted repeats and the segments necessary for passive mobility are well conserved. Element bif2 has retained rudiments of the coding sequence of exon 0 and exon 3, but the reading frame is destroyed by insertions and deletions. The comparison of theD. bifasciata P-elements with P-elements ofDrosophila melanogaster andDrosophila nebulosa reveals that the two latter sequences are more similar to each other than either of them is to theD. bifasciata elements. This finding contradicts the phylogenetic relationship of the species and can be taken as an indirect but unequivocal evidence for recent horizontal gene transfer from a relative ofD. nebulosa to the gene pool ofD. melanogaster. The P-elements ofD. bifasciata are phylogenetically ancient and have evolved independently for about 50 million years. A higher substitution rate at the third codon position as well as a predominance of conservative replacements at the amino acid level indicates that the P-elements ofD. bifasciata have been under selective constraint over a long period and that immobilization has occurred only recently.     

On the evolution of Tn21-like multiresistance transposons: sequence analysis of the gene (aacC1) for gentamicin acetyltransferase-3-I(AAC(3)-I), another member of the Tn21-based expression cassette

Summary The aminoglycoside-3-O-acetyltransferase-I gene (aacC1) from R plasmids of two incompatibility groups (R1033 [Tn1696], and R135) was cloned and sequenced. In the case of R1033, it was shown that theaacC gene is coded by a precise insertion of 833 bp between theaadA promoter and its structural gene in a Tn21 related transposon (Tn1696). This insertion occurs at the same target sequence as that of the OXA-1 β-lactamase gene insertion in Tn2603. Upstream of theaacC gene, we found an open reading frame (ORF) which is probably implicated in the site-specific recombinational events involved in the evolution of this family of genetic elements. These results provide additional confirmation of the role of Tn21 elements as naturally occurring interspecific transposition and expression casssettes.     

Invasion of thehobo transposable element studied byin situ hybridization on polytene chromosomes ofDrosophila melanogaster

The invasion kinetics ofhobo transposable element in theDrosophila melanogaster genome was studied byin situ hybridization on the polytene chromosomes. Six independent lines ofDrosophila melanogaster flies that had been previously transformed by microinjection of the pHFL1 plasmid containing a completehobo element were followed over 50 generations. We observed thathobo elements were scattered on each of the chromosome arms, with more insertion sites on the 3R arm. The total number of insertion sites remains quite small, between four and six, at generation 52. On the 2R arm, a short inversion appeared once at generation 52. Most of the integration sites reported here were already described for several transposons but some of them appear to be hotspots forhobo elements.     

IS<Emphasis Type="Italic">Rm31</Emphasis>, a new insertion sequence of the IS<Emphasis Type="Italic">66</Emphasis> family in<Emphasis Type="Italic"> Sinorhizobium meliloti</Emphasis>

Sinorhizobium meliloti natural populations show a high level of genetic polymorphism possibly due to the presence of mobile genetic elements such as insertion sequences (IS), transposons, and bacterial mobile introns. The analysis of the DNA sequence polymorphism of the nod region of S. meliloti pSymA megaplasmid in an Italian isolate led to the discovery of a new insertion sequence, ISRm31. ISRm31 is 2,803 bp long and has 22-bp-long terminal inverted repeat sequences, 8-bp direct repeat sequences generated by transposition, and three ORFs (A, B, C) coding for proteins of 124, 115, and 541 amino acids, respectively. ORF A and ORF C are significantly similar to members of the transposase family. Amino acid and nucleotide sequences indicate that ISRm31 is a member of the IS66 family. ISRm31 sequences were found in 30.5% of the Italian strains analyzed, and were also present in several collection strains of the Rhizobiaceae family, including S. meliloti strain 1021. Alignment of targets sites in the genome of strains carrying ISRm31 suggested that ISRm31 inserts randomly into S. meliloti genomes. Moreover, analysis of ISRm31 insertion sites revealed DNA sequences not present in the recently sequenced S. meliloti strain 1021 genome. In fact, ISRm31 was in some cases linked to DNA fragments homologous to sequences found in other rhizobia species.     

The use of a non-LTR element to date the formation of the Sdic gene cluster

Transposable elements comprise a considerable part of eukaryotic genomes, and there is increasing evidence for their role in the evolution of genomes. The number of active transposable elements present in the host genome at any given time is probably small relative to the number of elements that no longer transpose. The elements that have lost the ability to transpose tend to evolve neutrally. For example, non-LTR retrotransposons often become 5′ truncated due to their own transposition mechanism and hence lose their ability to transpose. The resulting transposons can be characterized as “dead-on-arrival” (DOA) elements. Because they are abundant and ubiquitous, and evolve neutrally in the location where they were inserted, these DOA non-LTR elements make a useful tool to date molecular events. There are four copies of a “dead-on-arrival” RT1C element on the recently formed Sdic gene cluster of Drosophila melanogaster, that are not present in the equivalent region of the other species of the melanogaster subgroup. The life history of the RT1C elements in the genome of D. melanogaster was used to determine the insertion chronology of the elements in the cluster and to date the duplication events that originated this cluster.     

Drosophila virilis histone gene clusters lacking H1 coding segments

Summary Approximately 30–40% ofDrosophila virilis DNA complementary to clonedDrosophila histone genes is reduced to 3.4-kilobase-pair (kbp) segments by Bgl I or Bgl II digestion. The core histone genes of a 3.4-kbp Bgl II segment cloned in the plasmid pDv3/3.4 have the same order as theD. melanogaster core histone genes in the plasmid cDm500: . Nonetheless, pDv3/3.4 and cDm500 have different histone gene configurations: In pDv3/3.4, the region between the H2B and H3 genes contains 0.35 kbp and cannot encode histone H1; in cDm500, the region contains 2.0 kbp and encodes histone H1. The lack of an H1 gene between the H2B and H3 genes in 30–40% ofD. virilis histone gene clusters suggests that changes in histone gene arrays have occurred during the evolution ofDrosophila. The ancestors of modernDrosophila may have possessed multiple varieties of histone gene clusters, which were subsequently lost differentially in thevirilis andmelanogaster lineages. Alternatively, they may have possessed a single variety, which was rearranged during evolution. The H1 genes ofD. virilis andD. melanogaster did not cross-hybridize in vitro under conditions that maintain stable duplexes between DNAs that are 75% homologous. Consequently,D. virilis H1 genes could not be visualized by hybridization to an H1-specific probe and thus remain unidentified. Our observations suggest that the coding segments in the H1 genes ofD. virilis andD. melanogaster are >25% divergent. Our estimate of sequence divergence in the H1 genes ofD. virilis andD. melanogaster seems high until one considers that the coding sequences of cloned H1 genes from the closely related speciesD. melanogaster andD. simulans are 5% divergent.     

The bz-rcy allele of the Cy transposable element system of Zea mays contains a Mu-like element insertion

Summary The receptive component of theCy transposable element system (rcy: Mu7) at theBz locus ofZea mays L. is 2.2 kb and has long terminal inverted repeats. The insertion is flanked by a 9 bp duplication. In the presence of an autonomousCy element in the genome,rcy: Mu7 is excised frombz-rcy in a manner consistent with a model suggested previously. The termini ofrcy: Mu7 have 85% sequence similarity with theMu1 element ofZ. mays. This is consistent with the observation thatMu1 can behave genetically like a receptive component of theCy system.     

rtpA, a gene encoding a bacterial two-component sensor kinase, determines pathogenic traits ofPseudomonas tolaasii, the causal agent of brown blotch disease of a cultivated mushroom,Pleurotus ostreatus

Pseudomonas tolaasii strain PT814 produces extracellular toxins, tolaasins, and a volatile toxin, tovsin, that are responsible for the induction of brown blotch and rotting, respectively, in a cultivated mushroom,Pleurotus ostreatus. Insertions of single transposon mini-Tn5Km 1 into the chromosome ofP. tolaasii strain PT814 generated mutants that are pleiotropically defective in tolaasin and protease production, and altered in colony morphology. The mutants, however, produce tovsin at the level of wild-type. Variants phenotypically similar to the pleiotropic mutants ofP. tolaasii strain PT814 spontaneously occurred inP. tolaasii strain S8501 at 22–30°C in vitro. The occurrence of variants was significantly reduced in the presence of extracts ofP ostreatus or at a temperature of 15–20°C. ThertpA gene (rtpA=regulator gene of tolaasin production and other pleiotropic traits) isolated from aP. tolaasii strain PT814 gene library restored the wild-type phenotype in both the mini-Tn5km 1 insertion and spontaneous mutants. mini-Tn5km 1 insertions were also located in the allele ofrtpA. Nucleotide sequencing of thertpA DNA revealed an open reading frame of 2,751 bp predicted to encode a protein consisting of 917 amino acid residues with a molecular mass of 100.6 kDa and displaying the conserved amino acid sequence of both sensor, and receiver domains of “bacterial two-component regulators”. The data suggest that the machinery responding to environmental stimuli is essential for the pathogenic interaction ofP. tolaasii with the mushroom.     

Organization and possible origin of the Bari-1 cluster in the heterochromatic h39 region of Drosophila melanogaster

The molecular organization of the heterochromatic h39 region of the Drosophila melanogaster second chromosome has been investigated by studying two BAC clones identified both by Southern blotting and by FISH experiments as containing tandem arrays of Bari1, a transposable element present only in this region. Such BAC clones appear to contain different portions of the h39 region since they differ in the DNA sequences flanking the Bari1 repeats on both sides. Thus, the 80 Bari1 copies estimated to be present in the h39 region are split into at least two separated subregions. On the basis of the analysis of the flanking sequences a possible mechanism depending on an aberrant activity of the Bari1 transposase is proposed for the genesis of the heterochromatic tandem arrays of the element.