FUNNEL-GSEA: FUNctioNal ELastic-net Regression in Time-course Gene Set Enrichment Analysis

Overview

Journal Bioinformatics

Publisher Oxford University Press

Specialty Biology

Date 2017 Mar 24

PMID 28334094

Citations 13

Authors

Yun Zhang

David J Topham

Juilee Thakar

Xing Qiu

Affiliations

Soon will be listed here.

Abstract

Motivation: Gene set enrichment analyses (GSEAs) are widely used in genomic research to identify underlying biological mechanisms (defined by the gene sets), such as Gene Ontology terms and molecular pathways. There are two caveats in the currently available methods: (i) they are typically designed for group comparisons or regression analyses, which do not utilize temporal information efficiently in time-series of transcriptomics measurements; and (ii) genes overlapping in multiple molecular pathways are considered multiple times in hypothesis testing.

Results: We propose an inferential framework for GSEA based on functional data analysis, which utilizes the temporal information based on functional principal component analysis, and disentangles the effects of overlapping genes by a functional extension of the elastic-net regression. Furthermore, the hypothesis testing for the gene sets is performed by an extension of Mann-Whitney U test which is based on weighted rank sums computed from correlated observations. By using both simulated datasets and a large-scale time-course gene expression data on human influenza infection, we demonstrate that our method has uniformly better receiver operating characteristic curves, and identifies more pathways relevant to immune-response to human influenza infection than the competing approaches.

Availability And Implementation: The methods are implemented in R package FUNNEL, freely and publicly available at: https://github.com/yunzhang813/FUNNEL-GSEA-R-Package .

Contact: xing_qiu@urmc.rochester.edu or juilee_thakar@urmc.rochester.edu.

Supplementary Information: Supplementary data are available at Bioinformatics online.

Citing Articles

Alzheimer's disease: an integrative bioinformatics and machine learning analysis reveals glutamine metabolism-associated gene biomarkers.

Xing N, Yan J, Gao R, Zhang A, He H, Zheng M BMC Pharmacol Toxicol. 2025; 26(1):19.

PMID: 39875978 PMC: 11776333. DOI: 10.1186/s40360-025-00852-z.

Purine metabolism-related genes and immunization in thyroid eye disease were validated using bioinformatics and machine learning.

Wu Z, Gao Y, Cao L, Peng Q, Yao X Sci Rep. 2023; 13(1):18391.

PMID: 37884559 PMC: 10603126. DOI: 10.1038/s41598-023-45048-9.

SMYD3 associates with the NuRD (MTA1/2) complex to regulate transcription and promote proliferation and invasiveness in hepatocellular carcinoma cells.

Yang Y, Qiu R, Zhao S, Shen L, Tang B, Weng Q BMC Biol. 2022; 20(1):294.

PMID: 36575438 PMC: 9795622. DOI: 10.1186/s12915-022-01499-6.

FastMix: a versatile data integration pipeline for cell type-specific biomarker inference.

Zhang Y, Sun H, Mandava A, Aevermann B, Kollmann T, Scheuermann R Bioinformatics. 2022; 38(20):4735-4744.

PMID: 36018232 PMC: 9801972. DOI: 10.1093/bioinformatics/btac585.

QPRT Acts as an Independent Prognostic Factor in Invasive Breast Cancer.

Zhou L, Mu L, Jiang W, Yang Q J Oncol. 2022; 2022:6548644.

PMID: 35251175 PMC: 8894024. DOI: 10.1155/2022/6548644.

References

Cheng S . Thyroid hormone receptor mutations and disease: beyond thyroid hormone resistance. Trends Endocrinol Metab. 2005; 16(4):176-82. DOI: 10.1016/j.tem.2005.03.008. View

Shimomura Y, Wajid M, Ishii Y, Shapiro L, Petukhova L, Gordon D . Disruption of P2RY5, an orphan G protein-coupled receptor, underlies autosomal recessive woolly hair. Nat Genet. 2008; 40(3):335-9. DOI: 10.1038/ng.100. View

