Protein interactions involved in tRNA gene-specific integration of Dictyostelium discoideum non-long terminal repeat retrotransposon TRE5-A

Nanopore sequencing for mapping of retrotransposon integration sites in the Dictyostelium discoideum genome

The unicellular eukaryote Dictyostelium discoideum has a gene-dense haploid genome. This configuration presents mobile elements with the particular challenge of replicating without causing excessive damage to the host through insertional mutagenesis or recombination between repetitive sequences. D. discoideum harbors an active population of the retrotransposon TRE5-A that integrates in a narrow window of ~50 bp upstream of tRNA genes. We assume that this integration preference was developed to avoid the disruption of protein-coding genes. Therefore, we recently mapped new integrations of a genetically tagged TRE5-A element at tRNA genes using PCR-based enrichment of integration junctions. However, the PCR-based enrichment produced several artificial DNA fusions that prevented the mapping of integration sites in unknown places of the genome. Here, we reanalyzed the previous experiment using nanopore sequencing. We summarize the advantages and limitations of direct genome resequencing for the mapping of mobile element integrations.

Retrotransposon targeting to RNA polymerase III-transcribed genes

Retrotransposons are genetic elements that are similar in structure and life cycle to retroviruses by replicating via an RNA intermediate and inserting into a host genome. The Saccharomyces cerevisiae (S. cerevisiae) Ty1-5 elements are long terminal repeat (LTR) retrotransposons that are members of the Ty1-copia (Pseudoviridae) or Ty3-gypsy (Metaviridae) families. Four of the five S. cerevisiae Ty elements are inserted into the genome upstream of RNA Polymerase (Pol) III-transcribed genes such as transfer RNA (tRNA) genes. This particular genomic locus provides a safe environment for Ty element insertion without disruption of the host genome and is a targeting strategy used by retrotransposons that insert into compact genomes of hosts such as S. cerevisiae and the social amoeba Dictyostelium. The mechanism by which Ty1 targeting is achieved has been recently solved due to the discovery of an interaction between Ty1 Integrase (IN) and RNA Pol III subunits. We describe the methods used to identify the Ty1-IN interaction with Pol III and the Ty1 targeting consequences if the interaction is perturbed. The details of Ty1 targeting are just beginning to emerge and many unexplored areas remain including consideration of the 3-dimensional shape of genome. We present a variety of other retrotransposon families that insert adjacent to Pol III-transcribed genes and the mechanism by which the host machinery has been hijacked to accomplish this targeting strategy. Finally, we discuss why retrotransposons selected Pol III-transcribed genes as a target during evolution and how retrotransposons have shaped genome architecture.

Retrotransposon Domestication and Control in Dictyostelium discoideum

Transposable elements, identified in all eukaryotes, are mobile genetic units that can change their genomic position. Transposons usually employ an excision and reintegration mechanism, by which they change position, but not copy number. In contrast, retrotransposons amplify via RNA intermediates, increasing their genomic copy number. Hence, they represent a particular threat to the structural and informational integrity of the invaded genome. The social amoeba Dictyostelium discoideum, model organism of the evolutionary Amoebozoa supergroup, features a haploid, gene-dense genome that offers limited space for damage-free transposition. Several of its contemporary retrotransposons display intrinsic integration preferences, for example by inserting next to transfer RNA genes or other retroelements. Likely, any retrotransposons that invaded the genome of the amoeba in a non-directed manner were lost during evolution, as this would result in decreased fitness of the organism. Thus, the positional preference of the Dictyostelium retroelements might represent a domestication of the selfish elements. Likewise, the reduced danger of such domesticated transposable elements led to their accumulation, and they represent about 10% of the current genome of D. discoideum. To prevent the uncontrolled spreading of retrotransposons, the amoeba employs control mechanisms including RNA interference and heterochromatization. Here, we review TRE5-A, DIRS-1 and Skipper-1, as representatives of the three retrotransposon classes in D. discoideum, which make up 5.7% of the Dictyostelium genome. We compile open questions with respect to their mobility and cellular regulation, and suggest strategies, how these questions might be addressed experimentally.

TRE5-A retrotransposition profiling reveals putative RNA polymerase III transcription complex binding sites on the Dictyostelium extrachromosomal rDNA element

The amoeba Dictyostelium discoideum has a haploid genome in which two thirds of the DNA encodes proteins. Consequently, the space available for selfish mobile elements to expand without excess damage to the host genome is limited. The non-long terminal repeat retrotransposon TRE5-A maintains an active population in the D. discoideum genome and apparently adapted to this gene-dense environment by targeting positions ~47 bp upstream of tRNA genes that are devoid of protein-coding regions. Because only ~24% of tRNA genes are associated with a TRE5-A element in the reference genome, we evaluated whether TRE5-A retrotransposition is limited to this subset of tRNA genes. We determined that a tagged TRE5-A element (TRE5-A^bsr) integrated at 384 of 405 tRNA genes, suggesting that expansion of the current natural TRE5-A population is not limited by the availability of targets. We further observed that TRE5-A^bsr targets the ribosomal 5S gene on the multicopy extrachromosomal DNA element that carries the ribosomal RNA genes, indicating that TRE5-A integration may extend to the entire RNA polymerase III (Pol III) transcriptome. We determined that both natural TRE5-A and cloned TRE5-A^bsr retrotranspose to locations on the extrachromosomal rDNA element that contain tRNA gene-typical A/B box promoter motifs without displaying any other tRNA gene context. Based on previous data suggesting that TRE5-A targets tRNA genes by locating Pol III transcription complexes, we propose that A/B box loci reflect Pol III transcription complex assembly sites that possess a function in the biology of the extrachromosomal rDNA element.

Integration site selection by retroviruses and transposable elements in eukaryotes

Transposable elements and retroviruses are found in most genomes, can be pathogenic and are widely used as gene-delivery and functional genomics tools. Exploring whether these genetic elements target specific genomic sites for integration and how this preference is achieved is crucial to our understanding of genome evolution, somatic genome plasticity in cancer and ageing, host-parasite interactions and genome engineering applications. High-throughput profiling of integration sites by next-generation sequencing, combined with large-scale genomic data mining and cellular or biochemical approaches, has revealed that the insertions are usually non-random. The DNA sequence, chromatin and nuclear context, and cellular proteins cooperate in guiding integration in eukaryotic genomes, leading to a remarkable diversity of insertion site distribution and evolutionary strategies.

Convergent evolution of tRNA gene targeting preferences in compact genomes

Background

In gene-dense genomes, mobile elements are confronted with highly selective pressure to amplify without causing excessive damage to the host. The targeting of tRNA genes as potentially safe integration sites has been developed by retrotransposons in various organisms such as the social amoeba Dictyostelium discoideum and the yeast Saccharomyces cerevisiae. In D. discoideum, tRNA gene-targeting retrotransposons have expanded to approximately 3 % of the genome. Recently obtained genome sequences of species representing the evolutionary history of social amoebae enabled us to determine whether the targeting of tRNA genes is a generally successful strategy for mobile elements to colonize compact genomes.

Results

During the evolution of dictyostelids, different retrotransposon types independently developed the targeting of tRNA genes at least six times. DGLT-A elements are long terminal repeat (LTR) retrotransposons that display integration preferences ~15 bp upstream of tRNA gene-coding regions reminiscent of the yeast Ty3 element. Skipper elements are chromoviruses that have developed two subgroups: one has canonical chromo domains that may favor integration in centromeric regions, whereas the other has diverged chromo domains and is found ~100 bp downstream of tRNA genes. The integration of D. discoideum non-LTR retrotransposons ~50 bp upstream (TRE5 elements) and ~100 bp downstream (TRE3 elements) of tRNA genes, respectively, likely emerged at the root of dictyostelid evolution. We identified two novel non-LTR retrotransposons unrelated to TREs: one with a TRE5-like integration behavior and the other with preference ~4 bp upstream of tRNA genes.

Conclusions

Dictyostelid retrotransposons demonstrate convergent evolution of tRNA gene targeting as a probable means to colonize the compact genomes of their hosts without being excessively mutagenic. However, high copy numbers of tRNA gene-associated retrotransposons, such as those observed in D. discoideum, are an exception, suggesting that the targeting of tRNA genes does not necessarily favor the amplification of position-specific integrating elements to high copy numbers under the repressive conditions that prevail in most host cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0073-9) contains supplementary material, which is available to authorized users.

Retrotransposons. An RNA polymerase III subunit determines sites of retrotransposon integration

