{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,11,11]],"date-time":"2025-11-11T13:29:45Z","timestamp":1762867785716,"version":"3.37.3"},"reference-count":35,"publisher":"Oxford University Press (OUP)","issue":"23","license":[{"start":{"date-parts":[[2019,5,22]],"date-time":"2019-05-22T00:00:00Z","timestamp":1558483200000},"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\/501100012166","name":"National Key R&D Program of China","doi-asserted-by":"crossref","award":["2017YFC0908600"],"award-info":[{"award-number":["2017YFC0908600"]}],"id":[{"id":"10.13039\/501100012166","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/501100001809","name":"Natural Science Foundation of China","doi-asserted-by":"crossref","award":["31571355"],"award-info":[{"award-number":["31571355"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"crossref"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2019,12,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Evolution of gene across species is usually subject to the stabilizing selection to maintain the optimal expression level. While it is generally accepted that the resulting expression conservation may vary considerably among genes, statistically reliable estimation remains challenging, due to few species included in current comparative RNA-seq data with high number of unknown parameters.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>In this paper, we develop a gamma distribution model to describe how the strength of expression conservation (denoted by W) varies among genes. Given the high throughput RNA-seq datasets from multiple species, we then formulate an empirical Bayesian procedure to estimate W for each gene. Our case studies showed that those W-estimates are useful to study the evolutionary pattern of expression conservation.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>Our method has been implemented in the R-package software, TreeExp, which is publically available at Github develop site https:\/\/github.com\/hr1912\/TreeExp. It involves three functions: estParaGamma, estParaQ and estParaWBayesian. The manual for software TreeExp is available at https:\/\/github.com\/hr1912\/TreeExp\/tree\/master\/vignettes. For any question, one may contact Dr Hang Ruan (Hang.Ruan@uth.tmc.edu).<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btz405","type":"journal-article","created":{"date-parts":[[2019,5,16]],"date-time":"2019-05-16T11:46:16Z","timestamp":1558007176000},"page":"5030-5038","source":"Crossref","is-referenced-by-count":10,"title":["Estimating the strength of expression conservation from high throughput RNA-seq data"],"prefix":"10.1093","volume":"35","author":[{"given":"Xun","family":"Gu","sequence":"first","affiliation":[{"name":"Department of Genetics, Development, and Cell Biology, Iowa State University , Ames, IA, USA"}]},{"given":"Hang","family":"Ruan","sequence":"additional","affiliation":[{"name":"MOE Key Laboratory of Contemporary Anthropology, Fudan University , Shanghai, China"},{"name":"Department of Biochemistry and Molecular Biology, McGovern Medical School , The University of Texas Health Science Center at Houston, Houston, TX, USA"}]},{"given":"Jingwen","family":"Yang","sequence":"additional","affiliation":[{"name":"MOE Key Laboratory of Contemporary Anthropology, Fudan University , Shanghai, China"},{"name":"Human Phenome Institute, Fudan University , Shanghai, China"}]}],"member":"286","published-online":{"date-parts":[[2019,5,22]]},"reference":[{"key":"2023013108321827400_btz405-B1","doi-asserted-by":"crossref","first-page":"1133","DOI":"10.1073\/pnas.0812009106","article-title":"Optimization of gene expression by natural selection","volume":"106","author":"Bedford","year":"2009","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013108321827400_btz405-B2","first-page":"85","article-title":"Similarities and differences in genome-wide expression data of six organisms","volume":"2","author":"Bergmann","year":"2004","journal-title":"PLoS Biol"},{"key":"2023013108321827400_btz405-B3","first-page":"380","article-title":"No evidence that protein noise-induced epigenetic epistasis constrains gene expression evolution","volume":"34","author":"Boross","year":"2017","journal-title":"Mol. Biol. Evol"},{"key":"2023013108321827400_btz405-B4","doi-asserted-by":"crossref","first-page":"343","DOI":"10.1038\/nature10532","article-title":"The evolution of gene expression levels in mammalian organs","volume":"478","author":"Brawand","year":"2011","journal-title":"Nature"},{"key":"2023013108321827400_btz405-B5","doi-asserted-by":"crossref","first-page":"6907","DOI":"10.1073\/pnas.0406882102","article-title":"The THAP domain of THAP1 is a large C2CH module with zinc-dependent sequence-specific DNA-binding activity","volume":"102","author":"Clouaire","year":"2005","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013108321827400_btz405-B6","doi-asserted-by":"crossref","first-page":"372","DOI":"10.1016\/j.tig.2007.04.003","article-title":"MicroRNA regulation and interspecific variation of gene expression","volume":"23","author":"Cui","year":"2007","journal-title":"Trends Genet"},{"key":"2023013108321827400_btz405-B7","doi-asserted-by":"crossref","first-page":"588","DOI":"10.1038\/nature03842","article-title":"Optimality and evolutionary tuning of the expression level of a protein","volume":"436","author":"Dekel","year":"2005","journal-title":"Nature"},{"key":"2023013108321827400_btz405-B8","doi-asserted-by":"crossref","first-page":"1915","DOI":"10.1534\/genetics.104.036871","article-title":"Mechanisms of haploinsufficiency revealed by genome-wide profiling in yeast","volume":"169","author":"Deutschbauer","year":"2005","journal-title":"Genetics"},{"key":"2023013108321827400_btz405-B9","doi-asserted-by":"crossref","first-page":"340","DOI":"10.1126\/science.1068996","article-title":"Intra- and interspecific variation in primate gene expression patterns","volume":"296","author":"Enard","year":"2002","journal-title":"Science"},{"key":"2023013108321827400_btz405-B10","doi-asserted-by":"crossref","first-page":"242","DOI":"10.1038\/nature04559","article-title":"Expression profiling in primates reveals a rapid evolution of human transcription factors","volume":"440","author":"Gilad","year":"2006","journal-title":"Nature"},{"key":"2023013108321827400_btz405-B11","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1016\/S0168-9525(02)00040-9","article-title":"Induced gene expression in human brain after the split from chimpanzee","volume":"19","author":"Gu","year":"2003","journal-title":"Trends Genet"},{"key":"2023013108321827400_btz405-B12","doi-asserted-by":"crossref","first-page":"243","DOI":"10.1093\/bib\/bbv035","article-title":"Statistical detection of differentially expressed genes based on RNA-seq: from biological to phylogenetic replicates","volume":"17","author":"Gu","year":"2016","journal-title":"Brief. Bioinform"},{"key":"2023013108321827400_btz405-B13","doi-asserted-by":"crossref","first-page":"249","DOI":"10.1093\/bib\/bbv041","article-title":"Understanding tissue expression evolution: from expression phylogeny to phylogenetic network","volume":"17","author":"Gu","year":"2016","journal-title":"Brief. Bioinform"},{"key":"2023013108321827400_btz405-B14","doi-asserted-by":"crossref","first-page":"2779","DOI":"10.1073\/pnas.0610797104","article-title":"Tissue-driven hypothesis of genomic evolution and sequence-expression correlations","volume":"104","author":"Gu","year":"2007","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023013108321827400_btz405-B15","doi-asserted-by":"crossref","first-page":"1404","DOI":"10.2307\/2410878","article-title":"Translating between microevolutionary process and macroevolutionary patterns: the correlation structure of interspecific data","volume":"50","author":"Hansen","year":"1996","journal-title":"Evolution"},{"key":"2023013108321827400_btz405-B16","doi-asserted-by":"crossref","first-page":"222","DOI":"10.1016\/j.semcdb.2011.12.004","article-title":"The evolution of gene expression and the transcriptome-phenotype relationship","volume":"23","author":"Harrison","year":"2012","journal-title":"Semin. Cell Dev. Biol"},{"key":"2023013108321827400_btz405-B17","doi-asserted-by":"crossref","first-page":"682","DOI":"10.1371\/journal.pbio.0020132","article-title":"A neutral model of transcriptome evolution","volume":"2","author":"Khaitovich","year":"2004","journal-title":"PLoS Biol"},{"key":"2023013108321827400_btz405-B18","doi-asserted-by":"crossref","first-page":"1850","DOI":"10.1126\/science.1108296","article-title":"Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees","volume":"309","author":"Khaitovich","year":"2005","journal-title":"Science"},{"key":"2023013108321827400_btz405-B19","doi-asserted-by":"crossref","first-page":"107","DOI":"10.1126\/science.1090005","article-title":"Evolution at 2 Levels in humans and chimpanzees","volume":"188","author":"King","year":"1975","journal-title":"Science"},{"key":"2023013108321827400_btz405-B20","doi-asserted-by":"crossref","first-page":"314","DOI":"10.2307\/2407703","article-title":"Natural selection and random genetic drift in phenotypic evolution","volume":"30","author":"Lande","year":"1976","journal-title":"Evolution"},{"key":"2023013108321827400_btz405-B21","doi-asserted-by":"crossref","first-page":"168","DOI":"10.1038\/nrg3404","article-title":"Genotype to phenotype: lessons from model organisms for human genetics","volume":"14","author":"Lehner","year":"2013","journal-title":"Nat. Rev. Genet"},{"key":"2023013108321827400_btz405-B22","first-page":"126","article-title":"Rates of divergence in gene expression profiles of primates, mice, and flies: stabilizing selection and variability among functional categories","volume":"59","author":"Lemos","year":"2005","journal-title":"Evolution"},{"key":"2023013108321827400_btz405-B23","doi-asserted-by":"crossref","first-page":"194","DOI":"10.1038\/nature01771","article-title":"Dosage sensitivity and the evolution of gene families in yeast","volume":"424","author":"Papp","year":"2003","journal-title":"Nature"},{"key":"2023013108321827400_btz405-B24","doi-asserted-by":"crossref","first-page":"645","DOI":"10.1038\/msb.2013.2","article-title":"Epigenetic epistatic interactions constrain the evolution of gene expression","volume":"9","author":"Park","year":"2013","journal-title":"Mol. Syst. Biol"},{"key":"2023013108321827400_btz405-B25","doi-asserted-by":"crossref","first-page":"138","DOI":"10.1038\/ng1086","article-title":"Evolution of gene expression in the Drosophila melanogaster subgroup","volume":"33","author":"Rifkin","year":"2003","journal-title":"Nat. Genet"},{"key":"2023013108321827400_btz405-B26","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"edgeR: a Bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2023013108321827400_btz405-B27","doi-asserted-by":"crossref","first-page":"201","DOI":"10.1093\/molbev\/mst190","article-title":"Modeling gene expression evolution with an extended Ornstein\u2013Uhlenbeck process accounting for within-species variation","volume":"31","author":"Rohlfs","year":"2014","journal-title":"Mol. Biol. Evol"},{"key":"2023013108321827400_btz405-B28","doi-asserted-by":"crossref","first-page":"695","DOI":"10.1093\/sysbio\/syv042","article-title":"Phylogenetic ANOVA: the expression variance and evolution model for quantitative trait evolution","volume":"64","author":"Rohlfs","year":"2015","journal-title":"Syst. Biol"},{"key":"2023013108321827400_btz405-B29","doi-asserted-by":"crossref","first-page":"394","DOI":"10.1002\/jez.b.22707","article-title":"TREEEXP1.0: R package for analyzing expression evolution based on RNA-seq data","volume":"326","author":"Ruan","year":"2016","journal-title":"J. Exp. Zool. Part B"},{"key":"2023013108321827400_btz405-B30","doi-asserted-by":"crossref","first-page":"319","DOI":"10.1016\/j.molcel.2005.12.011","article-title":"Mapping pathways and phenotypes by systematic gene overexpression","volume":"21","author":"Sopko","year":"2006","journal-title":"Mol. Cell"},{"key":"2023013108321827400_btz405-B31","doi-asserted-by":"crossref","first-page":"109","DOI":"10.1016\/j.tig.2007.12.004","article-title":"Evolution of gene sequence and gene expression are not correlated in yeast","volume":"24","author":"Tirosh","year":"2008","journal-title":"Trends Genet"},{"key":"2023013108321827400_btz405-B32","doi-asserted-by":"crossref","first-page":"830","DOI":"10.1038\/ng1819","article-title":"A genetic signature of interspecies variations in gene expression","volume":"38","author":"Tirosh","year":"2006","journal-title":"Nat. Genet"},{"key":"2023013108321827400_btz405-B33","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1038\/nrg2484","article-title":"RNA-seq: a revolutionary tool for transcriptomics","volume":"10","author":"Wang","year":"2009","journal-title":"Nat. Rev. Genet"},{"key":"2023013108321827400_btz405-B34","doi-asserted-by":"crossref","first-page":"e34935","DOI":"10.1371\/journal.pone.0034935","article-title":"A selection index for gene expression evolution and its application to the divergence between humans and chimpanzees","volume":"7","author":"Warnefors","year":"2012","journal-title":"PLoS One"},{"key":"2023013108321827400_btz405-B35","doi-asserted-by":"crossref","first-page":"2893","DOI":"10.1093\/molbev\/msr116","article-title":"Predominant gain of promoter TATA box after gene duplication associated with stress responses","volume":"28","author":"Zou","year":"2011","journal-title":"Mol. Biol. Evol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btz405\/28928401\/btz405.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/23\/5030\/48978261\/bioinformatics_35_23_5030.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/35\/23\/5030\/48978261\/bioinformatics_35_23_5030.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,1,31]],"date-time":"2023-01-31T17:38:51Z","timestamp":1675186731000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/35\/23\/5030\/5494699"}},"subtitle":[],"editor":[{"given":"Inanc","family":"Birol","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2019,5,22]]},"references-count":35,"journal-issue":{"issue":"23","published-print":{"date-parts":[[2019,12,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btz405","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"type":"print","value":"1367-4803"},{"type":"electronic","value":"1367-4811"}],"subject":[],"published-other":{"date-parts":[[2019,12,1]]},"published":{"date-parts":[[2019,5,22]]}}}