High frequency cDNA recombination of the saccharomyces retrotransposon Ty5: The LTR mediates formation of tandem elements

Tandem LTR-retrotransposon structures are common and highly polymorphic in plant genomes

BACKGROUND

LTR-retrotransposons (LTR-RT) are a major component of plant genomes and important drivers of genome evolution. Most LTR-RT copies in plant genomes are defective elements found as truncated copies, nested insertions or as part of more complex structures. The recent availability of highly contiguous plant genome assemblies based on long-read sequences now allows to perform detailed characterization of these complex structures and to evaluate their importance for plant genome evolution.

RESULTS

The detailed analysis of two rice loci containing complex LTR-RT structures showed that they consist of tandem arrays of LTR copies sharing internal LTRs. Our analyses suggests that these LTR-RT tandems are the result of a single insertion and not of the recombination of two independent LTR-RT elements. Our results also suggest that gypsy elements may be more prone to form these structures. We show that these structures are highly polymorphic in rice and therefore have the potential to generate genetic variability. We have developed a computational pipeline (IDENTAM) that scans genome sequences and identifies tandem LTR-RT candidates. Using this tool, we have detected 266 tandems in a pangenome built from the genomes of 76 accessions of cultivated and wild rice, showing that tandem LTR-RT structures are frequent and highly polymorphic in rice. Running IDENTAM in the Arabidopsis, almond and cotton genomes showed that LTR-RT tandems are frequent in plant genomes of different size, complexity and ploidy level. The complexity of differentiating intra-element variations at the nucleotide level among haplotypes is very high, and we found that graph-based pangenomic methodologies are appropriate to resolve these structures.

CONCLUSIONS

Our results show that LTR-RT elements can form tandem arrays. These structures are relatively abundant and highly polymorphic in rice and are widespread in the plant kingdom. Future studies will contribute to understanding how these structures originate and whether the variability that they generate has a functional impact.

Marsupial satellite DNA as faithful reflections of long terminal repeat (LTR) retroelement structure

Long terminal repeat (LTR) retroelements, including endogenous retroviruses, are one of the origins of satellite DNAs. However, the vast majority of satellite DNAs originating from LTR retroelements consist of parts of the element. In addition, they frequently contain sequences unrelated to that element. Here we report a novel marsupial satellite DNA (named walbRep) that contains, and consists solely of, the entire sequence of an LTR retroelement (the walb element). As is common with LTR retroelements, walb copies exhibit length variation. We focused on the abundance of copies of a specific length (2.7 kb) in the genome of the red-necked wallaby. Cloning and analyses of long genomic DNA fragments revealed a satellite DNA in which the LTR sequence (0.4 kb) and the sequence of the internal region of a nonautonomous walb copy (2.3 kb) were repeated alternately. The junctions between these two components exhibited the same end-to-end arrangements as those in the walb element. This satellite organization could be accounted for by a simple formation model that includes slippage during chromosome pairing followed by homologous recombination but does not invoke any other types of rearrangements. We discuss the possible reasons why satellite DNAs having such structures are rarely found in mammals.

Identification of an integrase-independent pathway of retrotransposition

Retroviruses and long terminal repeat retrotransposons rely on integrase (IN) to insert their complementary DNA (cDNA) into the genome of host cells. Nevertheless, in the absence of IN, retroelements can retain notable levels of insertion activity. We have characterized the IN-independent pathway of Tf1 and found that insertion sites had homology to the primers of reverse transcription, which form single-stranded DNAs at the termini of the cDNA. In the absence of IN activity, a similar bias was observed with HIV-1. Our studies showed that the Tf1 insertions result from single-strand annealing, a noncanonical form of homologous recombination mediated by Rad52. By expanding our analysis of insertions to include repeat sequences, we found most formed tandem elements by inserting at preexisting copies of a related transposable element. Unexpectedly, we found that wild-type Tf1 uses the IN-independent pathway as an alternative mode of insertion.

Rapid evolution at the Drosophila telomere: transposable element dynamics at an intrinsically unstable locus

Drosophila telomeres have been maintained by three families of active transposable elements (TEs), HeT-A, TAHRE, and TART, collectively referred to as HTTs, for tens of millions of years, which contrasts with an unusually high degree of HTT interspecific variation. While the impacts of conflict and domestication are often invoked to explain HTT variation, the telomeres are unstable structures such that neutral mutational processes and evolutionary tradeoffs may also drive HTT evolution. We leveraged population genomic data to analyze nearly 10,000 HTT insertions in 85  Drosophila melanogaster genomes and compared their variation to other more typical TE families. We observe that occasional large-scale copy number expansions of both HTTs and other TE families occur, highlighting that the HTTs are, like their feral cousins, typically repressed but primed to take over given the opportunity. However, large expansions of HTTs are not caused by the runaway activity of any particular HTT subfamilies or even associated with telomere-specific TE activity, as might be expected if HTTs are in strong genetic conflict with their hosts. Rather than conflict, we instead suggest that distinctive aspects of HTT copy number variation and sequence diversity largely reflect telomere instability, with HTT insertions being lost at much higher rates than other TEs elsewhere in the genome. We extend previous observations that telomere deletions occur at a high rate, and surprisingly discover that more than one-third do not appear to have been healed with an HTT insertion. We also report that some HTT families may be preferentially activated by the erosion of whole telomeres, implying the existence of HTT-specific host control mechanisms. We further suggest that the persistent telomere localization of HTTs may reflect a highly successful evolutionary strategy that trades away a stable insertion site in order to have reduced impact on the host genome. We propose that HTT evolution is driven by multiple processes, with niche specialization and telomere instability being previously underappreciated and likely predominant.

Double insertion of transposable elements provides a substrate for the evolution of satellite DNA

Eukaryotic genomes are replete with repeated sequences in the form of transposable elements (TEs) dispersed across the genome or as satellite arrays, large stretches of tandemly repeated sequences. Many satellites clearly originated as TEs, but it is unclear how mobile genetic parasites can transform into megabase-sized tandem arrays. Comprehensive population genomic sampling is needed to determine the frequency and generative mechanisms of tandem TEs, at all stages from their initial formation to their subsequent expansion and maintenance as satellites. The best available population resources, short-read DNA sequences, are often considered to be of limited utility for analyzing repetitive DNA due to the challenge of mapping individual repeats to unique genomic locations. Here we develop a new pipeline called ConTExt that demonstrates that paired-end Illumina data can be successfully leveraged to identify a wide range of structural variation within repetitive sequence, including tandem elements. By analyzing 85 genomes from five populations of Drosophila melanogaster, we discover that TEs commonly form tandem dimers. Our results further suggest that insertion site preference is the major mechanism by which dimers arise and that, consequently, dimers form rapidly during periods of active transposition. This abundance of TE dimers has the potential to provide source material for future expansion into satellite arrays, and we discover one such copy number expansion of the DNA transposon hobo to approximately 16 tandem copies in a single line. The very process that defines TEs—transposition—thus regularly generates sequences from which new satellites can arise.

Phosphorylation regulates integration of the yeast Ty5 retrotransposon into heterochromatin

The yeast Ty5 retrotransposon preferentially integrates into heterochromatin at the telomeres and silent mating loci. Target specificity is mediated by a small domain of Ty5 integrase (the targeting domain, TD), which interacts with the heterochromatin protein Sir4 and tethers the integration complex to target sites. Here we demonstrate that TD is phosphorylated and that phosphorylation is required for interaction with Sir4. The yeast cell, therefore, through posttranslational modification, controls Ty5's mutagenic potential: when TD is phosphorylated, insertions occur in gene-poor heterochromatin, thereby minimizing deleterious consequences of transposition; however, in the absence of phosphorylation, Ty5 integrates throughout the genome, frequently causing mutations. TD phosphorylation is reduced under stress conditions, specifically starvation for amino acids, nitrogen, or fermentable carbon. This suggests that Ty5 target specificity changes in response to nutrient availability and is consistent with McClintock's hypothesis that mobile elements restructure host genomes as an adaptive response to environmental challenge.

SREBP controls oxygen-dependent mobilization of retrotransposons in fission yeast

