MetaTranscriptomics amount of viable reads
1
0
Entering edit mode
3 months ago
zjhennin • 0

Hey all, I have a more scientific question in nature. Essentially, my group has done meta-transcriptomics for the first time in our lab and we encountered some issues with extraction. Under the guidance of the sequencing facility we then used a whole lane to sequence everything (no rRNA depletion) and recieved ~ 40million reads per sample. Assuming we have 90% rRNA we'll end up with ~4million reads. Is this enough for differential expression analysis? Has this type of analysis been done on communities with this low level of sequencing in the past? I'm trying to find technical validity in moving forward, as well as trying to learn what limitations analysis wise we should be concerned about. Any help, references, discussion, etc. would be greatly appreciated!

metagenomics rna-seq transcriptomics • 488 views
ADD COMMENT
0
Entering edit mode

If you are working with a synthetic community composed by few strains, 4 million reads are not great not terrible.

With natural community however, you will probably get some (differential expressed?) transcripts representing highly transcribed genes from the most abundant and transcriptionally active microbes. How much informative this could be? it totally depends on what are you looking for.

ADD REPLY
0
Entering edit mode

This is a natural community derived from marine sampling. One of our main ideas was to look at differentially expressed transcripts for specific metabolisms. So in essence, we may be able to determine some active pathways in the most abundant microbes with these data?

ADD REPLY
1
Entering edit mode

So in essence, we may be able to determine some active pathways in the most abundant microbes with these data?

If the community isn't to complex, yes... but keep your expectations low

ADD REPLY
1
Entering edit mode
3 months ago

You should try Krakenuniq or Kraken 2 in your RNA-Seq samples...

ADD COMMENT
0
Entering edit mode

This will determine the taxonomy of who is active, but not the genes present correct? Good idea, but as far as my knowledge goes these are mainly used with 16s databases (which we should have plenty of rRNA present), but can gene (to find mRNA) databases be used?

ADD REPLY

Login before adding your answer.

Traffic: 1968 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6