Model-based identification of conditionally-essential genes from transposon-insertion sequencing data

The understanding of bacterial gene function has been greatly enhanced by recent advancements in the deep sequencing of microbial genomes. Transposon insertion sequencing methods combines next-generation sequencing techniques with transposon mutagenesis for the exploration of the essentiality of genes under different environmental conditions. We propose a model-based method that uses regularized negative binomial regression to estimate the change in transposon insertions attributable to gene-environment changes in this genetic interaction study without transformations or uniform normalization. An empirical Bayes model for estimating the local false discovery rate combines unique and total count information to test for genes that show a statistically significant change in transposon counts. When applied to RB-TnSeq (randomized barcode transposon sequencing) and Tn-seq (transposon sequencing) libraries made in strains of Caulobacter crescentus using both total and unique count data the model was able to identify a set of conditionally beneficial or conditionally detrimental genes for each target condition that shed light on their functions and roles during various stress conditions.     

Genome-Scale Metabolic Network Validation of Shewanella oneidensis Using Transposon Insertion Frequency Analysis

Transposon mutagenesis, in combination with parallel sequencing, is becoming a powerful tool for en-masse mutant analysis. A probability generating function was used to explain observed miniHimar transposon insertion patterns, and gene essentiality calls were made by transposon insertion frequency analysis (TIFA). TIFA incorporated the observed genome and sequence motif bias of the miniHimar transposon. The gene essentiality calls were compared to: 1) previous genome-wide direct gene-essentiality assignments; and, 2) flux balance analysis (FBA) predictions from an existing genome-scale metabolic model of Shewanella oneidensis MR-1. A three-way comparison between FBA, TIFA, and the direct essentiality calls was made to validate the TIFA approach. The refinement in the interpretation of observed transposon insertions demonstrated that genes without insertions are not necessarily essential, and that genes that contain insertions are not always nonessential. The TIFA calls were in reasonable agreement with direct essentiality calls for S. oneidensis, but agreed more closely with E. coli essentiality calls for orthologs. The TIFA gene essentiality calls were in good agreement with the MR-1 FBA essentiality predictions, and the agreement between TIFA and FBA predictions was substantially better than between the FBA and the direct gene essentiality predictions.     

A collection of sequenced and mapped Ds transposon insertion sites in Arabidopsis thaliana

Insertional mutagenesis is a powerful tool for generating knockout mutations that facilitate associating biological functions with as yet uncharacterized open reading frames (ORFs) identified by genomic sequencing or represented in EST databases. We have generated a collection of Dissociation (Ds) transposon lines with insertions on all 5 Arabidopsis chromosomes. Here we report the insertion sites in 260 independent single-transposon lines, derived from four different Ds donor sites. We amplified and determined the genomic sequence flanking each transposon, then mapped its insertion site by identity of the flanking sequences to the corresponding sequence in the Arabidopsis genome database. This constitutes the largest collection of sequence-mapped Ds insertion sites unbiased by selection against the donor site. Insertion site clusters have been identified around three of the four donor sites on chromosomes 1 and 5, as well as near the nucleolus organizers on chromosomes 2 and 4. The distribution of insertions between ORFs and intergenic sequences is roughly proportional to the ratio of genic to intergenic sequence. Within ORFs, insertions cluster near the translational start codon, although we have not detected insertion site selectivity at the nucleotide sequence level. A searchable database of insertion site sequences for the 260 transposon insertion sites is available at http://sgio2.biotec.psu.edu/sr. This and other collections of Arabidopsis lines with sequence-identified transposon insertion sites are a valuable genetic resource for functional genomics studies because the transposon location is precisely known, the transposon can be remobilized to generate revertants, and the Ds insertion can be used to initiate further local mutagenesis.     

A collection of sequenced and mapped Ds transposon insertion sites in Arabidopsis thaliana

