-
Notifications
You must be signed in to change notification settings - Fork 597
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Added a '--prefer-mane-transcripts' mode that enforces MANE_Select tagged Gencode transcripts where possible #9012
Conversation
…gged Gencode transcripts where possible
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks pretty good. Some requests for more tests.
Also, can you describe what you changed and / or added to the git LFS / large files?
...in/java/org/broadinstitute/hellbender/tools/funcotator/BaseFuncotatorArgumentCollection.java
Outdated
Show resolved
Hide resolved
...roadinstitute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactory.java
Outdated
Show resolved
Hide resolved
src/main/java/org/broadinstitute/hellbender/tools/funcotator/dataSources/DataSourceUtils.java
Outdated
Show resolved
Hide resolved
...roadinstitute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactory.java
Outdated
Show resolved
Hide resolved
@@ -863,6 +892,11 @@ private List<GencodeFuncotation> createFuncotationsHelper(final VariantContext v | |||
transcriptList = gtfFeature.getTranscripts(); | |||
} | |||
|
|||
// Filter out the non-MANE_Select/Mane_Plus_Clinical transcripts if we're only using MANE transcripts: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Can you rename the input variable in the following method to reflect that the transcripts no longer need to be only basic
?
private List<GencodeFuncotation> createFuncotationsHelper(final VariantContext variant, final Allele altAllele, final ReferenceContext reference, final List<GencodeGtfTranscriptFeature> basicTranscripts)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
that is not accurate though, the transcripts still have to be basic. This simply applies a layer on-top of the basic transcript filtering that enforces that they are MANE select as well (which my cursory survey of gtf files seemed to indicated travels together anyway)
src/test/java/org/broadinstitute/hellbender/tools/funcotator/FuncotatorIntegrationTest.java
Show resolved
Hide resolved
...roadinstitute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactory.java
Outdated
Show resolved
Hide resolved
src/test/java/org/broadinstitute/hellbender/tools/funcotator/FuncotatorIntegrationTest.java
Show resolved
Hide resolved
@jonn-smith i added a readme in the LFS files. Looks like github won't display it though which is slightly annoying. Here is what it says (WARNING.txt: NOTE: this gencode.v43.PIK3CA.gtf file is slightly edited for the sake of testing MANE_SELECT and MANE_PLUS_CLINICAL. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
A few more quick changes requested.
Also, can you rename any of the resource files checked into the large
files repo that you modified to have descriptive names? Something like "MANE_PLUS_CLINICAL_test_data.gtf" would be fine.
...itute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactoryUnitTest.java
Outdated
Show resolved
Hide resolved
...itute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactoryUnitTest.java
Outdated
Show resolved
Hide resolved
src/test/java/org/broadinstitute/hellbender/tools/funcotator/FuncotatorIntegrationTest.java
Outdated
Show resolved
Hide resolved
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good.
👍
@jonn-smith