Retrotransposons are mobile genetic elements that proliferate through an RNA intermediate. Transposons do not encode transcription factors and thus rely on host factors for mRNA expression and survival. Despite information regarding conditions under which elements are upregulated, much remains to be learned about the regulatory mechanisms or factors controlling retrotransposon expression. Here, we report that low oxygen activates the fission yeast Tf2 family of retrotransposons. Sre1, the yeast ortholog of the mammalian membrane-bound transcription factor sterol regulatory element binding protein (SREBP), directly induces the expression and mobilization of Tf2 retrotransposons under low oxygen. Sre1 binds to DNA sequences in the Tf2 long terminal repeat that functions as an oxygen-dependent promoter. We find that Tf2 solo long terminal repeats throughout the genome direct oxygen-dependent expression of adjacent coding and noncoding sequences, providing a potential mechanism for the generation of oxygen-dependent gene expression.

Molecular characterization and genomic distribution of Isis: a new retrotransposon of Drosophila buzzatii

A new transposable element, Isis, is identified as a LTR retrotransposon in Drosophila buzzatii. DNA sequence analysis shows that Isis contains three long ORFs similar to gag, pol and env genes of retroviruses. The ORF1 exhibits sequence homology to matrix, capsid and nucleocapsid gag proteins and ORF2 encodes a putative protease (PR), a reverse transcriptase (RT), an Rnase H (RH) and an integrase (IN) region. The analysis of a putative env product, encoded by the env ORF3, shows a degenerated protein containing several stop codons. The molecular study of the putative proteins coded by this new element shows striking similarities to both Ulysses and Osvaldo elements, two LTR retrotransposons, present in D. virilis and D. buzzatii, respectively. Comparisons of the predicted Isis RT to several known retrotransposons show strong phylogenetic relationships to gypsy-like elements, particulary to Ulysses retrotransposon. Studies of Isis chromosomal distribution show a strong hybridization signal in centromeric and pericentromeric regions, and a scattered distribution along all chromosomal arms. The existence of insertional polymorphisms between different strains and high molecular weight bands by Southern blot suggests the existence of full-sized copies that have been active recently. The presence of euchromatic insertion sites coincident between Isis and Osvaldo could indicate preferential insertion sites of Osvaldo element into Isis sequence or vice versa. Moreover, the presence of Isis in different species of the buzzatii complex indicates the ancient origin of this element.

Happy together: the life and times of Ty retrotransposons and their hosts

The aim of this review is to describe the level of intimacy between Ty retrotransposons (Ty1-Ty5) and their host the yeast Saccharomyces cerevisiae. The effects of Ty location in the genome and of host proteins on the expression and mobility of Ty elements are highlighted. After a brief overview of Ty diversity and evolution, we describe the factors that dictate Ty target-site preference and the impact of targeting on Ty and adjacent gene expression. Studies on Ty3 and Ty5 have been especially informative in unraveling the role of host factors (Pol III machinery and silencing proteins, respectively) and integrase in controlling the specificity of integration. In contrast, not much is known regarding Ty1, Ty2 and Ty4, except that their insertion depends on the transcriptional competence of the adjacent Pol III gene and might be influenced by some chromatin components. This review also brings together recent findings on the regulation of Ty1 retrotransposition. A large number of host proteins (over 30) involved in a wide range of cellular processes controls either directly or indirectly Ty1 mobility, primarily at post-transcriptional steps. We focus on several genes for which more detailed analyses have permitted the elaboration of regulatory models. In addition, this review describes new data revealing that repression of Ty1 mobility also involves two forms of copy number control that act at both the trancriptional and post-transcriptional levels. Since S. cerevisiae lacks the conserved pathways for copy number control via transcriptional and post-transcriptional gene silencing found in other eukaryotes, Ty1 copy number control must be via another mechanism whose features are outlined. Ty1 response to stress also implicates activation at both transcriptional and postranscriptional steps of Ty1. Finally, we provide several insights in the role of Ty elements in chromosome evolution and yeast adaptation and discuss the factors that might limit Ty ectopic recombination.

Variability, Recombination, and Mosaic Evolution of the Barley BARE-1 Retrotransposon

