BAMBAM: parallel comparative analysis of high-throughput sequencing data
A differential sequence object is constructed on the basis of alignment of sub-strings via incremental synchronization of sequence strings using known positions of the sub-strings relative to a reference genome sequence. An output file is then generated that comprises only relevant changes with respect to the reference genome.
1. A genomic difference sequence analysis system comprising:
a first database storing a tumor genomic sequence of a patient's tumor and a germline genomic sequence of a patient's healthy tissue; and
a sequence analysis engine coupled with the first database and comprising a tangible, non-transitory computer readable memory storing software instructions and at least one processor, wherein the processor is configured to execute the software instructions to:
obtain, from the first database, a tumor string representing at least portion of the tumor genomic sequence and a germline string representing at least a portion of the germline genomic sequence, the tumor string comprising tumor reads and the germline string comprising germline reads, wherein the tumor reads of the tumor string and the germline reads of the germline string overlap with each other at at least one genomic position;
generate a local differential string that includes one or more differences between the tumor string and the germline string at the at least one genomic position by (1) incrementally obtaining a set of tumor reads and a set of germline reads that overlap a respective genomic position of the at least one genomic position and (2) comparing the set of tumor reads to the set of germline reads to determine the one or more differences at the respective genomic position, wherein the local differential string comprises a protein-coding sequence; and
update, in a second database, a differential sequence object based on the local differential string such that the differential sequence object includes one or more of the differences, where the differential sequence object is associated with the patient in the second database.
2. The system of claim 1 , wherein the tumor genomic sequence comprises a polypeptide-coding sequence of the patient's tumor.
3. The system of claim 2 , wherein the tumor string comprises the polypeptide-coding sequence of the patient's tumor.
4. The system of claim 1 , wherein the germline genomic sequence comprises a polypeptide-coding sequence of the patient's healthy tissue.
5. The system of claim 4 , wherein the germline string comprises the polypeptide-coding sequence of the patient's tumor.
6. The system of claim 1 , wherein the differential sequence object comprises a polypeptide-coding sequence.
7. The system of claim 1 , wherein the differential sequence object comprises differences in protein-coding sequences between the tumor string and the germline string.
8. The system of claim 1 , wherein the differential sequence object comprises differences in nucleic acid sequences between the tumor string and the germline string.
9. The system of claim 8 , wherein the differences in nucleic acid sequences comprise DNA sequence or RNA sequences.
10. The system of claim 1 , wherein the processor of the sequence analysis engine is further configured to execute the software instructions to calculate a patient-specific deviation by comparing the differential sequence object to a reference genome.
11. The system of claim 10 , where the patient-specific deviation is part of a patient-specific deviation profile.
12. The system of claim 1 , wherein the first database is configured to store the tumor genomic sequence and the germline genomic sequence according to at least one of the following file formats: BAM format and SAM format.
13. The system of claim 1 , wherein the differential sequence object comprises a constellation of local differential strings between the tumor genomic sequence and the germline genomic sequence.
14. The system of claim 1 , wherein the at least one genomic position is with respect to a known reference genomic sequence.
15. The system of claim 1 , wherein the processor of the sequence analysis engine is further configured to execute the software instructions to generate a patient-specific instruction based on a comparison of the differential sequence object to other differential sequence objects of subjects that have known conditions.
16. The system of claim 15 , wherein the patient-specific instruction includes at least one of the following: a treatment, a diagnosis, a prognosis, a predicted treatment outcome, a risk assessment, and a prescription.
17. The system of claim 1 , wherein the processor of the sequence analysis engine is further configured to execute the software instructions to update the differential sequence object based on the local differential string according to at least one of the following management functionalities:
create the differential sequence object and include the local differential string;
modify the differential sequence object according to the local differential string;
append the differential sequence object with the local differential string;
delete at least a portion of the differential sequence object;
copy at least a portion of the differential sequence object; and
analyze at least a portion of the differential sequence object.
18. The system of claim 1 , wherein the differential sequence object comprises metadata attributes.
19. The system of claim 18 , wherein the metadata attributes at least one of the following: a time stamp, a sample time, a patient name, a tissue type, a tissue state, a neoplastic growth, a ploidy, a gene copy number, a repeat copy number, an inversion, a deletion, an insertion, a viral insertion, a somatic mutation, a germline mutation, a rearrangement, a transposition, and a loss of heterozygosity.
20. The system of claim 18 , wherein differential sequence object is queryable via the metadata attributes.
21. The system of claim 1 , the germline genomic sequence of the patient's healthy tissue is obtained from a blood sample from the patient.