RNA-seq analysis between 2 closely related strains of the same species

Question: RNA-seq analysis between 2 closely related strains of the same species

2.0 years ago by

Hi, I am working on an RNA-seq analysis of a wild type and a gamma-irradiated mutant of a non-model organism. The aim is to identify differentially expressed genes between them. So, what I have done is generating 2 separate de novo assemblies, identifying common sequences between the 2 with at least 90% identity using cd-hit-est-2d, and using it as a mapping reference to do DE.

My question is whether my current workflow is fine to be continued or there is any generally accepted workflow to apply in my case? What do you think? Thanks.

rna-seq next-gen assembly • 659 views

ADD COMMENT • link •

modified 2.0 years ago by kristoffer.vittingseerup • 2.0k • written 2.0 years ago by GiantSilverSoy • 0

2.0 years ago by

kristoffer.vittingseerup • 2.0k

European Union

kristoffer.vittingseerup • 2.0k wrote:

I'm assuming you dont have a reference genome due to the non-model organisme comment (if you have you should take a different approach).

The other possible approach you could do would be to do a de-novo assembly based on the pooled data and then quantify that in each of your samples.

Both approaches have problems: The drawback of your solutions is that you rely on the % identical cutoff and you might get a 1:many or many:many relationships that are hard to untangle. The drawback of my suggestion is that you might assemble something that is not pressent in the actual samples.

I think i like mine a little better because you can be certain that it is the same transcript/gene that you quantify in both samles.

ADD COMMENT • link written 2.0 years ago by kristoffer.vittingseerup • 2.0k

Comparing expression when you map different samples to different assemblies is just too hard, so I'd second making a single coassembly. But an alternative approach to assembling all the reads together would be:

1) Assemble the control.
2) Map irradiated sample to control assembly, with fairly loose tolerance for mapping to account for polymorphisms.
3) Assemble unmapped reads.
4) Combine the assemblies and use that as your reference.

That might reduce the number of spurious assembled sequences that are just due to polymorphisms between the samples.

ADD REPLY • link modified 2.0 years ago • written 2.0 years ago by Brian Bushnell ♦ 16k

Similar posts • Search »