Skip to content

viktormiok/tigaR

Repository files navigation

version GitHub code size in bytes GitHub issues

dependencies commit GitHub

Edit with Gitpod

tigaR

Overview

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 the estimates of the DNA copy number more stable, model parameters are also multivariate using triplets of features, imposing a spatial prior for the copy number effect.

Application

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.
  • Identify miRNA targets for 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.

Installation

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")

Docker

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

Data

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). They 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_*

Tutorials

Please see the following tutorials for detailed examples of how to use tigaR:

tigaR walkthrough:

License

tigaR is distributed under the MIT license. Please read the license before using tigaR, distributed in the LICENSE file.

References

Publications related to tigaR include:

Please cite the relevant publications if you use tigaR.