Thank you for your reply! I read that section and it brings up Cufflinks as well as iReckon, SLIDE, and StringTie. However, Montebello seems the most applicable since it couples "isoform discovery and quantification". Do you have any preference or experience with these methods? It also seems that most of these methods require sam files which require much more rigorous processing, which could prove to be a problem.
Identifying isoforms de novo
Hello biostars. I recently ran rsem on my RNA-seq data and came back with unusual isoform results. In the IsoPct (isoform percentage) column I only have 100 or 0, meaning that there was no identifiable gene/isoform or that there was only one isoform of the gene. However, I find this highly unlikely and I believe that this is due to the fact that isoforms were not directly annotated in my reference genome or actual data. How would I go about finding the isoforms for my genes using fasta files (that have individual genes not divided by chromosome) and bam files for my varying conditions? I do not have isoforms annotated so this will have to be de novo.
I have tried various programs available such as flipflop which requires sam files, I have bam files that are much too big to convert (>10 gb). Also, I have tried GESS which requires fasta files for each chromosome in the reference genome (I only have a reference genome with all of the individual genes not divided by chromosomes). I used hisat and HTseq to retrieve my bam files and gene counts.
Much appreciated.
• 3,155 views
•
link
updated
by
Istvan Albert
•
written
by
jjrin •
•
written
by
jjrin •
1 answer
In this case you may not need a full de-novo transcript assembly since it appears that the genome is already known. This makes the process a little easier. Traditionally people used Cufflinks to do this, but the tool has been falling out of favor. For more current methods consult the literature.
See for example the Transcript Discovery section in
A survey of best practices for RNA-seq data analysis, Genome Biol. 2016
• 0 views
•
link
written
by
Istvan Albert
Aligning against known transcripts can't tell you much about unknown transcripts.
In general, to discover new isoforms you will need alignments against the genome. If those are not available then you would need to perform a de-novo transcriptome assembly of the reads.
• 0 views
•
link
written
by
Istvan Albert
What if i still have the whole genome fasta file available? GESS doesn't ask for the gtf annotation file so would it be beneficial to work with the whole genome data to find isoforms or not? I'm currently working on slide which isn't working as well as my references are invalid '6L'.
Log in to answer this question.
More posts like this
-
Different annotations for different isoforms from the same gene
written by Chikae •
-
Confusion about StringTie
written by DGTool •
-
Inputs required for CNVnator
written by sjlee4019 •
-
RSEM isoform.results outputs for edgeR
written by weichengz •
-
Unknown Error running GESS
written by jjrin •
-
Invalid reference '6L' even though naming of BAM files, fasta files, and gtf files are the same?
written by jjrin •
-
Why does cuffdiff test differential isoform expression for genes with only one annotated transcript?
written by unsupervised_learner •
-
mapping RNAseq reads to a genbank isoforme that is not in Reference Annotation
written by fta.mirzadeh •
-
AUGUSTUS principal isoform
written by jomaco
-
How does MISO quantify genes with single isoform
written by Anil Kesarwani
Was this a de novo assembly?
Not entirely, we had a reference genome (Taejoon lab) for assembly and for identifying exons, CDS, etc. However, isoforms weren't included and we would like to find a way to identify them.