Short tandem repeat number estimation from paired-end reads for multiple individuals by considering coalescent tree

BMC Genomics

Research

Naoki Nariai¹  Takahiro Mimori²  Kaname Kojima²  Yosuke Kawai²  Masao Nagasaki²  Takanori Hasegawa² 

[1] Institute for Genomic Medicine, University of California, San Diego, 9500 Gilman Drive #0761, 92093-0761, San Diego, USA;Tohoku Medical Megabank Organization, Tohoku University, 2-1, Seiryo-machi, Aoba-ku, 980-8573, Sendai, Japan;

关键词: High-throughput sequencing;    Short tandem repeat;    Coalescent theory;   

DOI  :  10.1186/s12864-016-2821-0

来源: Springer

【 摘 要 】

BackgroundTwo types of approaches are mainly considered for the repeat number estimation in short tandem repeat (STR) regions from high-throughput sequencing data: approaches directly counting repeat patterns included in sequence reads spanning the region and approaches based on detecting the difference between the insert size inferred from aligned paired-end reads and the actual insert size. Although the accuracy of repeat numbers estimated with the former approaches is high, the size of target STR regions is limited to the length of sequence reads. On the other hand, the latter approaches can handle STR regions longer than the length of sequence reads. However, repeat numbers estimated with the latter approaches is less accurate than those with the former approaches.ResultsWe proposed a new statistical model named coalescentSTR that estimates repeat numbers from paired-end read distances for multiple individuals simultaneously by connecting the read generative model for each individual with their genealogy. In the model, the genealogy is represented by handling coalescent trees as hidden variables, and the summation of the hidden variables is taken on coalescent trees sampled based on phased genotypes located around a target STR region with Markov chain Monte Carlo. In the sampled coalescent trees, repeat number information from insert size data is propagated, and more accurate estimation of repeat numbers is expected for STR regions longer than the length of sequence reads.For finding the repeat numbers maximizing the likelihood of the model on the estimation of repeat numbers, we proposed a state-of-the-art belief propagation algorithm on sampled coalescent trees.ConclusionsWe verified the effectiveness of the proposed approach from the comparison with existing methods by using simulation datasets and real whole genome and whole exome data for HapMap individuals analyzed in the 1000 Genomes Project.

【 授权许可】

CC BY   
© The Author(s) 2016

【 图 表 】

12951_2016_177_Article_IEq1.gif

12951_2015_155_Article_IEq65.gif

Fig. 4

Fig. 5

Fig. 3

Fig. 1

Fig. 2

Fig. 1

12936_2015_836_Article_IEq14.gif

12951_2017_255_Article_IEq41.gif

12936_2015_836_Article_IEq13.gif

Fig. 2

12936_2017_2014_Article_IEq29.gif

Fig. 2

Fig. 4

12936_2017_2045_Article_IEq18.gif

Fig. 2

Fig. 2

Fig. 1

Fig. 1

Fig. 3

Fig. 2

12936_2016_1182_Article_IEq39.gif

12951_2017_255_Article_IEq40.gif

Fig. 1

Fig. 4

Fig. 2

Fig. 6

Fig. 2

Fig. 1

Fig. 1

Fig. 4

Fig. 5

Fig. 6

Fig. 2

Fig. 5

Fig. 1

Fig. 1

Fig. 2

12951_2015_155_Article_IEq62.gif

12951_2017_292_Article_IEq1.gif

Fig. 6

Fig. 3

Fig. 5

Fig. 6

Fig. 4

Fig. 4

12937_2016_133_Article_IEq1.gif

Fig. 3

Fig. 3

Fig. 1

Fig. 2

Fig. 10