DNA sequence bias during Tn5 transposition

A Targeted Genome-scale Overexpression Platform for Proteobacteria

Targeted, genome-scale gene perturbation screens using Clustered Regularly Interspaced Short Palindromic Repeats interference (CRISPRi) and activation (CRISPRa) have revolutionized eukaryotic genetics, advancing medical, industrial, and basic research. Although CRISPRi knockdowns have been broadly applied in bacteria, options for genome-scale overexpression face key limitations. Here, we develop a facile approach for genome-scale gene overexpression in bacteria we call, "CRISPRtOE" (CRISPR transposition and OverExpression). We create a platform for comprehensive gene targeting using CRISPR-associated transposition (CAST) and show that transposition occurs at a higher frequency in non-transcribed DNA. We then demonstrate that CRISPRtOE can upregulate gene expression in Proteobacteria with medical and industrial relevance by integrating synthetic promoters of varying strength upstream of target genes. Finally, we employ CRISPRtOE screening at the genome-scale in Escherichia coli, recovering known antibiotic targets and genes with unexplored roles in antibiotic function. We envision that CRISPRtOE will be a valuable overexpression tool for antibiotic mode of action, industrial strain optimization, and gene function discovery in bacteria.

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.

High-throughput split-protein profiling by combining transposon mutagenesis and regulated protein-protein interactions with deep sequencing

Splitting a protein at a position may lead to self- or assisted-complementary fragments depending on whether two resulting fragments can reconstitute to maintain the native function spontaneously or require assistance from two interacting molecules. Assisted complementary fragments with high contrast are an important tool for probing biological interactions. However, only a small number of assisted-complementary split-variants have been identified due to manual, labour-intensive optimization of a candidate gene. Here, we introduce a technique for high-throughput split-protein profiling (HiTS) that allows fast identification of self- and assisted complementary positions by transposon mutagenesis, a rapamycin-regulated FRB-FKBP protein interaction pair, and deep sequencing. We test this technique by profiling three antibiotic-resistant genes (fosfomycin-resistant gene, fosA3, erythromycin-resistant gene, ermB, and chloramphenicol-resistant gene, catI). Self- and assisted complementary fragments discovered by the high-throughput technique were subsequently confirmed by low-throughput testing of individual split positions. Thus, the HiTS technique provides a quicker alternative for discovering the proteins with suitable self- and assisted-complementary split positions when combining with a readout such as fluorescence, bioluminescence, cell survival, gene transcription or genome editing.

Analysis of Gene Essentiality from TnSeq Data Using Transit

TnSeq, or sequencing of transposon insertion libraries, has proven to be a valuable method for probing the functions of genes in a wide range of bacteria. TnSeq has found many applications for studying genes involved in core functions (such as cell division or metabolism), stress response, virulence, etc., as well as to identify potential drug targets. Two of the most commonly used transposons in practice are Himar1, which inserts randomly at TA dinucleotides, and Tn5, which can insert more broadly throughout the genome. These insertions cause putative gene function disruption, and clones with insertions in genes that cannot tolerate disruption (in a given condition) are eliminated from the population. Deep sequencing can be used to efficiently profile the surviving members, with insertions in genes that can be inferred to be non-essential. Data from TnSeq experiments (i.e. transposon insertion counts at specific genomic locations) is inherently noisy, making rigorous statistical analysis (e.g. quantifying significance) challenging. In this chapter, we describe Transit, a Python-based software package for analyzing TnSeq data that combines a variety of data processing tools, quality assessment methods, and analytical algorithms for identifying essential (or conditionally essential) genes.

Fast-Seq: A Simple Method for Rapid and Inexpensive Validation of Packaged Single-Stranded Adeno-Associated Viral Genomes in Academic Settings

Adeno-associated viral (AAV) vectors have shown great promise in gene delivery as evidenced by recent FDA approvals. Despite efforts to optimize manufacturing for good manufacturing practice (GMP) productions, few academic laboratories have the resources to assess vector composition. One critical component of vector quality is packaged genome fidelity. Errors in viral genome replication and packaging can result in the incorporation of faulty genomes with mutations, truncations, or rearrangements, compromising vector potency. Thus, sequence validation of packaged genome composition is an important quality control (QC), even in academic settings. We developed Fast-Seq, an end-to-end method for extraction, purification, sequencing, and data analysis of packaged single-stranded AAV (ssAAV) genomes intended for non-GMP preclinical environments. We validated Fast-Seq on ssAAV vectors with three different genome compositions (CAG-GFP, CAG-tdTomato, EF1α-FLuc), three different genome sizes (2.9, 3.6, 4.4 kb), packaged in four different capsid serotypes (AAV1, AAV2, AAV5, and AAV8), and produced using the two most common production methods (Baculovirus-Sf9 and human HEK293), from both common commercial vendors and academic core facilities supplying academic laboratories. We achieved an average genome coverage of >1,400 × and an average inverted terminal repeat coverage of >280 × , despite the many differences in composition of each ssAAV sample. When compared with other ssAAV next-generation sequencing (NGS) methods for GMP settings, Fast-Seq has several unique advantages: Tn5 transposase-based fragmentation rather than sonication, 125 × less input DNA, simpler adapter ligation, compatibility with commonly available inexpensive sequencing instruments, and free open-source data analysis code in a preassembled customizable Docker container designed for novices. Fast-Seq can be completed in 18 h, is more cost-effective than other NGS methods, and is more accurate than Sanger sequencing, which is generally only applied at 1-2 × sequencing depth. Fast-Seq is a rapid, simple, and inexpensive methodology to validate packaged ssAAV genomes in academic settings.

A Novel Next-Generation Sequencing and Analysis Platform to Assess the Identity of Recombinant Adeno-Associated Viral Preparations from Viral DNA Extracts

Recombinant adeno-associated virus (rAAV) vectors are increasingly popular gene delivery tools in biological systems. They are safe and lead to high-level, long-term transgene expression. rAAV are available in multiple serotypes, natural or engineered, which enable targeting to a wide array of tissues and cell types. In addition, rAAVs are relatively easily produced in a well-equipped lab or obtained from a viral vector core facility. Unfortunately, there is no standardization of quality control assays beyond titering and purity assessments. Next-generation sequencing (NGS) can be used to identify rAAV preparations. Because the rAAV genome is single stranded, previous studies have assumed that rAAV genomes must be converted to double strands before NGS. We demonstrate that rAAV DNA extracts exist primarily as double-stranded species. We hypothesize that these molecules form from the natural base pairing of complementary [+] and [−] strands after DNA extraction and show that rAAV DNA extracts are sufficient templates for downstream NGS without the labor-intensive double-stranding step. Here, we provide a detailed protocol for the simple and rapid NGS of rAAV genomes from DNA extracts. With this protocol, users can quickly confirm the identity of an rAAV preparation and detect the presence of contaminating rAAV DNA. In addition, we share custom Python scripts that allow users to accurately determine the serotype and detect Cre-independent DNA recombination events in rAAV containing Lox sites within minutes. We have used these scripts to analyze more than 100 rAAV preparations. Although we focused on the detection of cross-contaminating rAAV DNA and recombination events, our Python scripts can be customized to detect other sequences or events, such as reverse packaging of plasmid backbone or DNA from the packaging cell line. We find that the NGS of rAAV DNA extracts, termed viral genome sequencing, is a simple and powerful method for rAAV validation.

Targeted insertional mutagenesis libraries for deep domain insertion profiling

