feat: StringTie novel transcript discovery + hybrid GTF#182
Draft
pinin4fjords wants to merge 4 commits into
Draft
feat: StringTie novel transcript discovery + hybrid GTF#182pinin4fjords wants to merge 4 commits into
pinin4fjords wants to merge 4 commits into
Conversation
Add reference-guided novel-transcript assembly from RNA-seq BAMs (or Ribo-seq fallback with tightened defaults), or accept a user-supplied novel GTF via --novel_gtf. The two sources feed a common chain: gffcompare class-code filter + optional strand-aware rRNA/repeat blacklist + concat with the canonical backbone, via the upstream gtf_hybridmerge_gffcompare subworkflow. The hybrid GTF is published at <outdir>/stringtie/hybrid_reference.gtf and exposed on the `hybrid_gtf` workflow channel. With default --skip_stringtie true and no --novel_gtf, the channel equals ch_canonical_gtf so downstream wiring is uniform. Wiring of the hybrid GTF into ORF callers lives in follow-on issue #165 (--extended_orf_analysis). New params: --skip_stringtie (default true) --novel_gtf --stringtie_class_codes (default 'u', intergenic only) --rrna_blacklist --extra_stringtie_args --extra_stringtie_merge_args --stringtie_ribo_fallback_args (default '-m 100 -c 5 -j 3 -f 0.05 -g 100') Modules and subworkflows added: stringtie/stringtie, stringtie/merge, gffcompare, bedtools/intersect, bam_stringtie_merge, gtf_hybridmerge_gffcompare
Member
|
Warning Newer version of the nf-core template is available. Your pipeline is using an old version of the nf-core template: 3.5.1. For more documentation on how to update your pipeline, please see the nf-core documentation and Synchronisation documentation. |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Summary
Adds reference-guided novel-transcript assembly via StringTie (preferring RNA-seq BAMs, falling back to Ribo-seq with tightened defaults) or a user-supplied annotation via
--novel_gtf. Either source flows through a common filtering chain (gffcompare class-code filter, optional strand-aware rRNA/repeat blacklist, concat with the canonical backbone) and produces a hybrid annotation at<outdir>/stringtie/hybrid_reference.gtf.This is the foundation for #165 / #171 / extended ORF discovery. With the default
--skip_stringtie trueand no--novel_gtf, the newhybrid_gtfworkflow channel equalsch_canonical_gtffrom #161, so downstream wiring is uniform whether novel discovery is on or off.Changes
NOVEL_TRANSCRIPT_DISCOVERYchainingbam_stringtie_merge(per-sample StringTie + merge) and the upstreamgtf_hybridmerge_gffcomparesubworkflow (gffcompare + class-code filter + optional blacklist intersect + concat with canonical).stringtie/stringtie,stringtie/merge,gffcompare,bedtools/intersect.bam_stringtie_merge,gtf_hybridmerge_gffcompare(PR feat(subworkflows): gtf_hybridmerge_gffcompare modules#11729, merged).--skip_stringtie(defaulttrue),--novel_gtf,--stringtie_class_codes(default'u', intergenic only),--rrna_blacklist,--extra_stringtie_args,--extra_stringtie_merge_args,--stringtie_ribo_fallback_args(default'-m 100 -c 5 -j 3 -f 0.05 -g 100'to tighten the Ribo-seq fallback).hybrid_gtffrom the riboseq workflow.Caveats
--stringtie_class_codes = ukeeps only unambiguously intergenic transcripts. Stranded users can opt into'u,x'to recover antisense; unsafe for non-stranded protocols.--rrna_blacklistis silently skipped when not supplied.Stacked PR notes
Sixth in the stack splitting #174. Targets #181 (
feat/162-ribotish-qc-canonical-doc).Closes #164
🤖 Generated with Claude Code