Tn5/IS50 target recognition

Variations in the persistence of 5'-end genomic and subgenomic SARS-CoV-2 RNAs in wastewater from aircraft, airports and wastewater treatment plants

Wastewater sequencing has become a powerful supplement to clinical testing in monitoring SARS-CoV-2 infections in the post-COVID-19 pandemic era. While its applications in measuring the viral burden and main circulating lineages in the community have proved their efficacy, the variations in sequencing quality and coverage across the different regions of the SARS-CoV-2 genome are not well understood. Furthermore, it is unclear how different sample origins, viral extraction and concentration methods and environmental factors impact the reads sequenced from wastewater. Using high-coverage, amplicon-based, paired-end read sequencing of viral RNA extracted from wastewater collected directly from aircraft, pooled from different aircraft and airport buildings or from regular wastewater plants, we assessed the genome coverage across the sample groups with a focus on the 5'-end region covering the leader sequence and investigated whether it was possible to detect subgenomic RNA from viral material recovered from wastewater. We identified distinct patterns in the persistence of the different genomic regions across the different types of wastewaters and the existence of chimeric reads mapping to non-amplified regions. Our findings suggest that preservation of the 5'-end of the genome and the ability to detect subgenomic RNA reads, though highly susceptible to environment and sample processing conditions, may be indicative of the quality and amount of the viral RNA present in wastewater.

Bacterial genome engineering using CRISPR-associated transposases

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated transposases have the potential to transform the technology landscape for kilobase-scale genome engineering, by virtue of their ability to integrate large genetic payloads with high accuracy, easy programmability and no requirement for homologous recombination machinery. These transposons encode efficient, CRISPR RNA-guided transposases that execute genomic insertions in Escherichia coli at efficiencies approaching ~100%. Moreover, they generate multiplexed edits when programmed with multiple guides, and function robustly in diverse Gram-negative bacterial species. Here we present a detailed protocol for engineering bacterial genomes using CRISPR-associated transposase (CAST) systems, including guidelines on the available vectors, customization of guide RNAs and DNA payloads, selection of common delivery methods, and genotypic analysis of integration events. We further describe a computational CRISPR RNA design algorithm to avoid potential off-targets, and a CRISPR array cloning pipeline for performing multiplexed DNA insertions. The method presented here allows the isolation of clonal strains containing a novel genomic integration event of interest within 1-2 weeks using available plasmid constructs and standard molecular biology techniques.

Bacterial genome engineering using CRISPR RNA-guided transposases

CRISPR-associated transposons (CASTs) have the potential to transform the technology landscape for kilobase-scale genome engineering, by virtue of their ability to integrate large genetic payloads with high accuracy, easy programmability, and no requirement for homologous recombination machinery. These transposons encode efficient, CRISPR RNA-guided transposases that execute genomic insertions in E. coli at efficiencies approaching ~100%, generate multiplexed edits when programmed with multiple guides, and function robustly in diverse Gram-negative bacterial species. Here we present a detailed protocol for engineering bacterial genomes using CAST systems, including guidelines on the available homologs and vectors, customization of guide RNAs and DNA payloads, selection of common delivery methods, and genotypic analysis of integration events. We further describe a computational crRNA design algorithm to avoid potential off-targets and CRISPR array cloning pipeline for DNA insertion multiplexing. Starting from available plasmid constructs, the isolation of clonal strains containing a novel genomic integration event-of-interest can be achieved in 1 week using standard molecular biology techniques.

Epigenetic Application of ATAC-Seq Based on Tn5 Transposase Purification Technology

Background

Assays of transposase accessible chromatin sequencing (ATAC-seq) is an efficient assay to investigate chromatin accessibility, which depends on the activity of a robust Tn5 transposase to fragment the genome while cutting in the sequencing adapters.

Methods

We propose reliable approaches for purifying hyperactive Tn5 transposase by chitin magnetic bead sorting. Double-stranded DNA of J76 cells and 293T cells were digested and subjected to tagmentation as test samples with Tn5 transposase, and libraries were established and sequenced. Sequencing data was then analyzed for peak calling, GO enrichment, and motif analysis.

Results

We report a set of rapid, efficient, and low-cost methods for ATAC-seq library construction and data analysis, through large-scale and rapid sequencing. These methods can provide a reference for the study of epigenetic regulation of gene expression.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing

Background

Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 μg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits.

Results

We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits’ fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection.

Conclusions

The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08316-y.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Modeling Site-Specific Nucleotide Biases Affecting Himar1 Transposon Insertion Frequencies in TnSeq Data Sets

TnSeq is a widely used methodology for determining gene essentiality, conditional fitness, and genetic interactions in bacteria. The Himar1 transposon is restricted to insertions at TA dinucleotides, but otherwise, few site-specific biases have been identified. As a result, most analytical approaches assume that insertions are expected to be randomly distributed among TA sites in nonessential regions. However, through analysis of Himar1 transposon libraries in Mycobacterium tuberculosis, we demonstrate that there are site-specific biases that affect the frequency of insertion of the Himar1 transposon at different TA sites. We use machine learning and statistical models to characterize patterns in the nucleotides surrounding TA sites that correlate with high or low insertion counts. We then develop a quantitative model based on these patterns that can be used to predict the expected counts at each TA site based on nucleotide context, which can explain up to half of the variance in insertion counts. We show that these insertion preferences exist in Himar1 TnSeq data sets from other mycobacterial and nonmycobacterial species. We present an improved method for identification of essential genes, called TTN-Fitness, that can better distinguish true biological fitness effects by comparing observed counts to expected counts based on our site-specific model of insertion preferences. Compared to previous essentiality methods, TTN-Fitness can make finer distinctions among genes whose disruption causes a fitness defect (or advantage), separating them out from the large pool of nonessentials, and is able to classify many smaller genes (with few TA sites) that were previously characterized as uncertain. IMPORTANCE When using the Himar1 transposon to create transposon insertion mutant libraries, it is known that the transposon is restricted to insertions at TA dinucleotide sites throughout the genome, and the absence of insertions is used to infer which genes are essential (or conditionally essential) in a bacterial organism. It is widely assumed that insertions in nonessential regions are otherwise random, and this assumption is used as the basis of several methods for statistical analysis of TnSeq data. In this paper, we show that the nucleotide sequence surrounding TA sites influences the magnitude of insertions, and these Himar1 insertion preferences (sequence biases) can partially explain why some sites have higher counts than others. We use this predictive model to make improved estimates of the fitness effects of genes, which help make finer distinctions of the phenotype and biological consequences of disruption of nonessential genes.

SorTn-seq: a high-throughput functional genomics approach to discovering regulators of bacterial gene expression

We recently developed a high-throughput functional genomics approach, named 'SorTn-seq', to identify factors affecting expression of any gene of interest in bacteria. Our approach facilitates high-throughput screening of complex mutant pools, a task previously hindered by a lack of suitable techniques. SorTn-seq combines high-density, Tn5-like transposon mutagenesis with fluorescence-activated cell sorting of a strain harboring a promoter-fluorescent reporter fusion, to isolate mutants with altered gene expression. The transposon mutant pool is sorted into different bins on the basis of fluorescence, and mutants are deep-sequenced to identify transposon insertions. DNA is prepared for sequencing by using commercial kits augmented with custom primers, enhancing ease of use and reproducibility. Putative regulators are identified by comparing the number of insertions per genomic feature in the different sort bins, by using existing bioinformatic pipelines and software packages. SorTn-seq can be completed in 1-2 weeks and requires general microbiology skills and basic flow cytometry experience.

The histone chaperone HIR maintains chromatin states to control nitrogen assimilation and fungal virulence

Adaptation to changing environments and immune evasion is pivotal for fitness of pathogens. Yet, the underlying mechanisms remain largely unknown. Adaptation is governed by dynamic transcriptional re-programming, which is tightly connected to chromatin architecture. Here, we report a pivotal role for the HIR histone chaperone complex in modulating virulence of the human fungal pathogen Candida albicans. Genetic ablation of HIR function alters chromatin accessibility linked to aberrant transcriptional responses to protein as nitrogen source. This accelerates metabolic adaptation and increases the release of extracellular proteases, which enables scavenging of alternative nitrogen sources. Furthermore, HIR controls fungal virulence, as HIR1 deletion leads to differential recognition by immune cells and hypervirulence in a mouse model of systemic infection. This work provides mechanistic insights into chromatin-coupled regulatory mechanisms that fine-tune pathogen gene expression and virulence. Furthermore, the data point toward the requirement of refined screening approaches to exploit chromatin modifications as antifungal strategies.

Fast and inexpensive whole-genome sequencing library preparation from intact yeast cells

Through the increase in the capacity of sequencing machines massively parallel sequencing of thousands of samples in a single run is now possible. With the improved throughput and resulting drop in the price of sequencing, the cost and time for preparation of sequencing libraries have become the major bottleneck in large-scale experiments. Methods using a hyperactive variant of the Tn5 transposase efficiently generate libraries starting from cDNA or genomic DNA in a few hours and are highly scalable. For genome sequencing, however, the time and effort spent on genomic DNA isolation limit the practicability of sequencing large numbers of samples. Here, we describe a highly scalable method for preparing high-quality whole-genome sequencing libraries directly from Saccharomyces cerevisiae cultures in less than 3 h at 34 cents per sample. We skip the rate-limiting step of genomic DNA extraction by directly tagmenting lysed yeast spheroplasts and add a nucleosome release step prior to enrichment PCR to improve the evenness of genomic coverage. Resulting libraries do not show any GC bias and are comparable in quality to libraries processed from genomic DNA with a commercially available Tn5-based kit. We use our protocol to investigate CRISPR/Cas9 on- and off-target edits and reliably detect edited variants and shared polymorphisms between strains. Our protocol enables rapid preparation of unbiased and high-quality, sequencing-ready indexed libraries for hundreds of yeast strains in a single day at a low price. By adjusting individual steps of our workflow, we expect that our protocol can be adapted to other organisms.

Evaluating coverage bias in next-generation sequencing of Escherichia coli

Whole-genome sequencing is essential to many facets of infectious disease research. However, technical limitations such as bias in coverage and tagmentation, and difficulties characterising genomic regions with extreme GC content have created significant obstacles in its use. Illumina has claimed that the recently released DNA Prep library preparation kit, formerly known as Nextera Flex, overcomes some of these limitations. This study aimed to assess bias in coverage, tagmentation, GC content, average fragment size distribution, and de novo assembly quality using both the Nextera XT and DNA Prep kits from Illumina. When performing whole-genome sequencing on Escherichia coli and where coverage bias is the main concern, the DNA Prep kit may provide higher quality results; though de novo assembly quality, tagmentation bias and GC content related bias are unlikely to improve. Based on these results, laboratories with existing workflows based on Nextera XT would see minor benefits in transitioning to the DNA Prep kit if they were primarily studying organisms with neutral GC content.

The landscape of chromatin accessibility in skeletal muscle during embryonic development in pigs

BACKGROUND

The development of skeletal muscle in pigs during the embryonic stage is precisely regulated by transcriptional mechanisms, which depend on chromatin accessibility. However, how chromatin accessibility plays a regulatory role during embryonic skeletal muscle development in pigs has not been reported. To gain insight into the landscape of chromatin accessibility and the associated genome-wide transcriptome during embryonic muscle development, we performed ATAC-seq and RNA-seq analyses of skeletal muscle from pig embryos at 45, 70 and 100 days post coitus (dpc).

RESULTS

In total, 21,638, 35,447 and 60,181 unique regions (or peaks) were found across the embryos at 45 dpc (LW45), 70 dpc (LW70) and 100 dpc (LW100), respectively. More than 91% of the peaks were annotated within - 1 kb to 100 bp of transcription start sites (TSSs). First, widespread increases in specific accessible chromatin regions (ACRs) from embryos at 45 to 100 dpc suggested that the regulatory mechanisms became increasingly complicated during embryonic development. Second, the findings from integrated ATAC-seq and RNA-seq analyses showed that not only the numbers but also the intensities of ACRs could control the expression of associated genes. Moreover, the motif screening of stage-specific ACRs revealed some transcription factors that regulate muscle development-related genes, such as MyoG, Mef2c, and Mef2d. Several potential transcriptional repressors, including E2F6, OTX2 and CTCF, were identified among the genes that exhibited different regulation trends between the ATAC-seq and RNA-seq data.

CONCLUSIONS

This work indicates that chromatin accessibility plays an important regulatory role in the embryonic muscle development of pigs and regulates the temporal and spatial expression patterns of key genes in muscle development by influencing the binding of transcription factors. Our results contribute to a better understanding of the regulatory dynamics of genes involved in pig embryonic skeletal muscle development.

Infection trains the host for microbiota-enhanced resistance to pathogens

The microbiota shields the host against infections in a process known as colonization resistance. How infections themselves shape this fundamental process remains largely unknown. Here, we show that gut microbiota from previously infected hosts display enhanced resistance to infection. This long-term functional remodeling is associated with altered bile acid metabolism leading to the expansion of taxa that utilize the sulfonic acid taurine. Notably, supplying exogenous taurine alone is sufficient to induce this alteration in microbiota function and enhance resistance. Mechanistically, taurine potentiates the microbiota's production of sulfide, an inhibitor of cellular respiration, which is key to host invasion by numerous pathogens. As such, pharmaceutical sequestration of sulfide perturbs the microbiota's composition and promotes pathogen invasion. Together, this work reveals a process by which the host, triggered by infection, can deploy taurine as a nutrient to nourish and train the microbiota, promoting its resistance to subsequent infection.

Recent Advances in Genetic Tools for Acinetobacter baumannii

Acinetobacter baumannii is classified as a top priority pathogen by the World Health Organization (WHO) because of its widespread resistance to all classes of antibiotics. This makes the need for understanding the mechanisms of resistance and virulence critical. Therefore, tools that allow genetic manipulations are vital to unravel the mechanisms of multidrug resistance (MDR) and virulence in A. baumannii. A host of current strategies are available for genetic manipulations of A. baumannii laboratory-strains, including ATCC^® 17978^TM and ATCC^® 19606^T, but depending on susceptibility profiles, these strategies may not be sufficient when targeting strains newly obtained from clinic, primarily due to the latter's high resistance to antibiotics that are commonly used for selection during genetic manipulations. This review highlights the most recent methods for genetic manipulation of A. baumannii including CRISPR based approaches, transposon mutagenesis, homologous recombination strategies, reporter systems and complementation techniques with the spotlight on those that can be applied to MDR clinical isolates.

