Question: Bio-statistics for metatranscriptomic
gravatar for amelie.rouger
3.3 years ago by
amelie.rouger0 wrote:

Hi everybody, I am trying to identified differentially expressed (DE) genes of several bacterial metatranscriptomes.

To sum up the experimental plan. We have 2 bacterial communities (named E and U) inoculated under 3 different conditions A, B and C (reference condition). I collected both DNA and RNA from time 7 and 9 and also the experiment was performed in duplicate. We sequenced RNA (each 24 samples) and metagenome (from pooled of DNA of several samples).

I assembled and annotated metagenome (141 921 contigs) and use it as reference to map RNA reads with bowtie2.

I performed DESeq2 analysis on those result and I obtained a large part of DE gene. So I used edgeR to change the normalization way (I had tested TMM, RLE and upperquartile methods) and I obtained around 10 000 DE genes for each comparison of condition (AvsC BvsC and AvsB).

The problem is that I don't know what is the next step. With a small list of DE gene I would blast the sequence on NCBI and try to identified from which bacterial species, this gene came from and try to identified metabolism pathwas by hand. But with a too large number of DE gene, I don't know how to do.

Do you think that I have to find another way of normalization to reduce the number of DE gene? Or do you have any idea of software I could use to group those gene by "category" or "class" to go throught the metabolic pathway?

Thank's for your help.


rna-seq stats pathway • 1.1k views
ADD COMMENTlink modified 3.2 years ago by Biostar ♦♦ 20 • written 3.3 years ago by amelie.rouger0

Did you map the RNA to the DNA? Maybe there are whole genomes which are DE, which is interesting by itself but in a different context, this might narrow down the list of genes.

ADD REPLYlink written 3.3 years ago by Asaf7.6k

How did you annotate the metagenome assembly?

ADD REPLYlink written 3.3 years ago by h.mon29k

Thank's for answering. Yes I used metagenome to map the RNA reads. I assembled metagenome with IDBA and annotated it with PROKKA.

When I used available genome on ncbi to map the RNA reads I aligned only 30% of reads, that why we used metagenomes to map RNA reads.

ADD REPLYlink written 3.3 years ago by amelie.rouger0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2102 users visited in the last hour