Qiu X, Klebanov L, Yakovlev A . Correlation between gene expression levels and limitations of the empirical bayes methodology for finding differentially expressed genes. Stat Appl Genet Mol Biol. 2006; 4:Article34. DOI: 10.2202/1544-6115.1157. View

Wang Z, Gerstein M, Snyder M . RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2008; 10(1):57-63. PMC: 2949280. DOI: 10.1038/nrg2484. View

Qiu X, Yakovlev A . Some comments on instability of false discovery rate estimation. J Bioinform Comput Biol. 2006; 4(5):1057-68. DOI: 10.1142/s0219720006002338. 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

Ruparelia N, Godec J, Lee R, Chai J, DallArmellina E, McAndrew D . Acute myocardial infarction activates distinct inflammation and proliferation pathways in circulating monocytes, prior to recruitment, and identified through conserved transcriptional responses in mice and humans. Eur Heart J. 2015; 36(29):1923-34. PMC: 4571177. DOI: 10.1093/eurheartj/ehv195. View

Garber M, Grabherr M, Guttman M, Trapnell C . Computational methods for transcriptome annotation and quantification using RNA-seq. Nat Methods. 2011; 8(6):469-77. DOI: 10.1038/nmeth.1613. View

Nam D . De-correlating expression in gene-set analysis. Bioinformatics. 2010; 26(18):i511-6. PMC: 2935420. DOI: 10.1093/bioinformatics/btq380. View

10.

Kanehisa M, Goto S . KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999; 28(1):27-30. PMC: 102409. DOI: 10.1093/nar/28.1.27. View

11.

Sohn I, Owzar K, George S, Kim S, Jung S . A permutation-based multiple testing method for time-course microarray experiments. BMC Bioinformatics. 2009; 10:336. PMC: 2772858. DOI: 10.1186/1471-2105-10-336. View

12.

Nueda M, Sebastian P, Tarazona S, Garcia-Garcia F, Dopazo J, Ferrer A . Functional assessment of time course microarray data. BMC Bioinformatics. 2009; 10 Suppl 6:S9. PMC: 2697656. DOI: 10.1186/1471-2105-10-S6-S9. View

13.

Luo W, Friedman M, Shedden K, Hankenson K, Woolf P . GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009; 10:161. PMC: 2696452. DOI: 10.1186/1471-2105-10-161. View

14.

Wu D, Smyth G . Camera: a competitive gene set test accounting for inter-gene correlation. Nucleic Acids Res. 2012; 40(17):e133. PMC: 3458527. DOI: 10.1093/nar/gks461. View

15.

Zhang K, Wang H, Bathke A, Harrar S, Piepho H, Deng Y . Gene set analysis for longitudinal gene expression data. BMC Bioinformatics. 2011; 12:273. PMC: 3142525. DOI: 10.1186/1471-2105-12-273. View

16.

Robinson M, McCarthy D, Smyth G . edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009; 26(1):139-40. PMC: 2796818. DOI: 10.1093/bioinformatics/btp616. View

17.

Thakar J, Mohanty S, Phillip West A, Joshi S, Ueda I, Wilson J . Aging-dependent alterations in gene expression and a mitochondrial signature of responsiveness to human influenza vaccination. Aging (Albany NY). 2015; 7(1):38-52. PMC: 4356402. DOI: 10.18632/aging.100720. View

18.

Wang L, Chen G, Li H . Group SCAD regression analysis for microarray time course gene expression data. Bioinformatics. 2007; 23(12):1486-94. DOI: 10.1093/bioinformatics/btm125. View

19.

Oron A, Jiang Z, Gentleman R . Gene set enrichment analysis using linear models and diagnostics. Bioinformatics. 2008; 24(22):2586-91. PMC: 2579710. DOI: 10.1093/bioinformatics/btn465. View

20.

Lee S, Li Z, Wu Z, Aau M, Guan P, Karuturi R . Context-specific regulation of NF-κB target gene expression by EZH2 in breast cancers. Mol Cell. 2011; 43(5):798-810. DOI: 10.1016/j.molcel.2011.08.011. View