Difference between revisions of "GTRD"
Ivan Yevshin (Talk | contribs) (→Web interface) |
Ivan Yevshin (Talk | contribs) (→General info) |
||
| Line 12: | Line 12: | ||
* [http://genome.ucsc.edu/ENCODE/ ENCODE]. | * [http://genome.ucsc.edu/ENCODE/ ENCODE]. | ||
| − | Initial | + | Initial raw data were systematically collected and uniformly processed using specially developed workflow (pipeline) for BioUML platform: |
- sequenced reads were aligned to reference genome using Bowtie; | - sequenced reads were aligned to reference genome using Bowtie; | ||
- peaks were identified using MACS and SISSR algorithms | - peaks were identified using MACS and SISSR algorithms | ||
Revision as of 12:02, 31 August 2015
GTRD (Gene Transcription Regulation Database) is a database of transcription factor binding sites identified from ChIP-seq experiments that were systematically collected and uniformly processed using special workflow (pipeline) for BioUML platform.
Contents |
General info
GTRD (Gene Transcription Regulation Database) is a database of transcription factor binding sites identified from ChIP-seq experiments from:
Initial raw data were systematically collected and uniformly processed using specially developed workflow (pipeline) for BioUML platform: - sequenced reads were aligned to reference genome using Bowtie; - peaks were identified using MACS and SISSR algorithms - further refinement of obtained peaks - position weight matrices (PWM) were constructed by different methods (ChIPMunk, our own methods) - ROC curves were calculated to estimate and compare built PWM - site models (PWMs + thresholds) were constucted for recognition TF binding sites.
GTRD database is freely available for non-commercial organizations via web interface.
Database statistics
GTRD uses 2417 ChIP-seq experiments for 470 distinct sequence specific transcription factors.
General statistics:
| Object type | Total count | Per ChIP-seq experiment |
|---|---|---|
| ChIP-seq reads | 80.8 × 109 | 34.9 × 106 |
| Reads aligned | 58.8 × 109 | 25.6 × 106 |
| ChIP-seq peaks | 59.5 × 106 | 32 899 |
In average each transcription factor is measured in 4.07 ChIP-seq experiments, but 284 (60%) transcription factors measured only in one experiment.
The ten most studied transcription factors are listed bellow:
| Transcription Factor | Number of ChIP-seq experiments |
|---|---|
| CTCF | 195 |
| c-Myc | 45 |
| ERα | 44 |
| NRSF | 37 |
| C/EBPβ | 37 |
| GATA-1 | 33 |
| NF-κB p65 | 30 |
| Max | 30 |
| PU.1 | 29 |
| GR | 24 |
Database structure
The metadata concerning GTRD is stored in MySQL tables.
Each ChIP-seq experiment has a row in 'chip_experiments' table, which assigns id and stores basic information about experiment. 'chip_experiments' table has following structure:
| Column | Description | Example value |
|---|---|---|
| id | Unique experiment identifier | EXP000489 |
| antibody | Antibody used in chromatin immunoprecipitation | sc-345 |
| tfClassId | Id in TFClass[1] database of target transcription factor, NULL for control experiments | 6.2.1.0.1 |
| cell_line | Studied cell line | HeLa S3 |
| specie | Species latin name | Homo sapiens |
| treatment | Cell treatment or conditions | IFN gamma |
| control_id | Id of control experiment, NULL for control experiments or experiments without control | EXP000490 |
The links to external databases stored in 'external_refs' table:
| Column | Description | Example values |
|---|---|---|
| id | Experiment identifier | EXP000489 |
| external_db | External database name | GEO or PUBMED or ENCODE or SRA |
| external_db_id | Identifier in external database | GSM320736 |
GTRD uses following object identifiers:
| Template | Object type | Example |
|---|---|---|
| EXPXXXXXX | ChIP-seq experiment | EXP000489 |
| READSXXXXXX | Collection of ChIP-seq reads | READS000770 |
| ALIGNSXXXXXX | Collection of read alignments | ALIGNS010001 |
| PEAKSXXXXXX | Collection of ChIP-seq peaks | PEAKS010000 |
The relationship between these objects is provided by 'hub' table:
| Column | Description | Example values |
|---|---|---|
| input | Input object identifier | READS000770 |
| input_type | Type of input object | ReadsGTRDType |
| output | Output object identifier | EXP000489 |
| output_type | Type of output object | ExperimentGTRDType |
ChIP-seq reads, alignments and peaks links to experiments with hub table in the following way:
| input | input_type | output | output_type |
|---|---|---|---|
| READS000770 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| READS000771 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| READS000772 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| READS000773 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| READS000774 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| READS000775 | ReadsGTRDType | EXP000489 | ExperimentGTRDType |
| ALIGNS010001 | AlignmentsGTRDType | EXP000489 | ExperimentGTRDType |
| PEAKS010000 | PeaksGTRDType | EXP000489 | ExperimentGTRDType |
Web interface
| This page or section is a stub. Please add screenshots here! |
Web interface to GTRD is available here. It provides capabilities for searching and browsing GTRD.
Start page
Start page contains search box and the links to browse all experiments in the repository tree, explore databases statistics and transcription factor classification tree.
Search capabilities
The ChIP-seq experiments contained in GTRD can be queried from the search box on the Start page. GTRD uses Lucene engine for indexing and quering ChIP-seq experiments that provides rich syntax for searching. The search can be performed by transcription factor(name or class), cell line, antibody or treatment/conditions.
For example, to search for STAT transcription factors enter stat* in the search field and press enter. You can restrict query to HeLa cells treated with interferon with following query:
tfTitle:stat* AND cellLine:hela AND treatment:IFN
The list of matching ChIP-seq experiments will appear in the 'Search result' tab. Select one of them to view detailed ChIP-seq experiment information in the information box. The information box also provides links to the experiment data: reads, alignments and peaks.
To view peaks or alignments click the link in the ChIP-seq experiment information box, the 
track will be opened as table. The track can be exported by pressing 'Export' button or opened in the genome browser by pressing 'Open as track' button in the general control panel.
Repository structure
GTRD is organized in hierarchical Repository.
The GTRD/Data folder contains following items:
- experiments - ChIP-seq experiments metainformation
- sequences - Raw ChIP-seq reads in fastq.gz format
- alignments - ChIP-seq read alignments
- peaks - ChIP-seq peaks identified by MACS and SISSRs peak callers
- matrices - Position weight matrices for transcription factor binding sites
- site models - Models for recognition of transcription factor binding sites
- statistics - Summary GTRD statistics
- views - TFClass centric view of GTRD data avalable as
tree-tables
GTRD/Dictionaries/classification is TFClass classification tree used by GTRD to reference transcription factors.
