Key pathways involved in prostate cancer based on gene set enrichment analysis and meta analysis

Prostate cancer is one of the most common male malignant neoplasms; however, its causes are not completely understood. A few recent studies have used gene expression profiling of prostate cancer to identify differentially expressed genes and possible relevant pathways. However, few studies have examined the genetic mechanics of prostate cancer at the pathway level to search for such pathways. We used gene set enrichment analysis and a meta-analysis of six independent studies after standardized microarray preprocessing, which increased concordance between these gene datasets. Based on gene set enrichment analysis, there were 12 down- and 25 up-regulated mixing pathways in more than two tissue datasets, while there were two down- and two up-regulated mixing pathways in three cell datasets. Based on the meta-analysis, there were 46 and nine common pathways in the tissue and cell datasets, respectively. Three up- and 10 down-regulated crossing pathways were detected with combined gene set enrichment analysis and meta-analysis. We found that genes with small changes are difficult to detect by classic univariate statistics; they can more easily be identified by pathway analysis. After standardized microarray preprocessing, we applied gene set enrichment analysis and a meta-analysis to increase the concordance in identifying biological mechanisms involved in prostate cancer. The gene pathways that we identified could provide insight concerning the development of prostate cancer.     

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.     

A network-based gene expression signature informs prognosis and treatment for colorectal cancer patients

Background

Several studies have reported gene expression signatures that predict recurrence risk in stage II and III colorectal cancer (CRC) patients with minimal gene membership overlap and undefined biological relevance. The goal of this study was to investigate biological themes underlying these signatures, to infer genes of potential mechanistic importance to the CRC recurrence phenotype and to test whether accurate prognostic models can be developed using mechanistically important genes.

Methods and Findings

We investigated eight published CRC gene expression signatures and found no functional convergence in Gene Ontology enrichment analysis. Using a random walk-based approach, we integrated these signatures and publicly available somatic mutation data on a protein-protein interaction network and inferred 487 genes that were plausible candidate molecular underpinnings for the CRC recurrence phenotype. We named the list of 487 genes a NEM signature because it integrated information from Network, Expression, and Mutation. The signature showed significant enrichment in four biological processes closely related to cancer pathophysiology and provided good coverage of known oncogenes, tumor suppressors, and CRC-related signaling pathways. A NEM signature-based Survival Support Vector Machine prognostic model was trained using a microarray gene expression dataset and tested on an independent dataset. The model-based scores showed a 75.7% concordance with the real survival data and separated patients into two groups with significantly different relapse-free survival (p = 0.002). Similar results were obtained with reversed training and testing datasets (p = 0.007). Furthermore, adjuvant chemotherapy was significantly associated with prolonged survival of the high-risk patients (p = 0.006), but not beneficial to the low-risk patients (p = 0.491).

Conclusions

The NEM signature not only reflects CRC biology but also informs patient prognosis and treatment response. Thus, the network-based data integration method provides a convergence between biological relevance and clinical usefulness in gene signature development.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.     

A cross-study gene set enrichment analysis identifies critical pathways in endometriosis

Background  

Endometriosis is an enigmatic disease. Gene expression profiling of endometriosis has been used in several studies, but few studies went further to classify subtypes of endometriosis based on expression patterns and to identify possible pathways involved in endometriosis. Some of the observed pathways are more inconsistent between the studies, and these candidate pathways presumably only represent a fraction of the pathways involved in endometriosis.     

GEPIS--quantitative gene expression profiling in normal and cancer tissues

MOTIVATION: Expression profiling in diverse tissues is fundamental to understanding gene function as well as therapeutic target identification. The vast collection of expressed sequence tags (ESTs) and the associated tissue source information provides an attractive opportunity for studying gene expression. RESULTS: To facilitate EST-based expression analysis, we developed GEPIS (gene expression profiling in silico), a tool that integrates EST and tissue source information to compute gene expression patterns in a large panel of normal and tumor samples. We found EST-based expression patterns to be consistent with published papers as well as our own experimental results. We also built a GEPIS Regional Atlas that depicts expression characteristics of all genes in a selected genomic region. This program can be adapted for large-scale screening for genes with desirable expression patterns, as illustrated by our large-scale mining for tissue- and tumor-specific genes. AVAILABILITY: The email server version of the GEPIS application is freely available at http://share.gene.com/share/gepis. An interactive version of GEPIS will soon be freely available at http://www.cgl.ucsf.edu/Research/genentech/gepis/. The source code, modules, data and gene lists can be downloaded at http://share.gene.com/share/gepis.     

