Question: CNV calling using ExomeDepth
gravatar for Nandini
21 months ago by
Nandini840 wrote:


I'm working to establish a pipeline for germline CNV calling. Our target consists of 150 genes with approx 1500 exons. Out data is HiSeq data (96 samples per run)

To establish a good reference set, I'm following this paper to calculate the inter-sample variation in coverage, using the rpkmCV for surveyed exons across reference samples selected by ExomeDepth. I know the basic formula for this is SD/Mean; but just not sure how to implement this in my data.

Does anyone know how to do this ? Or any suggestions on QC for the same ?

Thank you

cnv reference exomedepth • 1.4k views
ADD COMMENTlink modified 21 months ago by andrew.j.skelton735.9k • written 21 months ago by Nandini840
gravatar for andrew.j.skelton73
21 months ago by
andrew.j.skelton735.9k wrote:

ExomeDepth works on coverage, and there's functionality within ExomeDepth to determine an optimal reference / set of references to use. See section 5 here in the vignette. ExomeDepth typically goes with highly correlated samples, hence it's essential that your samples are from the same batch of sequencing. One thing to consider is if you have trios or related individuals in a pedigree, or a mix of probands and healthy samples, then which samples should / should not be included in your reference for a given test.

If you want to use the inter/ intra sample variation, then you'd apply this to the imported normalised counts. You could do something like an rLog or VST transformation before hand, but this will take some experimentation as that paper is quite ambiguous in its methods.

ADD COMMENTlink written 21 months ago by andrew.j.skelton735.9k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 660 users visited in the last hour