{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,8]],"date-time":"2026-10-08T19:24:46Z","timestamp":1791487486609,"version":"4.3.4"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1007784","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,3,3]],"date-time":"2021-03-03T00:00:00Z","timestamp":1614729600000}}],"reference-count":44,"publisher":"Public Library of Science (PLoS)","issue":"2","license":[{"start":{"date-parts":[[2021,2,19]],"date-time":"2021-02-19T00:00:00Z","timestamp":1613692800000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>Rare variants are thought to play an important role in the etiology of complex diseases and may explain a significant fraction of the missing heritability in genetic disease studies. Next-generation sequencing facilitates the association of rare variants in coding or regulatory regions with complex diseases in large cohorts at genome-wide scale. However, rare variant association studies (RVAS) still lack power when cohorts are small to medium-sized and if genetic variation explains a small fraction of phenotypic variance. Here we present a novel Bayesian rare variant Association Test using Integrated Nested Laplace Approximation (BATI). Unlike existing RVAS tests, BATI allows integration of individual or variant-specific features as covariates, while efficiently performing inference based on full model estimation. We demonstrate that BATI outperforms established RVAS methods on realistic, semi-synthetic whole-exome sequencing cohorts, especially when using meaningful biological context, such as functional annotation. We show that BATI achieves power above 70% in scenarios in which competing tests fail to identify risk genes, e.g. when risk variants in sum explain less than 0.5% of phenotypic variance. We have integrated BATI, together with five existing RVAS tests in the \u2018Rare Variant Genome Wide Association Study\u2019 (rvGWAS) framework for data analyzed by whole-exome or whole genome sequencing. rvGWAS supports rare variant association for genes or any other biological unit such as promoters, while allowing the analysis of essential functionalities like quality control or filtering. Applying rvGWAS to a Chronic Lymphocytic Leukemia study we identified eight candidate predisposition genes, including EHMT2 and COPS7A.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1007784","type":"journal-article","created":{"date-parts":[[2021,2,20]],"date-time":"2021-02-20T00:11:28Z","timestamp":1613779888000},"page":"e1007784","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":14,"title":["Efficient and flexible Integration of variant characteristics in rare variant association studies using integrated nested Laplace approximation"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-7965-7675","authenticated-orcid":true,"given":"Hana","family":"Susak","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8835-6890","authenticated-orcid":true,"given":"Laura","family":"Serra-Saurina","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9075-4276","authenticated-orcid":true,"given":"German","family":"Demidov","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5006-8140","authenticated-orcid":true,"given":"Raquel","family":"Rabionet","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2265-0824","authenticated-orcid":true,"given":"Laura","family":"Dom\u00e8nech","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8989-4220","authenticated-orcid":true,"given":"Mattia","family":"Bosio","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7857-0623","authenticated-orcid":true,"given":"Francesc","family":"Muyas","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xavier","family":"Estivill","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5416-3177","authenticated-orcid":true,"given":"Ge\u00f2rgia","family":"Escaram\u00eds","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7416-9568","authenticated-orcid":true,"given":"Stephan","family":"Ossowski","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2021,2,19]]},"reference":[{"key":"pcbi.1007784.ref001","article-title":"Multiple rare alleles contribute to low plasma levels of HDL cholesterol","author":"JC Cohen","year":"2004","journal-title":"Science (80-)"},{"key":"pcbi.1007784.ref002","article-title":"Targeted resequencing identifies PTCH1 as a major contributor to ocular developmental anomalies and extends the SOX2 regulatory network","author":"N Chassaing","year":"2016","journal-title":"Genome Res"},{"key":"pcbi.1007784.ref003","article-title":"De Novo and Rare Variants at Multiple Loci Support the Oligogenic Origins of Atrioventricular Septal Heart Defects","author":"JR Priest","year":"2016","journal-title":"PLoS Genet"},{"key":"pcbi.1007784.ref004","article-title":"Systematic functional testing of rare variants: Contributions of CFI to age-related macular degeneration","author":"PL Tan","year":"2017","journal-title":"Investig Ophthalmol Vis Sci"},{"key":"pcbi.1007784.ref005","article-title":"Evolution and Functional Impact of Rare Coding Variation from Deep Sequencing of Human Exomes Broad GO, Seattle GO, on behalf of the NHLBI Exome Sequencing Project","author":"JA Tennessen","year":"2012","journal-title":"Science (80-)"},{"key":"pcbi.1007784.ref006","article-title":"An abundance of rare functional variants in 202 drug target genes sequenced in 14,002 people","author":"MR Nelson","year":"2012","journal-title":"Science (80-)"},{"key":"pcbi.1007784.ref007","article-title":"Searching for missing heritability: Designing rare variant association studies","author":"O Zuk","year":"2014","journal-title":"Proc Natl Acad Sci"},{"key":"pcbi.1007784.ref008","article-title":"Methods for Detecting Associations with Rare Variants for Common Diseases: Application to Analysis of Sequence Data","author":"B Li","year":"2008","journal-title":"Am J Hum Genet"},{"key":"pcbi.1007784.ref009","article-title":"Pooled Association Tests for Rare Variants in Exon-Resequencing Studies","author":"AL Price","year":"2010","journal-title":"Am J Hum Genet"},{"key":"pcbi.1007784.ref010","article-title":"A groupwise association test for rare mutations using a weighted sum statistic","author":"BE Madsen","year":"2009","journal-title":"PLoS Genet"},{"key":"pcbi.1007784.ref011","article-title":"A novel adaptive method for the analysis of next-generation sequencing data to detect complex trait associations with rare variants due to gene main effects and interactions","author":"DJ Liu","year":"2010","journal-title":"PLoS Genet"},{"key":"pcbi.1007784.ref012","article-title":"Rare-variant association testing for sequencing data with the sequence kernel association test","author":"MC Wu","year":"2011","journal-title":"Am J Hum Genet"},{"key":"pcbi.1007784.ref013","article-title":"Testing for an unusual distribution of rare variants","author":"BM Neale","year":"2011","journal-title":"PLoS Genet"},{"key":"pcbi.1007784.ref014","article-title":"The Power of Gene-Based Rare Variant Methods to Detect Disease-Associated Variation and Test Hypotheses About Complex Disease","author":"L Moutsianas","year":"2015","journal-title":"PLoS Genet"},{"key":"pcbi.1007784.ref015","article-title":"Optimal tests for rare variant effects in sequencing association studies.","author":"S Lee","year":"2012","journal-title":"Biostatistics"},{"key":"pcbi.1007784.ref016","article-title":"Hierarchical bayesian model for rare variant association analysis integrating genotype uncertainty in human sequence data","author":"L He","year":"2015","journal-title":"Genet Epidemiol"},{"key":"pcbi.1007784.ref017","article-title":"A Unified Mixed-Effects Model for Rare-Variant Association in Sequencing Studies","author":"J Sun","year":"2013","journal-title":"Genet Epidemiol"},{"key":"pcbi.1007784.ref018","article-title":"Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations","author":"H Rue","year":"2009","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"pcbi.1007784.ref019","article-title":"Markov Chain Monte Carlo Convergence Diagnostics: A Comparative Review","author":"MK Cowles","year":"1996","journal-title":"J Am Stat Assoc"},{"key":"pcbi.1007784.ref020","article-title":"Approximate Bayesian inference for hierarchical Gaussian Markov random field models","author":"H\u00e5 Rue","year":"2007","journal-title":"J Stat Plan Inference"},{"key":"pcbi.1007784.ref021","article-title":"Bayesian computing with INLA: New features","author":"TG Martins","year":"2013","journal-title":"Comput Stat Data Anal"},{"key":"pcbi.1007784.ref022","article-title":"Spatial and spatio-temporal models with R-INLA","author":"M Blangiardo","year":"2013","journal-title":"Spatial and Spatio-temporal Epidemiology"},{"key":"pcbi.1007784.ref023","article-title":"Accurate approximations for posterior moments and marginal densities","author":"L Tierney","year":"1986","journal-title":"J Am Stat Assoc"},{"key":"pcbi.1007784.ref024","author":"DJ Spiegelhalter","year":"2002","journal-title":"Bayesian measures of model complexity and fitJ R Stat Soc Ser B Stat Methodol"},{"key":"pcbi.1007784.ref025","article-title":"Animal Models and Integrated Nested Laplace Approximations.","author":"AM Holand","year":"2013","journal-title":"G3&#58; Genes|Genomes|Genetics"},{"key":"pcbi.1007784.ref026","article-title":"The Human Gene Mutation Database: Building a comprehensive mutation repository for clinical and molecular genetics, diagnostic testing and personalized genomic medicine","author":"PD Stenson","year":"2014","journal-title":"Human Genetics"},{"key":"pcbi.1007784.ref027","article-title":"A map of human genome variation from population scale sequencing","author":"G Abecasis","year":"2010","journal-title":"Nature"},{"key":"pcbi.1007784.ref028","article-title":"A global reference for human genetic variation","author":"1000 Genomes Project Consortium","year":"2015","journal-title":"Nature"},{"key":"pcbi.1007784.ref029","article-title":"Epidemiology in practice: Case-control studies","author":"S Lewallen","year":"1998","journal-title":"Community Eye Heal J"},{"key":"pcbi.1007784.ref030","article-title":"Evaluating the heritability explained by known susceptibility variants: A survey of ten complex diseases","author":"HC So","year":"2011","journal-title":"Genet Epidemiol"},{"key":"pcbi.1007784.ref031","volume-title":"Global Cancer Observatory: Cancer Today.","author":"J Ferlay","year":"2018"},{"key":"pcbi.1007784.ref032","article-title":"International network of cancer genome projects","author":"TJ Hudson","year":"2010","journal-title":"Nature"},{"key":"pcbi.1007784.ref033","article-title":"Non-coding recurrent mutations in chronic lymphocytic leukaemia","author":"XS Puente","year":"2015","journal-title":"Nature"},{"key":"pcbi.1007784.ref034","article-title":"Allele balance bias identifies systematic genotyping errors and false disease associations","author":"F Muyas","year":"2019","journal-title":"Hum Mutat"},{"issue":"5","key":"pcbi.1007784.ref035","doi-asserted-by":"crossref","first-page":"955","DOI":"10.1007\/s10637-018-0613-x","article-title":"GLP overexpression is associated with poor prognosis in Chronic Lymphocytic Leukemia and its inhibition induces leukemic cell death","volume":"36","author":"JC Alves-Silva","year":"2018","journal-title":"Invest New Drugs"},{"key":"pcbi.1007784.ref036","article-title":"An integrative approach to predicting the functional effects of small indels in non-coding regions of the human genome","author":"M Ferlaino","year":"2017","journal-title":"BMC Bioinformatics"},{"key":"pcbi.1007784.ref037","article-title":"M-CAP eliminates a majority of variants of uncertain significance in clinical exomes at high sensitivity","author":"KA Jagadeesh","year":"2016","journal-title":"Nat Genet"},{"key":"pcbi.1007784.ref038","article-title":"REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants","author":"NM Ioannidis","year":"2016","journal-title":"Am J Hum Genet"},{"key":"pcbi.1007784.ref039","article-title":"Spatial modelling of lupus incidence over 40 years with changes in census areas","author":"Y Li","year":"2012","journal-title":"J R Stat Soc Ser C Appl Stat"},{"key":"pcbi.1007784.ref040","article-title":"Direct fitting of dynamic models using integrated nested Laplace approximations\u2014INLA","author":"R Ruiz-C\u00e1rdenas","year":"2012","journal-title":"Comput Stat Data Anal"},{"key":"pcbi.1007784.ref041","article-title":"Estimating stochastic volatility models using integrated nested laplace approximations","author":"S Martino","year":"2011","journal-title":"Eur J Financ"},{"key":"pcbi.1007784.ref042","article-title":"Sensitivity analysis in Bayesian generalized linear mixed models for binary data","author":"M Roos","year":"2011","journal-title":"Bayesian Anal"},{"key":"pcbi.1007784.ref043","article-title":"Using integrated nested Laplace approximations for the evaluation of veterinary surveillance data from Switzerland: A case-study","author":"B Schr\u00f6dle","year":"2011","journal-title":"J R Stat Soc Ser C Appl Stat"},{"key":"pcbi.1007784.ref044","article-title":"Bayesian bivariate meta-analysis of diagnostic test studies using integrated nested Laplace approximations","author":"M Paul","year":"2010","journal-title":"Stat Med"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1007784","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,3,3]],"date-time":"2021-03-03T00:00:00Z","timestamp":1614729600000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1007784","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,3,3]],"date-time":"2021-03-03T14:48:58Z","timestamp":1614782938000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1007784"}},"subtitle":[],"editor":[{"given":"Yue","family":"Li","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2021,2,19]]},"references-count":44,"journal-issue":{"issue":"2","published-online":{"date-parts":[[2021,2,19]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1007784","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2020.03.12.988584","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,2,19]]}}}