SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Start analysis from DESeq2-normalized read counts? friducha Bioinformatics 1 07-04-2015 05:40 AM
DESeq and non-raw counts (but also not rpkm normalized) wingless Bioinformatics 0 06-26-2015 11:11 AM
Getting normalized counts IsBeth Bioinformatics 2 05-31-2014 09:20 AM
edgeR mean normalized reads counts for each gene antoza Introductions 2 02-11-2014 08:22 AM
converting bam files to non-normalized read counts lpn Bioinformatics 4 10-09-2012 07:52 PM

Reply
 
Thread Tools
Old 06-28-2018, 09:04 AM   #1
seqprone
Junior Member
 
Location: USA

Join Date: Oct 2015
Posts: 8
Default What can you learn about your samples just by looking at normalized read counts?

Hi everyone,

My study involves multiple cell lineages which were isolated from the mammary gland and sorted based on surface markers. The FACS data 100% agrees with published literature and what we expected to see.

I isolated RNA from all these cell lineages and treatment groups and sent it for sequencing.

I recently finished the RNASeq analysis where I modified a script from another lab and made it work. I showed my script and outputs to several people who know better and they all agreed that the script is working flawlessly. I followed up with a GSEA analysis which looks ok to me.

The problem is, that the DEG list and GSEA analysis doesn't show all of the genes/pathways that we expected to see.

The immediate question was if my sorted population was contaminated or and/or the samples got mixed up/flipped at the sequencing core.

The sequencing core obviously denied the chances of samples flipping or being mislabeled and I want to believe them.

I suggested that, if the purity of cell lineage is a concern, then let me run a qRTPCR for the lineage-specific genes. Looking at fold change values of the KRT5, KRT14, KRT18, ESR1 etc would tell me if they are basal or luminal population.

To this, some of the people responded that why don't I look at the normalized gene counts and see the expression level of my lineage-specific markers?

My concern is, can I really do that? Can I look at the normalized read counts and say that these samples are basal cells because they have a higher count or reads of KRT5 and these are luminal cells because they have a higher count of KRT18?

Obviously, I know very little about the statistics that goes into creating all these data files but I was under the impression that the normalized read counts don't give you enough information unless you put it through an appropriate statistical analysis (like DESeq or EdgeR)

I'd appreciate any feedback or a link to any text that might help.
seqprone is offline   Reply With Quote
Old 07-05-2018, 09:59 PM   #2
Dario1984
Senior Member
 
Location: Sydney, Australia

Join Date: Jun 2011
Posts: 163
Default

You can't simply use GSEA on short read RNA-seq data counts. The number of alignments counted for a gene is proportional to its abundance and its length. You need to normalise for gene length.
Dario1984 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:48 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO