DNA transposon Hermes inserts into DNA in nucleosome-free regions in vivo

Tn-seq of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin

C. glabrata is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable Tn-seq methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131 kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1 . Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the drugging of which by SDZ 90-215 or myriocin enhanced the potency of micafungin in vitro . These findings provide insights into complexity of the C. glabrata species as well as strategies for improving antifungal efficacy.

Summary

Candida glabrata is an emerging pathogen with large genetic diversity and genome plasticity. The type strain CBS138 and a laboratory derivative were mutagenized with the Hermes transposon and profiled using Tn-seq. Numerous genes that regulate innate and acquired resistance to an important clinical antifungal were uncovered, including a pleiotropic drug resistance gene (PDR1) and a duplication of part of one chromosome. Compounds that target PDR1 and other genes may augment the potency of existing antifungals.

Minimal requirements for the epigenetic inheritance of engineered silent chromatin domains

Mechanisms enabling genetically identical cells to differentially regulate gene expression are complex and central to organismal development and evolution. While gene silencing pathways involving DNA sequence-specific recruitment of histone-modifying enzymes are prevalent in nature, examples of sequence-independent heritable gene silencing are scarce. Studies of the fission yeast Schizosaccharomyces pombe indicate that sequence-independent propagation of heterochromatin can occur but requires numerous multisubunit protein complexes and their diverse activities. Such complexity has so far precluded a coherent articulation of the minimal requirements for heritable gene silencing by conventional in vitro reconstitution approaches. Here, we take an unconventional approach to defining these requirements by engineering sequence-independent silent chromatin inheritance in budding yeast Saccharomyces cerevisiae cells. The mechanism conferring memory upon these cells is remarkably simple and requires only two proteins, one that recognizes histone H3 lysine 9 methylation (H3K9me) and catalyzes the deacetylation of histone H4 lysine 16 (H4K16), and another that recognizes deacetylated H4K16 and catalyzes H3K9me. Together, these bilingual "read-write" proteins form an interdependent positive feedback loop that is sufficient for the transmission of DNA sequence-independent silent information over multiple generations.

Unraveling the palindromic and nonpalindromic motifs of retroviral integration site sequences by statistical mixture models

A weak palindromic nucleotide motif is the hallmark of retroviral integration site alignments. Given that the majority of target sequences are not palindromic, the current model explains the symmetry by an overlap of the nonpalindromic motif present on one of the half-sites of the sequences. Here, we show that the implementation of multicomponent mixture models allows for different interpretations consistent with the existence of both palindromic and nonpalindromic submotifs in the sets of integration site sequences. We further show that the weak palindromic motifs result from freely combined site-specific submotifs restricted to only a few positions proximal to the site of integration. The submotifs are formed by either palindrome-forming nucleotide preference or nucleotide exclusion. Using the mixture models, we also identify HIV-1-favored palindromic sequences in Alu repeats serving as local hotspots for integration. The application of the novel statistical approach provides deeper insight into the selection of retroviral integration sites and may prove to be a valuable tool in the analysis of any type of DNA motifs.

Copy number variation alters local and global mutational tolerance

Copy number variants (CNVs), duplications and deletions of genomic sequences, contribute to evolutionary adaptation but can also confer deleterious effects and cause disease. Whereas the effects of amplifying individual genes or whole chromosomes (i.e., aneuploidy) have been studied extensively, much less is known about the genetic and functional effects of CNVs of differing sizes and structures. Here, we investigated Saccharomyces cerevisiae (yeast) strains that acquired adaptive CNVs of variable structures and copy numbers following experimental evolution in glutamine-limited chemostats. Although beneficial in the selective environment, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genome-wide genetic interactions in CNV strains. We find that CNVs increase mutational target size, confer increased mutational tolerance in amplified essential genes, and result in novel genetic interactions with unlinked genes. We validated a novel genetic interaction between different CNVs and BMH1 that was common to multiple strains. We also analyzed global gene expression and found that transcriptional dosage compensation does not affect most genes amplified by CNVs, although gene-specific transcriptional dosage compensation does occur for ∼12% of amplified genes. Furthermore, we find that CNV strains do not show previously described transcriptional signatures of aneuploidy. Our study reveals the extent to which local and global mutational tolerance is modified by CNVs with implications for genome evolution and CNV-associated diseases, such as cancer.

Zinc-finger BED domains drive the formation of the active Hermes transpososome by asymmetric DNA binding

The Hermes DNA transposon is a member of the eukaryotic hAT superfamily, and its transposase forms a ring-shaped tetramer of dimers. Our investigation, combining biochemical, crystallography and cryo-electron microscopy, and in-cell assays, shows that the full-length Hermes octamer extensively interacts with its transposon left-end through multiple BED domains of three Hermes protomers contributed by three dimers explaining the role of the unusual higher-order assembly. By contrast, the right-end is bound to no BED domains at all. Thus, this work supports a model in which Hermes multimerizes to gather enough BED domains to find its left-end among the abundant genomic DNA, facilitating the subsequent interaction with the right-end.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

Genome-wide analysis of heat stress-stimulated transposon mobility in the human fungal pathogen Cryptococcus deneoformans

We recently reported transposon mutagenesis as a significant driver of spontaneous mutations in the human fungal pathogen Cryptococcus deneoformans during murine infection. Mutations caused by transposable element (TE) insertion into reporter genes were dramatically elevated at high temperatures (37° vs. 30°) in vitro, suggesting that heat stress stimulates TE mobility in the Cryptococcus genome. To explore the genome-wide impact of TE mobilization, we generated transposon accumulation lines by in vitro passage of C. deneoformans strain XL280α for multiple generations at both 30° and at the host-relevant temperature of 37°. Utilizing whole-genome sequencing, we identified native TE copies and mapped multiple de novo TE insertions in these lines. Movements of the T1 DNA transposon occurred at both temperatures with a strong bias for insertion between gene-coding regions. By contrast, the Tcn12 retrotransposon integrated primarily within genes and movement occurred exclusively at 37°. In addition, we observed a dramatic amplification in copy number of the Cnl1 (Cryptococcus neoformans LINE-1) retrotransposon in subtelomeric regions under heat-stress conditions. Comparing TE mutations to other sequence variations detected in passaged lines, the increase in genomic changes at elevated temperatures was primarily due to mobilization of the retroelements Tcn12 and Cnl1. Finally, we found multiple TE movements (T1, Tcn12, and Cnl1) in the genomes of single C. deneoformans isolates recovered from infected mice, providing evidence that mobile elements are likely to facilitate microevolution and rapid adaptation during infection.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

A multiplexed, three-dimensional pooling and next-generation sequencing strategy for creating barcoded mutant arrays: construction of a Schizosaccharomyces pombe transposon insertion library

Arrayed libraries of defined mutants have been used to elucidate gene function in the post-genomic era. Yeast haploid gene deletion libraries have pioneered this effort, but are costly to construct, do not reveal phenotypes that may occur with partial gene function and lack essential genes required for growth. We therefore devised an efficient method to construct a library of barcoded insertion mutants with a wider range of phenotypes that can be generalized to other organisms or collections of DNA samples. We developed a novel but simple three-dimensional pooling and multiplexed sequencing approach that leveraged sequence information to reduce the number of required sequencing reactions by orders of magnitude, and were able to identify the barcode sequences and DNA insertion sites of 4391 Schizosaccharomyces pombe insertion mutations with only 40 sequencing preparations. The insertion mutations are in the genes and untranslated regions of nonessential, essential and noncoding RNA genes, and produced a wider range of phenotypes compared to the cognate deletion mutants, including novel phenotypes. This mutant library represents both a proof of principle for an efficient method to produce novel mutant libraries and a valuable resource for the S. pombe research community.

