Gene expression signature of human HepG2 cell line

Hepatocellular carcinoma (HCC) is the fifth most common cancer worldwide and is associated with various clinico-pathological characteristics such as genetic mutations and viral infections. Therefore, numerous laboratories look out for identifying always new putative markers for the improvement of HCC diagnosis/prognosis. Many molecular profiling studies investigated gene expression changes related to HCC. HepG2 represents a pure cell line of human liver carcinoma, often used as HCC model due to the absence of viral infection. In this study we compare gene expression profiles associated with HepG2 (as HCC model) and normal hepatocyte cells by microarray technology. Hierarchical cluster analysis of genes evidenced that 2646 genes significantly down-regulated in HepG2 cells compared to hepatocytes whereas a further 3586 genes significantly up-regulated. By using the Ingenuity Pathway Analysis (IPA) program, we have classified the genes that were differently expressed and studied the functional networks correlating these genes in the complete human interactome. Moreover, to confirm the differentially expressed genes as well as the reliability of our microarray data, we performed a quantitative Real time RT-PCR analysis on 9 up-regulated and 11 down-regulated genes, respectively. In conclusion this work i) provides a gene signature of human hepatoma cells showing genes that change their expression as a consequence of liver cancer in the absence of any genetic mutations or viral infection, ii) evidences new differently expressed genes found in our signature compared to previous published studies and iii) suggests some genes on which to focus future studies to understand if they can be used to improve the HCC prognosis/diagnosis.     

Feature genes of hepatitis B virus-positive hepatocellular carcinoma, established by its molecular discrimination approach using prediction analysis of microarray

Recent introduction of a learning algorithm for cDNA microarray analysis has permitted to select feature set to accurately distinguish human cancers according to their pathological judgments. Here, we demonstrate that hepatitis B virus-positive hepatocellular carcinoma (HCC) could successfully be identified from non-tumor liver tissues by supervised learning analysis of gene expression profiling. Through learning and cross-validating HCC sample set, we could identify an optimized set of 44 genes to discriminate the status of HCC from non-tumor liver tissues. In an analysis of other blind-tested HCC sample sets, this feature set was found to be statistically significant, indicating the reproducibility of our molecular discrimination approach with the defined genes. One prominent finding was an asymmetrical distribution pattern of expression profiling in HCC, in which the number of down-regulated genes was greater than that of up-regulated genes. In conclusion, the present findings indicate that application of learning algorithm to HCC may establish a reliable feature set of genes to be useful for therapeutic target of HCC, and that the asymmetric expression pattern may emphasize the importance of suppressed genes in HCC.     

Gene expression profile analysis in human hepatocellular carcinoma by cDNA microarray

We performed gene expression profiling of normal and hepatocellular carcinoma (HCC) liver tissues using a high-density microarray that contained 3,063 human cDNA. The results of a microarray hybridization experiment from eight different HCC tissues were analyzed and classified by the Cluster program. Among these differentially-expressed genes, the galectin-3, serine/threonine kinase SGK, translation factor eIF-4A, -4B, -3, fibroblast growth factor receptor, and ribosomal protein L35A were up-regulated; the mRNAs of Nip3, decorin, and the insulin-like growth factor binding protein-3 were down-regulated in HCC. The differential expression of these genes was further confirmed by an RT-PCR analysis. In addition, our data suggest that the gene expression profile of HCC varies according to the histological types.     

Integrative genomics based identification of potential human hepatocarcinogenesis-associated cell cycle regulators: RHAMM as an example

DNA microarray has been widely used to examine gene expression profile of different human tumors. The information generated from microarray analysis usually represents the overall range of cancer-associated abnormality associated with gene regulation. In order to identify key regulatory genes involved in carcinogenesis of human cancer, hypothesis driven data mining of the microarray data plus experimental validation becomes a critical approach in the post-genome era. Here, we present an integrative genomic analysis of published microarray data and homolog gene database. Over 20,000 genes were examined to reveal 16 genes specific to vertebrates, cell cycle G2/M regulated, and overexpressed in human HCC. Using Affymetrix microarray analysis, we found that all 16 genes were up-regulated in human HCC. Among these 16 genes, we experimentally validated the up-regulation of receptor for hyaluronan-mediated motility (RHAMM) in different cell model systems. We first confirmed elevation of RHAMM in the G2/M phase of synchronized HeLa cells. We also found that RHAMM had an elevated level of expression in all the HCC samples we examined and it was induced during the G2/M phase of regenerating mouse hepatocytes after partial hepatectomy. Thus, the expression of RHAMM appears to be tightly regulated during mammalian cell cycle G2/M progression. The ectopic overexpression of RHAMM in 293T cells resulted in the accumulation of cells at G2/M phase. RHAMM-induced mitotic arrest of cells was predominantly in the prophase. Taken together, using an integrated functional genomic approach, we have uncovered a set of genes that may play specific roles in cell cycle progression and in HCC development. To elucidate the function of these genes in cell cycle regulation may shed light on the control mechanism of human HCC in the future.     

