Summary: | Abstract
Motivation: Sequence models based on deep neural networks have achieved state-of-the-art performance on regulatory genomics prediction tasks, such as chromatin accessibility and transcription factor binding. But despite their high
accuracy, their contributions to a mechanistic understanding of the biology of regulatory elements is often hindered
by the complexity of the predictive model and thus poor interpretability of its decision boundaries. To address this, we
introduce seqgra, a deep learning pipeline that incorporates the rule-based simulation of biological sequence data and
the training and evaluation of models, whose decision boundaries mirror the rules from the simulation process.
Results: We show that seqgra can be used to (i) generate data under the assumption of a hypothesized model of
genome regulation, (ii) identify neural network architectures capable of recovering the rules of said model and (iii)
analyze a model’s predictive performance as a function of training set size and the complexity of the rules behind
the simulated data.
Availability and implementation: The source code of the seqgra package is hosted on GitHub (https://github.com/gif
ford-lab/seqgra). seqgra is a pip-installable Python package. Extensive documentation can be found at https://
kkrismer.github.io/seqgra.
|