Fold-Change-Specific Enrichment Analysis (FSEA): Quantification of Transcriptional Response Magnitude for Functional Gene Groups

Gene expression profiling data contains more information than is routinely extracted with standard approaches. Here we present Fold-Change-Specific Enrichment Analysis (FSEA), a new method for functional annotation of differentially expressed genes from transcriptome data with respect to their fold changes. FSEA identifies Gene Ontology (GO) terms, which are shared by the group of genes with a similar magnitude of response, and assesses these changes. GO terms found by FSEA are fold-change-specifically (e.g., weakly, moderately, or strongly) affected by a stimulus under investigation. We demonstrate that many responses to abiotic factors, mutations, treatments, and diseases occur in a fold-change-specific manner. FSEA analyses suggest that there are two prevailing responses of functionally-related gene groups, either weak or strong. Notably, some of the fold-change-specific GO terms are invisible by classical algorithms for functional gene enrichment, Singular Enrichment Analysis (SEA), and Gene Set Enrichment Analysis (GSEA). These are GO terms not enriched compared to the genome background but strictly regulated by a factor within specific fold-change intervals. FSEA analysis of a cancer-related transcriptome suggested that the gene groups with a tightly coordinated response can be the valuable source to search for possible regulators, markers, and therapeutic targets in oncogenic processes. Availability and Implementation: FSEA is implemented as the FoldGO Bioconductor R package and a web-server.

[1]  Mukesh Jain,et al.  Differential quantitative regulation of specific gene groups and pathways under drought stress in rice. , 2019, Genomics.

[2]  Tanya Barrett,et al.  The Gene Expression Omnibus Database , 2016, Statistical Genomics.

[3]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[4]  E. Giannoni,et al.  Targeting stromal-induced pyruvate kinase M2 nuclear translocation impairs OXPHOS and prostate cancer metastatic spread , 2015, Oncotarget.

[5]  Kara Dolinski,et al.  The BioGRID interaction database: 2019 update , 2018, Nucleic Acids Res..

[6]  R. Tibshirani,et al.  Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Gennady Korotkevich,et al.  Fast gene set enrichment analysis , 2021 .

[8]  Y. Rustum,et al.  Overexpression of Bax enhances antitumor activity of chemotherapeutic agents in human head and neck squamous cell carcinoma. , 2000, Clinical cancer research : an official journal of the American Association for Cancer Research.

[9]  G. Dallmann,et al.  Role of coenzymes in cancer metabolism. , 2020, Seminars in cell & developmental biology.

[10]  Silvio C. E. Tosatto,et al.  InterPro in 2019: improving coverage, classification and access to protein sequence annotations , 2018, Nucleic Acids Res..

[11]  J. Neoptolemos,et al.  Protein kinase C isoenzyme patterns characteristically modulated in early prostate cancer. , 1999, The American journal of pathology.

[12]  Teruhiko Fujii,et al.  Protein Kinase C Promotes Apoptosis in LNCaP Prostate Cancer Cells through Activation of p38 MAPK and Inhibition of the Akt Survival Pathway* , 2003, Journal of Biological Chemistry.

[13]  Zalmiyah Zakaria,et al.  A Review on Bioinformatics Enrichment Analysis Tools Towards Functional Analysis of High Throughput Gene Set Data , 2015 .

[14]  The Impact of Gene Expression Patterns in Breast Cancer. , 2016, Clinical chemistry.

[15]  Santhilata Kuppili Venkata,et al.  The Network of Cancer Genes (NCG): a comprehensive catalogue of known and candidate cancer genes from cancer sequencing screens , 2018, Genome Biology.

[16]  V. Mironova,et al.  Auxin regulates functional gene groups in a fold-change-specific manner in Arabidopsis thaliana roots , 2017, Scientific Reports.

[17]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[18]  J. Nickels,et al.  Cancerous inhibitor of protein phosphatase 2A promotes premature chromosome segregation and aneuploidy in prostate cancer cells through association with shugoshin , 2015, Tumor Biology.

[19]  Xin Lu,et al.  Integration of Metabolomics and Transcriptomics Reveals Major Metabolic Pathways and Potential Biomarker Involved in Prostate Cancer* , 2015, Molecular & Cellular Proteomics.

[20]  Alexey Sergushichev,et al.  An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation , 2016 .

[21]  P. Kapranov,et al.  On the importance of small changes in RNA expression. , 2013, Methods.

[22]  Minoru Kanehisa,et al.  KEGG: new perspectives on genomes, pathways, diseases and drugs , 2016, Nucleic Acids Res..

[23]  A. Flisser,et al.  Calreticulin in phagocytosis and cancer: opposite roles in immune response outcomes , 2019, Apoptosis.

[24]  Na Ye,et al.  Direct Activation of Bax Protein for Cancer Therapy , 2016, Medicinal research reviews.

[25]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[26]  The transcription factor ERG increases expression of neurotransmitter receptors on prostate cancer cells , 2015, BMC Cancer.

[27]  Gordon K Smyth,et al.  The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads , 2018, bioRxiv.

[28]  Chien‐Hung Yeh,et al.  FBXW7: a critical tumor suppressor of human cancers , 2018, Molecular Cancer.

[29]  J. Céraline,et al.  Dual effects of constitutively active androgen receptor and full-length androgen receptor for N-cadherin regulation in prostate cancer , 2017, Oncotarget.

[30]  Eva Bezak,et al.  A review of the development of tumor vasculature and its effects on the tumor microenvironment , 2017, Hypoxia.

[31]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[32]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..