System, method and computer accessible-medium for multiplexing base calling and/or alignment
An exemplary system, method and computer-accessible medium for multiplexing base-calling of a plurality of nucleic acid molecules in the same flow cell is provided. When multiplexing for just two nucleic acid molecules, it can operate by selecting a first base call for a first nucleic acid molecule and a second base call for a second nucleic acid molecule, after having placed them in the same flow cell and obtaining the combined raw intensity output. It can use as prior the appropriate reference genome sequences from which the nucleic acid molecules can be derived in order to create a score function, which can additionally be constrained by various penalty functions. It can derive accuracy and speed by using a branch and bound strategy as well as by performing alignment and base-calling in one step per cycle.
1. A nucleic acid sequencer comprising a single flow cell and a computer arrangement configured to perform procedures comprising:
obtain a single raw intensity output derived from a first nucleic acid sequence and a second nucleic acid sequence, wherein the first nucleic acid sequence is in a same flow cell as the second nucleic acid sequence;
obtain at least one reference sequence for a genome related to the first and second nucleic acid sequence;
generate and optimize a score function for a plurality of intensities where the score function is configured to determine: (a) that a map path is consistent with one of the first nucleic acid sequence or the second nucleic acid sequence in view of a noise component from another one of the first nucleic acid sequence or the second nucleic acid sequence, and (b) an alignment to the at least one reference sequence, subject to at least one constraint that penalizes the score function based on a fraction of bases that are not consistent with the map path;
apply the optimized score function to the combined raw intensity output to deconvolve the single raw intensity output into a plurality of constituent intensities;
map the deconvolved single raw intensity output to a set of possible alignments, the at least one reference sequence and a set of possible base calls; and
select the first base call and the second base call from the set of possible base calls.
2. The nucleic acid sequencer of claim 1 , wherein the combined raw intensity output includes raw intensity outputs received simultaneously from the first and second nucleic acid sequences.
3. The nucleic acid sequencer of claim 1 , wherein the combined raw intensity output includes a first raw intensity output received from the first nucleic acid sequence at a first time and a second raw intensity output received from the second nucleic acid sequence at a second time, and wherein the second time is later than the first time.
4. The nucleic acid sequencer of claim 1 , wherein the first nucleic acid sequence is the same as the second nucleic acid sequence.
5. The nucleic acid sequencer of claim 1 , wherein the first nucleic acid sequence is different than the second nucleic acid sequence.
6. The nucleic acid sequencer of claim 1 , wherein the at least one constraint includes at least one penalty function.
7. The nucleic acid sequencer of claim 1 , wherein the computer arrangement is further configured to generate the combined raw intensity output by introducing a plurality of primers at different cycles in a Sequencing-By-Synthesis process.
8. The nucleic acid sequencer of claim 1 , wherein the computer arrangement is further configured to obtain the at least one reference sequence independently from the combined raw intensity output.
9. The nucleic acid sequencer of claim 8 , wherein the set of possible base call includes a first set of possible base calls and a second set of possible base calls, and wherein the computer arrangement is further configured to map the combined raw intensity output to the first set of possible base calls for the first nucleic acid sequence and to the second set of possible base calls for the second nucleic acid sequence.
10. The nucleic acid sequencer of claim 1 , wherein the first and second nucleic acid sequences are deoxyribonucleic acid molecules.
11. The nucleic acid sequencer of claim 1 , wherein the first and second nucleic acid sequences are ribonucleic acid molecules.
12. The nucleic acid sequencer of claim 1 , wherein the at least one reference sequence includes a single reference sequence related to the first and second nucleic acid sequences.
13. The nucleic acid sequencer of claim 1 , wherein the at least one reference sequence includes at least two references, and wherein one of the reference sequences is related to the first nucleic acid sequence and another one of the reference sequences is related to the second nucleic acid sequence.
14. The nucleic acid sequencer of claim 1 , wherein the computer arrangement is configured to obtain the combined raw intensity from a single channel on the nucleic acid sequencer.
15. The nucleic acid sequencer of claim 1 , wherein the computer arrangement is configured to obtain (i) a first raw intensity output of the combined raw intensity output from a first channel on the nucleic acid sequencer, and (ii) a second raw intensity output of the combined raw intensity output from a second channel on the nucleic acid sequencer, and wherein the first channel is different than the second channel.
16. The nucleic acid sequencer of claim 1 , wherein the combined raw intensity output includes a first raw intensity output from the first nucleic acid sequence received simultaneously with a second raw intensity output from the second nucleic acid sequence.
17. The nucleic acid sequencer of claim 1 , wherein the first and second nucleic acid sequences are from at least one virus.
18. A system for selecting a first base call for a first nucleic acid sequence and a second base call for a second nucleic acid sequence, comprising:
a nucleic acid sequencer comprising a single flow cell, wherein the nucleic acid sequencer is configured to:
obtain a single raw intensity output derived from the first nucleic acid sequence and the second nucleic acid sequence, wherein the first nucleic acid sequence is in a same flow cell as the second nucleic acid sequence;
obtain at least one reference sequence for a genome related to the first and second nucleic acid sequence;
with the nucleic acid sequencer, generate and optimize a score function, subject to at least one constraint, based on a plurality of intensities and the at least one reference sequence, wherein the score function is configured to determine: (a) that a map path is consistent with one of the first nucleic acid sequence or the second nucleic acid sequence in view of a noise component from another one of the first nucleic acid sequence or the second nucleic acid sequence, and (b) an alignment to the at least one reference sequence, that penalizes the score function based on a fraction of bases that are not consistent with the map;
apply the optimized score function to the combined raw intensity output;
map, by the nucleic acid sequencer, the deconvolved single raw intensity output to a set of possible alignments, the at least one reference sequence and a set of possible base calls; and
select the first base call and the second base call from the set of possible base calls.
19. A method for selecting a first base call for a first nucleic acid sequence and a second base call for a second nucleic acid sequence, comprising:
obtaining, via a nucleic acid sequencer with a single flow cell, a single raw intensity output derived from the first nucleic acid sequence and the second nucleic acid sequence, wherein the first nucleic acid sequence is in a same flow cell as the second nucleic acid sequence;
obtaining at least one reference sequence for a genome related to the first and second nucleic acid sequence;
with the nucleic acid sequencer, generating and optimizing a score function, for a plurality of intensities where the score function is configured to determine: (a) that a map path is consistent with one of the first nucleic acid sequence or the second nucleic acid sequence in view of a noise component from another one of the first nucleic acid sequence or the second nucleic acid sequence, and (b) an alignment to the at least one reference sequence, that penalizes the score function based on a fraction of bases that are not consistent with the map path;
applying the optimized score function to the combined raw intensity output to deconvolve the single raw intensity output into a plurality of constituent intensities;
mapping, by the nucleic acid sequencer, the deconvolved single raw intensity output to a set of possible alignments, the at least one reference sequence and a set of possible base calls; and
selecting, by the nucleic acid sequencer, the first base call and the second base call from the set of possible base calls.
20. The method of claim 19 , further comprising:
placing the first and second nucleic acid molecules on a single flow cell; and
inserting the single flow cell into the nucleic acid sequencer with a single flow cell.