Evolution of the Gypsy endogenous retrovirus in the Drosophila melanogaster subgroup

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

Tirant Stealthily Invaded Natural Drosophila melanogaster Populations during the Last Century

It was long thought that solely three different transposable elements (TEs)-the I-element, the P-element, and hobo-invaded natural Drosophila melanogaster populations within the last century. By sequencing the "living fossils" of Drosophila research, that is, D. melanogaster strains sampled from natural populations at different time points, we show that a fourth TE, Tirant, invaded D. melanogaster populations during the past century. Tirant likely spread in D. melanogaster populations around 1938, followed by the I-element, hobo, and, lastly, the P-element. In addition to the recent insertions of the canonical Tirant, D. melanogaster strains harbor degraded Tirant sequences in the heterochromatin which are likely due to an ancient invasion, likely predating the split of D. melanogaster and D. simulans. These degraded insertions produce distinct piRNAs that were unable to prevent the novel Tirant invasion. In contrast to the I-element, P-element, and hobo, we did not find that Tirant induces any hybrid dysgenesis symptoms. This absence of apparent phenotypic effects may explain the late discovery of the Tirant invasion. Recent Tirant insertions were found in all investigated natural populations. Populations from Tasmania carry distinct Tirant sequences, likely due to a founder effect. By investigating the TE composition of natural populations and strains sampled at different time points, insertion site polymorphisms, piRNAs, and phenotypic effects, we provide a comprehensive study of a natural TE invasion.

A new approach for detecting adventitious viruses shows Sf-rhabdovirus-negative Sf-RVN cells are suitable for safe biologicals production

Background

Adventitious viral contamination in cell substrates used for biologicals production is a major safety concern. A powerful new approach that can be used to identify adventitious viruses is a combination of bioinformatics tools with massively parallel sequencing technology. Typically, this involves mapping or BLASTN searching individual reads against viral nucleotide databases. Although extremely sensitive for known viruses, this approach can easily miss viruses that are too dissimilar to viruses in the database. Moreover, it is computationally intensive and requires reference cell genome databases. To avoid these drawbacks, we set out to develop an alternative approach. We reasoned that searching genome and transcriptome assemblies for adventitious viral contaminants using TBLASTN with a compact viral protein database covering extant viral diversity as the query could be fast and sensitive without a requirement for high performance computing hardware.

Results

We tested our approach on Spodoptera frugiperda Sf-RVN, a recently isolated insect cell line, to determine if it was contaminated with one or more adventitious viruses. We used Illumina reads to assemble the Sf-RVN genome and transcriptome and searched them for adventitious viral contaminants using TBLASTN with our viral protein database. We found no evidence of viral contamination, which was substantiated by the fact that our searches otherwise identified diverse sequences encoding virus-like proteins. These sequences included Maverick, R1 LINE, and errantivirus transposons, all of which are common in insect genomes. We also identified previously described as well as novel endogenous viral elements similar to ORFs encoded by diverse insect viruses.

Conclusions

Our results demonstrate TBLASTN searching massively parallel sequencing (MPS) assemblies with a compact, manually curated viral protein database is more sensitive for adventitious virus detection than BLASTN, as we identified various sequences that encoded virus-like proteins, but had no similarity to viral sequences at the nucleotide level. Moreover, searches were fast without requiring high performance computing hardware. Our study also documents the enhanced biosafety profile of Sf-RVN as compared to other Sf cell lines, and supports the notion that Sf-RVN is highly suitable for the production of safe biologicals.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12896-017-0412-z) contains supplementary material, which is available to authorized users.

Wolbachia influences the maternal transmission of the gypsy endogenous retrovirus in Drosophila melanogaster

