Question: Using Trinity with case/control data
1
gravatar for bdeonovic
3.2 years ago by
bdeonovic120
United States
bdeonovic120 wrote:

I would like to do RNA-seq differential expression analysis on an experiment run on 8 pigs (4 cases and 4 controls, not matched). I've already done the tophat/cuffdiff approach using the pig genome. My PI also wants me to try out Trinity (without using the pig genome sequence).

I noticed that trinity can take in multiple fastq files at a time. How should I run my Trinity run(s)? Do I do 8 separate Trinity runs? Two runs for the two conditions? 

 

Thanks

sequencing rna-seq trinity • 1.1k views
ADD COMMENTlink modified 3.2 years ago by st.ph.n1.9k • written 3.2 years ago by bdeonovic120
2
gravatar for st.ph.n
3.2 years ago by
st.ph.n1.9k
Philadelphia, PA
st.ph.n1.9k wrote:

If you're looking to do a DE analysis on your Trinity results, then provide all of your fastq's at the same time. Take note if you have PE reads or SE reads. They are supplied differently to Trinity on the command line. This will create a single de novo assembly. If you assemble each sample separately, you'll have trouble with the RSEM step when aligning the raw reads back to the assembly to get abundance estimates. You can do individual runs, if you simply just want a de novo assembly for that sample. But for the DE analysis, you'll need to create a single assembly, combining all samples as inputs.

Example:

if PE: cat all your read 1's, and cat all your read 2's

Trinity --seqType fq --JM 50G --left reads_1.fq  --right reads_2.fq --CPU 6

For the DE analysis downstream you'll be able to specify replicates, and sample names.

If all fastq's are SE there is a --single option to use instead of specifying read orientation.

Check out the Trinity site: Trinity (http://trinityrnaseq.sourceforge.net/#typical_options); and for full usage there is a command line parameter: --show_full_usage

Hope this helps.

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by st.ph.n1.9k

This is what I ended up doing, except instead of cat-ing the reads together I provided a list of them (separated by a space) since Trinity says I can provide 1 or more fastq files. 

ADD REPLYlink written 3.2 years ago by bdeonovic120

You can provide a list of file names for the --left and --right command line parameters but that list of names should be a comma separated list, not space separated. I fear that in this case Trinity may have taken the first file in each list and ignored the rest.

ADD REPLYlink written 3.2 years ago by kmcarr00270

The Trinity documentation is very shoddy...This is what it says:

If paired reads:
#      --left  <string>    :left reads, one or more (separated by space)
#      --right <string>    :right reads, one or more (separated by space)
#
ADD REPLYlink written 3.2 years ago by bdeonovic120

I've always just concatenated the files in respect to read orientation.

ADD REPLYlink written 3.2 years ago by st.ph.n1.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 2.3.0
Traffic: 1003 users visited in the last hour