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.

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

Low increase in discovered isoforms with RNA004 data #443

Open
dietvin opened this issue Aug 23, 2024 · 1 comment
Open

Low increase in discovered isoforms with RNA004 data #443

dietvin opened this issue Aug 23, 2024 · 1 comment

Comments

@dietvin
Copy link

dietvin commented Aug 23, 2024

Hi,
we used Bambu for transcript discovery on RNA002 and RNA004 data of the same sample. As expected the new chemistry produced a lot more reads (roughly 4x depth), but Bambu only identified roughly 10 novel isoforms that were not identified with the RNA002 data. We were expecting larger differences and wonder whether this is due to wrong usage of Bambu or if it's just due to the underlying data?

I used the default parameters to run Bambu:

se = bambu(reads = bam.paths, annotations = annotations, genome = fa.path, quant = TRUE, ncore = 16)

Any help would be appreciated. Thank you.
Vincent

@andredsim
Copy link
Collaborator

Hi Vincent,

Thanks for your report. This is actually what we would expect, as we designed Bambu to be robust to changes in read depth. Higher read depth can lead to the identification of new transcripts, it disproportionately also increases the amount of false positives. If you are confident in the quality of your sample and would like to discover more transcripts, I would recommend increasing the NDR parameter which will increase the number of discovered isoforms. Check the output of your last run for the NDR threshold it applied (it would be around 0.1). A first attempt would be increasing it by 0.1 to 0.2.

I hope this helps,
Andre Sim

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants