Methods for accurate sequence data and modified base position determination

• US 9,747,414 B2
• Filed: 05/23/2014
• Issued: 08/29/2017
• Est. Priority Date: 11/07/2008
• Status: Active Grant

• Current Assignee: Industrial Technology Research Institute
• Original Assignee: Industrial Technology Research Institute
• Inventors: Pan, Chao-Chi, Fann, Jenn-Yeh, Chiou, Chung-Fan, Chien, Hung-Chi, Chen, Hui-Ling
• Primary Examiner(s): Martinell, James
• Application Number:
  US14/286,904
• Publication Number:
  US 20140256570A1
• Time in Patent Office: 1,194 Days
• Field of Search: None
• CPC Class Codes: C12Q 1/6869 Methods for sequencing
• CPC Class Codes: C12Q 2531/125 Rolling circle
• CPC Class Codes: C12Q 2537/165 Mathematical modelling, e.g. logarithm, ratio
• CPC Class Codes: G16B 30/00 ICT specially adapted for sequence analysis involving nucleotides or amino acids
• CPC Class Codes: G16B 30/10 Sequence alignment; Homology search

Abstract

Disclosed herein are methods of determining the sequence and/or positions of modified bases in a nucleic acid sample present in a circular molecule with a nucleic acid insert of known sequence comprising obtaining sequence data of at least two insert-sample units. In some embodiments, the methods comprise obtaining sequence data using circular pair-locked molecules. In some embodiments, the methods comprise calculating scores of sequences of the nucleic acid inserts by comparing the sequences to the known sequence of the nucleic acid insert, and accepting or rejecting repeats of the sequence of the nucleic acid sample according to the scores of one or both of the sequences of the inserts immediately upstream or downstream of the repeats of the sequence of the nucleic acid sample.

Claims

1. A method of determining the sequence of a nucleic acid sample comprising:

a. generating a nucleic acid molecule from a linear pair-locked molecule to obtain sequence data from the linear pair-locked molecule,

the linear pair-locked molecule and the nucleic acid molecule each comprising at least two insert-sample units each comprising a nucleic acid insert and the nucleic acid sample, wherein the nucleic acid insert has a known sequence, and the nucleic acid insert is immediately upstream or downstream to the nucleic acid sample,

and the sequence data comprising sequence data of at least two insert-sample units, wherein the sequence data of at least two insert-sample units each comprise a repeat of the sequence of the nucleic acid sample;

b. calculating scores of the sequences of at least two inserts of the sequence data of step (a) by comparing the sequences to the known sequence of the insert;

c. accepting or rejecting at least two of the repeats of the sequence of the nucleic acid sample of the sequence data of step (a) according to the scores of one or both of the sequences of the inserts;

d. compiling an accepted sequence set comprising at least one repeat of the sequence of the nucleic acid sample accepted in step (c); and

e. determining the sequence of the nucleic acid sample using the accepted sequence set.

2. The method of claim 1, wherein obtaining sequence data comprises single molecule sequencing.

3. The method of claim 2, wherein the single molecule sequencing comprises sequencing by a method chosen from single molecule sequencing by synthesis and ligation sequencing.

4. The method of claim 3, wherein the single molecule sequencing comprises real-time single molecule sequencing by synthesis.

5. The method of claim 3, wherein the single molecule sequencing comprises single molecule sequencing by synthesis by a method chosen from pyrosequencing, reversible terminator sequencing, and third-generation sequencing.

6. The method of claim 3, wherein the single molecule sequencing comprises nanopore sequencing.

7. The method of claim 1, further comprising providing the linear pair-locked molecule, wherein the nucleic acid sample is ligated to the nucleic acid insert to form an insert-sample unit of the linear pair-locked molecule.

8. The method of claim 1, wherein the nucleic acid sample is obtained from an RNA sample.

9. The method of claim 1, wherein the nucleic acid sample is obtained from a genomic DNA sample.

