STR-realigner: a realignment method for short tandem repeat regions

BMC Genomics

Research Article

Takahiro Mimori¹  Kaname Kojima¹  Yosuke Kawai¹  Masao Nagasaki¹  Kazuharu Misawa¹ 

[1] Tohoku Medical Megabank Organization, Tohoku University, 2-1, Seiryo-machi, Aoba-ku, 980-8573, Sendai, Japan;

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

DOI  :  10.1186/s12864-016-3294-x

 received in 2016-05-13, accepted in 2016-11-15,  发布年份 2016

来源: Springer

【 摘 要 】

BackgroundIn the estimation of repeat numbers in a short tandem repeat (STR) region from high-throughput sequencing data, two types of strategies are mainly taken: a strategy based on counting repeat patterns included in sequence reads spanning the region and a strategy based on estimating the difference between the actual insert size and the insert size inferred from paired-end reads. The quality of sequence alignment is crucial, especially in the former approaches although usual alignment methods have difficulty in STR regions due to insertions and deletions caused by the variations of repeat numbers.ResultsWe proposed a new dynamic programming based realignment method named STR-realigner that considers repeat patterns in STR regions as prior knowledge. By allowing the size change of repeat patterns with low penalty in STR regions, accurate realignment is expected. For the performance evaluation, publicly available STR variant calling tools were applied to three types of aligned reads: synthetically generated sequencing reads aligned with BWA-MEM, those realigned with STR-realigner, those realigned with ReviSTER, and those realigned with GATK IndelRealigner. From the comparison of root mean squared errors between estimated and true STR region size, the results for the dataset realigned with STR-realigner are better than those for other cases. For real data analysis, we used a real sequencing dataset from Illumina HiSeq 2000 for a parent-offspring trio. RepeatSeq and lobSTR were applied to the sequence reads for these individuals aligned with BWA-MEM, those realigned with STR-realigner, ReviSTER, and GATK IndelRealigner. STR-realigner shows the best performance in terms of consistency of the size of estimated STR regions in Mendelian inheritance. Root mean squared error values were also calculated from the comparison of these estimated results with STR region sizes obtained from high coverage PacBio sequencing data, and the results from the realigned sequencing data with STR-realigner showed the least (the best) root mean squared error value.ConclusionsThe effectiveness of the proposed realignment method for STR regions was verified from the comparison with an existing method on both simulation datasets and real whole genome sequencing dataset.

【 授权许可】

CC BY   
© The Author(s) 2016

【 图 表 】

12864_2017_4269_Article_IEq3.gif

12864_2017_4269_Article_IEq1.gif

12864_2017_3655_Article_IEq6.gif

12864_2017_4133_Article_IEq36.gif

12864_2015_2198_Article_IEq11.gif

12864_2017_4132_Article_IEq34.gif

12864_2017_4132_Article_IEq33.gif

12864_2017_3777_Article_IEq5.gif

12864_2015_2242_Article_IEq1.gif

12864_2017_3733_Article_IEq43.gif

12864_2016_3425_Article_IEq6.gif