Automatic Peak Selection by a Benjamini-Hochberg-Based Algorithm

Abstract

A common issue in bioinformatics is that computational methods often generate a large number of predictions sorted according to certain confidence scores. A key problem is then determining how many predictions must be selected to include most of the true predictions while maintaining reasonably high precision. In nuclear magnetic resonance (NMR)-based protein structure determination, for instance, computational peak picking methods are becoming more and more common, although expert-knowledge remains the method of choice to determine how many peaks among thousands of candidate peaks should be taken into consideration to capture the true peaks. Here, we propose a Benjamini-Hochberg (B-H)-based approach that automatically selects the number of peaks. We formulate the peak selection problem as a multiple testing problem. Given a candidate peak list sorted by either volumes or intensities, we first convert the peaks into -values and then apply the B-H-based algorithm to automatically select the number of peaks. The proposed approach is tested on the state-of-the-art peak picking methods, including WaVPeak [1] and PICKY [2]. Compared with the traditional fixed number-based approach, our approach returns significantly more true peaks. For instance, by combining WaVPeak or PICKY with the proposed method, the missing peak rates are on average reduced by 20% and 26%, respectively, in a benchmark set of 32 spectra extracted from eight proteins. The consensus of the B-H-selected peaks from both WaVPeak and PICKY achieves 88% recall and 83% precision, which significantly outperforms each individual method and the consensus method without using the B-H algorithm. The proposed method can be used as a standard procedure for any peak picking method and straightforwardly applied to some other prediction selection problems in bioinformatics. The source code, documentation and example data of the proposed method is available at http://sfb.kaust.edu.sa/pages/software.aspx.

Revealing the Potential Pathogenesis of Glioma by Utilizing a Glioma Associated Protein-Protein Interaction Network

“…Only the genes with adjusted p value, namely, false discovery rate (FDR)<0.0005 were screened out, and defined as DEGs. FDR<0.0005 implies that only less than 0.05 % of significant DEGs are false positive [15], indicating that the results of our DEGs screening are reliable.…”

“…Then the raw p values of t-test were adjusted by Benjamini-Hochberg algorithm [15]. Only the genes with adjusted p value, namely, false discovery rate (FDR)<0.0005 were screened out, and defined as DEGs.…”

Identification of MicroRNAs Associated with Nitrogen Use Efficiency and Fertility in Rice

“…It was further reported that miR169 is involved in nitrogen-starvation responses in Arabidopsis (Zhao et al, 2011). GmNFYA3 which is a target gene of miR169, has been reported to be a positive regulator of plant tolerance to drought stress (Abbas et al, 2013). MiR156k is responsive to salt and drought in rice and regulates flowering time and shoot development (Kruszka et al, 2012).…”

A Sequence-Based Dynamic Ensemble Learning System for Protein Ligand-Binding Site Prediction

“…So far, nuclear magnetic resonance (NMR) spectroscopy [2], [3], [4], [5], [6], [7], [8], [9] and X-ray crystallography [10] have been used to determine protein structures. Pintacuda et al employed lanthanide ions for the determination of protein-ligand binding sites [2].…”