10. The method of claim 1, wherein the nucleic acid insert comprises a promoter and synthesizing the product nucleic acid molecule comprises contacting the promoter with an RNA polymerase that recognizes the promoter followed by synthesizing a product nucleic acid molecule comprising ribonucleotide residues.

11. The method of claim 1, wherein the nucleic acid insert has a melting temperature ranging from 30°

C. to 90°

C.

12. The method of claim 1, wherein the nucleic acid insert has a length ranging from 14 to 200 nucleotide residues.

13. The method of claim 1, wherein the accepted sequence set is in a form chosen from a multi-element variable and a single data string comprising the sequence data of step (a) which has been processed to delete, overwrite, or omit repeats of the sequence of the nucleic acid sample rejected in step (d).

14. The method of claim 1, wherein the accepted sequence set is in the form of the multi-element variable of a type chosen from a list, array, hash, and matrix.

15. The method of claim 1, wherein at least two repeats of the sequence of the nucleic acid sample are accepted in step (c), and determining the sequence of the nucleic acid sample comprises determining a consensus sequence based on the at least two repeats of the sequence of the nucleic acid sample accepted in step (c).

16. The method of claim 15, wherein the consensus sequence comprises probabilistically represented bases at at least one position at which the at least two repeats of the sequence of the nucleic acid sample accepted in step (d) differ.

17. The method of claim 15, wherein at least three repeats of the sequence of the nucleic acid sample are accepted in step (d), and determining a consensus sequence comprises determining the majority vote of the at least three repeats of the sequence of the nucleic acid sample accepted in step (d).

18. The method of claim 15, wherein the consensus sequence is a position weight matrix.

19. The method of claim 15, wherein the consensus sequence is a flat sequence.

20. The method of claim 19, wherein the flat sequence comprises at least one ambiguity code.

21. The method of claim 15, wherein the consensus sequence comprises confidence levels.

22. The method of claim 21, wherein the confidence levels are expressed in a form chosen from base frequency, information content, and Phred quality score.

23. The method of claim 21, wherein steps (b)-(f) of claim 1 are performed in real time, and the consensus sequence and confidence levels are updated in real time.

24. The method of claim 23, wherein the method is performed until a set minimum level of confidence at a preselected percentage of positions of the consensus sequence is achieved.

25. The method of claim 24, further comprising generating an alert when a preselected percentage of positions reach the set minimum level of confidence.

26. The method of claim 24, wherein the set minimum level of confidence is chosen from 90%, 95%, 99%, 99.5%, 99.9%, 99.95%, and 99.99% base call accuracy.

27. The method of claim 1, further comprising repeating the steps of claim 1 with at least one other nucleic acid sample from the same source, species, or strain as the nucleic acid sample of claim 1 that has a sequence that partially overlaps the sequence of the nucleic acid sample of claim 1, thereby determining at least one other sequence, and assembling the at least one other sequence with the sequence of step (f) to form a contig.

28. The method of claim 1, wherein the scores of step (c) are used to estimate a confidence level of the sequence data of step (b) as a whole.

29. The method of claim 1, wherein calculating scores comprises determining the number of mismatches between the at least two inserts of the sequence data and the known sequence of the insert.

30. The method of claim 1, wherein calculating scores comprises determining the percent identity of the at least two inserts of the sequence data to the known sequence of the insert.

31. The method of claim 1, wherein calculating scores comprises performing an alignment between the at least two inserts of the sequence data and the known sequence of the insert.

32. The method of claim 31, wherein performing an alignment comprises using an algorithm chosen from BLAST, MEGABLAST, Smith-Waterman alignment, and Needleman-Wunsch alignment.

33. The method of claim 1, wherein the scores are generated on a basis chosen from a count basis and a rate basis.

34. The method of claim 1, wherein accepting or rejecting at least two of the repeats of the sequence of the nucleic acid sample of the sequence data of step (b) comprises accepting those of the at least two repeats of the sequence of the nucleic acid sample that are immediately upstream or downstream of a sample insert sequence with a score greater than or equal to a pre-determined threshold, and rejecting those that are not.