Source-linked AI summary

Comprehensive Evaluation of Deep Learning Architectures for Prediction of DNA/RNA Sequence Binding Specificities

Ameni Trabelsi, Mohamed Chaabane, Asa Ben Hur

arXiv:1901.10526v1cs.LGq-bio.QMstat.ML

TL;DR

The best-performing architecture for DNA/RNA binding prediction was unclear, and comparisons were complicated by differences in model selection. The authors developed deepRAM to compare architectures with a shared automatic procedure; hybrid CNN/RNN models and k-mer embeddings performed well, with accuracy gains accompanied by interpretability and training-time costs.

  • Problem

    It was unclear which deep-learning architecture best predicts DNA/RNA binding, while model-selection differences made fair comparisons difficult.

    Source: §1 · p. 1, §2.1 · p. 2
  • Method

    The authors used deepRAM to systematically compare CNN, RNN, and hybrid architectures with uniform automatic model selection across DNA- and RNA-binding datasets.

    Source: §2.1 · p. 2, §2.2 · p. 2
  • Results

    Across the benchmarks, ECBLSTM was the strongest model; its median AUC was 0.930 versus 0.902 for DeepBind on ChIP-seq data.

    Source: §3.1 · p. 5, §4 · p. 7
  • Takeaways & Limitations

    Hybrid models and k-mer embeddings improved performance, but the accuracy gains came with reduced interpretability and longer training times.

    Source: §2.2 · p. 2, §4 · p. 7

    Hybrid CNN/RNN models require sufficient training data and are harder to interpret, which may limit their usefulness for motif discovery.

    Source: §4 · p. 7

Abstract

from arXiv · show

Motivation: Deep learning architectures have recently demonstrated their power in predicting DNA- and RNA-binding specificities. Existing methods fall into three classes: Some are based on Convolutional Neural Networks (CNNs), others use Recurrent Neural Networks (RNNs), and others rely on hybrid architectures combining CNNs and RNNs. However, based on existing studies it is still unclear which deep learning architecture is achieving the best performance. Thus an in-depth analysis and evaluation of the different methods is needed to fully evaluate their relative. Results: In this study, We present a systematic exploration of various deep learning architectures for predicting DNA- and RNA-binding specificities. For this purpose, we present deepRAM, an end-to-end deep learning tool that provides an implementation of novel and previously proposed architectures; its fully automatic model selection procedure allows us to perform a fair and unbiased comparison of deep learning architectures. We find that an architecture that uses k-mer embedding to represent the sequence, a convolutional layer and a recurrent layer, outperforms all other methods in terms of model accuracy. Our work provides guidelines that will assist the practitioner in choosing the best architecture for the task at hand, and provides some insights on the differences between the models learned by convolutional and recurrent networks. In particular, we find that although recurrent networks improve model accuracy, this comes at the expense of a loss in the interpretability of the features learned by the model. Availability and implementation: The source code for deepRAM is available at https://github.com/MedChaabane/deepRAM

1 Introduction

The paper addresses the lack of a clear, fair comparison among deep-learning architectures for DNA/RNA binding prediction, where binding data are costly to obtain and simple sequence models may miss specificity patterns. It introduces a systematic comparison across benchmark datasets and sequence representations, reporting improved performance from hybrid models and k-mer embeddings.

1 Introduction

  • Only some DNA- and RNA-binding sites have been experimentally identified, and experiments are expensive and time-consuming; PWMs are common but more complex models may capture specificity better.The paper notes prior evidence that sequence specificity can be better captured with more complex models.Source: §1 · p. 1

2.1 Datasets

  • The evaluation uses 83 ENCODE ChIP-seq experiments and 31 human CLIP-seq experiments covering 19 proteins, with defined training and test examples for both benchmarks.ChIP-seq positives are 101-bp peak-centered regions with dinucleotide-matched shuffled negatives; CLIP-seq experiments each contain 30,000 training and 10,000 test examples.Source: §2.1 · p. 2

1 Introduction

  • A fair architecture comparison requires uniform model selection because deep-learning performance is sensitive to hyperparameter tuning; the study therefore uses a fully automatic, shared selection strategy.The architecture variants differ in network components, depth, and input representation.Source: §2.1 · p. 2

2.2 Model architectures

  • The study compares one-hot encoding with word2vec k-mer embeddings, alongside CNN, RNN, and hybrid CNN/RNN architectures for DNA/RNA sequences.The k-mer representation maps overlapping k-mers to vectors whose contexts are preserved.Source: §2.2 · p. 2, §2.1 · p. 2
  • Hybrid CNN/RNN models improved performance over vanilla CNNs, with the advantage increasing as more training examples were available; k-mer embeddings also helped, especially for RBP prediction.The authors note that improved accuracy comes with reduced interpretability and longer training times.Source: §2.2 · p. 2

2 Methods

The methods describe convolutional, recurrent, and hybrid sequence models, then compare architecture variants under a shared automatic model-selection procedure. deepRAM packages these models and supports motif extraction and visualization.

Convolutional Networks.

  • CNN filters scan local sequence patterns, while deeper convolution and pooling layers can capture features over larger spatial ranges, including motif interactions.The first convolutional layer can be viewed as a motif scanner; pooling aggregates nearby filter outputs to provide consistency to small shifts.Source: Convolutional Networks · pp. 2–3

