SYSTEMS AND METHODS FOR HYBRID ASSEMBLY OF NUCLEIC ACID SEQUENCES

• US 20120330559A1
• Filed: 06/20/2012
• Published: 12/27/2012
• Est. Priority Date: 06/21/2011
• Status: Abandoned Application

• Current Assignee: Life Technologies Corporation (Thermo Fisher Scientific Incorporated)
• Original Assignee: Life Technologies Corporation (Thermo Fisher Scientific Incorporated)
• Inventors: JIANG, Hongshan, XU, Zhao, INGMAN, Max
• Application Number:
  US13/528,470
• Publication Number:
  US 20120330559A1
• Time in Patent Office: Days
• US Class Current: 702/19
• 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
• CPC Class Codes: G16B 30/20 Sequence assembly

Abstract

Systems and methods for assembling a nucleic acid sequence are disclosed. A plurality of single fragment sequence reads and a plurality of paired fragment sequence reads are received. Each paired fragment sequence read comprises at least two sequence reads separated by an insert. Single fragment sequence reads are assembled into a plurality of contigs, and the paired fragment sequence reads are mapped to the contigs. Further, gap regions comprising a portion of the partially assembled nucleic acid sequence for which the single fragment sequence reads do not map are identified, and hanging pairwise sequence reads of the mapped paired fragment sequence reads are used to fill in the gap region.

Claims

1. A computer implemented method for assembling a nucleic acid sequence, comprising:

receiving, into a memory, a plurality of single fragment sequence reads and a plurality of paired fragment sequence reads, each paired fragment sequence read comprising at least two sequence reads separated by an insert;

assembling the single fragment sequence reads into a plurality of contigs;

mapping the paired fragment sequence reads to the contigs;

identifying a gap region comprising a portion of the partially assembled nucleic acid sequence for which the single fragment sequence reads do not map, and

utilizing hanging pairwise sequence reads of the mapped paired fragment sequence reads to fill in the gap region using a processor.

2. The computer implemented method of claim 1, further comprising estimating a size of the gap region based on a size of the contigs, an insert size distribution, and mapped locations of paired fragment sequence reads spanning the gap region.

3. The computer implemented method of claim 2, further comprising determining the insert size distribution from paired fragment sequence reads having both sequence reads mapped to a same contig.

4. The computer implemented method of claim 1, further comprising identifying first and second contigs as adjacent when a first sequence read of a paired fragment sequence read is mapped to the first contig and a second sequence read of a paired fragment sequence read is mapped to the second contig.

5. The computer implemented method of claim 1, further comprising classifying contigs of the plurality of contigs as unique contigs or repeat contigs.

6. The computer implemented method of claim 1, further comprising determining a mismatch error rate.

7. The computer implemented method of claim 1, further comprising determining a rearrangement error frequency.

8. The computer implemented method of claim 1, further comprising using a directed node graph to represent the relationships between the plurality of contigs.

9. A system for assembling a nucleic acid sequence, comprising:

a computing device, including;

a contig assembly engine configured to assemble single fragment sequence reads into one or more contigs;

a mapping engine configured to map a plurality of paired fragment sequence reads to the assembled contigs, each paired fragment sequence read comprising at least two sequence reads separated by an insert;

a scaffolding engine configured to form a sequence scaffold from the mapped paired fragment sequence reads and contigs; and

a gap-filling engine configured to utilize hanging pairwise sequences of the mapped paired fragment sequence reads to fill in gap regions in the sequence scaffold.

10. The system of claim 9, wherein the single fragment sequence reads have a length of greater than about 100 bases.

11. The system of claim 9, wherein the scaffolding engine is further configured to estimate the size of the gap region based on a size of the contigs, an insert size distribution, and mapped locations of paired fragment sequence reads spanning the gap region.

12. The system of claim 11, wherein the scaffolding engine is further configured to determine the insert size distribution from paired fragment sequence reads having both sequence reads mapped to a same contig.

13. The system of claim 9, wherein the scaffolding engine is further configured to identify first and second contigs as adjacent when a first sequence read of a paired fragment sequence read is mapped to the first contig and a second sequence read of a paired fragment sequence read is mapped to the second contig.

14. The system of claim 9, wherein the contig assembly engine is further configured to classify contigs of the plurality of contigs as unique contigs or repeat contigs.

15. The system of claim 9, wherein the scaffolding engine is further configured to use a directed node graph to represent the relationships between the plurality of contigs.

16. A non-transitory computer readable media having a computer readable program code embodied therein, the computer readable program code adapted to be executed by a processor to implement a method for annotating called variants in a sample genome, comprising:

receiving a plurality of single fragment sequence reads and a plurality of paired fragment sequence reads, each paired fragment sequence read comprising at least two sequence reads separated by an insert;

assembling the single fragment sequence reads into a plurality of contigs;

mapping the paired fragment sequence reads to the contigs;

identifying a gap region comprising a portion of the partially assembled nucleic acid sequence for which the single fragment sequence reads do not map; and

utilizing hanging pairwise sequence of the mapped paired fragment sequence reads to fill in the gap region.

17. The non-transitory computer readable media of claim 16, further comprising estimating a size of the gap region based on a size of the contigs, an insert size distribution, and mapped locations of paired fragment sequence reads spanning the gap region.

18. The non-transitory computer readable media of claim 17, further comprising determining the insert size distribution from paired fragment sequence reads having both sequence reads mapped to a same contig.

19. The non-transitory computer readable media of claim 16, further comprising using a directed node graph to represent the relationships between the plurality of contigs.

20. The non-transitory computer readable media of claim 16, further comprising identifying first and second contigs as adjacent when a first sequence read of a paired fragment sequence read is mapped to the first contig and a second sequence read of a paired fragment sequence read is mapped to the second contig.