HISAT2/StringTie/Ballgown for Novel Transcript Isoform Differential Expression
Entering edit mode
ra711 • 0
Last seen 3.8 years ago

I ran some mm9 data through the HISAT2/StringTie/Ballgown RNA-seq pipeline recently and went through the listed protocol HISAT2/StringTie/Ballgown Protocol. When I used Ballgown, I was not able to identify any novel transcript isoforms as compared with Table 3 from data analyzed in the linked paper. The first thing I tried was to use no filter, but that did not work. I am wondering if there was an issue in how I fed the data into Ballgown in terms of my experimental conditions?

My experiment consists of two experimental levels. The first being a genotype of SMAD4 KO or WT. The second being a treatment of either TGF-beta, BMP, or Ctrl. In total, I have six experimental groups.

Right now I have the following lines of code when trying to analyze the differential expression of the genes and transcripts.

DE by transcript

resultstranscripts <- stattest(bgmm9, feature='transcript', covariate='genotype', adjustvars=c('treatment'), getFC=TRUE, meas='FPKM')

DE by gene

resultsgenes <- stattest(bgmm9, feature='gene', covariate='genotype', adjustvars=c('treatment'), getFC=TRUE, meas='FPKM')

I appreciate any advice or guidance. Please let me know if you think there might be any issues besides the experimental grouping/leveling that is preventing me from getting novel transcript isoforms as found in Table 3 of the linked Nature protocol paper. Thank you so much.

isoform RNA-seq hisat2 stringtie ballgown • 1.4k views

Login before adding your answer.

Traffic: 971 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6