The endosymbiotic bacteria of the genus Wolbachia are present in most insects and are maternally transmitted through the germline. Moreover, these intracellular bacteria exert antiviral activity against insect RNA viruses, as in Drosophila melanogaster, which could explain the prevalence of Wolbachia bacteria in natural populations. Wolbachia is maternally transmitted in D. melanogaster through a mechanism that involves distribution at the posterior pole of mature oocytes and then incorporation into the pole cells of the embryos. In parallel, maternal transmission of several endogenous retroviruses is well documented in D. melanogaster. Notably, gypsy retrovirus is expressed in permissive follicle cells and transferred to the oocyte and then to the offspring by integrating into their genomes. Here, we show that the presence of Wolbachia wMel reduces the rate of gypsy insertion into the ovo gene. However, the presence of Wolbachia does not modify the expression levels of gypsy RNA and envelope glycoprotein from either permissive or restrictive ovaries. Moreover, Wolbachia affects the pattern of distribution of the retroviral particles and the gypsy envelope protein in permissive follicle cells. Altogether, our results enlarge the knowledge of the antiviral activity of Wolbachia to include reducing the maternal transmission of endogenous retroviruses in D. melanogaster.

Animals have established complex relationships with bacteria and viruses that spread horizontally among individuals or are vertically transmitted, i.e., from parents to offspring. It is well established that members of the genus Wolbachia, maternally inherited symbiotic bacteria present mainly in arthropods, reduce the replication of several RNA viruses transmitted horizontally. Here, we demonstrate for the first time that Wolbachia diminishes the maternal transmission of gypsy, an endogenous retrovirus in Drosophila melanogaster. We hypothesize that gypsy cannot efficiently integrate into the germ cells of offspring during embryonic development in the presence of Wolbachia because both are competitors for localization to the posterior pole of the egg. More generally, it would be of interest to analyze the influence of Wolbachia on vertically transmitted exogenous viruses, such as some arboviruses.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

A genomic perspective to assessing quality of mass-reared SIT flies used in Mediterranean fruit fly (Ceratitis capitata) eradication in California

Background

Temperature sensitive lethal (tsl) mutants of the tephritid C. capitata are used extensively in control programs involving sterile insect technique in California. These flies are artificially reared and treated with ionizing radiation to render males sterile for further release en masse into the field to compete with wild males and disrupt establishment of invasive populations. Recent research suggests establishment of C. capitata in California, despite the fact that over 250 million sterile flies are released weekly as part of the state’s preventative program. In this project, genome-level quality assessment was performed, measured as expression differences between the Vienna-7 tsl mutants used in SIT programs and wild flies. RNA-seq was performed to provide a genome-wide map of the messenger RNA populations in C. capitata, and to investigate significant expression changes in Vienna-7 mass reared flies.

Results

Flies from the Vienna-7 colony showed a markedly reduced abundance of transcripts related to visual and chemical responses, including light stimuli, neural development and signaling pathways when compared to wild flies. In addition, genes associated with muscle development and locomotion were shown to be reduced. This suggests that the Vienna-7 line may be less competitive in mating and host plant finding where these stimuli are utilized. Irradiated flies showed several transcripts representing stress associated with irradiation.

Conclusions

There are significant changes at the transcriptome level that likely alter the competitiveness of mass reared flies and provide justification for pursuing methods for strain improvement, increasing competitiveness of mass-reared flies, or exploring alternative SIT approaches to increase the efficiency of eradication programs.

The rhizome of life: what about metazoa?

The increase in huge number of genomic sequences in recent years has contributed to various genetic events such as horizontal gene transfer (HGT), gene duplication and hybridization of species. Among them HGT has played an important role in the genome evolution and was believed to occur only in Bacterial and Archaeal genomes. As a result, genomes were found to be chimeric and the evolution of life was represented in different forms such as forests, networks and species evolution was described more like a rhizome, rather than a tree. However, in the last few years, HGT has also been evidenced in other group such as metazoa (for example in root-knot nematodes, bdelloid rotifers and mammals). In addition to HGT, other genetic events such as transfer by retrotransposons and hybridization between more closely related lineages are also well established. Therefore, in the light of such genetic events, whether the evolution of metazoa exists in the form of a tree, network or rhizome is highly questionable and needs to be determined. In the current review, we will focus on the role of HGT, retrotransposons and hybridization in the metazoan evolution.

Structure-based modeling of the functional HIV-1 intasome and its inhibition

The intasome is the basic recombination unit of retroviral integration, comprising the integrase protein and the ends of the viral DNA made by reverse transcription. Clinical inhibitors preferentially target the DNA-bound form of integrase as compared with the free protein, highlighting the critical requirement for detailed understanding of HIV-1 intasome structure and function. Although previous biochemical studies identified integrase residues that contact the DNA, structural details of protein-protein and protein-DNA interactions within the functional intasome were lacking. The recent crystal structure of the prototype foamy virus (PFV) integrase-viral DNA complex revealed numerous details of this related integration machine. Structures of drug-bound PFV intasomes moreover elucidated the mechanism of inhibitor action. Herein we present a model for the HIV-1 intasome assembled using the PFV structure as template. Our results pinpoint previously identified protein-DNA contacts within the quaternary structure and reveal hitherto unknown roles for Arg20 and Lys266 in DNA binding and integrase function. Models for clinical inhibitors bound at the HIV-1 integrase active site were also constructed and compared with previous studies. Our findings highlight the structural basis for HIV-1 integration and define the mechanism of its inhibition, which should help in formulating new drugs to inhibit viruses resistant to first-in-class compounds.

Evolution of Tom, 297, 17.6 and rover retrotransposons in Drosophilidae species

LTR retrotransposons are the most abundant transposable elements in Drosophila and are believed to have contributed significantly to genome evolution. Different reports have shown that many LTR retrotransposon families in Drosophila melanogaster emerged from recent evolutionary episodes of transpositional activity. To contribute to the knowledge of the evolutionary history of Drosophila LTR retrotransposons and the mechanisms that control their abundance, distribution and diversity, we conducted analyses of four related families of LTR retrotransposons, 297, 17.6, rover and Tom. Our results show that these elements seem to be restricted to species from the D. melanogaster group, except for 17.6, which is also present in D. virilis and D. mojavensis. Genetic divergences and phylogenetic analyses of a 1-kb fragment region of the pol gene illustrate that the evolutionary dynamics of Tom, 297, 17.6 and rover retrotransposons are similar in several aspects, such as low codon bias, the action of purifying selection and phylogenies that are incongruent with those of the host species. We found an extremely complex association among the retrotransposon sequences, indicating that different processes shaped the evolutionary history of these elements, and we detected a very high number of possible horizontal transfer events, corroborating the importance of lateral transmission in the evolution and maintenance of LTR retrotransposons.

Drosophila viruses and the study of antiviral host-defense

The fruit fly Drosophila melanogaster is a powerful model to study host-pathogen interactions. Most studies so far have focused on extracellular pathogens such as bacteria and fungi. More recently, viruses have come to the front, and RNA interference was shown to play a critical role in the control of viral infections in drosophila. We review here our current knowledge on drosophila viruses. A diverse set of RNA viruses belonging to several families (Rhabdoviridae, Dicistroviridae, Birnaviridae, Reoviridae, Errantiviridae) has been reported in D. melanogaster. By contrast, no DNA virus has been recovered up to now. The drosophila viruses represent powerful tools to study virus-cell interactions in vivo. Analysis of the literature however reveals that for many of them, important gaps exist in our understanding of their replication cycle, genome organization, morphology or pathogenesis. The data obtained in the past few years on antiviral defense mechanisms in drosophila, which point to evolutionary conserved pathways, highlight the potential of the D. melanogaster model to study antiviral innate immunity and to better understand the complex interaction between arthropod-borne viruses and their insect vectors.

Novel clades of chromodomain-containing Gypsy LTR retrotransposons from mosses (Bryophyta)

