Access to DNA establishes a secondary target site bias for the yeast retrotransposon Ty5

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

The local integration preference of the Tf1 retrotransposon in Schizosaccharomyces pombe

Transposons are mobile DNAs that can move to different locations in host genomes. The integration site selection of transposons is critical for both themselves and host cells. Studies on the integration of retrotransposons and retroviruses have focused more on the global preference than on the local preference. The local preferences of retrotransposons are usually weak and of large diversity. Here, we analyzed hundreds of thousands of independent integration events of the Tf1 retrotransposon in Schizosaccharomyces pombe. The consensus sequence at the Tf1 integration sites shows a palindromic pattern, which can be divided into four sections, each of them contains one or more CGnTA units with a period of 10 base pairs, indicating interaction with subunits of the integrase oligomer in the pre-integration complex. Moreover, the analysis on the nucleosome occupancy flanking Tf1 target sites shows that Tf1 integration favors regions with one entire nucleosome depletion.

Light and shadow on the mechanisms of integration site selection in yeast Ty retrotransposon families

Transposable elements are ubiquitous in genomes. Their successful expansion depends in part on their sites of integration in their host genome. In Saccharomyces cerevisiae, evolution has selected various strategies to target the five Ty LTR-retrotransposon families into gene-poor regions in a genome, where coding sequences occupy 70% of the DNA. The integration of Ty1/Ty2/Ty4 and Ty3 occurs upstream and at the transcription start site of the genes transcribed by RNA polymerase III, respectively. Ty5 has completely different integration site preferences, targeting heterochromatin regions. Here, we review the history that led to the identification of the cellular tethering factors that play a major role in anchoring Ty retrotransposons to their preferred sites. We also question the involvement of additional factors in the fine-tuning of the integration site selection, with several studies converging towards an importance of the structure and organization of the chromatin.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

A small targeting domain in Ty1 integrase is sufficient to direct retrotransposon integration upstream of tRNA genes

Integration of transposable elements into the genome is mutagenic. Mechanisms targeting integrations into relatively safe locations, hence minimizing deleterious consequences for cell fitness, have emerged during evolution. In budding yeast, integration of the Ty1 LTR retrotransposon upstream of RNA polymerase III (Pol III)-transcribed genes requires interaction between Ty1 integrase (IN1) and AC40, a subunit common to Pol I and Pol III. Here, we identify the Ty1 targeting domain of IN1 that ensures (i) IN1 binding to Pol I and Pol III through AC40, (ii) IN1 genome-wide recruitment to Pol I- and Pol III-transcribed genes, and (iii) Ty1 integration only at Pol III-transcribed genes, while IN1 recruitment by AC40 is insufficient to target Ty1 integration into Pol I-transcribed genes. Swapping the targeting domains between Ty5 and Ty1 integrases causes Ty5 integration at Pol III-transcribed genes, indicating that the targeting domain of IN1 alone confers Ty1 integration site specificity.

Diverse transposable element landscapes in pathogenic and nonpathogenic yeast models: the value of a comparative perspective

Genomics and other large-scale analyses have drawn increasing attention to the potential impacts of transposable elements (TEs) on their host genomes. However, it remains challenging to transition from identifying potential roles to clearly demonstrating the level of impact TEs have on genome evolution and possible functions that they contribute to their host organisms. I summarize TE content and distribution in four well-characterized yeast model systems in this review: the pathogens Candida albicans and Cryptococcus neoformans, and the nonpathogenic species Saccharomyces cerevisiae and Schizosaccharomyces pombe. I compare and contrast their TE landscapes to their lifecycles, genomic features, as well as the presence and nature of RNA interference pathways in each species to highlight the valuable diversity represented by these models for functional studies of TEs. I then review the regulation and impacts of the Ty1 and Ty3 retrotransposons from Saccharomyces cerevisiae and Tf1 and Tf2 retrotransposons from Schizosaccharomyces pombe to emphasize parallels and distinctions between these well-studied elements. I propose that further characterization of TEs in the pathogenic yeasts would enable this set of four yeast species to become an excellent set of models for comparative functional studies to address outstanding questions about TE-host relationships.

Genome-Wide Mapping of Yeast Retrotransposon Integration Target Sites

Transposable elements (TEs) are present in virtually all organisms. TE integration into genomes contributes to their structure and evolution, but can also be harmful in some cases. Deciphering where and how TE integration is targeted is fundamental to understand their intricate relationship with their host and explore the outcome of TE mobility on genome evolution and cell fitness. In general, TEs display integration site preference, which differs between elements. High-throughput mapping of de novo insertions by deep sequencing has recently allowed identifying genome-wide integration preferences of several TEs. These studies have provided invaluable clues to address the molecular determinants of integration site preference. Here, we provide a step-by-step methodology to generate massive de novo insertion events and prepare a library of genomic DNA for next-generation sequencing. We also describe a primary bioinformatic procedure to map these insertions in the genome. The whole procedure comes from our recent work on the integration of Ty1 in Saccharomyces cerevisiae, but could be easily adapted to the study of other TEs.