IncC conjugative plasmids and SXT/R391 elements repair double-strand breaks caused by CRISPR-Cas during conjugation

Bacteria have evolved defence mechanisms against bacteriophages. Restriction-modification systems provide innate immunity by degrading invading DNAs that lack proper methylation. CRISPR-Cas systems provide adaptive immunity by sampling the genome of past invaders and cutting the DNA of closely related DNA molecules. These barriers also restrict horizontal gene transfer mediated by conjugative plasmids. IncC conjugative plasmids are important contributors to the global dissemination of multidrug resistance among pathogenic bacteria infecting animals and humans. Here, we show that IncC conjugative plasmids are highly resilient to host defence systems during entry into a new host by conjugation. Using a TnSeq strategy, we uncover a conserved operon containing five genes (vcrx089-vcrx093) that confer a novel host defence evasion (hde) phenotype. We show that vcrx089-vcrx090 promote resistance against type I restriction-modification, whereas vcrx091-vcxr093 promote CRISPR-Cas evasion by repairing double-strand DNA breaks via recombination between short sequence repeats. vcrx091, vcrx092 and vcrx093 encode a single-strand binding protein, and a single-strand annealing recombinase and double-strand exonuclease related to Redβ and λExo of bacteriophage λ, respectively. Homologous genes of the integrative and conjugative element R391 also provide CRISPR-Cas evasion. Hence, the conserved hde operon considerably broadens the host range of large families of mobile elements spreading multidrug resistance.

FASTQINS and ANUBIS: two bioinformatic tools to explore facts and artifacts in transposon sequencing and essentiality studies

Transposon sequencing is commonly applied for identifying the minimal set of genes required for cellular life; a major challenge in fields such as evolutionary or synthetic biology. However, the scientific community has no standards at the level of processing, treatment, curation and analysis of this kind data. In addition, we lack knowledge about artifactual signals and the requirements a dataset has to satisfy to allow accurate prediction. Here, we have developed FASTQINS, a pipeline for the detection of transposon insertions, and ANUBIS, a library of functions to evaluate and correct deviating factors known and uncharacterized until now. ANUBIS implements previously defined essentiality estimate models in addition to new approaches with advantages like not requiring a training set of genes to predict general essentiality. To highlight the applicability of these tools, and provide a set of recommendations on how to analyze transposon sequencing data, we performed a comprehensive study on artifacts corrections and essentiality estimation at a 1.5-bp resolution, in the genome-reduced bacterium Mycoplasma pneumoniae. We envision FASTQINS and ANUBIS to aid in the analysis of Tn-seq procedures and lead to the development of accurate genome essentiality estimates to guide applications such as designing live vaccines or growth optimization.

Highly efficient gene transfer in the mouse gut microbiota is enabled by the Incl2 conjugative plasmid TP114

The gut microbiota is a suspected hotspot for bacterial conjugation due to its high density and diversity of microorganisms. However, the contribution of different conjugative plasmid families to horizontal gene transfer in this environment remains poorly characterized. Here, we systematically quantified the transfer rates in the mouse intestinal tract for 13 conjugative plasmids encompassing 10 major incompatibility groups. The vast majority of these plasmids were unable to perform conjugation in situ or only reached relatively low transfer rates. Surprisingly, IncI₂ conjugative plasmid TP114 was identified as a proficient DNA delivery system in this environment, with the ability to transfer to virtually 100% of the probed recipient bacteria. We also show that a type IV pilus present in I-complex conjugative plasmids plays a crucial role for the transfer of TP114 in the mouse intestinal microbiota, most likely by contributing to mating pair stabilization. These results provide new insights on the mobility of genes in the gut microbiota and highlights TP114 as a very efficient DNA delivery system of interest for microbiome editing tools.

ATAC-Seq Identifies Chromatin Landscapes Linked to the Regulation of Oxidative Stress in the Human Fungal Pathogen Candida albicans

Human fungal pathogens often encounter fungicidal stress upon host invasion, but they can swiftly adapt by transcriptional reprogramming that enables pathogen survival. Fungal immune evasion is tightly connected to chromatin regulation. Hence, fungal chromatin modifiers pose alternative treatment options to combat fungal infections. Here, we present an assay for transposase-accessible chromatin using sequencing (ATAC-seq) protocol adapted for the opportunistic pathogen Candida albicans to gain further insight into the interplay of chromatin accessibility and gene expression mounted during fungal adaptation to oxidative stress. The ATAC-seq workflow not only facilitates the robust detection of genomic regions with accessible chromatin but also allows for the precise modeling of nucleosome positions in C. albicans. Importantly, the data reveal genes with altered chromatin accessibility in upstream regulatory regions, which correlate with transcriptional regulation during oxidative stress. Interestingly, many genes show increased chromatin accessibility without change in gene expression upon stress exposure. Such chromatin signatures could predict yet unknown regulatory factors under highly dynamic transcriptional control. Additionally, de novo motif analysis in genomic regions with increased chromatin accessibility upon H₂O₂ treatment shows significant enrichment for Cap1 binding sites, a major factor of oxidative stress responses in C. albicans. Taken together, the ATAC-seq workflow enables the identification of chromatin signatures and highlights the dynamics of regulatory mechanisms mediating environmental adaptation of C. albicans.

Transposase-assisted tagmentation of RNA/DNA hybrid duplexes

Tn5-mediated transposition of double-strand DNA has been widely utilized in various high-throughput sequencing applications. Here, we report that the Tn5 transposase is also capable of direct tagmentation of RNA/DNA hybrids in vitro. As a proof-of-concept application, we utilized this activity to replace the traditional library construction procedure of RNA sequencing, which contains many laborious and time-consuming processes. Results of Transposase-assisted RNA/DNA hybrids Co-tagmEntation (termed ‘TRACE-seq’) are compared to traditional RNA-seq methods in terms of detected gene number, gene body coverage, gene expression measurement, library complexity, and differential expression analysis. At the meantime, TRACE-seq enables a cost-effective one-tube library construction protocol and hence is more rapid (within 6 hr) and convenient. We expect this tagmentation activity on RNA/DNA hybrids to have broad potentials on RNA biology and chromatin research.

Exploring the fitness benefits of genome reduction in Escherichia coli by a selection-driven approach

Artificial simplification of bacterial genomes is thought to have the potential to yield cells with reduced complexity, enhanced genetic stability, and improved cellular economy. Of these goals, economical gains, supposedly due to the elimination of superfluous genetic material, and manifested in elevated growth parameters in selected niches, have not yet been convincingly achieved. This failure might stem from limitations of the targeted genome reduction approach that assumes full knowledge of gene functions and interactions, and allows only a limited number of reduction trajectories to interrogate. To explore the potential fitness benefits of genome reduction, we generated successive random deletions in E. coli by a novel, selection-driven, iterative streamlining process. The approach allows the exploration of multiple streamlining trajectories, and growth periods inherent in the procedure ensure selection of the fittest variants of the population. By generating single- and multiple-deletion strains and reconstructing the deletions in the parental genetic background, we showed that favourable deletions can be obtained and accumulated by the procedure. The most reduced multiple-deletion strain, obtained in five deletion cycles (2.5% genome reduction), outcompeted the wild-type, and showed elevated biomass yield. The spectrum of advantageous deletions, however, affecting only a few genomic regions, appears to be limited.

The fitness landscape of the African Salmonella Typhimurium ST313 strain D23580 reveals unique properties of the pBT1 plasmid

We have used a transposon insertion sequencing (TIS) approach to establish the fitness landscape of the African Salmonella enterica serovar Typhimurium ST313 strain D23580, to complement our previous comparative genomic and functional transcriptomic studies. We used a genome-wide transposon library with insertions every 10 nucleotides to identify genes required for survival and growth in vitro and during infection of murine macrophages. The analysis revealed genomic regions important for fitness under two in vitro growth conditions. Overall, 724 coding genes were required for optimal growth in LB medium, and 851 coding genes were required for growth in SPI-2-inducing minimal medium. These findings were consistent with the essentiality analyses of other S. Typhimurium ST19 and S. Typhi strains. The global mutagenesis approach also identified 60 sRNAs and 413 intergenic regions required for growth in at least one in vitro growth condition. By infecting murine macrophages with the transposon library, we identified 68 genes that were required for intra-macrophage replication but did not impact fitness in vitro. None of these genes were unique to S. Typhimurium D23580, consistent with a high conservation of gene function between S. Typhimurium ST313 and ST19 and suggesting that novel virulence factors are not involved in the interaction of strain D23580 with murine macrophages. We discovered that transposon insertions rarely occurred in many pBT1 plasmid-encoded genes (36), compared with genes carried by the pSLT-BT virulence plasmid and other bacterial plasmids. The key essential protein encoded by pBT1 is a cysteinyl-tRNA synthetase, and our enzymological analysis revealed that the plasmid-encoded CysRSpBT1 had a lower ability to charge tRNA than the chromosomally-encoded CysRSchr enzyme. The presence of aminoacyl-tRNA synthetases in plasmids from a range of Gram-negative and Gram-positive bacteria suggests that plasmid-encoded essential genes are more common than had been appreciated.

EGFP insertional mutagenesis reveals multiple FXR2P fibrillar states with differing ribosome association in neurons

RNA-binding proteins (RBPs) function in higher-order assemblages such as RNA granules to regulate RNA localization and translation. The Fragile X homolog FXR2P is an RBP essential for formation of neuronal Fragile X granules that associate with axonal mRNA and ribosomes in the intact brain. However, the FXR2P domains important for assemblage formation in a cellular system are unknown. Here we used an EGFP insertional mutagenesis approach to probe for FXR2P intrinsic features that influence its structural states. We tested 18 different in-frame FXR2P^EGFP fusions in neurons and found that the majority did not impact assemblage formation. However, EGFP insertion within a 23 amino acid region of the low complexity (LC) domain induced FXR2P^EGFP assembly into two distinct fibril states that were observed in isolation or in highly-ordered bundles. FXR2P^EGFP fibrils exhibited different developmental timelines, ultrastructures and ribosome associations. Formation of both fibril types was dependent on an intact RNA-binding domain. These results suggest that restricted regions of the LC domain, together with the RNA-binding domain, may be important for FXR2P structural state organization in neurons.

Summary: A mutagenesis study reveals that the higher-order structural states of the RBP FXR2P in neurons can be regulated by manipulation of the LC and RNA-binding domains.

Biochemical analysis of nucleosome targeting by Tn5 transposase

Tn5 transposase is a bacterial enzyme that integrates a DNA fragment into genomic DNA, and is used as a tool for detecting nucleosome-free regions of genomic DNA in eukaryotes. However, in chromatin, the DNA targeting by Tn5 transposase has remained unclear. In the present study, we reconstituted well-positioned 601 dinucleosomes, in which two nucleosomes are connected with a linker DNA, and studied the DNA integration sites in the dinucleosomes by Tn5 transposase in vitro. We found that Tn5 transposase preferentially targets near the entry-exit DNA regions within the nucleosome. Tn5 transposase minimally cleaved the dinucleosome without a linker DNA, indicating that the linker DNA between two nucleosomes is important for the Tn5 transposase activity. In the presence of a 30 base-pair linker DNA, Tn5 transposase targets the middle of the linker DNA, in addition to the entry-exit sites of the nucleosome. Intriguingly, this Tn5-targeting characteristic is conserved in a dinucleosome substrate with a different DNA sequence from the 601 sequence. Therefore, the Tn5-targeting preference in the nucleosomal templates reported here provides important information for the interpretation of Tn5 transposase-based genomics methods, such as ATAC-seq.

Mutagenesis of Vibrio fischeri and Other Marine Bacteria Using Hyperactive Mini-Tn5 Derivatives

Mutagenizing bacterial genomes with selectable transposon insertions is an effective approach for identifying the genes underlying important phenotypes. Specific bacteria may require different tools and methods for effective transposon mutagenesis, and here we describe methods to mutagenize Vibrio fischeri using an engineered mini-Tn5 transposon with synthetic "mosaic" transposon ends. The transposon is delivered by conjugation on a plasmid that cannot replicate in V. fischeri and that encodes a hyperactive transposase outside the transposon itself. The chromosomal location of insertions can be readily identified by cloning and/or PCR-based methods described here. Although developed in V. fischeri, these tools and methods have proven effective in some other bacteria as well.

Genome-wide analysis of chromatin accessibility using ATAC-seq

Programs of gene transcription are controlled by cis-acting DNA elements, including enhancers, silencers, and promoters. Local accessibility of chromatin has proven to be a highly informative structural feature for identifying such regulatory elements, which tend to be relatively open due to their interactions with proteins. Recently, ATAC-seq (assay for transposase-accessible chromatin using sequencing) has emerged as one of the most powerful approaches for genome-wide chromatin accessibility profiling. This method assesses DNA accessibility using hyperactive Tn5 transposase, which simultaneously cuts DNA and inserts sequencing adaptors, preferentially in regions of open chromatin. ATAC-seq is a relatively simple procedure which can be applied to only a few thousand cells. It is well-suited to developing embryos of sea urchins and other echinoderms, which are a prominent experimental model for understanding the genomic control of animal development. In this chapter, we present a protocol for applying ATAC-seq to embryonic cells of sea urchins.

Efficient transposon mutagenesis mediated by an IPTG-controlled conditional suicide plasmid

Background

Transposon mutagenesis is highly valuable for bacterial genetic and genomic studies. The transposons are usually delivered into host cells through conjugation or electroporation of a suicide plasmid. However, many bacterial species cannot be efficiently conjugated or transformed for transposon saturation mutagenesis. For this reason, temperature-sensitive (ts) plasmids have also been developed for transposon mutagenesis, but prolonged incubation at high temperatures to induce ts plasmid loss can be harmful to the hosts and lead to enrichment of mutants with adaptive genetic changes. In addition, the ts phenotype of a plasmid is often strain- or species-specific, as it may become non-ts or suicidal in different bacterial species.

Results

We have engineered several conditional suicide plasmids that have a broad host range and whose loss is IPTG-controlled. One construct, which has the highest stability in the absence of IPTG induction, was then used as a curable vector to deliver hyperactive miniTn5 transposons for insertional mutagenesis. Our analyses show that these new tools can be used for efficient and regulatable transposon mutagenesis in Escherichia coli, Acinetobacter baylyi and Pseudomonas aeruginosa. In P. aeruginosa PAO1, we have used this method to generate a Tn5 insertion library with an estimated diversity of ~ 10⁸, which is ~ 2 logs larger than the best transposon insertional library of PAO1 and related Pseudomonas strains previously reported.

Conclusion

