Reconstructing cancer karyotypes from short read data: the half empty and half full glass

BMC Bioinformatics

Methodology Article

Rami Eitan¹  Ron Shamir¹ 

[1] Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv-Yafo, Israel;

关键词: Cancer;    Karyotypes;    Genome rearrangements;    Structural and numerical variations;    Deep sequencing;    Reconstruction;    Graph theory;    Integer linear programming;   

DOI  :  10.1186/s12859-017-1929-9

 received in 2017-06-20, accepted in 2017-11-06,  发布年份 2017

来源: Springer

【 摘 要 】

BackgroundDuring cancer progression genomes undergo point mutations as well as larger segmental changes. The latter include, among others, segmental deletions duplications, translocations and inversions.The result is a highly complex, patient-specific cancer karyotype. Using high-throughput technologies of deep sequencing and microarrays it is possible to interrogate a cancer genome and produce chromosomal copy number profiles and a list of breakpoints (“jumps”) relative to the normal genome. This information is very detailed but local, and does not give the overall picture of the cancer genome. One of the basic challenges in cancer genome research is to use such information to infer the cancer karyotype.We present here an algorithmic approach, based on graph theory and integer linear programming, that receives segmental copy number and breakpoint data as input and produces a cancer karyotype that is most concordant with them. We used simulations to evaluate the utility of our approach, and applied it to real data.ResultsBy using a simulation model, we were able to estimate the correctness and robustness of the algorithm in a spectrum of scenarios. Under our base scenario, designed according to observations in real data, the algorithm correctly inferred 69% of the karyotypes. However, when using less stringent correctness metrics that account for incomplete and noisy data, 87% of the reconstructed karyotypes were correct. Furthermore, in scenarios where the data were very clean and complete, accuracy rose to 90%–100%. Some examples of analysis of real data, and the reconstructed karyotypes suggested by our algorithm, are also presented.ConclusionWhile reconstruction of complete, perfect karyotype based on short read data is very hard, a large fraction of the reconstruction will still be correct and can provide useful information.

【 授权许可】

CC BY   
© The Author(s). 2017

【 图 表 】

12864_2017_3781_Article_IEq3.gif

12864_2017_3781_Article_IEq2.gif

12864_2017_3990_Article_IEq12.gif

12864_2016_3425_Article_IEq8.gif

12864_2016_3425_Article_IEq7.gif

12864_2017_3733_Article_IEq41.gif

12864_2016_3256_Article_IEq5.gif

12864_2017_3670_Article_IEq14.gif

12864_2017_4071_Article_IEq5.gif

12864_2017_4071_Article_IEq4.gif

12864_2017_3920_Article_IEq1.gif

12864_2016_2580_Article_IEq1.gif

12894_2016_184_Article_IEq3.gif

12864_2017_4071_Article_IEq2.gif

12888_2016_877_Article_IEq22.gif

12864_2017_4132_Article_IEq18.gif

12864_2017_4132_Article_IEq16.gif

12888_2016_877_Article_IEq20.gif

12864_2015_2296_Article_IEq115.gif

12711_2017_365_Article_IEq33.gif

12864_2016_3353_Article_IEq26.gif