Retrotransposons are the major component of plant genomes. Chromodomain-containing Gypsy long terminal repeat (LTR) retrotransposons are widely distributed in eukaryotes. Four distinct clades of chromodomain-containing Gypsy retroelements are known from the vascular plants: Reina, CRM, Galadriel and Tekay. At the same time, almost nothing is known about the repertoire of LTR retrotransposons in bryophyte genomes. We have combined a search of chromodomain-containing Gypsy retroelements in Physcomitrella genomic sequences and an experimental investigation of diverse moss species. The computer-based mining of the chromodomain-containing LTR retrotransposons allowed us to describe four different elements from Physcomitrella. Four novel clades were identified that are evolutionarily distinct from the chromodomain-containing Gypsy LTR retrotransposons of other plants.

Gypsy endogenous retrovirus maintains potential infectivity in several species of Drosophilids

Background

Sequences homologous to the gypsy retroelement from Drosophila melanogaster are widely distributed among drosophilids. The structure of gypsy includes an open reading frame resembling the retroviral gene env, which is responsible for the infectious properties of retroviruses.

Results

In this study we report molecular and phylogeny analysis of the complete env gene from ten species of the obscura group of the genus Drosophila and one species from the genus Scaptomyza.

Conclusion

The results indicate that in most cases env sequences could produce a functional Env protein and therefore maintain the infectious capability of gypsy in these species.

Transcriptome Analysis of ESTs from a Chaetognath Reveals a Deep-Branching Clade of Retrovirus-Like Retrotransposons

Chaetognaths constitute a small marine phylum exhibiting several characteristic which are highly unusual in animal genomes, including two classes of both rRNA and protein ribosomal genes. As in this phylum presence of retrovirus-like elements has never been documented, analysis of a published expressed sequence tag (EST) collection of the chaetognath Spadella cephaloptera has been made. Twelve sequences representing transcript sections of reverse transcriptase domain of active retrotransposons were isolated from~11,000 ESTs. Five of them are originated from Gypsy retrovirus-like elements, whereas the other are transcripts from a Bel-Pao LTR-retrotransposon, a Penelope-like element and LINE retrotransposons. Moreover, a part of a putative integrase has also been found. Phylogenetic analyses suggest a deep-branching clade of the retrovirus-like elements, which is in agreement with the probably Cambrian origin of the phylum. Moreover, retrotransposons have not been found in telomeric-like transcripts which are probably constituted by both vertebrate and arthropod canonical repeats.

Multiple invasions of Errantivirus in the genus Drosophila

Aiming to contribute to the knowledge of the evolutionary history of Errantivirus, a phylogenetic analysis of the env gene sequences of Errantivirus gypsy, gtwin, gypsy2, gypsy3, gypsy4 and gypsy6 was carried out in 33 Drosophilidae species. Most sequences were obtained from in silico searches in the Drosophila genomes. The complex evolutionary pattern reported by other authors for the gypsy retroelement was also observed in the present study, including vertical transmission, ancestral polymorphism, stochastic loss and horizontal transfer. Moreover, the elements gypsy2, gypsy3, gypsy4 and gypsy6 were shown to have followed an evolutionary model that is similar to gypsy. Fifteen new possible cases of horizontal transfer were suggested. The infectious potential of these elements may help elucidate the evolutionary scenario described in the present study.

Polymorphism of canonical and noncanonical gypsy sequences in different species of Drosophila melanogaster subgroup: possible evolutionary relations

