Skip to content
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.

Sign up for GitHub

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

Jump to bottom

Added a '--prefer-mane-transcripts' mode that enforces MANE_Select tagged Gencode transcripts where possible #9012

Merged
merged 3 commits into from
Oct 23, 2024
Merged

Added a '--prefer-mane-transcripts' mode that enforces MANE_Select tagged Gencode transcripts where possible #9012

merged 3 commits into from
Oct 23, 2024

Conversation

jamesemery
Copy link
Collaborator

@jamesemery jamesemery mentioned this pull request Oct 18, 2024
Open
jonn-smith
jonn-smith requested changes Oct 22, 2024
View reviewed changes
Copy link
Collaborator

@jonn-smith jonn-smith left a 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
...roadinstitute/hellbender/tools/funcotator/dataSources/gencode/GencodeFuncotationFactory.java Outdated
@@ -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:
Copy link
Collaborator

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)

Copy link
Collaborator Author

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
@jamesemery
Copy link
Collaborator Author

@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.
transcript "ENST00000643187.1" is not actually classified as MANE_Plus_Clinical but it is included in the file for testing purposes.)

jonn-smith
jonn-smith requested changes Oct 22, 2024
View reviewed changes
Copy link
Collaborator

@jonn-smith jonn-smith left a 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
jonn-smith
jonn-smith approved these changes Oct 22, 2024
View reviewed changes
Copy link
Collaborator

@jonn-smith jonn-smith left a 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 jonn-smith merged commit c4860d4 into master Oct 23, 2024
20 checks passed
@jonn-smith jonn-smith deleted the je_ManeTranscriptModeFunctotator branch October 23, 2024 12:49
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Reviewers

@jonn-smith jonn-smith jonn-smith approved these changes

Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging this pull request may close these issues.
2 participants
@jamesemery @jonn-smith