Loss-of-function mutation survey revealed that genes with background-dependent fitness are rare and functionally related in yeast

In different individuals, the same mutation can lead to different phenotypes due to genetic background effects. This is commonly observed in various systems, including many human diseases. While isolated examples of such background effects have been observed, a systematic view across a large number of individuals is still lacking. Here, we surveyed genetic background effects associated with gene loss-of-function mutations across a population of natural isolates of the yeast Saccharomyces cerevisiae. We found that ∼15% of genes can display a background-dependent fitness change. Genes related to mitochondrial functions are significantly overrepresented, and showed reversed patterns of fitness gain or loss with genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport, suggesting a potential functional rewiring.

In natural populations, the same mutation can lead to different phenotypic outcomes due to the genetic variation that exists among individuals. Such genetic background effects are commonly observed, including in the context of many human diseases. However, systematic characterization of these effects at the species level is still lacking to date. Here, we sought to comprehensively survey background-dependent traits associated with gene loss-of-function (LoF) mutations in 39 natural isolates of Saccharomyces cerevisiae using a transposon saturation strategy. By analyzing the modeled fitness variability of a total of 4,469 genes, we found that 15% of them, when impacted by a LoF mutation, exhibited a significant gain- or loss-of-fitness phenotype in certain natural isolates compared with the reference strain S288C. Out of these 632 genes with predicted background-dependent fitness effects, around 2/3 impact multiple backgrounds with a gradient of predicted fitness change while 1/3 are specific to a single genetic background. Genes related to mitochondrial function are significantly overrepresented in the set of genes showing a continuous variation and display a potential functional rewiring with other genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport. Such rewiring effects are likely modulated by both the genetic background and the environment. While background-specific cases are rare and span diverse cellular processes, they can be functionally related at the individual level. All genes with background-dependent fitness effects tend to have an intermediate connectivity in the global genetic interaction network and have shown relaxed selection pressure at the population level, highlighting their potential evolutionary characteristics.

A phylogenetic and proteomic reconstruction of eukaryotic chromatin evolution

Histones and associated chromatin proteins have essential functions in eukaryotic genome organization and regulation. Despite this fundamental role in eukaryotic cell biology, we lack a phylogenetically-comprehensive understanding of chromatin evolution. Here, we combine comparative proteomics and genomics analysis of chromatin in eukaryotes and archaea. Proteomics uncovers the existence of histone post-translational modifications in Archaea. However, archaeal histone modifications are scarce, in contrast with the highly conserved and abundant marks we identify across eukaryotes. Phylogenetic analysis reveals that chromatin-associated catalytic functions (e.g., methyltransferases) have pre-eukaryotic origins, whereas histone mark readers and chaperones are eukaryotic innovations. We show that further chromatin evolution is characterized by expansion of readers, including capture by transposable elements and viruses. Overall, our study infers detailed evolutionary history of eukaryotic chromatin: from its archaeal roots, through the emergence of nucleosome-based regulation in the eukaryotic ancestor, to the diversification of chromatin regulators and their hijacking by genomic parasites.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

Mapping of Nonhomologous End Joining-Mediated Integration Facilitates Genome-Scale Trackable Mutagenesis in Yarrowia lipolytica

Genome-scale mutagenesis, phenotypic screening, and tracking the causal mutations is a powerful approach for genetic analysis. However, classic mutagenesis approaches require extensive effort to identify causal mutations. It is desirable to demonstrate a powerful approach for rapid trackable mutagenesis. Here, we mapped the distribution of nonhomologous end joining (NHEJ)-mediated integration for the first time and demonstrated that it can be used for constructing the genome-scale trackable mutagenesis library in Yarrowia lipolytica. The sequencing of 9.15 × 10⁵ insertions showed that NHEJ-mediated integration inserted DNA randomly across the chromosomes, and the transcriptional regulatory regions exhibited integration preference. The insertions were located in both nucleosome-occupancy regions and nucleosome-free regions. Using NHEJ-mediated integration to construct the genome-scale mutagenesis library, the new targets that improved β-carotene biosynthesis and acetic acid tolerance were identified rapidly. This mutagenesis approach is readily applicable to other organisms with strong NHEJ preference and will contribute to cell factory construction.

Insertion Specificity of the hATx-6 Transposase of Hydra magnipapillata

Transposable elements (TE) are mobile genetic elements, present in all domains of life. They commonly encode a single transposase enzyme, that performs the excision and reintegration reactions, and these enzymes have been used in mutagenesis and creation of next-generation sequencing libraries. All transposases have some bias in the DNA sequence they bind to when reintegrating the TE DNA. We sought to identify a transposase that showed minimal sequence bias and could be produced recombinantly, using information from the literature and a novel bioinformatic analysis, resulting in the selection of the hATx-6 transposase from Hydra vulgaris (aka Hydra magnipapillata) for further study. This transposase was tested and shown to be active both in vitro and in vivo, and we were able to demonstrate very low sequence bias in its integration preference. This transposase could be an excellent candidate for use in biotechnology, such as the creation of next-generation sequencing libraries.

Mapping open chromatin by ATAC-seq in bread wheat

Gene transcription is largely regulated by cis-regulatory elements. Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is an emerging technology that can accurately map cis-regulatory elements in animals and plants. However, the presence of cell walls and chloroplasts in plants hinders the extraction of high-quality nuclei, thereby affects the quality of ATAC-seq data. Meanwhile, it is tricky to perform ATAC-seq with different tissue types, especially for those with limited size and amount. Moreover, with rapid growth of ATAC-seq datasets from plants, powerful and easy-to-use data analysis pipelines for ATAC-seq, especially for wheat is lacking. Here, we provided an all-in-one solution for mapping open chromatin in wheat including both experimental and data analysis procedure. We efficiently obtained nuclei with less cell debris from various wheat tissues. High-quality ATAC-seq data from young spike and ovary, which are hard to harvest were generated. We determined that the saturation sequencing depth of wheat ATAC-seq is about 16 Gb. Particularly, we developed a powerful and easy-to-use online pipeline to analyze the wheat ATAC-seq data and this pipeline can be easily extended to other plant species. The method developed here will facilitate plant regulatory genome study not only for wheat but also for other plant species.

OUP accepted manuscript

The Sleeping Beauty (SB) transposon system is a popular tool for genome engineering, but random integration into the genome carries a certain genotoxic risk in therapeutic applications. Here we investigate the role of amino acids H187, P247 and K248 in target site selection of the SB transposase. Structural modeling implicates these three amino acids located in positions analogous to amino acids with established functions in target site selection in retroviral integrases and transposases. Saturation mutagenesis of these residues in the SB transposase yielded variants with altered target site selection properties. Transposon integration profiling of several mutants reveals increased specificity of integrations into palindromic AT repeat target sequences in genomic regions characterized by high DNA bendability. The H187V and K248R mutants redirect integrations away from exons, transcriptional regulatory elements and nucleosomal DNA in the human genome, suggesting enhanced safety and thus utility of these SB variants in gene therapy applications.

