A multilocus approach for accurate variant calling in low-copy repeats using whole-genome sequencing.

Autor:	Prodanov T; Bioinformatics and Systems Biology Graduate Program, University of California San Diego, La Jolla, CA 92093, United States.; Institute for Medical Biometry and Bioinformatics, Medical Faculty, Heinrich Heine University, Düsseldorf 40225, Germany.; Center for Digital Medicine, Heinrich Heine University, Düsseldorf 40225, Germany., Bansal V; School of Medicine, University of California San Diego, La Jolla, CA 92093, United States.
Jazyk:	angličtina
Zdroj:	Bioinformatics (Oxford, England) [Bioinformatics] 2023 Jun 30; Vol. 39 (39 Suppl 1), pp. i279-i287.
DOI:	10.1093/bioinformatics/btad268
Abstrakt:	Motivation: Low-copy repeats (LCRs) or segmental duplications are long segments of duplicated DNA that cover > 5% of the human genome. Existing tools for variant calling using short reads exhibit low accuracy in LCRs due to ambiguity in read mapping and extensive copy number variation. Variants in more than 150 genes overlapping LCRs are associated with risk for human diseases. Methods: We describe a short-read variant calling method, ParascopyVC, that performs variant calling jointly across all repeat copies and utilizes reads independent of mapping quality in LCRs. To identify candidate variants, ParascopyVC aggregates reads mapped to different repeat copies and performs polyploid variant calling. Subsequently, paralogous sequence variants that can differentiate repeat copies are identified using population data and used for estimating the genotype of variants for each repeat copy. Results: On simulated whole-genome sequence data, ParascopyVC achieved higher precision (0.997) and recall (0.807) than three state-of-the-art variant callers (best precision = 0.956 for DeepVariant and best recall = 0.738 for GATK) in 167 LCR regions. Benchmarking of ParascopyVC using the genome-in-a-bottle high-confidence variant calls for HG002 genome showed that it achieved a very high precision of 0.991 and a high recall of 0.909 across LCR regions, significantly better than FreeBayes (precision = 0.954 and recall = 0.822), GATK (precision = 0.888 and recall = 0.873) and DeepVariant (precision = 0.983 and recall = 0.861). ParascopyVC demonstrated a consistently higher accuracy (mean F1 = 0.947) than other callers (best F1 = 0.908) across seven human genomes. Availability and Implementation: ParascopyVC is implemented in Python and is freely available at https://github.com/tprodanov/ParascopyVC. (© The Author(s) 2023. Published by Oxford University Press.)
Databáze:	MEDLINE
Externí odkaz:	Zobrazit plný text záznamu