Merging CAGE experiments
Merging CAGE experiments. This page describes the problem of merging independent CAGE-seq experiments and approaches to solving it.
Problem statement
Transcription of genes begins at genomic positions called transcription start sites (TSS). CAGE is a high-throughput transcriptome analysis technique that can identify active TSSs with one base resolution and their relative activities. It was shown by CAGE method that different sets of TSSs can operate under different conditions, and that transcription can start from several closely spaced TSSs within the promoter. All this complicates the comparative analysis of CAGE experiments carried out in different conditions. We have developed a method that allow us to combine independent CAGE experiments and to obtain a pooled set of TSSs with accurately defined boundaries. Iterative application of this method to a large set of CAGE experiments allows the construction of a reference TSS set. The presence of such a reference set makes it easy to compare TSS activities in different experiments, as well as to identify previously unknown TSS in the incoming data.
Algorithm overview
The method accepts two data sets (Reference and NewData) as input. Each of the sets consists of CAGE peaks and a corresponding full genome profile of the 5' ends of CAGE reads. The result is a set of non-overlapping NewReference peaks that reflect all TSSs from the input sets. If we intersect two sets of CAGE peaks (Reference and NewData) at genomic coordinates, the following types of peaks can be identified (Fig. 1): 1. Previously unknown - NewData peaks do not intersect with Reference peaks. 2. Not active in NewData - Reference peaks do not intersect with NewData. 3. Previously known, active in NewData - intersecting peaks Reference and NewData.
Results on rn6 data
Further directions