Hijacking Transposable Elements for Saturation Mutagenesis in Fungi

Transposable elements are present in almost all known genomes, these endogenous transposons have recently been referred to as the mobilome. They are now increasingly used in research in order to make extensive mutant libraries in different organisms. Fungi are an essential part of our lives on earth, they influence the availability of our food and they live inside our own bodies both as commensals and pathogenic organisms. Only few fungal species have been studied extensively, mainly due to the lack of appropriate molecular genetic tools. The use of transposon insertion libraries can however help to rapidly advance our knowledge of (conditional) essential genes, compensatory mutations and drug target identification in fungi. Here we give an overview of some recent developments in the use of different transposons for saturation mutagenesis in different fungi.

Proteasome Control of [URE3] Prion Propagation by Degradation of Anti-Prion Proteins Cur1 and Btn2 in Saccharomyces cerevisiae

[URE3] is a prion of the nitrogen catabolism controller, Ure2p, and [PSI+] is a prion of the translation termination factor Sup35p in S. cerevisiae. Btn2p cures [URE3] by sequestration of Ure2p amyloid filaments. Cur1p, paralogous to Btn2p, also cures [URE3], but by a different (unknown) mechanism. We find that an array of mutations impairing proteasome assembly or MG132 inhibition of proteasome activity result in loss of [URE3]. In proportion to their prion-curing effects, each mutation affecting proteasomes elevates the cellular concentration of the anti-prion proteins Btn2 and Cur1. Of >4,600 proteins detected by SILAC, Btn2p was easily the most overexpressed in a pre9Δ (α3 core subunit) strain. Indeed, deletion of BTN2 and CUR1 prevents the prion-curing effects of proteasome impairment. Surprisingly, the 15 most unstable yeast proteins are not increased in pre9Δ cells suggesting altered proteasome specificity rather than simple inactivation. Hsp42, a chaperone that cooperates with Btn2 and Cur1 in curing [URE3], is also necessary for the curing produced by proteasome defects, although Hsp42p levels are not substantially altered by a proteasome defect. We find that pre9Δ and proteasome chaperone mutants that most efficiently lose [URE3], do not destabilize [PSI+] or alter cellular levels of Sup35p. A tof2 mutation or deletion likewise destabilizes [URE3], and elevates Btn2p, suggesting that Tof2p deficiency inactivates proteasomes. We suggest that when proteasomes are saturated with denatured/misfolded proteins, their reduced degradation of Btn2p and Cur1p automatically upregulates these aggregate-handling systems to assist in the clean-up.

The Integration Preference of Sleeping Beauty at Non-TA Site Is Related to the Transposon End Sequences

Recently, we proved that Sleeping Beauty (SB) transposon integrates into non-TA sites at a lower frequency. Here, we performed a further study on the non-TA integration of SB and showed that (1) SB can integrate into non-TA sites in HEK293T cells as well as in mouse cell lines; (2) Both the hyperactive transposase SB100X and the traditional SB11 catalyze integrations at non-TA sites; (3) The consensus sequence of the non-TA target sites only occurs at the opposite side of the sequenced junction between the transposon end and the genomic sequences, indicating that the integrations at non-TA sites are mainly aberrant integrations; and (4) The consensus sequence of the non-TA target sites is corresponding to the transposon end sequence. The consensus sequences changed following the changes of the transposon ends. This result indicated that the interaction between the SB transposon end and genomic DNA (gDNA) may be involved in the target site selection of the SB integrations at non-TA sites.

Iterative epigenomic analyses in the same single cell

Gene expression in individual cells is epigenetically regulated by DNA modifications, histone modifications, transcription factors, and other DNA-binding proteins. It has been shown that multiple histone modifications can predict gene expression and reflect future responses of bulk cells to extracellular cues. However, the predictive ability of epigenomic analysis is still limited for mechanistic research at a single cell level. To overcome this limitation, it would be useful to acquire reliable signals from multiple epigenetic marks in the same single cell. Here, we propose a new approach and a new method for analysis of several components of the epigenome in the same single cell. The new method allows reanalysis of the same single cell. We found that reanalysis of the same single cell is feasible, provides confirmation of the epigenetic signals, and allows application of statistical analysis to identify reproduced reads using data sets generated only from the single cell. Reanalysis of the same single cell is also useful to acquire multiple epigenetic marks from the same single cells. The method can acquire at least five epigenetic marks: H3K27ac, H3K27me3, mediator complex subunit 1, a DNA modification, and a DNA-interacting protein. We can predict active signaling pathways in K562 single cells using the epigenetic data and confirm that the predicted results strongly correlate with actual active signaling pathways identified by RNA-seq results. These results suggest that the new method provides mechanistic insights for cellular phenotypes through multilayered epigenome analysis in the same single cells.

Innate immunity to yeast prions: Btn2p and Cur1p curing of the [URE3] prion is prevented by 60S ribosomal protein deficiency or ubiquitin/proteasome system overactivity

[URE3] is an amyloid-based prion of Ure2p, a negative regulator of poor nitrogen source catabolism in Saccharomyces cerevisiae. Overproduced Btn2p or its paralog Cur1p, in processes requiring Hsp42, cure the [URE3] prion. Btn2p cures by collecting Ure2p amyloid filaments at one place in the cell. We find that rpl4aΔ, rpl21aΔ, rpl21bΔ, rpl11bΔ, and rpl16bΔ (large ribosomal subunit proteins) or ubr2Δ (ubiquitin ligase targeting Rpn4p, an activator of proteasome genes) reduce curing by overproduced Btn2p or Cur1p. Impaired curing in ubr2Δ or rpl21bΔ is restored by an rpn4Δ mutation. No effect of rps14aΔ or rps30bΔ on curing was observed, indicating that 60S subunit deficiency specifically impairs curing. Levels of Hsp42p, Sis1p, or Btn3p are unchanged in rpl4aΔ, rpl21bΔ, or ubr2Δ mutants. Overproduction of Cur1p or Btn2p was enhanced in rpn4Δ and hsp42Δ mutants, lower in ubr2Δ strains, and restored to above wild-type levels in rpn4Δ ubr2Δ strains. As in the wild-type, Ure2N-GFP colocalizes with Btn2-RFP in rpl4aΔ, rpl21bΔ, or ubr2Δ strains, but not in hsp42Δ. Btn2p/Cur1p overproduction cures [URE3] variants with low seed number, but seed number is not increased in rpl4aΔ, rpl21bΔ or ubr2Δ mutants. Knockouts of genes required for the protein sorting function of Btn2p did not affect curing of [URE3], nor did inactivation of the Hsp104 prion-curing activity. Overactivity of the ubiquitin/proteasome system, resulting from 60S subunit deficiency or ubr2Δ, may impair Cur1p and Btn2p curing of [URE3] by degrading Cur1p, Btn2p or another component of these curing systems.

Identification of Essential Genes and Fluconazole Susceptibility Genes in Candida glabrata by Profiling Hermes Transposon Insertions

Within the budding yeasts, the opportunistic pathogen Candida glabrata and other members of the Nakaseomyces clade have developed virulence traits independently from C. albicans and C. auris. To begin exploring the genetic basis of C. glabrata virulence and its innate resistance to antifungals, we launched the Hermes transposon from a plasmid and sequenced more than 500,000 different semi-random insertions throughout the genome. With machine learning, we identified 1278 protein-encoding genes (25% of total) that could not tolerate transposon insertions and are likely essential for C. glabrata fitness in vitro. Interestingly, genes involved in mRNA splicing were less likely to be essential in C. glabrata than their orthologs in S. cerevisiae, whereas the opposite is true for genes involved in kinetochore function and chromosome segregation. When a pool of insertion mutants was challenged with the first-line antifungal fluconazole, insertions in several known resistance genes (e.g., PDR1, CDR1, PDR16, PDR17, UPC2A, DAP1, STV1) and 15 additional genes (including KGD1, KGD2, YHR045W) became hypersensitive to fluconazole. Insertions in 200 other genes conferred significant resistance to fluconazole, two-thirds of which function in mitochondria and likely down-regulate Pdr1 expression or function. Knockout mutants of KGD2 and IDH2, which consume and generate alpha-ketoglutarate in mitochondria, exhibited increased and decreased resistance to fluconazole through a process that depended on Pdr1. These findings establish the utility of transposon insertion profiling in forward genetic investigations of this important pathogen of humans.

Leveraging the Hermes Transposon to Accelerate the Development of Nonconventional Yeast-based Microbial Cell Factories

We broadened the usage of DNA transposon technology by demonstrating its capacity for the rapid creation of expression libraries for long biochemical pathways, which is beyond the classical application of building genome-scale knockout libraries in yeasts. This strategy efficiently leverages the readily available fine-tuning impact provided by the diverse transcriptional environment surrounding each random integration locus. We benchmark the transposon-mediated integration against the nonhomologous end joining-mediated strategy. The latter strategy was demonstrated for achieving pathway random integration in other yeasts but is associated with a high false-positive rate in the absence of a high-throughput screening method. Our key innovation of a nonreplicable circular DNA platform increased the possibility of identifying top-producing variants to 97%. Compared to the classical DNA transposition protocol, the design of a nonreplicable circular DNA skipped the step of counter-selection for plasmid removal and thus not only reduced the time required for the step of library creation from 10 to 5 d but also efficiently removed the "transposition escapers", which undesirably represented almost 80% of the entire population as false positives. Using one endogenous product (i.e., shikimate) and one heterologous product (i.e., (S)-norcoclaurine) as examples, we presented a streamlined procedure to rapidly identify high-producing variants with titers significantly higher than the reported data in the literature. We selected Scheffersomyces stipitis, a representative nonconventional yeast, as a demo, but the strategy can be generalized to other nonconventional yeasts. This new exploration of transposon technology, therefore, adds a highly versatile tool to accelerate the development of novel species as microbial cell factories for producing value-added chemicals.

Comparing the utility of in vivo transposon mutagenesis approaches in yeast species to infer gene essentiality

In vivo transposon mutagenesis, coupled with deep sequencing, enables large-scale genome-wide mutant screens for genes essential in different growth conditions. We analyzed six large-scale studies performed on haploid strains of three yeast species (Saccharomyces cerevisiae, Schizosaccaromyces pombe, and Candida albicans), each mutagenized with two of three different heterologous transposons (AcDs, Hermes, and PiggyBac). Using a machine-learning approach, we evaluated the ability of the data to predict gene essentiality. Important data features included sufficient numbers and distribution of independent insertion events. All transposons showed some bias in insertion site preference because of jackpot events, and preferences for specific insertion sequences and short-distance vs long-distance insertions. For PiggyBac, a stringent target sequence limited the ability to predict essentiality in genes with few or no target sequences. The machine learning approach also robustly predicted gene function in less well-studied species by leveraging cross-species orthologs. Finally, comparisons of isogenic diploid versus haploid S. cerevisiae isolates identified several genes that are haplo-insufficient, while most essential genes, as expected, were recessive. We provide recommendations for the choice of transposons and the inference of gene essentiality in genome-wide studies of eukaryotic haploid microbes such as yeasts, including species that have been less amenable to classical genetic studies.

Transposable elements employ distinct integration strategies with respect to transcriptional landscapes in eukaryotic genomes

Transposable elements (TEs) are ubiquitous DNA segments capable of moving from one site to another within host genomes. The extant distributions of TEs in eukaryotic genomes have been shaped by both bona fide TE integration preferences in eukaryotic genomes and by selection following integration. Here, we compare TE target site distribution in host genomes using multiple de novo transposon insertion datasets in both plants and animals and compare them in the context of genome-wide transcriptional landscapes. We showcase two distinct types of transcription-associated TE targeting strategies that suggest a process of convergent evolution among eukaryotic TE families. The integration of two precision-targeting elements are specifically associated with initiation of RNA Polymerase II transcription of highly expressed genes, suggesting the existence of novel mechanisms of precision TE targeting in addition to passive targeting of open chromatin. We also highlight two features that can facilitate TE survival and rapid proliferation: tissue-specific transposition and minimization of negative impacts on nearby gene function due to precision targeting.

Structural Insights on Retroviral DNA Integration: Learning from Foamy Viruses

Foamy viruses (FV) are retroviruses belonging to the Spumaretrovirinae subfamily. They are non-pathogenic viruses endemic in several mammalian hosts like non-human primates, felines, bovines, and equines. Retroviral DNA integration is a mandatory step and constitutes a prime target for antiretroviral therapy. This activity, conserved among retroviruses and long terminal repeat (LTR) retrotransposons, involves a viral nucleoprotein complex called intasome. In the last decade, a plethora of structural insights on retroviral DNA integration arose from the study of FV. Here, we review the biochemistry and the structural features of the FV integration apparatus and will also discuss the mechanism of action of strand transfer inhibitors.

Fitness Landscape of the Fission Yeast Genome

The relationship between DNA sequence, biochemical function, and molecular evolution is relatively well-described for protein-coding regions of genomes, but far less clear in noncoding regions, particularly, in eukaryote genomes. In part, this is because we lack a complete description of the essential noncoding elements in a eukaryote genome. To contribute to this challenge, we used saturating transposon mutagenesis to interrogate the Schizosaccharomyces pombe genome. We generated 31 million transposon insertions, a theoretical coverage of 2.4 insertions per genomic site. We applied a five-state hidden Markov model (HMM) to distinguish insertion-depleted regions from insertion biases. Both raw insertion-density and HMM-defined fitness estimates showed significant quantitative relationships to gene knockout fitness, genetic diversity, divergence, and expected functional regions based on transcription and gene annotations. Through several analyses, we conclude that transposon insertions produced fitness effects in 66-90% of the genome, including substantial portions of the noncoding regions. Based on the HMM, we estimate that 10% of the insertion depleted sites in the genome showed no signal of conservation between species and were weakly transcribed, demonstrating limitations of comparative genomics and transcriptomics to detect functional units. In this species, 3'- and 5'-untranslated regions were the most prominent insertion-depleted regions that were not represented in measures of constraint from comparative genomics. We conclude that the combination of transposon mutagenesis, evolutionary, and biochemical data can provide new insights into the relationship between genome function and molecular evolution.

