Dataset Information

Correcting for experiment-specific variability in expression compendia can remove underlying signals.

ABSTRACT: MOTIVATION:In the past two decades, scientists in different laboratories have assayed gene expression from millions of samples. These experiments can be combined into compendia and analyzed collectively to extract novel biological patterns. Technical variability, or "batch effects," may result from combining samples collected and processed at different times and in different settings. Such variability may distort our ability to extract true underlying biological patterns. As more integrative analysis methods arise and data collections get bigger, we must determine how technical variability affects our ability to detect desired patterns when many experiments are combined. OBJECTIVE:We sought to determine the extent to which an underlying signal was masked by technical variability by simulating compendia comprising data aggregated across multiple experiments. METHOD:We developed a generative multi-layer neural network to simulate compendia of gene expression experiments from large-scale microbial and human datasets. We compared simulated compendia before and after introducing varying numbers of sources of undesired variability. RESULTS:The signal from a baseline compendium was obscured when the number of added sources of variability was small. Applying statistical correction methods rescued the underlying signal in these cases. However, as the number of sources of variability increased, it became easier to detect the original signal even without correction. In fact, statistical correction reduced our power to detect the underlying signal. CONCLUSION:When combining a modest number of experiments, it is best to correct for experiment-specific noise. However, when many experiments are combined, statistical correction reduces our ability to extract underlying patterns.

SUBMITTER: Lee AJ

PROVIDER: S-EPMC7607552 | biostudies-literature | 2020 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Correcting for experiment-specific variability in expression compendia can remove underlying signals.

Lee Alexandra J AJ Park YoSon Y Doing Georgia G Hogan Deborah A DA Greene Casey S CS

GigaScience 20201101 11

<h4>Motivation</h4>In the past two decades, scientists in different laboratories have assayed gene expression from millions of samples. These experiments can be combined into compendia and analyzed collectively to extract novel biological patterns. Technical variability, or "batch effects," may result from combining samples collected and processed at different times and in different settings. Such variability may distort our ability to extract true underlying biological patterns. As more integra ...[more]

PMID: 33140829

Similar Datasets

Project description:BACKGROUND: Compendia of gene expression profiles under chemical and genetic perturbations constitute an invaluable resource from a systems biology perspective. However, the perturbational nature of such data imposes specific challenges on the computational methods used to analyze them. In particular, traditional clustering algorithms have difficulties in handling one of the prominent features of perturbational compendia, namely partial coexpression relationships between genes. Biclustering methods on the other hand are specifically designed to capture such partial coexpression patterns, but they show a variety of other drawbacks. For instance, some biclustering methods are less suited to identify overlapping biclusters, while others generate highly redundant biclusters. Also, none of the existing biclustering tools takes advantage of the staple of perturbational expression data analysis: the identification of differentially expressed genes. RESULTS: We introduce a novel method, called ENIGMA, that addresses some of these issues. ENIGMA leverages differential expression analysis results to extract expression modules from perturbational gene expression data. The core parameters of the ENIGMA clustering procedure are automatically optimized to reduce the redundancy between modules. In contrast to the biclusters produced by most other methods, ENIGMA modules may show internal substructure, i.e. subsets of genes with distinct but significantly related expression patterns. The grouping of these (often functionally) related patterns in one module greatly aids in the biological interpretation of the data. We show that ENIGMA outperforms other methods on artificial datasets, using a quality criterion that, unlike other criteria, can be used for algorithms that generate overlapping clusters and that can be modified to take redundancy between clusters into account. Finally, we apply ENIGMA to the Rosetta compendium of expression profiles for Saccharomyces cerevisiae and we analyze one pheromone response-related module in more detail, demonstrating the potential of ENIGMA to generate detailed predictions. CONCLUSION: It is increasingly recognized that perturbational expression compendia are essential to identify the gene networks underlying cellular function, and efforts to build these for different organisms are currently underway. We show that ENIGMA constitutes a valuable addition to the repertoire of methods to analyze such data.

