sigven / pcgr

Personal Cancer Genome Reporter (PCGR)
https://sigven.github.io/pcgr
MIT License
251 stars 47 forks source link

Regarding CNV analysis on germline VCF #219

Closed AidanGCr closed 9 months ago

AidanGCr commented 10 months ago

Hello,

I understand that CPSR does not include an analysis framework as of now for CNV / SV analysis. However, CNV analysis is done within PCGR. I was wondering, would running a germline sample in tumor-only mode allow for effective prioritization of CNV calls? Of course, this would be with relevant settings disabled such as the following (as mentioned in the readme):

--exclude_dbsnp_nonsomatic
--exclude_likely_het_germline
--exclude_likely_hom_germline
--exclude_nonexonic

Thanks, Aidan

sigven commented 10 months ago

Dear Aidan,

Thanks for your question. Note that germline variants in the form of copy number is not a particularly frequent phenomenon when it comes to cancer predisposition, this is why this is currently omitted in the CPSR input. Also not that the above filtering strategy only works for SNVs/InDels, and is thus it will not really be applicable for filtering a germline sample. PCGR only accepts somatic CNA calls - no filtering is performed on those.

kind regards, Sigve

AidanGCr commented 10 months ago

Hi Sigve,

I see, thank you for clarifying about the filtering. Considering, that the SV/CNV filtering is done beforehand on my end, with data then being transformed into the tsv format, PCGR would then classify all CNV calls it could, with the assumption of them all being somatic, correct?

Thanks, Aidan

sigven commented 10 months ago

Hi Aidan,

Yes, exactly. This is correctly interpreted. I am also now working on a redesign for the CNV input and reporting (considering allele-specific copy numbers (major + minor allele)), rather than the log2 ratios as it is currently handled. I am also working on providing better visualizations for these data, e.g. so that the user can view the copy number segments better: https://caravagnalab.github.io/CNAqc/articles/a2_Plots.html#segments-plots

So be prepared for some changes :-), will keep you posted.

best, Sigve

AidanGCr commented 9 months ago

Great! I look forward to seeing the updates!

Thanks, Aidan