Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
progressiveMauve Orthologs Output Adjuvant Bioinformatics 0 02-15-2012 02:18 PM
RNA-Seq: ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count da Newsbot! Literature Watch 0 11-18-2011 03:20 AM
RNA-Seq: Comparative Analysis of RNA-Seq Alignment Algorithms and the RNA-Seq Unified Newsbot! Literature Watch 3 07-31-2011 08:08 PM
RNA-Seq: RNA-Seq Analysis of Gene Expression and Alternative Splicing by Double-Rando Newsbot! Literature Watch 0 03-03-2011 03:00 AM
Suggestions for RNA-seq quantitation? jjvalenz5 Bioinformatics 4 03-01-2011 01:28 PM

Thread Tools
Old 05-09-2012, 03:25 PM   #1
Junior Member
Location: Monterey, CA

Join Date: Mar 2010
Posts: 1
Default suggestions for RNA-Seq DE analysis of orthologs from 2 different references/taxa

Hi All,
I'm currently trying to do differential expression analysis comparing expression differences at orthologous loci between two closely related taxa that have been aligned to separate references. The issue is that I am hoping to directly compare two count values (# of reads uniquely aligned) that were obtained from alignment to separate reference sequences of different lengths.

-sample 1_taxa1_contigX_ortho1_contiglength1257 has a count of 1200
-sample 2_taxa2_contigY_ortho1_contiglength1989 has a count of 57
can I do differential expression analysis to see if these two expression levels are constitutively different?

So far there are two separate issues:
1-Normalizing for sequencing depth-
a-Most DEG analysis programs (e.g. DESeq, edgeR, etc.) do a within experiment normalization for sequencing depth, however, I have two separate references with a subset of orthologous loci so normalizing different counts at the subset orthologs only would ignore the counts at non-orthologs. This doesn't seem appropriate for other reasons as well.

2-Accounting for differences in contig length-
a-Since contigX for ortholog1 from taxa1 is a different length than contigY for ortholog1 from taxa2, there could be some length bias in alignment efficiency that could mess up the differential expression analysis.

My proposed solution would be to resurrect RPKM normalization for this particular case because it independently normalizes each sample based on total number of reads and the different contig lengths. I would do this:
1- RPKM normalization for all reference contigs for taxa1 and separately for all reference contigs for taxa2
2- Pull out the counts for the ortholog subset
3- Combine the data from each taxa at the orthologs
4- Do some type of differential expression analysis comparing counts from taxa1 at ortholog1 to taxa2 at ortholog1

However, I am not aware of any DE package that will take RPKM normalized data. I've managed to trick DESeq into doing it by using a sizeFactor of 1 after RPKMing my data and rounding to integers but this seems less than rigorous and/or true to the statistical framework.

Any suggestions for other methods to account for the different sequencing depths and different contig lengths, as well as possible statistical frameworks to analyze custom normalized data would be much appreciated.

This is for a microbe system where we are trying to investigate expression differences among different genotypic populations of organisms that are divergent enough to warrant separate reference assemblies but have a well matched subset of orthologous loci.

Any assistance or links to other posts that might be helpful would be much appreciated.
Many thanks,
HopkinsDan is offline   Reply With Quote

differential expression, gene expression pattern, orthologs, rna-seq, rpkm

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 02:49 AM.

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