Project description:Blood Oxygen Level Dependent (BOLD) functional MRI is a complex neurovascular signal whose magnitude depends on baseline physiological factors such as cerebral blood flow (CBF). Because baseline CBF varies across the brain and is altered with aging, the interpretation of stand-alone aging-related BOLD changes can be misleading. The primary objective of this study was to develop a methodology that combines task fMRI and arterial spin labeling (ASL) techniques to sensitize task-induced BOLD activity by covarying out the baseline physiology (i.e., CBF) in an aging model. We recruited 11 younger and 13 older healthy participants who underwent ASL and an overt language fMRI task (semantic category member generation). We measured in-scanner language performance to investigate the effect of BOLD sensitization on BOLD-behavior relationships. The results demonstrate that our correction approach is effective at enhancing the specificity and sensitivity of the BOLD signal in both groups. In addition, the correction strengthens the statistical association between task BOLD activity and behavioral performance. Although CBF has inherent age dependence, our results show that retaining the age factor within CBF aides in greater sensitization of task fMRI signals. From a cognitive standpoint, compared to young adults, the older participants showed a delayed domain-general language-related task activity possibly due to compromised vessel compliance. Further, assessment of functional evolution of corrected BOLD activity revealed biphasic BOLD dynamics in both groups where BOLD deactivation may reflect greater semantic demand or increased premium on domain general executive functioning in response to task difficulty. Although it was promising to note that the predictability of behavior using the proposed methodology outperforms other methodologies (i.e., no correction and normalization by division), and provides moderate stability and adequate power, further work with a larger cohort and other task designs is necessary to improve the stability of predicting associated behavior. In summary, we recommend correction of task fMRI signals by covarying out baseline CBF especially when comparing groups with different neurovascular properties. Given that ASL and BOLD fMRI are well established and widely employed techniques, our proposed multi-modal methodology can be readily implemented into data processing pipelines to obtain more accurate BOLD activation maps.

Project description:BackgroundMicroarrays are the main technology for large-scale transcriptional gene expression profiling, but the large bodies of data available in public databases are not useful due to the large heterogeneity. There are several initiatives that attempt to bundle these data into expression compendia, but such resources for bacterial organisms are scarce and limited to integration of experiments from the same platform or to indirect integration of per experiment analysis results.Methodology/principal findingsWe have constructed comprehensive organism-specific cross-platform expression compendia for three bacterial model organisms (Escherichia coli, Bacillus subtilis, and Salmonella enterica serovar Typhimurium) together with an access portal, dubbed COLOMBOS, that not only provides easy access to the compendia, but also includes a suite of tools for exploring, analyzing, and visualizing the data within these compendia. It is freely available at http://bioi.biw.kuleuven.be/colombos. The compendia are unique in directly combining expression information from different microarray platforms and experiments, and we illustrate the potential benefits of this direct integration with a case study: extending the known regulon of the Fur transcription factor of E. coli. The compendia also incorporate extensive annotations for both genes and experimental conditions; these heterogeneous data are functionally integrated in the COLOMBOS analysis tools to interactively browse and query the compendia not only for specific genes or experiments, but also metabolic pathways, transcriptional regulation mechanisms, experimental conditions, biological processes, etc.Conclusions/significanceWe have created cross-platform expression compendia for several bacterial organisms and developed a complementary access port COLOMBOS, that also serves as a convenient expression analysis tool to extract useful biological information. This work is relevant to a large community of microbiologists by facilitating the use of publicly available microarray experiments to support their research.

Project description:Thousands of Pseudomonas aeruginosa RNA sequencing (RNA-seq) gene expression profiles are publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). In this work, the transcriptional profiles from hundreds of studies performed by over 75 research groups were reanalyzed in aggregate to create a powerful tool for hypothesis generation and testing. Raw sequence data were uniformly processed using the Salmon pseudoaligner, and this read mapping method was validated by comparison to a direct alignment method. We developed filtering criteria to exclude samples with aberrant levels of housekeeping gene expression or an unexpected number of genes with no reported values and normalized the filtered compendia using the ratio-of-medians method. The filtering and normalization steps greatly improved gene expression correlations for genes within the same operon or regulon across the 2,333 samples. Since the RNA-seq data were generated using diverse strains, we report the effects of mapping samples to noncognate reference genomes by separately analyzing all samples mapped to cDNA reference genomes for strains PAO1 and PA14, two divergent strains that were used to generate most of the samples. Finally, we developed an algorithm to incorporate new data as they are deposited into the SRA. Our processing and quality control methods provide a scalable framework for taking advantage of the troves of biological information hibernating in the depths of microbial gene expression data and yield useful tools for P. aeruginosa RNA-seq data to be leveraged for diverse research goals. IMPORTANCE Pseudomonas aeruginosa is a causative agent of a wide range of infections, including chronic infections associated with cystic fibrosis. These P. aeruginosa infections are difficult to treat and often have negative outcomes. To aid in the study of this problematic pathogen, we mapped, filtered for quality, and normalized thousands of P. aeruginosa RNA-seq gene expression profiles that were publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The resulting compendia facilitate analyses across experiments, strains, and conditions. Ultimately, the workflow that we present could be applied to analyses of other microbial species.

Dataset Information

Correcting for experiment-specific variability in expression compendia can remove underlying signals.

Publications

Correcting for experiment-specific variability in expression compendia can remove underlying signals.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets