Comments (3)
Hey Vlad,
The gnomAD threshold gives you the opportunity to essentially rule out (from the report) any variant that is relatively common in the population. Most pathogenic variants in cancer susceptibility genes (here I am thinking of the ones with high penetrance, i.e. with an P/LP status in ClinVar) have very low allele frequency, or are even absent from gnomAD. This is supported by standard pathogenicity scoring criteria (ACMG), which gives you a score towards "pathogenic" if the allele frequency is very low, whereas it gives you a score towards "benign" if the allele frequency is high (likely benign). Basically, setting this threshold to a low level in the configuration file will tend to "enrich" your report for the variants that are more likely of pathogenic nature. My understanding is also that regulatory concerns regarding the eligibility to report on relatively frequent variants may exist (at least here in Norway).
Note that the situation for variants found through GWAS, the story is different, they have a much higher population frequency, but with a lower penetrance. The gnomAD allele frequency filtering is not applied to any variant that overlaps with known GWAS hit variants (at least in the dev version (:-)), that should probably be communicated more clearly).
Either way, setting it to a high level will keep all variants in the predisposition genes, but the accumulation of benign variants in the report will consequently tend to rise. I guess it boils down to finding a balance here, and at which threshold this can be achieved.
best,
Sigve
from cpsr.
That totally makes sense, really appreciate the detailed explanation! I was trying to hunt several missing MHS variants in hg38 compared to GRCh37 runs (e.g. rs1650697). For some reason the gnomAD frequencies are reported differently, e.g. GRCh37:
GLOBAL_AF_GNOMAD 0
NON_CANCER_AF_GLOBAL 0.8565
And hg38:
GLOBAL_AF_GNOMAD 0.857
NON_CANCER_AF_GLOBAL 0.8565
And because of the missing GLOBAL_AF_GNOMAD, in 37 the variants are kept. Even though the representation of the variant in gnomad_cpsr/gnomad_cpsr.vcf
looks identical. Probably not a big deal since in the worst case it just means a bit more clutter, without losing anything important.
Thanks again!
from cpsr.
Cutoff's are now only applied for unclassified variants (non-ClinVar).
from cpsr.
Related Issues (20)
- Custom panel HOT 11
- ERROR while running HOT 3
- cpsr_validate_input.py: error: argument vcf_validation: invalid int value: '/home/zdjzyx01/biosoft/cpsr/cpsr.toml' HOT 11
- Missed variant due to transcript pick? HOT 9
- cpsr_configuration_default.toml missing from PCGR data HOT 2
- Error: Problem with `mutate()` input `ACMG_PVS1_1`. HOT 4
- Error in loadNamespace(name) : there is no package called 'reactable' HOT 2
- x object 'GWAS_CITATION' not found HOT 3
- ImportError: No module named toml HOT 1
- minor typo in src/cpsr_functions.R HOT 1
- Error: `x` and `y` must share the same src HOT 5
- genes being excluded from CPSR exploratory track HOT 5
- Not able to download the data bundles (404 Error) HOT 2
- ACMG guideline updated to version 3.0. New genes were added to the secondary finding gene list. HOT 2
- make infromation like ckb to cpsr HOT 4
- Encounter the error in cpsr-report-generation() HOT 3
- Installation Instructions Only Describe Conda HOT 2
- Add sif file for running on an HPC HOT 1
- Biomarker variants not in Pathogenic/Likely Pathogenic section of report HOT 1
- adding clones IDs + child clone emerging outside parent clone HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from cpsr.