Dataset Information

TigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities.

ABSTRACT:

Background

To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies.

Results

Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect.

Conclusion

With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set.

SUBMITTER: Miok V

PROVIDER: S-EPMC4288633 | biostudies-literature | 2014 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities.

Miok Viktorian V Wilting Saskia M SM van de Wiel Mark A MA Jaspers Annelieke A van Noort Paula I PI Brakenhoff Ruud H RH Snijders Peter J F PJ Steenbergen Renske D M RD van Wieringen Wessel N WN

BMC bioinformatics 20141002

<h4>Background</h4>To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative ge ...[more]

PMID: 25278371

Similar Datasets

Project description:This study investigated the expression pattern of retinoblastoma binding protein 4 (RBBP4) gene in glioma and explored its associations with clinicopathologic characteristics and prognosis of patients. Data retrieved from the GEPIA, CGGA, HPA and TIMER databases were processed to analyze RBBP4 expression in glioma and investigate its relationship with clinicopathologic characteristics, tumor immune infiltration and prognosis in glioma patients. Immunohistochemistry was applied to determine the expression of RBBP4 protein in glioma tissue. Additionally, the Coexpedia database was visited to identify co-expressed genes for RBBP4 gene, while the Cytoscape software was run to visualize the enriched GO entries and KEGG pathways of these co-expressed genes. The expression levels of RBBP4 in lower-grade glioma (LGG) and glioblastoma (GBM) tissues were markedly elevated when compared to normal tissues (both p < 0.05). The up-regulation of RBBP4 expression was associated with an increase in WHO grade (II-IV), wild-type IDH, and 1p/19q non-codeletion (all p < 0.05). Multi-variate Cox regression analysis showed that both increased abundance of infiltrating macrophages and up-regulated RBBP4 expression independently predicted poor survival outcomes in LGG patients (both p < 0.05). Furthermore, RBBP4 expression exhibited significant positive correlations with the abundance of infiltrating B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil, and dendritic cell in LGG (all p < 0.05). Functional enrichment analyses indicated that the co-expressed genes associated with RBBP4 were highly involved in pathways such as the p53 signaling pathway, cell cycle, DNA replication, glutathione metabolism, as well as biological processes including cell cycle process, DNA replication, and DNA repair. High levels of RBBP4 are predictive for the poor survival outcome of LGG patients. RBBP4 gene, therefore, is expected to be a potential biomarker for prognosis of LGG and a target for immunotherapy. Highlights • The function and value of RBBP4 gene in prognosis of LGG have been less well studied by now.• Our study first used bioinformatics to reveal the expression pattern and prognostic value of RBBP4 in LGG.• Our study demonstrated that RBBP4 gene is promising as a target for clinical research and treatment of LGG.

Project description:BackgroundComparative Genomic Hybridization (CGH) is a molecular approach for detecting DNA Copy Number Alterations (CNAs) in tumor, which are among the key causes of tumorigenesis. However in the post-genomic era, most studies in cancer biology have been focusing on Gene Expression Profiling (GEP) but not CGH, and as a result, an enormous amount of GEP data had been accumulated in public databases for a wide variety of tumor types. We exploited this resource of GEP data to define possible recurrent CNAs in tumor. In addition, the CNAs identified by GEP would be more functionally relevant CNAs in the disease pathogenesis since the functional effects of CNAs can be reflected by altered gene expression.MethodsWe proposed a novel computational approach, coined virtual CGH (vCGH), which employs hidden Markov models (HMMs) to predict DNA CNAs from their corresponding GEP data. vCGH was first trained on the paired GEP and CGH data generated from a sufficient number of tumor samples, and then applied to the GEP data of a new tumor sample to predict its CNAs.ResultsUsing cross-validation on 190 Diffuse Large B-Cell Lymphomas (DLBCL), vCGH achieved 80% sensitivity, 90% specificity and 90% accuracy for CNA prediction. The majority of the recurrent regions defined by vCGH are concordant with the experimental CGH, including gains of 1q, 2p16-p14, 3q27-q29, 6p25-p21, 7, 11q, 12 and 18q21, and losses of 6q, 8p23-p21, 9p24-p21 and 17p13 in DLBCL. In addition, vCGH predicted some recurrent functional abnormalities which were not observed in CGH, including gains of 1p, 2q and 6q and losses of 1q, 6p and 8q. Among those novel loci, 1q, 6q and 8q were significantly associated with the clinical outcomes in the DLBCL patients (p < 0.05).ConclusionsWe developed a novel computational approach, vCGH, to predict genome-wide genetic abnormalities from GEP data in lymphomas. vCGH can be generally applied to other types of tumors and may significantly enhance the detection of functionally important genetic abnormalities in cancer research.