We have developed a number of IPTG-controlled conditional suicide plasmids. By exploiting one of them for transposon delivery, a highly efficient and broadly useful mutagenesis system has been developed. As the assay condition is mild, we believe that our methodology will have broad applications in microbiology research.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1319-0) contains supplementary material, which is available to authorized users.

Transposase-Mediated Chromosomal Integration of Exogenous Genes in Acidithiobacillus ferrooxidans

The development of Acidithiobacillus ferrooxidans as a non-model host organism for synthetic biology is hampered by a lack of genetic tools and techniques. New plating and liquid-based selection methods were developed to improve the identification of transformed cell lines. Enabled by these methods, a hyperactive transposase was used to generate mutants with integrated genes for the expression of the superfolder green fluorescent protein (sfGFP) gene or a 2-keto decarboxylase (KDC) gene, which enabled the production and secretion of isobutyric acid (IBA). An inverse PCR method was used to identify the insertion sites of the KDC gene in several mutants, leading to the identification of a region on the chromosome that may be suitable for future genetic insertions. These results demonstrate that functional exogenous metabolic genes have been chromosomally integrated into A. ferrooxidans, and this advance will facilitate the future development of these cells for new biotechnology applications.IMPORTANCE Acidithiobacillus ferrooxidans is an iron- and sulfur-oxidizing chemolithoautotroph and is a key member of the microbial consortia used in industrial biomining applications. There is interest in exploiting these cells for other metal recovery applications as well as in developing them as unique nonmodel microbial cell factories. Plasmid-driven expression of exogenous genes has been reported, and homologous recombination has been used to knock out some gene expression. Here, new selection protocols facilitated the development of a transposition method for chromosomal integration of exogenous genes into A. ferrooxidans This greatly expands the available genetic toolbox, which will open the door to greater metabolic engineering efforts for these cells.

Transposons: the agents of antibiotic resistance in bacteria

Transposons are a group of mobile genetic elements that are defined as a DNA sequence. Transposons can jump into different places of the genome; for this reason, they are called jumping genes. However, some transposons are always kept at the insertion site in the genome. Most transposons are inactivated and as a result, cannot move. Transposons are divided into two main groups: retrotransposons (class І) and DNA transposons (class ІІ). Retrotransposons are often found in eukaryotes. DNA transposons can be found in both eukaryotes and prokaryotes. The bacterial transposons belong to the DNA transposons and the Tn family, which are usually the carrier of additional genes for antibiotic resistance. Transposons can transfer from a plasmid to other plasmids or from a DNA chromosome to plasmid and vice versa that cause the transmission of antibiotic resistance genes in bacteria. The treatment of bacterial infectious diseases is difficult because of existing antibiotic resistance that part of this antibiotic resistance is caused by transposons. Bacterial infectious diseases are responsible for the increasing rise in world mortality rate. In this review, transposons and their roles have been studied in bacterial antibiotic resistance, in detail.

Simplified ChIP-exo assays

ChIP-seq and ChIP-exo identify where proteins bind along any genome in vivo. Although ChIP-seq is widely adopted in academic research, it has inherently high noise. In contrast, ChIP-exo has relatively low noise and achieves near-base pair resolution. Consequently, and unlike other genomic assays, ChIP-exo provides structural information on genome-wide binding proteins. Construction of ChIP-exo libraries is technically difficult. Here we describe greatly simplified ChIP-exo methods, each with use-specific advantages. This is achieved through assay optimization and use of Tn5 tagmentation and/or single-stranded DNA ligation. Greater library yields, lower processing time, and lower costs are achieved. In comparing assays, we reveal substantial limitations in other ChIP-based assays. Importantly, the new ChIP-exo assays allow high-resolution detection of some protein-DNA interactions in organs and in as few as 27,000 cells. It is suitable for high-throughput parallelization. The simplicity of ChIP-exo now makes it a highly appropriate substitute for ChIP-seq, and for broader adoption.

While ChIP-exo is low noise and highly informative regarding genome-wide binding proteins, libraries are difficult to construct. Here the authors present a simplified ChIP-exo method for high-resolution detection of interactions.

Robustness encoded across essential and accessory replicons of the ecologically versatile bacterium Sinorhizobium meliloti

Bacterial genome evolution is characterized by gains, losses, and rearrangements of functional genetic segments. The extent to which large-scale genomic alterations influence genotype-phenotype relationships has not been investigated in a high-throughput manner. In the symbiotic soil bacterium Sinorhizobium meliloti, the genome is composed of a chromosome and two large extrachromosomal replicons (pSymA and pSymB, which together constitute 45% of the genome). Massively parallel transposon insertion sequencing (Tn-seq) was employed to evaluate the contributions of chromosomal genes to growth fitness in both the presence and absence of these extrachromosomal replicons. Ten percent of chromosomal genes from diverse functional categories are shown to genetically interact with pSymA and pSymB. These results demonstrate the pervasive robustness provided by the extrachromosomal replicons, which is further supported by constraint-based metabolic modeling. A comprehensive picture of core S. meliloti metabolism was generated through a Tn-seq-guided in silico metabolic network reconstruction, producing a core network encompassing 726 genes. This integrated approach facilitated functional assignments for previously uncharacterized genes, while also revealing that Tn-seq alone missed over a quarter of wild-type metabolism. This work highlights the many functional dependencies and epistatic relationships that may arise between bacterial replicons and across a genome, while also demonstrating how Tn-seq and metabolic modeling can be used together to yield insights not obtainable by either method alone.

S. meliloti, which has traditionally facilitated ground-breaking insights into symbiotic communication, is also emerging as an excellent model for studying the evolution of functional relationships between bacterial chromosomes and anciently acquired accessory replicons. Multi-replicon genome architecture is present in ~ 10% of presently sequenced bacterial genomes. The S. meliloti genome is composed of three circular replicons, two of which are dispensable even though they encompass nearly half of the protein-coding genes in this organism. The construction of strains lacking these replicons has enabled a straightforward, genome-wide analysis of interactions between the chromosome and the non-essential replicons, revealing extensive functional cooperation between these genomic components. This analysis enabled a substantial refinement of a metabolic network model for S. meliloti. The integration of massively parallel genotype-phenotype screening with in silico metabolic reconstruction has enhanced our understanding of metabolic network structure as it relates to genome evolution in S. meliloti, and exemplifies an approach that may be productively applied to other taxa. The combined experimental and computational approach employed here further provides unique insights into the pervasive genetic interactions that may exist within large bacterial genomes.

Chromatin Accessibility-Based Characterization of the Gene Regulatory Network Underlying Plasmodium falciparum Blood-Stage Development

Underlying the development of malaria parasites within erythrocytes and the resulting pathogenicity is a hardwired program that secures proper timing of gene transcription and production of functionally relevant proteins. How stage-specific gene expression is orchestrated in vivo remains unclear. Here, using the assay for transposase accessible chromatin sequencing (ATAC-seq), we identified ∼4,000 regulatory regions in P. falciparum intraerythrocytic stages. The vast majority of these sites are located within 2 kb upstream of transcribed genes and their chromatin accessibility pattern correlates positively with abundance of the respective mRNA transcript. Importantly, these regions are sufficient to drive stage-specific reporter gene expression and DNA motifs enriched in stage-specific sets of regulatory regions interact with members of the P. falciparum AP2 transcription factor family. Collectively, this study provides initial insights into the in vivo gene regulatory network of P. falciparum intraerythrocytic stages and should serve as a valuable resource for future studies.

Large-Scale Low-Cost NGS Library Preparation Using a Robust Tn5 Purification and Tagmentation Protocol

Efficient preparation of high-quality sequencing libraries that well represent the biological sample is a key step for using next-generation sequencing in research. Tn5 enables fast, robust, and highly efficient processing of limited input material while scaling to the parallel processing of hundreds of samples. Here, we present a robust Tn5 transposase purification strategy based on an N-terminal His₆-Sumo3 tag. We demonstrate that libraries prepared with our in-house Tn5 are of the same quality as those processed with a commercially available kit (Nextera XT), while they dramatically reduce the cost of large-scale experiments. We introduce improved purification strategies for two versions of the Tn5 enzyme. The first version carries the previously reported point mutations E54K and L372P, and stably produces libraries of constant fragment size distribution, even if the Tn5-to-input molecule ratio varies. The second Tn5 construct carries an additional point mutation (R27S) in the DNA-binding domain. This construct allows for adjustment of the fragment size distribution based on enzyme concentration during tagmentation, a feature that opens new opportunities for use of Tn5 in customized experimental designs. We demonstrate the versatility of our Tn5 enzymes in different experimental settings, including a novel single-cell polyadenylation site mapping protocol as well as ultralow input DNA sequencing.

Assay for Transposase Accessible Chromatin (ATAC-Seq) to Chart the Open Chromatin Landscape of Human Pancreatic Islets

The regulatory mechanisms that ensure an accurate control of gene transcription are central to cellular function, development and disease. Such mechanisms rely largely on noncoding regulatory sequences that allow the establishment and maintenance of cell identity and tissue-specific cellular functions.The study of chromatin structure and nucleosome positioning allowed revealing transcription factor accessible genomic sites with regulatory potential, facilitating the comprehension of tissue-specific cis-regulatory networks. Recently a new technique coupled with high-throughput sequencing named Assay for Transposase Accessible Chromatin (ATAC-seq) emerged as an efficient method to chart open chromatin genome wide. The application of such technique to different cell types allowed unmasking tissue-specific regulatory elements and characterizing cis-regulatory networks. Herein we describe the implementation of the ATAC-seq method to human pancreatic islets, a tissue playing a central role in the control of glucose metabolism.

Statistical analysis of genetic interactions in Tn-Seq data

Tn-Seq is an experimental method for probing the functions of genes through construction of complex random transposon insertion libraries and quantification of each mutant's abundance using next-generation sequencing. An important emerging application of Tn-Seq is for identifying genetic interactions, which involves comparing Tn mutant libraries generated in different genetic backgrounds (e.g. wild-type strain versus knockout strain). Several analytical methods have been proposed for analyzing Tn-Seq data to identify genetic interactions, including estimating relative fitness ratios and fitting a generalized linear model. However, these have limitations which necessitate an improved approach. We present a hierarchical Bayesian method for identifying genetic interactions through quantifying the statistical significance of changes in enrichment. The analysis involves a four-way comparison of insertion counts across datasets to identify transposon mutants that differentially affect bacterial fitness depending on genetic background. Our approach was applied to Tn-Seq libraries made in isogenic strains of Mycobacterium tuberculosis lacking three different genes of unknown function previously shown to be necessary for optimal fitness during infection. By analyzing the libraries subjected to selection in mice, we were able to distinguish several distinct classes of genetic interactions for each target gene that shed light on their functions and roles during infection.

Himar1 Transposon for Efficient Random Mutagenesis in Aggregatibacter actinomycetemcomitans

Aggregatibacter actinomycetemcomitans is the primary etiological agent of aggressive periodontal disease. Identification of novel virulence factors at the genome-wide level is hindered by lack of efficient genetic tools to perform mutagenesis in this organism. The Himar1 mariner transposon is known to yield a random distribution of insertions in an organism’s genome with requirement for only a TA dinucleotide target and is independent of host-specific factors. However, the utility of this system in A. actinomycetemcomitans is unknown. In this study, we found that Himar1 transposon mutagenesis occurs at a high frequency (×10^-4), and can be universally applied to wild-type A. actinomycetemcomitans strains of serotypes a, b, and c. The Himar1 transposon inserts were stably inherited in A. actinomycetemcomitans transconjugants in the absence of antibiotics. A library of 16,000 mutant colonies of A. actinomycetemcomitans was screened for reduced biofilm formation. Mutants with transposon inserts in genes encoding pilus, putative ion transporters, multidrug resistant proteins, transcription regulators and enzymes involved in the synthesis of extracellular polymeric substance, bacterial metabolism and stress response were discovered in this screen. Our results demonstrated the utility of the Himar1 mutagenesis system as a novel genetic tool for functional genomic analysis in A. actinomycetemcomitans.

Correcting nucleotide-specific biases in high-throughput sequencing data

Background

High-throughput sequence (HTS) data exhibit position-specific nucleotide biases that obscure the intended signal and reduce the effectiveness of these data for downstream analyses. These biases are particularly evident in HTS assays for identifying regulatory regions in DNA (DNase-seq, ChIP-seq, FAIRE-seq, ATAC-seq). Biases may result from many experiment-specific factors, including selectivity of DNA restriction enzymes and fragmentation method, as well as sequencing technology-specific factors, such as choice of adapters/primers and sample amplification methods.

Results

We present a novel method to detect and correct position-specific nucleotide biases in HTS short read data. Our method calculates read-specific weights based on aligned reads to correct the over- or underrepresentation of position-specific nucleotide subsequences, both within and adjacent to the aligned read, relative to a baseline calculated in assay-specific enriched regions. Using HTS data from a variety of ChIP-seq, DNase-seq, FAIRE-seq, and ATAC-seq experiments, we show that our weight-adjusted reads reduce the position-specific nucleotide imbalance across reads and improve the utility of these data for downstream analyses, including identification and characterization of open chromatin peaks and transcription-factor binding sites.

Conclusions

A general-purpose method to characterize and correct position-specific nucleotide sequence biases fills the need to recognize and deal with, in a systematic manner, binding-site preference for the growing number of HTS-based epigenetic assays. As the breadth and impact of these biases are better understood, the availability of a standard toolkit to correct them will be important.

Electronic supplementary material

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

Assessing Tn5 and Sleeping Beauty for transpositional transgenesis by cytoplasmic injection into bovine and ovine zygotes

Transgenic domestic animals represent an alternative to bioreactors for large-scale production of biopharmaceuticals and could also provide more accurate biomedical models than rodents. However, their generation remains inefficient. Recently, DNA transposons allowed improved transgenesis efficiencies in mice and pigs. In this work, Tn5 and Sleeping Beauty (SB) transposon systems were evaluated for transgenesis by simple cytoplasmic injection in livestock zygotes. In the case of Tn5, the transposome complex of transposon nucleic acid and Tn5 protein was injected. In the case of SB, the supercoiled plasmids encoding a transposon and the SB transposase were co-injected. In vitro produced bovine zygotes were used to establish the cytoplasmic injection conditions. The in vitro cultured blastocysts were evaluated for reporter gene expression and genotyped. Subsequently, both transposon systems were injected in seasonally available ovine zygotes, employing transposons carrying the recombinant human factor IX driven by the beta-lactoglobulin promoter. The Tn5 approach did not result in transgenic lambs. In contrast, the Sleeping Beauty injection resulted in 2 lambs (29%) carrying the transgene. Both animals exhibited cellular mosaicism of the transgene. The extraembryonic tissues (placenta or umbilical cord) of three additional animals were also transgenic. These results show that transpositional transgenesis by cytoplasmic injection of SB transposon components can be applied for the production of transgenic lambs of pharmaceutical interest.

