Comparing Gene Set Analysis Methods on Single-nucleotide Polymorphism Data from Genetic Analysis Workshop 16

Overview

Journal BMC Proc

Publisher Biomed Central

Specialty Biology

Date 2009 Dec 19

PMID 20018093

Citations 33

Authors

Nathan L Tintle

Bryce Borchers

Marshall Brown

Airat Bekmetjev

Affiliations

Soon will be listed here.

Abstract

Recently, gene set analysis (GSA) has been extended from use on gene expression data to use on single-nucleotide polymorphism (SNP) data in genome-wide association studies. When GSA has been demonstrated on SNP data, two popular statistics from gene expression data analysis (gene set enrichment analysis [GSEA] and Fisher's exact test [FET]) have been used. However, GSEA and FET have shown a lack of power and robustness in the analysis of gene expression data. The purpose of this work is to investigate whether the same issues are also true for the analysis of SNP data. Ultimately, we conclude that GSEA and FET are not optimal for the analysis of SNP data when compared with the SUMSTAT method. In analysis of real SNP data from the Framingham Heart Study, we find that SUMSTAT finds many more gene sets to be significant when compared with other methods. In an analysis of simulated data, SUMSTAT demonstrates high power and better control of the type I error rate. GSA is a promising approach to the analysis of SNP data in GWAS and use of the SUMSTAT statistic instead of GSEA or FET may increase power and robustness.

Citing Articles

Live-bearing cockroach genome reveals convergent evolutionary mechanisms linked to viviparity in insects and beyond.

Fouks B, Harrison M, Mikhailova A, Marchal E, English S, Carruthers M iScience. 2023; 26(10):107832.

PMID: 37829199 PMC: 10565785. DOI: 10.1016/j.isci.2023.107832.

Integrative pathway and network analysis provide insights on flooding-tolerance genes in soybean.

Jhan L, Yang C, Huang C, Lai M, Huang Y, Baiya S Sci Rep. 2023; 13(1):1980.

PMID: 36737640 PMC: 9898312. DOI: 10.1038/s41598-023-28593-1.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation.

Treaster S, Deelen J, Daane J, Murabito J, Karasik D, Harris M Sci Adv. 2023; 9(2):eadd2743.

PMID: 36630509 PMC: 9833670. DOI: 10.1126/sciadv.add2743.

Signatures of Positive Selection in the Genome of : A Subspecies of European Honeybees.

Huang Q, Zhu Y, Fouks B, He X, Niu Q, Zheng H Life (Basel). 2022; 12(10).

PMID: 36295077 PMC: 9605442. DOI: 10.3390/life12101642.

Understanding signatures of positive natural selection in human zinc transporter genes.

Roca-Umbert A, Caro-Consuegra R, Londono-Correa D, Rodriguez-Lozano G, Vicente R, Bosch E Sci Rep. 2022; 12(1):4320.

PMID: 35279701 PMC: 8918337. DOI: 10.1038/s41598-022-08439-y.

References

Nam D, Kim S . Gene-set approach for expression pattern analysis. Brief Bioinform. 2008; 9(3):189-97. DOI: 10.1093/bib/bbn001. View

Dinu I, Potter J, Mueller T, Liu Q, Adewale A, Jhangri G . Improving gene set analysis of microarray data by SAM-GS. BMC Bioinformatics. 2007; 8:242. PMC: 1931607. DOI: 10.1186/1471-2105-8-242. View

Nettleton D, Recknor J, Reecy J . Identification of differentially expressed gene categories in microarray studies using nonparametric multivariate analysis. Bioinformatics. 2007; 24(2):192-201. DOI: 10.1093/bioinformatics/btm583. View

Wang K, Li M, Bucan M . Pathway-based approaches for analysis of genomewide association studies. Am J Hum Genet. 2007; 81(6):1278-83. PMC: 2276352. DOI: 10.1086/522374. View

Allison D, Cui X, Page G, Sabripour M . Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2005; 7(1):55-65. DOI: 10.1038/nrg1749. View

Subramanian A, Tamayo P, Mootha V, Mukherjee S, Ebert B, Gillette M . Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102(43):15545-50. PMC: 1239896. DOI: 10.1073/pnas.0506580102. View

Goeman J, Buhlmann P . Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007; 23(8):980-7. DOI: 10.1093/bioinformatics/btm051. View

Flicek P, Aken B, Beal K, Ballester B, Caccamo M, Chen Y . Ensembl 2008. Nucleic Acids Res. 2007; 36(Database issue):D707-14. PMC: 2238821. DOI: 10.1093/nar/gkm988. View

Tian L, Greenberg S, Kong S, Altschuler J, Kohane I, Park P . Discovering statistically significant pathways in expression profiling studies. Proc Natl Acad Sci U S A. 2005; 102(38):13544-9. PMC: 1200092. DOI: 10.1073/pnas.0506577102. View

10.

Amos C . Successful design and conduct of genome-wide association studies. Hum Mol Genet. 2007; 16 Spec No. 2:R220-5. PMC: 2691963. DOI: 10.1093/hmg/ddm161. View

11.

Khatri P, Draghici S . Ontological analysis of gene expression data: current tools, limitations, and open problems. Bioinformatics. 2005; 21(18):3587-95. PMC: 2435250. DOI: 10.1093/bioinformatics/bti565. View

12.

Chasman D . On the utility of gene set methods in genomewide association studies of quantitative traits. Genet Epidemiol. 2008; 32(7):658-68. DOI: 10.1002/gepi.20334. View

13.

Mootha V, Lindgren C, Eriksson K, Subramanian A, Sihag S, Lehar J . PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003; 34(3):267-73. DOI: 10.1038/ng1180. View

14.

Tintle N, Best A, DeJongh M, Van Bruggen D, Heffron F, Porwollik S . Gene set analyses for interpreting microarray experiments on prokaryotic organisms. BMC Bioinformatics. 2008; 9:469. PMC: 2587482. DOI: 10.1186/1471-2105-9-469. View