The R-package tigaR performs temporal integrative differential expression analysis of omics data, both sequencings count RNAseq and continuous microarray data.
A generalized linear mixed model employing low-rank thin-plate splines describes the temporal variation in gene expression. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, the 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.
With the proposed R-package tigaR for analysis of time-course multilevel molecular continuous (microarray) and count (RNAseq) data, more profound insight may be gained through:
- identification of temporal differential gene expression, where the method yields improvements in sensitivity, specificity, and reproducibility compared to existing methods.
- employing the same spline in modeling up/down-regulated genes may be identified over the cell lines while using different splines more flexibility may be allowed in capturing the temporal variation over time.
- Identify temporal differential expression induced by DNA copy number abnormalities or/and miRNA expression levels.
- identification of miRNA targets to the mRNA gene expression in a time-course fashion.
Note: If you can choose to use Windows or Unix/Linux, opt for the latter. tigaR runs more efficiently under Unix/Linux than Windows. NOTE: When running tigaR you may see *** WARNINGS *** from INLA (e.g. on eigenvalues, or convergence, or even something like 18500 Aborted...). They can currently not be suppressed, because C-code produces them. Please ignore them.
The package tigaR depends on ShrinkBayes and on R >= 3.0.0 and is available from GitHub. This requires the package devtools:
devtools::install_github("viktormiok/tigaR", build_vignettes=TRUE)Please restart R before loading the package and its documentation:
library(tigaR)
utils::help(tigaR)
utils::vignette("tigaR")If your system configuration makes installing tigaR natively difficult, a docker container is an alternative way to get tigaR running.
Note: Docker Machine has Memory and CPU limits on Mac OS X. To control it, please check instructions either for CLI or for Docker Desktop.
For building a Docker image from the Dockerfile, download the Dockerfile (available in this repo) and run the following command to create it:
docker build -t tigaR .
This will create a tigaR docker image on your system (please be patient, as the build could take approximately 30-50 minutes to finish).
You can then run it using the following command:
docker run -d -p 8787:8787 -e PASSWORD=pass --name tigaR -it tigaR
All the data required for performing temporal integrative genomics analysis and published in the reference articles have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO) and are accessible through the GEO Series accession numbers:
| Data type | GEO number |
|---|---|
| CGH Arrays | GSE138724 |
| mRNA Arrays | GSE138079 |
| miRNA Arrays | GSE78279 |
To access one of the data sets for instance GSE78279 you need to run the code below. Unpacking the data requires tar and gunzip, which should already be available on most systems.
cd ../ #To get to the main GitHub repo folder
mkdir -p data/tigaR_data_analysis/
cd data/tigaR_data_analysis/
wget ftp://ftp.ncbi.nlm.nih.gov/geo/series/GSE78nnn/GSE78279/suppl/GSE78279_RAW.tar
mkdir GSE78279_RAW
tar -C GSE78279_RAW -xvf GSE78279_RAW.tar
gunzip GSE78279_RAW/*_Regional_*
Please see the following tutorials for detailed examples of how to use tigaR:
tigaR is distributed under the MIT license. Please read the license before using tigaR, distributed in the LICENSE file.
Publications related to tigaR include:
-
Miok, V., Wilting, S.M., Van de Wiel, M.A., Jaspers, A., van Noort, P.I., Brakenhoff, R.H., Snijders, P.J.F., Steenbergen, R.D.M., Van Wieringen, W.N. (2014), "tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities", BMC Bioinformatics, 15, 327.
-
Babion, I., Miok, V., Jaspers, A., Huseinovic, A., Steenbergen, R.D., van Wieringen, W.N., Wilting, S.M. (2020), "Identification of Deregulated Pathways, Key Regulators, and Novel miRNA-mRNA Interactions in HPV-Mediated Transformation", Cancers, 12(3), 700.
-
Babion, I., Miok, V., Jaspers, A., van Wieringen, W.N., Steenbergen, R.D., Wilting, S.M. (2018), "Comprehensive molecular profiling of HPV-induced transformation over time", Cancer Research, 22(3), 200.
Please cite the relevant publications if you use tigaR.