Structural insights into the mechanism of double strand break formation by Hermes, a hAT family eukaryotic DNA transposase

Some DNA transposons relocate from one genomic location to another using a mechanism that involves generating double-strand breaks at their transposon ends by forming hairpins on flanking DNA. The same double-strand break mode is employed by the V(D)J recombinase at signal-end/coding-end junctions during the generation of antibody diversity. How flanking hairpins are formed during DNA transposition has remained elusive. Here, we describe several co-crystal structures of the Hermes transposase bound to DNA that mimics the reaction step immediately prior to hairpin formation. Our results reveal a large DNA conformational change between the initial cleavage step and subsequent hairpin formation that changes which strand is acted upon by a single active site. We observed that two factors affect the conformational change: the complement of divalent metal ions bound by the catalytically essential DDE residues, and the identity of the –2 flanking base pair. Our data also provides a mechanistic link between the efficiency of hairpin formation (an A:T basepair is favored at the –2 position) and Hermes' strong target site preference. Furthermore, we have established that the histidine residue within a conserved C/DxxH motif present in many transposase families interacts directly with the scissile phosphate, suggesting a crucial role in catalysis.

Old cogs, new tricks: the evolution of gene expression in a chromatin context

Sophisticated gene-regulatory mechanisms probably evolved in prokaryotes billions of years before the emergence of modern eukaryotes, which inherited the same basic enzymatic machineries. However, the epigenomic landscapes of eukaryotes are dominated by nucleosomes, which have acquired roles in genome packaging, mitotic condensation and silencing parasitic genomic elements. Although the molecular mechanisms by which nucleosomes are displaced and modified have been described, just how transcription factors, histone variants and modifications and chromatin regulators act on nucleosomes to regulate transcription is the subject of considerable ongoing study. We explore the extent to which these transcriptional regulatory components function in the context of the evolutionarily ancient role of chromatin as a barrier to processes acting on DNA and how chromatin proteins have diversified to carry out evolutionarily recent functions that accompanied the emergence of differentiation and development in multicellular eukaryotes.

Validating genome-wide CRISPR-Cas9 function improves screening in the oleaginous yeast Yarrowia lipolytica

Genome-wide mutational screens are central to understanding the genetic underpinnings of evolved and engineered phenotypes. The widespread adoption of CRISPR-Cas9 genome editing has enabled such screens in many organisms, but identifying functional sgRNAs still remains a challenge. Here, we developed a methodology to quantify the cutting efficiency of each sgRNA in a genome-scale library, and in doing so improve screens in the biotechnologically important yeast Yarrowia lipolytica. Screening in the presence and absence of native DNA repair enabled high-throughput quantification of sgRNA function leading to the identification of high efficiency sgRNAs that cover 94% of genes. Library validation enhanced the classification of essential genes by identifying inactive guides that create false negatives and mask the effects of successful disruptions. Quantification of guide effectiveness also creates a dataset from which determinants of CRISPR-Cas9 can be identified. Finally, application of the library identified novel mutations for metabolic engineering of high lipid accumulation.

Historical Meiotic Crossover Hotspots Fueled Patterns of Evolutionary Divergence in Rice

Recombination plays an integral role in the creation of novel genetic variation in sexually reproducing species. Despite this important role, the determinants and evolution of crossover hotspots have remained poorly understood in plants. Here, we present a comparative analysis of two rice (Oryza sativa) historical recombination maps from two subspecies (indica and japonica) using 150 resequenced genomes. Fine-scale recombination rates and crossover hotspots were validated by comparison with a consensus genetic map and empirically derived crossovers, respectively. Strikingly, nearly 80% of crossover hotspots were unique to each subspecies, despite their relatively recent divergence and broad-scale correlated recombination rates. Crossover hotspots were enriched with Stowaway and P instability factor (PIF)/Harbinger transposons and overlapped accessible chromatin regions. Increased nucleotide diversity and signatures of population differentiation augmented by Stowaway and PIF/Harbinger transposons were prevalent at subspecies-specific crossover hotspots. Motifs derived from lineage-specific indica and japonica crossover hotspots were nearly identical in the two subspecies, implicating a core set of crossover motifs in rice. Finally, Stowaway and PIF/Harbinger transposons were associated with stabilized G/C bias within highly active hotspots, suggesting that hotspot activity can be fueled by de novo variation. These results provide evolutionary insight into historical crossover hotspots as potentially powerful drivers of sequence and subspecies evolution in plants.

Transposon insertional mutagenesis in Saccharomyces uvarum reveals trans-acting effects influencing species-dependent essential genes

To understand how complex genetic networks perform and regulate diverse cellular processes, the function of each individual component must be defined. Comprehensive phenotypic studies of mutant alleles have been successful in model organisms in determining what processes depend on the normal function of a gene. These results are often ported to newly sequenced genomes by using sequence homology. However, sequence similarity does not always mean identical function or phenotype, suggesting that new methods are required to functionally annotate newly sequenced species. We have implemented comparative analysis by high-throughput experimental testing of gene dispensability in Saccharomyces uvarum, a sister species of Saccharomyces cerevisiae. We created haploid and heterozygous diploid Tn7 insertional mutagenesis libraries in S. uvarum to identify species-dependent essential genes, with the goal of detecting genes with divergent functions and/or different genetic interactions. Comprehensive gene dispensability comparisons with S. cerevisiae predicted diverged dispensability at 12% of conserved orthologs, and validation experiments confirmed 22 differentially essential genes. Despite their differences in essentiality, these genes were capable of cross-species complementation, demonstrating that trans-acting factors that are background-dependent contribute to differential gene essentiality. This study shows that direct experimental testing of gene disruption phenotypes across species can inform comparative genomic analyses and improve gene annotations. Our method can be widely applied in microorganisms to further our understanding of genome evolution.

Gene Essentiality Analyzed by In Vivo Transposon Mutagenesis and Machine Learning in a Stable Haploid Isolate of Candida albicans

Comprehensive understanding of an organism requires that we understand the contributions of most, if not all, of its genes. Classical genetic approaches to this issue have involved systematic deletion of each gene in the genome, with comprehensive sets of mutants available only for very-well-studied model organisms. We took a different approach, harnessing the power of in vivo transposition coupled with deep sequencing to identify >500,000 different mutations, one per cell, in the prevalent human fungal pathogen Candida albicans and to map their positions across the genome. The transposition approach is efficient and less labor-intensive than classic approaches. Here, we describe the production and analysis (aided by machine learning) of a large collection of mutants and the comprehensive identification of 1,610 C. albicans genes that are essential for growth under standard laboratory conditions. Among these C. albicans essential genes, we identify those that are also essential in two distantly related model yeasts as well as those that are conserved in all four major human fungal pathogens and that are not conserved in the human genome. This list of genes with functions important for the survival of the pathogen provides a good starting point for the development of new antifungal drugs, which are greatly needed because of the emergence of fungal pathogens with elevated resistance and/or tolerance of the currently limited set of available antifungal drugs.