McClintock: An Integrated Pipeline for Detecting Transposable Element Insertions in Whole-Genome Shotgun Sequencing Data

Transposable element (TE) insertions are among the most challenging types of variants to detect in genomic data because of their repetitive nature and complex mechanisms of replication . Nevertheless, the recent availability of large resequencing data sets has spurred the development of many new methods to detect TE insertions in whole-genome shotgun sequences. Here we report an integrated bioinformatics pipeline for the detection of TE insertions in whole-genome shotgun data, called McClintock (https://github.com/bergmanlab/mcclintock), which automatically runs and standardizes output for multiple TE detection methods. We demonstrate the utility of McClintock by evaluating six TE detection methods using simulated and real genome data from the model microbial eukaryote, Saccharomyces cerevisiae We find substantial variation among McClintock component methods in their ability to detect nonreference TEs in the yeast genome, but show that nonreference TEs at nearly all biologically realistic locations can be detected in simulated data by combining multiple methods that use split-read and read-pair evidence. In general, our results reveal that split-read methods detect fewer nonreference TE insertions than read-pair methods, but generally have much higher positional accuracy. Analysis of a large sample of real yeast genomes reveals that most McClintock component methods can recover known aspects of TE biology in yeast such as the transpositional activity status of families, target preferences, and target site duplication structure, albeit with varying levels of accuracy. Our work provides a general framework for integrating and analyzing results from multiple TE detection methods, as well as useful guidance for researchers studying TEs in yeast resequencing data.

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.

Transcription factor regulation and chromosome dynamics during pseudohyphal growth

A multiplexed analysis of the transcriptional regulation of yeast pseudohyphal growth recorded the binding of 28 different transcription factors with barcoded transposons. A core set of target genes is identified, and a process of DNA looping at the FLO11 locus that provides transcriptional memory for expression of the gene is described.

Pseudohyphal growth is a developmental pathway seen in some strains of yeast in which cells form multicellular filaments in response to environmental stresses. We used multiplexed transposon “Calling Cards” to record the genome-wide binding patterns of 28 transcription factors (TFs) in nitrogen-starved yeast. We identified TF targets relevant for pseudohyphal growth, producing a detailed map of its regulatory network. Using tools from graph theory, we identified 14 TFs that lie at the center of this network, including Flo8, Mss11, and Mfg1, which bind as a complex. Surprisingly, the DNA-binding preferences for these key TFs were unknown. Using Calling Card data, we predicted the in vivo DNA-binding motif for the Flo8-Mss11-Mfg1 complex and validated it using a reporter assay. We found that this complex binds several important targets, including FLO11, at both their promoter and termination sequences. We demonstrated that this binding pattern is the result of DNA looping, which regulates the transcription of these targets and is stabilized by an interaction with the nuclear pore complex. This looping provides yeast cells with a transcriptional memory, enabling them more rapidly to execute the filamentous growth program when nitrogen starved if they had been previously exposed to this condition.

Epigenetic control of mobile DNA as an interface between experience and genome change

Mobile DNA in the genome is subject to RNA-targeted epigenetic control. This control regulates the activity of transposons, retrotransposons and genomic proviruses. Many different life history experiences alter the activities of mobile DNA and the expression of genetic loci regulated by nearby insertions. The same experiences induce alterations in epigenetic formatting and lead to trans-generational modifications of genome expression and stability. These observations lead to the hypothesis that epigenetic formatting directed by non-coding RNA provides a molecular interface between life history events and genome alteration.

Expansion of microsatellites on evolutionary young Y chromosome

Sex chromosomes are an ideal system to study processes connected with suppressed recombination. We found evidence of microsatellite expansion, on the relatively young Y chromosome of the dioecious plant sorrel (Rumex acetosa, XY1Y2 system), but no such expansion on the more ancient Y chromosomes of liverwort (Marchantia polymorpha) and human. The most expanding motifs were AC and AAC, which also showed periodicity of array length, indicating the importance of beginnings and ends of arrays. Our data indicate that abundance of microsatellites in genomes depends on the inherent expansion potential of specific motifs, which could be related to their stability and ability to adopt unusual DNA conformations. We also found that the abundance of microsatellites is higher in the neighborhood of transposable elements (TEs) suggesting that microsatellites are probably targets for TE insertions. This evidence suggests that microsatellite expansion is an early event shaping the Y chromosome where this process is not opposed by recombination, while accumulation of TEs and chromosome shrinkage predominate later.

Two large-scale analyses of Ty1 LTR-retrotransposon de novo insertion events indicate that Ty1 targets nucleosomal DNA near the H2A/H2B interface

BACKGROUND

Over the years, a number of reports have revealed that Ty1 integration occurs in a 1-kb window upstream of Pol III-transcribed genes with an approximate 80-bp periodicity between each integration hotspot and that this targeting requires active Pol III transcription at the site of integration. However, the molecular bases of Ty1 targeting are still not understood.

FINDINGS

The publications by Baller et al. and Mularoni et al. in the April issue of Genome Res. report the first high-throughput sequencing analysis of Ty1 de novo insertion events. Their observations converge to the same conclusion, that Ty1 targets a specific surface of the nucleosome at he H2A/H2B interface.

CONCLUSION

This discovery is important, and should help identifying factor(s) involved in Ty1 targeting. Recent data on transposable elements and retroviruses integration site choice obtained by large-scale analyses indicate that transcription and chromatin structure play an important role in this process. The studies reported in this commentary add a new evidence of the importance of chromatin in integration selectivity that should be of interest for everyone interested in transposable elements integration.

Evolutionary genomics of transposable elements in Saccharomyces cerevisiae

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.

Nucleosomes affect local transformation efficiency

Genetic transformation is a natural process during which foreign DNA enters a cell and integrates into the genome. Apart from its relevance for horizontal gene transfer in nature, transformation is also the cornerstone of today's recombinant gene technology. Despite its importance, relatively little is known about the factors that determine transformation efficiency. We hypothesize that differences in DNA accessibility associated with nucleosome positioning may affect local transformation efficiency. We investigated the landscape of transformation efficiency at various positions in the Saccharomyces cerevisiae genome and correlated these measurements with nucleosome positioning. We find that transformation efficiency shows a highly significant inverse correlation with relative nucleosome density. This correlation was lost when the nucleosome pattern, but not the underlying sequence was changed. Together, our results demonstrate a novel role for nucleosomes and also allow researchers to predict transformation efficiency of a target region and select spots in the genome that are likely to yield higher transformation efficiency.

Retrotransposon profiling of RNA polymerase III initiation sites

Although retroviruses are relatively promiscuous in choice of integration sites, retrotransposons can display marked integration specificity. In yeast and slime mold, some retrotransposons are associated with tRNA genes (tDNAs). In the Saccharomyces cerevisiae genome, the long terminal repeat retrotransposon Ty3 is found at RNA polymerase III (Pol III) transcription start sites of tDNAs. Ty1, 2, and 4 elements also cluster in the upstream regions of these genes. To determine the extent to which other Pol III-transcribed genes serve as genomic targets for Ty3, a set of 10,000 Ty3 genomic retrotranspositions were mapped using high-throughput DNA sequencing. Integrations occurred at all known tDNAs, two tDNA relics (iYGR033c and ZOD1), and six non-tDNA, Pol III-transcribed types of genes (RDN5, SNR6, SNR52, RPR1, RNA170, and SCR1). Previous work in vitro demonstrated that the Pol III transcription factor (TF) IIIB is important for Ty3 targeting. However, seven loci that bind the TFIIIB loader, TFIIIC, were not targeted, underscoring the unexplained absence of TFIIIB at those sites. Ty3 integrations also occurred in two open reading frames not previously associated with Pol III transcription, suggesting the existence of a small number of additional sites in the yeast genome that interact with Pol III transcription complexes.

A nucleosomal surface defines an integration hotspot for the Saccharomyces cerevisiae Ty1 retrotransposon

Ty1, the most abundant retrotransposon in Saccharomyces cerevisiae, integrates preferentially upstream of genes transcribed by RNA polymerase III (Pol III). Targeting is likely due to interactions between the Ty1 integration complex and a feature of chromatin characteristic of sites of Pol III transcription. To better understand Ty1 targeting determinants, >150,000 Ty1 insertions were mapped onto the S. cerevisiae genome sequence. Logistic regression was used to assess relationships between patterns of Ty1 integration and various genomic features, including genome-wide data sets of histone modifications and transcription-factor binding sites. Nucleosomes were positively associated with Ty1 insertions, and fine-scale mapping of insertions upstream of genes transcribed by Pol III indicated that Ty1 preferentially integrates into nucleosome-bound DNA near the H2A/H2B interface. Outside of genes transcribed by Pol III, Ty1 avoids coding sequences, a pattern that is not due to selection, but rather reflects a preference for nucleosome-rich sites flanking genes. Ty1 insertion sites were also mapped in four mutant lines that affect Ty1 transposition frequency or integration specificity (rrm3Δ, hos2Δ, rtt109Δ, and rad6Δ). Patterns of integration were largely preserved in the mutants, although significantly more insertions into coding sequences were observed in the rad6Δ strain, suggesting a loosening of target specificity in this mutant that lacks an enzyme involved in ubiquitinating H2A. Overall, our data suggest that nucleosomes are necessary for Ty1 integration, and that a secondary factor, likely a histone modification or nucleosome-bound factor enriched at sites of Pol III transcription, determines preferred target sites.