Resolving a Systematic Error in STARR-seq for QuantitativeEnhancer Activity Mapping

Abstract

STARR-seq assesses millions of fragments in parallel measuring enhancer activity quantitatively. Here we show that STARR-seq is critically flawed with a systematic error in the cells of Arabidopsis thaliana (A. thaliana). Large amount of self-transcripts (STs) is lost during reverse transcription because these STs are polyadenylated after alternative polyadenylation sites (APAS) inside the test sequences. We solved this problem by using specially designed primer and recovered self-transcribed sequences independent from the PAS usage. In A. thaliana, we identified active enhancers and also enhancers quiescent in their endogenous genomic loci. Different from traditional STARR-seq identified enhancers, enhancers identified by new method are highly enriched in sequences proximal to the 5 and 3 ends of genes, and their epigenetic states correlate with gene expression levels. Our solution applies to methods based on self-transcript quantification. In addition, our results provide an invaluable functional enhancer activity map and insights into the functional complexity of enhancers in A. thaliana.