1. Profiling of Short-Tandem-Repeat Disease Alleles in 12,632 Human Whole Genomes
- Author
-
Haibao Tang, Ewen F. Kirkness, Victor Lavrenko, Claire Hou, William H. Biggs, Smriti Ramakrishnan, Amalio Telenti, C. Thomas Caskey, Martin M. Fabani, Boyko Kakaradov, David Heckerman, Christoph Lippert, Barry W. Hicks, Ernesto Guzman, Franz Josef Och, and J. Craig Venter
- Subjects
Adult ,Male ,0301 basic medicine ,Adolescent ,STR multiplex system ,Computational biology ,Biology ,Genome ,Article ,microsatellites ,DNA sequencing ,03 medical and health sciences ,Genetics ,Humans ,Child ,Genotyping ,Alleles ,Genetics (clinical) ,Polymorphism, Genetic ,Genome, Human ,High-Throughput Nucleotide Sequencing ,population genetics ,Sequence Analysis, DNA ,Middle Aged ,short tandem repeats ,genome sequencing ,Genetics, Population ,030104 developmental biology ,genotyping ,STR analysis ,Microsatellite ,Female ,genetic disorder ,Human genome ,Nanopore sequencing ,trinucleotide repeat expansion ,Software ,Microsatellite Repeats - Abstract
Short tandem repeats (STRs) are hyper-mutable sequences in the human genome. They are often used in forensics and population genetics and are also the underlying cause of many genetic diseases. There are challenges associated with accurately determining the length polymorphism of STR loci in the genome by next-generation sequencing (NGS). In particular, accurate detection of pathological STR expansion is limited by the sequence read length during whole-genome analysis. We developed TREDPARSE, a software package that incorporates various cues from read alignment and paired-end distance distribution, as well as a sequence stutter model, in a probabilistic framework to infer repeat sizes for genetic loci, and we used this software to infer repeat sizes for 30 known disease loci. Using simulated data, we show that TREDPARSE outperforms other available software. We sampled the full genome sequences of 12,632 individuals to an average read depth of approximately 30× to 40× with Illumina HiSeq X. We identified 138 individuals with risk alleles at 15 STR disease loci. We validated a representative subset of the samples (n = 19) by Sanger and by Oxford Nanopore sequencing. Additionally, we validated the STR calls against known allele sizes in a set of GeT-RM reference cell-line materials (n = 6). Several STR loci that are entirely guanine or cytosines (G or C) have insufficient read evidence for inference and therefore could not be assayed precisely by TREDPARSE. TREDPARSE extends the limit of STR size detection beyond the physical sequence read length. This extension is critical because many of the disease risk cutoffs are close to or beyond the short sequence read length of 100 to 150 bases.
- Published
- 2017
- Full Text
- View/download PDF