This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Need tools recommendation for plant RNAseq data

Hi,

I have RNA data from 10 samples of Cannabis plant. I am planning to do the following analyses with these samples

-) De novo Assembly (ex: Trinity, Velvet, SOAPdenovo, etc. Which is recommended for plants?)

-) Variant calling ( does GATK support plant variant calling?

-) Functional annotation like COG, GO, KEGG, swissprot, NR database (ex: BLAST2GO tool, any other recommendation)

-) Differential gene expression (DESEQ2 or EDGER or LIMMA, which one is recommended for plant dataset?)

rna-seq plants assembly variant calling
• 2,335 views • link
updated by apeltzer apeltzer • written by bioinforesearchquestions bioinforesearchquestions

There is already a (partially) similar paper, look at its methods to draw some ideas. If you follow the paper methodology, upgrade the tools used (hisat - stringtie - ballgown onstead of tophat - cufflinks). There is a Cannabis genome, unfortunately, it seems it is not annotated.

• 0 views • link
written by h.mon h.mon

1 answer

First of all: What kind of data do you have? Is it short read data? Is it long reads (e.g. PacBio SMRT/Nanopore?).

  • If you have short read data, try SOAPDenovo for assembly. In case you have long read data, I had pretty good experiences with "Canu".

For Variant calling I assume that GATK should work - but didn't use it on plant data yet.

  • DESEQ2 should work fine, just make sure to use a proper reference genome when mapping reads and then annotate properly before doing the expression analysis.
• 0 views • link
written by apeltzer apeltzer

Log in to answer this question.

More posts like this