Dataset Information

Bayesian analysis of gene essentiality based on sequencing of transposon insertion libraries.

ABSTRACT:

Motivation

Next-generation sequencing affords an efficient analysis of transposon insertion libraries, which can be used to identify essential genes in bacteria. To analyse this high-resolution data, we present a formal Bayesian framework for estimating the posterior probability of essentiality for each gene, using the extreme-value distribution to characterize the statistical significance of the longest region lacking insertions within a gene. We describe a sampling procedure based on the Metropolis-Hastings algorithm to calculate posterior probabilities of essentiality while simultaneously integrating over unknown internal parameters.

Results

Using a sequence dataset from a transposon library for Mycobacterium tuberculosis, we show that this Bayesian approach predicts essential genes that correspond well with genes shown to be essential in previous studies. Furthermore, we show that by using the extreme-value distribution to characterize genomic regions lacking transposon insertions, this method is capable of identifying essential domains within genes. This approach can be used for analysing transposon libraries in other organisms and augmenting essentiality predictions with statistical confidence scores.

SUBMITTER: DeJesus MA

PROVIDER: S-EPMC3597147 | biostudies-literature | 2013 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bayesian analysis of gene essentiality based on sequencing of transposon insertion libraries.

DeJesus Michael A MA Zhang Yanjia J YJ Sassetti Christopher M CM Rubin Eric J EJ Sacchettini James C JC Ioerger Thomas R TR

Bioinformatics (Oxford, England) 20130129 6

<h4>Motivation</h4>Next-generation sequencing affords an efficient analysis of transposon insertion libraries, which can be used to identify essential genes in bacteria. To analyse this high-resolution data, we present a formal Bayesian framework for estimating the posterior probability of essentiality for each gene, using the extreme-value distribution to characterize the statistical significance of the longest region lacking insertions within a gene. We describe a sampling procedure based on t ...[more]

PMID: 23361328

Similar Datasets

Project description:The massive sequencing of transposon insertion mutant libraries (Tn-Seq) represents a commonly used method to determine essential genes in bacteria. Using a hypersaturated transposon mutant library consisting of 400,096 unique Tn insertions, 523 genes were classified as essential in Escherichia coli K-12 MG1655. This provided a useful genome-wide gene essentiality landscape for rapidly identifying 233 of 301 essential genes previously validated by a knockout study. However, there was a discrepancy in essential gene sets determined by conventional gene deletion methods and Tn-Seq, although different Tn-Seq studies reported different extents of discrepancy. We have elucidated two causes of this discrepancy. First, 68 essential genes not detected by Tn-Seq contain nonessential subgenic domains that are tolerant to transposon insertion, which leads to the false assignment of an essential gene as a nonessential or dispensable gene. These genes exhibited a high level of transposon insertion in their subgenic nonessential domains. In contrast, 290 genes were additionally categorized as essential by Tn-Seq, although their knockout mutants were available. The comparative analysis of Tn-Seq and high-resolution footprinting of nucleoid-associated proteins (NAPs) revealed that a protein-DNA interaction hinders transposon insertion. We identified 213 false-positive genes caused by NAP-genome interactions. These two limitations have to be considered when addressing essential bacterial genes using Tn-Seq. Furthermore, a comparative analysis of high-resolution Tn-Seq with other data sets is required for a more accurate determination of essential genes in bacteria. IMPORTANCE Transposon mutagenesis is an efficient way to explore gene essentiality of a bacterial genome. However, there was a discrepancy between the essential gene set determined by transposon mutagenesis and that determined using single-gene knockout strains. In this study, we generated a hypersaturated Escherichia coli transposon mutant library comprising approximately 400,000 different mutants. Determination of transposon insertion sites using next-generation sequencing provided a high-resolution essentiality landscape of the E. coli genome. We identified false negatives of essential gene discovery due to the permissive insertion of transposons in the C-terminal region. Comparisons between the transposon insertion landscape with binding profiles of DNA-binding proteins revealed interference of nucleoid-associated proteins to transposon insertion, generating false positives of essential gene discovery. Consideration of these findings is required to avoid the misinterpretation of transposon mutagenesis results.

Project description:Mycobacterium abscessus is an emerging opportunistic human pathogen that naturally resists most major classes of antibiotics, making infections difficult to treat. Thus far, little is known about M. abscessus physiology, pathogenesis, and drug resistance. Genome-wide analyses have comprehensively catalogued genes with essential functions in Mycobacterium tuberculosis and Mycobacterium avium subsp. hominissuis (here, M. avium) but not in M. abscessus. By optimizing transduction conditions, we achieved full saturation of TA insertion sites with Himar1 transposon mutagenesis in the M. abscessus ATCC 19977T genome, as confirmed by deep sequencing prior to essentiality analyses of annotated genes and other genomic features. The overall densities of inserted TA sites (85.7%), unoccupied TA sites (14.3%), and nonpermissive TA sites (8.1%) were similar to results in M. tuberculosis and M. avium. Of the 4,920 annotated genes, 326 were identified as essential, 269 (83%) of which have mutual homology with essential M. tuberculosis genes, while 39 (12%) are homologous to genes that are not essential in M. tuberculosis and M. avium, and 11 (3.4%) only have homologs in M. avium. Interestingly, 7 (2.1%) essential M. abscessus genes have no homologs in either M. tuberculosis or M. avium, two of which were found in phage-like elements. Most essential genes are involved in DNA replication, RNA transcription and translation, and posttranslational events to synthesize important macromolecules. Some essential genes may be involved in M. abscessus pathogenesis and antibiotics response, including certain essential tRNAs and new short open reading frames. Our findings will help to pave the way for better understanding of M. abscessus and benefit development of novel bactericidal drugs against M. abscessus. IMPORTANCE Limited knowledge regarding Mycobacterium abscessus pathogenesis and intrinsic resistance to most classes of antibiotics is a major obstacle to developing more effective strategies to prevent and mitigate disease. Using optimized procedures for Himar1 transposon mutagenesis and deep sequencing, we performed a comprehensive analysis to identify M. abscessus genetic elements essential for in vitro growth and compare them to similar data sets for M. tuberculosis and M. avium subsp. hominissuis. Most essential M. abscessus genes have mutual homology with essential M. tuberculosis genes, providing a foundation for leveraging available knowledge from M. tuberculosis to develop more effective drugs and other interventions against M. abscessus. A small number of essential genes unique to M. abscessus deserve further attention to gain insights into what makes M. abscessus different from other mycobacteria. The essential genes and other genomic features such as short open reading frames and noncoding RNA identified here will provide useful information for future study of M. abscessus pathogenicity and new drug development.

Dataset Information

Bayesian analysis of gene essentiality based on sequencing of transposon insertion libraries.

Motivation

Results

Publications

Bayesian analysis of gene essentiality based on sequencing of transposon insertion libraries.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets