SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Species abundance and rarefaction in metagenomics data OliverDeusch Metagenomics 1 12-13-2012 09:17 AM
cross species RNAseq sidderb Bioinformatics 0 04-11-2012 03:12 AM
View Cross Machine Image Processing Data problem in V2.5.3 AndyChung 454 Pyrosequencing 0 05-26-2011 07:47 PM
PubMed: NGSQC: cross-platform quality analysis pipeline for deep sequencing data. Newsbot! Literature Watch 0 03-30-2011 03:20 AM
Best Cross species aligner Khanjan Bioinformatics 0 02-07-2011 07:59 AM

Reply
 
Thread Tools
Old 05-23-2013, 07:14 AM   #1
trelek2
Junior Member
 
Location: Bristol

Join Date: May 2013
Posts: 2
Default cross-species data - questions about normalization

Hi,

I have some data form various samples (cell types) in different species.
I want to compare and analyze gene expression variability across the different species.

I've plotted the average expression (tags per million 'tpm' data) for each sample. I found that the average expression is more or less the same for all samples from a given species but the average expression values vary across species greatly (for example all human samples are about 20 and all mouse samples are about 25 but all dog samples are at about 100 tmp).

I'm guessing that this is because the data is not normalized and therefore the data is not comparable before performing a normalization.

I tried using different recently developed methods for data normalization (RLE, TMM) implemented in EdgeR as well as DeSeq and normalized the data from different species (by entrez id) separately for each cell type. This however still does not give much more similar average expressions for the different species. The averages of samples from the same species have now become much more dissimilar.

The only samples that went to a similar average expression level across the different species is Universal RNA - which is a mix of different tissues rather than a specific cell type.


I'm really confused, does the above mean the normalization is fine and I shouldn't worry about the fact that average expression values between different species are dissimilar, or does this rather imply something is wrong?

Maybe I should normalize by ignoring the gene ids and just looking at the whole expression profile (negative binomial) by sorting the genes according to their expression values and then having a different normalization factor for every expression value, whichever gene it is in every sample, so that I end up with 2 negative binomial curves that look the same but the genes on the x axis will be differently ordered depending on the species?

thanks.
trelek2 is offline   Reply With Quote
Old 05-23-2013, 08:29 AM   #2
mbblack
Senior Member
 
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245
Default

Does each species have it's own controls? In other words, can each species be analyzed as an independent experiment?

If so, you could normalize each independently, analyze differential expression in each independently, and then simply compare homologous genes amongst the significantly differentially expressed genes in each species. I would suggest basing significance in each of the three independent analyses by simultaniously applying a statistical and a fold change cutoff to get the most robust differential gene lists for each.

That's how I've dealt with comparisons between, for example, rat and human hepatocytes exposed to dioxin. Set up the experiments separately and analyzed each species as an independent differential gene expression experiment, then compare the analyzed gene lists for significant homologous genes in each.

Otherwise, you have a very complex normalization situation where you may have to perform independent normalizations for each species followed by some form of scaling correction/normalization for the meta-analysis of all three.
__________________
Michael Black, Ph.D.
ScitoVation LLC. RTP, N.C.

Last edited by mbblack; 05-23-2013 at 08:32 AM.
mbblack is offline   Reply With Quote
Reply

Tags
edger, normalization, rle, species, tmm

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 10:02 PM.


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