Mobile genetic elements constitute a substantial part of eukaryotic genome and play an important role in its organization and functioning. Co-evolution of retrotransposons and their hosts resulted in the establishment of control systems employing mechanisms of RNA interference that seem to be impossible to evade. However, "active" copies of endogenous retrovirus gypsy escape cellular control in some cases, while its evolutionary elder "inactive" variants do not. To clarify the evolutionary relationship between "active" and "inactive" gypsy we combined two approaches: the analysis of gypsy sequences, isolated from G32 Drosophila melanogaster strain and from different Drosophila species of the melanogaster subgroup, as well as the study of databases, available on the Internet. No signs of "intermediate" (between "active" and "inactive") gypsy form were found in GenBank, and four full-size G32 gypsy copies demonstrated a convergence that presumably involves gene conversion. No "active" gypsy were revealed among PCR generated gypsy ORF3 sequences from the various Drosophila species indicating that "active" gypsy appeared in some population of D. melanogaster and then started to spread out. Analysis of sequences flanking gypsy variants in G32 revealed their predominantly heterochromatic location. Discrepancy between the structure of actual gypsy sites in G32 and corresponding sequences in database might indicate significant inter-strain heterochromatin diversity.

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.

CR1 clade of non-LTR retrotransposons from Maculinea butterflies (Lepidoptera: Lycaenidae): evidence for recent horizontal transmission

Background

Non-long terminal repeat (non-LTR) retrotransposons are mobile genetic elements that propagate themselves by reverse transcription of an RNA intermediate. Non-LTR retrotransposons are known to evolve mainly via vertical transmission and random loss. Horizontal transmission is believed to be a very rare event in non-LTR retrotransposons. Our knowledge of distribution and diversity of insect non-LTR retrotransposons is limited to a few species – mainly model organisms such as dipteran genera Drosophila, Anopheles, and Aedes. However, diversity of non-LTR retroelements in arthropods seems to be much richer. The present study extends the analysis of non-LTR retroelements to CR1 clade from four butterfly species of genus Maculinea (Lepidoptera: Lycaenidae).

The lycaenid genus Maculinea, the object of interest for evolutionary biologists and also a model group for European biodiversity studies, possesses a unique, specialized myrmecophilous lifestyle at larval stage. Their caterpillars, after three weeks of phytophagous life on specific food plants drop to the ground where they are adopted to the ant nest by Myrmica foraging workers.

Results

We found that the genome of Maculinea butterflies contains multiple CR1 lineages of non-LTR retrotransposons, including those from MacCR1A, MacCR1B and T1Q families. A comparative analysis of RT nucleotide sequences demonstrated an extremely high similarity among elements both in interspecific and intraspecific comparisons. CR1A-like elements were found only in family Lycaenidae. In contrast, MacCR1B lineage clones were extremely similar to CR1B non-LTR retrotransposons from Bombycidae moths: silkworm Bombyx mori and Oberthueria caeca.

Conclusion

The degree of coding sequence similarity of the studied elements, their discontinuous distribution, and results of divergence-versus-age analysis make it highly unlikely that these sequences diverged at the same time as their host taxa. The only reasonable alternative explanation is horizontal transfer. In addition, phylogenetic markers for population analysis of Maculinea could be developed based on the described non-LTR retrotransposons.

Evolutionary pathways of the tirant LTR retrotransposon in the Drosophila melanogaster subgroup of species

Tirant, a LTR retrotransposon with copies scattered over the chromosome arms of Drosophila melanogaster, is in the process of being lost from the chromosome arms of most natural populations of the sister species D. simulans. In an attempt to clarify the dynamics and evolution of tirant, we have studied the regulatory and reverse transcriptase regions in copies of the nine closely related species of the D. melanogaster subgroup. We show that tirant is mainly vertically transmitted in these species, with the exception of a horizontal transfer event from an ancestor of D. melanogaster to D. teissieri. We propose that, in four of the species (D. melanogaster, D. simulans, D. sechellia, and D. mauritiana), the observed patterns of evolution of the regulatory region vary with genome constraints and with the history and biogeography of the species.

Analyzing patterns of microbial evolution using the mauve genome alignment system

During the course of evolution, genomes can undergo large-scale mutation events such as rearrangement and lateral transfer. Such mutations can result in significant variations in gene order and gene content among otherwise closely related organisms. The Mauve genome alignment system can successfully identify such rearrangement and lateral transfer events in comparisons of multiple microbial genomes even under high levels of recombination. This chapter outlines the main features of Mauve and provides examples that describe how to use Mauve to conduct a rigorous multiple genome comparison and study evolutionary patterns.

Evolutionary pattern of the gtwin retrotransposon in the Drosophila melanogaster subgroup

The gtwin retrotransposon was recently discovered in the Drosophila melanogaster genome and it is evolutionarily closer to gypsy endogenous retrovirus. This study has identified gtwin homologous sequences in the genome of D. simulans, D. sechellia, D. erecta and D. yakuba by performing homology searches against the public genome database of Drosophila species. The phylogenetic analyses of the gtwin env gene sequences of these species have shown some incongruities with the host species phylogeny, suggesting some horizontal transfer events for this retroelement. Moreover, we reported the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster. The results suggest that the gtwin element may be an infectious retrovirus able to invade the genome of new species, supporting the gtwin evolutionary picture shown in this work.

Molecular phylogeny of the Drosophila melanogaster species subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba-D. teissieri and D. erecta-D. orena species pairs relative to the D. melanogaster-D. simulans-D. mauritiana-D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed "vectorette PCR" to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba-D. teissieri and D. erecta-D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica-D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models.

Complex evolution of gypsy in Drosophilid species

In an endeavor to contribute to the comprehension of the evolution of transposable elements (TEs) in the genome of host species, we investigated the phylogenetic relationships of sequences homologous to the retrotransposon gypsy of Drosophila melanogaster in 19 species of Drosophila, in Scaptodrosophila latifasciaeformis, and in Zaprionus indianus. This phylogenetic study was based on approximately 500 base pairs of the env gene. Our analyses showed considerable discrepancy between the phylogeny of gypsy elements and the relationship of their host species, and they allow us to infer a complex evolutionary pattern that could include ancestral polymorphism, vertical transmission, and several cases of horizontal transmission.

Sequence divergence within transposable element families in the Drosophila melanogaster genome

The availability of the sequenced Drosophila melanogaster genome provides an opportunity to study sequence variation between copies within transposable element families. In this study,we analyzed the 624 copies of 22 transposable element (TE) families (14 LTR retrotransposons, five non-LTR retrotransposons, and three transposons). LTR and non-LTR retrotransposons possessed far fewer divergent elements than the transposons,suggesting that the difference depends on the transposition mechanism. However,there was not a continuous range of divergence of the copies in each class,which were either very similar to the canonical elements,or very divergent from them. This sequence homogeneity among TE family copies matches the theoretical models of the dynamics of these repeated sequences. The sequenced Drosophila genome thus appears to be composed of a mixture of TEs that are still active and of ancient relics that have degenerated and the distribution of which along the chromosomes results from natural selection. This clearly demonstrates that the TEs are highly active within the genome,suggesting that the genetic variability of the Drosophila genome is still being renewed by the action of TEs.

Drosophila germline invasion by the endogenous retrovirus gypsy: involvement of the viral env gene

The endogenous retrovirus gypsy is expressed at high levels in mutant flamenco female flies. Gypsy viral particles extracted from such flies can infect naive flamenco individuals raised in the presence of these extracts mixed into their food. This results in the integration of new proviruses into the germline genome. These proviruses can then increase their copy number by (1) expression in the flamenco female somatic cells, (2) transfer into the oocyte and (3) integration into the genome of the progeny. Surprisingly, unlike the infection observed in the feeding experiments, this strategy of endogenous proviral multiplication does not seem to involve the expression of the viral env gene.

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.

Abundance, distribution, and transcriptional activity of repetitive elements in the maize genome

