BAMBAM: parallel comparative analysis of high-throughput sequencing data
The present invention relates to methods for evaluating and/or predicting the outcome of a clinical condition, such as cancer, metastasis, AIDS, autism, Alzheimer's, and/or Parkinson's disorder. The methods can also be used to monitor and track changes in a patient's DNA and/or RNA during and following a clinical treatment regime. The methods may also be used to evaluate protein and/or metabolite levels that correlate with such clinical conditions. The methods are also of use to ascertain the probability outcome for a patient's particular prognosis.
1. A computer-based sequence analysis system comprising:
a computer readable memory configured to store at least a first and second digital sequence datasets associated with a healthy tissue and a suspected diseased tissue, respectively; and
a sequence analysis engine having at least one processor coupled with the computer readable memory and configured, upon execution of software instructions, to:
determine a common sequence location within the first and second digital sequence datasets;
generate at least one pair of sequence pileups from the first and second digital sequence datasets and that overlap at the common sequence location;
infer genotypes for the healthy tissue and suspected diseased tissue at the common sequence location as a function of the at least one pair of sequence pileups;
identify a sequence difference at the common sequence location and between genotypes;
filter false positives based on reads in the at least one pair of sequence pileups; and
store the sequence difference as a differential sequence object in a device memory.
2. The system of claim 1 , wherein the first and second digital sequences are derived from tissues of the same person.
3. The system of claim 1 , wherein at least one of the first and second digital sequences is derived from blood.
4. The system of claim 1 , wherein the second digital sequence is derived from a tumor tissue.
5. The system of claim 1 , wherein at least one of the first and second digital sequence datasets comprises DNA sequences.
6. The system of claim 1 , wherein at least one of the first and second digital sequence datasets comprises RNA sequences.
7. The system of claim 6 , wherein the RNA sequences represent at least one of the following types of RNA: mRNA, rRNA, tRNA, miRNA, and asRNA.
8. The system of claim 1 , wherein at least one of the first and second digital sequence datasets comprises genomic sequences.
9. The system of claim 1 , wherein at least one of the first and second digital sequence datasets comprises a proteomic sequence.
10. The system of claim 9 , wherein the proteomic sequence comprises a polypetide sequence.
11. The system of claim 1 , wherein the first and the second digital sequence datasets comprise short reads.
12. The system of claim 1 , wherein the at least one pair of sequence pileups includes sequence reads having up to 100 positions.
13. The system of claim 1 , wherein the at least one pair of sequence pileups includes more than 30 reads.
14. The system of claim 1 , wherein one of the at least one of first and second digital sequence datasets includes at least a billion reads.
15. The system of claim 1 , wherein the false positives are filtered based on read information of reads within the first and second digital sequence datasets.
16. The system of claim 15 , wherein the read information includes at least one of the following: a strand to which a read maps, an allele position in the reads, and an average quality of alleles.
17. The system of claim 1 , wherein the false positives are filtered based on allele distribution.
18. The system of claim 17 , wherein the allele distribution indicates variant alleles are found near a tail of a read.
19. The system of claim 17 , wherein the false positives are filtered when the allele distribution is skewed from an expected distribution.
20. The system of claim 19 , wherein the expected distribution is a random distribution.
21. The system of claim 1 , wherein the differential sequence object includes at least one variant call.