zeeev
commented
4 years ago Hi @djakubosky,
Thanks for reaching out, these are really good questions, especially for large cohorts. I've made some bullet points that might help guide you. If you find this helpful or come up with your own tricks I'd like to add them to a wiki.
- If you have related individuals batch them together. WhamG becomes more sensitive when you run a trio or quad together. However, there are performance considerations, mostly ram. WhamG stores the discordance read-pair mappings in memory.
- WhamG accepts and exclusion list, you'll want to use this. If I remember correctly it's just a list of sequences (not coordinates). ChromosomeY, alt haplotypes, and bait sequences are smart to ignore, for WhamG. It will decrease your runtimes.
- After you generate the calls, take a look at the distributions of the reported info fields. For example I like to look at the support, or depth of the calls. I designed WhamG to be sensitive so you'll have false positives. For example, segmental duplication regions are highly suspicious.
- Merging helps. If you can run at least 2-3 different callers you gain specificity, and sensitivity. Lumpy/WhamG/Manta seem to be a good combo, although WhamG and Lumpy tend to have a fair amount of overlap (see 1kg phase3 paper). Merging different callers can be a gain pain. There are several tools out there. I wrote one, but I don't know if it's best in class, that was a while ago.
- For the 1KG project, I found, that genotyping the calls helped a lot; I used SVTYPER. You can do a little mendelian or Hardy–Weinberg filtering. I will caution SVTYPER does have a high FN rate (don't have anything formal to show you).
djakubosky
Hi, I was curious about the suggested workflow for WHAMG on ~1K genomes. My assumption is that this is what you would want to do:
1) Run WHAMG individually on all genomes -> many vcfs 2) Filtering? - (not sure if you'd suggest filtering at this stage of the process on each individual VCF) 3) Run mergeSVcallers on these VCFs to create a set of positions 4) Genotype the putative SVs at these positions with something- (eg SVTyper) 5) Merge genotyped variants into one VCF
Does this sound reasonable? If this is the proposed approach might be helpful to add a little more detail in the wiki!
Thanks for a nice tool!