Fragmentation-based methods and systems for de novo sequencing

• US 20050009053A1
• Filed: 04/22/2004
• Published: 01/13/2005
• Est. Priority Date: 04/25/2003
• Status: Abandoned Application

• Current Assignee: Sequenom Incorporated (Laboratory Corporation Of America Holdings)
• Original Assignee: Sequenom Incorporated (Laboratory Corporation Of America Holdings)
• Inventors: Boom, Dirk van den, Boecker, Sebastian
• Application Number:
  US10/830,943
• Publication Number:
  US 20050009053A1
• Time in Patent Office: Days
• US Class Current: 435/6
• CPC Class Codes: C12Q 1/6872 involving mass spectrometry
• CPC Class Codes: C12Q 2521/531 Glycosylase
• CPC Class Codes: C12Q 2523/107 Chemical cleaving agents
• CPC Class Codes: C12Q 2523/313 Irradiation, e.g. UV irradiation
• CPC Class Codes: C12Q 2525/186 incorporating a non-extendable or blocking moiety
• CPC Class Codes: G01N 33/6848 Methods of protein analysis involving mass spectrometry
• CPC Class Codes: G16B 30/00 ICT specially adapted for sequence analysis involving nucleotides or amino acids
• CPC Class Codes: Y02A 50/30 Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change
• CPC Class Codes: Y02A 90/10 Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

Methods and systems, particularly mass spectrometric methods and systems, for the analysis and sequencing of biomolecules, particularly nucleic acids, by fragmentation are provided.

Claims

1. A method of obtaining sequence information from a target biomolecule, comprising:

fragmenting the target biomolecule into a plurality of fragments by partial cleavage;

performing mass spectrometry on the plurality of fragments to produce mass spectra of the fragments;

extracting peak information from the produced mass spectra;

constructing sequencing graphs using the extracted peak information; and

traversing the sequencing graphs to reconstruct the sequence information of the target biomolecule.

2. The method of claim 1, wherein constructing sequencing graphs includes generating a plurality of graphs having vertices and edges, each sequencing graph of the plurality of graphs representing a sequencing graph with a distinct cleavage reaction different from cleavage reactions used in other sequencing graphs of the plurality of graphs.

3. The method of claim 1, wherein each fragment of the plurality of fragments comprises a compomer.

4. The method of claim 3, wherein traversing the sequencing graphs includes tracing through each sequencing graph in the plurality of graphs, starting at a source vertex.

5. The method of claim 4, wherein traversing the sequencing graphs further includes setting the source vertex as a current vertex.

6. The method of claim 5, wherein traversing the sequencing graphs further includes setting a current sequence with the compomer of the current vertex.

7. The method of claim 6, wherein traversing the sequencing graphs further includes proceeding to the current vertex of the sequencing graph of an untested cleavage reaction.

8. The method of claim 7, wherein traversing the sequencing graphs further includes moving to a connecting vertex to the current vertex through an edge.

9. The method of claim 8, wherein traversing the sequencing graph further includes processing the connecting vertex.

10. The method of claim 9, wherein traversing the sequencing graphs further includes producing a candidate sequence by combining the traversed edge and vertex to the current sequence.

11. The method of claim 10, wherein traversing the sequencing graphs further includes determining whether the current vertex is an ending vertex.

12. The method of claim 11, wherein traversing the sequencing graphs further includes determining whether a length of the reconstructed sequence has reached a predetermined threshold.

13. The method of claim 12, wherein traversing the sequencing graphs further includes outputting the current sequence as a candidate sequence if the current vertex is the ending vertex and the length of the reconstructed sequence has reached the predetermined threshold.

14. The method of claim 12, wherein traversing the sequencing graphs further includes performing recursion after edge traversion if the current vertex is not the ending vertex.

15. The method of claim 12, wherein traversing the sequencing graphs further includes performing recursion after edge traversion if the length of the reconstructed sequence has not reached the predetermined threshold.

