Identifying and quantifying transcript variants in RNA-seq data
Hi,
I've processed FASTQ files with nf-core/rnaseq pipeline (GRCh38, hisat2 aligner, StringTie transcript assembly & quantification). After DEG & GSEA analysis, I'm interested in looking into transcript variants. Specifically, I want to look for truncated MDK transcript (if present) - https://www.ncbi.nlm.nih.gov/nuccore/DQ323888.1 . How do I do that? I'd really appreciate for pointing me in the right direction, which tools to use etc. Is stringtie --merge followed by gffcompare a good way to do it in this case?
Thank you! Mieszko
• 1,257 views
•
link
written
by
MieszkoL •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Transcript level quantification and gene length scaling for rRNA depleted libraries
written by marc.zimmerli •
-
Current best practice for alignment and quantification for downstream analysis in DESeq2
written by yura.grabovska •
-
How to to lncRNA analysis from WTS data ?
written by Sowmya Pulapet •
-
RNA seq data analysis without control samples
written by vasanthmandla •
-
Distinguishing protein transcript encoded within lincRNA in RNA-seq DE analysis
written by sysboolean
-
How do you find novel transcripts using GFFcompare?
written by c_u
-
Stringtie with Gencode gtf file doesn't give any gene_ids or transcript_ids?
written by Biologist
-
(Modern, mid-2016) RNA-seq software pipeline
written by jp
try using one of the pseudo-aligners (kallisto/salmon) https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-021-04198-1
Otherwise you need to try like DEXseq or something and recreate your "isoforms" manually
If I'm understanding correctly, salmon/kallisto approach is limited by transcriptome used for annotation. Linked MDK transcript (DQ323888.1) is not in GRCh38.p13 and therefore it can't be detected this way, right?