Domain recombination is a key principle in protein evolution and protein engineering, but inserting a donor domain into every position of a target protein is not easily experimentally accessible. Most contemporary domain insertion profiling approaches rely on DNA transposons, which are constrained by sequence bias. Here, we establish Saturated Programmable Insertion Engineering (SPINE), an unbiased, comprehensive, and targeted domain insertion library generation technique using oligo library synthesis and multi-step Golden Gate cloning. Through benchmarking to MuA transposon-mediated library generation on four ion channel genes, we demonstrate that SPINE-generated libraries are enriched for in-frame insertions, have drastically reduced sequence bias as well as near-complete and highly-redundant coverage. Unlike transposon-mediated domain insertion that was severely biased and sparse for some genes, SPINE generated high-quality libraries for all genes tested. Using the Inward Rectifier K⁺ channel Kir2.1, we validate the practical utility of SPINE by constructing and comparing domain insertion permissibility maps. SPINE is the first technology to enable saturated domain insertion profiling. SPINE could help explore the relationship between domain insertions and protein function, and how this relationship is shaped by evolutionary forces and can be engineered for biomedical applications.

Statistical analysis of variability in TnSeq data across conditions using zero-inflated negative binomial regression

Background

Deep sequencing of transposon mutant libraries (or TnSeq) is a powerful method for probing essentiality of genomic loci under different environmental conditions. Various analytical methods have been described for identifying conditionally essential genes whose tolerance for insertions varies between two conditions. However, for large-scale experiments involving many conditions, a method is needed for identifying genes that exhibit significant variability in insertions across multiple conditions.

Results

In this paper, we introduce a novel statistical method for identifying genes with significant variability of insertion counts across multiple conditions based on Zero-Inflated Negative Binomial (ZINB) regression. Using likelihood ratio tests, we show that the ZINB distribution fits TnSeq data better than either ANOVA or a Negative Binomial (in a generalized linear model). We use ZINB regression to identify genes required for infection of M. tuberculosis H37Rv in C57BL/6 mice. We also use ZINB to perform a analysis of genes conditionally essential in H37Rv cultures exposed to multiple antibiotics.

Conclusions

Our results show that, not only does ZINB generally identify most of the genes found by pairwise resampling (and vastly out-performs ANOVA), but it also identifies additional genes where variability is detectable only when the magnitudes of insertion counts are treated separately from local differences in saturation, as in the ZINB model.

Defining the core essential genome of Pseudomonas aeruginosa

Genomics offered the promise of transforming antibiotic discovery by revealing many new essential genes as good targets, but the results fell short of the promise. While numerous factors contributed to the disappointing yield, one factor was that essential genes for a bacterial species were often defined based on a single or limited number of strains grown under a single or limited number of in vitro laboratory conditions. In fact, the essentiality of a gene can depend on both the genetic background and growth condition. We thus developed a strategy for more rigorously defining the core essential genome of a bacterial species by studying many pathogen strains and growth conditions. We assessed how many strains must be examined to converge on a set of core essential genes for a species. We used transposon insertion sequencing (Tn-Seq) to define essential genes in nine strains of Pseudomonas aeruginosa on five different media and developed a statistical model, FiTnEss, to classify genes as essential versus nonessential across all strain-medium combinations. We defined a set of 321 core essential genes, representing 6.6% of the genome. We determined that analysis of four strains was typically sufficient in P. aeruginosa to converge on a set of core essential genes likely to be essential across the species across a wide range of conditions relevant to in vivo infection, and thus to represent attractive targets for novel drug discovery.

Genetic and epigenetic determinants of inter-individual variability in responses to toxicants

It is well established that genetic variability has a major impact on susceptibility to common diseases, responses to drugs and toxicants, and influences disease-related outcomes. The appreciation that epigenetic marks also vary across the population is growing with more data becoming available from studies in humans and model organisms. In addition, the links between genetic variability, toxicity outcomes and epigenetics are being actively explored. Recent studies demonstrate that gene-by-environment interactions involve both chromatin states and transcriptional regulation, and that epigenetics provides important mechanistic clues to connect expression-related quantitative trait loci (QTL) and disease outcomes. However, studies of Gene×Environment×Epigenetics further extend the complexity of the experimental designs and create a challenge for selecting the most informative epigenetic readouts that can be feasibly performed to interrogate multiple individuals, exposures, tissue types and toxicity phenotypes. We propose that among the many possible epigenetic experimental methodologies, assessment of chromatin accessibility coupled with total RNA levels provides a cost-effective and comprehensive option to sufficiently characterize the complexity of epigenetic and regulatory activity in the context of understanding the inter-individual variability in responses to toxicants.

Haplotype phasing of whole human genomes using bead-based barcode partitioning in a single tube

Haplotype-resolved genome sequencing promises to unlock a wealth of information in population and medical genetics. However, for the vast majority of genomes sequenced to date, haplotypes have not been determined because of cumbersome haplotyping workflows that require fractions of the genome to be sequenced in a large number of compartments. Here we demonstrate barcode partitioning of long DNA molecules in a single compartment using "on-bead" barcoded tagmentation. The key to the method that we call "contiguity preserving transposition" sequencing on beads (CPTv2-seq) is transposon-mediated transfer of homogenous populations of barcodes from beads to individual long DNA molecules that get fragmented at the same time (tagmentation). These are then processed to sequencing libraries wherein all sequencing reads originating from each long DNA molecule share a common barcode. Single-tube, bulk processing of long DNA molecules with ∼150,000 different barcoded bead types provides a barcode-linked read structure that reveals long-range molecular contiguity. This technology provides a simple, rapid, plate-scalable and automatable route to accurate, haplotype-resolved sequencing, and phasing of structural variants of the genome.

Illumina error profiles: resolving fine-scale variation in metagenomic sequencing data

Background

Illumina’s sequencing platforms are currently the most utilised sequencing systems worldwide. The technology has rapidly evolved over recent years and provides high throughput at low costs with increasing read-lengths and true paired-end reads. However, data from any sequencing technology contains noise and our understanding of the peculiarities and sequencing errors encountered in Illumina data has lagged behind this rapid development.

Results

We conducted a systematic investigation of errors and biases in Illumina data based on the largest collection of in vitro metagenomic data sets to date. We evaluated the Genome Analyzer II, HiSeq and MiSeq and tested state-of-the-art low input library preparation methods. Analysing in vitro metagenomic sequencing data allowed us to determine biases directly associated with the actual sequencing process. The position- and nucleotide-specific analysis revealed a substantial bias related to motifs (3mers preceding errors) ending in “GG”. On average the top three motifs were linked to 16 % of all substitution errors. Furthermore, a preferential incorporation of ddGTPs was recorded. We hypothesise that all of these biases are related to the engineered polymerase and ddNTPs which are intrinsic to any sequencing-by-synthesis method. We show that quality-score-based error removal strategies can on average remove 69 % of the substitution errors - however, the motif-bias remains.

Conclusion

Single-nucleotide polymorphism changes in bacterial genomes can cause significant changes in phenotype, including antibiotic resistance and virulence, detecting them within metagenomes is therefore vital. Current error removal techniques are not designed to target the peculiarities encountered in Illumina sequencing data and other sequencing-by-synthesis methods, causing biases to persist and potentially affect any conclusions drawn from the data. In order to develop effective diagnostic and therapeutic approaches we need to be able to identify systematic sequencing errors and distinguish these errors from true genetic variation.

Electronic supplementary material

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

Low-Cost, High-Throughput Sequencing of DNA Assemblies Using a Highly Multiplexed Nextera Process

In recent years, next-generation sequencing (NGS) technology has greatly reduced the cost of sequencing whole genomes, whereas the cost of sequence verification of plasmids via Sanger sequencing has remained high. Consequently, industrial-scale strain engineers either limit the number of designs or take short cuts in quality control. Here, we show that over 4000 plasmids can be completely sequenced in one Illumina MiSeq run for less than $3 each (15× coverage), which is a 20-fold reduction over using Sanger sequencing (2× coverage). We reduced the volume of the Nextera tagmentation reaction by 100-fold and developed an automated workflow to prepare thousands of samples for sequencing. We also developed software to track the samples and associated sequence data and to rapidly identify correctly assembled constructs having the fewest defects. As DNA synthesis and assembly become a centralized commodity, this NGS quality control (QC) process will be essential to groups operating high-throughput pipelines for DNA construction.

Impact of three Illumina library construction methods on GC bias and HLA genotype calling

Next-generation sequencing (NGS) is increasingly recognized for its ability to overcome allele ambiguity and deliver high-resolution typing in the HLA system. Using this technology, non-uniform read distribution can impede the reliability of variant detection, which renders high-confidence genotype calling particularly difficult to achieve in the polymorphic HLA complex. Recently, library construction has been implicated as the dominant factor in instigating coverage bias. To study the impact of this phenomenon on HLA genotyping, we performed long-range PCR on 12 samples to amplify HLA-A, -B, -C, -DRB1, and -DQB1, and compared the relative contribution of three Illumina library construction methods (TruSeq Nano, Nextera, Nextera XT) in generating downstream bias. Here, we show high GC% to be a good predictor of low sequencing depth. Compared to standard TruSeq Nano, GC bias was more prominent in transposase-based protocols, particularly Nextera XT, likely through a combination of transposase insertion bias being coupled with a high number of PCR enrichment cycles. Importantly, our findings demonstrate non-uniform read depth can have a direct and negative impact on the robustness of HLA genotyping, which has clinical implications for users when choosing a library construction strategy that aims to balance cost and throughput with data quality.

Characterization of plasmid pPO1 from the hyperacidophile Picrophilus oshimae

Picrophilus oshimae and Picrophilus torridus are free-living, moderately thermophilic and acidophilic organisms from the lineage of Euryarchaeota. With a pH optimum of growth at pH 0.7 and the ability to even withstand molar concentrations of sulphuric acid, these organisms represent the most extreme acidophiles known. So far, nothing is known about plasmid biology in these hyperacidophiles. Also, there are no genetic tools available for this genus. We have mobilized the 7.6 Kbp plasmid from P. oshimae in E. coli by introducing origin-containing transposons and described the plasmid in terms of its nucleotide sequence, copy number in the native host, mode of replication, and transcriptional start sites of the encoded ORFs. Plasmid pPO1 may encode a restriction/modification system in addition to its replication functions. The information gained from the pPO1 plasmid may prove useful in developing a cloning system for this group of extreme acidophiles.

Transposon Tn5

Tn5 was one of the first transposons to be identified ( 10 ). As a result of Tn5's early discovery and its simple macromolecular requirements for transposition, the Tn5 system has been a very productive tool for studying the molecular mechanism of DNA transposition. These studies are of broad value because they offer insights into DNA transposition in general, because DNA transposition is a useful model with which to understand other types of protein-DNA interactions such as retroviral DNA integration and the DNA cleavage events involved in immunoglobulin gene formation, and because Tn5-derived tools are useful adjuncts in genetic experimentation.

Functional fusion proteins by random transposon-based GFP insertion

Fusions with fluorescent proteins are usually created by fusing the ends of two coding sequences. Appending the coding region of a fluorescent protein to the N- or C-terminus of another protein is typically the easiest way of creating a functional, fluorescent fusion protein. Another strategy involves placing the fluorescent protein in the middle of another protein. Such sandwich fusions are feasible, and there are many reasons for creating these fusion proteins. For example, sandwich fusions can be used to place two fluorescent proteins close to one another for optimization of a biosensor based on fluorescence resonance energy transfer, or they can be used to place the fluorescent protein in a region that moves during conformational changes of the host protein. Designing a sandwich fusion that produces a functional, fluorescent fusion protein is often challenging. This protocol describes an alternative approach. A simple, in vitro, transposon reaction is used to randomly insert the sequence encoding a fluorescent fusion protein into a target protein. This random labeling strategy makes it possible to create a small library of sandwich fusion proteins that can then be screened for activity. The approach makes it possible to test many possible solutions to the complex problem of building new biosensors.

A bifunctional DNA binding region in Tn5 transposase