16. The method of claim 1, wherein traversing the sequencing graphs further includes backtracking to search for unexplored branching possibilities in the plurality of graphs.

17. A method for producing a candidate sequence of a biomolecule, comprising:

receiving a plurality of sequencing graphs, each sequencing graph having a plurality of vertices and edges, where each vertex represents a compomer of the biomolecule, and each edge represents a cut base of the sequencing graph; and

generating the candidate sequence by traversing the plurality of sequencing graphs.

18. The method of claim 17, further comprising:

traversing the plurality of sequencing graphs by tracing through each sequencing graph, starting at a source vertex.

19. The method of claim 18, wherein traversing the plurality of sequencing graphs includes setting the source vertex as a current vertex.

20. The method of claim 19, wherein traversing the plurality of sequencing graphs further includes setting the candidate sequence of the biomolecule as a compomer of the current vertex.

21. The method of claim 20, wherein traversing the plurality of sequencing graphs further includes proceeding to the current vertex of the sequencing graph of an untested cut base.

22. The method of claim 21, wherein traversing the plurality of sequencing graphs further includes moving to a connecting vertex from the current vertex through an edge.

23. The method of claim 22, wherein traversing the plurality of sequencing graphs further includes resetting the candidate sequence by appending compomers of the traversed edge and the connecting vertex to the previous-candidate sequence.

24. A program product for use in a computer that executes program instructions recorded in a computer-readable media to produce a candidate sequence of a biomolecule, the program product comprising:

a recordable medium; and

a plurality of computer-readable program instructions on the recordable media that are executable by the computer to perform a method comprising;

receiving a plurality of sequencing graphs, each sequencing graph having a plurality of vertices and edges, where each vertex represents a compomer of the biomolecule, and each edge represents a cut base of the sequencing graph; and

generating the candidate sequence by traversing the plurality of sequencing graphs.

25. The program product of claim 24, further comprising:

traversing the plurality of sequencing graphs by tracing through each sequencing graph, starting at a source vertex.

26. The program product of claim 25, wherein traversing the plurality of sequencing graphs includes setting the source vertex as a current vertex.

27. The program product of claim 26, wherein traversing the plurality of sequencing graphs further includes setting the candidate sequence of the biomolecule as a compomer of the current vertex.

28. The program product of claim 27, wherein traversing the plurality of sequencing graphs further includes proceeding to the current vertex of the sequencing graph of an untested cut base.

29. The program product of claim 28, wherein traversing the plurality of sequencing graphs further includes moving to a connecting vertex from the current vertex through an edge.

30. The program product of claim 29, wherein traversing the plurality of sequencing graphs further includes the candidate sequence by appending compomers of the traversed edge and the connecting vertex to the candidate sequence.

83. A combination of the program product of claim 24 and one or more specific cleavage reagents.

84. A system, comprising a computer, the program product of claim 24, and one or more specific cleavage reagents.

31. A sequencing system for obtaining sequence information from a target biomolecule, comprising:

a biomolecule workstation configured to process the target biomolecule into a plurality fragments and to produce mass spectra; and

an analysis computer configured to construct sequencing graphs using the mass spectra of the target biomelcule.

32. The system of claim 31, wherein the biomolecule workstation includes a processing station configured to receive and prepare one or more molecular samples for analysis.

33. The system of claim 32, wherein the processing station includes a cleaving element configured to provide for cleavage reactions on the one or more molecular samples to produce partially cleaved fragments.

34. The system of claim 33, wherein the biomolecule workstation includes a mass measuring station to perform mass spectrometry on the cleaved fragments.

35. The system of claim 34, wherein the biomolecule workstation includes a robotic device configured to move the molecular sample from the processing station to the mass measuring station.

36. The system of claim 35, wherein the robotic device includes a plurality of subsystems that ensure movement between the processing station and the mass measuring station to preserve the integrity of the samples.

