Dataset Information

Cluster oligonucleotide signatures for rapid identification by sequencing.

ABSTRACT: BACKGROUND:Oligonucleotide signatures (signatures) have been widely used for studying microbial diversity and function in wet-lab settings, but using them for accurate in silico identification of organisms from high-throughput sequencing (HTS) data is only a proof of concept. Existing signature design programs for sequence signatures (signatures matching exactly one sequence) or clade signatures (signatures matching every sequence in a phylogenetic clade) are not able to identify all possible polymorphic sites for sequences with high similarity and perform poorly when handling large genome sequencing datasets. RESULTS:We introduce cluster signatures: subsequences that match perfectly and exclusively any group of sequences in a data set. Cluster signatures provide complete recall for primer/probe design and increased discrimination between sequences beyond that of clade signatures. Using cluster signatures for in silico identification of HTS targets achieves good precision/recall and running time performance. This method has been implemented into an open source tool, the Automated Oligonucleotide Design Pipeline (adop), included in supplementary material and available at: https://bitbucket.org/wenchen_aafc/aodp_v2.0_release . CONCLUSIONS:Cluster signatures provide a rapid and universal analysis tool to identify all possible short diagnostic DNA markers and variants from any DNA sequencing dataset. They are particularly useful in discriminating genetic material from closely related organisms and in detecting deleterious mutations in highly or perfectly conserved genomic sites.

SUBMITTER: Zahariev M

PROVIDER: S-EPMC6284311 | biostudies-literature | 2018 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Cluster oligonucleotide signatures for rapid identification by sequencing.

Zahariev Manuel M Chen Wen W Visagie Cobus M CM Lévesque C André CA

BMC bioinformatics 20181029 1

<h4>Background</h4>Oligonucleotide signatures (signatures) have been widely used for studying microbial diversity and function in wet-lab settings, but using them for accurate in silico identification of organisms from high-throughput sequencing (HTS) data is only a proof of concept. Existing signature design programs for sequence signatures (signatures matching exactly one sequence) or clade signatures (signatures matching every sequence in a phylogenetic clade) are not able to identify all pos ...[more]

PMID: 30522439

Similar Datasets

Project description:Third-generation DNA sequencing has enabled sequencing of long, unamplified DNA fragments with minimal steps. Direct sequencing of ssDNA or RNA gives valuable insights like base-level modifications, phosphoramidite synthesis yield estimates and strand quality analysis, without the need to add the complimentary strand. Direct sequencing of single-stranded nucleic acid species is challenging as they are non-compatible to the double-stranded sequencing adapters used by manufacturers. The MinION platform from Oxford Nanopore Technologies performs sequencing by passing single-strands of DNA through a layer of biological nanopore sensors; although sequencing is performed on single-strands, the recommended template by the manufacturer is double-stranded. We have identified that the MinION platform can perform sequencing of short, single-strand oligonucleotides directly without amplification or second-strand synthesis by performing a single annealing step before library preparation. Short 5' phosphorylated oligos when annealed to an adapter sequence can be directly sequenced in the 5' to 3' direction via nanopores. Adapter sequences were designed to bind to the 5' end of the oligos and to leave a 3' adenosine overhang after binding to their target. The 3' adenosine overhang of the adapter and the terminal phosphate makes the 5' end of the oligo analogous to an end-prepared dsDNA, rendering it compatible with ligation-based library preparation for sequencing. An oligo-pool containing 42,000, 120 nt orthogonal sequences was phosphorylated and sequenced using this method and ~90% of these sequences were recovered with high accuracy using BLAST. In the nanopore raw data, we have identified that empty signals can be wrongly identified as a valid read by the MinION platform and sometimes multiple signals containing several strands can be fused into a single raw sequence file due to segmentation faults in the software. This direct oligonucleotide sequencing method enables novel applications in DNA data storage systems where short oligonucleotides are the primary information carriers.

