Role of RNA polymerase III transcription factors in the selection of integration sites by the dictyostelium non-long terminal repeat retrotransposon TRE5-A

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.

A host factor supports retrotransposition of the TRE5-A population in Dictyostelium cells by suppressing an Argonaute protein

Background

In the compact and haploid genome of Dictyostelium discoideum control of transposon activity is of particular importance to maintain viability. The non-long terminal repeat retrotransposon TRE5-A amplifies continuously in D. discoideum cells even though it produces considerable amounts of minus-strand (antisense) RNA in the presence of an active RNA interference machinery. Removal of the host-encoded C-module-binding factor (CbfA) from D. discoideum cells resulted in a more than 90 % reduction of both plus- and minus-strand RNA of TRE5-A and a strong decrease of the retrotransposition activity of the cellular TRE5-A population. Transcriptome analysis revealed an approximately 230-fold overexpression of the gene coding for the Argonaute-like protein AgnC in a CbfA-depleted mutant.

Results

The D. discoideum genome contains orthologs of RNA-dependent RNA polymerases, Dicer-like proteins, and Argonaute proteins that are supposed to represent RNA interference pathways. We analyzed available mutants in these genes for altered expression of TRE5-A. We found that the retrotransposon was overexpressed in mutants lacking the Argonaute proteins AgnC and AgnE. Because the agnC gene is barely expressed in wild-type cells, probably due to repression by CbfA, we employed a new method of promoter-swapping to overexpress agnC in a CbfA-independent manner. In these strains we established an in vivo retrotransposition assay that determines the retrotransposition frequency of the cellular TRE5-A population. We observed that both the TRE5-A steady-state RNA level and retrotransposition rate dropped to less than 10 % of wild-type in the agnC overexpressor strains.

Conclusions

The data suggest that TRE5-A amplification is controlled by a distinct pathway of the Dictyostelium RNA interference machinery that does not require RNA-dependent RNA polymerases but involves AgnC. This control is at least partially overcome by the activity of CbfA, a factor derived from the retrotransposon’s host. This unusual regulation of mobile element activity most likely had a profound effect on genome evolution in D. discoideum.

Electronic supplementary material

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

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.

Experimental annotation of the human pathogen Candida albicans coding and noncoding transcribed regions using high-resolution tiling arrays

BACKGROUND

Compared to other model organisms and despite the clinical relevance of the pathogenic yeast Candida albicans, no comprehensive analysis has been done to provide experimental support of its in silico-based genome annotation.

RESULTS

We have undertaken a genome-wide experimental annotation to accurately uncover the transcriptional landscape of the pathogenic yeast C. albicans using strand-specific high-density tiling arrays. RNAs were purified from cells growing under conditions relevant to C. albicans pathogenicity, including biofilm, lab-grown yeast and serum-induced hyphae, as well as cells isolated from the mouse caecum. This work provides a genome-wide experimental validation for a large number of predicted ORFs for which transcription had not been detected by other approaches. Additionally, we identified more than 2,000 novel transcriptional segments, including new ORFs and exons, non-coding RNAs (ncRNAs) as well as convincing cases of antisense gene transcription. We also characterized the 5' and 3' UTRs of expressed ORFs, and established that genes with long 5' UTRs are significantly enriched in regulatory functions controlling filamentous growth. Furthermore, we found that genomic regions adjacent to telomeres harbor a cluster of expressed ncRNAs. To validate and confirm new ncRNA candidates, we adapted an iterative strategy combining both genome-wide occupancy of the different subunits of RNA polymerases I, II and III and expression data. This comprehensive approach allowed the identification of different families of ncRNAs.

CONCLUSIONS

In summary, we provide a comprehensive expression atlas that covers relevant C. albicans pathogenic developmental stages in addition to the discovery of new ORF and non-coding genetic elements.

A sequence similar to tRNA 3 Lys gene is embedded in HIV-1 U3-R and promotes minus-strand transfer

We identified a sequence embedded in the U3/R region of HIV-1 RNA that is highly complementary to human tRNA₃^Lys. The free energy of annealing to tRNA₃^Lys is significantly lower for this sequence and the primer-binding site than for other similar length viral sequences. The only interruption in complementarity is a 29-nucleotide segment inserted where a tRNA intron would be expected. The insert contains the TATA box for viral RNA transcription. The embedded sequence includes a nine-nucleotide segment previously reported to aid minus strand transfer by binding the primer tRNA₃^Lys. Reconstituting transfer in vitro, we show that including segments from the embedded sequence in the acceptor template, beyond the nine nucleotides, further increases transfer efficiency. We propose that a tRNA₃^Lys gene was incorporated during HIV-1 evolution and retained largely intact because of its roles in transcription and strand transfer.

TFIIIC binding sites function as both heterochromatin barriers and chromatin insulators in Saccharomyces cerevisiae

Chromosomal sites of RNA polymerase III (Pol III) transcription have been demonstrated to have "extratranscriptional" functions, as the assembled Pol III complex can act as chromatin boundaries or pause sites for replication forks, can alter nucleosome positioning or affect transcription of neighboring genes, and can play a role in sister chromatid cohesion. Several studies have demonstrated that assembled Pol III complexes block the propagation of heterochromatin-mediated gene repression. Here we show that in Saccharomyces cerevisiae tRNA genes (tDNAs) and even partially assembled Pol III complexes containing only the transcription factor TFIIIC can exhibit chromatin boundary functions both as heterochromatin barriers and as insulators to gene activation. Both the TRT2 tDNA and the ETC4 site which binds only the TFIIIC complex prevented an upstream activation sequence from activating the GAL promoters in our assay system, effectively acting as chromatin insulators. Additionally, when placed downstream from the heterochromatic HMR locus, ETC4 blocked the ectopic spread of Sir protein-mediated silencing, thus functioning as a barrier to repression. Finally, we show that TRT2 and the ETC6 site upstream of TFC6 in their natural contexts display potential insulator-like functions, and ETC6 may represent a novel case of a Pol III factor directly regulating a Pol II promoter. The results are discussed in the context of how the TFIIIC transcription factor complex may function to demarcate chromosomal domains in yeast and possibly in other eukaryotes.

Patterns of evolution in the unique tRNA gene arrays of the genus Entamoeba

Genome sequencing of the protistan parasite Entamoeba histolytica HM-1:IMSS revealed that almost all the tRNA genes are organized into tandem arrays that make up over 10% of the genome. The 25 distinct array units contain up to 5 tRNA genes each and some also encode the 5S RNA. Between adjacent genes in array units are complex short tandem repeats (STRs) resembling microsatellites. To investigate the origins and evolution of this unique gene organization, we have undertaken a genome survey to determine the array unit organization in 4 other species of Entamoeba-Entamoeba dispar, Entamoeba moshkovskii, Entamoeba terrapinae, and Entamoeba invadens-and have explored the STR structure in other isolates of E. histolytica. The genome surveys revealed that E. dispar has the same array unit organization as E. histolytica, including the presence and numerical variation of STRs between adjacent genes. However, the individual repeat sequences are completely different to those in E. histolytica. All other species of Entamoeba studied also have tandem arrays of clustered tRNA genes, but the gene composition of the array units often differs from that in E. histolytica/E. dispar. None of the other species' arrays exhibit the complex STRs between adjacent genes although simple tandem duplications are occasionally seen. The degree of similarity in organization reflects the phylogenetic relationships among the species studied. Within individual isolates of E. histolytica most copies of the array unit are uniform in sequence with only minor variation in the number and organization of the STRs. Between isolates, however, substantial differences in STR number and organization can exist although the individual repeat sequences tend to be conserved. The origin of this unique gene organization in the genus Entamoeba clearly predates the common ancestor of the species investigated to date and their function remains unclear.

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

Mobile genetic elements that reside in gene-dense genomes face the problem of avoiding devastating insertional mutagenesis of genes in their host cell genomes. To meet this challenge, some Saccharomyces cerevisiae long terminal repeat (LTR) retrotransposons have evolved targeted integration at safe sites in the immediate vicinity of tRNA genes. Integration of yeast Ty3 is mediated by interactions of retrotransposon protein with the tRNA gene-specific transcription factor IIIB (TFIIIB). In the genome of the social amoeba Dictyostelium discoideum, the non-LTR retrotransposon TRE5-A integrates approximately 48 bp upstream of tRNA genes, yet little is known about how the retrotransposon identifies integration sites. Here, we show direct protein interactions of the TRE5-A ORF1 protein with subunits of TFIIIB, suggesting that ORF1p is a component of the TRE5-A preintegration complex that determines integration sites. Our results demonstrate that evolution has put forth similar solutions to prevent damage of diverse, compact genomes by different classes of mobile elements.

Suffix-specific RNAi leads to silencing of F element in Drosophila melanogaster

Separate conserved copies of suffix, a short interspersed Drosophila retroelement (SINE), and also divergent copies in the 3′ untranslated regions of the three genes, have already been described. Suffix has also been identified on the 3′ end of the Drosophila non-LTR F element, where it forms the last conserved domain of the reverse transcriptase (RT). In our current study, we show that the separate copies of suffix are far more actively transcribed than their counterparts on the F element. Transcripts from both strands of suffix are present in RNA preparations during all stages of Drosophila development, providing the potential for the formation of double-stranded RNA and the initiation of RNA interference (RNAi). Using in situ RNA hybridization analysis, we have detected the expression of both sense and antisense suffix transcripts in germinal cells. These sense and antisense transcripts are colocalized in the primary spermatocytes and in the cytoplasm of the nurse cells, suggesting that they form double-stranded RNA. We performed further analyses of suffix-specific small RNAs using northern blotting and SI nuclease protection assays. Among the total RNA preparations isolated from embryos, larvae, pupae and flies, suffix-specific small interfering RNAs (siRNAs) were detected only in pupae. In wild type ovaries, both the siRNAs and longer suffix-specific Piwi-interacting RNAs (piRNAs) were observed, whereas in ovaries of the Dicer-2 mutant, only piRNAs were detected. We further found by 3′ RACE that in pupae and ovaries, F element transcripts lacking the suffix sequence are also present. Our data provide direct evidence that suffix-specific RNAi leads to the silencing of the relative LINE (long interspersed element), F element, and suggests that SINE-specific RNA interference could potentially downregulate a set of genes possessing SINE stretches in their 5′ or 3′ non-coding regions. These data also suggest that double stranded RNAs possessing suffix are processed by both RNAi and an additional silencing mechanism.