Statistical alignment based on fragment insertion and deletion models

Motivation: The topic of this paper is the estimation of alignments and mutation rates based on stochastic sequence-evolution models that allow insertions and deletions of subsequences ('fragments') and not just single bases. The model we propose is a variant of a model introduced by Thorne et al., (J. Mol. Evol., 34, 3-16, 1992). The computational tractability of the model depends on certain restrictions in the insertion/deletion process; possible effects we discuss. Results: The process of fragment insertion and deletion in the sequence-evolution model induces a hidden Markov structure at the level of alignments and thus makes possible efficient statistical alignment algorithms. As an example we apply a sampling procedure to assess the variability in alignment and mutation parameter estimates for HVR1 sequences of human and orangutan, improving results of previous work. Simulation studies give evidence that estimation methods based on the proposed model also give satisfactory results when applied to data for which the restrictions in the insertion/deletion process do not hold. Availability: The source code of the software for sampling alignments and mutation rates for a pair of DNA sequences according to the fragment insertion and deletion model is freely available from http://www.math.uni-frankfurt.de/~stoch/software/mcmcsalut under the terms of the GNU public license (GPL, 2000).     

ngLOC: software and web server for predicting protein subcellular localization in prokaryotes and eukaryotes

ABSTRACT: BACKGROUND: Understanding protein subcellular localization is a necessary component toward understanding the overall function of a protein. Numerous computational methods have been published over the past decade, with varying degrees of success. Despite the large number of published methods in this area, only a small fraction of them are available for researchers to use in their own studies. Of those that are available, many are limited by predicting only a small number of major organelles in the cell. Additionally, the majority of methods predict only a single location, even though it is known that a large fraction of the proteins in eukaryotic species shuttle between locations to carry out their function. FINDINGS: We present a software package and a web server for predicting subcellular localization of protein sequences based on the ngLOC method. ngLOC is an n-gram-based Bayesian classifier that predicts subcellular localization of proteins both in prokaryotes and eukaryotes. The overall prediction accuracy varies from 89.8% to 91.4% across species. This program can predict 11 distinct locations each in plant and animal species. ngLOC also predicts 4 and 5 distinct locations on gram-positive and gram-negative bacterial datasets, respectively. CONCLUSIONS: ngLOC is a generic method that can be trained by data from a variety of species or classes for predicting protein subcellular localization. The standalone software is freely available for academic use under GNU GPL, and the ngLOC web server is also accessible at http://ngloc.unmc.edu.     

REMA: A computer-based mapping tool for analysis of restriction sites in multiple DNA sequences

REMA is an interactive web-based program which predicts endonuclease cut sites in DNA sequences. It analyses multiple sequences simultaneously and predicts the number and size of fragments as well as provides restriction maps. The users can select single or paired combinations of all commercially available enzymes. Additionally, REMA permits prediction of multiple sequence terminal fragment sizes and suggests suitable restriction enzymes for maximally discriminatory results. REMA is an easy to use, web based program which will have a wide application in molecular biology research. Availability: REMA is written in Perl and is freely available for non-commercial use. Detailed information on installation can be obtained from Jan Szubert (jan.szubert@gmail.com) and the web based application is accessible on the internet at the URL http://www.macaulay.ac.uk/rema Contact: b.singh@macaulay.ac.uk.     

A simple method using PyrosequencingTM to identify de novo SNPs in pooled DNA samples

A practical way to reduce the cost of surveying single-nucleotide polymorphism (SNP) in a large number of individuals is to measure the allele frequencies in pooled DNA samples. Pyrosequencing(TM) has been frequently used for this application because signals generated by this approach are proportional to the amount of DNA templates. The Pyrosequencing(TM) pyrogram is determined by the dispensing order of dNTPs, which is usually designed based on the known SNPs to avoid asynchronistic extensions of heterozygous sequences. Therefore, utilizing the pyrogram signals to identify de novo SNPs in DNA pools has never been undertook. Here, in this study we developed an algorithm to address this issue. With the sequence and pyrogram of the wild-type allele known in advance, we could use the pyrogram obtained from the pooled DNA sample to predict the sequence of the unknown mutant allele (de novo SNP) and estimate its allele frequency. Both computational simulation and experimental Pyrosequencing(TM) test results suggested that our method performs well. The web interface of our method is available at http://life.nctu.edu.tw/～yslin/PSM/.     

Mass analysis peptide sequence prediction (MAPSP)