BARE-1 is a highly abundant, copia-like, LTR (long terminal repeat) retrotransposon in the genus Hordeum. The LTRs provide the promoter, terminator, and polyadenylation signals necessary for the replicational life cycle of retrotransposons. We have examined the variability and evolution of BARE-1-like elements, focusing on the LTRs. Three groups were found, corresponding to each of the Hordeum genome types analyzed, which predate the divergence of these types. The most variable LTR regions are tandem repeats near the 3' end and the promoter. In barley (H. vulgare L.), two main classes of LTR promoters were defined, corresponding to BARE-1 and to a new class we call BARE-2. These can be considered as families within the group I BARE elements. Although less abundant in cultivated barley than is BARE-1, BARE-2 is transcriptionally active in leaves and calli. A sequenced BARE-2 has more than 99% similar LTRs and perfect terminal direct repeats (TDRs), indicating it is a recent insertion, but the coding region, especially gag, is disrupted by frameshifts and stop codons. BARE-2 appears to be a chimeric element resulting from retrotransposon recombination by strand switching during replication, with LTRs and 5'UTR more similar to BARE-1 and the rest more similar to Wis-2. We provide evidence as well for another form of recombination, where LTR-LTR recombination has generated tandem multimeric BARE-1 elements in which internal coding domains are interspersed with shared LTRs. The data indicate that recombination contributes to the complexity and plasticity of retroelement evolution in plant genomes.

Post-transcriptional Cosuppression of Ty1 Retrotransposition

To determine whether homology-dependent gene silencing or cosuppression mechanisms underlie copy number control (CNC) of Ty1 retrotransposition, we introduced an active Ty1 element into a naïve strain. Single Ty1 element retrotransposition was elevated in a Ty1-less background, but decreased dramatically when additional elements were present. Transcription from the suppressing Ty1 elements enhanced CNC but translation or reverse transcription was not required. Ty1 CNC occurred with a transcriptionally active Ty2 element, but not with Ty3 or Ty5 elements. CNC also occurred when the suppressing Ty1 elements were transcriptionally silenced, fused to the constitutive PGK1 promoter, or contained a minimal segment of mostly TYA1-gag sequence. Ty1 transcription of a multicopy element expressed from the GAL1 promoter abolished CNC, even when the suppressing element was defective for transposition. Although Ty1 RNA and TyA1-gag protein levels increased with the copy number of expressible elements, a given element's transcript level varied less than twofold regardless of whether the suppressing elements were transcriptionally active or repressed. Furthermore, a decrease in the synthesis of Ty1 cDNA is strongly associated with Ty1 CNC. Together our results suggest that Ty1 cosuppression can occur post-transcriptionally, either prior to or during reverse transcription.

Controlling integration specificity of a yeast retrotransposon

Retrotransposons and retroviruses integrate nonrandomly into eukaryotic genomes. For the yeast retrotransposon Ty5, integration preferentially occurs within domains of heterochromatin. Targeting to these locations is determined by interactions between an amino acid sequence motif at the C terminus of Ty5 integrase (IN) called the targeting domain, and the heterochromatin protein Sir4p. Here we show that new Ty5 integration hot spots are created when Sir4p is tethered to ectopic DNA sites. Targeting to sites of tethered Sir4p is abrogated by single amino acid substitutions in either IN or Sir4p that prevent their interaction. Ty5 target specificity can be altered by replacing the IN-targeting domain with other peptide motifs that interact with known protein partners. Integration occurs at high efficiency and in close proximity to DNA sites where the protein partners are tethered. These findings define a mechanism by which retrotransposons shape their host genomes and suggest ways in which retroviral integration can be controlled.

Ty5 gag mutations increase retrotransposition and suggest a role for hydrogen bonding in the function of the nucleocapsid zinc finger