Insertional mutagenesis is a powerful tool for generating knockout mutations that facilitate associating biological functions with as yet uncharacterized open reading frames (ORFs) identified by genomic sequencing or represented in EST databases. We have generated a collection of Dissociation(Ds) transposon lines with insertions on all 5 Arabidopsischromosomes. Here we report the insertion sites in 260 independent single-transposon lines, derived from four different Ds donor sites. We amplified and determined the genomic sequence flanking each transposon, then mapped its insertion site by identity of the flanking sequences to the corresponding sequence in the Arabidopsisgenome database. This constitutes the largest collection of sequence-mapped Ds insertion sites unbiased by selection against the donor site. Insertion site clusters have been identified around three of the four donor sites on chromosomes 1 and 5, as well as near the nucleolus organizers on chromosomes 2 and 4. The distribution of insertions between ORFs and intergenic sequences is roughly proportional to the ratio of genic to intergenic sequence. Within ORFs, insertions cluster near the translational start codon, although we have not detected insertion site selectivity at the nucleotide sequence level. A searchable database of insertion site sequences for the 260 transposon insertion sites is available at http://sgio2.biotec.psu.edu/sr. This and other collections of Arabidopsislines with sequence-identified transposon insertion sites are a valuable genetic resource for functional genomics studies because the transposon location is precisely known, the transposon can be remobilized to generate revertants, and the Ds insertion can be used to initiate further local mutagenesis.     

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.     

Clonal Expansion Analysis of Transposon Insertions by High-Throughput Sequencing Identifies Candidate Cancer Genes in a PiggyBac Mutagenesis Screen

Somatic transposon mutagenesis in mice is an efficient strategy to investigate the genetic mechanisms of tumorigenesis. The identification of tumor driving transposon insertions traditionally requires the generation of large tumor cohorts to obtain information about common insertion sites. Tumor driving insertions are also characterized by their clonal expansion in tumor tissue, a phenomenon that is facilitated by the slow and evolving transformation process of transposon mutagenesis. We describe here an improved approach for the detection of tumor driving insertions that assesses the clonal expansion of insertions by quantifying the relative proportion of sequence reads obtained in individual tumors. To this end, we have developed a protocol for insertion site sequencing that utilizes acoustic shearing of tumor DNA and Illumina sequencing. We analyzed various solid tumors generated by PiggyBac mutagenesis and for each tumor >10⁶ reads corresponding to >10⁴ insertion sites were obtained. In each tumor, 9 to 25 insertions stood out by their enriched sequence read frequencies when compared to frequencies obtained from tail DNA controls. These enriched insertions are potential clonally expanded tumor driving insertions, and thus identify candidate cancer genes. The candidate cancer genes of our study comprised many established cancer genes, but also novel candidate genes such as Mastermind-like1 (Mamld1) and Diacylglycerolkinase delta (Dgkd). We show that clonal expansion analysis by high-throughput sequencing is a robust approach for the identification of candidate cancer genes in insertional mutagenesis screens on the level of individual tumors.     

A benchmark and an algorithm for detecting germline transposon insertions and measuring de novo transposon insertion frequencies

Transposons are genomic parasites, and their new insertions can cause instability and spur the evolution of their host genomes. Rapid accumulation of short-read whole-genome sequencing data provides a great opportunity for studying new transposon insertions and their impacts on the host genome. Although many algorithms are available for detecting transposon insertions, the task remains challenging and existing tools are not designed for identifying de novo insertions. Here, we present a new benchmark fly dataset based on PacBio long-read sequencing and a new method TEMP2 for detecting germline insertions and measuring de novo ‘singleton’ insertion frequencies in eukaryotic genomes. TEMP2 achieves high sensitivity and precision for detecting germline insertions when compared with existing tools using both simulated data in fly and experimental data in fly and human. Furthermore, TEMP2 can accurately assess the frequencies of de novo transposon insertions even with high levels of chimeric reads in simulated datasets; such chimeric reads often occur during the construction of short-read sequencing libraries. By applying TEMP2 to published data on hybrid dysgenic flies inflicted by de-repressed P-elements, we confirmed the continuous new insertions of P-elements in dysgenic offspring before they regain piRNAs for P-element repression. TEMP2 is freely available at Github: https://github.com/weng-lab/TEMP2.     

Use of Illumina sequencing to identify transposon insertions underlying mutant phenotypes in high‐copy Mutator lines of maize

