Question: How to use PASA to add novel gene from RNA-seq alignment to existing gene feature files ( GFF )?
0
gravatar for Maimai
2.3 years ago by
Maimai0
Maimai0 wrote:

Hi,

I'm new for the PASA software. I've been dealing with gene annotation in drosophila recently and I faced a problem when dealing transcriptome annotation with PASA. The question is about the TRANSDECODER parameter. When I added this parameter in the step "Performing an annotation comparison and generating an updated gene set" of the script ../scripts/Launch_PASA_pipeline.pl, it didn't add novel gene comparing to the one without adding TRANSDECODER parameter. The thing is, I am wondering if it can help me to find new gene of from transcript alignment that the existing annotation don't have and add those genes to updated GFF files.

Best, Maimai

rna-seq annotation • 1.1k views
ADD COMMENTlink modified 2.3 years ago by Juke-342.1k • written 2.3 years ago by Maimai0
1
gravatar for Juke-34
2.3 years ago by
Juke-342.1k
Sweden
Juke-342.1k wrote:

Yes I confirm you it can help you to find new genes from transcript alignments as well as new isoforms for gene already annotated. Could you share your command line ?

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Juke-342.1k

Here's the command line: nohup $tool/PASApipeline-2.0.2/scripts/Launch_PASA_pipeline.pl -c annotCompare.config -A -g genome.fasta -t existing_annotation.gff3 --TRANSDECODER &

Thank you so much for ur reply~~

ADD REPLYlink written 2.3 years ago by Maimai0

You have to provide the annotation you want to update with the option --annots_gff3 (i.e the help page). You used instead the -t that is to provide the transcripts file to use to perform the new predictions (in fasta format). So obviously is missing this latest one.

ADD REPLYlink written 2.3 years ago by Juke-342.1k

I changed my command: nohup $tool/PASApipeline-2.0.2/scripts/Launch_PASA_pipeline.pl -c annotCompare.config -A -L --annots_gff3 existing.annotation.gff3 -g genome.fasta -t transcript.fasta --TRANSDECODER& But the output still have no novel gene been added, there are some revised on adding UTR though.

The option --annots_gff3, should it be added in the step of transcript alignment or annotation comparison ?

Thank you so much !

ADD REPLYlink modified 2.3 years ago • written 2.3 years ago by Maimai0

I guess it's because you have to use on top of the parameters already used, the "-R" and "--ALIGNERS blat,map" to launch the alignment pipeline to pre-process your transcript.fasta. Otherwise you transcript.fasta file is useless...

ADD REPLYlink written 2.3 years ago by Juke-342.1k
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: 1159 users visited in the last hour