Mobile genetic elements are ubiquitous. Their integration site influences genome stability and gene expression. The Ty1 retrotransposon of the yeast Saccharomyces cerevisiae integrates upstream of RNA polymerase III (Pol III)-transcribed genes, yet the primary determinant of target specificity has remained elusive. Here we describe an interaction between Ty1 integrase and the AC40 subunit of Pol III and demonstrate that AC40 is the predominant determinant targeting Ty1 integration upstream of Pol III-transcribed genes. Lack of an integrase-AC40 interaction dramatically alters target site choice, leading to a redistribution of Ty1 insertions in the genome, mainly to chromosome ends. The mechanism of target specificity allows Ty1 to proliferate and yet minimizes genetic damage to its host.

A superfamily of DNA transposons targeting multicopy small RNA genes

Target-specific integration of transposable elements for multicopy genes, such as ribosomal RNA and small nuclear RNA (snRNA) genes, is of great interest because of the relatively harmless nature, stable inheritance and possible application for targeted gene delivery of target-specific transposable elements. To date, such strict target specificity has been observed only among non-LTR retrotransposons. We here report a new superfamily of sequence-specific DNA transposons, designated Dada. Dada encodes a DDE-type transposase that shows a distant similarity to transposases encoded by eukaryotic MuDR, hAT, P and Kolobok transposons, as well as the prokaryotic IS256 insertion element. Dada generates 6-7 bp target site duplications upon insertion. One family of Dada DNA transposons targets a specific site inside the U6 snRNA genes and are found in various fish species, water flea, oyster and polycheate worm. Other target sequences of the Dada transposons are U1 snRNA genes and different tRNA genes. The targets are well conserved in multicopy genes, indicating that copy number and sequence conservation are the primary constraints on the target choice of Dada transposons. Dada also opens a new frontier for target-specific gene delivery application.

In vitro targeting of strand transfer by the Ty3 retroelement integrase

The Saccharomyces cerevisiae long terminal repeat retrotransposon Ty3 integrates within one or two nucleotides of the transcription initiation sites of genes transcribed by RNA polymerase III. In this study the minimal components required to re-constitute position-specific strand transfer by Ty3 integrase are defined. Ty3 integrase targeted by a synthetic fusion of RNA polymerase III transcription factor IIIB subunits, Brf1 and TBP, mediated position-specific strand transfer of duplex oligonucleotides representing the ends of the Ty3 cDNA. These results further delimit the TFIIIB domains targeted by the Ty3 element and show that IN is the Ty3 component sufficient in vitro to target integration. These results underscore the commonality of protein interactions that mediate transcription and retrotransposon targeting. Surprisingly, in the presence of MnCl(2), strand transfer was TFIIIB-independent and targeted sequences resembling the Ty3 terminal inverted repeat.

Dynamic interactions between transposable elements and their hosts

Transposable elements (TEs) have a unique ability to mobilize to new genomic locations, and the major advance of second-generation DNA sequencing has provided insights into the dynamic relationship between TEs and their hosts. It now is clear that TEs have adopted diverse strategies - such as specific integration sites or patterns of activity - to thrive in host environments that are replete with mechanisms, such as small RNAs or epigenetic marks, that combat TE amplification. Emerging evidence suggests that TE mobilization might sometimes benefit host genomes by enhancing genetic diversity, although TEs are also implicated in diseases such as cancer. Here, we discuss recent findings about how, where and when TEs insert in diverse organisms.

Genetically tagged TRE5-A retrotransposons reveal high amplification rates and authentic target site preference in the Dictyostelium discoideum genome

Retrotransposons contribute significantly to the evolution of eukaryotic genomes. They replicate by producing DNA copies of their own RNA, which are integrated at new locations in the host cell genome. In the gene-dense genome of the social amoeba Dictyostelium discoideum, retrotransposon TRE5-A avoids insertional mutagenesis by targeting the transcription factor (TF) IIIC/IIIB complex and integrating ∼ 50 bp upstream of tRNA genes. We generated synthetic TRE5-A retrotransposons (TRE5-A(bsr)) that were tagged with a selection marker that conferred resistance to blasticidin after a complete retrotransposition cycle. We found that the TRE5-A(bsr) elements were efficiently mobilized in trans by proteins expressed from the endogenous TRE5-A population found in D. discoideum cells. ORF1 protein translated from TRE5-A(bsr) elements significantly enhanced retrotransposition. We observed that the 5' untranslated region of TRE5-A could be replaced by an unrelated promoter, whereas the 3' untranslated region of TRE5-A was essential for retrotransposition. A predicted secondary structure in the RNA of the 3' untranslated region of TRE5-A may be involved in the retrotransposition process. The TRE5-A(bsr) elements were capable of identifying authentic integration targets in vivo, including formerly unnoticed, putative binding sites for TFIIIC on the extrachromosomal DNA element that carries the ribosomal RNA genes.

Dictyostelium transfer RNA gene-targeting retrotransposons: Studying mobile element-host interactions in a compact genome

The model species of social amoebae, Dictyostelium discoideum, has a compact genome consisting of about two thirds protein-coding regions, with intergenic regions that are rarely larger than 1,000 bp. We hypothesize that the haploid state of D. discoideum cells provides defense against the amplification of mobile elements whose transposition activities would otherwise lead to the accumulation of heterozygous, potentially lethal mutations in diploid populations. We further speculate that complex transposon clusters found on D. discoideum chromosomes do not a priori result from integration preferences of these transposons, but that the clusters instead result from negative selection against cells harboring insertional mutations in genes. D. discoideum cells contain a fraction of retrotransposons that are found in the close vicinity of tRNA genes. Growing evidence suggests that these retrotransposons use active recognition mechanisms to determine suitable integration sites. However, the question remains whether these retrotransposons also cause insertional mutagenesis of genes, resulting in their enrichment at tRNA genes, which are relatively safe sites in euchromatic regions. Recently developed in vivo retrotransposition assays will allow a detailed, genome-wide analysis of de novo integration events in the D. discoideum genome.

The C-module-binding factor supports amplification of TRE5-A retrotransposons in the Dictyostelium discoideum genome

Retrotransposable elements are molecular parasites that have invaded the genomes of virtually all organisms. Although retrotransposons encode essential proteins to mediate their amplification, they also require assistance by host cell-encoded machineries that perform functions such as DNA transcription and repair. The retrotransposon TRE5-A of the social amoeba Dictyostelium discoideum generates a notable amount of both sense and antisense RNAs, which are generated from element-internal promoters, located in the A module and the C module, respectively. We observed that TRE5-A retrotransposons depend on the C-module-binding factor (CbfA) to maintain high steady-state levels of TRE5-A transcripts and that CbfA supports the retrotransposition activity of TRE5-A elements. The carboxy-terminal domain of CbfA was found to be required and sufficient to mediate the accumulation of TRE5-A transcripts, but it did not support productive retrotransposition of TRE5-A. This result suggests different roles for CbfA protein domains in the regulation of TRE5-A retrotransposition frequency in D. discoideum cells. Although CbfA binds to the C module in vitro, the factor regulates neither C-module nor A-module promoter activity in vivo. We speculate that CbfA supports the amplification of TRE5-A retrotransposons by suppressing the expression of an as yet unidentified component of the cellular posttranscriptional gene silencing machinery.

Genomic organization of eukaryotic tRNAs

BACKGROUND

Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities.

RESULTS

In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size.

CONCLUSIONS

The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.

Targeted gene insertion for molecular medicine

Genomic insertion of a functional gene together with suitable transcriptional regulatory elements is often required for long-term therapeutical benefit in gene therapy for several genetic diseases. A variety of integrating vectors for gene delivery exist. Some of them exhibit random genomic integration, whereas others have integration preferences based on attributes of the targeted site, such as primary DNA sequence and physical structure of the DNA, or through tethering to certain DNA sequences by host-encoded cellular factors. Uncontrolled genomic insertion bears the risk of the transgene being silenced due to chromosomal position effects, and can lead to genotoxic effects due to mutagenesis of cellular genes. None of the vector systems currently used in either preclinical experiments or clinical trials displays sufficient preferences for target DNA sequences that would ensure appropriate and reliable expression of the transgene and simultaneously prevent hazardous side effects. We review in this paper the advantages and disadvantages of both viral and non-viral gene delivery technologies, discuss mechanisms of target site selection of integrating genetic elements (viruses and transposons), and suggest distinct molecular strategies for targeted gene delivery.