Knowing the full set of essential genes for a given organism provides important information about ways to promote, and to limit, its growth and survival. For many non-model organisms, the lack of a stable haploid state and low transformation efficiencies impede the use of conventional approaches to generate a genome-wide comprehensive set of mutant strains and the identification of the genes essential for growth. Here we report on the isolation and utilization of a highly stable haploid derivative of the human pathogenic fungus Candida albicans, together with a modified heterologous transposon and machine learning (ML) analysis method, to predict the degree to which all of the open reading frames are required for growth under standard laboratory conditions. We identified 1,610 C. albicans essential genes, including 1,195 with high “essentiality confidence” scores, thereby increasing the number of essential genes (currently 66 in the Candida Genome Database) by >20-fold and providing an unbiased approach to determine the degree of confidence in the determination of essentiality. Among the genes essential in C. albicans were 602 genes also essential in the model budding and fission yeasts analyzed by both deletion and transposon mutagenesis. We also identified essential genes conserved among the four major human pathogens C. albicans, Aspergillus fumigatus, Cryptococcus neoformans, and Histoplasma capsulatum and highlight those that lack homologs in humans and that thus could serve as potential targets for the design of antifungal therapies.

Trac-looping measures genome structure and chromatin accessibility

Long-range chromatin interactions play critical roles in genome organization and regulation of transcription. We now report transposase-mediated analysis of chromatin looping (Trac-looping) for simultaneous detection of multiscale genome-wide chromatin interactions among regulatory elements and chromatin accessibility. With this technique, a bivalent oligonucleotide linker is inserted between two interacting regions such that the chromatin interactions are captured without prior chromatin fragmentation and proximity-based ligation. Application of Trac-looping to human CD4⁺ T cells revealed substantial reorganization of enhancer-promoter interactions associated with changes in gene expression after T cell receptor stimulation.

Genome-Wide Determination of Gene Essentiality by Transposon Insertion Sequencing in Yeast Pichia pastoris

In many prokaryotes but limited eukaryotic species, the combination of transposon mutagenesis and high-throughput sequencing has greatly accelerated the identification of essential genes. Here we successfully applied this technique to the methylotrophic yeast Pichia pastoris and classified its conditionally essential/non-essential gene sets. Firstly, we showed that two DNA transposons, TcBuster and Sleeping beauty, had high transposition activities in P. pastoris. By merging their insertion libraries and performing Tn-seq, we identified a total of 202,858 unique insertions under glucose supported growth condition. We then developed a machine learning method to classify the 5,040 annotated genes into putatively essential, putatively non-essential, ambig1 and ambig2 groups, and validated the accuracy of this classification model. Besides, Tn-seq was also performed under methanol supported growth condition and methanol specific essential genes were identified. The comparison of conditionally essential genes between glucose and methanol supported growth conditions helped to reveal potential novel targets involved in methanol metabolism and signaling. Our findings suggest that transposon mutagenesis and Tn-seq could be applied in the methylotrophic yeast Pichia pastoris to classify conditionally essential/non-essential gene sets. Our work also shows that determining gene essentiality under different culture conditions could help to screen for novel functional components specifically involved in methanol metabolism.

Functional genomics for the oleaginous yeast Yarrowia lipolytica

Oleaginous yeasts are valuable systems for biosustainable production of hydrocarbon-based chemicals. Yarrowia lipolytica is one of the best characterized of these yeast with respect to genome annotation and flux analysis of metabolic processes. Nonetheless, progress is hampered by a dearth of genome-wide tools enabling functional genomics. In order to remedy this deficiency, we developed a library of Y. lipolytica insertion mutants via transposon mutagenesis. The Hermes DNA transposon was expressed to achieve saturation mutagenesis of the genome. Over 534,000 independent insertions were identified by next-generation sequencing. Poisson analysis of insertion density classified ~ 22% of genes as essential. As expected, most essential genes have homologs in Saccharomyces cerevisiae and Schizosaccharomyces pombe, and the majority of those are also essential. As an obligate aerobe, Y. lipolytica has significantly more respiration - related genes that are classified as essential than do S. cerevisiae and S. pombe. Contributions of non-essential genes to growth in glucose and glycerol carbon sources were assessed and used to evaluate two recent genome-scale models of Y. lipolytica metabolism. Fluorescence-activated cell sorting identified mutants in which lipid accumulation is increased. Our findings provide insights into biosynthetic pathways, compartmentalization of enzymes, and distinct functions of paralogs. This functional genomic analysis of the oleaginous yeast Y. lipolytica provides an important resource for modeling, bioengineering, and design of synthetic minimalized strains of respiratory yeasts.

Hermes Transposon Mutagenesis Shows [URE3] Prion Pathology Prevented by a Ubiquitin-Targeting Protein: Evidence for Carbon/Nitrogen Assimilation Cross Talk and a Second Function for Ure2p in Saccharomyces cerevisiae

[URE3] is an amyloid-based prion of Ure2p, a regulator of nitrogen catabolism. While most "variants" of the [URE3] prion are toxic, mild variants that only slightly slow growth are more widely studied. The existence of several antiprion systems suggests that some components may be protecting cells from potential detrimental effects of mild [URE3] variants. Our extensive Hermes transposon mutagenesis showed that disruption of YLR352W dramatically slows the growth of [URE3-1] strains. Ylr352wp is an F-box protein, directing selection of substrates for ubiquitination by a "cullin"-containing E₃ ligase. For efficient ubiquitylation, cullin-dependent E₃ ubiquitin ligases must be NEDDylated, modified by a ubiquitin-related peptide called NEDD8 (Rub1p in yeast). Indeed, we find that disruption of NEDDylation-related genes RUB1, ULA1, UBA3, and UBC12 is also counterselected in our screen. We find that like ylr352wΔ [URE3] strains, ylr352wΔ ure2Δ strains do not grow on nonfermentable carbon sources. Overexpression of Hap4p, a transcription factor stimulating expression of mitochondrial proteins, or mutation of GLN1, encoding glutamine synthetase, allows growth of ylr352w∆ [URE3] strains on glycerol media. Supplying proline as a nitrogen source shuts off the nitrogen catabolite repression (NCR) function of Ure2p, but does not slow growth of ylr352wΔ strains, suggesting a distinct function of Ure2p in carbon catabolism. Also, gln1 mutations impair NCR, but actually relieve the growth defect of ylr352wΔ [URE3] and ylr352wΔ ure2Δ strains, again showing that loss of NCR is not producing the growth defect and suggesting that Ure2p has another function. YLR352W largely protects cells from the deleterious effects of otherwise mild [URE3] variants or of a ure2 mutation (the latter a rarer event), and we name it LUG1 (lets [URE3]/ure2 grow).

Sleeping Beauty transposon integrates into non-TA dinucleotides

Background

Sleeping Beauty transposon (SB) has become an increasingly important genetic tool for generating mutations in vertebrate cells. It is widely thought that SB exclusively integrates into TA dinucleotides. However, this strict TA-preference has not been rigorously tested in large numbers of insertion sites that now can be detected with next generation sequencing. Li et al. found 71 SB insertions in non-TA dinucleotides in 2013, suggesting that TA dinucleotides are not the only sites of SB integration, yet further studies on this topic have not been carried out.

Results

