Question: When to compare independent trials during differential gene expression analysis
0
gravatar for GB
7 days ago by
GB0
GB0 wrote:

Hi, I’m trying to analyze the differential gene expression between wild type and a particular mutation. In total, I did 3 independent trials using the TruSeq stranded kit with ribosomal depletion.

After obtaining the raw reads per trial, I then performed quality control using fastp. I noticed that the polyA tail might have been sequenced because A/T is overrepresented in the fastp report. I assembled both polyA-trimmed and non-polyA-trimmed sets using Trinity for de-novo transcriptome assembly, after which I assessed the quality using read representation, transrate and metric statistics in hopes of identifying which assembly is better. I am now trying to quantify the transcripts per set-up to be able to do differential expression analysis. Of note, I am doing all these steps independently of each trial, that is, I’m doing all the steps in one trial before repeating it in another.

My questions are:

  1. At which step do I compare the 3 independent trials? For example, do I create one Trinity assembly by using the reads from all three trials? or do I do differential gene expression analysis separately for all the trials and then perform some sort of statistical analysis/ PCA at the end?

  2. Is there anything I can do with a low transrate assembly score (around 0.00012) despite good bowtie2 mapping when I mapped the reads back to the assembly?

Thank you very much in advance!

ADD COMMENTlink modified 2 days ago • written 7 days ago by GB0

What organism? Are you sure you need to assemble transcripts?

ADD REPLYlink written 2 days ago by Asaf8.0k

This is from a mouse cell line but I'm looking for non-coding RNAs (a lot of which are not yet annotated) as well as novel transcripts.

ADD REPLYlink modified 2 days ago • written 2 days ago by GB0
Please log in to add an answer.

Help
Access

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