SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
using NGS, what is the best miRNA expression normalization method? Giorgio C Bioinformatics 4 12-07-2011 06:32 AM
miRNA's expression normalization Giorgio C Bioinformatics 0 12-06-2011 10:54 AM
ChIP-Seq: A regression analysis of gene expression in ES cells reveals two gene class Newsbot! Literature Watch 0 03-05-2011 02:01 AM
Differential gene expression of gene clusters anjana.vr RNA Sequencing 1 10-28-2010 10:33 AM
A scaling normalization method for differential expression analysis of RNA-seq data severin Literature Watch 1 09-09-2010 11:09 PM

Reply
 
Thread Tools
Old 08-04-2011, 01:31 AM   #1
cedance
Senior Member
 
Location: Germany

Join Date: Feb 2011
Posts: 108
Default Gene expression normalization

Hi,
I am working on alternative splicing (AS) events on 4 different tomato species. I am trying to find "interesting" AS events. What do I mean by that?

Lets say for example, the AS event is "exon-skipping" (ES). After mapping RNA-Seq reads to tomato genome, I look for "junctions" (or intron coordinates) where they are "spliced" normally - normal junctions (NJ) and where the 3prime end Exon is skipped. So, for every junction, I have a count of reads that map to the junction normally (exactly where the intron is and supposed splicing should occur) and count of reads where at the same junction an ES event had occurred (the 3prime exon is skipped). At the end I have a table like this for each junction (I already remove where there is no ES event in ALL 4 species).

Junction 1:
S1 S2 S3 S4
ES 10 0 27 0
NJ 95 20 50 380

Then I do a fisher-test on this 2*4 table and correct for multiple-testing using Benjamini-Hochberg method (from R multtest package) to obtain those events that are significantly different across species.

Now of course the question is, what if 1) the gene where this junction (or intron) belongs is over- (or under-) expressed between these species. For ex: S2 has only a total of 20 reads mapped. 2) How about the number of reads for these species that was sequenced? 3) what about gene length? ( as the transcript abundance is also found to be positively correlated with gene length).

So, I have to somehow normalize this data. So far, with the exception of RPKM (which I am not convinced as an appropriate measure), all other methods were about finding differential expression of genes (and demand 2 or more samples), for ex: quantile normalization, TMM, the edgeR package etc. However, I would like to normalize gene expression in each of these samples.

Does anyone have an idea how to go about it? I would be very grateful for any ideas.

Thank you!

Last edited by cedance; 08-04-2011 at 01:50 AM.
cedance is offline   Reply With Quote
Reply

Tags
fisher test, gene expression, multiple testing, normalization, rna-seq

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 12:55 PM.


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