SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Inquiry: Educational Resources for Data Simulation in R for DE Analysis hyates Bioinformatics 2 02-11-2014 11:44 AM
Samtools multisample SNP caller amcrisan Bioinformatics 3 01-06-2014 05:57 PM
Which GATK step should be multisample? dr_sson Bioinformatics 0 05-10-2013 03:12 AM
Filtering multisample vcf files on DP DavyK Bioinformatics 3 01-17-2013 09:39 AM
Gatk multiSample realignment and recalibration seq_GA Bioinformatics 5 06-15-2011 01:02 AM

Reply
 
Thread Tools
Old 03-19-2014, 05:47 AM   #1
edge
Senior Member
 
Location: China

Join Date: Sep 2009
Posts: 199
Default Inquiry about using Cuffdiff for multisample analysis

Hi,

I have 10 RNA-seq samples now.
Samples Times (hours after inoculation)
Susceptible NIL 0
Susceptible NIL 3
Susceptible NIL 6
Susceptible NIL 9
Susceptible NIL 12
Resistant NIL 0
Resistant NIL 3
Resistant NIL 6
Resistant NIL 9
Resistant NIL 12

My main objective is to identify and characterize resistant genes associated to disease/bacteria-infection in my plant samples.
I have a reference genome and reference annotation of my plant samples.
From my RNA-seq data, I don't have any biological replicate or technical replicate.

As shown in cufflink link, http://cufflinks.cbcb.umd.edu/tutorial.html
I have run the cuffdiff as shown in the "Differential analysis with gene and transcript discovery" section.
cuffdiff merged_asm/merged.gtf S_NIL_0.bam S_NIL_3.bam S_NIL_6.bam S_NIL_9.bam R_NIL_0.bam R_NIL_3.bam R_NIL_6.bam R_NIL_9.bam

From the above command, it generates a list of output file :
var_model.info
tss_groups.read_group_tracking
tss_groups.fpkm_tracking
tss_groups.count_tracking
tss_group_exp.diff
splicing.diff
run.info
read_groups.info
promoters.diff
isoforms.read_group_tracking
isoforms.fpkm_tracking
isoforms.count_tracking
isoform_exp.diff
genes.read_group_tracking
genes.fpkm_tracking
genes.count_tracking
gene_exp.diff
cds.read_group_tracking
cds.fpkm_tracking
cds_exp.diff
cds.diff
cds.count_tracking
bias_params.info

Do anybody has similar experience to use Cuffdiff for multisample analysis ?
I hope can plot some PCA plot, heat-map, MA-plot, etc in order to identify and characterize resistant genes associated to disease/bacteria-infection in my plant samples.

Below is one of the journal I refer so far, http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3223079/
I have a similar experiment design as the journal above.
But I have no much idea how to use Cuffdiff output for further analysis.

Thanks for any advice.
edge is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:57 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO