{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,2,21]],"date-time":"2025-02-21T10:52:33Z","timestamp":1740135153782,"version":"3.37.3"},"reference-count":29,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2021,9,7]],"date-time":"2021-09-07T00:00:00Z","timestamp":1630972800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2021,9,7]],"date-time":"2021-09-07T00:00:00Z","timestamp":1630972800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"name":"Data Science and System Complexity Centre","award":["University of Groningen"],"award-info":[{"award-number":["University of Groningen"]}]},{"name":"Netherlands X-omics Initiative","award":["184.034.019"],"award-info":[{"award-number":["184.034.019"]}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"published-print":{"date-parts":[[2021,12]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec>\n                <jats:title>Background<\/jats:title>\n                <jats:p>In systems biology, it is important to reconstruct regulatory networks from quantitative molecular profiles. Gaussian graphical models (GGMs) are one of the most popular methods to this end. A GGM consists of nodes (representing the transcripts, metabolites or proteins) inter-connected by edges (reflecting their partial correlations). Learning the edges from quantitative molecular profiles is statistically challenging, as there are usually fewer samples than nodes (\u2018high dimensional problem\u2019). Shrinkage methods address this issue by learning a regularized GGM. However, it remains open to study how the shrinkage affects the final result and its interpretation.<\/jats:p>\n              <\/jats:sec><jats:sec>\n                <jats:title>Results<\/jats:title>\n                <jats:p>We show that the shrinkage biases the partial correlation in a non-linear way. This bias does not only change the magnitudes of the partial correlations but also affects their order. Furthermore, it makes networks obtained from different experiments incomparable and hinders their biological interpretation. We propose a method, referred to as \u2018un-shrinking\u2019 the partial correlation, which corrects for this non-linear bias. Unlike traditional methods, which use a fixed shrinkage value, the new approach provides partial correlations that are closer to the actual (population) values and that are easier to interpret. This is demonstrated on two gene expression datasets from <jats:italic>Escherichia coli<\/jats:italic> and <jats:italic>Mus musculus.<\/jats:italic><\/jats:p>\n              <\/jats:sec><jats:sec>\n                <jats:title>Conclusions<\/jats:title>\n                <jats:p>GGMs are popular undirected graphical models based on partial correlations. The application of GGMs to reconstruct regulatory networks is commonly performed using shrinkage to overcome the \u2018high-dimensional problem\u2019. Besides it advantages, we have identified that the shrinkage introduces a non-linear bias in the partial correlations. Ignoring this type of effects caused by the shrinkage can obscure the interpretation of the network, and impede the validation of earlier reported results.<\/jats:p>\n              <\/jats:sec>","DOI":"10.1186\/s12859-021-04313-2","type":"journal-article","created":{"date-parts":[[2021,9,7]],"date-time":"2021-09-07T19:02:21Z","timestamp":1631041341000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["The \u2018un-shrunk\u2019 partial correlation in Gaussian graphical models"],"prefix":"10.1186","volume":"22","author":[{"given":"Victor","family":"Bernal","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Rainer","family":"Bischoff","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2218-1140","authenticated-orcid":false,"given":"Peter","family":"Horvatovich","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Victor","family":"Guryev","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Marco","family":"Grzegorczyk","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2021,9,7]]},"reference":[{"key":"4313_CR1","doi-asserted-by":"publisher","first-page":"428","DOI":"10.1007\/0-387-21679-0_19","volume-title":"The analysis of gene expression data","author":"AJ Butte","year":"2003","unstructured":"Butte AJ, Kohane IS. Relevance networks: a first step toward finding genetic regulatory networks within microarray data. In: Parmigiani G, Garrett ES, Irizarry RAZSL, editors. The analysis of gene expression data. New York: Springer; 2003. p. 428\u201346."},{"key":"4313_CR2","doi-asserted-by":"publisher","first-page":"601","DOI":"10.1089\/106652700750050961","volume":"7","author":"N Friedman","year":"2000","unstructured":"Friedman N, Linial M, Nachman I, Pe\u2019er D, . Using Bayesian networks to analyze expression data. J Comput Biol. 2000;7:601\u201320.","journal-title":"J Comput Biol"},{"key":"4313_CR3","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4612-0493-0","volume-title":"Introduction to graphical modelling, 2nd edition","author":"D Edwards","year":"2000","unstructured":"Edwards D. Introduction to graphical modelling, 2nd edition. 2nd ed. New York: Springer-Verlag; 2000.","edition":"2"},{"key":"4313_CR4","first-page":"1","volume":"2019","author":"VA Huynh-Thu","year":"2019","unstructured":"Huynh-Thu VA, Sanguinetti G. Gene regulatory network inference: an introductory survey. Gene Regul Netw. 2019;2019:1\u201323.","journal-title":"Gene Regul Netw"},{"key":"4313_CR5","doi-asserted-by":"publisher","first-page":"2523","DOI":"10.1093\/bioinformatics\/btl391","volume":"22","author":"AV Werhli","year":"2006","unstructured":"Werhli AV, Grzegorczyk M, Husmeier D. Comparative evaluation of reverse engineering gene regulatory networks with relevance networks, graphical gaussian models and bayesian networks. Bioinformatics. 2006;22:2523\u201331.","journal-title":"Bioinformatics"},{"key":"4313_CR6","doi-asserted-by":"publisher","first-page":"S3","DOI":"10.1186\/1471-2105-8-S2-S3","volume":"8","author":"R Opgen-Rhein","year":"2007","unstructured":"Opgen-Rhein R, Strimmer K. Learning causal networks from systems biology time course data: An effective model selection procedure for the vector autoregressive process. BMC Bioinformatics. 2007;8:S3.","journal-title":"BMC Bioinformatics"},{"key":"4313_CR7","doi-asserted-by":"publisher","first-page":"37","DOI":"10.1186\/1752-0509-1-37","volume":"1","author":"R Opgen-Rhein","year":"2007","unstructured":"Opgen-Rhein R, Strimmer K. From correlation to causation networks: A simple approximate learning algorithm and its application to high-dimensional plant gene expression data. BMC Syst Biol. 2007;1:37.","journal-title":"BMC Syst Biol"},{"key":"4313_CR8","doi-asserted-by":"publisher","first-page":"e225","DOI":"10.1371\/journal.pcbi.0030225","volume":"3","author":"N Beerenwinkel","year":"2007","unstructured":"Beerenwinkel N, Antal T, Dingli D, Traulsen A, Kinzler KW, Velculescu VE, et al. Genetic progression and the waiting time to cancer. PLoS Comput Biol. 2007;3:e225.","journal-title":"PLoS Comput Biol"},{"key":"4313_CR9","doi-asserted-by":"publisher","first-page":"1483","DOI":"10.1038\/s41467-017-01525-0","volume":"8","author":"E Benedetti","year":"2017","unstructured":"Benedetti E, Pu\u010di\u0107-Bakovi\u0107 M, Keser T, Wahl A, Hassinen A, Yang JY, et al. Network inference from glycoproteomics data reveals new reactions in the IgG glycosylation pathway. Nat Commun. 2017;8:1483.","journal-title":"Nat Commun"},{"key":"4313_CR10","doi-asserted-by":"publisher","first-page":"gr074914107","DOI":"10.1101\/gr.074914.107","volume":"18","author":"MP Keller","year":"2008","unstructured":"Keller MP, Choi Y, Wang P, Davis DB, Rabaglia ME, Oler AT, et al. A gene expression network model of type 2 diabetes links cell cycle regulation in islets with diabetes susceptibility. Genome Res. 2008;18:gr074914107.","journal-title":"Genome Res"},{"key":"4313_CR11","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1038\/s41598-019-50963-x","volume":"15835","author":"K Imkamp","year":"2019","unstructured":"Imkamp K, Bernal V, Grzegorzcyk M, Horvatovich P, Vermeulen CJ, Heijink IH, et al. Gene network approach reveals co-expression patterns in nasal and bronchial epithelium. Sci Rep. 2019;15835:1\u201313. https:\/\/doi.org\/10.1038\/s41598-019-50963-x.","journal-title":"Sci Rep"},{"key":"4313_CR12","doi-asserted-by":"publisher","first-page":"603","DOI":"10.1162\/NECO_a_00936","volume":"29","author":"A Das","year":"2017","unstructured":"Das A, Sampson AL, Lainscsek C, Muller L, Lin W, Doyle JC, et al. Interpretation of the precision matrix and its application in estimating sparse brain connectivity during sleep spindles from human electrocorticography recordings. Neural Comput. 2017;29:603.","journal-title":"Neural Comput"},{"key":"4313_CR13","doi-asserted-by":"crossref","unstructured":"McNally RJ, Robinaugh DJ, Wu GWY, Wang L, Deserno MK, Borsboom D. Mental disorders as causal systems: a network approach to posttraumatic stress disorder. Clin Psychol Sci. 2015.","DOI":"10.1177\/2167702614553230"},{"key":"4313_CR14","doi-asserted-by":"publisher","first-page":"365","DOI":"10.1016\/S0047-259X(03)00096-4","volume":"88","author":"O Ledoit","year":"2004","unstructured":"Ledoit O, Wolf M. A well-conditioned estimator for large-dimensional covariance matrices. J Multivar Anal. 2004;88:365\u2013411.","journal-title":"J Multivar Anal"},{"key":"4313_CR15","doi-asserted-by":"publisher","first-page":"432","DOI":"10.1093\/biostatistics\/kxm045","volume":"9","author":"J Friedman","year":"2008","unstructured":"Friedman J, Hastie T, Tibshirani R. Sparse inverse covariance estimation with the graphical lasso. Biostatistics. 2008;9:432.","journal-title":"Biostatistics"},{"key":"4313_CR16","doi-asserted-by":"publisher","first-page":"603","DOI":"10.1016\/S0927-5398(03)00007-0","volume":"10","author":"O Ledoit","year":"2003","unstructured":"Ledoit O, Wolf M. Improved estimation of the covariance matrix of stock returns with an application to portfolio selection. J Empir Financ. 2003;10:603\u201321.","journal-title":"J Empir Financ"},{"key":"4313_CR17","doi-asserted-by":"publisher","first-page":"1175","DOI":"10.2202\/1544-6115.1175","volume":"4","author":"J Sch\u00e4fer","year":"2005","unstructured":"Sch\u00e4fer J, Strimmer K. A shrinkage approach to large-scale covariance matrix estimation and implications for functional genomics. Stat Appl Genet Mol Biol. 2005;4:1175\u201389. https:\/\/doi.org\/10.2202\/1544-6115.1175.","journal-title":"Stat Appl Genet Mol Biol"},{"key":"4313_CR18","doi-asserted-by":"publisher","first-page":"5011","DOI":"10.1093\/bioinformatics\/btz357","volume":"35","author":"V Bernal","year":"2019","unstructured":"Bernal V, Bischoff R, Guryev V, Grzegorczyk M, Horvatovich P. Exact hypothesis testing for shrinkage-based Gaussian graphical models. Bioinformatics. 2019;35:5011.","journal-title":"Bioinformatics"},{"key":"4313_CR19","first-page":"1059","volume":"13","author":"T Zhao","year":"2012","unstructured":"Zhao T, Liu H, Roeder K, Lafferty J, Wasserman L. The huge package for high-dimensional undirected graph estimation in R. J Mach Learn Res. 2012;13:1059.","journal-title":"J Mach Learn Res"},{"key":"4313_CR20","unstructured":"Schmidt-Heck W, Guthke R, Toepfer S, Reischer H, Duerrschmid K, Bayer K, et al. Reverse engineering of the stress response during expression of a recombinant protein. Proc EUNITE Symp. 2004; 10\u201312. http:\/\/citeseerx.ist.psu.edu\/viewdoc\/download?doi=10.1.1.129.539&rep=rep1&type=pdf."},{"key":"4313_CR21","doi-asserted-by":"publisher","first-page":"e17820","DOI":"10.1371\/journal.pone.0017820","volume":"6","author":"D Bottomly","year":"2011","unstructured":"Bottomly D, Walter NAR, Hunter JE, Darakjian P, Kawane S, Buck KJ, et al. Evaluating gene expression in C57BL\/6J and DBA\/2J mouse striatum using RNA-Seq and microarrays. PLoS ONE. 2011;6:e17820.","journal-title":"PLoS ONE"},{"key":"4313_CR22","doi-asserted-by":"publisher","first-page":"e47","DOI":"10.1093\/nar\/gkv007","volume":"43","author":"ME Ritchie","year":"2015","unstructured":"Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.","journal-title":"Nucleic Acids Res"},{"key":"4313_CR23","doi-asserted-by":"publisher","first-page":"289","DOI":"10.2307\/2346101","volume":"57","author":"Y Benjamini","year":"1995","unstructured":"Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57:289\u2013300. https:\/\/doi.org\/10.2307\/2346101.","journal-title":"J R Stat Soc"},{"key":"4313_CR24","volume-title":"Statistical power analysis for the behavioural sciences","author":"J Cohen","year":"1988","unstructured":"Cohen J. Statistical power analysis for the behavioural sciences. NJ Lawrence Earlbaum Assoc: Hillside; 1988."},{"key":"4313_CR25","first-page":"gkw937","volume":"2017","author":"D Szklarczyk","year":"2017","unstructured":"Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;2017:gkw937.","journal-title":"Nucleic Acids Res"},{"key":"4313_CR26","first-page":"28","volume":"43","author":"A Franceschini","year":"2013","unstructured":"Franceschini A, Simonovic M, Roth A, Von Mering C, Szklarczyk D, Frankild S, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;43:28.","journal-title":"Nucleic Acids Res"},{"key":"4313_CR27","doi-asserted-by":"publisher","first-page":"122","DOI":"10.1090\/S0002-9947-1941-0003498-3","volume":"49","author":"AC Berry","year":"1941","unstructured":"Berry AC. The accuracy of the Gaussian approximation to the sum of independent variates. Trans Am Math Soc. 1941;49:122.","journal-title":"Trans Am Math Soc"},{"key":"4313_CR28","first-page":"24","volume":"2003","author":"DW Zimmerman","year":"2003","unstructured":"Zimmerman DW, Zumbo BD, Williams RH. Bias in estimation and hypothesis testing of correlation. Psicol\u00f3gica. 2003;2003:24.","journal-title":"Psicol\u00f3gica"},{"key":"4313_CR29","unstructured":"Whittaker J. Graphical models in applied multivariate statistics. Wiley Publishing. 1990."}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-021-04313-2.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1186\/s12859-021-04313-2\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-021-04313-2.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,9,7]],"date-time":"2021-09-07T19:02:37Z","timestamp":1631041357000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-021-04313-2"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,9,7]]},"references-count":29,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2021,12]]}},"alternative-id":["4313"],"URL":"https:\/\/doi.org\/10.1186\/s12859-021-04313-2","relation":{},"ISSN":["1471-2105"],"issn-type":[{"type":"electronic","value":"1471-2105"}],"subject":[],"published":{"date-parts":[[2021,9,7]]},"assertion":[{"value":"11 September 2020","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"2 August 2021","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"7 September 2021","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Declarations"}},{"value":"Not applicable.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Ethics approval and consent to participate"}},{"value":"Not applicable.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Consent for publication"}},{"value":"None.","order":4,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing interests"}}],"article-number":"424"}}