Question: Extract Transcript sequence from fastq file (Plant RNA-Seq data)
0
gravatar for hafiz.talhamalik
3 months ago by
Pakistan
hafiz.talhamalik230 wrote:

I have a plant RNA-Seq data and I want to extract the transcript sequence of the some transcripts. How can I do that ??

rna-seq • 130 views
ADD COMMENTlink modified 3 months ago by kristoffer.vittingseerup3.0k • written 3 months ago by hafiz.talhamalik230
1

@ hafiz.talhamalik

Is this plant genome already sequenced and annotated? If so, go for reference based assembly (with reference sequence) and post alignment, you can extract transcripts (sequence) of interest using annotation file. For aligning RNAseq data, use splice aware aligners such as hisat2, star.

ADD REPLYlink modified 3 months ago • written 3 months ago by cpad011212k

yes it's known one Brassica Napus. I tried reference based assembly using trinity but assembly results were not good..!! any other good assembler ??

ADD REPLYlink written 3 months ago by hafiz.talhamalik230
2

I guess Trinity is for de novo assembly. Since reference genome and annotation are available, use reference based assembly. Authors in this paper (on B. napus transcriptome) https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-015-2062-7 used tophat. I would suggest hisat2 or star. You can use reference genome and annotation from here: https://plants.ensembl.org/Brassica_napus/Info/Index @ hafiz.talhamalik

ADD REPLYlink modified 3 months ago • written 3 months ago by cpad011212k

You can extract transcript sequences directly from ensembl plants biomart, you don't need RNA-seq data for that. In case it is transgenic with fused transcripts or something like that, then you can use blast to extract reads from fastq file and then assemble them.

ADD REPLYlink written 3 months ago by ashish300

ok let me try that.. and any idea of tool for annotation of plant assembly ?

ADD REPLYlink modified 3 months ago • written 3 months ago by hafiz.talhamalik230

MAKER-P is good.

ADD REPLYlink written 3 months ago by ashish300
0
gravatar for kristoffer.vittingseerup
3 months ago by
European Union
kristoffer.vittingseerup3.0k wrote:

If you have RNAseq you can also do a guided transcriptome assembly - basically look for transcripts in your data that does not appear in the annotation database - you can read more about considerations and tools here.

ADD COMMENTlink written 3 months ago by kristoffer.vittingseerup3.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1371 users visited in the last hour