In this study, we re-analyzed 600 million pairs of Illumina sequence reads from a high-throughput SB mutagenesis screen and identified 28 thousand SB insertions in non-TA sites. We recovered some of these non-TA sites using PCR and confirmed that at least a subset of the insertions at non-TA sites are real integrations. The consensus sequence of these non-TA sites shows an asymmetric pattern distinct from the symmetric pattern of the canonical TA sites. Perfect similarity between the downstream flanking sequence and SB transposon ends indicates there may be interaction between the transposon DNA binding domain of transposase and the target DNA.

Conclusion

The TA-preference of SB transposon is not as strict as what people had thought. And the SB integrations at non-TA sites might be guided by the interaction between the transposon DNA binding domain of SB transposase and the target DNA.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0113-8) contains supplementary material, which is available to authorized users.

Profiling Developmentally and Environmentally Controlled Chromatin Reprogramming

Dynamic reshuffling of the chromatin landscape is a recurrent theme orchestrated in many, if not all, plant developmental transitions and adaptive responses. Spatiotemporal variations of the chromatin properties on regulatory genes and on structural genomic elements trigger the establishment of distinct transcriptional contexts, which in some instances can epigenetically be inherited. Studies on plant cell plasticity during the differentiation of stem cells, including gametogenesis, or the specialization of vegetative cells in various organs, as well as the investigation of allele-specific gene regulation have long been impaired by technical challenges in generating specific chromatin profiles in complex or hardly accessible cell populations. Recent advances in increasing the sensitivity of genome-enabled technologies and in the isolation of specific cell types have allowed for overcoming such limitations. These developments hint at multilevel regulatory events ranging from nucleosome accessibility and composition to higher order chromatin organization and genome topology. Uncovering the large extent to which chromatin dynamics and epigenetic processes influence gene expression is therefore not surprisingly revolutionizing current views on plant molecular genetics and (epi)genomics as well as their perspectives in eco-evolutionary biology. Here, we introduce current methodologies to probe genome-wide chromatin variations for which protocols are detailed in this book chapter, with an emphasis on the plant model species Arabidopsis.

Functional mapping of yeast genomes by saturated transposition

Yeast is a powerful model for systems genetics. We present a versatile, time- and labor-efficient method to functionally explore the Saccharomyces cerevisiae genome using saturated transposon mutagenesis coupled to high-throughput sequencing. SAturated Transposon Analysis in Yeast (SATAY) allows one-step mapping of all genetic loci in which transposons can insert without disrupting essential functions. SATAY is particularly suited to discover loci important for growth under various conditions. SATAY (1) reveals positive and negative genetic interactions in single and multiple mutant strains, (2) can identify drug targets, (3) detects not only essential genes, but also essential protein domains, (4) generates both null and other informative alleles. In a SATAY screen for rapamycin-resistant mutants, we identify Pib2 (PhosphoInositide-Binding 2) as a master regulator of TORC1. We describe two antagonistic TORC1-activating and -inhibiting activities located on opposite ends of Pib2. Thus, SATAY allows to easily explore the yeast genome at unprecedented resolution and throughput.

DOI:http://dx.doi.org/10.7554/eLife.23570.001

Genes are stretches of DNA that carry the instructions to build and maintain cells. Many studies in genetics involve inactivating one or more genes and observing the consequences. If the loss of a gene kills the cell, that gene is likely to be vital for life. If it does not, the gene may not be essential, or a similar gene may be able to take over its role.

Baker’s yeast is a simple organism that shares many characteristics with human cells. Many yeast genes have a counterpart among human genes, and so studying baker’s yeast can reveal clues about our own genetics. Michel et al. report an adaptation for baker’s yeast of a technique called “Transposon sequencing”, which had been used in other single-celled organisms to study the effects of interrupting genes. Briefly, a virus-like piece of DNA, called a transposon, inserts randomly into the genetic material and switches off individual genes. The DNA is then sequenced to reveal every gene that can be disrupted without killing the cell, and remaining genes are inferred to be essential for life.

The approach, named SATAY (which is short for “saturated transposon analysis in yeast”), uses this strategy to create millions of baker’s yeast cells, each with a different gene switched off. Because the number of cells generated this way vastly exceeds the number of genes, every gene will be switched off by several independent transposons. Therefore the technique allows all yeast genes to be inactivated several times in one single experiment. The cells can be grown in varying conditions during the experiment, revealing the genes needed for survival in different situations. Non-essential genes can also be inactivated beforehand to uncover if any genes might be compensating for their absence.

In the future, this technique may be used to better understand human diseases, such as cancer, since many disease-causing genes in humans have counterparts in yeast.

DOI:http://dx.doi.org/10.7554/eLife.23570.002

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.

Improved genome sequencing using an engineered transposase

Background

Next-generation sequencing (NGS) has transformed genomic research by reducing turnaround time and cost. However, no major breakthrough has been made in the upstream library preparation methods until the transposase-based Nextera method was invented. Nextera combines DNA fragmentation and barcoding in a single tube reaction and therefore enables a very fast workflow to sequencing-ready DNA libraries within a couple of hours. When compared to the traditional ligation-based methods, transposed-based Nextera has a slight insertion bias.

Results

Here we present the discovery of a mutant transposase (Tn5-059) with a lowered GC insertion bias through protein engineering. We demonstrate Tn5-059 reduces AT dropout and increases uniformity of genome coverage in both bacterial genomes and human genome. We also observe higher library diversity generated by Tn5-059 when compared to Nextera v2 for human exomes, which leads to less sequencing and lower cost per genome. In addition, when used for human exomes, Tn5-059 delivers consistent library insert size over a range of input DNA, allowing up to a tenfold variance from the 50 ng input recommendation.

Conclusions

Enhanced DNA input tolerance of Tn5-059 can translate to flexibility and robustness of workflow. DNA input tolerance together with superior uniformity of coverage and lower AT dropouts extend the applications of transposase based library preps. We discuss possible mechanisms of improvements in Tn5-059, and potential advantages of using the new mutant in varieties of applications including microbiome sequencing and chromatin profiling.

Electronic supplementary material

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

Retroviruses integrate into a shared, non-palindromic DNA motif

Many DNA-binding factors, such as transcription factors, form oligomeric complexes with structural symmetry that bind to palindromic DNA sequences¹. Palindromic consensus nucleotide sequences are also found at the genomic integration sites of retroviruses^2-6 and other transposable elements^7-9, and it has been suggested that this palindromic consensus arises as a consequence of the structural symmetry in the integrase complex^2,3. However, we show here that the palindromic consensus sequence is not present in individual integration sites of human T-cell lymphotropic virus type 1 (HTLV-1) and human immunodeficiency virus type 1 (HIV-1), but arises in the population average as a consequence of the existence of a non-palindromic nucleotide motif that occurs in approximately equal proportions on the plus strand and the minus strand of the host genome. We develop a generally applicable algorithm to sort the individual integration site sequences into plus-strand and minus-strand subpopulations, and use this to identify the integration site nucleotide motifs of five retroviruses of different genera: HTLV-1, HIV-1, murine leukaemia virus (MLV), avian sarcoma leucosis virus (ASLV) and prototype foamy virus (PFV). The results reveal a non-palindromic motif that is shared between these retroviruses.

Mechanism for DNA transposons to generate introns on genomic scales

The discovery of introns four decades ago was one of the most unexpected findings in molecular biology. Introns are sequences interrupting genes that must be removed as part of messenger RNA production. Genome sequencing projects have shown that most eukaryotic genes contain at least one intron, and frequently many. Comparison of these genomes reveals a history of long evolutionary periods during which few introns were gained, punctuated by episodes of rapid, extensive gain. However, although several detailed mechanisms for such episodic intron generation have been proposed, none has been empirically supported on a genomic scale. Here we show how short, non-autonomous DNA transposons independently generated hundreds to thousands of introns in the prasinophyte Micromonas pusilla and the pelagophyte Aureococcus anophagefferens. Each transposon carries one splice site. The other splice site is co-opted from the gene sequence that is duplicated upon transposon insertion, allowing perfect splicing out of the RNA. The distributions of sequences that can be co-opted are biased with respect to codons, and phasing of transposon-generated introns is similarly biased. These transposons insert between pre-existing nucleosomes, so that multiple nearby insertions generate nucleosome-sized intervening segments. Thus, transposon insertion and sequence co-option may explain the intron phase biases and prevalence of nucleosome-sized exons observed in eukaryotes. Overall, the two independent examples of proliferating elements illustrate a general DNA transposon mechanism that can plausibly account for episodes of rapid, extensive intron gain during eukaryotic evolution.

Retroviral DNA Integration

The integration of a DNA copy of the viral RNA genome into host chromatin is the defining step of retroviral replication. This enzymatic process is catalyzed by the virus-encoded integrase protein, which is conserved among retroviruses and LTR-retrotransposons. Retroviral integration proceeds via two integrase activities: 3'-processing of the viral DNA ends, followed by the strand transfer of the processed ends into host cell chromosomal DNA. Herein we review the molecular mechanism of retroviral DNA integration, with an emphasis on reaction chemistries and architectures of the nucleoprotein complexes involved. We additionally discuss the latest advances on anti-integrase drug development for the treatment of AIDS and the utility of integrating retroviral vectors in gene therapy applications.

hAT Transposable Elements

hAT transposons are ancient in their origin and they are widespread across eukaryote kingdoms. They can be present in large numbers in many genomes. However, only a few active forms of these elements have so far been discovered indicating that, like all transposable elements, there is selective pressure to inactivate them. Nonetheless, there have been sufficient numbers of active hAT elements and their transposases characterized that permit an analysis of their structure and function. This review analyzes these and provides a comparison with the several domesticated hAT genes discovered in eukaryote genomes. Active hAT transposons have also been developed as genetic tools and understanding how these may be optimally utilized in new hosts will depend, in part, on understanding the basis of their function in genomes.

Continuous Influx of Genetic Material from Host to Virus Populations

Many genes of large double-stranded DNA viruses have a cellular origin, suggesting that host-to-virus horizontal transfer (HT) of DNA is recurrent. Yet, the frequency of these transfers has never been assessed in viral populations. Here we used ultra-deep DNA sequencing of 21 baculovirus populations extracted from two moth species to show that a large diversity of moth DNA sequences (n = 86) can integrate into viral genomes during the course of a viral infection. The majority of the 86 different moth DNA sequences are transposable elements (TEs, n = 69) belonging to 10 superfamilies of DNA transposons and three superfamilies of retrotransposons. The remaining 17 sequences are moth sequences of unknown nature. In addition to bona fide DNA transposition, we uncover microhomology-mediated recombination as a mechanism explaining integration of moth sequences into viral genomes. Many sequences integrated multiple times at multiple positions along the viral genome. We detected a total of 27,504 insertions of moth sequences in the 21 viral populations and we calculate that on average, 4.8% of viruses harbor at least one moth sequence in these populations. Despite this substantial proportion, no insertion of moth DNA was maintained in any viral population after 10 successive infection cycles. Hence, there is a constant turnover of host DNA inserted into viral genomes each time the virus infects a moth. Finally, we found that at least 21 of the moth TEs integrated into viral genomes underwent repeated horizontal transfers between various insect species, including some lepidopterans susceptible to baculoviruses. Our results identify host DNA influx as a potent source of genetic diversity in viral populations. They also support a role for baculoviruses as vectors of DNA HT between insects, and call for an evaluation of possible gene or TE spread when using viruses as biopesticides or gene delivery vectors.

While gene exchange is known to occur between viruses and their hosts, this phenomenon has never been studied at the level of the viral population. Here we report that each time a virus from the Baculoviridae family infects a moth, a large number (dozens to hundreds) and high diversity of moth DNA sequences (86 different sequences) can integrate into replicating viral genomes. These findings show that viral populations carry a measurable load of host DNA sequences, further supporting the role of viruses as vectors of horizontal transfer of DNA between insect species. The potential uncontrolled gene spread associated with the use of viruses produced in insect cells as gene delivery vectors and/or biopesticides should therefore be evaluated.

Genome-wide Profiling Reveals Remarkable Parallels Between Insertion Site Selection Properties of the MLV Retrovirus and the piggyBac Transposon in Primary Human CD4(+) T Cells

The inherent risks associated with vector insertion in gene therapy need to be carefully assessed. We analyzed the genome-wide distributions of Sleeping Beauty (SB) and piggyBac (PB) transposon insertions as well as MLV retrovirus and HIV lentivirus insertions in human CD4⁺ T cells with respect to a panel of 40 chromatin states. The distribution of SB transposon insertions displayed the least deviation from random, while the PB transposon and the MLV retrovirus showed unexpected parallels across all chromatin states. Both MLV and PB insertions are enriched at transcriptional start sites (TSSs) and co-localize with BRD4-associated sites. We demonstrate physical interaction between the PB transposase and bromodomain and extraterminal domain proteins (including BRD4), suggesting convergent evolution of a tethering mechanism that directs integrating genetic elements into TSSs. We detect unequal biases across the four systems with respect to targeting genes whose deregulation has been previously linked to serious adverse events in gene therapy clinical trials. The SB transposon has the highest theoretical chance of targeting a safe harbor locus in the human genome. The data underscore the significance of vector choice to reduce the mutagenic load on cells in clinical applications.

Transposable element detection from whole genome sequence data

The number of software tools available for detecting transposable element insertions from whole genome sequence data has been increasing steadily throughout the last ~5 years. Some of these methods have unique features suiting them for particular use cases, but in general they follow one or more of a common set of approaches. Here, detection and filtering approaches are reviewed in the light of transposable element biology and the current state of whole genome sequencing. We demonstrate that the current state-of-the-art methods still do not produce highly concordant results and provide resources to assist future development in transposable element detection methods.

Cell-based screens and phenomics with fission yeast

Next-generation sequencing approaches have considerably advanced our understanding of genome function and regulation. However, the knowledge of gene function and complex cellular processes remains a challenge and bottleneck in biological research. Phenomics is a rapidly emerging area, which seeks to rigorously characterize all phenotypes associated with genes or gene variants. Such high-throughput phenotyping under different conditions can be a potent approach toward gene function. The fission yeast Schizosaccharomyces pombe (S. pombe) is a proven eukaryotic model organism that is increasingly used for genomewide screens and phenomic assays. In this review, we highlight current large-scale, cell-based approaches used with S. pombe, including computational colony-growth measurements, genetic interaction screens, parallel profiling using barcodes, microscopy-based cell profiling, metabolomic methods and transposon mutagenesis. These diverse methods are starting to offer rich insights into the relationship between genotypes and phenotypes.