GC bias between samples tolerance in RNAseq?
Does anyone have any info on the level of GC content disparity that is tolerable between samples? For example, if I have one group that has an average GC content of 41% (+/- 3) that I want to compare to one with an average GC content of 44% (+/- 2)- does this need to be accounted for in my analysis or is this considered equal? Any info is appreciated.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to check gene properties (eg GC, repetitiveness, average expression level)?
written by bioinfo •
-
Wobbly average GC profile spanning entire R1 read length; what the cause and is it bad
written by William
-
Cutting fasta file into chunks?
written by zhousun21 •
-
Generating genesets matched for GC content and length
written by gradstudentNew •
-
GC bias within a gene?
written by annen •
-
GC content and average fold coverage of contigs
written by BAGeno
-
Difficult to Sequence Regions for WGS
written by gmcinnes •
-
Batch effects in sequencing data
written by Andrew •
-
Why Cannot Limma Package Do Differential Expression Between Two Samples Without Replication?
written by shangzhong0619 •
-
Rscript To Heatmap The Overall Gc/At Content Of Several Sequences
written by daniel.soronellas
You've got bigger problems than GC content if you're trying to compare RNASeq between two different genomes. What is the experimental design?
Sorry, I guess I was unclear. They are samples from the same species, just two groups- a control group and a treated group. If I look at the reads from the control group, they have a GC content of 44% (+/- 2) compared to my treated group that is 41% (+/- 3). Is this within tolerance for looking at differential expression between treated and control or something that I need to account for?
That looks like a pretty big difference to me - which I wouldn't expect for a library between members of the same species. How many replicates do you have per group?
That's what I thought... I have 5 replicates per group.
Have you used the
cqnplot()function from the cqn package? That'd be more useful than comparing averages.I have not.. I will try that and see what it looks like. Thanks.