Long terminal repeat (LTR) retrotransposons have been shown to make up much of the maize genome. Although these elements are known to be prevalent in plant genomes of a middle-to-large size, little information is available on the relative proportions composed by specific families of elements in a single genome. We sequenced a library of randomly sheared genomic DNA from maize to characterize this genome. BLAST analysis of these sequences demonstrated that the maize genome is composed of diverse sequences that represent numerous families of retrotransposons. The largest families contain the previously described elements Huck, Ji, and Opie. Approximately 5% of the sequences are predicted to encode proteins. The genomic abundance of 16 families of elements was estimated by hybridization to an array of 10,752 maize bacterial artificial chromosome (BAC) clones. Comparisons of the number of elements present on individual BACs indicated that retrotransposons are in general randomly distributed across the maize genome. A second library was constructed that was selected to contain sequences hypomethylated in the maize genome. Sequence analysis of this library indicated that retroelements abundant in the genome are poorly represented in hypomethylated regions. Fifty-six retroelement sequences corresponding to the integrase and reverse transcriptase domains were isolated from approximately 407,000 maize expressed sequence tags (ESTs). Phylogenetic analysis of these and the genomic retroelement sequences indicated that elements most abundant in the genome are less abundant at the transcript level than are more rare retrotransposons. Additional phylogenies also demonstrated that rice and maize retrotransposon families are frequently more closely related to each other than to families within the same species. An analysis of the GC content of the maize genomic library and that of maize ESTs did not support recently published data that the gene space in maize is found within a narrow GC range, but does indicate that genic sequences have a higher GC content than intergenic sequences (52% vs. 47% GC).

Evolutionary dynamics and evolutionary history in the RTE clade of non-LTR retrotransposons

This study examined the evolutionary dynamics of Bov-B LINEs in vertebrates and the evolution of the RTE clade of non-LTR retrotransposons. The first full-length reptilian Bov-B LINE element is described; it is 3.2 kb in length, with a structural organization typical of the RTE clade of non-LTR retrotransposons. The long-term evolution of Bov-B LINEs was studied in 10 species of Squamata by analysis of a PCR-amplified 1.8-kb fragment encoding part of apurinic/apyrimidinic endonuclease, the intervening domain, and the palm/fingers subdomain of reverse transcriptase. A very high level of conservation in Squamata Bov-B long interspersed nuclear elements has been found, reaching 86% identity in the nearly 600 amino acids of ORF2. The same level of conservation exists between the ancestral snake lineage and Ruminantia. Such a high level is exceptional when compared with the level of conservation observed in nuclear and mitochondrial proteins and in other transposable elements. The RTE clade has been found to be much more widely distributed than previously thought, and novel representatives have been discovered in plants, brown algae, annelids, crustaceans, mollusks, echinoderms, and teleost fishes. Evolutionary relationships in the RTE clade were deduced at the amino acid level from three separate regions of ORF2. By using different independent methods, including the divergence-versus-age analysis, several examples of horizontal transfer in the RTE clade were recognized, with important implications for the existence of HT in non-LTR retrotransposons.

Drosophila euchromatic LTR retrotransposons are much younger than the host species in which they reside

The recent release of the complete euchromatic genome sequence of Drosophila melanogaster offers a unique opportunity to explore the evolutionary history of transposable elements (TEs) within the genome of a higher eukaryote. In this report, we describe the annotation and phylogenetic comparison of 178 full-length long terminal repeat (LTR) retrotransposons from the sequenced component of the D. melanogaster genome. We report the characterization of 17 LTR retrotransposon families described previously and five newly discovered element families. Phylogenetically, these families can be divided into three distinct lineages that consist of members from the canonical Copia and Gypsy groups as well as a newly discovered third group containing BEL, mazi, and roo elements. Each family consists of members with average pairwise identities > or =99% at the nucleotide level, indicating they may be the products of recent transposition events. Consistent with the recent transposition hypothesis, we found that 70% (125/178) of the elements (across all families) have identical intra-element LTRs. Using the synonymous substitution rate that has been calculated previously for Drosophila (.016 substitutions per site per million years) and the intra-element LTR divergence calculated here, the average age of the remaining 30% (53/178) of the elements was found to be 137,000 +/-89,000 yr. Collectively, these results indicate that many full-length LTR retrotransposons present in the D. melanogaster genome have transposed well after this species diverged from its closest relative Drosophila simulans, 2.3 +/-.3 million years ago.