The software tool MAPSP allows the combinatorial prediction of novel short peptides such as hormones with common sequence features. In addition, it assists in de novo sequencing in general. The tool was designed for use in conjunction with the analytical identification method of mass spectrometry (MS) and it can considerably speed-up the analysis of unknowns. Availability: The web interface is freely available at http://mapsp.ifg.uni-muenster.de/     

Efficient Identification of Arabidopsis Knock-Out Mutants Using DNA-Arrays of Transposon Flanking Sequences

Abstract: Mutants obtained by insertional mutagenesis are widely used for determining gene-phenotype relationships. In Arabidopsis thaliana, several populations mutagenized either by T-DNA or transposon insertion are available for screening for knockout mutants in genes of interest. We have so far screened our Arabidopsis population mutagenized with the Zea mays transposon En-1/Spm for insertion mutations in 718 genes, using PCR on DNA pools. Although successful, this common approach is too time consuming for use in systematic screening of all 25 498 predicted genes of the Arabidopsis genome. We therefore investigated the use of DNA arrays for the direct identification of mutants in our population. All transposon-flanking regions from individual plants are amplified by PCR and subsequently spotted at high density onto nylon membranes. A single hybridization experiment with a gene-specific probe then allows one to identify candidate mutant plants. The efficiency of each separate step was determined and optimized. Screening of filters representing 2880 plants for insertions in 144 genes and subsequent investigation of some of the potential insertion mutants suggest that an overall screening efficiency of 50 % is attained.     

HaploSearch: a tool for haplotype-sequence two-way transformation

Comparison of available mitochondrial DNA data is some times hindered by the data presentation format. HaploSearch is a simple tool for transforming DNA sequences into haplotype data and vice versa, speeding up the manipulation of large datasets. Although designed for mitochondrial DNA, HaploSearch could be used with any kind of DNA type. HaploSearch program, detailed software instructions and example files are freely available on the web http://www.haplosite.com/haplosearch/.     

Structural features of 3C6/C7 interband chromatin organization in Drosophila melanogaster polytene chromosomes

Mechanisms of chromatin decompaction in interbands of Drosophila polytene chromosomes have been studied. Using the example of interband 3C6/C7 of the X chromosome, we investigated the ability of different DNA segments to form an interband in a new genetic environment. We applied site-specific FLP recombination between two transposons with FRT-sites that allows introducing the DNA fragments from the interband 3C6/C7 into pICon(dv) transposon located in cytologically well-characterized 84F region of chromosome 3 followed by electron microscopic analysis of changes in the region caused by insertion of the DNA fragments into the transposon. It was shown that the insertion of a 276-bp DNA fragment from the 3C6/C7 region into the pICon(dv) transposon leads to the formation of a new interband between two thin bands represented by the transposon material. This DNA fragment is the known minimal sequence that is necessary and sufficient for interband generation. In addition, the sequence containing three copies repeated in tandem of 0.9 kb DNA from the interband 3C6/C7, including the 276-bp fragment, were integrated in the transposon. The presence of introduced DNA fragments did not change the morphology of the resulting interband. It was shown that the sites of DNase I hypersensitivity were saved in transposon sequences introduced into a new genetic environment. The data obtained allow analysis to be started of specific factors (proteins, DNA motifs, etc.) that determine the formation of decompacted chromatin in a certain interband region and chromomeric organization of interphase chromosomes in Drosophila as a whole.     

MINS2: revisiting the molecular code for transmembrane-helix recognition by the Sec61 translocon

To be fully functional, membrane proteins should not only fold, but also get inserted into the membrane, which is mediated by the Sec61 translocon. Recent experimental studies have attempted to elucidate how the Sec61 translocon accomplishes this delicate task by measuring the translocon-mediated membrane insertion free energies of 357 systematically designed peptides. On the basis of this data set, we have developed MINS2, a novel sequence-based computational method for predicting the membrane insertion free energies of protein sequences. A benchmark analysis of MINS2 shows that MINS2 signi.cantly outperforms previously proposed methods. Importantly, the application of MINS2 to known membrane protein structures shows that a better prediction of membrane insertion free energies does not lead to a better prediction of transmembrane segments of polytopic membrane proteins. AVAILABILITY: A web server for MINS2 is publicly available at http://service.bioinformatik.uni-saarland.de/mins. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Efficient insertion mutagenesis strategy for bacterial genomes involving electroporation of in vitro-assembled DNA transposition complexes of bacteriophage mu

An efficient insertion mutagenesis strategy for bacterial genomes based on the phage Mu DNA transposition reaction was developed. Incubation of MuA transposase protein with artificial mini-Mu transposon DNA in the absence of divalent cations in vitro resulted in stable but inactive Mu DNA transposition complexes, or transpososomes. Following delivery into bacterial cells by electroporation, the complexes were activated for DNA transposition chemistry after encountering divalent metal ions within the cells. Mini-Mu transposons were integrated into bacterial chromosomes with efficiencies ranging from 10(4) to 10(6) CFU/microg of input transposon DNA in the four species tested, i.e., Escherichia coli, Salmonella enterica serovar Typhimurium, Erwinia carotovora, and Yersinia enterocolitica. Efficiency of integration was influenced mostly by the competence status of a given strain or batch of bacteria. An accurate 5-bp target site duplication flanking the transposon, a hallmark of Mu transposition, was generated upon mini-Mu integration into the genome, indicating that a genuine DNA transposition reaction was reproduced within the cells of the bacteria studied. This insertion mutagenesis strategy for microbial genomes may be applicable to a variety of organisms provided that a means to introduce DNA into their cells is available.     

In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction

We have developed a website, www.in-silico.com, which runs a software program that performs three basic tasks in completely sequenced bacterial genomes by in silico analysis: PCR amplification, amplified fragment length polymorphism (AFLP-PCR) and endonuclease restriction. For PCR, after selection of the genome and introduction of primers, fragment size, DNA sequence and corresponding open reading frame (ORF) identity of the resulting PCR product is computed. Plasmids of sequenced species may be included in the analysis. Theoretical AFLP-PCR analyzes similar parameters, and includes a suggestion tool providing a list of commercial restriction enzyme pairs yielding up to 50 amplicons in the selected genome. Endonuclease restriction analysis of complete genomes and plasmids calculates the number of restriction sites for endonucleases in a given genome. If the number of fragments is 50 or fewer, pulsed field gel electrophoresis image and restriction maps are illustrated. Other tools that have been included in this site are ORF search by name and DNA to protein translation as well as restriction digestion of user-defined DNA sequences. AVAILABILITY: This is a new molecular biology resource freely available over the Internet at http://www.in-silico.com     

Natural genetic variation caused by transposable elements in humans

Transposons and transposon-like repetitive elements collectively occupy 44% of the human genome sequence. In an effort to measure the levels of genetic variation that are caused by human transposons, we have developed a new method to broadly detect transposon insertion polymorphisms of all kinds in humans. We began by identifying 606,093 insertion and deletion (indel) polymorphisms in the genomes of diverse humans. We then screened these polymorphisms to detect indels that were caused by de novo transposon insertions. Our method was highly efficient and led to the identification of 605 nonredundant transposon insertion polymorphisms in 36 diverse humans. We estimate that this represents 25-35% of approximately 2075 common transposon polymorphisms in human populations. Because we identified all transposon insertion polymorphisms with a single method, we could evaluate the relative levels of variation that were caused by each transposon class. The average human in our study was estimated to harbor 1283 Alu insertion polymorphisms, 180 L1 polymorphisms, 56 SVA polymorphisms, and 17 polymorphisms related to other forms of mobilized DNA. Overall, our study provides significant steps toward (i) measuring the genetic variation that is caused by transposon insertions in humans and (ii) identifying the transposon copies that produce this variation.     

Efficient Insertion Mutagenesis Strategy for Bacterial Genomes Involving Electroporation of In Vitro-Assembled DNA Transposition Complexes of Bacteriophage Mu

An efficient insertion mutagenesis strategy for bacterial genomes based on the phage Mu DNA transposition reaction was developed. Incubation of MuA transposase protein with artificial mini-Mu transposon DNA in the absence of divalent cations in vitro resulted in stable but inactive Mu DNA transposition complexes, or transpososomes. Following delivery into bacterial cells by electroporation, the complexes were activated for DNA transposition chemistry after encountering divalent metal ions within the cells. Mini-Mu transposons were integrated into bacterial chromosomes with efficiencies ranging from 10⁴ to 10⁶ CFU/μg of input transposon DNA in the four species tested, i.e., Escherichia coli, Salmonella enterica serovar Typhimurium, Erwinia carotovora, and Yersinia enterocolitica. Efficiency of integration was influenced mostly by the competence status of a given strain or batch of bacteria. An accurate 5-bp target site duplication flanking the transposon, a hallmark of Mu transposition, was generated upon mini-Mu integration into the genome, indicating that a genuine DNA transposition reaction was reproduced within the cells of the bacteria studied. This insertion mutagenesis strategy for microbial genomes may be applicable to a variety of organisms provided that a means to introduce DNA into their cells is available.     

GeneTrack--a genomic data processing and visualization framework

MOTIVATION: High-throughput 'ChIP-chip' and 'ChIP-seq' methodologies generate sufficiently large data sets that analysis poses significant informatics challenges, particularly for research groups with modest computational support. To address this challenge, we devised a software platform for storing, analyzing and visualizing high resolution genome-wide binding data. GeneTrack automates several steps of a typical data processing pipeline, including smoothing and peak detection, and facilitates dissemination of the results via the web. Our software is freely available via the Google Project Hosting environment at http://genetrack.googlecode.com     

PILER: identification and classification of genomic repeats

SUMMARY: Repeated elements such as satellites and transposons are ubiquitous in eukaryotic genomes. De novo computational identification and classification of such elements is a challenging problem. Therefore, repeat annotation of sequenced genomes has historically largely relied on sequence similarity to hand-curated libraries of known repeat families. We present a new approach to de novo repeat annotation that exploits characteristic patterns of local alignments induced by certain classes of repeats. We describe PILER, a package of efficient search algorithms for identifying such patterns. Novel repeats found using PILER are reported for Homo sapiens, Arabidopsis thalania and Drosophila melanogaster. AVAILABILITY: The PILER software is freely available at http://www.drive5.com/piler.     

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

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

Detecting hierarchical structure in molecular characteristics of disease using transitive approximations of directed graphs

MOTIVATION: Molecular diagnostics aims at classifying diseases into clinically relevant sub-entities based on molecular characteristics. Typically, the entities are split into subgroups, which might contain several variants yielding a hierarchical model of the disease. Recent years have introduced a plethora of new molecular screening technologies to molecular diagnostics. As a result molecular profiles of patients became complex and the classification task more difficult. RESULTS: We present a novel tool for detecting hierarchical structure in binary datasets. We aim for identifying molecular characteristics, which are stochastically implying other characteristics. The final hierarchical structure is encoded in a directed transitive graph where nodes represent molecular characteristics and a directed edge from a node A to a node B denotes that almost all cases with characteristic B also display characteristic A. Naturally, these graphs need to be transitive. In the core of our modeling approach lies the problem of calculating good transitive approximations of given directed but not necessarily transitive graphs. By good transitive approximation we understand transitive graphs, which differ from the reference graph in only a small number of edges. It is known that the problem of finding optimal transitive approximation is NP-complete. Here we develop an efficient heuristic for generating good transitive approximations. We evaluate the computational efficiency of the algorithm in simulations, and demonstrate its use in the context of a large genome-wide study on mature aggressive lymphomas. AVAILABILITY: The software used in our analysis is freely available from http://compdiag.uni-regensburg.de/software/transApproxs.shtml.     

Characterization of three Agrobacterium tumefaciens avirulent mutants with chromosomal mutations that affect induction of vir genes.

Three Agrobacterium tumefaciens mutants with chromosomal mutations that affect bacterial virulence were isolated by transposon mutagenesis. Two of the mutants were avirulent on all hosts tested. The third mutant, Ivr-211, was a host range mutant which was avirulent on Bryophyllum diagremontiana, Nicotiana tabacum, N. debneyi, N. glauca, and Daucus carota but was virulent on Zinnia elegans and Lycopersicon esculentum (tomato). That the mutant phenotype was due to the transposon insertion was determined by cloning the DNA containing the transposon insertion and using the cloned DNA to replace the wild-type DNA in the parent bacterial strain by marker exchange. The transposon insertions in the three mutants mapped at three widely separated locations on the bacterial chromosome. The effects of the mutations on various steps in tumor formation were examined. All three mutants showed no alteration in binding to carrot cells. However, none of the mutants showed any induction of vir genes by acetosyringone under conditions in which the parent strain showed vir gene induction. When the mutant bacteria were examined for changes in surface components, it was found that all three of the mutants showed a similar alteration in lipopolysaccharide (LPS). LPS from the mutants was larger in size and more heavily saccharide substituted than LPS from the parent strain. Two of the mutants showed no detectable alteration in outer membrane and periplasmic space proteins. The third mutant, Ivr-225, was missing a 79-kDa surface peptide. The reason(s) for the failure of vir gene induction in these mutants and its relationship, if any, to the observed alteration in LPS are unknown.