Tn5 transposition is a complicated process that requires the formation of a highly ordered protein-DNA structure, a synaptic complex, to catalyse the movement of a sequence of DNA (transposon) into a target DNA. Much is known about the structure of the synaptic complex and the positioning of protein-DNA contacts, although many protein-DNA contacts remain largely unstudied. In particular, there is little evidence for the positioning of donor DNA and target DNA. In this communication, we describe the isolation and analysis of mutant transposases that have, for the first time, provided genetic and biochemical evidence for the stage-specific positioning of both donor and target DNAs within the synaptic complex. Furthermore, we have provided evidence that some of the amino acids that contact donor DNA also contact target DNA, and therefore suggest that these amino acids help define a bifunctional DNA binding region responsible for these two transposase-DNA binding events.

Transposon mutagenesis identifies sites upstream of the Neisseria gonorrhoeae pilE gene that modulate pilin antigenic variation

Gene conversion mediates the variation of virulence-associated surface structures on pathogenic microorganisms, which prevents host humoral immune responses from being effective. One of the best-studied gene conversion systems is antigenic variation (Av) of the pilin subunit of the Neisseria gonorrhoeae type IV pilus. To identify cis-acting DNA sequences that facilitate Av, the 700-bp region upstream of the pilin gene pilE was targeted for transposon mutagenesis. Four classes of transposon-associated mutations were isolated, distinguishable by their pilus-associated phenotypes: (i) insertions that did not alter Av or piliation, (ii) insertions that blocked Av, (iii) insertions that interfered with Av, and (iv) insertions that interfered with pilus expression and Av. Mutagenesis of the pilE promoter did not affect the frequency of Av, directly demonstrating that pilin Av is independent of pilE transcription. Two stretches of sequence upstream of pilE were devoid of transposon insertions, and some deletions in these regions were not recoverable, suggesting that they are essential for gonococcal viability. Insertions that blocked pilin Av were located downstream of the RS1 repeat sequence, and deletion of the region surrounding these insertions completely abrogated pilin Av, confirming that specific sequences 5' to pilE are essential for the recombination events underlying pilin Av.

Formation, characterization and partial purification of a Tn5 strand transfer complex

DNA transposition reactions typically involve a strand transfer step wherein the transposon ends are covalently joined by the transposase protein to a short target site. There is very little known about the transposase-DNA interactions that direct this process, and thus our overall understanding of the dynamics of DNA transposition reactions is limited. Tn5 presents an attractive system for defining such interactions because it has been possible to solve the structure of at least one Tn5 transposition intermediate: a transpososome formed with pre-cleaved ends. However, insertion specificity in the Tn5 system is low and this has hampered progress in generating target-containing transpososomes that are homogeneous in structure (i.e. where a single target site is engaged) and therefore suitable for biochemical and structural analysis. We have developed a system where the Tn5 transpososome integrates almost exclusively into a single target site within a short DNA fragment. The key to establishing this high degree of insertion specificity was to use a target DNA with tandem repeats of a previously characterized Tn5 insertion hotspot. The target DNA requirements to form this strand transfer complex are evaluated. In addition, we show that target DNAs missing single phosphate groups at specific positions are better substrates for strand transfer complex formation relative to the corresponding unmodified DNA fragments. Moreover, utilization of missing phosphate substrates can increase the degree of target site selection. A method for concentrating and partially purifying the Tn5 strand transfer complex is described.

Targeting Tn5 transposase identifies human immunodeficiency virus type 1 inhibitors

Human immunodeficiency virus (HIV) type 1 (HIV-1) integrase is an underutilized drug target for the treatment of HIV infection. One limiting factor is the lack of costructural data for use in the rational design or modification of integrase inhibitors. Tn5 transposase is a structurally well characterized, related protein that may serve as a useful surrogate. However, little data exist on inhibitor cross-reactivity. Here we screened 16,000 compounds using Tn5 transposase as the target and identified 20 compounds that appear to specifically inhibit complex assembly. Six were found to also inhibit HIV-1 integrase. These compounds likely interact with a highly conserved region presumably within the catalytic core. Most promising, several cinnamoyl derivatives were found to inhibit HIV transduction in cells. The identification of integrase inhibitors from a screen using Tn5 transposase as the target illustrates the utility of Tn5 as a surrogate for HIV-1 integration even though the relationship between the two systems is limited to the active site architecture and catalytic mechanism.