Significance analysis of groups of genes in expression profiling studies

MOTIVATION: Gene class testing (GCT) is a statistical approach to determine whether some functionally predefined classes of genes express differently under two experimental conditions. GCT computes the P-value of each gene class based on the null distribution and the gene classes are ranked for importance in accordance with their P-values. Currently, two null hypotheses have been considered: the Q1 hypothesis tests the relative strength of association with the phenotypes among the gene classes, and the Q2 hypothesis assesses the statistical significance. These two hypotheses are related but not equivalent. METHOD: We investigate three one-sided and two two-sided test statistics under Q1 and Q2. The null distributions of gene classes under Q1 are generated by permuting gene labels and the null distributions under Q2 are generated by permuting samples. RESULTS: We applied the five statistics to a diabetes dataset with 143 gene classes and to a breast cancer dataset with 508 GO (Gene Ontology) terms. In each statistic, the null distributions of the gene classes under Q1 are different from those under Q2 in both datasets, and their rankings can be different too. We clarify the one-sided and two-sided hypotheses, and discuss some issues regarding the Q1 and Q2 hypotheses for gene class ranking in the GCT. Because Q1 does not deal with correlations among genes, we prefer test based on Q2. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Comprehensive gene expression analysis reveals multiple signal pathways associated with prostate cancer

Prostate cancer (PC) depends on androgenic signaling for growth and survival. To data, the exact molecular mechanism of hormone controlling proliferation and tumorigenesis in the PC remains unclear. Therefore, in this study, we explored the differentially expressed genes (DEGs) and identified featured genes related to hormone stimulus from PC. Two sets of gene expression data, including PC and normal control sample, were downloaded from Gene Expression Omnibus (GEO) database. The t-test was used to identify DEGs between PC and controls. Gene ontology (GO) functional annotation was applied to analyze the function of DEGs and screen hormone-related DEGs. Then these hormone-related DEGs were further analyzed in constructed cancer network and Human Protein Reference Database to screen important signaling pathways they participated in. A total of 912 DEGs were obtained which included 326 up-regulated genes and 586 down-regulated genes. GO functional enrichment analysis identified 50 hormone-related DEGs associated with PC. After pathway and PPI network analysis, we found these hormone-related DEGs participated in several important signaling pathways including TGF-β (TGFB2, TGFB3 and TGFBR2), MAPK (TGFB2, TGFB3 and TGFBR2), insulin (PIK3R3, SHC1 and EIF4EBP1), and p53 signaling pathways (CCND2 and CDKN1A). In addition, a total of five hormone-related DEGs (SHC1, CAV1, RXRA, CDKN1A and SRF) were located in the center of PPI network and 12 hormone-related DEGs formed six protein modules. These important signal pathways and hormone-related DEGs may provide potential therapeutic targets for PC.     

Prognostic classification of breast cancer and gene expression profiling

Clinical and pathological heterogeneity of breast cancer, partly responsible of therapeutic failures, reflects complex and combinatory molecular alterations until now poorly documented by classical investigation tools. Thorough molecular typing is crucial. The advent of DNA microarray-based gene expression profiling allowed consistent progresses in this direction. A novel molecular taxonomy of breast cancer has been defined, signatures that predict clinical outcome or therapeutic response have been identified, some of them being tested in ongoing prospective clinical trials. In this review, we present the main results and their potential clinical applications. We also discuss their current limits and future hopes in the therapeutic management of patients.     

Detection of gene pathways with predictive power for breast cancer prognosis

Background  

Prognosis is of critical interest in breast cancer research. Biomedical studies suggest that genomic measurements may have independent predictive power for prognosis. Gene profiling studies have been conducted to search for predictive genomic measurements. Genes have the inherent pathway structure, where pathways are composed of multiple genes with coordinated functions. The goal of this study is to identify gene pathways with predictive power for breast cancer prognosis. Since our goal is fundamentally different from that of existing studies, a new pathway analysis method is proposed.     

Gene expression profiles and pathway enrichment analysis to identification of differentially expressed gene and signaling pathways in epithelial ovarian cancer based on high-throughput RNA-seq data

