Bambam: parallel comparative analysis of high-throughput sequencing data

• US 9,824,181 B2
• Filed: 05/27/2016
• Issued: 11/21/2017
• Est. Priority Date: 05/25/2010
• Status: Active Grant

• Current Assignee: Regents of the University of California (University of California)
• Original Assignee: Regents of the University of California (University of California)
• Inventors: Sanborn, John Zachary, Haussler, David
• Primary Examiner(s): Smith, Paulinho E
• Application Number:
  US15/167,512
• Publication Number:
  US 20160275256A1
• Time in Patent Office: 543 Days
• Field of Search: None
• CPC Class Codes: C12Q 1/6886 for cancer immunoassay for cancer G01N33/574
• CPC Class Codes: C12Q 2600/106 Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
• CPC Class Codes: C12Q 2600/118 Prognosis of disease development
• CPC Class Codes: C12Q 2600/156 Polymorphic or mutational markers
• CPC Class Codes: G06F 2203/04806 Zoom, i.e. interaction techniques or interactors for controlling the zooming operation
• CPC Class Codes: G06F 3/04845 for image manipulation, e.g. dragging, rotation, expansion or change of colour
• CPC Class Codes: G06F 40/169 Annotation, e.g. comment data or footnotes
• CPC Class Codes: G06N 7/01 Probabilistic graphical models, e.g. probabilistic networks
• CPC Class Codes: G06T 11/206 Drawing of charts or graphs
• CPC Class Codes: G16B 20/20 Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
• 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 40/00 ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
• CPC Class Codes: G16H 10/40 for data related to laboratory analysis, e.g. patient specimen analysis
• CPC Class Codes: G16H 10/60 for patient-specific data, e.g. for electronic patient records
• CPC Class Codes: G16H 50/20 for computer-aided diagnosis, e.g. based on medical expert systems
• CPC Class Codes: G16H 50/30 for calculating health indices; for individual health risk assessment
• CPC Class Codes: G16H 70/20 relating to practices or guidelines
• CPC Class Codes: Y02A 90/10 Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
• CPC Class Codes: Y02A 90/30 Assessment of water resources

Abstract

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.

Claims

1. A parallel genomic comparative analysis system comprising:

a memory; and

an sequence analysis engine coupled with the memory and configured to;

identify a genomic position within a reference genome;

access a first file storing tumor sequence data including short reads associated with a tumor tissue;

access a second file storing match normal sequence data short reads associated with a matched normal tissue;

store in the memory a tumor dataset having tumor short read sequences from the first file where the tumor short read sequences overlap the genomic position;

store in the memory a matched normal dataset having matched normal short read sequences from the second file and that overlap the genomic position;

select a tumor genotype and a matched normal genotype that maximize a joint probability as a function of the tumor short read sequences and the match normal short read sequences at the genomic position, wherein the joint probability depends on one of a probability calculated as a multinomial operating as a function of the matched normal genotype or as a probability calculated as a multinomial operating as a function of the tumor genotype; and

store a difference between the tumor genotype and the matched normal genotype in a device memory.

2. The system of claim 1, wherein the memory is configured to store the tumor dataset and the matched normal dataset simultaneously.

3. The system of claim 1, wherein the sequence analysis engine is further configured to synchronize the first and the second file.

4. The system of claim 1, wherein the sequence analysis engine is further configured to read the first and the second files at the same time to access the files.

5. The system of claim 1, wherein the genomic position comprises a common genomic location between the first and the second file relative to the reference genome.

6. The system of claim 1, wherein the difference is selected from the group consisting of:

a somatic variant, a germline variant, a single nucleotide polymorphism, an allele-specific copy number, a loss of heterozygosity, a structural rearrangement, a chromosomal fusion, and a breakpoint.

7. The system of claim 1, wherein the analysis engine is further configured to calculate a confidence score of the tumor genotype and the match normal genotype pair.

