Weighted set enrichment of gene expression data
© Qureshi and Sacan; licensee BioMed Central Ltd. 2013
Published: 23 October 2013
Sets of genes that are known to be associated with each other can be used to interpret microarray data. This gene set approach to microarray data analysis can illustrate patterns of gene expression which may be more informative than analyzing the expression of individual genes. Various statistical approaches exist for the analysis of gene sets. There are three main classes of these methods: over-representation analysis, functional class scoring, and pathway topology based methods.
We propose weighted hypergeometric and weighted chi-squared methods in order to assign a rank to the degree to which each gene participates in the enrichment. Each gene is assigned a weight determined by the absolute value of its log fold change, which is then raised to a certain power. The power value can be adjusted as needed. Datasets from the Gene Expression Omnibus are used to test the method. The significantly enriched pathways are validated through searching the literature in order to determine their relevance to the dataset.
Although these methods detect fewer significantly enriched pathways, they can potentially produce more relevant results. Furthermore, we compare the results of different enrichment methods on a set of microarray studies all containing data from various rodent neuropathic pain models.
Our method is able to produce more consistent results than other methods when evaluated on similar datasets. It can also potentially detect relevant pathways that are not identified by the standard methods. However, the lack of biological ground truth makes validating the method difficult.
Due to their ability to provide comprehensive snapshots of cellular activity, microarrays have become a widely utilized tool in bio-medical sciences. Microarray-based gene expression detection has been used for biomarker discovery as well as diagnostic and prognostic purposes [1–4]. Online microarray experiment repositories such as Gene Expression Omnibus (GEO) [5, 6], ArrayExpress , and Stanford Microarray Database (SMD)  are invaluable resources containing gene expression profiles that span multiple developmental stages, experimental conditions, and model organisms [9, 10]. There are numerous challenges presented by the expanding availability of microarray data. The difficulty of interpreting the lists of significant genes produced by microarray experiments is a major challenge. The staggering number and diversity of the differentially expressed genes can be hard to interpret in a biologically meaningful way. As a result several statistical methods for gene set enrichment have been developed. The set of differentially expressed genes is compared to gene sets from various databases including Gene Ontology (GO)  or the Kyoto Encyclopedia of Genes and Genomes (KEGG) .
Huang et al. reviewed and classified 68 available tools for the statistical analysis of gene sets . Huang classified the available pathway enrichment methods into three categories: over-representation analysis (ORA), functional class scoring (FCS), and pathway topology (PT) based methods [13, 14]. In ORA, a list of genes is compiled by selecting genes based on their significance, fold change, or both. ORA techniques seek to identify whether the gene list is over-represented in a gene set or pathway. In ORA approaches, if k genes from the list are found in a pathway then the probability of finding k or more genes is calculated. The resulting p-values are used to determine whether or not a pathway or gene set is significantly enriched. The probability can be calculated using the chi-squared distribution, Fisher's Exact Test, the binomial probability distribution, or the hypergeometric distribution .
In functional class scoring approaches, such as Gene Set Enrichment Analysis (GSEA) , all genes are considered when calculating enrichment instead of a pre-selected list [14, 15]. This can deliver improved statistical power . In the FCS approaches, genes are assigned ranks. In GSEA, a gene's rank is determined by its correlation with the experimental sample classifications. When calculating the significance of a gene set, the null hypothesis is that the genes in a set are randomly distributed throughout the ranked list of genes from the microarray experiment. GSEA creates a null distribution by randomly permuting the labels of the samples and producing lists of genes ranked by their correlation with the newly shuffled sample labels. Using this null distribution to estimate the significance is analogous to a weighted Kolmogorov-Smirnov-like statistic . In contrast, Parametric Analysis of Gene Set Enrichment (PAGE) determines a z-score for a set and uses normal distribution to determine significance .
Both ORA and FCS approaches ignore the connections between genes in a pathway, however PT-based approaches integrate the information contain in the edges of a pathway when determining the enrichment. The disadvantage of PT-based approaches is that they cannot be applied to the Gene Ontology . ScorePAGE computes a score that represents the similarity between pairs of genes, and then divides this score by the number of edges between the two genes [14, 17]. Another approach is that of Signaling Pathway Impact Analysis (SPIA), which computes a "perturbation factor" for each gene in a pathway. This is given by the change in expression of the gene and by a linear function of the perturbation factors of all the other genes in the pathway. The "impact factor" of the pathway is a statistic calculated by taking the sum of the perturbation factors of the genes in the pathway [14, 18].
We have previously proposed a method for enriching gene sets that is a hybrid of over-representation and functional class scoring . Our method requires the contribution of all the genes in the dataset. Each gene contributes to the enrichment score in proportion to its fold change. Like ORA we calculate significance using the hypergeometric or chi-square distribution; however our method weighs the probability calculation by the fold change of the genes. In our method each gene is assigned a score based on its fold change, and we create a pseudo pathway, which is proportionally larger than the original pathway. We then calculate the significance of sampling the sum of the scores of the genes from the larger pseudo pathway.
We applied this pathway enrichment methodology in order to perform a meta-analysis of rodent neuropathic pain microarray experiments. Neuropathic pain is a chronic condition resulting from damage to any part of the nervous system or from diseases affecting an area of the nervous system. Neuropathic pain is typically accompanied by inflammation  and sensory and motor dysfunction . Up to eight percent of the general population is affected by neuropathic pain [22, 23]. While there is no clear etiology for neuropathic pain, spinal cord injury, diabetes, alcoholism, chemotherapy, chronic viral infection, transverse myelitis, and strokes are common causes. Due to the complex etiology and symptoms of neuropathic pain and its poorly understood mechanisms, the classification of chronic pain syndromes has remained largely subjective. Common treatments are able to produce better than moderate pain relief in only one third of patients . Treatments such as opiates, tricyclic antidepressants, anti-convulsants, anti-epileptics, topical analgesics, and NMDA-antagonists are used despite their limited efficacy and harmful side-effects [25, 26]. There are several rodent models of neuropathic pain such as nerve ligation, chronic constriction, and spared nerve injury [27, 28].
Gene ID mapping
Before we could begin performing enrichment analysis, we needed to construct a back-end database containing relevant information from various databases. Towards this end, we stored all the KEGG pathways and the genes involved in each pathway in a database. We further created a database to map the correspondence of Entrez genes with Affymetrix probe identifiers to enable gene identifier conversion. The correspondence between Entrez gene identifiers and KEGG gene identifiers was also mapped. Ultimately, a database for the KEGG pathway information and a database for gene identifier conversion were created in SQLite. Only genes from Homo Sapiens, Rattus Norvegicus, and Mus Musculus were included in the database. We mapped Affymetrix gene identifiers to Entrez Gene identifiers for Affymetrix microarray datasets with binary classifications obtained from the Gene Expression Omnibus (GEO) . Because multiple Affymetrix probes can map to a single Entrez Gene, we took the mean of the fold-change of the corresponding probes and the minimum of their p-values.
where N is the number of genes on the array, m is the number of significant genes, n is the number of genes in the particular KEGG pathway, and k is the number of genes that are both significant and present in the particular KEGG pathway. Thus we were able to calculate the significance of the enrichment of the KEGG pathways and rank them by their significance.
where all variables represent the same quantities that they do in equation 1, and all quantities are rounded to the nearest whole number. We ranked the pathways using this p-value.
The 2x2 table used to calculate the chi-squared statistic.
Genes on Array
N1r = n11 + n12
Not in Pathway
N2r = n21 + n22
N1c = n11+n21
N2c = n12+n22
N = n11+n12+n21+n22
where r is the number of rows in the table and c is the number of columns. We compute a weighted chi-squared statistic by constructing a table similar to Table 1, but in the place of the significant genes column we use the pathway score calculated by Equation 5 and the sum of the scores of all the genes on the array. Unlike the hypergeometric probability distribution, the chi-squared probability distribution is continuous. We did not need to discretize our data.
Experiments and results
The results of hypergeometric enrichment of the genes that are significant at the 0
# of significant genes
Natural killer cell mediated cytotoxicity
Cytokine-cytokine receptor interaction
ErbB signaling pathway
Non-small cell lung cancer
Vibrio cholerae infection
The results of weighted hypergeometric enrichment of the C.Pneumonia infection dataset
Glycosphingolipid biosynthesis - globo series
Glycosphingolipid biosynthesis - ganglio series
Pantothenate and CoA biosynthesis
D-Arginine and D-ornithine metabolism
Vitamin digestion and absorption
Primary bile acid biosynthesis
Ether lipid metabolism
The results of weighted chi-squared enrichment of the C. Pneumonia infection dataset
Drug metabolism - cytochrome P450
Amyotrophic lateral sclerosis (ALS)
Pathways in cancer
Cell adhesion molecules (CAMs)
Calcium signaling pathway
NOD-like receptor signaling pathway
Glycosphingolipid biosynthesis - ganglio series
Glyoxylate and dicarboxylate metabolism
Glycosphingolipid biosynthesis - globo series
The p-values and fold changes of the genes in the glycosphingolipid biosynthesis-globo pathway
Entrez Gene ID
We performed a literature search to assess the relevance of the top-10 pathways produced by the various methods. Since the dataset utilized involved the infection of cells, pathways related to the immune system should be enriched. Hypergeometric enrichment identified two potentially relevant pathways: natural killer cell mediated cytotoxicity and V. Cholerae infection. The top-2 pathways produced by this method were nitrogen metabolism and biotin metabolism. Weighted hypergeometric enrichment identified both types of glycosphingolipid biosynthesis as the top-2 pathways, with a glycosaminoglycan degradation related pathway as the third ranked pathway. Glycans and glycosylation are essential components of the antigen-presenting function of dendritic cells . Glycosphingolipids are proteins present in the plasma membrane that are known to be involved in immune function. They can act as cell-surface antigens [32, 33]. The standard method failed to detect these pathways; where as the weighted hypergeometric method uncovered the action of these pathways and helped elucidate mechanisms of the infection of the cells. In addition, despite finding fewer significant pathways the weighted chi-squared method also detected the glycosphingolipid synthesis pathways among its top-10 pathways, although at lower ranks than weighted hypergeometric enrichment. These pathways, which are the top-ranked pathways by the weighted hypergeometric method, are more biologically relevant than the top-ranked pathways generated by the standard hypergeometric method.
Weighted hypergeometric and chi-squared enrichment extend over-representation analysis to include change in expression of the genes and include all genes instead of a pre-selected list. These approaches enable every gene to contribute to the enrichment in proportion to their fold-change. Changing the power parameter enables one to adjust how much the expression change of the genes contributes to the enrichment score of the pathway. Our approach combines ORA and FCS methodologies. Unlike GSEA  our methods can detect pathways comprised of both up and down regulated genes by means of the score calculated for each gene. This is because we consider only the magnitude of expression change and not its direction with our score. There already exists a modification of GSEA that allows enrichment of pathways with bidirectional gene expression .
Weighted enrichment methods are much more conservative than unweighted methods. Because the weighted hypergeometric enrichment methods are so conservative, they produce no significant results when corrected for multiple comparisons. The Benjamini-Hochberg false discovery correction  was applied to the weighted enrichment, and the results are depicted in Tables 2, 3, 4. Table 2 shows that after false discovery rate (FDR) correction there are no significant pathways (p < 0.05), and that each pathway has an FDR-corrected p-value of 1. However, the FDR correction is not suitable for application to enrichment analysis because FDR has a high variability and should be applied to a larger number of p-values than those generated by enrichment . Furthermore, it has been shown that most multiple comparison corrections decrease the power of the analysis and are also too conservative [13, 42]. Additionally, the p-values resulting from enrichment analyses can be fragile and sensitive to non-statistical aspects of their calculation such as the data sources or the mapping of gene names between different conventions; these issues cannot be resolved by correction for multiple comparisons . Huang et al. advise using prior biological knowledge to assess the enriched pathways, and that the results of enrichment should only be guidelines for an investigator . Thus, FDR values were only included to be thorough when describing the results of these methods. We advise considering only the top-10 pathways instead of multiple comparison correction. Furthermore, we evaluated the consistency of our methods by considering the top-10 pathways enriched in data from several similar experiments. We were able to demonstrate that weighted hypergeometric enrichment produced the most consistent results.
Validating the weighted enrichment methods has proved to be challenging because there is no ground truth to compare the enriched pathways against. As a result, validation of the methods was performed based on literature search, which is not a complete or objective analysis. Literature search-based validation is biased towards already known pathways. There is no way of knowing whether pathways enriched by the dataset that have not been previously identified in the literature are actually associated with the disease or are falsely identified as enriched. Furthermore, our method is still sensitive to the handling of gene identifier mapping. Another drawback of this methodology is that it ignores the topology of the pathways. It is possible, for example, that an increase in the expression of a gene could be canceled out by a decrease in the expression of a downstream gene that is up regulated by the first gene. There is no way to address this situation when using our method. However, this drawback does not apply to the enrichment of Gene Ontology terms, which are arranged hierarchically.
We have proposed weighted hypergeometric and chi-squared methods to enrich gene sets. These methods can produce more biologically relevant results for KEGG pathway enrichment than the standard hypergeometric approach, despite the fact that the problem of Type II errors is inadequately addressed by correcting for multiple comparisons. We also showed that our method tends to produce more consistent results when using data from similar experiments. Despite only showing the results of KEGG pathway enrichment, these methods can also be applied to the Gene Ontology classifications as well as any other set of genes.
The publication costs for this article were funded by the corresponding author.
This article has been published as part of BMC Systems Biology Volume 7 Supplement 4, 2013: Selected articles from the IEEE International Conference on Bioinformatics and Biomedicine 2012: Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/7/S4.
- Ntzani EE, Ioannidis JPA: Predictive ability of DNA microarrays for cancer outcomes and correlates: an empirical assessment. Lancet. 2003, 362 (9394): 1439-1444. 10.1016/S0140-6736(03)14686-7.View ArticlePubMedGoogle Scholar
- van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AAM, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.View ArticlePubMedGoogle Scholar
- Pomeroy SL, Tamayo P, Gaasenbeek M, Sturla LM, Angelo M, McLaughlin ME, Kim JYH, Goumnerova LC, Black PM, Lau C, et al: Prediction of central nervous system embryonal tumour outcome based on gene expression. Nature. 2002, 415 (6870): 436-442. 10.1038/415436a.View ArticlePubMedGoogle Scholar
- Simon R: Diagnostic and prognostic prediction using gene expression profiles in high-dimensional microarray data. Br J Cancer. 2003, 89 (9): 1599-1604. 10.1038/sj.bjc.6601326.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, et al: NCBI GEO: archive for functional genomics data sets--10 years on. Nucleic Acids Research. 2011, 39 (suppl 1): D1005-D1010.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Research. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M, et al: ArrayExpress--a public database of microarray experiments and gene expression profiles. Nucleic Acids Research. 2007, 35 (suppl 1): D747-D750.PubMed CentralView ArticlePubMedGoogle Scholar
- Sherlock G, Hernandez-Boussard T, Kasarskis A, Binkley G, Matese JC, Dwight SS, Kaloper M, Weng S, Jin H, Ball CA, et al: The Stanford Microarray Database. Nucleic Acids Research. 2001, 29 (1): 152-155. 10.1093/nar/29.1.152.PubMed CentralView ArticlePubMedGoogle Scholar
- Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101 (25): 9309-9314. 10.1073/pnas.0401994101.PubMed CentralView ArticlePubMedGoogle Scholar
- Dawany N, Tozeren A: Asymmetric microarray data produces gene lists highly predictive of research literature on multiple cancer types. BMC Bioinformatics. 2010, 11 (1): 483-10.1186/1471-2105-11-483.PubMed CentralView ArticlePubMedGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, et al: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Research. 2004, 32 (Database): D258-261.PubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Research. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.PubMed CentralView ArticleGoogle Scholar
- Khatri P, Sirota M, Butte AJ: Ten Years of Pathway Analysis: Current Approaches and Outstanding Challenges. PLoS Comput Biol. 2012, 8 (2): e1002375-10.1371/journal.pcbi.1002375.PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim S-Y, Volsky D: PAGE: Parametric Analysis of Gene Set Enrichment. BMC Bioinformatics. 2005, 6 (1): 144-10.1186/1471-2105-6-144.PubMed CentralView ArticlePubMedGoogle Scholar
- Rahnenfuhrer J, Domingues FS, Maydt J, Lengauer T: Calculating the statistical significance of changes in pathway activity from gene expression data. Stat Appl Genet Mol Biol. 2004, 3: Article16-PubMedGoogle Scholar
- Tarca AL, Draghici S, Khatri P, Hassan SS, Mittal P, Kim J-s, Kim CJ, Kusanovic JP, Romero R: A novel signaling pathway impact analysis. Bioinformatics. 2009, 25 (1): 75-82. 10.1093/bioinformatics/btn577.PubMed CentralView ArticlePubMedGoogle Scholar
- Qureshi R, Sacan A: A weighted hypergeometric statistic for the enrichment of gene sets. Bioinformatics and Biomedicine (BIBM), 2012 IEEE International Conference on: 4-7 Oct 2012. 2012, 1-6. 10.1109/BIBM.2012.6392685.View ArticleGoogle Scholar
- Merskey H, Bogduk N: International Association for the Study of Pain. Task Force on Taxonomy Classification of chronic pain : descriptions of chronic pain syndromes and definitions of pain terms, 2nd edn. 1994, Seattle: IASP PressGoogle Scholar
- Rasmussen PV, Sindrup SH, Jensen TS, Bach FW: Symptoms and signs in patients with suspected neuropathic pain. Pain. 2004, 110 (1-2): 461-469. 10.1016/j.pain.2004.04.034.View ArticlePubMedGoogle Scholar
- Torrance N, Smith BH, Bennett MI, Lee AJ: The Epidemiology of Chronic Pain of Predominantly Neuropathic Origin. Results From a General Population Survey. The Journal of Pain. 2006, 7 (4): 281-289. 10.1016/j.jpain.2005.11.008.View ArticlePubMedGoogle Scholar
- Bouhassira D, Lantéri-Minet M, Attal N, Laurent B, Touboul C: Prevalence of chronic pain with neuropathic characteristics in the general population. PAIN. 2008, 136 (3): 380-387. 10.1016/j.pain.2007.08.013.View ArticlePubMedGoogle Scholar
- Sindrup SH, Jensen TS: Efficacy of pharmacological treatments of neuropathic pain: an update and effect related to mechanism of drug action. Pain. 1999, 83 (3): 389-400. 10.1016/S0304-3959(99)00154-2.View ArticlePubMedGoogle Scholar
- Barclay J, Clark AK, Ganju P, Gentry C, Patel S, Wotherspoon G, Buxton F, Song C, Ullah J, Winter J, et al: Role of the cysteine protease cathepsin S in neuropathic hyperalgesia. PAIN. 2007, 130 (3): 225-234. 10.1016/j.pain.2006.11.017.View ArticlePubMedGoogle Scholar
- Costigan M, Belfer I, Griffin RS, Dai F, Barrett LB, Coppola G, Wu T, Kiselycznyk C, Poddar M, Lu Y, et al: Multiple chronic pain states are associated with a common amino acid-changing allele in KCNS1. Brain. 2010, 133 (9): 2519-2527. 10.1093/brain/awq195.PubMed CentralView ArticlePubMedGoogle Scholar
- Decosterd I, Woolf CJ: Spared nerve injury: an animal model of persistent peripheral neuropathic pain. Pain. 2000, 87 (2): 149-158. 10.1016/S0304-3959(00)00276-1.View ArticlePubMedGoogle Scholar
- Kim SH, Chung JM: An experimental model for peripheral neuropathy produced by segmental spinal nerve ligation in the rat. Pain. 1992, 50 (3): 355-363. 10.1016/0304-3959(92)90041-9.View ArticlePubMedGoogle Scholar
- Drǎghici S, Khatri P, Martins RP, Ostermeier GC, Krawetz SA: Global functional profiling of gene expression. Genomics. 2003, 81 (2): 98-104. 10.1016/S0888-7543(02)00021-6.View ArticlePubMedGoogle Scholar
- Njau F, Geffers R, Thalmann J, Haller H, Wagner AD: Restriction of Chlamydia pneumoniae replication in human dendritic cell by activation of indoleamine 2,3-dioxygenase. Microbes and Infection. 2009, 11 (13): 1002-1010. 10.1016/j.micinf.2009.07.006.View ArticlePubMedGoogle Scholar
- Erbacher A, Gieseke F, Handgretinger R, Müller I: Dendritic cells: Functional aspects of glycosylation and lectins. Human Immunology. 2009, 70 (5): 308-312. 10.1016/j.humimm.2009.02.005.View ArticlePubMedGoogle Scholar
- Ichikawa S, Hirabayashi Y: Glucosylceramide synthase and glycosphingolipid synthesis. Trends in cell biology. 1998, 8 (5): 198-202. 10.1016/S0962-8924(98)01249-5.View ArticlePubMedGoogle Scholar
- Uemura A, Watarai S, Iwasaki T, Kodama H: Induction of Immune Responses against Glycosphingolipid Antigens: Comparison of Antibody Responses in Mice Immunized with Antigen Associated with Liposomes Prepared from Various Phospholipids. Journal of Veterinary Medical Science. 2005, 67 (12): 1197-1201. 10.1292/jvms.67.1197.View ArticlePubMedGoogle Scholar
- Levin ME, Jin JG, Ji R-R, Tong J, Pomonis JD, Lavery DJ, Miller SW, Chiang LW: Complement activation in the peripheral nervous system following the spinal nerve ligation model of neuropathic pain. PAIN. 2008, 137 (1): 182-201. 10.1016/j.pain.2007.11.005.View ArticlePubMedGoogle Scholar
- von Schack D, Agostino MJ, Murray BS, Li Y, Reddy PS, Chen J, Choe SE, Strassle BW, Li C, Bates B, et al: Dynamic Changes in the MicroRNA Expression Profile Reveal Multiple Regulatory Mechanisms in the Spinal Nerve Ligation Model of Neuropathic Pain. PLoS ONE. 2011, 6 (3): e17670-10.1371/journal.pone.0017670.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou Y, Zhou Z-S, Zhao Z-Q: Neomycin blocks capsaicin-evoked responses in rat dorsal root ganglion neurons. Neuroscience Letters. 2001, 315 (1-2): 98-102. 10.1016/S0304-3940(01)02356-4.View ArticlePubMedGoogle Scholar
- Van Den Wijngaard RM, Welting O, Bulmer DC, Wouters MM, Lee K, De Jonge WJ, Boeckxstaens GE: Possible role for TRPV1 in neomycin-induced inhibition of visceral hypersensitivity in rat. Neurogastroenterology & Motility. 2009, 21 (8): 863-e860. 10.1111/j.1365-2982.2009.01287.x.View ArticleGoogle Scholar
- Materazzi S, Nassini R, Andrè E, Campi B, Amadesi S, Trevisani M, Bunnett NW, Patacchini R, Geppetti P: Cox-dependent fatty acid metabolites cause pain through activation of the irritant receptor TRPA1. Proceedings of the National Academy of Sciences. 2008, 105 (33): 12045-12050. 10.1073/pnas.0802354105.View ArticleGoogle Scholar
- Saxena V, Orgill D, Kohane I: Absolute enrichment: gene set enrichment analysis for homeostatic systems. Nucleic Acids Research. 2006, 34 (22): e151-10.1093/nar/gkl766.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995, 57 (1): 289-300.Google Scholar
- Gold DL, Miecznikowski JC, Liu S: Error control variability in pathway-based microarray analysis. Bioinformatics. 2009, 25 (17): 2216-2221. 10.1093/bioinformatics/btp385.PubMed CentralView ArticlePubMedGoogle Scholar
- Bluthgen N, Brand K, Cajavec B, Swat M, Herzel H, Beule D: Biological profiling of gene groups utilizing Gene Ontology. Genome informatics International Conference on Genome Informatics. 2005, 16 (1): 106-115.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.