Diversity, origin, and distribution of retrotransposons (gypsy and copia) in conifers

We examined the diversity, evolution, and genomic organization of retroelements in a wide range of gymnosperms. In total, 165 fragments of the reverse transcriptase (RT) gene domain were sequenced from PCR products using newly designed primers for gypsy-like retrotransposons and well-known primers for copia-like retrotransposons; representatives of long interspersed nuclear element (LINE) retroposons were also found. Gypsy and copia-like retroelements are a major component of the gymnosperm genome, and in situ hybridization showed that individual element families were widespread across the chromosomes, consistent with dispersion and amplification via an RNA intermediate. Most of the retroelement families were widely distributed among the gymnosperms, including species with wide taxonomic separation from the Northern and Southern Hemispheres. When the gymnosperm sequences were analyzed together with retroelements from other species, the monophyletic origin of plant copia, gypsy, and LINE groups was well supported, with an additional clade including badnaviral and other, probably virus-related, plant sequences as well as animal and fungal gypsy elements. Plant retroelements showed high diversity within the phylogenetic trees of both copia and gypsy RT domains, with, for example, retroelement sequences from Arabidopsis thaliana being present in many supported groupings. No primary branches divided major taxonomic clades such as angiosperms, monocotyledons, gymnosperms, or conifers or (based on smaller samples) ferns, Gnetales, or Sphenopsida (Equisetum), suggesting that much of the existing diversity was present early in plant evolution, or perhaps that horizontal transfer of sequences has occurred. Within the phylogenetic trees for both gypsy and copia, two clearly monophyletic gymnosperm/conifer clades were revealed, providing evidence against recent horizontal transfer. The results put the evolution of the large and relatively conserved genome structure of gymnosperms into the context of the diversity of other groups of plants.

Jule from the fish Xiphophorus is the first complete vertebrate Ty3/Gypsy retrotransposon from the Mag family

Jule is the second complete long-terminal-repeat (LTR) Ty3/Gypsy retrotransposon identified to date in vertebrates. Jule, first isolated from the poeciliid fish Xiphophorus maculatus, is 4.8 kb in length, is flanked by two 202-bp LTRs, and encodes Gag (structural core protein) and Pol (protease, reverse transcriptase, RNase H, and integrase, in that order) but no envelope. There are three to four copies of Jule per haploid genome in X. maculatus. Two of them are located in a subtelomeric region of the sex chromosomes, where they are associated with the Xmrk receptor tyrosine kinase genes, of which oncogenic versions are responsible for the formation of hereditary melanoma in Xiphophorus. One almost intact copy of Jule was found in the first intron of the X-chromosomal allele of the Xmrk proto-oncogene, and a second, more corrupted copy is present only 56 nt downstream of the polyadenylation signal of the Xmrk oncogene. Jule-related elements were detected by Southern blot hybridization with less than 10 copies per haploid genome in numerous other poeciliids, as well as in more divergent fishes, including the medakafish Oryzias latipes and the tilapia Oreochromis niloticus. Database searches also identified Jule-related sequences in the zebrafish Danio rerio and in both genome project pufferfishes, Fugu rubripes and Tetraodon nigroviridis. Phylogenetic analysis revealed that Jule is the first member of the Mag family of Ty3/Gypsy retrotransposons described to date in vertebrates. This family includes the silkworm Mag and sea urchin SURL retrotransposons, as well as sequences from the nematode Caenorhabditis elegans. Additional related elements were identified in the genomes of the malaria mosquito Anopheles gambiae and the nematode Ascaris lumbricoides. Phylogeny of Mag-related elements suggested that the Mag family of retrotransposons is polyphyletic and is constituted of several ancient lineages that diverged before their host genomes more than 600 MYA.