ImpulseDE: detection of differentially expressed genes in time series data using impulse models.
Ontology highlight
ABSTRACT: Summary
Perturbations in the environment lead to distinctive gene expression changes within a cell. Observed over time, those variations can be characterized by single impulse-like progression patterns. ImpulseDE is an R package suited to capture these patterns in high throughput time series datasets. By fitting a representative impulse model to each gene, it reports differentially expressed genes across time points from a single or between two time courses from two experiments. To optimize running time, the code uses clustering and multi-threading. By applying ImpulseDE , we demonstrate its power to represent underlying biology of gene expression in microarray and RNA-Seq data.Availability and implementation
ImpulseDE is available on Bioconductor ( https://bioconductor.org/packages/ImpulseDE/ ).Contact
niryosef@berkeley.edu.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Sander J
PROVIDER: S-EPMC5859984 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Among multiple subtypes of tissue or cell, subtype-specific differentially-expressed genes (SDEGs) are defined as being most-upregulated in only one subtype but not in any other. Detecting SDEGs plays a critical role in the molecular characterization and deconvolution of multicellular complex tissues. Classic differential analysis assumes a null hypothesis whose test statistic is not subtype-specific, thus can produce a high false positive rate and/or lower detection power. Here we first introduce a One-Versus-Everyone Fold Change (OVE-FC) test for detecting SDEGs. We then propose a scaled test statistic (OVE-sFC) for assessing the statistical significance of SDEGs that applies a mixture null distribution model and a tailored permutation test. The OVE-FC/sFC test was validated on both type 1 error rate and detection power using extensive simulation data sets generated from real gene expression profiles of purified subtype samples. The OVE-FC/sFC test was then applied to two benchmark gene expression data sets of purified subtype samples and detected many known or previously unknown SDEGs. Subsequent supervised deconvolution results on synthesized bulk expression data, obtained using the SDEGs detected from the independent purified expression data by the OVE-FC/sFC test, showed superior performance in deconvolution accuracy when compared with popular peer methods.
| S-EPMC7801594 | biostudies-literature
Project description:We present a general algorithm to detect genes differentially expressed between two nonhomogeneous time-series data sets. As increasing amounts of high-throughput biological data become available, a major challenge in genomic and computational biology is to develop methods for comparing data from different experimental sources. Time-series whole-genome expression data are a particularly valuable source of information because they can describe an unfolding biological process such as the cell cycle or immune response. However, comparisons of time-series expression data sets are hindered by biological and experimental inconsistencies such as differences in sampling rate, variations in the timing of biological processes, and the lack of repeats. Our algorithm overcomes these difficulties by using a continuous representation for time-series data and combining a noise model for individual samples with a global difference measure. We introduce a corresponding statistical method for computing the significance of this differential expression measure. We used our algorithm to compare cell-cycle-dependent gene expression in wild-type and knockout yeast strains. Our algorithm identified a set of 56 differentially expressed genes, and these results were validated by using independent protein-DNA-binding data. Unlike previous methods, our algorithm was also able to identify 22 non-cell-cycle-regulated genes as differentially expressed. This set of genes is significantly correlated in a set of independent expression experiments, suggesting additional roles for the transcription factors Fkh1 and Fkh2 in controlling cellular activity in yeast.
| S-EPMC193530 | biostudies-literature
Project description:BackgroundRNA-seq has become a standard technology to quantify mRNA. The measured values usually vary by several orders of magnitude, and while the detection of differences at high values is statistically well grounded, the significance of the differences for rare mRNAs can be weakened by the presence of biological and technical noise.ResultsWe have developed a method for cleaning RNA-seq data, which improves the detection of differentially expressed genes and specifically genes with low to moderate transcription. Using a data modeling approach, parameters of randomly distributed mRNA counts are identified and reads, most probably originating from technical noise, are removed. We demonstrate that the removal of this random component leads to the significant increase in the number of detected differentially expressed genes, more significant pvalues and no bias towards low-count genes.ConclusionApplication of RNAdeNoise to our RNA-seq data on polysome profiling and several published RNA-seq datasets reveals its suitability for different organisms and sequencing technologies such as Illumina and BGI, shows improved detection of differentially expressed genes, and excludes the subjective setting of thresholds for minimal RNA counts. The program, RNA-seq data, resulted gene lists and examples of use are in the supplementary data and at https://github.com/Deyneko/RNAdeNoise .
| S-EPMC9670425 | biostudies-literature
Project description:MotivationModern high-throughput biotechnologies such as microarray are capable of producing a massive amount of information for each sample. However, in a typical high-throughput experiment, only limited number of samples were assayed, thus the classical 'large p, small n' problem. On the other hand, rapid propagation of these high-throughput technologies has resulted in a substantial collection of data, often carried out on the same platform and using the same protocol. It is highly desirable to utilize the existing data when performing analysis and inference on a new dataset.ResultsUtilizing existing data can be carried out in a straightforward fashion under the Bayesian framework in which the repository of historical data can be exploited to build informative priors and used in new data analysis. In this work, using microarray data, we investigate the feasibility and effectiveness of deriving informative priors from historical data and using them in the problem of detecting differentially expressed genes. Through simulation and real data analysis, we show that the proposed strategy significantly outperforms existing methods including the popular and state-of-the-art Bayesian hierarchical model-based approaches. Our work illustrates the feasibility and benefits of exploiting the increasingly available genomics big data in statistical inference and presents a promising practical strategy for dealing with the 'large p, small n' problem.Availability and implementationOur method is implemented in R package IPBT, which is freely available from https://github.com/benliemory/IPBT CONTACT: yuzhu@purdue.edu; zhaohui.qin@emory.eduSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC4907396 | biostudies-literature
Project description:BackgroundTechnological advances involving RNA-Seq and Bioinformatics allow quantifying the transcriptional levels of genes in cells, tissues, and cell lines, permitting the identification of Differentially Expressed Genes (DEGs). DESeq2 and edgeR are well-established computational tools used for this purpose and they are based upon generalized linear models (GLMs) that consider only fixed effects in modeling. However, the inclusion of random effects reduces the risk of missing potential DEGs that may be essential in the context of the biological phenomenon under investigation. The generalized linear mixed models (GLMM) can be used to include both effects.MethodsWe present DEGRE (Differentially Expressed Genes with Random Effects), a user-friendly tool capable of inferring DEGs where fixed and random effects on individuals are considered in the experimental design of RNA-Seq research. DEGRE preprocesses the raw matrices before fitting GLMMs on the genes and the derived regression coefficients are analyzed using the Wald statistical test. DEGRE offers the Benjamini-Hochberg or Bonferroni techniques for P-value adjustment.ResultsThe datasets used for DEGRE assessment were simulated with known identification of DEGs. These have fixed effects, and the random effects were estimated and inserted to measure the impact of experimental designs with high biological variability. For DEGs' inference, preprocessing effectively prepares the data and retains overdispersed genes. The biological coefficient of variation is inferred from the counting matrices to assess variability before and after the preprocessing. The DEGRE is computationally validated through its performance by the simulation of counting matrices, which have biological variability related to fixed and random effects. DEGRE also provides improved assessment measures for detecting DEGs in cases with higher biological variability. We show that the preprocessing established here effectively removes technical variation from those matrices. This tool also detects new potential candidate DEGs in the transcriptome data of patients with bipolar disorder, presenting a promising tool to detect more relevant genes.ConclusionsDEGRE provides data preprocessing and applies GLMMs for DEGs' inference. The preprocessing allows efficient remotion of genes that could impact the inference. Also, the computational and biological validation of DEGRE has shown to be promising in identifying possible DEGs in experiments derived from complex experimental designs. This tool may help handle random effects on individuals in the inference of DEGs and presents a potential for discovering new interesting DEGs for further biological investigation.
| S-EPMC10078460 | biostudies-literature
Project description:BackgroundTime-course microarray experiments are being increasingly used to characterize dynamic biological processes. In these experiments, the goal is to identify genes differentially expressed in time-course data, measured between different biological conditions. These differentially expressed genes can reveal the changes in biological process due to the change in condition which is essential to understand differences in dynamics.ResultsIn this paper, we propose a novel method for finding differentially expressed genes in time-course data and across biological conditions (say C1 and C2). We model the expression at C1 using Principal Component Analysis and represent the expression profile of each gene as a linear combination of the dominant Principal Components (PCs). Then the expression data from C2 is projected on the developed PCA model and scores are extracted. The difference between the scores is evaluated using a hypothesis test to quantify the significance of differential expression. We evaluate the proposed method to understand differences in two case studies (1) the heat shock response of wild-type and HSF1 knockout mice, and (2) cell-cycle between wild-type and Fkh1/Fkh2 knockout Yeast strains.ConclusionIn both cases, the proposed method identified biologically significant genes.
| S-EPMC2435549 | biostudies-literature
Project description:BackgroundBoolean models of biological signalling-regulatory networks are increasingly used to formally describe and understand complex biological processes. These models may become inconsistent as new data become available and need to be repaired. In the past, the focus has been shed on the inference of (classes of) models given an interaction network and time-series data sets. However, repair of existing models against new data is still in its infancy, where the process is still manually performed and therefore slow and prone to errors.ResultsIn this work, we propose a method with an associated tool to suggest repairs over inconsistent Boolean models, based on a set of atomic repair operations. Answer Set Programming is used to encode the minimal repair problem as a combinatorial optimization problem. In particular, given an inconsistent model, the tool provides the minimal repairs that render the model capable of generating dynamics coherent with a (set of) time-series data set(s), considering either a synchronous or an asynchronous updating scheme.ConclusionsThe method was validated using known biological models from different species, as well as synthetic models obtained from randomly generated networks. We discuss the method's limitations regarding each of the updating schemes and the considered minimization algorithm.
| S-EPMC6434824 | biostudies-literature
Project description:BackgroundOligonucleotide microarrays measure the relative transcript abundance of thousands of mRNAs in parallel. A large number of procedures for normalization and detection of differentially expressed genes have been proposed. However, the relative impact of these methods on the detection of differentially expressed genes remains to be determined.ResultsWe have employed four different normalization methods and all possible combinations with three different statistical algorithms for detection of differentially expressed genes on a prototype dataset. The number of genes detected as differentially expressed differs by a factor of about three. Analysis of lists of genes detected as differentially expressed, and rank correlation coefficients for probability of differential expression shows that a high concordance between different methods can only be achieved by using the same normalization procedure.ConclusionsNormalization has a profound influence of detection of differentially expressed genes. This influence is higher than that of three subsequent statistical analysis procedures examined. Algorithms incorporating more array-derived information than gene-expression values alone are urgently needed.
| S-EPMC126238 | biostudies-literature
Project description:BackgroundOsteoarthritis (OA) is the most common chronic disorder of joints; however, the key genes and transcription factors (TFs) associated with OA are still unclear. Through bioinformatics tools, the study aimed to understand the mechanism of genes associated with the development of OA.MethodsFour gene expression profiling datasets were used to identify differentially expressed genes (DEGs) between OA and healthy control samples by a meta-analysis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed with Multifaceted Analysis Tool for Human Transcriptome (MATHT). Subsequently, a protein-protein interaction (PPI) network was constructed for these DEGs. Significant network modules were identified using ReactomeFIViz, and the pathway of each module was enriched using MATHT. In addition, TFs in the DEGs were identified.ResultsIn total, 690 DEGs were identified between OA and healthy control samples, including 449 upregulated and 241 downregulated DEGs. Additionally, 622 nodes and 2752 interactions constituted the PPI network, including 401 upregulated and 221 downregulated DEGs. Among them, FOS, TWIST1, POU2F1, SMARCA4, and CREBBP were also identified as TFs. RT-PCR results showed that the expression levels of Fos, Twist1, Pou2f1, Smarca4, and Crebbp decreased in mice with OA. In addition, FOS, TWIST1, SMARCA4, and CREBBP were involved in the positive regulation of transcription from the RNA polymerase II promoter.ConclusionsTWIST1, POU2F1, SMARCA4, and CREBBP may play an important role in OA pathology.
| S-EPMC5842536 | biostudies-literature
Project description:Colorectal cancer (CRC) is a prevalent cancer with high morbidity and mortality rates worldwide. Late diagnosis is a significant contributor to low survival rates in a minority of cases. The study aimed to perform a robust pipeline using integrated bioinformatics tools that will enable us to identify potential diagnostic and prognostic biomarkers for early detection of CRC by exploring differentially expressed genes (DEGs). In addition to, testing the capability of replacing chemotherapy with plant extract in CRC treatment by validating it using real-time PCR. RNA-seq data from cancerous and adjacent normal tissues were pre-processed and analyzed using various tools such as FastQC, Kallisto, DESeq@ R package, g:Profiler, GNEMANIA-CytoScape and CytoHubba, resulting in the identification of 1641 DEGs enriched in various signaling routes. MMP7, TCF21, and VEGFD were found to be promising diagnostic biomarkers for CRC. An in vitro experiment was conducted to examine the potential anticancer properties of 5-fluorouracile, Withania somnifera extract, and their combination. The extract was found to exhibit a positive trend in gene expression and potential therapeutic value by targeting the three genes; however, further trials are required to regulate the methylation promoter. Molecular docking tests supported the findings by revealing a stable ligand-receptor complex. In conclusion, the study's analysis workflow is precise and robust in identifying DEGs in CRC that may serve as biomarkers for diagnosis and treatment. Additionally, the identified DEGs can be used in future research with larger sample sizes to analyze CRC survival.
| S-EPMC11289010 | biostudies-literature