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

Metatranscriptomic Data #48

Open
sh1ro-1 opened this issue Dec 9, 2023 · 4 comments
Open

Metatranscriptomic Data #48

sh1ro-1 opened this issue Dec 9, 2023 · 4 comments

Comments

@sh1ro-1
Copy link

sh1ro-1 commented Dec 9, 2023

I have some metatranscriptomic data and I want to run it through ARGs-OAP, is this reasonable?
If not, any other suggestions would be greatly appreciated.
Thank you so much for taking the time to read this, looking forward to your reply!

@xinehc
Copy link
Owner

xinehc commented Dec 11, 2023

hi,

args-oap can only handle short reads so if your metatranscriptomic data are generated by e.g. Illumina data then ok, but the results may not be directly comparable with metagenomic data.

@sh1ro-1
Copy link
Author

sh1ro-1 commented Jan 7, 2024

My metatranscriptomic data are generated by Illumina data, so args-oap is available.
I have analyzed the RPKM results that I obtained, but I am very confused about why the results may not be directly comparable with metagenomic data, probably because I do not understand the principle well enough. I do not know how to deal with this issue. If possible, I hope you can clear up my confusion, thank you very much.

@xinehc
Copy link
Owner

xinehc commented Jan 9, 2024

RPKM requires normalising against read counts. For metagenomic, reads are randomly distributed along the genomes; for metatranscriptomic, only genes that express will be sequenced. So the total number of reads (normalising constant) has different meaning in metagenomic/metatranscriptomic, which will lead to a difference.

@jiaojiaoguan
Copy link

RPKM requires normalising against read counts. For metagenomic, reads are randomly distributed along the genomes; for metatranscriptomic, only genes that express will be sequenced. So the total number of reads (normalising constant) has different meaning in metagenomic/metatranscriptomic, which will lead to a difference.

Hello,
I have a question about the RPKM. I think the RPKM for metagenomic and metatranscriptomic are same. For example, in metagenomic,

ARG1 mapped reads number =10 length =1 ARG2 mapped reads number =20 length =2
Then, the RPKM for ARG1 is 101e9/130.

I think the number of the mapped reads for each arg is the same for metagenomic (30) and metatranscriptomic (30). So, the total number of the mapped reads is the same (the normalizing constant is the same).

what do you think of it? Thanks very much for your answer.

Best,
jiaojiao

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

3 participants