Empirical Bayes identification [correction of identication] of tumor progression genes from microarray data

The use of microarray data has become quite commonplace in medical and scientific experiments. We focus here on microarray data generated from cancer studies. It is potentially important for the discovery of biomarkers to identify genes whose expression levels correlate with tumor progression. In this article, we propose a simple procedure for the identification of such genes, which we term tumor progression genes. The first stage involves estimation based on the proportional odds model. At the second stage, we calculate two quantities: a q-value, and a shrinkage estimator of the test statistic is constructed to adjust for the multiple testing problem. The relationship between the proposed method with the false discovery rate is studied. The proposed methods are applied to data from a prostate cancer microarray study.     

A Bayesian mixture model for metaanalysis of microarray studies

The increased availability of microarray data has been calling for statistical methods to integrate findings across studies. A common goal of microarray analysis is to determine differentially expressed genes between two conditions, such as treatment vs control. A recent Bayesian metaanalysis model used a prior distribution for the mean log-expression ratios that was a mixture of two normal distributions. This model centered the prior distribution of differential expression at zero, and separated genes into two groups only: expressed and nonexpressed. Here, we introduce a Bayesian three-component truncated normal mixture prior model that more flexibly assigns prior distributions to the differentially expressed genes and produces three groups of genes: up and downregulated, and nonexpressed. We found in simulations of two and five studies that the three-component model outperformed the two-component model using three comparison measures. When analyzing biological data of Bacillus subtilis, we found that the three-component model discovered more genes and omitted fewer genes for the same levels of posterior probability of differential expression than the two-component model, and discovered more genes for fixed thresholds of Bayesian false discovery. We assumed that the data sets were produced from the same microarray platform and were prescaled.     

Discovery of differentially expressed genes related to histological subtype of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common human malignancies in the world. To identify the histological subtype-specific genes of HCC, we analyzed the gene expression profile of 10 HCC patients by means of cDNA microarray. We proposed a systematic approach for determining the discriminatory genes and revealing the biological phenomena of HCC with cDNA microarray data. First, normalization of cDNA microarray data was performed to reduce or minimize systematic variations. On the basis of the suitably normalized data, we identified specific genes involved in histological subtype of HCC. Two classification methods, Fisher's discriminant analysis (FDA) and support vector machine (SVM), were used to evaluate the reliability of the selected genes and discriminate the histological subtypes of HCC. This study may provide a clue for the needs of different chemotherapy and the reason for heterogeneity of the clinical responses according to histological subtypes.     

OpenDMAP: An open source, ontology-driven concept analysis engine, with applications to capturing knowledge regarding protein transport, protein interactions and cell-type-specific gene expression

Background

Microarray technology provides an efficient means for globally exploring physiological processes governed by the coordinated expression of multiple genes. However, identification of genes differentially expressed in microarray experiments is challenging because of their potentially high type I error rate. Methods for large-scale statistical analyses have been developed but most of them are applicable to two-sample or two-condition data.

Results

We developed a large-scale multiple-group F-test based method, named ranking analysis of F-statistics (RAF), which is an extension of ranking analysis of microarray data (RAM) for two-sample t-test. In this method, we proposed a novel random splitting approach to generate the null distribution instead of using permutation, which may not be appropriate for microarray data. We also implemented a two-simulation strategy to estimate the false discovery rate. Simulation results suggested that it has higher efficiency in finding differentially expressed genes among multiple classes at a lower false discovery rate than some commonly used methods. By applying our method to the experimental data, we found 107 genes having significantly differential expressions among 4 treatments at <0.7% FDR, of which 31 belong to the expressed sequence tags (ESTs), 76 are unique genes who have known functions in the brain or central nervous system and belong to six major functional groups.

Conclusion

Our method is suitable to identify differentially expressed genes among multiple groups, in particular, when sample size is small.     

Bayesian hierarchical modeling for time course microarray experiments

Time course microarray experiments designed to characterize the dynamic regulation of gene expression in biological systems are becoming increasingly important. One critical issue that arises when examining time course microarray data is the identification of genes that show different temporal expression patterns among biological conditions. Here we propose a Bayesian hierarchical model to incorporate important experimental factors and to account for correlated gene expression measurements over time and over different genes. A new gene selection algorithm is also presented with the model to simultaneously identify genes that show changes in expression among biological conditions, in response to time and other experimental factors of interest. The algorithm performs well in terms of the false positive and false negative rates in simulation studies. The methodology is applied to a mouse model time course experiment to correlate temporal changes in azoxymethane-induced gene expression profiles with colorectal cancer susceptibility.     

Sex genes for genomic analysis in human brain: internal controls for comparison of probe level data extraction.

Background

Genomic studies of complex tissues pose unique analytical challenges for assessment of data quality, performance of statistical methods used for data extraction, and detection of differentially expressed genes. Ideally, to assess the accuracy of gene expression analysis methods, one needs a set of genes which are known to be differentially expressed in the samples and which can be used as a "gold standard". We introduce the idea of using sex-chromosome genes as an alternative to spiked-in control genes or simulations for assessment of microarray data and analysis methods.

Results

Expression of sex-chromosome genes were used as true internal biological controls to compare alternate probe-level data extraction algorithms (Microarray Suite 5.0 [MAS5.0], Model Based Expression Index [MBEI] and Robust Multi-array Average [RMA]), to assess microarray data quality and to establish some statistical guidelines for analyzing large-scale gene expression. These approaches were implemented on a large new dataset of human brain samples. RMA-generated gene expression values were markedly less variable and more reliable than MAS5.0 and MBEI-derived values. A statistical technique controlling the false discovery rate was applied to adjust for multiple testing, as an alternative to the Bonferroni method, and showed no evidence of false negative results. Fourteen probesets, representing nine Y- and two X-chromosome linked genes, displayed significant sex differences in brain prefrontal cortex gene expression.

Conclusion

In this study, we have demonstrated the use of sex genes as true biological internal controls for genomic analysis of complex tissues, and suggested analytical guidelines for testing alternate oligonucleotide microarray data extraction protocols and for adjusting multiple statistical analysis of differentially expressed genes. Our results also provided evidence for sex differences in gene expression in the brain prefrontal cortex, supporting the notion of a putative direct role of sex-chromosome genes in differentiation and maintenance of sexual dimorphism of the central nervous system. Importantly, these analytical approaches are applicable to all microarray studies that include male and female human or animal subjects.

     

A systems biology-based classifier for hepatocellular carcinoma diagnosis

Aim

The diagnosis of hepatocellular carcinoma (HCC) in the early stage is crucial to the application of curative treatments which are the only hope for increasing the life expectancy of patients. Recently, several large-scale studies have shed light on this problem through analysis of gene expression profiles to identify markers correlated with HCC progression. However, those marker sets shared few genes in common and were poorly validated using independent data. Therefore, we developed a systems biology based classifier by combining the differential gene expression with topological features of human protein interaction networks to enhance the ability of HCC diagnosis.

Methods and Results

In the Oncomine platform, genes differentially expressed in HCC tissues relative to their corresponding normal tissues were filtered by a corrected Q value cut-off and Concept filters. The identified genes that are common to different microarray datasets were chosen as the candidate markers. Then, their networks were analyzed by GeneGO Meta-Core software and the hub genes were chosen. After that, an HCC diagnostic classifier was constructed by Partial Least Squares modeling based on the microarray gene expression data of the hub genes. Validations of diagnostic performance showed that this classifier had high predictive accuracy (85.88∼92.71%) and area under ROC curve (approximating 1.0), and that the network topological features integrated into this classifier contribute greatly to improving the predictive performance. Furthermore, it has been demonstrated that this modeling strategy is not only applicable to HCC, but also to other cancers.

Conclusion

Our analysis suggests that the systems biology-based classifier that combines the differential gene expression and topological features of human protein interaction network may enhance the diagnostic performance of HCC classifier.     

Application of the S-score algorithm for analysis of oligonucleotide microarrays

In the past several years, oligonucleotide microarrays have emerged as a widely used tool for the simultaneous, non-biased measurement of expression levels for thousands of genes. Several challenges exist in successfully utilizing this biotechnology; principal among these is analysis of microarray data. An experiment to measure differential gene expression can consist of a dozen microarrays, each consisting of over a hundred thousand data points. Previously, we have described the use of a novel algorithm for analyzing oligonucleotide microarrays and assessing changes in gene expression. This algorithm describes changes in expression in terms of the statistical significance (S-score) of change, which combines signals detected by multiple probe pairs according to an error model characteristic of oligonucleotide arrays. Software is available that simplifies the use of the application of this algorithm so that it may be applied to improving the analysis of oligonucleotide microarray data. The application of this method to problems of the central nervous system is discussed.     

Gene expression analyses of subchondral bone in early experimental osteoarthritis by microarray

Osteoarthritis (OA) is a degenerative joint disease that affects both cartilage and bone. A better understanding of the early molecular changes in subchondral bone may help elucidate the pathogenesis of OA. We used microarray technology to investigate the time course of molecular changes in the subchondral bone in the early stages of experimental osteoarthritis in a rat model. We identified 2,234 differentially expressed (DE) genes at 1 week, 1,944 at 2 weeks and 1,517 at 4 weeks post-surgery. Further analyses of the dysregulated genes indicated that the events underlying subchondral bone remodeling occurred sequentially and in a time-dependent manner at the gene expression level. Some of the identified dysregulated genes that were identified have suspected roles in bone development or remodeling; these genes include Alp, Igf1, Tgf β1, Postn, Mmp3, Tnfsf11, Acp5, Bmp5, Aspn and Ihh. The differences in the expression of these genes were confirmed by real-time PCR, and the results indicated that our microarray data accurately reflected gene expression patterns characteristic of early OA. To validate the results of our microarray analysis at the protein level, immunohistochemistry staining was used to investigate the expression of Mmp3 and Aspn protein in tissue sections. These analyses indicate that Mmp3 protein expression completely matched the results of both the microarray and real-time PCR analyses; however, Aspn protein expression was not observed to differ at any time. In summary, our study demonstrated a simple method of separation of subchondral bone sample from the knee joint of rat, which can effectively avoid bone RNA degradation. These findings also revealed the gene expression profiles of subchondral bone in the rat OA model at multiple time points post-surgery and identified important DE genes with known or suspected roles in bone development or remodeling. These genes may be novel diagnostic markers or therapeutic targets for OA.     

A Regression-Based Differential Expression Detection Algorithm for Microarray Studies with Ultra-Low Sample Size

Global gene expression analysis using microarrays and, more recently, RNA-seq, has allowed investigators to understand biological processes at a system level. However, the identification of differentially expressed genes in experiments with small sample size, high dimensionality, and high variance remains challenging, limiting the usability of these tens of thousands of publicly available, and possibly many more unpublished, gene expression datasets. We propose a novel variable selection algorithm for ultra-low-n microarray studies using generalized linear model-based variable selection with a penalized binomial regression algorithm called penalized Euclidean distance (PED). Our method uses PED to build a classifier on the experimental data to rank genes by importance. In place of cross-validation, which is required by most similar methods but not reliable for experiments with small sample size, we use a simulation-based approach to additively build a list of differentially expressed genes from the rank-ordered list. Our simulation-based approach maintains a low false discovery rate while maximizing the number of differentially expressed genes identified, a feature critical for downstream pathway analysis. We apply our method to microarray data from an experiment perturbing the Notch signaling pathway in Xenopus laevis embryos. This dataset was chosen because it showed very little differential expression according to limma, a powerful and widely-used method for microarray analysis. Our method was able to detect a significant number of differentially expressed genes in this dataset and suggest future directions for investigation. Our method is easily adaptable for analysis of data from RNA-seq and other global expression experiments with low sample size and high dimensionality.     

Identification of differentially expressed gene categories in microarray studies using nonparametric multivariate analysis

MOTIVATION: The field of microarray data analysis is shifting emphasis from methods for identifying differentially expressed genes to methods for identifying differentially expressed gene categories. The latter approaches utilize a priori information about genes to group genes into categories and enhance the interpretation of experiments aimed at identifying expression differences across treatments. While almost all of the existing approaches for identifying differentially expressed gene categories are practically useful, they suffer from a variety of drawbacks. Perhaps most notably, many popular tools are based exclusively on gene-specific statistics that cannot detect many types of multivariate expression change. RESULTS: We have developed a nonparametric multivariate method for identifying gene categories whose multivariate expression distribution differs across two or more conditions. We illustrate our approach and compare its performance to several existing procedures via the analysis of a real data set and a unique data-based simulation study designed to capture the challenges and complexities of practical data analysis. We show that our method has good power for differentiating between differentially expressed and non-differentially expressed gene categories, and we utilize a resampling based strategy for controlling the false discovery rate when testing multiple categories. AVAILABILITY: R code (www.r-project.org) for implementing our approach is available from the first author by request.