37. The system of claim 36, wherein the plurality of subsystems include a mechanical lifting device to pick up the sample from the processing station and move the sample to the mass measuring station.

38. The system of claim 34, wherein the mass measuring station and the analysis computer are interconnected over a network.

39. The system of claim 38, wherein the network includes a local area network (LAN).

40. The system of claim 38, wherein the network includes a wireless communication channel.

41. The system of claim 38, wherein the network includes a wide area network (WAN).

42. The system of claim 41, wherein the wide area network (WAN) is the Internet.

43. The system of claim 31, wherein the analysis computer includes a neural network element to learn an efficient way to process the cleavages to obtain the sequence information of the target biomolecule.

44. A method of obtaining sequence information from a target biomolecule, comprising:

fragmenting the target biomolecule into at least two fragments by partial cleavage at specific cleavage sites;

determining the molecular weights of the at least two fragments;

determining the possible compositions of the at least two fragments;

ordering the possible compositions of the at least two fragments according to the number of specific cleavage sites that are not cleaved in each fragment;

constructing at least one sequencing graph that is a graph theoretical representation of the ordered compositions for the at least two fragments; and

traversing the at least one sequencing graph to reconstruct one or more underlying sequence candidates of the target biomolecule.

45. The method of claim 44, further comprising scoring the one or more underlying sequence candidates and determining the rank order of fitness.

46. The method of claim 45, wherein the scoring is done by statistical analysis.

47. The method of claim 46, wherein the scoring is done by maximum likelihood statistical analysis.

48. The method of claim 44 wherein the target biomolecule is DNA, and the compositions of the at least two fragments are the base compositions.

49. The method of claim 44, wherein the target biomolecule is RNA, and the compositions of the at least two fragments are the base compositions.

50. The method of claim 44, wherein the target biomolecule is a protein, and the compositions of the at least two fragments are the amino acid compositions.

51. The method of claim 44, wherein the molecular weights of the fragments are determined by mass spectrometry.

52. The method of claim 44, wherein the sequencing graph is a subgraph of a de Bruijn graph.

53. The method of claim 44, wherein the sequencing graph is traversed in a subgraph that is a walk.

60. The method of claim 44, wherein the target biomolecule contains a sequence variation.

61. The method of claim 60, wherein the sequence variation is a mutation or a polymorphism.

62. The method of claim 61, wherein the mutation is an insertion, a deletion or a substitution.

63. The method of claim 61, wherein the polymorphism is a single nucleotide polymorphism.

64. The method of claim 44, wherein the target is a target nucleic acid molecule from an organism selected from the group consisting of eukaryotes, prokaryotes and viruses.

65. The method of claim 64, wherein the organism is a bacterium.

66. The method of claim 65, wherein the bacterium is selected from the group consisting of Helicobacter pyloris, Borelia burgdorferi, Legionella pneumophilia, Mycobacteria sp. (e.g. M. tuberculosis, M. avium, M. intracellulare, M. kansaii, M. gordonae), Staphylococcus aureus, Neisseria gonorrheae, Neisseria meningitidis, Listeria monocytogenes, Streptococcus pyogenes, Streptococcus agalactiae, Streptococcus sp., Streptococcus faecalis, Streptococcus bovis, Streptococcus pneumoniae, Campylobacter sp., Enterococcus sp., Haemophilus influenzae, Bacillus antracis, Corynebacterium diphtheriae, Corynebacterium sp., Erysipelothrix rhusiopathiae, Clostridium perfringens, Clostridium tetani, Enterobacter aerogenes, Klebsiella pneumoniae, Pasturella multocida, Bacteroides sp., Fusobacterium nucleatum, Streptobacillus moniliformis, Treponema palladium, Treponema pertenue, Leptospira and Actinomyces israelli.

67. The method of claim 44, wherein a specific cleavage reagent is an RNAse.

68. The method of claim 67, wherein a specific cleavage reagents are selected from among the RNase T₁, RNase U₂, the RNase PhyM, RNase A, chicken liver RNase (RNase CL3) and cusavitin.

69. The method of claim 44, wherein a specific cleavage reagent is a glycosylase.

70. The method of claim 44, wherein sequence variations in the target biomolecule permit genotyping a subject, forensic analysis, disease diagnosis or disease prognosis.

71. The method of claim 44, wherein the method determines epigenetic changes in a target nucleic acid molecule relative to a reference nucleic acid molecule.

54. A method of obtaining nucleic acid sequence information from a target nucleic acid molecule, comprising:

subjecting the nucleic acid molecule to partial cleavage reactions with one or more specific cleavage reagents, thereby generating two or more fragments that are specific cleavage products;

determining the molecular weights of the two or more fragments;

determining the possible base compositions of the two or more fragments;

ordering the possible base compositions of the two or more fragments according to the number of specific cleavage sites that are not cleaved in each fragment;

constructing one or more sequencing graphs that are graph theoretical representations of the ordered base compositions for the two or more fragments; and

traversing the one or more sequencing graphs to reconstruct one or more underlying sequence candidates, wherein each sequencing graph corresponds to the ordered base compositions derived from a partial cleavage reaction with one base-specific cleavage reagent.

55. The method of claim 54, wherein the one or more sequencing graphs are subgraphs of de Bruijn graphs that are traversed in a subgraph that is a walk.

56. The method of claim 54, wherein the nucleic acid molecule is subject to partial cleavage with two or more base-specific cleavage reagents and two or more sequencing graphs are constructed.

57. The method of claim 56, wherein the two or more sequencing graphs are traversed serially.

58. The method of claim 56, wherein the two or more sequencing graphs are traversed in parallel.

59. The method of claim 54, wherein the molecular weights of the two or more fragments are determined by mass spectrometry.

72. A program product for use in a computer that executes program instructions recorded in a computer-readable media to obtain sequence information in a target biomolecule, the program product comprising:

a recordable medium; and

a plurality of computer-readable program instructions on the recordable media that are executable by the computer to perform a method comprising;

a) determining mass signals of target biomolecule fragments produced from partially cleaving a target biomolecule into fragments by contacting the target biomolecule with one or more base-specific cleavage reagents;

b) determining the possible compositions of the at least two fragments;

c) ordering the possible compositions of the at least two fragments according to the number of specific cleavage sites that are not cleaved in each fragment;

d) constructing at least one sequencing graph that is a graph theoretical representation of the ordered compositions for the at least two fragments; and

e) traversing the at least one sequencing graph to reconstruct one or more underlying sequence candidates of the target biomolecule.

73. The program product of claim 72, wherein the computer executable method further comprises scoring the candidate sequences and determining a rank order of sequence fitness.

74. The program product of claim 73, wherein determining a rank order of sequence fitness further comprises subjecting each of the target biomolecule candidate sequences to one or more statistical algorithms.

75. The program product of claim 72, wherein the masses are determined by mass spectrometry.

76. The method of claim 72, wherein the target biomolecule is a nucleic acid.

77. A combination of the program product of claim 72 and one or more specific cleavage reagents.

78. A system, comprising a computer, the program product of claim 72, and one or more specific cleavage reagents.

79. The combination of claim 77, further comprising:

one or more reference nucleic acid molecules; and

/or one or more natural or modified nucleoside triphosphates.

80. A kit for determining de novo sequence information in one or more target nucleic acid molecules, comprising a combination of claim 77, and optionally instructions for determining de novo sequence information.

81. The kit of claim 80, wherein a specific cleavage reagent is an RNAse.

82. The kit of claim 81, wherein the RNAses are selected from among the RNase T₁, RNase U₂, the RNase PhyM, RNase A, chicken liver RNase (RNase CL3) and cusavitin.