The Ty5 retrotransposon of Saccharomyces paradoxus transposes in Saccharomyces cerevisiae at frequencies 1,000-fold lower than do the native Ty1 elements. The low transposition activity of Ty5 could be due to differences in cellular environments between these yeast species or to naturally occurring mutations in Ty5. By screening of a Ty5 mutant library, two single mutants (D252N and Y68C) were each found to increase transposition approximately sixfold. When combined, transposition increased 36-fold, implying that the two mutations act independently. Neither mutation affected Ty5 protein synthesis, processing, cDNA recombination, or target site choice. However, cDNA levels in both single mutants and the double mutant were significantly higher than in the wild type. The D252N mutation resides in the zinc finger of nucleocapsid and increases the potential for hydrogen bonding with nucleic acids. We generated other mutations that increase the hydrogen bonding potential (i.e., D252R and D252K) and found that they similarly increased transposition. This suggests that hydrogen bonding within the zinc finger motif is important for cDNA production and builds upon previous studies implicating basic amino acids flanking the zinc finger as important for zinc finger function. Although NCp zinc fingers differ from the zinc finger motifs of cellular enzymes, the requirement for efficient hydrogen bonding is likely universal.

Coupling of enhancer and insulator properties identified in two retrotransposons modulates their mutagenic impact on nearby genes

We recently reported a novel transposition system in which two retroelements from Drosophila melanogaster, ZAM and Idefix, are highly mobilized and preferentially insert within intergenic regions. Among the loci where new copies are detected, a hot spot for their insertion was identified at the white locus, where up to three elements occurred within a 3-kb fragment upstream of the transcriptional start site of white. We have used these insertions as molecular entry points to throw light on the mutagenic effect exerted by multiple insertions of retrotransposons within intergenic regions of a genome. Analysis of the molecular mechanisms by which ZAM and Idefix elements interfere with the regulation of the white gene has shown that ZAM bears cis-acting regulatory sequences able to enhance transcription of the white gene in the eyes of the flies. This activation may be counteracted by Idefix, which acts as an insulator able to isolate the white gene from the upstream ZAM enhancer. In addition to revealing a novel insulator sequence with its own specific features, our data clearly illustrate how retroelements can act as epigenetic factors able to interfere with the transcriptional regulation of their host.

Targeting of the yeast Ty5 retrotransposon to silent chromatin is mediated by interactions between integrase and Sir4p

The Ty5 retrotransposons of Saccharomyces cerevisiae integrate preferentially into regions of silent chromatin at the telomeres and silent mating loci (HMR and HML). We define a Ty5-encoded targeting domain that spans 6 amino acid residues near the C terminus of integrase (LXSSXP). The targeting domain establishes silent chromatin when it is tethered to a weakened HMR-E silencer, and it disrupts telomeric silencing when it is overexpressed. As determined by both yeast two-hybrid and in vitro binding assays, the targeting domain interacts with the C terminus of Sir4p, a structural component of silent chromatin. This interaction is abrogated by mutations in the targeting domain that disrupt integration into silent chromatin, suggesting that recognition of Sir4p by the targeting domain is the primary determinant in Ty5 target specificity.

Expression and processing of proteins encoded by the Saccharomyces retrotransposon Ty5

Retroelements (retrotransposons and retroviruses) have two genes in common: gag, which specifies structural proteins that form a virus or virus-like particle, and pol, which specifies catalytic proteins required for replication. For many retroelements, gag and pol are present on separate reading frames. Their expression is highly regulated, and the ratio of Gag to Pol is critical for retroelement replication. The Saccharomyces retrotransposon Ty5 contains a single open reading frame, and we characterized Gag and Pol expression by generating transpositionally active Ty5 elements with epitope tags at the N terminus or C terminus or within the integrase coding region. Immunoblot analysis identified two Gag species (Gag-p27 and Gag-p37), reverse transcriptase (Pol-p59), and integrase (Pol-p80), all of which are largely insoluble in the absence of urea or ionic detergent. These proteins result from proteolytic processing of a polyprotein, because elements with mutations in the presumed active site of Ty5 protease express a single tagged protein (Gag-Pol-p182). Protease mutants are also transpositionally inactive. In a time course experiment, we monitored protein expression, proteolytic processing, and transposition of a Ty5 element with identical epitope tags at its N and C termini. Both transposition and the abundance of Gag-p27 increased over time. In contrast, the levels of Gag-p37 and reverse transcriptase peaked after approximately 14 h of induction and then gradually decreased. This may be due to differences in stability of Gag-p27 relative to Gag-p37 and reverse transcriptase. The ratio of Ty5 Gag to Pol averaged 5:1 throughout the time course experiment, suggesting that differential protein stability regulates the amounts of these proteins.