8. The system of claim 7, wherein the analysis engine is configured to calculate the confidence scores as a posterior probability.

9. The system of claim 7, wherein the analysis engine is further configured to store the confidence score in the device memory with the difference.

10. The system of claim 1, wherein at least one of the first file and the second file comprises at least one of a BAM file and a SAM file.

11. The system of claim 1, wherein the tumor dataset comprises all tumor short read sequences in the first file that overlap the genomic position.

12. The system of claim 1, wherein the matched normal dataset comprises all matched normal short read sequences in the second file that overlap the genomic position.

13. The system of claim 1, wherein the tumor sequence data of the tumor tissue and the matched normal sequence data of the matched normal tissue originate from the same person.

14. A parallel genomic comparative analysis

system comprising:

a memory; and

a sequence analysis engine coupled with the memory and configured to;

access a first file storing tumor sequence data including short reads associated with a tumor tissue;

access a second file storing matched normal sequence data short reads associated with a matched normal tissue;

align, relative to a first genomic position within a reference genome, the short reads associated with the tumor tissue with the short reads associated with the matched normal tissue;

process at the same time all aligned short reads to determine a difference between the tumor sequence data and the matched normal sequence data;

store a difference between the tumor sequence data and the matched normal sequence data in a device memory;

align, relative to a second genomic position within the reference genome, the short reads associated with the tumor tissue with the short reads associated with the matched normal tissue;

process at the same time all aligned short reads to determine a second difference between the tumor sequence data and the matched normal sequence data; and

store the second difference between the tumor sequence data and the matched normal sequence data in the device memory.

15. The system of claim 14 wherein the sequence analysis engine is further configured to select a tumor genotype and a matched normal genotype that maximize a joint probability as a function of the tumor short read sequences and the match normal short read sequences at the genomic position.

16. The system of claim 15 wherein maximizing the joint probability is defined by

P(D_r,Ac,r)−

−

P(IVG,)P(G,Ir′

GA.;

),P 

 

(1)

D₄ , , , G;

_LIFG7Ar)*(DG₄,G₄ , , , LO(9₁₀ . . . ⁻

i6;

₄₈) 

 

(2)

where r is an observed reference allele, a is a fraction of normal contamination, where tumor and germline genotypes are defined by Gt=(ti, t₂) and Gg=(gi,g2), where t₁, t2, gi, g2 c {A, T, C, G}, and where tumor and germline sequence data are defined as a set of reads DM d_i¹, d_i² . . . , d_iⁿ} and D_g={d_gⁱ, d_g² . . . , d_gⁿ}, respectively, with the observed bases d_i′

, d_g′

E {A, T, C, G}.

17. The system of claim 15 wherein a probability of germline alleles given a germline genotype is modeled as a multinomial over four nucleotides:

where n is the total number of germline reads at the genomic position and n_A, n_G, nc, n_t are reads supporting each observed allele, and a probability of tumor alleles given a tumor genotype is modeled as a multinomial over four nucleotides;

18. A parallel genomic comparative analysis system comprising:

a memory; and

an sequence analysis engine coupled with the memory and configured to;

identify a genomic position within a reference genome;

access a first file storing tumor sequence data including short reads associated with a tumor tissue;

access a second file storing match normal sequence data short reads associated with a matched normal tissue;

store in the memory a tumor dataset having tumor short read sequences from the first file where the tumor short read sequences overlap the genomic position, wherein the tumor dataset comprises all tumor short read sequences in the first file that overlap the genomic position;

store in the memory a matched normal dataset having matched normal short read sequences from the second file and that overlap the genomic position;

select a tumor genotype and a matched normal genotype that maximize a joint probability as a function of the tumor short read sequences and the match normal short read sequences at the genomic position; and

store a difference between the tumor genotype and the matched normal genotype in a device memory.

19. The system of claim 18, wherein the matched normal dataset comprises all matched normal short read sequences in the second file that overlap the genomic position.