Dataset Information

Impact of normalization methods on high-throughput screening data with high hit rates and drug testing with dose-response data.

ABSTRACT: Most data analysis tools for high-throughput screening (HTS) seek to uncover interesting hits for further analysis. They typically assume a low hit rate per plate. Hit rates can be dramatically higher in secondary screening, RNAi screening and in drug sensitivity testing using biologically active drugs. In particular, drug sensitivity testing on primary cells is often based on dose-response experiments, which pose a more stringent requirement for data quality and for intra- and inter-plate variation. Here, we compared common plate normalization and noise-reduction methods, including the B-score and the Loess a local polynomial fit method under high hit-rate scenarios of drug sensitivity testing. We generated simulated 384-well plate HTS datasets, each with 71 plates having a range of 20 (5%) to 160 (42%) hits per plate, with controls placed either at the edge of the plates or in a scattered configuration.We identified 20% (77/384) as the critical hit-rate after which the normalizations started to perform poorly. Results from real drug testing experiments supported this estimation. In particular, the B-score resulted in incorrect normalization of high hit-rate plates, leading to poor data quality, which could be attributed to its dependency on the median polish algorithm. We conclude that a combination of a scattered layout of controls per plate and normalization using a polynomial least squares fit method, such as Loess helps to reduce column, row and edge effects in HTS experiments with high hit-rates and is optimal for generating accurate dose-response curves.john.mpindi@helsinki.fi.Supplementary information: R code and Supplementary data are available at Bioinformatics online.

SUBMITTER: Mpindi JP

PROVIDER: S-EPMC4653387 | biostudies-other | 2015 Dec

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

Impact of normalization methods on high-throughput screening data with high hit rates and drug testing with dose-response data.

Mpindi John-Patrick JP Swapnil Potdar P Dmitrii Bychkov B Jani Saarela S Saeed Khalid K Wennerberg Krister K Aittokallio Tero T Östling Päivi P Kallioniemi Olli O

Bioinformatics (Oxford, England) 20150807 23

<h4>Motivation</h4>Most data analysis tools for high-throughput screening (HTS) seek to uncover interesting hits for further analysis. They typically assume a low hit rate per plate. Hit rates can be dramatically higher in secondary screening, RNAi screening and in drug sensitivity testing using biologically active drugs. In particular, drug sensitivity testing on primary cells is often based on dose-response experiments, which pose a more stringent requirement for data quality and for intra- an ...[more]

PMID: 26254433

Similar Datasets

Project description:Recently, rapid improvements in technology and decrease in sequencing costs have made RNA-Seq a widely used technique to quantify gene expression levels. Various normalization approaches have been proposed, owing to the importance of normalization in the analysis of RNA-Seq data. A comparison of recently proposed normalization methods is required to generate suitable guidelines for the selection of the most appropriate approach for future experiments.In this paper, we compared eight non-abundance (RC, UQ, Med, TMM, DESeq, Q, RPKM, and ERPKM) and two abundance estimation normalization methods (RSEM and Sailfish). The experiments were based on real Illumina high-throughput RNA-Seq of 35- and 76-nucleotide sequences produced in the MAQC project and simulation reads. Reads were mapped with human genome obtained from UCSC Genome Browser Database. For precise evaluation, we investigated Spearman correlation between the normalization results from RNA-Seq and MAQC qRT-PCR values for 996 genes. Based on this work, we showed that out of the eight non-abundance estimation normalization methods, RC, UQ, Med, TMM, DESeq, and Q gave similar normalization results for all data sets. For RNA-Seq of a 35-nucleotide sequence, RPKM showed the highest correlation results, but for RNA-Seq of a 76-nucleotide sequence, least correlation was observed than the other methods. ERPKM did not improve results than RPKM. Between two abundance estimation normalization methods, for RNA-Seq of a 35-nucleotide sequence, higher correlation was obtained with Sailfish than that with RSEM, which was better than without using abundance estimation methods. However, for RNA-Seq of a 76-nucleotide sequence, the results achieved by RSEM were similar to without applying abundance estimation methods, and were much better than with Sailfish. Furthermore, we found that adding a poly-A tail increased alignment numbers, but did not improve normalization results.Spearman correlation analysis revealed that RC, UQ, Med, TMM, DESeq, and Q did not noticeably improve gene expression normalization, regardless of read length. Other normalization methods were more efficient when alignment accuracy was low; Sailfish with RPKM gave the best normalization results. When alignment accuracy was high, RC was sufficient for gene expression calculation. And we suggest ignoring poly-A tail during differential gene expression analysis.

Project description:BACKGROUND: High throughput methodologies such as microarrays, mass spectrometry and plate-based small molecule screens are increasingly used to facilitate discoveries from gene function to drug candidate identification. These large-scale experiments are typically carried out over the course of months and years, often without the controls needed to compare directly across the dataset. Few methods are available to facilitate comparisons of high throughput metabolic data generated in batches where explicit in-group controls for normalization are lacking. RESULTS: Here we describe MIPHENO (Mutant Identification by Probabilistic High throughput-Enabled Normalization), an approach for post-hoc normalization of quantitative first-pass screening data in the absence of explicit in-group controls. This approach includes a quality control step and facilitates cross-experiment comparisons that decrease the false non-discovery rates, while maintaining the high accuracy needed to limit false positives in first-pass screening. Results from simulation show an improvement in both accuracy and false non-discovery rate over a range of population parameters (p < 2.2 × 10(-16)) and a modest but significant (p < 2.2 × 10(-16)) improvement in area under the receiver operator characteristic curve of 0.955 for MIPHENO vs 0.923 for a group-based statistic (z-score). Analysis of the high throughput phenotypic data from the Arabidopsis Chloroplast 2010 Project (http://www.plastid.msu.edu/) showed ~ 4-fold increase in the ability to detect previously described or expected phenotypes over the group based statistic. CONCLUSIONS: Results demonstrate MIPHENO offers substantial benefit in improving the ability to detect putative mutant phenotypes from post-hoc analysis of large data sets. Additionally, it facilitates data interpretation and permits cross-dataset comparison where group-based controls are missing. MIPHENO is applicable to a wide range of high throughput screenings and the code is freely available as Additional file 1 as well as through an R package in CRAN.

Project description:The prevalence of lung disease due to infections with nontuberculous mycobacteria (NTM) has been increasing and surpassed tuberculosis (TB) in some countries. Treatment outcomes are often unsatisfactory, highlighting an urgent need for new anti-NTM medications. Although NTM in general do not respond well to TB specific drugs, the similarities between NTM and Mycobacterium tuberculosis at the molecular and cell structural level suggest that compound libraries active against TB could be leveraged for NTM drug discovery. Here we tested this hypothesis. The Pathogen Box from the Medicines for Malaria Venture (MMV) is a collection of 400 diverse drug-like compounds, among which 129 are known to be active against M. tuberculosis. By screening this compound collection against two NTM species, Mycobacterium abscessus and Mycobacterium avium, we showed that indeed the hit rates for NTM among TB active compounds is significantly higher compared to compounds that are not active against TB. MIC/dose response confirmation identified 10 top hits. Bactericidal activity determination demonstrated attractive potency for a subset of the confirmed hits. In vivo pharmacokinetic profiling showed that some of the compounds present reasonable starting points for medicinal chemistry programs. Three of the top hits were oxazolidinones, suggesting the potential for repositioning this class of protein synthesis inhibitors to replace linezolid which suffers from low potency. Two hits were inhibitors of the trehalose monomycolate transporter MmpL3, suggesting that this transmembrane protein may be an attractive target for NTM. Other hits are predicted to target a range of functions, including cell division (FtsZ), DNA gyrase (GyrB), dihydrofolate reductase, RNA polymerase and ABC transporters. In conclusion, our study showed that screening TB active compounds for activity against NTM resulted in high hit rates, suggesting that this may be an attractive approach to kick start NTM drug discovery projects. In addition, the work identified a series of novel high value NTM hits with associated candidate targets which can be followed up in hit-to-lead projects for the discovery of new NTM antibiotics.

Dataset Information

Impact of normalization methods on high-throughput screening data with high hit rates and drug testing with dose-response data.

Publications

Impact of normalization methods on high-throughput screening data with high hit rates and drug testing with dose-response data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure