A robust peak detection method for RNA structure inference by high-throughput contact mapping.
BACKGROUND
For high-throughput prediction of the helical arrangements of large RNA molecules, an innovative method termed multiplexed hydroxyl radical (*OH) cleavage analysis (MOHCA) has been proposed. A key step in this promising technique is to detect peaks accurately from noisy radioactivity profiles. Since manual peak finding is laborious and prone to error, an automated peak detection method to improve the accuracy and throughput of MOHCA is required. Existing methods were not applicable to MOHCA due to their high false positive rates.


RESULTS
We developed a two-step computational method that can detect peaks from MOHCA profiles in a robust manner. The first step exploits an ensemble of linear and non-linear signal processing techniques to find true peak candidates. In the second step, a binary classifier trained with the characteristics of true and false peaks is used to eliminate false peaks out of the peak candidates. We tested the proposed approach with 2002 MOHCA cleavage profiles and obtained the median recall, precision and F-measure values of 0.917, 0.750 and 0.830, respectively. Compared with the alternatives considered, the proposed method was able to handle false peaks substantially better, thus resulting in 51.0-71.8% higher median values of precision and F-measure.


BACKGROUND
The software and supplementary data are available at http://dna.korea.ac.kr/pub/mohca.
Kim
Jinkyu
J
School of Electrical Engineering, Korea University, Seoul 136-713, Korea.
Yu
Seunghak
S
Shim
Byonghyo
B
Kim
Hanjoo
H
Min
Hyeyoung
H
Chung
Eui-Young
EY
Das
Rhiju
R
Yoon
Sungroh
S
D001483:Base Sequence; D019295:Computational Biology / Q000379:methods*; D008958:Models, Molecular; D008969:Molecular Sequence Data; D009690:Nucleic Acid Conformation; D010363:Pattern Recognition, Automated; D012313:RNA / Q000737:chemistry*
D012313:RNA