{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,2]],"date-time":"2026-05-02T06:23:19Z","timestamp":1777702999502,"version":"3.51.4"},"reference-count":44,"publisher":"Oxford University Press (OUP)","issue":"3","license":[{"start":{"date-parts":[[2020,8,13]],"date-time":"2020-08-13T00:00:00Z","timestamp":1597276800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["R01HG009124"],"award-info":[{"award-number":["R01HG009124"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["R01GM126553"],"award-info":[{"award-number":["R01GM126553"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000001","name":"National Science Foundation","doi-asserted-by":"publisher","award":["DMS1712933"],"award-info":[{"award-number":["DMS1712933"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100004543","name":"China Scholarship Council","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100004543","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2021,4,20]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Identifying cis-acting genetic variants associated with gene expression levels\u2014an analysis commonly referred to as expression quantitative trait loci (eQTLs) mapping\u2014is an important first step toward understanding the genetic determinant of gene expression variation. Successful eQTL mapping requires effective control of confounding factors. A common method for confounding effects control in eQTL mapping studies is the probabilistic estimation of expression residual (PEER) analysis. PEER analysis extracts PEER factors to serve as surrogates for confounding factors, which is further included in the subsequent eQTL mapping analysis. However, it is computationally challenging to determine the optimal number of PEER factors used for eQTL mapping. In particular, the standard approach to determine the optimal number of PEER factors examines one number\u00a0at a time and chooses a number that optimizes eQTLs discovery. Unfortunately, this standard approach involves multiple repetitive eQTL mapping procedures that are computationally expensive, restricting its use in large-scale eQTL mapping studies that being collected today.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>Here, we present a simple and computationally scalable alternative, Effect size Correlation for COnfounding determination (ECCO), to determine the optimal number of PEER factors used for eQTL mapping studies. Instead of performing repetitive eQTL mapping, ECCO jointly applies differential expression analysis and Mendelian randomization analysis, leading to substantial computational savings. In simulations and real data applications, we show that ECCO identifies a similar number of PEER factors required for eQTL mapping analysis as the standard approach but is two orders of magnitude faster. The computational scalability of ECCO allows for optimized eQTL discovery across 48 GTEx tissues for the first time, yielding an overall 5.89% power gain on the number of eQTL harboring genes (eGenes) discovered as compared to the previous GTEx recommendation that does not attempt to determine tissue-specific optimal number of PEER factors.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availabilityand implementation<\/jats:title>\n                  <jats:p>Our method is implemented in the ECCO software, which, along with its GTEx mapping results, is freely available at www.xzlab.org\/software.html. All R scripts used in this study are also available at this site.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Supplementary information<\/jats:title>\n                  <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btaa715","type":"journal-article","created":{"date-parts":[[2020,8,6]],"date-time":"2020-08-06T16:37:22Z","timestamp":1596731842000},"page":"296-302","source":"Crossref","is-referenced-by-count":9,"title":["Efficient and effective control of confounding in eQTL mapping studies through joint differential expression and Mendelian randomization analyses"],"prefix":"10.1093","volume":"37","author":[{"given":"Yue","family":"Fan","sequence":"first","affiliation":[{"name":"Key Laboratory of Trace Elements and Endemic Diseases of National Health and Family Planning Commission, School of Public Health, Health Science Center, Xi\u2019an Jiaotong University , Xi\u2019an, Shaanxi 710061, China"},{"name":"Department of Biostatistics, University of Michigan , Ann Arbor, MI 48109, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Huanhuan","family":"Zhu","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of Michigan , Ann Arbor, MI 48109, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yanyi","family":"Song","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of Michigan , Ann Arbor, MI 48109, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Qinke","family":"Peng","sequence":"additional","affiliation":[{"name":"Systems Engineering Institute, Xi\u2019an Jiaotong University , Xi\u2019an, Shaanxi 710049, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xiang","family":"Zhou","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of Michigan , Ann Arbor, MI 48109, USA"},{"name":"Center for Statistical Genetics, University of Michigan , Ann Arbor, MI 48109, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2020,8,13]]},"reference":[{"key":"2023051701213706800_btaa715-B1","doi-asserted-by":"crossref","first-page":"1011","DOI":"10.1016\/j.cell.2015.10.025","article-title":"The molecular taxonomy of primary prostate cancer","volume":"163","author":"Abeshouse","year":"2015","journal-title":"Cell"},{"key":"2023051701213706800_btaa715-B2","doi-asserted-by":"crossref","DOI":"10.1101\/787903","article-title":"The GTEx Consortium atlas of genetic regulatory effects across human tissues","author":"Aguet","year":"2019"},{"key":"2023051701213706800_btaa715-B3","doi-asserted-by":"crossref","first-page":"14","DOI":"10.1101\/gr.155192.113","article-title":"Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals","volume":"24","author":"Battle","year":"2014","journal-title":"Genome Res"},{"key":"2023051701213706800_btaa715-B4","first-page":"443","article-title":"Problems with instrumental variables estimation when the correlation between the instruments and the endogenous explanatory variable is weak","volume":"90","author":"Bound","year":"1995","journal-title":"J. Am. Stat. Assoc"},{"key":"2023051701213706800_btaa715-B5","doi-asserted-by":"crossref","first-page":"512","DOI":"10.1093\/ije\/dyv080","article-title":"Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression","volume":"44","author":"Bowden","year":"2015","journal-title":"Int. J. Epidemiol"},{"key":"2023051701213706800_btaa715-B6","doi-asserted-by":"crossref","first-page":"e1004461","DOI":"10.1371\/journal.pgen.1004461","article-title":"Cis and trans effects of human genomic variants on gene expression","volume":"10","author":"Bryois","year":"2014","journal-title":"PLos Genet"},{"key":"2023051701213706800_btaa715-B7","doi-asserted-by":"crossref","first-page":"755","DOI":"10.1093\/ije\/dyr036","article-title":"Avoiding bias from weak instruments in Mendelian randomization studies","volume":"40","author":"Burgess","year":"2011","journal-title":"Int. J. Epidemiol"},{"key":"2023051701213706800_btaa715-B8","article-title":"Controlling for confounding effects in single Cell RNA sequencing studies using both control and target genes","volume":"7, 1-14","author":"Chen","year":"2017","journal-title":"Sci. Rep"},{"key":"2023051701213706800_btaa715-B9","article-title":"A survey of best practices for RNA-seq data analysis","volume":"17, 13","author":"Conesa","year":"2016","journal-title":"Genome Biol"},{"key":"2023051701213706800_btaa715-B10","doi-asserted-by":"crossref","first-page":"530","DOI":"10.1038\/nature25160","article-title":"Genetic effects on gene expression across human tissues","volume":"553","author":"Consortium","year":"2018","journal-title":"Nature"},{"key":"2023051701213706800_btaa715-B11","doi-asserted-by":"crossref","first-page":"184","DOI":"10.1038\/nrg2537","article-title":"Mapping complex disease traits with global gene expression","volume":"10","author":"Cookson","year":"2009","journal-title":"Nat. Rev. Genet"},{"key":"2023051701213706800_btaa715-B12","doi-asserted-by":"crossref","first-page":"e1003864","DOI":"10.1371\/journal.pgen.1003864","article-title":"Partitioning the heritability of tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture","volume":"9","author":"Davis","year":"2013","journal-title":"PLoS Genet"},{"key":"2023051701213706800_btaa715-B13","doi-asserted-by":"crossref","first-page":"390","DOI":"10.1038\/nature10808","article-title":"DNase I sensitivity QTLs are a major determinant of human expression variation","volume":"482","author":"Degner","year":"2012","journal-title":"Nature"},{"key":"2023051701213706800_btaa715-B14","doi-asserted-by":"crossref","first-page":"539","DOI":"10.1093\/biostatistics\/kxr034","article-title":"Using control genes to correct for unwanted variation in microarray data","volume":"13","author":"Gagnon-Bartsch","year":"2012","journal-title":"Biostatistics"},{"key":"2023051701213706800_btaa715-B15","doi-asserted-by":"crossref","first-page":"1091","DOI":"10.1038\/ng.3367","article-title":"A gene-based association method for mapping traits using reference transcriptome data","volume":"47","author":"Gamazon","year":"2015","journal-title":"Nat. Genet"},{"key":"2023051701213706800_btaa715-B16","doi-asserted-by":"crossref","first-page":"e1004383","DOI":"10.1371\/journal.pgen.1004383","article-title":"Bayesian test for colocalisation between pairs of genetic association studies using summary statistics","volume":"10","author":"Giambartolomei","year":"2014","journal-title":"PLoS Genet"},{"key":"2023051701213706800_btaa715-B17","doi-asserted-by":"crossref","first-page":"575","DOI":"10.1038\/nrg2383","article-title":"The environmental contribution to gene expression profiles","volume":"9","author":"Gibson","year":"2008","journal-title":"Nat. Rev. Genet"},{"key":"2023051701213706800_btaa715-B18","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1038\/ng.3506","article-title":"Integrative approaches for large-scale transcriptome-wide association studies","volume":"48","author":"Gusev","year":"2016","journal-title":"Nat. Genet"},{"key":"2023051701213706800_btaa715-B19","doi-asserted-by":"crossref","first-page":"1245","DOI":"10.1016\/j.ajhg.2016.10.003","article-title":"Colocalization of GWAS and eQTL signals detects target genes","volume":"99","author":"Hormozdiari","year":"2016","journal-title":"Am. J. Hum. Genet"},{"key":"2023051701213706800_btaa715-B20","doi-asserted-by":"crossref","first-page":"506","DOI":"10.1038\/nature12531","article-title":"Transcriptome and genome sequencing uncovers functional variation in humans","volume":"501","author":"Lappalainen","year":"2013","journal-title":"Nature"},{"key":"2023051701213706800_btaa715-B21","doi-asserted-by":"crossref","first-page":"e161","DOI":"10.1371\/journal.pgen.0030161","article-title":"Capturing heterogeneity in gene expression studies by surrogate variable analysis","volume":"3","author":"Leek","year":"2007","journal-title":"PLoS Genet"},{"key":"2023051701213706800_btaa715-B22","doi-asserted-by":"crossref","first-page":"716","DOI":"10.1101\/gr.142521.112","article-title":"A cross-platform analysis of 14,177 expression quantitative trait loci derived from lymphoblastoid cell lines","volume":"23","author":"Liang","year":"2013","journal-title":"Genome Res"},{"key":"2023051701213706800_btaa715-B23","doi-asserted-by":"crossref","DOI":"10.1101\/395970","article-title":"GBAT: a gene-based association method for robust trans-gene regulation detection","author":"Liu","year":"2018"},{"key":"2023051701213706800_btaa715-B24","doi-asserted-by":"crossref","first-page":"e1007755","DOI":"10.1371\/journal.pgen.1007755","article-title":"Functional regulatory mechanism of smooth muscle cell-restricted LMOD1 coronary artery disease locus","volume":"14","author":"Nanda","year":"2018","journal-title":"PLoS Genet"},{"key":"2023051701213706800_btaa715-B25","doi-asserted-by":"crossref","first-page":"1940","DOI":"10.1016\/j.celrep.2017.05.018","article-title":"Effect of human genetic variability on gene expression in dorsal root ganglia and association with pain phenotypes","volume":"19","author":"Parisien","year":"2017","journal-title":"Cell Rep"},{"key":"2023051701213706800_btaa715-B26","doi-asserted-by":"crossref","first-page":"e1001276","DOI":"10.1371\/journal.pgen.1001276","article-title":"Joint genetic analysis of gene expression data with inferred cellular phenotypes","volume":"7","author":"Parts","year":"2011","journal-title":"PLoS Genet"},{"key":"2023051701213706800_btaa715-B27","doi-asserted-by":"crossref","first-page":"768","DOI":"10.1038\/nature08872","article-title":"Understanding mechanisms underlying human gene expression variation with RNA sequencing","volume":"464","author":"Pickrell","year":"2010","journal-title":"Nature"},{"key":"2023051701213706800_btaa715-B28","doi-asserted-by":"crossref","DOI":"10.1038\/s41467-019-10936-0","article-title":"Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits","volume":"10","author":"Porcu","year":"2019","journal-title":"Nat. Commun"},{"key":"2023051701213706800_btaa715-B29","doi-asserted-by":"crossref","first-page":"1220","DOI":"10.1093\/ije\/dym225","article-title":"Cohort profile: the cardiovascular risk in Young Finns Study","volume":"37","author":"Raitakari","year":"2008","journal-title":"Int. J. Epidemiol"},{"key":"2023051701213706800_btaa715-B30","doi-asserted-by":"crossref","first-page":"519","DOI":"10.1126\/science.1249547","article-title":"Polarization of the effects of autoimmune and neurodegenerative risk alleles in leukocytes","volume":"344","author":"Raj","year":"2014","journal-title":"Science"},{"key":"2023051701213706800_btaa715-B31","article-title":"Modelling local gene networks increases power to detect trans-acting genetic effects on gene expression","volume":"17, 33","author":"Rakitsch","year":"2016","journal-title":"Genome Biol"},{"key":"2023051701213706800_btaa715-B32","doi-asserted-by":"crossref","first-page":"1860","DOI":"10.12688\/f1000research.17145.1","article-title":"False positives in trans-eQTL and co-expression analyses arising from RNA-sequencing alignment errors","volume":"7","author":"Saha","year":"2018","journal-title":"F1000Research"},{"key":"2023051701213706800_btaa715-B33","doi-asserted-by":"crossref","first-page":"1353","DOI":"10.1093\/bioinformatics\/bts163","article-title":"Matrix eQTL: ultra fast eQTL analysis via large matrix operations","volume":"28","author":"Shabalin","year":"2012","journal-title":"Bioinformatics"},{"key":"2023051701213706800_btaa715-B34","doi-asserted-by":"crossref","first-page":"557","DOI":"10.2307\/2171753","article-title":"Instrumental variables regression with weak instruments","volume":"65","author":"Staiger","year":"1997","journal-title":"Econometrica"},{"key":"2023051701213706800_btaa715-B35","doi-asserted-by":"crossref","first-page":"1895","DOI":"10.2337\/db11-1378","article-title":"Hyperglycemia and a common variant of GCKR are associated with the levels of eight amino acids in 9,369 Finnish men","volume":"61","author":"Stan\u010d\u00e1kov\u00e1","year":"2012","journal-title":"Diabetes"},{"key":"2023051701213706800_btaa715-B36","doi-asserted-by":"crossref","first-page":"e1000770","DOI":"10.1371\/journal.pcbi.1000770","article-title":"A Bayesian framework to account for complex non-genetic factors in gene expression levels greatly increases power in eQTL studies","volume":"6","author":"Stegle","year":"2010","journal-title":"PLoS Comput. Biol"},{"key":"2023051701213706800_btaa715-B37","doi-asserted-by":"crossref","first-page":"521","DOI":"10.1016\/j.ajhg.2014.10.001","article-title":"Cross-tissue and tissue-specific eQTLs: partitioning the heritability of a complex trait","volume":"95","author":"Torres","year":"2014","journal-title":"Am. J. Hum. Genet"},{"key":"2023051701213706800_btaa715-B38","article-title":"The genetic architecture of gene expression levels in wild baboons","volume":"4, e04729","author":"Tung","year":"2015","journal-title":"Elife"},{"key":"2023051701213706800_btaa715-B39","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1214\/aoms\/1177731868","article-title":"The fitting of straight lines if both variables are subject to error","volume":"11","author":"Wald","year":"1940","journal-title":"Ann. Math. Stat"},{"key":"2023051701213706800_btaa715-B40","doi-asserted-by":"crossref","first-page":"430","DOI":"10.1038\/ng.2951","article-title":"Heritability and genomics of gene expression in peripheral blood","volume":"46","author":"Wright","year":"2014","journal-title":"Nat. Genet"},{"key":"2023051701213706800_btaa715-B41","first-page":"1","article-title":"Testing and controlling for horizontal pleiotropy with the probabilistic Mendelian randomization in transcriptome-wide association studies","author":"Yuan","year":"2020"},{"key":"2023051701213706800_btaa715-B42","doi-asserted-by":"crossref","first-page":"225","DOI":"10.1186\/s12916-019-1448-9","article-title":"Causal association of type 2 diabetes with amyotrophic lateral sclerosis: new evidence from Mendelian randomization using GWAS summary statistics","volume":"17","author":"Zeng","year":"2019","journal-title":"BMC Med"},{"key":"2023051701213706800_btaa715-B43","doi-asserted-by":"crossref","DOI":"10.1038\/s41467-017-00470-2","article-title":"Non-parametric genetic prediction of complex traits with latent Dirichlet process regression models","volume":"8","author":"Zeng","year":"2017","journal-title":"Nat. Commun"},{"key":"2023051701213706800_btaa715-B44","doi-asserted-by":"crossref","first-page":"688","DOI":"10.1093\/hmg\/ddy384","article-title":"Causal effects of blood lipids on amyotrophic lateral sclerosis: a Mendelian randomization study","volume":"28","author":"Zeng","year":"2019","journal-title":"Hum. Mol. Genet"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btaa715\/34111040\/btaa715.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/37\/3\/296\/50326094\/btaa715.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/37\/3\/296\/50326094\/btaa715.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,5,17]],"date-time":"2023-05-17T01:22:24Z","timestamp":1684286544000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/37\/3\/296\/5892253"}},"subtitle":[],"editor":[{"given":"Peter","family":"Robinson","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2020,8,13]]},"references-count":44,"journal-issue":{"issue":"3","published-print":{"date-parts":[[2021,4,20]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btaa715","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2021,2,1]]},"published":{"date-parts":[[2020,8,13]]}}}