The yeast retrotransposon Ty5 uses the anticodon stem-loop of the initiator methionine tRNA as a primer for reverse transcription

Retrotransposons and retroviruses replicate by reverse transcription of an mRNA intermediate. Most retroelements initiate reverse transcription from a host-encoded tRNA primer. DNA synthesis typically extends from the 3'-OH of the acceptor stem, which is complementary to sequences on the retroelement mRNA (the primer binding site, PBS). However, for some retrotransposons, including the yeast Ty5 elements, sequences in the anticodon stem-loop of the initiator methionine tRNA (IMT) are complementary to the PBS. We took advantage of the genetic tractability of the yeast system to investigate the mechanism of Ty5 priming. We found that transposition frequencies decreased at least 800-fold for mutations in the Ty5 PBS that disrupt complementarity with the IMT. Similarly, transposition was reduced at least 200-fold for IMT mutations in the anticodon stem-loop. Base pairing between the Ty5 PBS and IMT is essential for transposition, as compensatory changes that restored base pairing between the two mutant RNAs restored transposition significantly. An analysis of 12 imt mutants with base changes outside of the region of complementarity failed to identify other tRNA residues important for transposition. In addition, assays carried out with heterologous IMTs from Schizosaccharomyces pombe and Arabidopsis thaliana indicated that residues outside of the anticodon stem-loop have at most a fivefold effect on transposition. Our genetic system should make it possible to further define the components required for priming and to understand the mechanism by which Ty5's novel primer is generated.

cDNA of the yeast retrotransposon Ty5 preferentially recombines with substrates in silent chromatin

The yeast retrotransposon Ty5 preferentially integrates into regions of silent chromatin. Ty5 cDNA also recombines with homologous sequences, generating tandem elements or elements that have exchanged markers between cDNA and substrate. In this study, we demonstrate that Ty5 integration depends upon the conserved DD(35)E domain of integrase and cis-acting sequences at the end of the long terminal repeat (LTR) implicated in integrase binding. cDNA recombination requires Rad52p, which is responsible for homologous recombination. Interestingly, Ty5 cDNA recombines at least three times more frequently with substrates in silent chromatin than with a control substrate at an internal chromosomal locus. This preference depends upon the Ty5 targeting domain that is responsible for integration specificity, suggesting that localization of cDNA to silent chromatin results in the enhanced recombination. Recombination with a telomeric substrate occasionally generates highly reiterated Ty5 arrays, and mechanisms for tandem element formation were explored by using a plasmid-based recombination assay. Point mutations were introduced into plasmid targets, and recombination products were characterized to determine recombination initiation sites. Despite our previous observation of the importance of the LTR in forming tandem elements, recombination cannot simply be explained by crossover events between the LTRs of substrate and cDNA. We propose an alternative model based on single-strand annealing, where single-stranded cDNA initiates tandem element formation and the LTR is required for strand displacement to form a looped intermediate. Retrotransposons are increasingly found associated with chromosome ends, and amplification of Ty5 by both integration and recombination exemplifies how retroelements can contribute to telomere dynamics.

Schizosaccharomyces pombe retrotransposon Tf2 mobilizes primarily through homologous cDNA recombination

The Tf2 retrotransposon, found in the fission yeast Schizosaccharomyces pombe, is nearly identical to its sister element, Tf1, in its reverse transcriptase-RNase H and integrase domains but is very divergent in the gag domain, the protease, the 5' untranslated region, and the U3 domain of the long terminal repeats. It has now been demonstrated that a neo-marked copy of Tf2 overexpressed from a heterologous promoter can mobilize into the S. pombe genome and produce true transposition events. However, the Tf2-neo mobilization frequency is 10- to 20-fold lower than that of Tf1-neo, and 70% of the Tf2-neo events are homologous recombination events generated independently of a functional Tf2 integrase. Thus, the Tf2 element is primarily dependent on homologous recombination with preexisting copies of Tf2 for its propagation. Finally, production of Tf2-neo proteins and cDNA was also analyzed; surprisingly, Tf2 was found to produce its reverse transcriptase as a single species in which it is fused to protease, unlike all other retroviruses and retrotransposons.