Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
crAss for comparative metagenomics using cross-assembly dutilh Bioinformatics 19 12-08-2015 08:07 AM
who coined RNAseq? RNAseq as an alignment first approach only brachysclereid Bioinformatics 3 01-10-2012 12:17 PM
Merging samples to get cross-experiment variation cbeck Bioinformatics 0 09-29-2011 07:56 AM
Best Cross species aligner Khanjan Bioinformatics 0 02-07-2011 06:59 AM

Thread Tools
Old 04-11-2012, 02:12 AM   #1
Location: Cambridge UK

Join Date: Aug 2011
Posts: 14
Default cross species RNAseq

Hi All,

I have paired-end RNAseq data from the same tissue in two species (human and rat), which I want to compare. It seems to me this is fraught with issues:
  1. The annotation quality for the two species is different and so the reads determined to align to orthologous features will likely be different. It is often the UTRs which are not well annotated (it seems), so perhaps excluding reads that bind to them is a way forward?
  2. You are wholly reliant on the quality of the orthologue mapping.
  3. Orthologous *genes* are mapped, not transcripts, meaning that a representative value for the gene needs to be derived from (possibly) multiple expressed variants. One could simply sum read counts from all the exons (with htseq-count) for this?
  4. Count based tests such as DEseq (after mapping orthologous gene pairs) seems inappropriate given the differences in gene/transcript length between species etc. However, if I calculate RPKM in each species what would an appropriate test for differential expression be?
  5. I suspect there are others, but this seems enough for now!

I wonder if anyone else has wrestled with these and has any useful advice?

sidderb is offline   Reply With Quote

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 01:28 AM.

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