Detecting mutations and ploidy in chromosomal segments

• US 10,179,937 B2
• Filed: 04/21/2015
• Issued: 01/15/2019
• Est. Priority Date: 04/21/2014
• Status: Active Grant

• Current Assignee: Natera Incorporated (Natera)
• Original Assignee: Natera Incorporated (Natera)
• Inventors: Babiarz, Joshua, Constantin, Tudor Pompiliu, Eubank, Lane A., Gemelos, George, Hill, Matthew Micah, Kirkizlar, Huseyin Eser, Rabinowitz, Matthew, Sakarya, Onur, Sigurjonsson, Styrmir, Zimmermann, Bernhard
• Primary Examiner(s): Riggs, II, Larry D
• Application Number:
  US14/692,703
• Publication Number:
  US 20170107576A1
• Time in Patent Office: 1,365 Days
• Field of Search: None
• CPC Class Codes: C12Q 1/6869 Methods for sequencing
• CPC Class Codes: C12Q 1/6886 for cancer immunoassay for cancer G01N33/574
• CPC Class Codes: C12Q 2539/10 The purpose being sequence identification by analysis of gene expression or genome comparison characterised by
• CPC Class Codes: C12Q 2600/156 Polymorphic or mutational markers
• CPC Class Codes: C12Q 2600/158 Expression markers
• CPC Class Codes: C12Q 2600/16 Primer sets for multiplex assays
• CPC Class Codes: C12Q 2600/172 Haplotypes
• CPC Class Codes: G06N 20/00 Machine learning
• CPC Class Codes: G06N 7/01 Probabilistic graphical models, e.g. probabilistic networks
• CPC Class Codes: G16B 15/00 ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
• CPC Class Codes: G16B 20/00 ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
• CPC Class Codes: G16B 20/10 Ploidy or copy number detection
• CPC Class Codes: G16B 20/20 Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
• CPC Class Codes: G16B 25/00 ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
• CPC Class Codes: G16B 25/20 Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
• 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: G16B 40/20 Supervised data analysis
• CPC Class Codes: G16H 10/40 for data related to laboratory analysis, e.g. patient specimen analysis
• CPC Class Codes: G16H 50/20 for computer-aided diagnosis, e.g. based on medical expert systems
• CPC Class Codes: G16Z 99/00 Subject matter not provided for in other main groups of this subclass

Abstract

The invention provides methods, systems, and computer readable medium for detecting ploidy of chromosome segments or entire chromosomes, for detecting single nucleotide variants and for detecting both ploidy of chromosome segments and single nucleotide variants. In some aspects, the invention provides methods, systems, and computer readable medium for detecting cancer or a chromosomal abnormality in a gestating fetus.

Claims

1. A method for determining ploidy of a chromosomal segment in a sample of an individual, the method comprising:

isolating circulating tumor DNA from a blood sample, serum sample, or plasma sample from an individual suspected of having cancer;

amplifying at least 1,000 polymorphic loci relating to cancer-associated copy number variation from the circulating tumor DNA to obtain amplicons;

sequencing the amplicons to obtain allele frequency data generated from the sample, comprising the amount of each allele present in the sample at each loci in a set of polymorphic loci on the chromosomal segment;

generating phased allelic information for the set of polymorphic loci by estimating the phase of the allele frequency data;

generating individual probabilities of allele frequencies for the polymorphic loci for different ploidy states using the allele frequency data, wherein the individual probabilities are generated using a set of models of both different ploidy states and allelic imbalance fractions for the set of polymorphic loci;

generating joint probabilities for the set of polymorphic loci using the individual probabilities and the phased allelic information and considering the linkage between polymorphic loci on the chromosome segment; and

selecting, based on the joint probabilities, a best fit model indicative of chromosomal ploidy, thereby determining ploidy of the chromosomal segment, and the method further comprises determining based on the selecting, whether copy number variation is present in cells of a tumor of the individual.

2. The method of claim 1, wherein the allele frequency data is generated from nucleic acid sequence data.

3. The method of claim 1, further comprising correcting for errors in the allele frequency data and using the corrected allele frequency data for the generating individual probabilities step, wherein the errors that are corrected are allele amplification efficiency bias errors.

18. The method of claim 1, further comprising receiving haplotype information of the chromosome segment for a tumor of the individual and using the haplotype information to generate the set of models of different ploidy states and allelic imbalance fractions of the set of polymorphic loci.

19. The method of claim 1, wherein an average allelic imbalance of greater than 0.45% indicates that copy number variation is present in cells of a tumor of the individual.

21. The method of claim 1, further comprising determining whether a single nucleotide variant is present in the sample in a set of single nucleotide variants, wherein detecting either the copy number variation or the single nucleotide variant or both indicates the presence of circulating tumor nucleic acids in the sample.

22. The method of claim 1, wherein individual probabilities are generated by modeling ploidy states of a first homolog of the chromosome segment and a second homolog of the chromosome segment, wherein the ploidy states are:

(1) all cells have no deletion or amplification of the first homolog or the second homolog of the chromosome segment;

(2) at least some cells have a deletion of the first homolog or an amplification of the second homolog of the chromosome segment; and

(3) at least some cells have a deletion of the second homolog or an amplification of the first homolog of the chromosome segment.

4. A method for detecting chromosomal ploidy in a sample of an individual, the method comprising:

isolating circulating tumor DNA from a blood sample, serum sample, or plasma sample from an individual suspected of having cancer;

amplifying at least 1,000 polymorphic loci relating to cancer-associated copy number variation from the circulating tumor DNA to obtain amplicons;

sequencing the amplicons to obtain nucleic acid sequence data generated from the sample, for alleles at a set of polymorphic loci on a chromosome segment in the individual;

detecting allele frequencies at the set of loci using the nucleic acid sequence data;

correcting for allele amplification efficiency bias in the detected allele frequencies to generate corrected allele frequencies for the set of polymorphic loci;

generating phased allelic information for the set of polymorphic loci by estimating the phase of the nucleic acid sequence data;

generating individual probabilities of allele frequencies for the polymorphic loci for different ploidy states by comparing the corrected allele frequencies to a set of models of different ploidy states and allelic imbalance fractions of the set of polymorphic loci;

generating joint probabilities for the set of polymorphic loci by combining the individual probabilities considering the linkage between polymorphic loci on the chromosome segment; and

selecting, based on the joint probabilities, the best fit model indicative of chromosomal ploidy, and the method further comprises determining based on the selecting, whether copy number variation is present in cells of a tumor of the individual.

5. A method according to claim 4, wherein individual probabilities are generated by modeling ploidy states of a first homolog of the chromosome segment and a second homolog of the chromosome segment.

6. A method according to claim 4, wherein the errors that are corrected comprise ambient contamination and genotype contamination and are determined for homozygous alleles in the sample.

7. A method according to claim 4, wherein the selecting is performed by analyzing a magnitude of a difference between the phased allelic information and estimated allelic frequencies generated for the models, and the model with the smallest magnitude of difference is selected.

8. A method according to claim 4, wherein the individual probabilities of allele frequencies are generated based on a beta binomial model of expected and observed allele frequencies at the set of polymorphic loci.

9. The method of claim 8, where the individual probabilities are generated using a Bayesian classifier.

10. The method of claim 4, wherein the nucleic acid sequence data is generated by performing high throughput DNA sequencing of a plurality of copies of a series of amplicons generated using a multiplex amplification reaction, wherein each amplicon of the series of amplicons spans at least one polymorphic loci of the set of polymorphic loci and wherein each of the polymorphic loci of the set is amplified.

11. The method of claim 10, wherein the multiplex amplification reaction is performed under limiting primer conditions.

12. The method of claim 10, wherein chromosomal ploidy is determined for a set of chromosome segments in the sample by performing the method on each chromosome segment of the set.

13. The method of claim 12, wherein the multiplex amplification reaction comprises between 1,000 and 100,000 multiplex reactions.

14. The method of claim 4, wherein an average allelic imbalance of greater than 0.45% indicates that copy number variation is present in cells of a tumor of the individual.

15. The method of claim 4, wherein the sample is a plasma sample from an individual previously known to have a solid tumor.

16. The method of claim 4, further comprising determining whether a single nucleotide variant from a set of single nucleotide variants is present in the sample, wherein detecting either the copy number variation or the single nucleotide variant or both, indicates the presence of circulating tumor nucleic acids in the sample.

17. The method of claim 15, further comprising receiving haplotype information of the chromosome segment for a tumor of the individual and using the haplotype information to generate the set of models of different ploidy states and allelic imbalance fractions of the set of polymorphic loci.

20. A method according to claim 5, wherein the ploidy states are:

(1) all cells have no deletion or amplification of the first homolog or the second homolog of the chromosome segment;

(2) at least some cells have a deletion of the first homolog or an amplification of the second homolog of the chromosome segment; and

(3) at least some cells have a deletion of the second homolog or an amplification of the first homolog of the chromosome segment.