Project description:Rapid molecular identification of carbapenemase genes in Gram-negative bacteria is crucial for infection control and prevention, surveillance and for epidemiological purposes. Furthermore, it may have a significant impact upon determining the appropriate initial treatment and greatly benefit for critically ill patients. A novel oligonucleotide microarray-based assay was developed to simultaneously detect genes encoding clinically important carbapenemases as well as selected extended (ESBL) and narrow spectrum (NSBL) beta-lactamases directly from clonal culture material within few hours. Additionally, a panel of species specific markers was included to identify Escherichia coli, Pseudomonas aeruginosa, Citrobacter freundii/braakii, Klebsiella pneumoniae and Acinetobacter baumannii. The assay was tested using a panel of 117 isolates collected from urinary, blood and stool samples. For these isolates, phenotypic identifications and susceptibility tests were available. An independent detection of carbapenemase, ESBL and NSBL genes was carried out by various external reference laboratories using PCR methods. In direct comparison, the microarray correctly identified 98.2% of the covered carbapenemase genes. This included blaVIM (13 out of 13), blaGIM (2/2), blaKPC (27/27), blaNDM (5/5), blaIMP-2/4/7/8/13/14/15/16/31 (10/10), blaOXA-23 (12/13), blaOXA-40-group (7/7), blaOXA-48-group (32/33), blaOXA-51 (1/1) and blaOXA-58 (1/1). Furthermore, the test correctly identified additional beta-lactamases [blaOXA-1 (16/16), blaOXA-2 (4/4), blaOXA-9 (33/33), OXA-10 (3/3), blaOXA-51 (25/25), blaOXA-58 (2/2), CTX-M1/M15 (17/17) and blaVIM (1/1)]. In direct comparison to phenotypical identification obtained by VITEK or MALDI-TOF systems, 114 of 117 (97.4%) isolates, including Acinetobacter baumannii (28/28), Enterobacter spec. (5/5), Escherichia coli (4/4), Klebsiella pneumoniae (62/63), Klebsiella oxytoca (0/2), Pseudomonas aeruginosa (12/12), Citrobacter freundii (1/1) and Citrobacter braakii (2/2), were correctly identified by a panel of species specific probes. This assay might be easily extended, adapted and transferred to point of care platforms enabling fast surveillance, rapid detection and appropriate early treatment of infections caused by multiresistant Gram-negative bacteria.

Project description:BACKGROUND:Gene signatures are important to represent the molecular changes in the disease genomes or the cells in specific conditions, and have been often used to separate samples into different groups for better research or clinical treatment. While many methods and applications have been available in literature, there still lack powerful ones that can take account of the complex data and detect the most informative signatures. METHODS:In this article, we present a new framework for identifying gene signatures using Pareto-optimal cluster size identification for RNA-seq data. We first performed pre-filtering steps and normalization, then utilized the empirical Bayes test in Limma package to identify the differentially expressed genes (DEGs). Next, we used a multi-objective optimization technique, "Multi-objective optimization for collecting cluster alternatives" (MOCCA in R package) on these DEGs to find Pareto-optimal cluster size, and then applied k-means clustering to the RNA-seq data based on the optimal cluster size. The best cluster was obtained through computing the average Spearman's Correlation Score among all the genes in pair-wise manner belonging to the module. The best cluster is treated as the signature for the respective disease or cellular condition. RESULTS:We applied our framework to a cervical cancer RNA-seq dataset, which included 253 squamous cell carcinoma (SCC) samples and 22 adenocarcinoma (ADENO) samples. We identified a total of 582 DEGs by Limma analysis of SCC versus ADENO samples. Among them, 260 are up-regulated genes and 322 are down-regulated genes. Using MOCCA, we obtained seven Pareto-optimal clusters. The best cluster has a total of 35 DEGs consisting of all-upregulated genes. For validation, we ran PAMR (prediction analysis for microarrays) classifier on the selected best cluster, and assessed the classification performance. Our evaluation, measured by sensitivity, specificity, precision, and accuracy, showed high confidence. CONCLUSIONS:Our framework identified a multi-objective based cluster that is treated as a signature that can classify the disease and control group of samples with higher classification performance (accuracy 0.935) for the corresponding disease. Our method is useful to find signature for any RNA-seq or microarray data.

Dataset Information

Cluster oligonucleotide signatures for rapid identification by sequencing.

Publications

Cluster oligonucleotide signatures for rapid identification by sequencing.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets