Difference between revisions of "Merging CAGE experiments"

Revision as of 16:19, 4 March 2021

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.

Figure 1

The first two types of peaks go into the NewReference set without changes, and for the third, it is necessary to clarify the boundaries, since the intersection can be partial. When there is a partial intersection of the Reference and NewData peaks, overhanging ends (Fig. 2) as well as multiple intersections (Fig. 3, Fig. 4) can be observed.

Figure 2. CAGE peak overhangs.

Figure 3. One reference peak to several new peaks.

Figure 4. One new peak to several reference peaks.

A preliminary analysis of the rat CAGE data (Reference = FANTOM5, NewData = UEXP) showed a significant length of the overhangs (the average overhang length is 55% of the peak length for FANTOM5 peaks and 12% for UEXP peaks) as well as a significant proportion of multiple intersections (7.5% FANTOM5 peaks intersect with more than one UEXP peak and 1.1% of UEXP peaks intersect with more than one FANTOM5 peak). The presence of overhangs and multiple intersections can be caused both by the actual TSS activity in these regions in only one of the datasets or by the insufficient accuracy of determining the boundaries with the peak-caller due to insufficient coverage by reads. It is important to be able to distinguish between these cases, since in the first case this area should be included in the NewReference and in the second should not. This requires an analysis of the read density profile in these regions. To do this, the entire genome is divided into segments defined by the boundaries of the peaks of both datasets (Fig. 5).

Figure 5. Segmentation of CAGE peaks

Results on rn6 data

Further directions