SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
qubit and picogreen discrepancies anjama Sample Prep / Library Generation 3 07-02-2017 07:07 PM
pico green discrepancies schizosaccharomyces Sample Prep / Library Generation 5 10-17-2013 09:37 AM
Qubit vs Nanodrop discrepancies?? amylee86 General 1 07-10-2013 08:50 AM
tophat / bowtie mapping discrepancies vineeth_s Bioinformatics 3 06-10-2011 05:44 AM

Reply
 
Thread Tools
Old 09-18-2017, 08:38 AM   #1
biodiverse
Junior Member
 
Location: Pennsylvania

Join Date: Aug 2017
Posts: 3
Question Discrepancies in Demultiplexing programs - help request

I have Illumina PE data from a MiSeq run and am finding large discrepancies in the number of reads recovered while demultiplexing.
Originally, the fastq files provided with the run had very few reads per individual and a large (>4Gb) Undetermined file. I investigated the DemultiplexSummary provided with the run and found that the top 30 indexes found were in fact my indexes but for some reason (??) were not de-multiplexed correctly.
Using the Fastx-barcode splitter (HammonLab), allowing a 1nucleotide mismatch, I was able to recover a small number of additional reads from the Undetermined file but not nearly the number stated in the summary.
I then inputted the raw data into Geneious and after trimming the adapters with the bbduk input I de-multiplexed and found significantly more reads but still only approximately half of the number that should be present.

Some Numbers:
Barcode CAAAAG/CTTTTG
DemultiplexSummary: 1,119,171 reads
Fastq file unaltered from run: 1, 428 reads
Fastx Barcode Splitter (on undetermined file): 118,031 reads
Geneious: 574,498 reads

My questions are:

1 - why are the different programs giving such disparate results?

2- Am I misunderstanding the orientation of the barcodes in the reads and thus perhaps searching for them incorrectly? It is my understanding that the "adapter - barcode - read" order should have the barcode as the first 6 bases in the R1 read after adapter trimming (in my example CAAAAG). R2 should not have(?) the barcode - or I should not have to search for a barcode in R2 anyways as I have paired the data? I recovered the reads in Geneious using CTTTTG as that was the index listed in the demultiplexsummary but my barcode as listed in my primer order was CAAAAG so I am concerned that I am misunderstanding a fundamental piece of the puzzle here.

3 - Most importantly - how do I recover the 1million+ reads?!!?
Attached Files
File Type: txt Barcodes.txt (475 Bytes, 9 views)
File Type: txt DemultiplexSummaryF1L1.txt (11.9 KB, 6 views)
biodiverse is offline   Reply With Quote
Old 09-19-2017, 06:08 AM   #2
kmcarr
Senior Member
 
Location: USA, Midwest

Join Date: May 2008
Posts: 1,135
Default

Quote:
Originally Posted by biodiverse View Post
I have Illumina PE data from a MiSeq run and am finding large discrepancies in the number of reads recovered while demultiplexing.
Originally, the fastq files provided with the run had very few reads per individual and a large (>4Gb) Undetermined file. I investigated the DemultiplexSummary provided with the run and found that the top 30 indexes found were in fact my indexes but for some reason (??) were not de-multiplexed correctly.
Using the Fastx-barcode splitter (HammonLab), allowing a 1nucleotide mismatch, I was able to recover a small number of additional reads from the Undetermined file but not nearly the number stated in the summary.
I then inputted the raw data into Geneious and after trimming the adapters with the bbduk input I de-multiplexed and found significantly more reads but still only approximately half of the number that should be present.

Some Numbers:
Barcode CAAAAG/CTTTTG
DemultiplexSummary: 1,119,171 reads
Fastq file unaltered from run: 1, 428 reads
Fastx Barcode Splitter (on undetermined file): 118,031 reads
Geneious: 574,498 reads

My questions are:

1 - why are the different programs giving such disparate results?

2- Am I misunderstanding the orientation of the barcodes in the reads and thus perhaps searching for them incorrectly? It is my understanding that the "adapter - barcode - read" order should have the barcode as the first 6 bases in the R1 read after adapter trimming (in my example CAAAAG). R2 should not have(?) the barcode - or I should not have to search for a barcode in R2 anyways as I have paired the data? I recovered the reads in Geneious using CTTTTG as that was the index listed in the demultiplexsummary but my barcode as listed in my primer order was CAAAAG so I am concerned that I am misunderstanding a fundamental piece of the puzzle here.

3 - Most importantly - how do I recover the 1million+ reads?!!?
BD,

First things first; assuming that your libraries are standard Illumina design your barcodes are NOT part of R1. In the standard Illumina library design and run configuration the index read is completely separate from R1 or R2. Unless your sequence provider also gave you the index read file (would have something like "I1" in place of R1 or R2 in the file name) then you don't have the information needed to demultiplex your reads.

Can you provide a copy of the SampleSheet.csv file that was used for the MiSeq run? That will show the indexes actually used by the MiSeq to demultiplex your data. It sounds very likely that the SampleSheet.csv had the indexes entered in the wrong orientation.

To determine the proper orientation we would need to see an example of your completed adapter sequences, or know what kit was used to construct your libraries.
kmcarr is offline   Reply With Quote
Old 09-19-2017, 08:12 AM   #3
biodiverse
Junior Member
 
Location: Pennsylvania

Join Date: Aug 2017
Posts: 3
Default

Thank you so much for your reply. I've attached both my sample sheet and my illumina adapter list with the barcodes. I didn't prepare the libraries with a kit but it was a standard prep.
Could you please explain how the reads and barcodes are oriented? I have searched the literature for this but am clearly misunderstanding something.
Attached Files
File Type: txt SampleSheet.txt (1.6 KB, 6 views)
File Type: txt Indexes_and_adapers.txt (2.1 KB, 3 views)
biodiverse is offline   Reply With Quote
Old 09-19-2017, 08:44 AM   #4
jdk787
josh kinman
 
Location: Austin

Join Date: Apr 2014
Posts: 57
Default

It looks like the index sequences in your sample sheet are in the wrong orientation. The demux report shows that the reverse complement of your sample sheet index sequences is what was found in your run.


Here is a link to a useful document that shows all of the reads used by Illumina sequencers.
https://support.illumina.com/content...5057455-03.pdf
__________________
Josh Kinman
jdk787 is offline   Reply With Quote
Old 09-19-2017, 08:51 AM   #5
kmcarr
Senior Member
 
Location: USA, Midwest

Join Date: May 2008
Posts: 1,135
Default

Quote:
Originally Posted by biodiverse View Post
Thank you so much for your reply. I've attached both my sample sheet and my illumina adapter list with the barcodes. I didn't prepare the libraries with a kit but it was a standard prep.
Could you please explain how the reads and barcodes are oriented? I have searched the literature for this but am clearly misunderstanding something.
As suspected the index sequences in the SampleSheet are in the wrong orientation.

Code:
Index1 p7 end primer showing index read (i7) primer annealed 

                                  Index (i7) Read Sequencing Primer
                               <-CACTGACCTCAAGTCTGCACACGAGAAGGCTAG-5
5'-CAAGCAGAAGACGGCATACGAGATatcacgGTGACTGGAGTTCAGACGTGT-3'

(Obviously your full library fragment will extend beyond the 3' end of 
Index1 primer. Truncated for clarity.)

Index1 will be read as CGTGAT
To fix this problem you will need to correct the SampleSheet.csv file by reverse complementing the indexes. Then repeat the Bcl2Fastq conversion with demultiplexing. I imagine you will need to get your sequencing provider involved since they have the original base call (.bcl) files with the index read data.
kmcarr is offline   Reply With Quote
Old 09-19-2017, 12:12 PM   #6
biodiverse
Junior Member
 
Location: Pennsylvania

Join Date: Aug 2017
Posts: 3
Default

Thank you so much for your help, I will do this!
biodiverse is offline   Reply With Quote
Reply

Tags
demultiplex, fastx toolkit barcode, geneious, illumina fastq

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 01:59 AM.


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