Comprehensive and relaxed search for oligonucleotide signatures in hierarchically clustered sequence datasets.

MOTIVATION: PCR, hybridization, DNA sequencing and other important methods in molecular diagnostics rely on both sequence-specific and sequence group-specific oligonucleotide primers and probes. Their design depends on the identification of oligonucleotide signatures in whole genome or marker gene sequences. Although genome and gene databases are generally available and regularly updated, collections of valuable signatures are rare. Even for single requests, the search for signatures becomes computationally expensive when working with large collections of target (and non-target) sequences. Moreover, with growing dataset sizes, the chance of finding exact group-matching signatures decreases, necessitating the application of relaxed search methods. The resultant substantial increase in complexity is exacerbated by the dearth of algorithms able to solve these problems efficiently. RESULTS: We have developed CaSSiS, a fast and scalable method for computing comprehensive collections of sequence- and sequence group-specific oligonucleotide signatures from large sets of hierarchically clustered nucleic acid sequence data. Based on the ARB Positional Tree (PT-)Server and a newly developed BGRT data structure, CaSSiS not only determines sequence-specific signatures and perfect group-covering signatures for every node within the cluster (i.e. target groups), but also signatures with maximal group coverage (sensitivity) within a user-defined range of non-target hits (specificity) for groups lacking a perfect common signature. An upper limit of tolerated mismatches within the target group, as well as the minimum number of mismatches with non-target sequences, can be predefined. Test runs with one of the largest phylogenetic gene sequence datasets available indicate good runtime and memory performance, and in silico spot tests have shown the usefulness of the resulting signature sequences as blueprints for group-specific oligonucleotide probes. AVAILABILITY: Software and Supplementary Material are available at http://cassis.in.tum.de/.

PTPan--overcoming memory limitations in oligonucleotide string matching for primer/probe design.

MOTIVATION: Nucleic acid diagnostics has high demands for non-heuristic exact and approximate oligonucleotide string matching concerning in silico primer/probe design in huge nucleic acid sequence collections. Unfortunately, public sequence repositories grow much faster than computer hardware performance and main memory capacity do. This growth imposes severe problems on existing oligonucleotide primer/probe design applications necessitating new approaches based on space-efficient indexing structures. RESULTS: We developed PTPan (spoken Peter Pan, 'PT' is for Position Tree, the earlier name of suffix trees), a space-efficient indexing structure for approximate oligonucleotide string matching in nucleic acid sequence data. Based on suffix trees, it combines partitioning, truncation and a new suffix tree stream compression to deal with large amounts of aligned and unaligned data. PTPan operates efficiently in main memory and on secondary storage, balancing between memory consumption and runtime during construction and application. Based on PTPan, applications supporting similarity search and primer/probe design have been implemented, namely FindFamily, ProbeMatch and ProbeDesign. All three use a weighted Levenshtein distance metric for approximative queries to find and rate matches with indels as well as substitutions. We integrated PTPan in the worldwide used software package ARB to demonstrate usability and performance. Comparing PTPan and the original ARB index for the very large ssu-rRNA database SILVA, we recognized a shorter construction time, extended functionality and dramatically reduced memory requirements at the price of expanded, but very reasonable query times. PTPan enables indexing of huge nucleic acid sequence collections at reasonable application response times. Not being limited by main memory, PTPan constitutes a major advancement regarding rapid oligonucleotide string matching in primer/probe design now and in the future facing the enormous growth of molecular sequence data. AVAILABILITY: Supplementary Material, PTPan stand-alone library and ARB-PTPan binary on http://ptpan.lrr.in.tum.de/. CONTACT: meierh@in.tum.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Deciphering the evolution and metabolism of an anammox bacterium from a community genome.

Anaerobic ammonium oxidation (anammox) has become a main focus in oceanography and wastewater treatment. It is also the nitrogen cycle's major remaining biochemical enigma. Among its features, the occurrence of hydrazine as a free intermediate of catabolism, the biosynthesis of ladderane lipids and the role of cytoplasm differentiation are unique in biology. Here we use environmental genomics--the reconstruction of genomic data directly from the environment--to assemble the genome of the uncultured anammox bacterium Kuenenia stuttgartiensis from a complex bioreactor community. The genome data illuminate the evolutionary history of the Planctomycetes and allow us to expose the genetic blueprint of the organism's special properties. Most significantly, we identified candidate genes responsible for ladderane biosynthesis and biological hydrazine metabolism, and discovered unexpected metabolic versatility.

PhylArray: phylogenetic probe design algorithm for microarray.