High‐copy transposons have been effectively exploited as mutagens in a variety of organisms. However, their utility for phenotype‐driven forward genetics has been hampered by the difficulty of identifying the specific insertions responsible for phenotypes of interest. We describe a new method that can substantially increase the throughput of linking a disrupted gene to a known phenotype in high‐copy Mutator (Mu) transposon lines in maize. The approach uses the Illumina platform to obtain sequences flanking Mu elements in pooled, bar‐coded DNA samples. Insertion sites are compared among individuals of suitable genotype to identify those that are linked to the mutation of interest. DNA is prepared for sequencing by mechanical shearing, adapter ligation, and selection of DNA fragments harboring Mu flanking sequences by hybridization to a biotinylated oligonucleotide corresponding to the Mu terminal inverted repeat. This method yields dense clusters of sequence reads that tile approximately 400 bp flanking each side of each heritable insertion. The utility of the approach is demonstrated by identifying the causal insertions in four genes whose disruption blocks chloroplast biogenesis at various steps: thylakoid protein targeting (cpSecE), chloroplast gene expression (polynucleotide phosphorylase and PTAC12), and prosthetic group attachment (HCF208/CCB2). This method adds to the tools available for phenotype‐driven Mu tagging in maize, and could be adapted for use with other high‐copy transposons. A by‐product of the approach is the identification of numerous heritable insertions that are unrelated to the targeted phenotype, which can contribute to community insertion resources.     

Genome scanning for conditionally essential genes in Salmonella enterica Serotype Typhimurium

As more whole-genome sequences become available, there is an increasing demand for high-throughput methods that link genes to phenotypes, facilitating discovery of new gene functions. In this study, we describe a new version of the Tn-seq method involving a modified EZ:Tn5 transposon for genome-wide and quantitative mapping of all insertions in a complex mutant library utilizing massively parallel Illumina sequencing. This Tn-seq method was applied to a genome-saturating Salmonella enterica serotype Typhimurium mutant library recovered from selection under 3 different in vitro growth conditions (diluted Luria-Bertani [LB] medium, LB medium plus bile acid, and LB medium at 42°C), mimicking some aspects of host stressors. We identified an overlapping set of 105 protein-coding genes in S. Typhimurium that are conditionally essential under at least one of the above selective conditions. Competition assays using 4 deletion mutants (pyrD, glnL, recD, and STM14_5307) confirmed the phenotypes predicted by Tn-seq data, validating the utility of this approach in discovering new gene functions. With continuously increasing sequencing capacity of next generation sequencing technologies, this robust Tn-seq method will aid in revealing unexplored genetic determinants and the underlying mechanisms of various biological processes in Salmonella and the other approximately 70 bacterial species for which EZ:Tn5 mutagenesis has been established.     

Rapid Transposon Liquid Enrichment Sequencing (TnLE-seq) for Gene Fitness Evaluation in Underdeveloped Bacterial Systems

Whole-genome fitness analysis in microbes that uses saturating transposon mutagenesis combined with massively parallel sequencing (Tn-seq) is providing a measure of the contribution of each gene to a given growth condition. With this technique, gene fitness profiles and essential genes are discovered by simultaneous analyses of whether the absence of each gene product alters the growth kinetics of the bacterium. Here we modify the standard Tn-seq procedure to simplify and shorten the process by including delivery of the transposon through conjugation and liquid culture enrichment of the mutant pool, creating transposon liquid enrichment sequencing (TnLE-seq). To illustrate the success of these modifications and the robustness of the procedure, analyses of gene fitness of two cultures of the strictly anaerobic bacterium Desulfovibrio vulgaris Hildenborough were performed, with growth on lactate as the electron donor and sulfate as the electron acceptor. These data demonstrate reproducibility and provide a base condition for analysis of fitness changes in deletion mutants and in various growth conditions. The procedural modifications will facilitate the application of this powerful genetic analysis to microbes lacking a facile genetic system. Pilot studies produced 2.5 × 10⁵ and 3.4 × 10⁵ unique insertion mutants in the anaerobe Desulfovibrio vulgaris Hildenborough grown under typical laboratory conditions in rich medium. These analyses provided two similar high-resolution maps of gene fitness across the genome, and the method was also applied to growth in minimal medium. These results were also compared to the coverage obtained with a ca. 13,000-member cataloged transposon library constructed by sequencing transposon insertion sites in individual mutants.     

