Discovering functional sequences with RELICS, an analysis method for tiling CRISPR screens

CRISPR screens are a powerful new technology for the identification of genome sequences that affect cellular phenotypes such as gene expression, survival, and proliferation. By tiling single-guide RNA (sgRNA) target sites across large genomic regions, CRISPR screens have the potential to systematically discovery novel functional sequences, however, a lack of purpose-built analysis tools limits the effectiveness of this approach. Here we describe RELICS, a Bayesian hierarchical model for the discovery of functional sequences from tiling CRISPR screens. RELICS considers the overlapping effects of multiple nearby functional sequences, accounts for the area of effect surrounding sgRNA target sites, models overdispersion in sgRNA counts, combines information across multiple pools, and estimates the number of functional sequences supported by the data. In simulations, RELICS outperforms existing methods and provides higher resolution predictions. We apply RELICS to published CRISPR interference and CRISPR activation screens and predict novel regulatory sequences, several of which we experimentally validate. In summary, RELICS is a powerful new analysis method for tiling CRISPR screens that enables the discovery of functional sequences with unprecedented resolution and accuracy.