RNN-based models.

  • RNNs update an internal state along the sequence to capture interactions between distant elements; the evaluated units include LSTM and GRU, with bidirectional variants processing both sequence directions.The models end with fully connected layer(s) and a sigmoid that estimates binding-site probability.Source: RNN-based models · p. 4, §2.3 · p. 4

Evaluated architectures

  • deepRAM compares existing models including DeepBind, DanQ, KEGRU, and dilated CNNs, alongside variants that alter convolution depth, embedding, and recurrent layers.The variants include DeepBind*, DanQ*, DeepBind-E*, ECLSTM, and ECBLSTM.Source: Evaluated architectures · p. 4, Table 1 · p. 4

2.3 Model training, selection, and evaluation

  • For each dataset, automatic selection samples 40 hyperparameter settings and scores each by average AUC across three-fold cross-validation before training five models with the best setting.The final model is selected using training performance among those five models, then evaluated on the test set.Source: §2.3 · p. 4

Model training.

  • deepRAM lets users configure input representation and CNN/RNN components, automatically extracts first-layer motifs, and provides an open-source package for model evaluation and analysis.Model-selection runs range from about one hour for a single-layer CNN to almost four hours for a convolutional and bidirectional-LSTM network on the stated setup.Source: Model training · p. 5

3 Results

Across DNA- and RNA-binding prediction tasks, ECBLSTM achieved higher median AUCs than DeepBind, and embedding, recurrent layers, and sufficient training data were associated with stronger performance. Dilated convolution also outperformed DeepBind*, while recurrent features were less interpretable than CNN features.

3.1 Deeper is better

  • Adding multiple convolutional layers, dilated convolution, and embedding improved performance over DeepBind, while recurrent layers further improved performance by capturing long-term dependencies.DanQ's bidirectional LSTM also improved its performance over DeepBind.Source: §3.1 · p. 5

2.5 deepRAM

  • deepRAM provides an end-to-end toolkit for binding-site and motif prediction; its DeepBind implementation nearly matched the original, and deep-learning methods outperformed iONMF.iDeepS reached a 0.917 median AUC on CLIP-seq, below all evaluated methods except DeepBind.Source: §3.2 · p. 5, §2.5 · p. 5

3.1 Deeper is better

  • Median AUC was 0.930 for ECBLSTM versus 0.902 for DeepBind on ChIP-seq, and 0.951 versus 0.914 on CLIP-seq.The figure summarizes AUC distributions and pairwise model comparisons across both tasks.Source: §3.1 · p. 5, Figure 2 · p. 6

3.3 Deeper is better with sufficient training data

  • Hybrid models ECBLSTM, ECLSTM, and DanQ* tended to exceed median AUC 0.983 on large datasets but fell behind DeepBind-E* on smaller datasets.The authors suggest hybrid models need sufficient training data; ECBLSTM still performed well on smaller datasets, which they attribute to effective regularization.Source: §3.2 · p. 6

3.2 k-mer embedding boosts model performance

  • On RBP datasets, embedding models all had median AUC above 0.94, versus below 0.935 for one-hot models; k-mer co-occurrence encodes contextual information.The authors used k = 3 and stride s = 1, which performed best in preliminary experiments.Source: §3.2 · pp. 5–6

3.3 Deeper is better with sufficient training data

  • For datasets with more than 10,000 positive training samples, median AUCs ranged from 0.967 (DeepBind) to 0.993 (ECBLSTM), versus 0.864 (DeepBind) to 0.879 (DeepBind-E*) below 10,000.The figure separates 38 experiments with fewer than 10,000 peaks from 45 with more than 10,000 peaks.Source: §3.2 · p. 6, Figure 3 · p. 7

3.4 Dilated convolution

  • DeepBind detected motifs representing binding events, whereas ECBLSTM's convolution stage extracted features spanning the whole sequence, consistent with reduced interpretability from RNNs.Motifs were extracted from first-layer filters by collecting strongly activating sequence fragments and forming position frequency matrices.Source: §3.1 · p. 5, §3.4 · p. 7
  • The three-layer dilated model, with dilations 1, 2, and 2, significantly outperformed DeepBind* on both tasks and had slightly higher median AUC than DanQ on RBP datasets.The latter comparison suggests dilated convolution can capture long-range relationships similarly to LSTMs.Source: §3.4 · p. 7

4 Conclusion

Hybrid CNN/RNN models perform accurately, but their learned features may be less useful for motif discovery because they are harder to interpret. deepRAM provides a tool for evaluating and analyzing DNA/RNA binding architectures.

  • Hybrid CNN/RNN accuracy requires sufficient training data, and their harder interpretation may limit usefulness in motif discovery.Source: §4 · p. 7
  • deepRAM packages the study’s architectures in a user-friendly tool for evaluating and analyzing DNA/RNA binding prediction models.Source: §4 · p. 7
  • ECBLSTM was the best-performing compared model, yet detected fewer informative motifs than DeepBind and DeepBind-E*.The authors hypothesize that CNN filters in the hybrid model learn information for the recurrent layer rather than localized motifs.Source: §4 · p. 7
Loading 1901.10526v1…