SUMMARY: Non-coding RNAs are often neglected during genome annotation due to their difficulty of detection relative to protein coding genes. FindNonCoding takes a pattern mining approach to capture the essential sequence motifs and hairpin loops representing a non-coding RNA family and quickly identify matches in genomes. FindNonCoding was designed for ease of use and accurately finds non-coding RNAs with a low false discovery rate. AVAILABILITY AND IMPLEMENTATION: FindNonCoding is implemented within the DECIPHER package (v2.19.3) for R (v4.1) available from Bioconductor. Pre-trained models of common non-coding RNA families are included for bacteria, archaea and eukarya. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Antibiotic resistance is an important public health problem. One potential solution is the development of synergistic antibiotic combinations, in which the combination is more effective than the component drugs. However, experimental progress in this direction is severely limited by the number of samples required to exhaustively test for synergy, which grows exponentially with the number of drugs combined. We introduce a new metric for antibiotic synergy, motivated by the popular Fractional Inhibitory Concentration Index and the Highest Single Agent model. We also propose a new experimental design that samples along all appropriately normalized diagonals in concentration space, and prove that this design identifies all synergies among a set of drugs while only sampling a small fraction of the possible combinations. We applied our method to screen two- through eight-way combinations of eight antibiotics at 10 concentrations each, which requires sampling only 2,560 unique combinations of antibiotic concentrations.

MOTIVATION: A core task of genomics is to identify the boundaries of protein coding genes, which may cover over 90% of a prokaryote's genome. Several programs are available for gene finding, yet it is currently unclear how well these programs perform and whether any offers superior accuracy. This is in part because there is no universal benchmark for gene finding and, therefore, most developers select their own benchmarking strategy. RESULTS: Here, we introduce AssessORF, a new approach for benchmarking prokaryotic gene predictions based on evidence from proteomics data and the evolutionary conservation of start and stop codons. We applied AssessORF to compare gene predictions offered by GenBank, GeneMarkS-2, Glimmer and Prodigal on genomes spanning the prokaryotic tree of life. Gene predictions were 88-95% in agreement with the available evidence, with Glimmer performing the worst but no clear winner. All programs were biased towards selecting start codons that were upstream of the actual start. Given these findings, there remains considerable room for improvement, especially in the detection of correct start sites. AVAILABILITY AND IMPLEMENTATION: AssessORF is available as an R package via the Bioconductor package repository. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

UNLABELLED: For numerous experimental applications, PCR primers must be designed to efficiently amplify a set of homologous DNA sequences while giving rise to amplicons with maximally diverse signatures. We developed DesignSignatures to automate the process of designing primers for high-resolution melting (HRM), fragment length polymorphism (FLP) and sequencing experiments. The program also finds the best restriction enzyme to further diversify HRM or FLP signatures. This enables efficient comparison across many experimental designs in order to maximize signature diversity. AVAILABILITY AND IMPLEMENTATION: DesignSignatures is accessible as a web tool at www.DECIPHER.cee.wisc.edu, or as part of the DECIPHER open source software package for R available from BioConductor. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BACKGROUND: Alignment of large and diverse sequence sets is a common task in biological investigations, yet there remains considerable room for improvement in alignment quality. Multiple sequence alignment programs tend to reach maximal accuracy when aligning only a few sequences, and then diminish steadily as more sequences are added. This drop in accuracy can be partly attributed to a build-up of error and ambiguity as more sequences are aligned. Most high-throughput sequence alignment algorithms do not use contextual information under the assumption that sites are independent. This study examines the extent to which local sequence context can be exploited to improve the quality of large multiple sequence alignments. RESULTS: Two predictors based on local sequence context were assessed: (i) single sequence secondary structure predictions, and (ii) modulation of gap costs according to the surrounding residues. The results indicate that context-based predictors have appreciable information content that can be utilized to create more accurate alignments. Furthermore, local context becomes more informative as the number of sequences increases, enabling more accurate protein alignments of large empirical benchmarks. These discoveries became the basis for DECIPHER, a new context-aware program for sequence alignment, which outperformed other programs on large sequence sets. CONCLUSIONS: Predicting secondary structure based on local sequence context is an efficient means of breaking the independence assumption in alignment. Since secondary structure is more conserved than primary sequence, it can be leveraged to improve the alignment of distantly related proteins. Moreover, secondary structure predictions increase in accuracy as more sequences are used in the prediction. This enables the scalable generation of large sequence alignments that maintain high accuracy even on diverse sequence sets. The DECIPHER R package and source code are freely available for download at DECIPHER.cee.wisc.edu and from the Bioconductor repository.