ITIS,a bioinformatics tool for accurate identification of transposon insertion sites using next-generation sequencing data

Background

Transposable elements constitute an important part of the genome and are essential in adaptive mechanisms. Transposition events associated with phenotypic changes occur naturally or are induced in insertional mutant populations. Transposon mutagenesis results in multiple random insertions and recovery of most/all the insertions is critical for forward genetics study. Using genome next-generation sequencing data and appropriate bioinformatics tool, it is plausible to accurately identify transposon insertion sites, which could provide candidate causal mutations for desired phenotypes for further functional validation.

Results

We developed a novel bioinformatics tool, ITIS (Identification of Transposon Insertion Sites), for localizing transposon insertion sites within a genome. It takes next-generation genome re-sequencing data (NGS data), transposon sequence, and reference genome sequence as input, and generates a list of highly reliable candidate insertion sites as well as zygosity information of each insertion. Using a simulated dataset and a case study based on an insertional mutant line from Medicago truncatula, we showed that ITIS performed better in terms of sensitivity and specificity than other similar algorithms such as RelocaTE, RetroSeq, TEMP and TIF. With the case study data, we demonstrated the efficiency of ITIS by validating the presence and zygosity of predicted insertion sites of the Tnt1 transposon within a complex plant system, M. truncatula.

Conclusion

This study showed that ITIS is a robust and powerful tool for forward genetic studies in identifying transposable element insertions causing phenotypes. ITIS is suitable in various systems such as cell culture, bacteria, yeast, insect, mammal and plant.

Electronic supplementary material

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

MTID: a database of Sleeping Beauty transposon insertions in mice

The Sleeping Beauty (SB) transposon system provides the first random insertional mutagen available for germline genetic screens in mice. In preparation for a large scale project to create, map and manage up to 5000 SB insertions, we have developed the Mouse Transposon Insertion Database (MTID; http://mouse.ccgb.umn.edu/transposon/). Each insertion's genomic position, as well as the distance between the insertion and the nearest annotated gene, are determined by a sequence analysis pipeline. Users can search the database using a specified nucleotide or genetic map position to identify the nearest insertion. Mouse reports describe insertions carried, strain, genotype and dates of birth and death. Insertion reports describes chromosome, nucleotide and genetic map positions, as well as nearest gene data from Ensembl, NCBI and Celera. The flanking sequence used to map the insertion is also provided. Researchers will be able to identify insertions of interest and request mice or frozen sperm that carry the insertion.     

Identification of cold shock gene loci in Sinorhizobium meliloti by using a luxAB reporter transposon

Using a luxAB reporter transposon, seven mutants of Sinorhizobium meliloti were identified as containing insertions in four cold shock loci. LuxAB activity was strongly induced (25- to 160-fold) after a temperature shift from 30 to 15 degrees C. The transposon and flanking host DNA from each mutant was cloned, and the nucleic acid sequence of the insertion site was determined. Unexpectedly, five of the seven luxAB mutants contained transposon insertions in the 16S and 23S rRNA genes of two of the three rrn operons of S. meliloti. Directed insertion of luxAB genes into each of the three rrn operons revealed that all three operons were similarly affected by cold shock. Two other insertions were found to be located downstream of a homolog of the major Escherichia coli cold shock gene, cspA. Although the cold shock loci were highly induced in response to a shift to low temperature, none of the insertions resulted in a statistically significant decrease in growth rate at 15 degrees C.     

Identification of the Maize Gravitropism Gene lazy plant1 by a Transposon-Tagging Genome Resequencing Strategy

Since their initial discovery, transposons have been widely used as mutagens for forward and reverse genetic screens in a range of organisms. The problems of high copy number and sequence divergence among related transposons have often limited the efficiency at which tagged genes can be identified. A method was developed to identity the locations of Mutator (Mu) transposons in the Zea mays genome using a simple enrichment method combined with genome resequencing to identify transposon junction fragments. The sequencing library was prepared from genomic DNA by digesting with a restriction enzyme that cuts within a perfectly conserved motif of the Mu terminal inverted repeats (TIR). Paired-end reads containing Mu TIR sequences were computationally identified and chromosomal sequences flanking the transposon were mapped to the maize reference genome. This method has been used to identify Mu insertions in a number of alleles and to isolate the previously unidentified lazy plant1 (la1) gene. The la1 gene is required for the negatively gravitropic response of shoots and mutant plants lack the ability to sense gravity. Using bioinformatic and fluorescence microscopy approaches, we show that the la1 gene encodes a cell membrane and nuclear localized protein. Our Mu-Taq method is readily adaptable to identify the genomic locations of any insertion of a known sequence in any organism using any sequencing platform.     

Sleeping Beauty transposition in the mouse genome is associated with changes in DNA methylation at the site of insertion

The Sleeping Beauty (SB) transposon (Tn) system is a nonviral gene delivery tool that has widespread application for transfer of therapeutic genes into the mammalian genome. To determine its utility as a gene delivery system, it was important to assess the epigenetic modifications associated with SB insertion into the genome and potential inactivation of the transgene. This study investigated the DNA methylation pattern of an SB Tn as well as the flanking genomic region at insertion sites in the mouse genome. The ubiquitous ROSA26 promoter and an initial part of the eGFP coding sequence in the SB Tn exhibited high levels of CpG methylation in transgenic mouse lines, irrespective of the chromosomal loci of the insertion sites. In contrast, no detectable CpG methylation in the endogenous mouse ROSA26 counterpart was observed in the same animals. Furthermore, significant hypomethylation was detected in neighboring chromosomal sequences of two unique SB Tn insertions compared to wild-type patterns. Taken together, these results suggest that SB Tn insertions into the mouse genome can be discriminated by DNA methylation machinery from an identical endogenous DNA sequence and can profoundly alter the DNA methylation status of the transgene cargo as well as flanking host genomic regions.     

Testing the palindromic target site model for DNA transposon insertion using the Drosophila melanogaster P-element

Understanding the molecular mechanisms that influence transposable element target site preferences is a fundamental challenge in functional and evolutionary genomics. Large-scale transposon insertion projects provide excellent material to study target site preferences in the absence of confounding effects of post-insertion evolutionary change. Growing evidence from a wide variety of prokaryotes and eukaryotes indicates that DNA transposons recognize staggered-cut palindromic target site motifs (TSMs). Here, we use over 10 000 accurately mapped P-element insertions in the Drosophila melanogaster genome to test predictions of the staggered-cut palindromic target site model for DNA transposon insertion. We provide evidence that the P-element targets a 14-bp palindromic motif that can be identified at the primary sequence level, which predicts the local spacing, hotspots and strand orientation of P-element insertions. Intriguingly, we find that the although P-element destroys the complete 14-bp target site upon insertion, the terminal three nucleotides of the P-element inverted repeats complement and restore the original TSM, suggesting a mechanistic link between transposon target sites and their terminal inverted repeats. Finally, we discuss how the staggered-cut palindromic target site model can be used to assess the accuracy of genome mappings for annotated P-element insertions.     

A pathway for the evolution of the plasmid NTP16 involving the novel kanamycin resistance transposon Tn4352

The kanamycin resistance determinant of the drug resistance plasmid NTP16 has been characterized by DNA sequencing and has been shown to possess all of the structural features of a transposable element. It is made up of a 1040-bp central region encoding a protein identical to the aminoglycoside 3'-phosphotransferase of Tn903, flanked by direct repeats of an element identical to IS26. This novel transposon has been designated Tn4352. Analysis of the host sequences flanking the transposon reveal that they are derived from a Tn3-like element, and contain no 8 base pair target size duplications which are normally created by the insertion of IS26-like elements. Comparison to the Tn3 sequence shows that the flanking sequences are noncontiguous within Tn3, with the clear implication that NTP16 has evolved from a similar plasmid encoding only ampicillin resistance (presumably NTP1) by the insertion of Tn4352 into the Tn3-like element, followed by a substantial deletion. The sequence analysis suggests that the initial insertion was into the tnpR gene of the ampicillin transposon, followed by a deletion extending to a specific site within tnpA.