Epithelial ovarian cancer (EOC) can be considered as a stressful and challenging disease among all women in the world, which has been associated with a poor prognosis and its molecular pathogenesis has remained unclear. In recent years, RNA Sequencing (RNA-seq) has become a functional and amazing technology for profiling gene expression. In the present study, RNA-seq raw data from Sequence Read Archive (SRA) of six tumor and normal ovarian sample was extracted, and then analysis and statistical interpretation was done with Linux and R Packages from the open-source Bioconductor. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of key genes and pathways involved in EOC. We identified 1091 Differential Expression Genes (DEGs) which have been reported in various studies of ovarian cancer as well as other types of cancer. Among them, 333 genes were up-regulated and 273 genes were down-regulated. In addition, Differentially Expressed Genes (DEGs) including RPL41, ALDH3A2, ERBB2, MIEN1, RBM25, ATF4, UPF2, DDIT3, HOXB8 and IL17D as well as Ribosome and Glycolysis/Gluconeogenesis pathway have had the potentiality to be used as targets for EOC diagnosis and treatment. In this study, unlike that of any other studies on various cancers, ALDH3A2 was most down-regulated gene in most KEGG pathways, and ATF4 was most up-regulated gene in leucine zipper domain binding term. In the other hand, RPL41 as a regulatory of cellular ATF4 level was up-regulated in many term and pathways and augmentation of ATF4 could justify the increase of RPL41 in the EOC. Pivotal pathways and significant genes, which were identified in the present study, can be used for adaptation of different EOC study. However, further molecular biological experiments and computational processes are required to confirm the function of the identified genes associated with EOC.     

Integrated gene expression profiling and linkage analysis in the rat

The combined application of genome-wide expression profiling from microarray experiments with genetic linkage analysis enables the mapping of expression quantitative trait loci (eQTLs) which are primary control points for gene expression across the genome. This approach allows for the dissection of primary and secondary genetic determinants of gene expression. The cis-acting eQTLs in practice are easier to investigate than the trans-regulated eQTLs because they are under simpler genetic control and are likely to be due to sequence variants within the gene itself or its neighboring regulatory elements. These genes are therefore candidates both for variation in gene expression and for contributions to whole-body phenotypes, particularly when these are located within known and relevant physiologic QTLs. Multiple trans-acting eQTLs tend to cluster to the same genetic location, implying shared regulatory control mechanisms that may be amenable to network analysis to identify gene clusters within the same metabolic pathway. Such clusters may ultimately underlie development of individual complex, whole-body phenotypes. The combined expression and linkage approach has been applied successfully in several mammalian species, including the rat which has specific features that demonstrate its value as a model for studying complex traits.     

Global gene expression profiling and cluster analysis in Xenopus laevis

We have undertaken a large-scale microarray gene expression analysis using cDNAs corresponding to 21,000 Xenopus laevis ESTs. mRNAs from 37 samples, including embryos and adult organs, were profiled. Cluster analysis of embryos of different stages was carried out and revealed expected affinities between gastrulae and neurulae, as well as between advanced neurulae and tadpoles, while egg and feeding larvae were clearly separated. Cluster analysis of adult organs showed some unexpected tissue-relatedness, e.g. kidney is more related to endodermal than to mesodermal tissues and the brain is separated from other neuroectodermal derivatives. Cluster analysis of genes revealed major phases of co-ordinate gene expression between egg and adult stages. During the maternal-early embryonic phase, genes maintaining a rapidly dividing cell state are predominantly expressed (cell cycle regulators, chromatin proteins). Genes involved in protein biosynthesis are progressively induced from mid-embryogenesis onwards. The larval-adult phase is characterised by expression of genes involved in metabolism and terminal differentiation. Thirteen potential synexpression groups were identified, which encompass components of diverse molecular processes or supra-molecular structures, including chromatin, RNA processing and nucleolar function, cell cycle, respiratory chain/Krebs cycle, protein biosynthesis, endoplasmic reticulum, vesicle transport, synaptic vesicle, microtubule, intermediate filament, epithelial proteins and collagen. Data filtering identified genes with potential stage-, region- and organ-specific expression. The dataset was assembled in the iChip microarray database, , which allows user-defined queries. The study provides insights into the higher order of vertebrate gene expression, identifies synexpression groups and marker genes, and makes predictions for the biological role of numerous uncharacterized genes.