This document provides a technical description of the transcription start site assembly pipeline that was used to generate EPDnew version 003 for D. melanogaster.
|Name||Genome Assembly||Promoters||Genes||PMID||Access data|
|ENSEMBL70||Apr 2006 BDGP R5/dm3||18409||13660||19420058||SOURCE||DOC||DATA|
|Hoskins et al., 2012||CAGE||1||17,979,809||21177961||SOURCE||DOC||DATA|
|Ni et al., 2010||CAGE||2||27,173,616||20495556||SOURCE||DOC||DATA|
Data was downloaded from BioMart selecting the following attributes:
Then, transcrips have been filtered according to the following rules:
Gene names were taken from the field 'Associated Gene Name'. Since the EPD format doesn't allow gene names longer than 18 characters, we checked whether the names repsected this limitation.
Transcripts with the same TSS position were merged under a common ID. As a conseguence of this, from the 23850 transcrips originally present in the ENSEMBL database, 5953 were merged, leaving 17897 uniquely mapped promoters in the input list.
The ENSEMBL TSS collection is stored as a tab-deliminated text file conforming to the SGA format under the name:
The six field contain the following kinds of information:
Note that the second and forth fields are invariant.
MachiBase data were generated with the oligo-capping technology. The source data were downloaded from:
According to the readme file included in the tar archive, the 5' end tags were mapped to the Drosophila genome using BLAT as alignment tool allowing for up to three mismatches.
The compressed version of this file is available from the MGA archive (see above) under the name:
Mapped sequence tags were extracted from Supplementary Data File 1 available from Genome Research at:
The downloaded source file is in SAM format and has been generated with the tag mapping program StatMap as described in the article cyted above. We extracted all tags with mapping quality scores greater or equal to 30.
The compressed version of this file is available from our ftp site (see above link) with the name:
BAM files for the SRA serie SRP001602 and SRX018832 were downloaded from SRA site and converted into SGA file using in house software.
Each sample in the collection (mRNA peaks and ENSEMBL TSS) was then processed in a pipeline aiming at validating transcription start sites with mRNA peaks. An Ensembl TSS was experimentally confirmed if a CAGE peak lied in a window of 200 bp around it and if it had a maximum high of at least 3 tags. The validated TSS was then shifted to the nearest base with the higher tag density.
Each sample in the dataset was used to generate a separate promoter collection. Potentially, the same transcript could be validated by multiple samples and it could have different start sites in different samples. To avoid redundancy, the individual collections were used as input for an additional step in the analysis (Assembly pipeline part B).
All promoter collections were merged into a unique file and further analysed. The promoter of a transcript was mantained in the list only if validated by at least two samples. Transcript validated by multiple samples could potentially have the TSS set on a broader region and not to single position. To avoid such inconsistency, for each transcript we selected the position that was validated by the larger number of samples as the true TSS.
Transcription Start Sites that mapped closed to other TSS that belonged to the same gene (200 bp window) were merged into a unique promoter following the same rule: the promoter that was validated by the higher number of samples was kept.
The experimentally validated promoter were stored in the EPDnew database that can be downloaded from our ftp site. Scientist are wellcome to use our other tools ChIP-Seq (for correlation analysis) and SSA (for motifs analysis around promoters) to analyse EPDnew database.