This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CNV analysis for targeted genomic region

I am currently trying to analyse the presence of CNVs in a targeted sequencing dataset generated from Ion Torrent platform. The data description is as follows.

  • Targeted region is a segment of the genome that spans multiple genes (not limited to exons).
  • Single end reads of variable length.
  • Reads are mapped to the hg19 assembly.
  • The sample are not paired. We have a set of test and control samples.

So far, I have tried analysis using various tools such as CNVkit, CNVSeq, freec, codex. None of them seem to work well with the kind of data that I have in hand.

Can anyone suggest an analysis strategy that I could try with such data?

Thanks.

cnv
• 4,343 views • link
updated by Ram Ram • written by Sandeep Sandeep

Also, I would like to know if I am better off using the whole genome approach of analysis by mapping the reads to the segment of the genome in question?

• 1 views • link
written by Sandeep Sandeep

You could try to calculate the coverage per gene/exon/target, median normalize over all/within sets, build a reference sample supposed to be copy number of two, calculate the ratio between target and reference sample and perform a cluster analysis for CNV detection per gene/exon/target. This is similar to the tools freely available but I obtain more reliable results using a individualized pipeline for targted NGS results.

• 1 views • link
written by Jimbou Jimbou

Thanks for the suggestion. Can you suggest tools that we could use to perform the above steps?

• 1 views • link
written by Sandeep Sandeep

2 answers

Are your using targeted amplicon capture? If so, there is a program designed specifically for this type of data called OncoCNV: https://oncocnv.curie.fr/

• 1 views • link
written by Eric T. Eric T.

Yes, I have used that already and seems like that is working fine for our kind of data.

• 1 views • link
written by Sandeep Sandeep

Hi,

I just met the same question as you. I got targeted sequencing data (~8M, a continous segment in chr22) for 20 cases and 2 controls, and wanted to detect CNV. I wanna try the oncocnv, but don't have a idea about how and where to get the bed file for amplicon coordinates. How do you prepare the amplicon file?

I also plan to try different softs for detecting CNV, however, no much softs were found to be designed for target sequencing data. What else appropriate software have you tried?

Thanks a lot.

• 0 views • link
written by cc • cc

Log in to answer this question.

More posts like this