Evolutionary dynamics of the Ty3/gypsy LTR retrotransposons in the genome of Anopheles gambiae

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis.     

About the origin of retroviruses and the co-evolution of the gypsy retrovirus with the Drosophila flamenco host gene

The gypsy element of Drosophila melanogaster is the first retrovirus identified so far in invertebrates. According to phylogenetic data, gypsy belongs to the same group as the Ty3 class of LTR-retrotransposons, which suggests that retroviruses evolved from this kind of retroelements before the radiation of vertebrates. There are other invertebrate retroelements that are also likely to be endogenous retroviruses because they share with gypsy some structural and functional retroviral-like characteristics. Gypsy is controlled by a Drosophila gene called flamenco, the restrictive alleles of which maintain the retrovirus in a repressed state. In permissive strains, functional gypsy elements transpose at high frequency and produce infective particles. Defective gypsy proviruses located in pericentromeric heterochromatin of all strains seem to be very old components of the genome of Drosophila melanogaster, which indicates that gypsy invaded this species, or an ancestor, a long time ago. At that time, Drosophila melanogaster presumably contained permissive alleles of the flamenco gene. One can imagine that the species survived to the increase of genetic load caused by the retroviral invasion because restrictive alleles of flamenco were selected. The characterization of a retrovirus in Drosophila, one of the most advanced model organisms for molecular genetics, provides us with an exceptional clue to study how a species can resist a retroviral invasion. This revised version was published online in August 2006 with corrections to the Cover Date.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

Pairing between gypsy insulators facilitates the enhancer action in trans throughout the Drosophila genome

The Suppressor of the Hairy wing [Su(Hw)] binding region within the gypsy retrotransposon is the best known chromatin insulator in Drosophila melanogaster. According to previous data, two copies of the gypsy insulator inserted between an enhancer and a promoter neutralize each other's actions, which is indicative of an interaction between the protein complexes bound to the insulators. We have investigated the role of pairing between the gypsy insulators located on homologous chromosomes in trans interaction between yellow enhancers and a promoter. It has been shown that trans activation of the yellow promoter strongly depends on the site of the transposon insertion, which is evidence for a role of surrounding chromatin in homologous pairing. The presence of the gypsy insulators in both homologous chromosomes even at a distance of 9 kb downstream from the promoter dramatically improves the trans activation of yellow. Moreover, the gypsy insulators have proved to stabilize trans activation between distantly located enhancers and a promoter. These data suggest that gypsy insulator pairing is involved in communication between loci in the Drosophila genome.     

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.     

Mobilization of the gypsy and copia retrotransposons in Drosophila melanogaster induces reversion of the ovo dominant female-sterile mutations: molecular analysis of revertant alleles

The ovo locus is required for the maintenance of the female germ line in Drosophila melanogaster. In the absence of an ovo⁺ gene, males are completely normal but females have no germ-line stem cells. Three dominant mutations at the ovo locus, called ovo^D, were observed to revert towards recessive alleles at high frequency when ovo^D males were crossed to females of the strain y v f mal. We have found that this strain contains an inordinately high number of gypsy transposable elements, and crossing it with the ovo^D strains results in the mobilization of both gypsy and copia, with high-frequency insertions into the ovo locus: of 16 revertants examined 12 have gypsy and four have copia inserted at 4E, the ovo cytological site. Using gypsy DNA as a tag we have cloned 32 kb of wild-type DNA sequences surrounding a gypsy insertion and characterized molecular rearrangements in several independent revertants: in 10 of them gypsy appears to be inserted into the same site. The orientation of gypsy is strictly correlated with whether the neighbouring lozenge-like mutation appears in the revertants. A distal limit of the ovo locus was molecularly determined from the breakpoint of a deletion affecting closely flanking regions.     

Analysis of genes associated with retrotransposons in the rice genome

Retrotransposons comprise a significant fraction of the rice genome. Despite their prevalence, the effects of retrotransposon insertions are not well understood, especially with regard to how they affect the expression of genes. In this study, we identified one-sixth of rice genes as being associated with retrotransposons, with insertions either in the gene itself or within its putative promoter region. Among genes with insertions in the promoter region, the likelihood of the gene being expressed was shown to be directly proportional to the distance of the retrotransposon from the translation start site. In addition, retrotransposon insertions in the transcribed region of the gene were found to be positively correlated with the presence of alternative splicing forms. Furthermore, preferential association of retrotransposon insertions with genes in several functional classes was identified. Some of the retrotransposons that are part of full-length cDNA (fl-cDNA) contribute splice sites and give rise to novel exons. Several interesting trends concerning the effects of retrotransposon insertions on gene expression were identified. Taken together, our data suggests that retrotransposon association with genes have a role in gene regulation. The data presented in this study provides a foundation for experimental studies to determine the role of retrotransposons in gene regulation.     

HB mobile element in the Drosophila melanogaster genome: structural and functional analyses

The structure was analyzed for 60 annotated copies of the mobile genetic element (MGE) HB from the Drosophila melanogaster genome. The genomic distribution of HB copies was studied, and preferential insertion sites (hot spots) were identified, which presumably amount to several kilobases. Structural analysis of the open reading frame (ORF) and terminal repeats of HB was performed. All 26 HB copies retaining the ORF sequence have a stop codon in the same position. Consequently, the HB ORF proved indeed to code for an enzyme of 148 amino acid residues, relatively small for Tc1-family transposases. The ORF consensus sequence was established. HB{}1185 was identified as the only HB copy potentially coding for a functional protein. All 37 repeat-containing HB copies were analyzed. Of these, only four had functional terminal sequences, lacking, however, a functional transposase gene. A new 7762-bp copy of MGE roo was found in the D. melanogaster genome; the copy was earlier unavailable from databases and represents an insert in the HB{}1605 sequence.