Fold-change threshold screening: a robust algorithm to unmask hidden gene expression patterns in noisy aggregated transcriptome data

Abstract

Transcriptomics is often used to investigate changes in an organism's genetic response to environmental contamination. Data noise can mask the effects of contaminants making it difficult to detect responding genes. Because the number of genes which are found differentially expressed in transcriptome data is often very large, algorithms are needed to reduce the number down to a few robust discriminative genes. We present an algorithm for aggregated analysis of transcriptome data which uses multiple fold-change thresholds (threshold screening) and p values from Bayesian generalized linear model in order to assess the robustness of a gene as a potential indicator for the treatments tested. The algorithm provides a robustness indicator (ROBI) as well as a significance profile, which can be used to assess the statistical significance of a given gene for different fold-change thresholds. Using ROBI, eight discriminative genes were identified from an exemplary dataset (Danio rerio FET treated with chlorpyrifos, methylmercury, and PCB) which could be potential indicators for a given substance. Significance profiles uncovered genetic effects and revealed appropriate fold-change thresholds for single genes or gene clusters. Fold-change threshold screening is a powerful tool for dimensionality reduction and feature selection in transcriptome data, as it effectively reduces the number of detected genes suitable for environmental monitoring. In addition, it is able to unmask patterns in altered genetic expression hidden by data noise and reduces the chance of type II errors, e.g., in environmental screening.

Multi-omics analysis of multiple missions to space reveal a theme of lipid dysregulation in mouse liver

“…Pathway analysis and subsequent predictions in each tissue were done using the statistically significant genes with a fold-change ≥ 1.2 (or ≤ −1.2) comparing flight conditions versus habitat ground controls (with either p-value < 0.05 or FDR < 0.05 as stated above). Low fold-change has become quite standard when trying to identify genes that are differentially expressed 67–69 , as low cutoffs are less affected by different data normalization schemes and they are less likely to eliminate key genes operating under very tight level regulation. Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems) was used to predict statistically significant activation or inhibition of upstream regulators using activation Z -score statistics ( ≥ 2, activated or ≤ −2, inhibited) 70 .…”

Multiple Pathways Involved in Palmitic Acid-Induced Toxicity: A System Biology Approach

“…Here, the function takes as argument two valid models for the "sybil" R package and a customizable threshold value to filter functions to be reported. In this aspect, we chose an arbitrary threshold value greater or equal to 2-fold times for reactions with an absolute change between the unconstrained and constrained metabolic scenarios, as reported in previous models (Hausen et al, 2015;Banos et al, 2017).…”

Danio rerio as a model in aquatic toxicology and sediment research

“…New and valuable knowledge concerning the application of such bioassays and the underlying cellular mechanisms have been gained. Moreover, microarray analyses were proven to be able to identify sets of genes, which represent various important cellular pathways and react to environmental contaminants (Hausen et al 2015).…”