Question: De novo transcriptome assembly using Trinity
0
gravatar for 1234anjalianjali1234
9 months ago by
India
1234anjalianjali123430 wrote:

Hi guys,

I have a total of 12 samples [6 infected sample (number- 7,8,9,10,11 and 12) of different time-point with their 6 control samples (number- 1,2,3,4,5 and 6)] from RNAseq, paired end (Illumina; GAIIx).

My Sample was distributed in 8 Lanes. Details:

LANE      SAMPLE
LANE_1: 7,8,9,10          [Control samples]
LANE_2: 1,2,3,4,5,6       [Infected samples]
LANE_3: 11,12,7,8         [Control samples]
LANE_4: PhiX control lane
LANE_5: 9,10,11,12        [Control samples]
LANE_6: 1,2,3,4           [Infected samples]
LANE_7: 5,6,1,2           [Infected samples]
LANE_8: 3,4,5,6           [Infected samples]

The data provided to me was in bcl format, which i have converted them into fastq format using bcltofastq software. Now, I have 16 fastq files, both forward and reverse:

Right Reads    Left Reads
LANE_1_R2.fq   LANE_1_R1.fq
LANE_2_R2.fq   LANE_2_R1.fq
LANE_3_R2.fq   LANE_3_R1.fq
LANE_4_R2.fq   LANE_4_R1.fq
LANE_5_R2.fq   LANE_5_R1.fq
LANE_6_R2.fq   LANE_6_R1.fq
LANE_7_R2.fq   LANE_7_R1.fq
LANE_8_R2.fq   LANE_8_R1.fq

My problems:

  1. For the denovo assembly do i need to feed all the R1 reads as left and R2 reads as right in TRINITY and do I need to remove PhiX Lane.

  2. How to align reads from each sample to the reference because i dont have seperate fastq file of each sample.

Many thanks in advance!

rna-seq illumina gaiix • 265 views
ADD COMMENTlink modified 9 months ago by JC9.1k • written 9 months ago by 1234anjalianjali123430
0
gravatar for JC
9 months ago by
JC9.1k
Mexico
JC9.1k wrote:
  1. In general yes, you can pass all the transcriptome to Trinity, check the "normalize" parameters to reduce the large assembly memory issue. Alternatively, you can assemble samples per treatment (I generally see better assembly because isoforms are less per sample) and then merge your assemblies.

  2. The Illumina Fastq IDs should indicate the index (https://en.wikipedia.org/wiki/FASTQ_format#Illumina_sequence_identifiers), and you need to know which index belongs to each library to demultiplex your files.

ADD COMMENTlink written 9 months ago by JC9.1k

Thankyou JC,

I will try this.

ADD REPLYlink written 9 months ago by 1234anjalianjali123430
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: 1755 users visited in the last hour