Same sequencing sample in multiple lanes. How to analyse it?
1
0
Entering edit mode
2.0 years ago
Federico • 0

Hi,

I have the following samples:

Lane 1:

  • 184631_S1_L001_trimmed_R1.fastq
  • 184631_S1_L001_trimmed_R2.fastq

Lane 2:

  • 184631_S1_L001_trimmed_R1.fastq
  • 184631_S1_L001_trimmed_R2.fastq

I would like to align these fastq files to the reference library. Bowtie2 is the tool I am using. As you can see, the same sample was spread across 2 lanes in order to obtain adequate sequencing depth. The expected output is a count matrix that will be further analysed using DESeq2.

How (and when) to proper treat and merge the same sample from two different lanes?

sequencing bowtie2 alignment samtools ngs • 1.8k views
ADD COMMENT
1
Entering edit mode

Are they technical replicates? If not, you can merge them before mapping with the simple cat command.

ADD REPLY
0
Entering edit mode

They are not. The shown samples are from the same multiplex, and the multiplex was then sequenced on 2 lanes. What are you suggesting is to merge them before the alignment? R1.fastq from both lanes and same for R2.fastq?

ADD REPLY
1
Entering edit mode

Yes, merge all forward and reverse read files separately.

cat Lane1/184631_S1_L001_trimmed_R1.fastq Lane2/184631_S1_L001_trimmed_R1.fastq > final_R1.fastq
cat Lane1/184631_S1_L001_trimmed_R2.fastq Lane2/184631_S1_L001_trimmed_R2.fastq > final_R2.fastq
ADD REPLY
0
Entering edit mode

If they same from the same library prep, they absolutely are technical replicates, no matter how many different lanes they are split over.

ADD REPLY
0
Entering edit mode

Thanks a lot for your suggestion!

ADD REPLY
2
Entering edit mode
2.0 years ago

You can merge the files before and after alignment.

Merging the FASTQ files perhaps makes life a little simpler in that fewer files need to be managed.

Merging after alignments gives you a bit of a better understanding of whether the two lanes exhibit different characteristics during alignment. You also get somewhat simpler parallel processing when generating alignments (in the sense that the files are already split).

ADD COMMENT
0
Entering edit mode

Hi Istvan! Would it be fine to add the counts of samples (same sample from different lanes) from count matrix instead of merging aligned files? Logically it should not make any difference.

Best Ekta

ADD REPLY
0
Entering edit mode

yes, I think adding up the counts is also a valid approach

ADD REPLY

Login before adding your answer.

Traffic: 1136 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6