Addressing the Mean-correlation Relationship in Co-expression Analysis

Overview

Journal PLoS Comput Biol

Specialty Biology

Date 2022 Mar 30

PMID 35353807

Authors

Yi Wang

Stephanie C Hicks

Kasper D Hansen

Affiliations

Soon will be listed here.

Abstract

Estimates of correlation between pairs of genes in co-expression analysis are commonly used to construct networks among genes using gene expression data. As previously noted, the distribution of such correlations depends on the observed expression level of the involved genes, which we refer to this as a mean-correlation relationship in RNA-seq data, both bulk and single-cell. This dependence introduces an unwanted technical bias in co-expression analysis whereby highly expressed genes are more likely to be highly correlated. Such a relationship is not observed in protein-protein interaction data, suggesting that it is not reflecting biology. Ignoring this bias can lead to missing potentially biologically relevant pairs of genes that are lowly expressed, such as transcription factors. To address this problem, we introduce spatial quantile normalization (SpQN), a method for normalizing local distributions in a correlation matrix. We show that spatial quantile normalization removes the mean-correlation relationship and corrects the expression bias in network reconstruction.

Citing Articles

Cell-type-specific mapping of enhancers and target genes from single-cell multimodal data.

Su C, Lee D, Jin P, Zhang J bioRxiv. 2024; .

PMID: 39386519 PMC: 11463474. DOI: 10.1101/2024.09.24.614814.

Massively integrated coexpression analysis reveals transcriptional regulation, evolution and cellular implications of the yeast noncanonical translatome.

Rich A, Acar O, Carvunis A Genome Biol. 2024; 25(1):183.

PMID: 38978079 PMC: 11232214. DOI: 10.1186/s13059-024-03287-7.

eQTLs identify regulatory networks and drivers of variation in the individual response to sepsis.

Burnham K, Milind N, Lee W, Kwok A, Cano-Gamez K, Mi Y Cell Genom. 2024; 4(7):100587.

PMID: 38897207 PMC: 11293594. DOI: 10.1016/j.xgen.2024.100587.

Network-based drug repurposing for schizophrenia.

Truong T, Liu Z, Panizzutti B, Kim J, Dean O, Berk M Neuropsychopharmacology. 2024; 49(6):983-992.

PMID: 38321095 PMC: 11039639. DOI: 10.1038/s41386-024-01805-6.

Cell-type-specific co-expression inference from single cell RNA-sequencing data.

Su C, Xu Z, Shan X, Cai B, Zhao H, Zhang J Nat Commun. 2023; 14(1):4846.

PMID: 37563115 PMC: 10415381. DOI: 10.1038/s41467-023-40503-7.

References

Leek J, Johnson W, Parker H, Jaffe A, Storey J . The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012; 28(6):882-3. PMC: 3307112. DOI: 10.1093/bioinformatics/bts034. View

Oldham M, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath S . Functional organization of the transcriptome in human brain. Nat Neurosci. 2008; 11(11):1271-82. PMC: 2756411. DOI: 10.1038/nn.2207. View

Battle A, Brown C, Engelhardt B, Montgomery S . Genetic effects on gene expression across human tissues. Nature. 2017; 550(7675):204-213. PMC: 5776756. DOI: 10.1038/nature24277. View

Marioni J, Mason C, Mane S, Stephens M, Gilad Y . RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008; 18(9):1509-17. PMC: 2527709. DOI: 10.1101/gr.079558.108. View

Wang Q, Armenia J, Zhang C, Penson A, Reznik E, Zhang L . Unifying cancer and normal RNA sequencing data from different sources. Sci Data. 2018; 5:180061. PMC: 5903355. DOI: 10.1038/sdata.2018.61. View

Leek J, Storey J . Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007; 3(9):1724-35. PMC: 1994707. DOI: 10.1371/journal.pgen.0030161. View

Ruprecht C, Proost S, Hernandez-Coronado M, Ortiz-Ramirez C, Lang D, Rensing S . Phylogenomic analysis of gene co-expression networks reveals the evolution of functional modules. Plant J. 2017; 90(3):447-465. DOI: 10.1111/tpj.13502. View

Leek J, Storey J . A general framework for multiple testing dependence. Proc Natl Acad Sci U S A. 2008; 105(48):18718-23. PMC: 2586646. DOI: 10.1073/pnas.0808709105. View

Boukas L, Havrilla J, Hickey P, Quinlan A, Bjornsson H, Hansen K . Coexpression patterns define epigenetic regulators associated with neurological dysfunction. Genome Res. 2019; 29(4):532-542. PMC: 6442390. DOI: 10.1101/gr.239442.118. View

10.

Johnson W, Li C, Rabinovic A . Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2006; 8(1):118-27. DOI: 10.1093/biostatistics/kxj037. View

11.

Parsana P, Ruberman C, Jaffe A, Schatz M, Battle A, Leek J . Addressing confounding artifacts in reconstruction of gene co-expression networks. Genome Biol. 2019; 20(1):94. PMC: 6521369. DOI: 10.1186/s13059-019-1700-9. View

12.

Ghazalpour A, Doss S, Zhang B, Wang S, Plaisier C, Castellanos R . Integrating genetic and network analysis to characterize genes related to mouse weight. PLoS Genet. 2006; 2(8):e130. PMC: 1550283. DOI: 10.1371/journal.pgen.0020130. View

13.

Bullard J, Purdom E, Hansen K, Dudoit S . Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics. 2010; 11:94. PMC: 2838869. DOI: 10.1186/1471-2105-11-94. View

14.

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

15.

Barrera L, Vedenko A, Kurland J, Rogers J, Gisselbrecht S, Rossin E . Survey of variation in human transcription factors reveals prevalent DNA binding changes. Science. 2016; 351(6280):1450-1454. PMC: 4825693. DOI: 10.1126/science.aad2257. View

16.

Friedman J, Hastie T, Tibshirani R . Sparse inverse covariance estimation with the graphical lasso. Biostatistics. 2007; 9(3):432-41. PMC: 3019769. DOI: 10.1093/biostatistics/kxm045. View

17.

Willsey A, Sanders S, Li M, Dong S, Tebbenkamp A, Muhle R . Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism. Cell. 2013; 155(5):997-1007. PMC: 3995413. DOI: 10.1016/j.cell.2013.10.020. View

18.

Workman C, Jensen L, Jarmer H, Berka R, Gautier L, Nielser H . A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biol. 2002; 3(9):research0048. PMC: 126873. DOI: 10.1186/gb-2002-3-9-research0048. View

19.

Freytag S, Gagnon-Bartsch J, Speed T, Bahlo M . Systematic noise degrades gene co-expression signals but can be corrected. BMC Bioinformatics. 2015; 16:309. PMC: 4583191. DOI: 10.1186/s12859-015-0745-3. View

20.

van Dam S, Vosa U, van der Graaf A, Franke L, de Magalhaes J . Gene co-expression analysis for functional classification and gene-disease predictions. Brief Bioinform. 2017; 19(4):575-592. PMC: 6054162. DOI: 10.1093/bib/bbw139. View