Question: RNA-Seq analysis for 2 genes
0
gravatar for zaid.barcelona
16 months ago by
zaid.barcelona0 wrote:

Hi, I did RNA-seq for 20 samples and I am interested in the differential expression of normal transcripts and novel transcripts within just 2 genes using galaxy. Do I have to map my reads to the whole reference genome and what tools can I use beside cuffdiff. in the end, i want to compare 10 of my normal samples transcripts with each of the 10 mutated samples (each has different mutations) thank you in advanced,

rna-seq • 464 views
ADD COMMENTlink modified 16 months ago by Devon Ryan1.9k • written 16 months ago by zaid.barcelona0
0
gravatar for Devon Ryan
16 months ago by
Devon Ryan1.9k
Germany
Devon Ryan1.9k wrote:

If you only care about 2 genes it would have been vastly cheaper and faster to do qPCR.

Anyway, yes, you must align the data to the entire genome (or transcriptome). You're better off with featureCounts and DEseq2 rather than cuff-anything. Analyse all of the genes but if you only care about 2 genes, then just take the unadjusted p-values for them.

ADD COMMENTlink written 16 months ago by Devon Ryan1.9k

Thanks Ryan for your answer. I am doing it for my thesis to prove that RNA-seq can be used as a method to detect Splice site mutations. Do you recommend any other tools that can detect the impact of the splicing mutations in BRCA1 and BRCA2? Thank you advanced

ADD REPLYlink written 15 months ago by zaid.barcelona0

It'd make sense to quantify transcripts with something like salmon and then look for isoform switching in DESeq2. That or using DEXseq should presumably work.

ADD REPLYlink written 15 months ago by Devon Ryan1.9k
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 16.09
Traffic: 175 users visited in the last hour