Large-Scale Transposition Mutagenesis of Streptomyces coelicolor Identifies Hundreds of Genes Influencing Antibiotic Biosynthesis

Gram-positive Streptomyces bacteria produce thousands of bioactive secondary metabolites, including antibiotics. To systematically investigate genes affecting secondary metabolism, we developed a hyperactive transposase-based Tn5 transposition system and employed it to mutagenize the model species Streptomyces coelicolor, leading to the identification of 51,443 transposition insertions. These insertions were distributed randomly along the chromosome except for some preferred regions associated with relatively low GC content in the chromosomal core. The base composition of the insertion site and its flanking sequences compiled from the 51,443 insertions implied a 19-bp expanded target site surrounding the insertion site, with a slight nucleic acid base preference in some positions, suggesting a relative randomness of Tn5 transposition targeting in the high-GC Streptomyces genome. From the mutagenesis library, 724 mutants involving 365 genes had altered levels of production of the tripyrrole antibiotic undecylprodigiosin (RED), including 17 genes in the RED biosynthetic gene cluster. Genetic complementation revealed that most of the insertions (more than two-thirds) were responsible for the changed antibiotic production. Genes associated with branched-chain amino acid biosynthesis, DNA metabolism, and protein modification affected RED production, and genes involved in signaling, stress, and transcriptional regulation were overrepresented. Some insertions caused dramatic changes in RED production, identifying future targets for strain improvement.IMPORTANCE High-GC Gram-positive streptomycetes and related actinomycetes have provided more than 100 clinical drugs used as antibiotics, immunosuppressants, and antitumor drugs. Their genomes harbor biosynthetic genes for many more unknown compounds with potential as future drugs. Here we developed a useful genome-wide mutagenesis tool based on the transposon Tn5 for the study of secondary metabolism and its regulation. Using Streptomyces coelicolor as a model strain, we found that chromosomal insertion was relatively random, except at some hot spots, though there was evidence of a slightly preferred 19-bp target site. We then used prodiginine production as a model to systematically survey genes affecting antibiotic biosynthesis, providing a global view of antibiotic regulation. The analysis revealed 348 genes that modulate antibiotic production, among which more than half act to reduce production. These might be valuable targets in future investigations of regulatory mechanisms, for strain improvement, and for the activation of silent biosynthetic gene clusters.

Sequencing thousands of single-cell genomes with combinatorial indexing

Single-cell genome sequencing has proven valuable for the detection of somatic variation, particularly in the context of tumor evolution. Current technologies suffer from high library construction costs, which restrict the number of cells that can be assessed and thus impose limitations on the ability to measure heterogeneity within a tissue. Here, we present single-cell combinatorial indexed sequencing (SCI-seq) as a means of simultaneously generating thousands of low-pass single-cell libraries for detection of somatic copy-number variants. We constructed libraries for 16,698 single cells from a combination of cultured cell lines, primate frontal cortex tissue and two human adenocarcinomas, and obtained a detailed assessment of subclonal variation within a pancreatic tumor.

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.

Comprehensive Essentiality Analysis of the Mycobacterium tuberculosis Genome via Saturating Transposon Mutagenesis

For decades, identifying the regions of a bacterial chromosome that are necessary for viability has relied on mapping integration sites in libraries of random transposon mutants to find loci that are unable to sustain insertion. To date, these studies have analyzed subsaturated libraries, necessitating the application of statistical methods to estimate the likelihood that a gap in transposon coverage is the result of biological selection and not the stochasticity of insertion. As a result, the essentiality of many genomic features, particularly small ones, could not be reliably assessed. We sought to overcome this limitation by creating a completely saturated transposon library in Mycobacterium tuberculosis In assessing the composition of this highly saturated library by deep sequencing, we discovered that a previously unknown sequence bias of the Himar1 element rendered approximately 9% of potential TA dinucleotide insertion sites less permissible for insertion. We used a hidden Markov model of essentiality that accounted for this unanticipated bias, allowing us to confidently evaluate the essentiality of features that contained as few as 2 TA sites, including open reading frames (ORF), experimentally identified noncoding RNAs, methylation sites, and promoters. In addition, several essential regions that did not correspond to known features were identified, suggesting uncharacterized functions that are necessary for growth. This work provides an authoritative catalog of essential regions of the M. tuberculosis genome and a statistical framework for applying saturating mutagenesis to other bacteria.

IMPORTANCE

Sequencing of transposon-insertion mutant libraries has become a widely used tool for probing the functions of genes under various conditions. The Himar1 transposon is generally believed to insert with equal probabilities at all TA dinucleotides, and therefore its absence in a mutant library is taken to indicate biological selection against the corresponding mutant. Through sequencing of a saturated Himar1 library, we found evidence that TA dinucleotides are not equally permissive for insertion. The insertion bias was observed in multiple prokaryotes and influences the statistical interpretation of transposon insertion (TnSeq) data and characterization of essential genomic regions. Using these insights, we analyzed a fully saturated TnSeq library for M. tuberculosis, enabling us to generate a comprehensive catalog of in vitro essentiality, including ORFs smaller than those found in any previous study, small (noncoding) RNAs (sRNAs), promoters, and other genomic features.

Validation of picogram- and femtogram-input DNA libraries for microscale metagenomics

High-throughput sequencing libraries are typically limited by the requirement for nanograms to micrograms of input DNA. This bottleneck impedes the microscale analysis of ecosystems and the exploration of low biomass samples. Current methods for amplifying environmental DNA to bypass this bottleneck introduce considerable bias into metagenomic profiles. Here we describe and validate a simple modification of the Illumina Nextera XT DNA library preparation kit which allows creation of shotgun libraries from sub-nanogram amounts of input DNA. Community composition was reproducible down to 100 fg of input DNA based on analysis of a mock community comprising 54 phylogenetically diverse Bacteria and Archaea. The main technical issues with the low input libraries were a greater potential for contamination, limited DNA complexity which has a direct effect on assembly and binning, and an associated higher percentage of read duplicates. We recommend a lower limit of 1 pg (∼100-1,000 microbial cells) to ensure community composition fidelity, and the inclusion of negative controls to identify reagent-specific contaminants. Applying the approach to marine surface water, pronounced differences were observed between bacterial community profiles of microliter volume samples, which we attribute to biological variation. This result is consistent with expected microscale patchiness in marine communities. We thus envision that our benchmarked, slightly modified low input DNA protocol will be beneficial for microscale and low biomass metagenomics.

A bend, flip and trap mechanism for transposon integration

Cut-and-paste DNA transposons of the mariner/Tc1 family are useful tools for genome engineering and are inserted specifically at TA target sites. A crystal structure of the mariner transposase Mos1 (derived from Drosophila mauritiana), in complex with transposon ends covalently joined to target DNA, portrays the transposition machinery after DNA integration. It reveals severe distortion of target DNA and flipping of the target adenines into extra-helical positions. Fluorescence experiments confirm dynamic base flipping in solution. Transposase residues W159, R186, F187 and K190 stabilise the target DNA distortions and are required for efficient transposon integration and transposition in vitro. Transposase recognises the flipped target adenines via base-specific interactions with backbone atoms, offering a molecular basis for TA target sequence selection. Our results will provide a template for re-designing mariner/Tc1 transposases with modified target specificities.

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