MOTIVATION: Microbial diversity is still largely unknown in most environments, such as soils. In order to get access to this microbial 'black-box', the development of powerful tools such as microarrays are necessary. However, the reliability of this approach relies on probe efficiency, in particular sensitivity, specificity and explorative power, in order to obtain an image of the microbial communities that is close to reality. RESULTS: We propose a new probe design algorithm that is able to select microarray probes targeting SSU rRNA at any phylogenetic level. This original approach, implemented in a program called 'PhylArray', designs a combination of degenerate and non-degenerate probes for each target taxon. Comparative experimental evaluations indicate that probes designed with PhylArray yield a higher sensitivity and specificity than those designed by conventional approaches. Applying the combined PhyArray/GoArrays strategy helps to optimize the hybridization performance of short probes. Finally, hybridizations with environmental targets have shown that the use of the PhylArray strategy can draw attention to even previously unknown bacteria.

Evaluation of sequence alignments and oligonucleotide probes with respect to three-dimensional structure of ribosomal RNA using ARB software package.

BACKGROUND: Availability of high-resolution RNA crystal structures for the 30S and 50S ribosomal subunits and the subsequent validation of comparative secondary structure models have prompted the biologists to use three-dimensional structure of ribosomal RNA (rRNA) for evaluating sequence alignments of rRNA genes. Furthermore, the secondary and tertiary structural features of rRNA are highly useful and successfully employed in designing rRNA targeted oligonucleotide probes intended for in situ hybridization experiments. RNA3D, a program to combine sequence alignment information with three-dimensional structure of rRNA was developed. Integration into ARB software package, which is used extensively by the scientific community for phylogenetic analysis and molecular probe designing, has substantially extended the functionality of ARB software suite with 3D environment. RESULTS: Three-dimensional structure of rRNA is visualized in OpenGL 3D environment with the abilities to change the display and overlay information onto the molecule, dynamically. Phylogenetic information derived from the multiple sequence alignments can be overlaid onto the molecule structure in a real time. Superimposition of both statistical and non-statistical sequence associated information onto the rRNA 3D structure can be done using customizable color scheme, which is also applied to a textual sequence alignment for reference. Oligonucleotide probes designed by ARB probe design tools can be mapped onto the 3D structure along with the probe accessibility models for evaluation with respect to secondary and tertiary structural conformations of rRNA. CONCLUSION: Visualization of three-dimensional structure of rRNA in an intuitive display provides the biologists with the greater possibilities to carry out structure based phylogenetic analysis. Coupled with secondary structure models of rRNA, RNA3D program aids in validating the sequence alignments of rRNA genes and evaluating probe target sites. Superimposition of the information derived from the multiple sequence alignment onto the molecule dynamically allows the researchers to observe any sequence inherited characteristics (phylogenetic information) in real-time environment. The extended ARB software package is made freely available for the scientific community via http://www.arb-home.de.

ARB: a software environment for sequence data.

The ARB (from Latin arbor, tree) project was initiated almost 10 years ago. The ARB program package comprises a variety of directly interacting software tools for sequence database maintenance and analysis which are controlled by a common graphical user interface. Although it was initially designed for ribosomal RNA data, it can be used for any nucleic and amino acid sequence data as well. A central database contains processed (aligned) primary structure data. Any additional descriptive data can be stored in database fields assigned to the individual sequences or linked via local or worldwide networks. A phylogenetic tree visualized in the main window can be used for data access and visualization. The package comprises additional tools for data import and export, sequence alignment, primary and secondary structure editing, profile and filter calculation, phylogenetic analyses, specific hybridization probe design and evaluation and other components for data analysis. Currently, the package is used by numerous working groups worldwide.

Graphical representation of ribosomal RNA probe accessibility data using ARB software package.

BACKGROUND: Taxon specific hybridization probes in combination with a variety of commonly used hybridization formats nowadays are standard tools in microbial identification. A frequently applied technology, fluorescence in situ hybridization (FISH), besides single cell identification, allows the localization and functional studies of the microbial community composition. Careful in silico design and evaluation of potential oligonucleotide probe targets is therefore crucial for performing successful hybridization experiments. RESULTS: The PROBE Design tools of the ARB software package take into consideration several criteria such as number, position and quality of diagnostic sequence differences while designing oligonucleotide probes. Additionally, new visualization tools were developed to enable the user to easily examine further sequence associated criteria such as higher order structure, conservation, G+C content, transition-transversion profiles and in situ target accessibility patterns. The different types of sequence associated information (SAI) can be visualized by user defined background colors within the ARB primary and secondary structure editors as well as in the PROBE Match tool. CONCLUSION: Using this tool, in silico probe design and evaluation can be performed with respect to in situ probe accessibility data. The evaluation of proposed probe targets with respect to higher-order rRNA structure is of importance for successful design and performance of in situ hybridization experiments. The entire ARB software package along with the probe accessibility data is available from the ARB home page http://www.arb-home.de

Development and application of oligonucleotide probes for in situ detection of thermotolerant Campylobacter in chicken faecal and liver samples.

Based on Campylobacter 16S- and 23S-rRNA sequence data oligonucleotide probes specific for thermotolerant campylobacters and for members of the genus Campylobacter have been developed. The 16S-rRNA-targeted probe CAMP 653, recommended for a comprehensive detection of members of the genus Campylobacter, specifically detected all Campylobacter strains used in this study. Detection of thermotolerant species has been achieved by the 23S-rRNA-targeted probe CAJECO1427. Optimal hybridisation conditions have been derived for both probes from melting profiles of fluorescence-labelled probe-target hybrids recorded in fluorescence in situ hybridisation experiments (FISH). The FISH assay was evaluated both by spiking poultry faecal samples with Campylobacter jejuni and by detecting Campylobacter in naturally colonized chickens. C. jejuni was reliably detected at levels of 10(6) cfu/g faeces after a 3- h enrichment step in Blood Preston Selective broth. Low level contaminations (

Diversity of bacteria growing in natural mineral water after bottling.

Bacterial growth occurs in noncarbonated natural mineral waters a few days after filling and storage at room temperature, a phenomenon known for more than 40 years. Using the full-cycle rRNA approach, we monitored the development of the planktonic bacterial community in a noncarbonated natural mineral water after bottling. Seven 16S rRNA gene libraries, comprising 108 clones in total, were constructed from water samples taken at various days after bottling and from two different bottle sizes. Sequence analyses identified 11 operational taxonomic units (OTUs), all but one affiliated with the betaproteobacterial order Burkholderiales (6 OTUs) or the class Alphaproteobacteria (4 OTUs). Fluorescence in situ hybridization (FISH) was applied in combination with DAPI (4',6'-diamidino-2-phenylindole) staining, viability staining, and microscopic counting to quantitatively monitor changes in bacterial community composition. A growth curve similar to that of a bacterium grown in a batch culture was recorded. In contrast to the current perception that Gammaproteobacteria are the most important bacterial components of natural mineral water in bottles, Betaproteobacteria dominated the growing bacterial community and accounted for 80 to 98% of all bacteria detected by FISH in the late-exponential and stationary-growth phases. Using previously published and newly designed genus-specific probes, members of the betaproteobacterial genera Hydrogenophaga, Aquabacterium, and Polaromonas were found to constitute a significant proportion of the bacterial flora (21 to 86% of all bacteria detected by FISH). For the first time, key genera responsible for bacterial growth in a natural mineral water were identified by applying molecular cultivation-independent techniques.

AxML: a fast program for sequential and parallel phylogenetic tree calculations based on the maximum likelihood method.

Heuristics for the NP-complete problem of calculating the optimal phylogenetic tree for a set of aligned rRNA sequences based on the maximum likelihood method are computationally expensive. In most existing algorithms the tree evaluation and branch length optimization functions, calculating the likelihood value for each tree topology examined in the search space, account for the greatest part of overall computation time. This paper introduces AxML, a program derived from fastDNAml, incorporating a fast topology evaluation function. The algorithmic optimizations introduced, represent a general approach for accelerating this function and are applicable to both sequential and parallel phylogeny programs, irrespective of their search space strategy. Therefore, their integration into three existing phylogeny programs rendered encouraging results. Experimental results on conventional processor architectures show a global run time improvement of 35% up to 47% for the various test sets and program versions we used.

Oligonucleotide microarray for 16S rRNA gene-based detection of all recognized lineages of sulfate-reducing prokaryotes in the environment.

For cultivation-independent detection of sulfate-reducing prokaryotes (SRPs) an oligonucleotide microarray consisting of 132 16S rRNA gene-targeted oligonucleotide probes (18-mers) having hierarchical and parallel (identical) specificity for the detection of all known lineages of sulfate-reducing prokaryotes (SRP-PhyloChip) was designed and subsequently evaluated with 41 suitable pure cultures of SRPs. The applicability of SRP-PhyloChip for diversity screening of SRPs in environmental and clinical samples was tested by using samples from periodontal tooth pockets and from the chemocline of a hypersaline cyanobacterial mat from Solar Lake (Sinai, Egypt). Consistent with previous studies, SRP-PhyloChip indicated the occurrence of Desulfomicrobium spp. in the tooth pockets and the presence of Desulfonema- and Desulfomonile-like SRPs (together with other SRPs) in the chemocline of the mat. The SRP-PhyloChip results were confirmed by several DNA microarray-independent techniques, including specific PCR amplification, cloning, and sequencing of SRP 16S rRNA genes and the genes encoding the dissimilatory (bi)sulfite reductase (dsrAB).