GREAT improves functional interpretation of cis-regulatory regions - PubMed (original) (raw)
. 2010 May;28(5):495-501.
doi: 10.1038/nbt.1630. Epub 2010 May 2.
Affiliations
- PMID: 20436461
- PMCID: PMC4840234
- DOI: 10.1038/nbt.1630
GREAT improves functional interpretation of cis-regulatory regions
Cory Y McLean et al. Nat Biotechnol. 2010 May.
Abstract
We developed the Genomic Regions Enrichment of Annotations Tool (GREAT) to analyze the functional significance of cis-regulatory regions identified by localized measurements of DNA binding events across an entire genome. Whereas previous methods took into account only binding proximal to genes, GREAT is able to properly incorporate distal binding sites and control for false positives using a binomial test over the input genomic regions. GREAT incorporates annotations from 20 ontologies and is available as a web application. Applying GREAT to data sets from chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-seq) of multiple transcription-associated factors, including SRF, NRSF, GABP, Stat3 and p300 in different developmental contexts, we recover many functions of these factors that are missed by existing gene-based tools, and we generate testable hypotheses. The utility of GREAT is not limited to ChIP-seq, as it could also be applied to open chromatin, localized epigenomic markers and similar functional data sets, as well as comparative genomics sets.
Figures
Figure 1
Enrichment analysis of a set of _cis_-regulatory regions. (a) The current prevailing methodology associates only proximal binding events with genes and performs a gene-list test of functional enrichments using tools originally designed for microarray analysis. (b) GREAT’s binomial approach over genomic regions uses the total fraction of the genome associated with a given ontology term (green bar) as the expected fraction of input regions associated with the term by chance.
Figure 2
Binding profiles and their effects on statistical tests. (a) ChIP-seq data sets of several regulatory proteins show that the majority of binding events lie well outside the proximal promoter, both for sequence-specific transcription factors (SRF and NRSF, ref. ; Stat3, ref. 43) and a general enhancer-associated protein (p300, refs. 33,43). Cell type is given in parentheses: H, human; M, mouse. (b) When not restricted to proximal promoters, the gene-based hypergeometric test (red) generates false positive enriched terms, especially at the size range of 1,000–50,000 input regions typical of a ChIP-seq set. Negligible false positive enrichment was observed for the region-based binomial test (blue). For each set size, we generated 1,000 random input sets in which each base pair in the human genome was equally likely to be included in each set, avoiding assembly gaps. We calculated all GO term enrichments for both hypergeometric and binomial tests using GREAT’s 5+1 kb basal promoter and up to 1 Mb extension association rule (see Results). Plotted is the average number of terms artificially significant at a threshold of 0.05 after application of the conservative Bonferroni correction. (c) GO enrichment P values using the genomic region-based binomial (x axis) and gene-based hypergeometric (y axis) tests on the SRF data8 with GREAT’s 5+1 kb basal promoter and up to 1 Mb extension association rule (see Results). b1 through b10 denote the top ten most enriched terms when we used the binomial test. h1 through h10 denote the top ten most enriched terms when we used the hypergeometric test. Terms significant by both tests (B ∩ H) provide specific and accurate annotations supported by multiple genes and binding events (Table 3). Terms significant by only the hypergeometric test (H\B) are general and often associated with genes of large regulatory domains, whereas terms significant by only the binomial test (B\H) cluster four to six genomic regions near only one or two genes annotated with the term (Supplementary Table 46).
Figure 3
Distal binding events contribute substantially to accurate functional enrichments of p300 limb peaks. We examined properties of the 2,105 p300 mouse embryonic limb peaks in the context of three known limb-related terms and a negative control term (GO cortical cytoskeleton). Three different association rules were used (see Results): a gene-based GREAT analysis using only peaks within 2 kb of the nearest transcription start site (labeled 2 kb), an analysis with 5+1 kb basal and up to 50 kb extension (50 kb), and an analysis with 5+1 kb basal and up to 1 Mb extension (1 Mb). For each term, we examined the relevance of distal binding peaks by comparing the experimental results (black bars) to the average values of 1,000 simulated data sets (gray bars) in which the 192 proximal ChIP-seq peaks within 2 kb of the nearest transcription start site were fixed and the 1,913 distal peaks were shuffled uniformly within the mouse genome, avoiding assembly gaps and proximal promoters. By design, simulation results for proximal, 2-kb GREAT are identical to the actual data and are thus omitted. (a) Lengthening a 2-kb proximal promoter to a 50-kb extension, expected to increase genome coverage per term (_p_π in Fig. 1b) by 25-fold, causes an actual increase of 19- to 24-fold; in contrast, lengthening a 50-kb extension rule to a 1-Mb extension rule, expected to raise genome coverage 20-fold, leads to an actual increase of only 2.5- to 6-fold because regulatory domains are not extended through neighboring genes. (b) As regulatory domains increase in length from only the proximal 2 kb up to 50 kb and 1 Mb, the number of relevant genes with a p300 limb peak in their regulatory domain increases. The added genes selected only by distal associations are typically enriched for limb functionality compared to simulated data. (c) As regulatory domains increase in length, the number of p300 limb peaks associated with a relevant gene in excess of the number expected by chance increases for all limb-related terms. (d) As in c, the inclusion of distal peaks markedly increases the statistical significance of the correct terms alone. *Statistical significance is measured using the hypergeometric test over genes for 2 kb to mimic current gene-based approaches, and using the binomial test over genomic regions for 50 kb and 1 Mb. Error bars indicate s.d.; NS, not significant at a threshold of 0.05 after false discovery rate multiple test correction; obs, observed; exp, expected. Note scale changes on x axes.
Similar articles
- Genome-wide identification of in vivo protein-DNA binding sites from ChIP-Seq data.
Jothi R, Cuddapah S, Barski A, Cui K, Zhao K. Jothi R, et al. Nucleic Acids Res. 2008 Sep;36(16):5221-31. doi: 10.1093/nar/gkn488. Epub 2008 Aug 6. Nucleic Acids Res. 2008. PMID: 18684996 Free PMC article. - Genome-wide analysis of transcription factor binding sites based on ChIP-Seq data.
Valouev A, Johnson DS, Sundquist A, Medina C, Anton E, Batzoglou S, Myers RM, Sidow A. Valouev A, et al. Nat Methods. 2008 Sep;5(9):829-34. doi: 10.1038/nmeth.1246. Nat Methods. 2008. PMID: 19160518 Free PMC article. - Annotation of gene promoters by integrative data-mining of ChIP-seq Pol-II enrichment data.
Gupta R, Wikramasinghe P, Bhattacharyya A, Perez FA, Pal S, Davuluri RV. Gupta R, et al. BMC Bioinformatics. 2010 Jan 18;11 Suppl 1(Suppl 1):S65. doi: 10.1186/1471-2105-11-S1-S65. BMC Bioinformatics. 2010. PMID: 20122241 Free PMC article. - ChIP-Seq Data Analysis to Define Transcriptional Regulatory Networks.
Pavesi G. Pavesi G. Adv Biochem Eng Biotechnol. 2017;160:1-14. doi: 10.1007/10_2016_43. Adv Biochem Eng Biotechnol. 2017. PMID: 28070596 Review. - Identification and computational analysis of gene regulatory elements.
Taher L, Narlikar L, Ovcharenko I. Taher L, et al. Cold Spring Harb Protoc. 2015 Jan 5;2015(1):pdb.top083642. doi: 10.1101/pdb.top083642. Cold Spring Harb Protoc. 2015. PMID: 25561628 Free PMC article. Review.
Cited by
- Droplet Hi-C enables scalable, single-cell profiling of chromatin architecture in heterogeneous tissues.
Chang L, Xie Y, Taylor B, Wang Z, Sun J, Armand EJ, Mishra S, Xu J, Tastemel M, Lie A, Gibbs ZA, Indralingam HS, Tan TM, Bejar R, Chen CC, Furnari FB, Hu M, Ren B. Chang L, et al. Nat Biotechnol. 2024 Oct 18. doi: 10.1038/s41587-024-02447-1. Online ahead of print. Nat Biotechnol. 2024. PMID: 39424717 - A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival.
Abassah-Oppong S, Zoia M, Mannion BJ, Rouco R, Tissières V, Spurrell CH, Roland V, Darbellay F, Itum A, Gamart J, Festa-Daroux TA, Sullivan CS, Kosicki M, Rodríguez-Carballo E, Fukuda-Yuzawa Y, Hunter RD, Novak CS, Plajzer-Frick I, Tran S, Akiyama JA, Dickel DE, Lopez-Rios J, Barozzi I, Andrey G, Visel A, Pennacchio LA, Cobb J, Osterwalder M. Abassah-Oppong S, et al. Nat Commun. 2024 Oct 10;15(1):8793. doi: 10.1038/s41467-024-53009-7. Nat Commun. 2024. PMID: 39389973 Free PMC article. - Dynamic landscape of m6A modifications and related post-transcriptional events in muscle-invasive bladder cancer.
Zhang L, Chen Z, Sun G, Li C, Wu P, Xu W, Zhu H, Zhang Z, Tang Y, Li Y, Li Y, Xu S, Li H, Chen M, Xiao F, Zhang Y, Zhang W. Zhang L, et al. J Transl Med. 2024 Oct 8;22(1):912. doi: 10.1186/s12967-024-05701-x. J Transl Med. 2024. PMID: 39380003 Free PMC article. - An epigenetic human cytomegalovirus infection score predicts viremia risk in seropositive lung transplant recipients.
Hsu FM, Mohanty RP, Rubbi L, Thompson M, Pickering H, Reed EF, Greenland JR, Schaenman JM, Pellegrini M. Hsu FM, et al. Epigenetics. 2024 Dec;19(1):2408843. doi: 10.1080/15592294.2024.2408843. Epub 2024 Oct 3. Epigenetics. 2024. PMID: 39360678 Free PMC article. - ZIC2 and ZIC3 promote SWI/SNF recruitment to safeguard progression towards human primed pluripotency.
Hossain I, Priam P, Reynoso SC, Sahni S, Zhang XX, Côté L, Doumat J, Chik C, Fu T, Lessard JA, Pastor WA. Hossain I, et al. Nat Commun. 2024 Oct 2;15(1):8539. doi: 10.1038/s41467-024-52431-1. Nat Commun. 2024. PMID: 39358345 Free PMC article.
References
- Johnson DS, Mortazavi A, Myers RM, Wold B. Genome-wide mapping of in vivo protein-DNA interactions. Science. 2007;316:1497–1502. - PubMed
- Mardis ER. ChIP-seq: welcome to the new frontier. Nat. Methods. 2007;4:613–614. - PubMed
Publication types
MeSH terms
Substances
LinkOut - more resources
Full Text Sources
Other Literature Sources
Miscellaneous