The complete set of DNA in a cell is referred to as its genome. Most genomes contain short fragments of DNA called transposons that can jump from one place to another. Transposons carry sections of DNA with them when they move, which creates diversity and can influence the evolution of a species. Transposons are also being exploited to develop tools for biotechnology and medical applications. One family of transposons – the Mariner/Tc1 family – has proved particularly useful in these endeavours because it is widespread in nature and can jump around the genomes of a broad range of species, including mammals.

DNA transposons are cut out of their position and then pasted at a new site by an enzyme called transposase, which is encoded by some of the DNA within the transposon. DNA is made up of strings of molecules called bases and Mariner/Tc1-family transposons can only insert into a new position in the genome at sites that have a specific sequence of two bases. However, it was not known how this target sequence is chosen and how the transposon inserts into it.

Morris et al. have now used a technique called X-ray crystallography to build a three-dimensional model of a Mariner/Tc1-family transposon as it inserts into a new position. The model shows that, as the transposon is pasted into its new site, the surrounding DNA bends. This causes two DNA bases in the surrounding DNA to flip out from their normal position in the DNA molecule, which enables them to be recognised by the transposase. Further experiments showed that this base-flipping is dynamic, that is, the two bases continuously flip in and out of position. Furthermore, Morris et al. identified which parts of the transposase enzyme are required for the transposon to be efficiently pasted into the genome.

Together these findings may help researchers to alter the transposase so that it can insert the transposon into different locations in a genome. This will hopefully lead to new tools for biotechnology and medical applications.

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

The design and analysis of transposon insertion sequencing experiments

Transposon insertion sequencing (TIS) is a powerful approach that can be extensively applied to the genome-wide definition of loci that are required for bacterial growth under diverse conditions. However, experimental design choices and stochastic biological processes can heavily influence the results of TIS experiments and affect downstream statistical analysis. In this Opinion article, we discuss TIS experimental parameters and how these factors relate to the benefits and limitations of the various statistical frameworks that can be applied to the computational analysis of TIS data.

Illumina MiSeq sequencing disfavours a sequence motif in the GFP reporter gene

Green fluorescent protein (GFP) is one of the most used reporter genes. We have used next-generation sequencing (NGS) to analyse the genetic diversity of a recombinant influenza A virus that expresses GFP and found a remarkable coverage dip in the GFP coding sequence. This coverage dip was present when virus-derived RT-PCR product or the parental plasmid DNA was used as starting material for NGS and regardless of whether Nextera XT transposase or Covaris shearing was used for DNA fragmentation. Therefore, the sequence coverage dip in the GFP coding sequence was not the result of emerging GFP mutant viruses or a bias introduced by Nextera XT fragmentation. Instead, we found that the Illumina MiSeq sequencing method disfavours the 'CCCGCC' motif in the GFP coding sequence.

A zero-inflated Poisson model for insertion tolerance analysis of genes based on Tn-seq data

MOTIVATION

Transposon insertion sequencing (Tn-seq) is an emerging technology that combines transposon mutagenesis with next-generation sequencing technologies for the identification of genes related to bacterial survival. The resulting data from Tn-seq experiments consist of sequence reads mapped to millions of potential transposon insertion sites and a large portion of insertion sites have zero mapped reads. Novel statistical method for Tn-seq data analysis is needed to infer functions of genes on bacterial growth.

RESULTS

In this article, we propose a zero-inflated Poisson model for analyzing the Tn-seq data that are high-dimensional and with an excess of zeros. Maximum likelihood estimates of model parameters are obtained using an expectation-maximization (EM) algorithm, and pseudogenes are utilized to construct appropriate statistical tests for the transposon insertion tolerance of normal genes of interest. We propose a multiple testing procedure that categorizes genes into each of the three states, hypo-tolerant, tolerant and hyper-tolerant, while controlling the false discovery rate. We evaluate the proposed method with simulation studies and apply the proposed method to a real Tn-seq data from an experiment that studied the bacterial pathogen, Campylobacter jejuniAvailability and implementation: We provide R code for implementing our proposed method at http://github.com/ffliu/TnSeq A user's guide with example data analysis is also available there.

CONTACT

pliu@iastate.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Tn-seq explorer: a tool for analysis of high-throughput sequencing data of transposon mutant libraries

Tn-seq is a high throughput technique for analysis of transposon mutant libraries. Tn-seq Explorer was developed as a convenient and easy-to-use package of tools for exploration of the Tn-seq data. In a typical application, the user will have obtained a collection of sequence reads adjacent to transposon insertions in a reference genome. The reads are first aligned to the reference genome using one of the tools available for this task. Tn-seq Explorer reads the alignment and the gene annotation, and provides the user with a set of tools to investigate the data and identify possibly essential or advantageous genes as those that contain significantly low counts of transposon insertions. Emphasis is placed on providing flexibility in selecting parameters and methodology most appropriate for each particular dataset. Tn-seq Explorer is written in Java as a menu-driven, stand-alone application. It was tested on Windows, Mac OS, and Linux operating systems. The source code is distributed under the terms of GNU General Public License. The program and the source code are available for download at http://www.cmbl.uga.edu/downloads/programs/Tn_seq_Explorer/ and https://github.com/sina-cb/Tn-seqExplorer.

Analysis of the genetic diversity of influenza A viruses using next-generation DNA sequencing

BACKGROUND

Influenza viruses exist as a large group of closely related viral genomes, also called quasispecies. The composition of this influenza viral quasispecies can be determined by an accurate and sensitive sequencing technique and data analysis pipeline. We compared the suitability of two benchtop next-generation sequencers for whole genome influenza A quasispecies analysis: the Illumina MiSeq sequencing-by-synthesis and the Ion Torrent PGM semiconductor sequencing technique.

RESULTS

We first compared the accuracy and sensitivity of both sequencers using plasmid DNA and different ratios of wild type and mutant plasmid. Illumina MiSeq sequencing reads were one and a half times more accurate than those of the Ion Torrent PGM. The majority of sequencing errors were substitutions on the Illumina MiSeq and insertions and deletions, mostly in homopolymer regions, on the Ion Torrent PGM. To evaluate the suitability of the two techniques for determining the genome diversity of influenza A virus, we generated plasmid-derived PR8 virus and grew this virus in vitro. We also optimized an RT-PCR protocol to obtain uniform coverage of all eight genomic RNA segments. The sequencing reads obtained with both sequencers could successfully be assembled de novo into the segmented influenza virus genome. After mapping of the reads to the reference genome, we found that the detection limit for reliable recognition of variants in the viral genome required a frequency of 0.5% or higher. This threshold exceeds the background error rate resulting from the RT-PCR reaction and the sequencing method. Most of the variants in the PR8 virus genome were present in hemagglutinin, and these mutations were detected by both sequencers.

CONCLUSIONS

Our approach underlines the power and limitations of two commonly used next-generation sequencers for the analysis of influenza virus gene diversity. We conclude that the Illumina MiSeq platform is better suited for detecting variant sequences whereas the Ion Torrent PGM platform has a shorter turnaround time. The data analysis pipeline that we propose here will also help to standardize variant calling in small RNA genomes based on next-generation sequencing data.