{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,25]],"date-time":"2026-08-25T22:15:03Z","timestamp":1787696103463,"version":"build-2784847793"},"reference-count":34,"publisher":"Public Library of Science (PLoS)","issue":"4","license":[{"start":{"date-parts":[[2021,4,16]],"date-time":"2021-04-16T00:00:00Z","timestamp":1618531200000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100014535","name":"Mayo Clinic Center for Individualized Medicine","doi-asserted-by":"crossref","id":[{"id":"10.13039\/100014535","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/100014535","name":"Mayo Clinic Center for Individualized Medicine","doi-asserted-by":"crossref","id":[{"id":"10.13039\/100014535","id-type":"DOI","asserted-by":"crossref"}]},{"name":"Illinois Alliance Fellowships for Technology-Based Healthcare Research program"},{"DOI":"10.13039\/501100008982","name":"National Science Foundation","doi-asserted-by":"publisher","award":["OCI-0725070 and ACI-1238993"],"award-info":[{"award-number":["OCI-0725070 and ACI-1238993"]}],"id":[{"id":"10.13039\/501100008982","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100008982","name":"National Science Foundation","doi-asserted-by":"publisher","award":["OCI-0725070 and ACI-1238993"],"award-info":[{"award-number":["OCI-0725070 and ACI-1238993"]}],"id":[{"id":"10.13039\/501100008982","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100008982","name":"National Science Foundation","doi-asserted-by":"publisher","award":["OCI-0725070 and ACI-1238993"],"award-info":[{"award-number":["OCI-0725070 and ACI-1238993"]}],"id":[{"id":"10.13039\/501100008982","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["www.plosone.org"],"crossmark-restriction":false},"short-container-title":["PLoS ONE"],"abstract":"<jats:p>Genetic studies have shifted to sequencing-based rare variants discovery after decades of success in identifying common disease variants by Genome-Wide Association Studies using Single Nucleotide Polymorphism chips. Sequencing-based studies require large sample sizes for statistical power and therefore often inadvertently introduce batch effects because samples are typically collected, processed, and sequenced at multiple centers. Conventionally, batch effects are first detected and visualized using Principal Components Analysis and then controlled by including batch covariates in the disease association models. For sequencing-based genetic studies, because all variants included in the association analyses have passed sequencing-related quality control measures, this conventional approach treats every variant as equal and ignores the substantial differences still remaining in variant qualities and characteristics such as genotype quality scores, alternative allele fractions (fraction of reads supporting alternative allele at a variant position) and sequencing depths. In the Alzheimer\u2019s Disease Sequencing Project (ADSP) exome dataset of 9,904 cases and controls, we discovered hidden variant-level differences between sample batches of three sequencing centers and two exome capture kits. Although sequencing centers were included as a covariate in our association models, we observed differences at the variant level in genotype quality and alternative allele fraction between samples processed by different exome capture kits that significantly impacted both the confidence of variant detection and the identification of disease-associated variants. Furthermore, we found that a subset of top disease-risk variants came exclusively from samples processed by one exome capture kit that was more effective at capturing the alternative alleles compared to the other kit. Our findings highlight the importance of additional variant-level quality control for large sequencing-based genetic studies. More importantly, we demonstrate that automatically filtering out variants with batch differences may lead to false negatives if the batch discordances come largely from quality differences and if the batch-specific variants have better quality.<\/jats:p>","DOI":"10.1371\/journal.pone.0249305","type":"journal-article","created":{"date-parts":[[2021,4,16]],"date-time":"2021-04-16T17:41:06Z","timestamp":1618594866000},"page":"e0249305","update-policy":"https:\/\/doi.org\/10.1371\/journal.pone.corrections_policy","source":"Crossref","is-referenced-by-count":9,"title":["Impact of variant-level batch effects on identification of genetic risk factors in large sequencing studies"],"prefix":"10.1371","volume":"16","author":[{"given":"Daniel P.","family":"Wickland","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yingxue","family":"Ren","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jason P.","family":"Sinnwell","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1783-4453","authenticated-orcid":true,"given":"Joseph S.","family":"Reddy","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Cyril","family":"Pottier","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Vivekananda","family":"Sarangi","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Minerva M.","family":"Carrasquillo","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Owen A.","family":"Ross","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Steven G.","family":"Younkin","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Nil\u00fcfer","family":"Ertekin-Taner","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Rosa","family":"Rademakers","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Matthew E.","family":"Hudson","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Liudmila Sergeevna","family":"Mainzer","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Joanna M.","family":"Biernacka","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-8896-2647","authenticated-orcid":true,"given":"Yan W.","family":"Asmann","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2021,4,16]]},"reference":[{"issue":"6","key":"pone.0249305.ref001","doi-asserted-by":"crossref","first-page":"498","DOI":"10.1016\/j.tibtech.2017.02.012","article-title":"Why batch effects matter in omics data, and how to avoid them","volume":"35","author":"WW Goh","year":"2017","journal-title":"Trends Biotechnol"},{"key":"pone.0249305.ref002","first-page":"207","article-title":"Repurposing germline exomes of the Cancer Genome Atlas demands a cautious approach and sample-specific variant filtering","volume":"21","author":"A Koire","year":"2016","journal-title":"Pacific Symp Biocomput"},{"issue":"1","key":"pone.0249305.ref003","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12885-019-5994-5","article-title":"Substantial batch effects in TCGA exome sequences undermine pan-cancer analysis of germline variants","volume":"19","author":"R Rasnic","year":"2019","journal-title":"BMC Cancer"},{"issue":"1","key":"pone.0249305.ref004","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12864-017-3770-y","article-title":"Pan-cancer analysis reveals technical artifacts in TCGA germline variant calls","volume":"18","author":"AR Buckley","year":"2017","journal-title":"BMC Genomics"},{"issue":"10","key":"pone.0249305.ref005","first-page":"1","article-title":"Whole-exome sequencing capture kit biases yield false negative mutation calls in TCGA cohorts","volume":"13","author":"VG Wang","year":"2018","journal-title":"PLoS One"},{"issue":"5","key":"pone.0249305.ref006","doi-asserted-by":"crossref","DOI":"10.1212\/NXG.0000000000000194","article-title":"Clinical\/Scientific Notes: The Alzheimer\u2019s disease sequencing project: Study design and sample selection","volume":"3","author":"GW Beecham","year":"2017","journal-title":"Neurol Genet"},{"issue":"12","key":"pone.0249305.ref007","doi-asserted-by":"crossref","first-page":"2074","DOI":"10.1371\/journal.pgen.0020190","article-title":"Population structure and eigenanalysis","volume":"2","author":"N Patterson","year":"2006","journal-title":"PLoS Genet"},{"issue":"8","key":"pone.0249305.ref008","doi-asserted-by":"crossref","first-page":"904","DOI":"10.1038\/ng1847","article-title":"Principal components analysis corrects for stratification in genome-wide association studies","volume":"38","author":"AL Price","year":"2006","journal-title":"Nat Genet"},{"key":"pone.0249305.ref009","doi-asserted-by":"crossref","first-page":"394","DOI":"10.1002\/gepi.20393","article-title":"On the adjustment for covariates in genetic association analysis: A novel, simple principle to infer direct causal effects","volume":"33","author":"S Vansteelandt","year":"2009","journal-title":"Genet Epidemiol"},{"issue":"6","key":"pone.0249305.ref010","doi-asserted-by":"crossref","DOI":"10.1515\/sagmb-2017-0054","article-title":"A practical approach to adjusting for population stratification in genome-wide association studies: Principal components and propensity scores (PCAPS)","volume":"17","author":"H Zhao","year":"2018","journal-title":"Stat Appl Genet Mol Biol"},{"key":"pone.0249305.ref011","first-page":"260","article-title":"Outgroup machine learning approach identifies single nucleotide variants in noncoding DNA associated with autism spectrum disorder","volume":"24","author":"M Varma","year":"2019","journal-title":"Pacific Symp Biocomput"},{"issue":"1","key":"pone.0249305.ref012","doi-asserted-by":"crossref","first-page":"29","DOI":"10.1093\/biostatistics\/kxv027","article-title":"Methods that remove batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses","volume":"17","author":"V Nygaard","year":"2016","journal-title":"Biostatistics"},{"key":"pone.0249305.ref013","doi-asserted-by":"crossref","first-page":"1297","DOI":"10.1101\/gr.107524.110","article-title":"The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data","volume":"20","author":"A McKenna","year":"2010","journal-title":"Genome Res"},{"issue":"3","key":"pone.0249305.ref014","doi-asserted-by":"crossref","first-page":"263","DOI":"10.1016\/j.jalz.2011.03.005","article-title":"The diagnosis of dementia due to Alzheimer\u2019s disease: Recommendations from the National Institute on Aging-Alzheimer\u2019s Association workgroups on diagnostic guidelines for Alzheimer\u2019s disease","volume":"7","author":"GM McKhann","year":"2011","journal-title":"Alzheimer\u2019s Dement"},{"key":"pone.0249305.ref015","doi-asserted-by":"crossref","first-page":"921","DOI":"10.1126\/science.8346443","article-title":"Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer\u2019s disease in late onset families","volume":"261","author":"E Corder","year":"1993","journal-title":"Science (80-)"},{"key":"pone.0249305.ref016","article-title":"Identification of missing variants by combining multiple analytic pipelines","volume":"19","author":"Y Ren","year":"2018","journal-title":"BMC Bioinformatics"},{"key":"pone.0249305.ref017","doi-asserted-by":"crossref","first-page":"839","DOI":"10.1016\/j.ajhg.2012.09.004","article-title":"Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data","volume":"91","author":"G Jun","year":"2012","journal-title":"Am J Hum Genet"},{"issue":"22","key":"pone.0249305.ref018","doi-asserted-by":"crossref","first-page":"2867","DOI":"10.1093\/bioinformatics\/btq559","article-title":"Robust relationship inference in genome-wide association studies","volume":"26","author":"A Manichaikul","year":"2010","journal-title":"Bioinformatics"},{"key":"pone.0249305.ref019","article-title":"Whole exome sequencing study identifies novel rare and common Alzheimer\u2019s-Associated variants involved in immune response and transcriptional regulation","author":"JC Bis","year":"2018","journal-title":"Mol Psychiatry"},{"key":"pone.0249305.ref020","doi-asserted-by":"crossref","first-page":"1655","DOI":"10.1101\/gr.094052.109","article-title":"Fast model-based estimation of ancestry in unrelated individuals","volume":"19","author":"DH Alexander","year":"2009","journal-title":"Genome Res"},{"key":"pone.0249305.ref021","doi-asserted-by":"crossref","first-page":"1549","DOI":"10.3233\/JAD-160953","article-title":"Localization and processing of the amyloid-\u03b2 protein precursor in mitochondria-associated membranes","volume":"55","author":"D Del Prete","year":"2017","journal-title":"J Alzheimer\u2019s Dis"},{"key":"pone.0249305.ref022","doi-asserted-by":"crossref","DOI":"10.3390\/molecules22101692","article-title":"\u03b2-amyloid and the pathomechanisms of Alzheimer\u2019s disease: A comprehensive view","volume":"22","author":"B Penke","year":"2017","journal-title":"Molecules"},{"key":"pone.0249305.ref023","doi-asserted-by":"crossref","first-page":"811","DOI":"10.1038\/s41593-018-0154-9","article-title":"A molecular network of the aging human brain provides insights into the pathology and cognitive decline of Alzheimer\u2019s disease","volume":"21","author":"S Mostafavi","year":"2018","journal-title":"Nat Neurosci"},{"issue":"22","key":"pone.0249305.ref024","doi-asserted-by":"crossref","first-page":"6097","DOI":"10.1158\/0008-5472.CAN-17-1018","article-title":"Cic loss promotes gliomagenesis via aberrant neural stem cell proliferation and differentiation","volume":"77","author":"R Yang","year":"2017","journal-title":"Cancer Res"},{"key":"pone.0249305.ref025","article-title":"TTBK2: A tau protein kinase beyond tau phosphorylation","author":"J Liao","year":"2015","journal-title":"Biomed Res Int"},{"issue":"7","key":"pone.0249305.ref026","article-title":"Pathological phosphorylation of tau and TDP-43 by TTBK1 and TTBK2 drives neurodegeneration","volume":"13","author":"LM Taylor","year":"2018","journal-title":"Mol Neurodegener"},{"issue":"3","key":"pone.0249305.ref027","doi-asserted-by":"crossref","first-page":"448","DOI":"10.1016\/j.ajhg.2018.07.019","article-title":"De novo mutations of CCNK cause a syndromic neurodevelopmental disorder with distinctive facial dysmorphism","volume":"103","author":"Y Fan","year":"2018","journal-title":"Am J Hum Genet"},{"key":"pone.0249305.ref028","unstructured":"Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv. 2013;1303.3997(00):http:\/\/arxiv.org\/abs\/1303.3997."},{"key":"pone.0249305.ref029","doi-asserted-by":"crossref","first-page":"186","DOI":"10.1016\/j.ygeno.2017.03.001","article-title":"Evaluation and assessment of read-mapping by multiple next-generation sequencing aligners based on genome-wide characteristics","volume":"109","author":"S Thankaswamy-Kosalai","year":"2017","journal-title":"Genomics"},{"key":"pone.0249305.ref030","doi-asserted-by":"crossref","first-page":"506","DOI":"10.1111\/nan.12452","article-title":"Whole-exome sequencing of the BDR cohort: evidence to support the role of the PILRA gene in Alzheimer\u2019s disease","volume":"44","author":"T Patel","year":"2018","journal-title":"Neuropathol Appl Neurobiol"},{"key":"pone.0249305.ref031","article-title":"From fastQ data to high-confidence variant calls: The genome analysis toolkit best practices pipeline","author":"GA Van der Auwera","year":"2013","journal-title":"Current Protocols in Bioinformatics"},{"issue":"2","key":"pone.0249305.ref032","doi-asserted-by":"crossref","first-page":"80","DOI":"10.4161\/fly.19695","article-title":"A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3","volume":"6","author":"P Cingolani","year":"2012","journal-title":"Fly (Austin)"},{"issue":"16","key":"pone.0249305.ref033","doi-asserted-by":"crossref","DOI":"10.1093\/nar\/gkq603","article-title":"ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data","volume":"38","author":"K Wang","year":"2010","journal-title":"Nucleic Acids Res"},{"key":"pone.0249305.ref034","doi-asserted-by":"crossref","DOI":"10.1186\/s13742-015-0047-8","article-title":"Second-generation PLINK: rising to the challenge of larger and richer datasets","volume":"4","author":"C Chang","year":"2015","journal-title":"Gigascience"}],"container-title":["PLOS ONE"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0249305","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,4,16]],"date-time":"2021-04-16T17:41:51Z","timestamp":1618594911000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0249305"}},"subtitle":[],"editor":[{"given":"Alvaro","family":"Galli","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2021,4,16]]},"references-count":34,"journal-issue":{"issue":"4","published-online":{"date-parts":[[2021,4,16]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pone.0249305","relation":{},"ISSN":["1932-6203"],"issn-type":[{"value":"1932-6203","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,4,16]]}}}