MASCARA: coexpression analysis in data from designed experiments

Abstract

Experiments in plant transcriptomics are usually designed to induce variation in a pathway of interest. Harsh experimental conditions can cause widespread transcriptional changes between groups. Discovering coexpression within a pathway of interest (here the strigolactone pathway) in this context is hampered by the dominant variance induced by the design. Minor changes in experimental conditions not controlled for may affect the plants, leading to small coordinated differences in genes within pathways of interest and related pathways between replicate plants in the same controlled experimental condition. These systematic differences are usually averaged out, but we argue here that they can be used to improve the detection of genes that co-express. We introduce a novel framework "MASCARA" which combines ANOVA simultaneous component analysis and partial least squares to remove the experimentally induced variance and investigate multivariate relationships in the non-designed variance. MASCARA is tested against a selection of competitors on simulated data, created to mimic a designed transcriptome study, where its benefit is demonstrated. In a coexpression analysis of a real dataset MASCARA detects several uncharacterised but relevant transcripts. Our results indicate that there is sufficient structure left in a typical dataset after correcting for experimental variance and that this residual information is useful to investigate coexpression. Author SummaryExperiments in the life sciences usually purposefully induce significant variance between different treatments, in order to activate or repress certain mechanisms of interest. Whilst this is necessary it can make it challenging to detect meaningful relationships within pathways of interest, particularly when the experimental conditions are drastically different. Instead of focusing on the drastic changes in response due to the different treatment, MASCARA uses the systematic synchronous variances between replicates to find related features within the pathway of interest. Through simulation studies and application to a real dataset, we demonstrate the effectiveness of MASCARA in detecting relevant transcripts and extracting coexpression patterns from gene expression data.