Methods and systems for detecting allelic imbalance in cell-free nucleic acid samples - Patents.com

Abstract

In this specification, it is recognized that the problem of distinguishing between allele imbalance samples and samples that contain contaminants or second genomes is encountered.When assaying the cell-free nucleic acid from samples that contain contaminants or second genomes, such samples may require additional manual inspection or additional sequencing runs.As a result, if failing to distinguish between allele imbalance samples and samples that contain contaminants or second genomes, the cost and time required for reliable assay of such samples may increase significantly.

Description

[Technical field]

[0001] cross reference This application claims the benefit of U.S. Provisional Patent Application No. 62 / 726,922, filed September 4, 2018, and U.S. Provisional Patent Application No. 62 / 810,625, filed February 26, 2019, each of which is incorporated by reference in its entirety. [Background technology]

[0002] background In a subject (e.g., a patient) with cancer, allelic imbalance can be caused by loss of heterozygosity and can result in a different mutant allele fraction (MAF) distribution in the assay of a cell-free nucleic acid sample from the subject compared to a sample without allelic imbalance. For example, a sample with allelic imbalance can contain germline variants with very low MAF. Germline variants with low MAF can also be observed when the sample is contaminated, for example, during processing for sequencing, or when the sample contains a second genome (other than the subject's genome), for example, from a transplant, blood transfusion, or fetus. Summary of the Invention [Means for solving the problem]

[0003] Abstract The present disclosure recognizes the problem faced in distinguishing between allelic imbalance samples and samples with contamination or second genome.When assaying cell-free nucleic acid from samples with contamination or second genome, such samples may require additional manual inspection or additional sequencing runs.As a result, failure to distinguish between allelic imbalance samples and samples with contamination or second genome samples may significantly increase the cost and time required to reliably assay such samples.The present disclosure provides methods and systems for identifying allelic imbalance or contamination in cell-free nucleic acid samples.These methods and systems may identify allelic imbalance or contamination by obtaining and analyzing quantitative measurements of small variants and copy number variations.

[0004] In one aspect, the disclosure provides a method for detecting the presence or absence of allelic imbalance in a sample from a subject, the method comprising: (a) sequencing a plurality of cell-free nucleic acid molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying germline variants present in the sample at a mutant allele fraction (MAF) for at least a portion of the plurality of aligned sequence reads, wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants identified in (c) that are between a plurality of discrete ranges of MAF values; and (e) detecting the presence or absence of the allelic imbalance in the sample based on a predetermined criterion by filtering the set of germline variants identified in (c) based on at least the quantitative measure of (d).

[0005] In one aspect, the disclosure provides a method for detecting the presence or absence of allelic imbalance in a sample from a subject, the method comprising: (a) sequencing a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying germline variants present in the sample at a mutant allele fraction (MAF) for at least a portion of the plurality of aligned sequence reads, where each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants identified in (c) that are between a plurality of discrete ranges of MAF values; and (e) detecting the presence or absence of the allelic imbalance in the sample based on a predetermined criterion by filtering the set of germline variants identified in (c) based on at least the quantitative measure of (d).

[0006] In some embodiments, the detecting in (e) comprises detecting one or more quantitative measurements from the plurality of aligned sequence reads indicative of a copy number variation (CNV) or a diploid gene (wherein the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid gene).

[0007] In some embodiments, the method further comprises detecting the presence or absence of contamination or a second genome in said sample if the absence of said allelic imbalance is detected in said sample.

[0008] In some embodiments, the set of germline variants comprises at least about 50, at least about 100, at least about 200, at least about 500, at least about 1,000, at least about 2,000, at least about 5,000, at least about 10,000, or more than about 10,000 different germline variants. In some embodiments, the set of genetic variants comprises genetic variants selected from the group consisting of single nucleotide variants (SNVs), insertions or deletions (indels), and fusions. In some embodiments, the sample is a bodily fluid sample selected from the group consisting of blood, plasma, serum, urine, saliva, mucosal secretions, sputum, stool, and tears. In some embodiments, the subject has a disease or disorder. In some embodiments, the disease is cancer.

[0009] In some embodiments, the method further comprises amplifying the cell-free DNA molecules prior to sequencing. In some embodiments, the method further comprises selectively enriching the cell-free DNA molecules for a set of loci prior to sequencing. In some embodiments, the method further comprises attaching one or more adapters comprising a barcode to the cell-free DNA molecules prior to sequencing. In some embodiments, the one or more adapters are randomly attached to both ends of the cell-free DNA molecules. In some embodiments, the cell-free DNA molecules are uniquely barcoded. In some embodiments, the cell-free DNA molecules are non-uniquely barcoded. In some embodiments, each barcode comprises a predefined or semi-random oligonucleotide sequence that, in combination with the diversity of molecules sequenced from a selected region, allows for the identification of unique cell-free DNA molecules. In some embodiments, the plurality of genomic regions comprises genetic variants found in COSMIC, The Cancer Genome Atlas (TCGA), or Exome Aggregation Consortium (ExAC). In some cases, the genetic variants may belong to a predefined set of clinically available variants. For example, such variants can be found in various databases of variants whose presence in a sample of a subject has been shown to be associated with or indicative of a disease or disorder (e.g., cancer) in the subject. Such databases of variants can include, for example, COSMIC (Catalogue of Somatic Mutations in Cancer), TCGA (The Cancer Genome Atlas), and ExAC (Exome Aggregation Consortium). In some embodiments, the multiple genomic regions include BRCA1 gene variants (e.g., BRCA1 P209L).A predefined set of such cataloged variants can be selected for further bioinformatics analysis because such variants are relevant to medical decisions (e.g., diagnosis, prognosis, treatment selection, targeted treatment, treatment monitoring, recurrence monitoring, etc.) Such a predefined set can be determined based on annotation information from public databases and clinical literature, as well as, for example, analysis of clinical samples (e.g., clinical samples from patient cohorts with known presence or absence of a disease or disorder).

[0010] In some embodiments, the plurality of discrete ranges of MAF values ​​include a first range of about 3% to about 40% and a second range of about 60% to about 97%. In some embodiments, the quantitative measure of (d) includes a plurality of sets of the genetic variants that are between the plurality of discrete ranges of MAF values. In some embodiments, the predetermined criteria includes the quantitative measure of (d) being greater than a predetermined germline variant threshold. In some embodiments, the predetermined germline variant threshold is about 21. In some embodiments, the one or more quantitative measures indicative of the CNV or the diploid gene are selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average. In some embodiments, the one or more quantitative measures indicative of the CNV or the diploid gene include two or more quantitative measures selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average. In some embodiments, the one or more quantitative measurements indicative of the CNV or diploid genes include three or more quantitative measurements selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average. In some embodiments, the predetermined criteria include one or more criteria selected from the group consisting of the following criteria: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold, a minimum CNV level across the samples is less than a predetermined minimum CNV threshold, a diploid gene fraction is less than a predetermined diploid fraction threshold, and an absolute value of the copy number average for the same germline variant is greater than a predetermined copy number average threshold and a MAF for the same germline variant is less than about 3%.In some embodiments, the predetermined criteria include two or more criteria selected from the group consisting of: the maximum CNV level across the samples is greater than a predetermined maximum CNV threshold, the minimum CNV level across the samples is less than a predetermined minimum CNV threshold, the diploid gene fraction is less than a predetermined diploid fraction threshold, and the absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and the MAF of the same germline variant is less than about 3%. In some embodiments, the predetermined criteria include three or more criteria selected from the group consisting of: the maximum CNV level across the samples is greater than a predetermined maximum CNV threshold, the minimum CNV level across the samples is less than a predetermined minimum CNV threshold, the diploid gene fraction is less than a predetermined diploid fraction threshold, and the absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and the MAF of the same germline variant is less than about 3%. In some embodiments, the predetermined criteria include the following criteria: the maximum CNV level across the samples is greater than a predetermined maximum CNV threshold, the minimum CNV level across the samples is less than a predetermined minimum CNV threshold, the diploid gene fraction is less than a predetermined diploid fraction threshold, and the absolute value of the copy number average for the same germline variant is greater than a predetermined copy number average threshold and the MAF for the same germline variant is less than about 3%. In some embodiments, the predetermined criteria include one or more thresholds selected from the group consisting of a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. In some embodiments, the predetermined criteria include two or more thresholds selected from the group consisting of: a maximum CNV threshold of about 0.20, about 0.21, or 0.22; a minimum CNV threshold of about -0.10, about -0.11, about -0.12, about -0.13, about -0.14, or about -0.15; a diploid fraction threshold of about 0.5, about 0.6, about 0.7, about 0.8, about 0.9, about 0.10; and a copy number average threshold of about 5, about 6, about 7, about 8, about 9, about 10, or about 15.In some embodiments, the predetermined criteria include three or more thresholds selected from the group consisting of the following thresholds: a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. In some embodiments, the predetermined criteria include a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10.

[0011] In some embodiments, the method further comprises detecting the presence of the contaminant or the second genome in the sample with a positive predictive value (PPV) of at least about 50%, at least about 55%, at least about 60%, at least about 65%, at least about 70%, at least about 75%, at least about 80%, at least about 85%, at least about 90%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, or at least about 99%. In some embodiments, the method further comprises detecting the absence of the contaminant or the second genome in the sample with a negative predictive value (NPV) of at least about 50%, at least about 55%, at least about 60%, at least about 65%, at least about 70%, at least about 75%, at least about 80%, at least about 85%, at least about 90%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, or at least about 99%. In some embodiments, the PPV and / or NPV are determined based on test data from a training set of samples with known contamination / allelic imbalance status (e.g., about 10 samples, about 20 samples, about 30 samples, about 40 samples, about 50 samples, about 100 samples, about 150 samples, about 200 samples, or about 250 samples).

[0012] In some embodiments, the method further comprises detecting the presence of the contamination or the second genome in the sample with a sensitivity of at least about 50%, at least about 55%, at least about 60%, at least about 65%, at least about 70%, at least about 75%, at least about 80%, at least about 85%, at least about 90%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, or at least about 99%.

[0013] In some embodiments, the method further comprises detecting the absence of the contamination or the second genome in the sample with a specificity of at least about 50%, at least about 55%, at least about 60%, at least about 65%, at least about 70%, at least about 75%, at least about 80%, at least about 85%, at least about 90%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, or at least about 99%.

[0014] In some embodiments, the method further comprises identifying the germline variant by: (i) determining a total allele count and a variant allele count for a nucleic acid variant from the cfDNA molecule; (ii) identifying an associated variable for the nucleic acid variant from the cfDNA molecule; (iii) determining a quantitative value for the associated variable for the nucleic acid variant; (iv) generating a statistical model for a predicted germline variant allele count at the genomic locus of the nucleic acid variant; (v) generating a probability value for the nucleic acid variant based at least in part on the statistical model for a predicted germline variant allele count, the quantitative value for the associated variable for the nucleic acid variant, and at least one of the total allele count and the variant allele count for the nucleic acid variant; and (vi) classifying the nucleic acid variant as (1) being of somatic origin if the p-value for the nucleic acid variant is less than a predetermined threshold, or (2) being of germline origin if the p-value for the nucleic acid variant is equal to or greater than a predetermined threshold.

[0015] In some embodiments, the method further comprises detecting an allele-specific loss in the sample based on at least one of the set of germline variants identified as present at a given MAF in (c). In some embodiments, the allele-specific loss in the sample is detected based on the presence of at least one of the set of germline variants in the sample from the subject at a MAF of less than 50%. In some embodiments, the allele-specific loss in the sample is detected based on the presence of at least one of the set of germline variants in the sample from the subject and in each of one or more additional samples from one or more subjects at a MAF of less than 50%. In some embodiments, the at least one of the set of germline variants is found in COSMIC, (The Cancer Genome Atlas; TGCA), or ExAC (Exome Aggregation Consortium). In some embodiments, the at least one of the set of germline variants is a BRCA1 gene variant. In some embodiments, the BRCA1 gene variant is BRCA1 P209L.

[0016] In another aspect, the disclosure provides a system that, when executed by at least one electronic processor, includes at least: (a) obtaining a plurality of sequence reads corresponding to a plurality of cell-free deoxyribonucleic acid (DNA) molecules from a sample of a subject; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; and (c) identifying, for at least a portion of the plurality of aligned sequence reads, a set of germline variants in the sample by identifying germline variants present in the sample at a mutant allele fraction (MAF), wherein each germline variant in the set of germline variants is identified as a mutant allele fraction (MAF). (d) determining a quantitative measure of the set of germline variants identified in (c) that are between a plurality of discrete ranges of MAF values; and (e) detecting the presence or absence of allelic imbalance in the sample based on predetermined criteria by filtering the set of germline variants identified in (c) based on at least the quantitative measure of (d).

[0017] In some embodiments, detecting in (e) further comprises detecting one or more quantitative measurements from the plurality of aligned sequence reads indicative of copy number variation (CNV) or diploid genes, wherein the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes. In some embodiments, the system further comprises a nucleic acid sequencer operably connected to the controller, wherein the nucleic acid sequencer is configured to process the plurality of cell-free DNA molecules from the sample to generate the plurality of sequence reads.

[0018] In some embodiments, the non-transitory computer executable instructions, when executed by at least one electronic processor, further perform: generating a report that optionally includes information about the presence or absence of the allelic imbalance in the sample and / or information about the presence or absence of the contamination or second genome in the sample. In some embodiments, the non-transitory computer executable instructions, when executed by at least one electronic processor, further perform: communicating the report to a third party (e.g., the subject from whom the sample originated, or a medical professional, etc.).

[0019] In one aspect, the disclosure provides a method for detecting the presence or absence of allelic imbalance in a sample from a subject, comprising: (a) accessing, by a computer system, a plurality of sequencing reads generated from a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample; (b) aligning, by the computer system, at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; and (c) identifying, by the computer system, germline variants present in the sample at a variant allele fraction (MAF) for at least a portion of the plurality of aligned sequence reads. (d) determining, by the computer system, a quantitative measure of the set of germline variants that is between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) detecting, by the computer system, the presence or absence of the allelic imbalance in the sample based on a predetermined criterion by filtering the set of germline variants identified in (c) based on at least the quantitative measure of (d).

[0020] In some embodiments, the detecting in (e) comprises (f) detecting, by the computer system, one or more quantitative measurements indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, wherein the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes.

[0021] In some embodiments, the method further comprises generating a report, optionally including information about the presence or absence of the allelic imbalance in the sample and / or information about the presence or absence of the contamination or second genome in the sample. In some embodiments, the method further comprises communicating the report to a third party (e.g., the subject from whom the sample originated, or a medical professional, etc.).

[0022] Another aspect of the present disclosure provides a non-transitory computer-readable medium that includes machine-executable code that, upon execution by one or more computer processors, performs any of the methods described above or elsewhere herein.

[0023] Another aspect of the present disclosure provides a system including one or more computer processors and a computer memory coupled thereto, the computer memory including machine executable code that, upon execution by the one or more computer processors, performs any of the methods described above or elsewhere herein.

[0024] In certain embodiments, for example, the following items are provided: (Item 1) 1. A method for detecting the presence or absence of allelic imbalance in a sample from a subject, comprising: (a) sequencing a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying, for at least a portion of the plurality of aligned sequence reads, germline variants present in the sample at a variant allele fraction (MAF), wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d) to detect the presence or absence of the allelic imbalance in the sample based on predetermined criteria. A method comprising: (Item 2) 2. The method of claim 1, wherein the detecting in (e) comprises detecting one or more quantitative measurements indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, and the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes. (Item 3) 3. The method of claim 1 or 2, further comprising detecting the presence or absence of contamination or a second genome in the sample if the absence of the allelic imbalance is detected in the sample. (Item 4) 4. The method of any one of items 1 to 3, wherein the set of germline variants comprises at least about 1,000 different germline variants. (Item 5) 5. The method of any one of items 1 to 4, wherein the set of genetic variants comprises genetic variants selected from the group consisting of single nucleotide variants (SNVs), insertions or deletions (indels), and fusions. (Item 6) 6. The method according to any one of items 1 to 5, wherein the sample is a body fluid sample selected from the group consisting of blood, plasma, serum, urine, saliva, mucosal secretions, sputum, stool, and tears. (Item 7) 7. The method according to any one of items 1 to 6, wherein the subject has a disease or disorder. (Item 8) 8. The method of claim 7, wherein the disease is cancer. (Item 9) 9. The method of any one of items 1 to 8, further comprising amplifying the cell-free DNA molecules prior to sequencing. (Item 10) 10. The method of any one of items 1 to 9, further comprising selectively enriching the cell-free DNA molecules, or the amplified cell-free DNA molecules, for a set of loci prior to sequencing. (Item 11) 11. The method of any one of items 1-10, further comprising attaching one or more adapters comprising molecular barcodes to the cell-free DNA molecules prior to sequencing. (Item 12) 12. The method of claim 11, wherein the one or more adaptors are randomly attached to both ends of the cell-free DNA molecule. (Item 13) 12. The method of claim 11, wherein the cell-free DNA molecules are uniquely barcoded with a molecular barcode. (Item 14) 12. The method of claim 11, wherein the cell-free DNA molecules are non-uniquely barcoded with a molecular barcode. (Item 15) 12. The method of claim 11, wherein each molecular barcode comprises a predetermined or semi-random oligonucleotide sequence that, in combination with the diversity of molecules sequenced from the selected region, enables the identification of a unique cell-free DNA molecule. (Item 16) 16. The method of any one of items 1 to 15, wherein the plurality of genomic regions comprises genetic variants found in COSMIC, TCGA (The Cancer Genome Atlas), or ExAC (Exome Aggregation Consortium). (Item 17) 17. The method of any one of items 1 to 16, wherein the plurality of discrete ranges of MAF values ​​comprises a first range of about 3% to about 40% and a second range of about 60% to about 97%. (Item 18) 20. The method of claim 17, wherein the quantitative measurements of (d) include a multiple set of the genetic variants that fall between a plurality of discrete ranges of MAF values. (Item 19) 20. The method of claim 18, wherein the predetermined criteria comprises the quantitative measure of (d) being greater than a predetermined germline variant threshold. (Item 20) 20. The method of claim 19, wherein the predetermined germline variant threshold is about 21. (Item 21) 21. The method of any one of items 2 or 17-20, wherein the one or more quantitative measures indicative of the CNV or diploid genes are selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average. (Item 22) 22. The method of claim 21, wherein the one or more quantitative measures indicative of the CNV or diploid genes comprise two or more quantitative measures selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average. (Item 23) 23. The method of claim 22, wherein the one or more quantitative measures indicative of the CNV or diploid genes comprise three or more quantitative measures selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, diploid gene fraction, and copy number average. (Item 24) 24. The method of any one of items 21 to 23, wherein the predetermined criteria include one or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than about 3%. (Item 25) 25. The method of claim 24, wherein the predetermined criteria comprises two or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than about 3%. (Item 26) 26. The method of claim 25, wherein the predetermined criteria comprises three or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than about 3%. (Item 27) 27. The method of claim 26, wherein the predetermined criteria include the following: a maximum CNV level across all samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across all samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average for the same germline variant is greater than a predetermined copy number average threshold and a MAF for the same germline variant is less than about 3%. (Item 28) 28. The method of any one of items 24 to 27, wherein the predetermined criteria comprises one or more thresholds selected from the group consisting of a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. (Item 29) 29. The method of claim 28, wherein the predetermined criteria comprises two or more thresholds selected from the group consisting of a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. (Item 30) 30. The method of claim 29, wherein the predetermined criteria comprises three or more thresholds selected from the group consisting of a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. (Item 31) 31. The method of claim 30, wherein the predetermined criteria include a maximum CNV threshold of about 0.22, a minimum CNV threshold of about -0.14, a diploid fraction threshold of about 0.7, and a copy number average threshold of about 10. (Item 32) 4. The method of claim 3, further comprising detecting the presence of the contamination or the second genome in the sample with a positive predictive value (PPV) of at least about 60%. (Item 33) 4. The method of claim 3, further comprising detecting the absence of said contaminant or said second genome in said sample with a negative predictive value (NPV) of at least about 90%. (Item 34) 4. The method of claim 3, further comprising detecting the presence of the contamination or the second genome in the sample with a sensitivity of at least about 90%. (Item 35) 35. The method of claim 34, further comprising detecting the presence of the contamination or the second genome in the sample with a sensitivity of at least about 99%. (Item 36) 4. The method of claim 3, further comprising detecting the absence of the contamination or the second genome in the sample with a specificity of at least about 35%. (Item 37) The germline variant is (i) determining total and mutant allele counts for nucleic acid variants from the cfDNA molecules; (ii) identifying associated variables of said nucleic acid variants from said cfDNA molecules; (iii) determining a quantitative value for said associated variable of said nucleic acid variant; (iv) generating a statistical model for the expected number of germline variant alleles at the genomic locus of the nucleic acid variant; (v) generating a probability value for the nucleic acid variant based at least in part on the statistical model for the predicted germline variant allele count, the quantitative value for the associated variable for the nucleic acid variant, and at least one of the total allele count and the variant allele count for the nucleic acid variant; and (vi) classifying the nucleic acid variant as (1) of somatic origin if the p-value for the nucleic acid variant is less than a predetermined threshold, or (2) of germline origin if the p-value for the nucleic acid variant is equal to or greater than a predetermined threshold. 37. The method according to any one of items 1 to 36, further comprising identifying the target gene by: (Item 38) 38. The method of any one of items 1 to 37, further comprising detecting an allele-specific loss in the sample based on at least one of the set of germline variants identified as present in a given MAF in (c). (Item 39) 39. The method of claim 38, wherein the allele-specific loss in the sample is detected based on the at least one of the set of germline variants being present in the sample from the subject at a MAF of less than 50%. (Item 40) 40. The method of claim 39, wherein the allele-specific loss in the sample is detected based on the at least one of the set of germline variants being present in the sample from the subject and in each of one or more additional samples from one or more subjects at a MAF below 50%. (Item 41) 41. The method of any one of items 38 to 40, wherein the at least one of the sets of germline variants is found in COSMIC, TCGA (The Cancer Genome Atlas), or ExAC (Exome Aggregation Consortium). (Item 42) 42. The method of claim 41, wherein the at least one of the set of germline variants is a BRCA1 gene variant. (Item 43) 43. The method of claim 42, wherein the BRCA1 gene variant is BRCA1 P209L. (Item 44) 44. The method according to any one of items 1 to 43, wherein at least a part of the method is executed by a computer system. (Item 45) A system, which when executed by at least one electronic processor, (a) obtaining a plurality of sequence reads corresponding to a plurality of cell-free deoxyribonucleic acid (DNA) molecules from a sample of a subject; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying, for at least a portion of the plurality of aligned sequence reads, germline variants present in the sample at a variant allele fraction (MAF), wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d) to detect the presence or absence of allelic imbalance in the sample based on predefined criteria. A system including a controller including a computer-readable medium including non-transitory computer-executable instructions for implementing the above-mentioned, or a controller that can access the computer-readable medium. (Item 46) 46. ​​The system of claim 45, wherein the detecting in (e) comprises detecting one or more quantitative measurements indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, and the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes. (Item 47) 47. The system of claim 45 or 46, further comprising a nucleic acid sequencer operably connected to the controller, the nucleic acid sequencer configured to process the plurality of cell-free DNA molecules from the sample to generate the plurality of sequence reads. (Item 48) 48. The system of any one of items 45 to 47, further comprising: generating a report, optionally including information about the presence or absence of the allelic imbalance in the sample and / or information about the presence or absence of the contamination or second genome in the sample, when the non-transitory computer executable instructions, when executed by at least one electronic processor, are configured to: (Item 49) 49. The system of claim 48, wherein the non-transitory computer-executable instructions, when executed by at least one electronic processor, further perform the following: communicating the report to a third party (e.g., the subject from whom the sample originated, or a medical professional, etc.). (Item 50) 1. A method for detecting the presence or absence of allelic imbalance in a sample from a subject, comprising: (a) accessing, by a computer system, a plurality of sequencing reads generated from a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample; (b) aligning, by the computer system, at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying, by the computer system, germline variants present in the sample at a variant allele fraction (MAF) for at least a portion of the plurality of aligned sequence reads, thereby identifying a set of germline variants in the sample, each germline variant in the set of germline variants having a corresponding MAF value; (d) determining, by the computer system, a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) detecting, by the computer system, the presence or absence of the allelic imbalance in the sample based on predetermined criteria by filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d). A method comprising: (Item 51) 51. The method of claim 50, wherein the detecting in (e) comprises detecting, by the computer system, one or more quantitative measurements indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, and the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes. (Item 52) 52. The system of any one of items 1 to 44 or 50 to 51, further comprising generating a report optionally comprising information about the presence or absence of the allelic imbalance in the sample and / or information about the presence or absence of the contamination or second genome in the sample. (Item 53) 53. The method of claim 52, further comprising communicating the report to the subject from whom the sample originated or to a third party, such as a medical professional. Additional aspects and advantages of the present disclosure will become readily apparent to those skilled in the art from the following detailed description, wherein only illustrative embodiments of the present disclosure are shown and described. As will be realized, the present disclosure is capable of other and different embodiments, and its several details are capable of modifications in various obvious respects, all without departing from the present disclosure. Accordingly, the drawings and description are to be regarded as illustrative in nature, and not as restrictive. [Brief description of the drawings]

[0025] [Figure 1] FIG. 1 shows an example of the methods provided herein.

[0026] [Diagram 2] Figure 2 shows an example workflow for detecting allelic imbalance or contamination in cell-free DNA samples.

[0027] [Diagram 3] FIG. 3 is a diagram illustrating a computer system programmed or otherwise configured to perform the methods provided herein. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0028] definition While various embodiments of the present disclosure have been shown and described herein, those skilled in the art will understand that such embodiments are presented by way of example only. Numerous variations, changes, and substitutions may occur to those skilled in the art without departing from the present disclosure. It should be understood that various alternatives to the embodiments of the present disclosure described herein may be employed.

[0029] Adapter: The term "adapter" refers to a short nucleic acid (e.g., less than 500 nucleotides, less than 100 nucleotides, or less than 50 nucleotides in length), usually at least partially double-stranded, for attachment to either or both ends of a sample nucleic acid molecule. The adapter may include primer binding sites that allow amplification of the nucleic acid molecule flanked by the adapter at both ends, and / or sequencing primer binding sites, including primer binding sites for next generation sequencing (NGS). The adapter may also include binding sites for capture probes, such as oligonucleotides attached to a flow cell support. The adapter may also include a tag, as described above. The tag is preferably positioned relative to the primer and sequencing primer binding sites such that the tag is included in the amplicon and sequencing reads of the nucleic acid molecule. The same or different adapters can be ligated to each end of the nucleic acid molecule. Optionally, the same adapters can be ligated to each end, except that the tags are different. A preferred adaptor is a Y-shaped adaptor, one end of which is blunt or overhanging, for binding to a nucleic acid molecule, which is also blunt or has one or more overhanging complementary nucleotides. Another preferred adaptor is a bell-shaped adaptor, which also has a blunt or overhanging end for binding to the nucleic acid to be analyzed.

[0030] Allelic imbalance: The term "allelic imbalance" generally refers to the difference in DNA levels between two alleles in a gene (e.g., as a result of loss of heterozygosity). Allelic imbalance can occur when the ratio of DNA levels between two alleles in a gene is not about 1. For example, allelic imbalance can occur as a result of genetic imprinting, in which epigenetics and environmental factors can affect the expression of one or both alleles in a given gene. As another example, cis-acting mutations can affect the regulation of one allele of a pair of alleles in a gene (e.g., by changes in promoter or enhancer regions (e.g., transcription factor binding sites) or changes to the 3'UTR region).

[0031] Allelic imbalance candidate: The term "allelic imbalance candidate" generally refers to a sample that is being analyzed (e.g., using the methods, systems, and media of the present disclosure) to detect the presence or absence of allelic imbalance or contamination.

[0032] Cell-free nucleic acid: The phrase "cell-free nucleic acid" may refer to nucleic acid that is not contained in or otherwise bound to a cell, in other words, nucleic acid remaining in a sample from which intact cells have been removed. Cell-free nucleic acid may refer to all unencapsulated nucleic acid originating from a body fluid (e.g., blood, urine, CSF, etc.) from a subject. Cell-free nucleic acid includes DNA (cfDNA), RNA (cfRNA), and hybrids thereof, including genomic DNA, mitochondrial DNA, circulating DNA, siRNA, miRNA, circulating RNA (cRNA), tRNA, rRNA, nucleolar RNA (snoRNA), Piwi-binding RNA (piRNA), long non-coding RNA (long ncRNA), or fragments of any of these. Cell-free nucleic acid may be double-stranded, single-stranded, or hybrids thereof. Cell-free nucleic acid may be released into body fluids through secretion or cell death processes (e.g., necrosis and apoptosis of cells). Cell-free nucleic acid may be found in exosomes. Some cell-free nucleic acids may be released from cancer cells into bodily fluids (e.g., circulating tumor DNA (ctDNA)). Other cell-free nucleic acids are released from healthy cells. ctDNA may be unencapsulated tumor-derived fragmented DNA. Cell-free fetal DNA (cffDNA) is fetal DNA that is freely circulating in the maternal bloodstream. Cell-free nucleic acids may have one or more epigenetic modifications. For example, cell-free nucleic acids may be acetylated, 5-methylated, ubiquitinated, phosphorylated, sumoylated, ribosylated, and / or citrullinated.

[0033] Contamination: The term "contamination" refers to any chemical or digital contamination of one sample by another sample. Contamination can come from a variety of sources, including, but not limited to, (1) assay-level contamination, such as physical carryover of liquid between samples (e.g., pipetting, automated liquid handling by sample preparation equipment or sequencer, handling of amplified materials); demultiplexing artifacts (e.g., base calling errors confounding sample indexes with poor pairwise Hamming distances; insertions / deletions confounding sample indexes with poor pairwise Hamming distances); reagent impurities (e.g., sample index oligos that are missing a level of oligos synthesized in the same batch; sample index oligos contaminated by oligos containing another sample index (either through carryover of synthesis errors); or (2) a sample containing a second genome.

[0034] Copy number variant: As used herein, "copy number variant," "CNV," or "copy number variation" refers to the phenomenon in which sections of the genome are repeated and the number of repeats in the genome varies between individuals within a population under consideration and between two conditions or states of an individual (e.g., CNVs can be different in an individual before and after treatment).

[0035] Deoxyribonucleic acid and ribonucleic acid: The term "DNA" refers to a natural or modified nucleotide that has a hydrogen group at the 2' position of the sugar moiety. DNA typically includes a chain of nucleotides that includes four types of nucleotide bases: adenine (A), thymine (T), cytosine (C), and guanine (G). As used herein, "ribonucleic acid" or "RNA" refers to a natural or modified nucleotide that has a hydroxyl group at the 2' position of the sugar moiety. RNA typically includes nucleotides that include four types of nucleotides: A, uracil (U), G, and C. As used herein, the term "nucleotide" refers to a natural or modified nucleotide. Certain pairs of nucleotides specifically bind to each other in a complementary manner (called complementary base pairing). In DNA, adenine (A) pairs with thymine (T) and cytosine (C) pairs with guanine (G). In RNA, adenine (A) pairs with uracil (U) and cytosine (C) pairs with guanine (G). When a first nucleic acid strand binds to a second nucleic acid strand that is composed of nucleotides complementary to the first strand, the two strands combine to form a duplex.

[0036] Germline variant: The terms "germline variant(s)" or "germline mutation(s)" are used interchangeably and refer to a heritable mutation (i.e., a mutation that does not occur after conception). A germline mutation may be the only mutation that can be inherited by offspring and may be present in every somatic and germline cell of the offspring.

[0037] Loss of heterozygosity: The term "loss of heterozygosity" (LOH) generally refers to a form of allelic imbalance in which one allele of an allele pair at a locus is completely lost. LOH can occur by many genetic mechanisms, such as physical deletion, chromosomal nondisjunction, mitotic nondisjunction followed by doubling of the remaining chromosome, mitotic recombination, and gene conversion. LOH can be detected based on the measurement of the mutant allele ratio or minor allele frequency at a locus. LOH can occur, for example, when a tumor suppressor gene is inactivated such that one allele of an allele pair at the tumor suppressor gene is mutated and the other allele is lost.

[0038] Minor allele frequency: As used herein, "minor allele frequency" refers to the frequency of a minor allele (e.g., not the most common allele) occurring in a given population of nucleic acids (e.g., a sample obtained from a subject). Genetic variants with low minor allele frequency are typically present at a relatively low frequency in samples.

[0039] Variant allele count: The term "variant allele count" refers to the number of nucleic acid molecules in a plurality of nucleic acid molecules (e.g., obtained or derived from a sample) that have a variant allele or allelic alteration at a particular genomic locus.

[0040] Mutant allele fraction: The phrase "mutant allele fraction", "mutation dose", or "MAF" refers to the fraction of nucleic acid molecules in a given sample that have an allelic alteration or mutation at a given genomic location. MAF is generally expressed as a fraction or percentage. For example, MAF is typically less than about 0.5, 0.1, 0.05, or 0.01 (i.e., less than about 50%, 10%, 5%, or 1%) of the total cellular variants or alleles present at a given locus.

[0041] Nucleic acid sequencing data: As used herein, "nucleic acid sequencing data", "nucleic acid sequencing information", "nucleic acid sequence", "nucleotide sequence", "genomic sequence", "gene sequence", "sequence information", or "fragment sequence", or "nucleic acid sequencing read" refers to any information or data that indicates the order of nucleotide bases (e.g., adenine, guanine, cytosine, and thymine or uracil) in a molecule of nucleic acid such as DNA or RNA (e.g., a whole genome, a whole transcriptome, an exome, an oligonucleotide, a polynucleotide, or a fragment). It should be understood that the present teachings contemplate sequence information obtained using any type of available technique, platform, or technology, including but not limited to capillary electrophoresis, microarrays, ligation-based systems, polymerase-based systems, hybridization-based systems, direct or indirect nucleotide discrimination systems, pyrosequencing, ion- or pH-based systems, and electronic signature-based systems.

[0042] Nucleic acid tag: As used herein, "nucleic acid tag" refers to a short nucleic acid (e.g., less than n nucleotides in length, where n is about 500 nucleotides, about 100 nucleotides, about 50 nucleotides, or about 10 nucleotides in length) that is used to distinguish nucleic acids from different samples (e.g., represents a sample index) or to distinguish different types or differently processed nucleic acid molecules in the same sample (e.g., represents a molecular barcode). Nucleic acid tags comprise a predetermined, predefined, non-random, random, or semi-random oligonucleotide sequence. Such nucleic acid tags can be used to label different nucleic acid molecules, or different nucleic acid samples or subsamples. Nucleic acid tags can be single-stranded, double-stranded, or at least partially double-stranded. Nucleic acid tags can have equal or different lengths, as appropriate. Nucleic acid tags may also include double-stranded molecules with one or more blunt ends, may include 5' or 3' single-stranded regions (e.g., overhangs), and / or may include one or more other single-stranded regions at other sites within a given molecule. Nucleic acid tags can be attached to one or both ends of other nucleic acids (e.g., sample nucleic acids to be amplified and / or sequenced). Nucleic acid tags can be decoded to reveal information such as the sample origin, form, or processing of a given nucleic acid. For example, nucleic acid tags can be used to enable the storage and / or parallel processing of multiple samples containing nucleic acids with different molecular barcodes and / or sample indexes, which are then analyzed by detecting (e.g., reading) the nucleic acid tags. Nucleic acid tags are also referred to as identifiers (e.g., molecular identifiers, sample identifiers). Additionally or alternatively, nucleic acid tags can also be used as molecular identifiers (e.g., to distinguish between different molecules or amplicons of different parent molecules in the same sample or subsample). This includes, for example, uniquely tagging different nucleic acid molecules in a given sample, or non-uniquely tagging such molecules.In the case of non-unique tagging amplification, a limited number of tags (i.e., molecular barcodes) may be used to tag each nucleic acid molecule such that different molecules, in combination with at least one molecular barcode, can be identified based on their endogenous sequence information (e.g., start and / or end positions where they are mapped to a selected reference genome, subsequences at one or both ends of the sequence, and / or length of the sequence). Typically, a sufficient number of different molecular barcodes are used such that the probability that any two molecules have the same endogenous sequence information (e.g., start and / or end positions, subsequences at one or both ends of the sequence, and / or length) and the same molecular barcode is low (e.g., less than about 10%, less than about 5%, less than about 1%, less than about 0.1%, less than about 0.01%, less than about 0.001%, or less than 0.0001% probability).

[0043] Polynucleotide: "Polynucleotide", "nucleic acid", "nucleic acid molecule", or "oligonucleotide" refers to a linear polymer of nucleosides (including deoxyribonucleosides, ribonucleosides, or analogs thereof) linked by internucleoside bonds. Typically, a polynucleotide contains at least three nucleosides. Oligonucleotides often range in size from a few monomeric units (e.g., 3-4) to hundreds of monomeric units. When a polynucleotide is represented by a sequence of letters (e.g., "ATGCCTG"), it is understood that unless otherwise indicated, the nucleotides are always in a 5'→3' orientation from left to right of the string, with "A" referring to deoxyadenosine, "C" referring to deoxycytidine, "G" referring to deoxyguanosine, and "T" referring to thymidine. The letters A, C, G, and T may be used to refer to the bases themselves or to the nucleosides or nucleotides that contain those bases, as is standard in the art.

[0044] Reference sequence: The phrase "reference sequence" refers to a known sequence that is used for the purpose of comparison with an experimentally determined sequence. For example, the known sequence can be a whole genome, a chromosome, or any fragment thereof. The reference typically comprises at least 20, 50, 100, 200, 250, 300, 350, 400, 450, 500, 1000, 10000, 50000, 100000, or more nucleotides. The reference sequence may be aligned with a single continuous sequence of a genome or chromosome, or may comprise discontinuous segments that align with different regions of a genome or chromosome. Reference human genomes include, for example, hG19 and hG38.

[0045] Second Genome: The term "second genome" refers to a nucleic acid sequence that is present in a subject but is associated with a genome that is not the subject's genome. Such genomes include, but are not limited to, genomes derived from grafts, viruses, therapeutic nucleic acid constructs, blood transfusions, fetuses, etc.

[0046] Sequencing: As used herein, the term "sequencing" or "sequencer" refers to any of a number of techniques used to determine the sequence of a biological molecule (e.g., a nucleic acid such as DNA or RNA). Exemplary sequencing methods include, but are not limited to, targeted sequencing, single molecule real-time sequencing, exon sequencing, electron microscope-based sequencing, panel sequencing, transistor-mediated sequencing, direct sequencing, random shotgun sequencing, Sanger dideoxy termination sequencing, whole genome sequencing, sequencing by hybridization, pyrosequencing, capillary electrophoresis, duplex sequencing, cycle sequencing, single base extension sequencing, solid-phase sequencing, and the like. Examples of the sequencing methods include PCR, high-throughput sequencing, massively parallel signature sequencing, emulsion PCR, co-amplification PCR at low denaturation temperature (COLD-PCR), multiplex PCR, reversible dye terminator sequencing, paired-end sequencing, near-term sequencing, exonuclease sequencing, ligation sequencing, short-read sequencing, single molecule sequencing, synthesis sequencing, real-time sequencing, reverse terminator sequencing, nanopore sequencing, 454 sequencing, Solexa Genome Analyzer sequencing, SOLiD™ sequencing, MS-PET sequencing, and combinations thereof. In some embodiments, sequencing can be performed by a commercially available genetic analyzer, such as, for example, Illumina or Applied Biosystems. The phrase "next-generation sequencing" or "NGS" refers to sequencing technologies that have increased throughput (e.g., the ability to generate hundreds of thousands of relatively small sequence reads at one time) compared to traditional Sanger sequencing or capillary electrophoresis-based techniques.Some examples of next generation sequencing techniques include, but are not limited to, sequencing by synthesis, sequencing by ligation, and sequencing by hybridization.

[0047] Subject: The term "subject" may refer to an animal, such as a mammalian species (preferably a human) or an avian (e.g., avian) species, or other organism (particularly a diploid organism). More specifically, the subject may be a mammal, such as a vertebrate, e.g., a mouse, a primate, a monkey, or a human. Animals include farm animals, sport animals, and pets. A subject may be a healthy individual, an individual having symptoms or signs, suspected of having a disease or a propensity for a disease, or an individual in need of treatment or suspected of being in need of treatment. [Mode for carrying out the invention]

[0048] I. Overview In cancer patients, allelic imbalance may be caused by loss of heterozygosity and may result in a different mutant allele fraction (MAF) distribution in the assay of cell-free nucleic acid samples from subjects compared to samples without allelic imbalance. For example, samples with allelic imbalance may contain germline variants with very low MAF. Germline variants with low MAF may also be observed when samples are contaminated, for example during processing for sequencing, or when samples contain a second genome (other than the subject's genome), for example from a transplant, blood transfusion, or fetus. Therefore, problems may be encountered when distinguishing between allelic imbalance samples and samples with contamination or containing a second genome.

[0049] When assaying cell-free nucleic acid from samples containing contamination or second genome, such samples may require additional manual review or additional sequencing runs. As a result, failure to distinguish between allelic imbalance samples and contaminated or second genome samples may significantly increase the cost and time required to reliably assay such samples. The present disclosure provides methods and systems for identifying allelic imbalance or contamination in cell-free nucleic acid samples. These methods and systems may identify allelic imbalance or contamination by obtaining and analyzing quantitative measurements of small variants and copy number variations.

[0050] The present disclosure provides a method and system for detecting allelic imbalance in a sample from a subject.In one embodiment, the present disclosure provides a method for detecting allelic imbalance in a sample from a subject, comprising: (a) sequencing a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) for at least a portion of the plurality of aligned sequence reads, identifying a set of germline variants in the sample by identifying germline variants present in the sample at a mutant allele fraction (MAF), where each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measurement value of the set of germline variants identified in (c) that is between a plurality of separate ranges of MAF values; and (e) filtering the set of germline variants identified in (c) based on at least the quantitative measurement value of (d), thereby detecting the allelic imbalance in the sample based on a predetermined criterion.

[0051] In some embodiments, the method further comprises (f) detecting one or more quantitative measurements from the plurality of aligned sequence reads indicative of copy number variation (CNV) or diploid genes, wherein the predetermined criteria comprises the one or more quantitative measurements indicative of the CNV or the diploid genes.

[0052] In some embodiments, the method further comprises detecting contamination in the sample if no allelic imbalance is detected in the sample.

[0053] FIG. 1 illustrates an example of a method 100 provided herein. The method 100 may include sequencing DNA molecules from a sample in which allelic imbalance or contamination is to be detected (as in operation 102) to generate sequence reads. The method 100 may then include aligning at least some of the sequence reads to a reference sequence to generate aligned sequence reads (as in operation 104). The method 100 may then include identifying a set of germline variants in the sample and their corresponding MAF values ​​for at least some of the aligned sequence reads (as in operation 106), or in certain embodiments, corresponding minor allele frequency values. The method 100 may then determine quantitative measurements of the germline variants that are between a plurality of discrete ranges of MAF values, or in certain embodiments, that are within discrete ranges of minor allele frequency values ​​(as in operation 108). Next, method 100 may include detecting allelic imbalance in the sample based on predetermined criteria by filtering the germline variants based on at least the quantitative measurements (as in operation 110).

[0054] The methods and systems provided herein may be particularly useful in the analysis of cell-free nucleic acid molecules (e.g., DNA or RNA molecules). In some cases, cell-free nucleic acid molecules may be extracted and isolated from a biological sample from a subject and may be readily available. Biological samples may include body fluid samples selected from the group including, but not limited to, blood, plasma, serum, urine, saliva, mucosal secretions, sputum, stool, and tears. Cell-free nucleic acid molecules may be extracted using a variety of methods, including, but not limited to, isopropanol precipitation and / or silica-based purification.

[0055] Biological samples may be collected from a number of subjects (e.g., disease-free subjects, subjects at risk for, showing symptoms of, or having a disease, such as cancer or a virus, or subjects at risk for, showing symptoms of, or having a genetic disorder. In some embodiments, the disease or disorder is selected from the group consisting of immunodeficiency disorders, hemophilia, thalassemia, sickle cell disease, blood disorders, chronic granulomatous disorders, congenital blindness, lysosomal storage disorders, muscular dystrophy, cancer, neurodegenerative disorders, viral infections, bacterial infections, epidermolysis bullosa, heart disease, lipid metabolism disorders, and diabetes, or a combination thereof.

[0056] After obtaining or preparing cell-free nucleic acid molecules, the cell-free nucleic acid molecules can be subjected to any of a number of different library preparation procedures to prepare nucleic acid molecules for sequencing.The cell-free nucleic acid molecules can be treated with one or more reagents (e.g., enzymes, adaptors, tags (e.g., barcodes), probes, etc.) before sequencing.The tagged molecules can then be used for downstream applications, such as sequencing reactions, which can track individual molecules.

[0057] In some embodiments, the method may further comprise an enrichment step prior to sequencing, whereby regions of the tagged molecules are selectively or non-selectively enriched.

[0058] Once sequencing data of the cell-free nucleic acid molecules has been collected, one or more bioinformatics processes may be applied to the sequence data to detect allelic imbalance or contamination of the cell-free nucleic acid sample.

[0059] In some cases, the sequence reads generated from the sequencing reaction can be aligned to a reference sequence to perform bioinformatics analysis. In various aspects of bioinformatics analysis, one or more thresholds can be set to ensure quality. For example, an alignment threshold can be set so that only highly homologous sequence reads (e.g., 10 or less mismatches between the reference sequence and the sequence read) are mapped to the reference sequence. In some cases, sequence reads that do not meet the quality threshold can be removed, for example, based on the chromatogram of the sequence reads. In some cases, the copy number or amount of a given sequence can be quantified based on the number of sequence reads that are mapped or aligned to the given sequence. In some cases, the overrepresentation of a sequence can be determined by comparing the copy number or amount of different sequences within all sequence reads.

[0060] In certain embodiments, the sample may be contacted with a sufficient number of adapters to make it unlikely (e.g., less than about 1%, less than about 0.1%, less than about 0.01%, less than about 0.001%, or less than about 0.0001%) that any two copies of the same nucleic acid will receive the same combination of adapter molecule barcodes or tags derived from the adapters attached to one or both ends. Using adapters in this manner allows sequence reads that have the same start and end points aligned (or mapped) to a reference sequence and are attached to the same combination of barcodes to be grouped into families of reads generated from the same original molecule. Such families may represent the sequences of the amplification products of the nucleic acids in the sample before amplification.

[0061] In some embodiments, the sequences of the family members modified by blunting and adaptor ligation may be compiled to derive a consensus nucleotide or a complete consensus sequence of the nucleic acid molecules in the original sample. In other words, a nucleotide occupying a particular position of the nucleic acid in the sample may be determined to be the consensus of the nucleotide occupying the corresponding position in the family member sequences. The consensus nucleotide may be determined by methods such as voting or confidence score, to name two non-limiting exemplary methods. A family may include sequences of one or both strands of a double-stranded nucleic acid. When a family member includes sequences of both strands from a double-stranded nucleic acid, the sequence of one strand is converted to its complementary sequence in order to compile the entire sequence and derive a consensus nucleotide or consensus sequence. Some families may include only a single member sequence. In this case, this sequence may be interpreted as the sequence of the nucleic acid in the sample before amplification. Alternatively, a family with only a single member sequence may be excluded from subsequent analysis.

[0062] The reference sequence may be one or more known sequences, for example, a known whole genome sequence or partial genome sequence from a subject, a whole genome sequence of a human subject. The reference sequence may be hG19. The sequenced nucleic acid may represent a sequence determined directly for the nucleic acid in the sample, or a consensus of sequences of an amplification product of such nucleic acid, as described above. The comparison may be performed at one or more designated positions of interest in the reference sequence. A subset of sequenced nucleic acids may be identified, including positions that correspond to the designated positions of the reference sequence when each sequence is maximally aligned. Within such a subset, it may be determined which, if any, sequenced nucleic acids contain a nucleotide mutation at the designated position, as well as, if necessary, which, if any, sequenced nucleic acids contain a reference nucleotide (i.e., the same as in the reference sequence). If the number of sequenced nucleic acids in the subset that contains a nucleotide variant exceeds a threshold, the mutated nucleotide may be called at the designated position. The threshold value may be a simple number, e.g., at least 1, 2, 3, 4, 5, 6, 7, 9, or 10 sequenced nucleic acids in the subset that contain nucleotide variants, among other possibilities, or may be a ratio, e.g., at least 0.5, 1, 2, 3, 4, 5, 10, 15, or 20 sequenced nucleic acids in the subset that contain nucleotide variants. The comparison may be repeated for any designated positions of interest in the reference sequence. Optionally, the comparison may be performed for designated positions that occupy at least 20, 100, 200, or 300 consecutive positions on the reference sequence, e.g., 20-500, or 50-300 consecutive positions.

[0063] The present disclosure also provides a system for carrying out or implementing the methods described herein. In certain embodiments, the system includes: (a) a nucleic acid sequencer that generates sequencing reads as signals from adaptor-tagged cfDNA molecules derived from one or more samples, where the adaptors include barcodes that identify redundant sequence reads derived from the same original cfDNA molecule together with start and end information from the cfDNA molecules; and (b) a computer that communicates with the nucleic acid sequencer through a communication network, where the computer accepts the signals into a computer memory, where the computer includes a computer processor and a computer readable medium, where the computer readable medium includes machine executable code that, when executed by the computer processor, carries out the method described below, and includes a method comprising: a) generating a sequence read from a plurality of cell-free data derived from the samples; a) sequencing an oxyribonucleic acid (DNA) molecule to generate a plurality of sequence reads; b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; c) for each of a plurality of genomic regions, determining from the plurality of aligned sequence reads a variant allele fraction (MAF) of the genomic region of the sample; d) for each of the plurality of genomic regions, determining from the plurality of aligned sequence reads whether the genomic region is a germline variant; e) determining a quantitative measure of the determined germline variant of the plurality of genomic regions that is between a plurality of discrete ranges of MAF values; and f) detecting allelic imbalance in the sample based on a predetermined criteria comprising the quantitative measure of the determined germline variant.

[0064] In some embodiments, the method performed by the computer processor further comprises grouping the sequence reads into families (each family comprises sequence reads that comprise the same barcode and have the same start and end positions), whereby each family comprises amplified sequence reads derived from the same original cfDNA molecule.

[0065] In some embodiments, the sequencer is a DNA sequencer. In some embodiments, the sequencer is designed for high-throughput sequencing, such as next-generation sequencing. In some embodiments, the system includes adaptor-tagged cfDNA molecules in the sequencer. In some embodiments, the adaptor-tagged cfDNA molecules are derived from one subject or multiple subjects. In some embodiments, the cfDNA molecules derived from the sample have unique or non-unique barcodes. II. General Features of the Method and System A. Sample

[0066] The sample can be any biological sample isolated from a subject. Samples can include body tissue, whole blood, platelets, serum, plasma, stool, red blood cells, white blood cells or leucocytes, endothelial cells, tissue biopsies (e.g., biopsies from known or suspected solid tumors), cerebrospinal fluid, synovial fluid, lymphatic fluid, peritoneal fluid, interstitial or extracellular fluid (e.g., interstitial fluid), gingival exudate, gingival crevicular fluid, bone marrow, pleural fluid, cerebrospinal fluid, saliva, mucus, sputum, semen, sweat, and urine. Samples are preferably body fluids, particularly blood and its fractions, and urine. Such samples include nucleic acids shed from tumors. Nucleic acids can include DNA and RNA, and can be in double-stranded and single-stranded forms. Sample may be in the form that is isolated from subject, or may be further processed to remove or add components such as cells, enrich one component with respect to another, or convert one form of nucleic acid into another (e.g., RNA into DNA, or single-stranded nucleic acid into double-stranded nucleic acid).Thus, for example, the body fluid for analysis is plasma or serum, which contains cell-free nucleic acid, such as cell-free DNA (cfDNA).

[0067] In some embodiments, the sample volume of bodily fluid taken from the subject depends on the desired read depth of the region to be sequenced. Exemplary volumes are about 0.4-40 ml, about 5-20 ml, about 10-20 ml. For example, the volume can be about 0.5 ml, about 1 ml, about 5 ml, about 10 ml, about 20 ml, about 30 ml, about 40 ml, or more (in milliliters). The volume of sampled plasma is typically between about 5 ml and about 20 ml.

[0068] A sample can contain varying amounts of nucleic acid. Typically, the amount of nucleic acid in a given sample is equal to the amount of multiple genome equivalents. For example, a sample of about 30 ng of DNA will contain about 10,000 (10 4 ) haploid human genome equivalents, whereas cfDNA contains approximately 200 billion (2 × 10 11Similarly, a sample of about 100 ng of DNA may contain about 30,000 haploid human genome equivalents, whereas in the case of cfDNA, it may contain about 600 billion individual molecules.

[0069] In some embodiments, the sample may contain nucleic acids from different sources, such as nucleic acids from cells and nucleic acids from acellular sources (such as blood samples). Typically, the sample contains nucleic acids with mutations. For example, the sample may contain DNA with germline mutations and / or somatic mutations. Typically, the sample contains DNA with cancer-associated mutations (e.g., cancer-associated somatic mutations).

[0070] Exemplary amounts of cell-free nucleic acid in a sample prior to amplification typically range from about 1 femtogram (fg) to about 1 microgram (μg), e.g., from about 1 picogram (pg) to about 200 nanograms (ng), from about 1 ng to about 100 ng, from about 10 ng to about 1000 ng. In some embodiments, the sample contains up to about 600 ng, up to about 500 ng, up to about 400 ng, up to about 300 ng, up to about 200 ng, up to about 100 ng, up to about 50 ng, or up to about 20 ng of cell-free nucleic acid molecules. Optionally, the amount is at least about 1 fg, at least about 10 fg, at least about 100 fg, at least about 1 pg, at least about 10 pg, at least about 100 pg, at least about 1 ng, at least about 10 ng, at least about 100 ng, at least about 150 ng, or at least about 200 ng of cell-free nucleic acid molecules. In certain embodiments, the amount is up to about 1 fg, about 10 fg, about 100 fg, about 1 pg, about 10 pg, about 100 pg, about 1 ng, about 10 ng, about 100 ng, about 150 ng, or about 200 ng of cell-free nucleic acid molecules. In some embodiments, the method includes obtaining between about 1 fg and about 200 ng of cell-free nucleic acid molecules from the sample.

[0071] Cell-free nucleic acids typically have a size distribution between about 100 and about 500 nucleotides in length, with about 90% of the molecules in a sample being between about 110 and about 230 nucleotides in length, with a mode of about 168 nucleotides in length, and a second minor peak within the length range of about 240 to about 440 nucleotides. In certain embodiments, the cell-free nucleic acids are about 160 to about 180 nucleotides in length, about 320 to about 360 nucleotides in length, or about 440 to about 480 nucleotides in length.

[0072] In some embodiments, the cell-free nucleic acid is separated from the body fluid by a partitioning step that separates the cell-free nucleic acid as found in solution from intact cells and other insoluble components of the body fluid. In some of these embodiments, partitioning includes techniques such as centrifugation or filtration. Alternatively, the cells in the body fluid are lysed and the cell-free and cellular nucleic acids are processed together. Generally, after addition of buffer and washing steps, the cell-free nucleic acid is precipitated, for example, with alcohol. In certain embodiments, additional purification steps, for example, silica-based columns, are used to remove contaminants or salts. For example, non-specific bulk carrier nucleic acid may be added throughout the reaction as needed to optimize certain aspects of the exemplary procedure, such as yield. After such processing, the sample typically contains various forms of nucleic acid, including double-stranded DNA, single-stranded DNA, and / or single-stranded RNA. If necessary, the single-stranded DNA and / or single-stranded RNA are converted to double-stranded form for inclusion in subsequent processing and analysis steps. B. Nucleic Acid Tags

[0073] In some embodiments, nucleic acid molecules (from a sample of polynucleotides) may be tagged with a sample index and / or a molecular barcode (commonly referred to as a "tag"). The tag may be incorporated or otherwise attached to an adapter by, among other methods, chemical synthesis, ligation (e.g., blunt-end ligation or sticky-end ligation), or overlap extension polymerase chain reaction (PCR). Such an adapter may ultimately be attached to the target nucleic acid molecule. In other embodiments, one or more iterations of an amplification cycle (e.g., PCR amplification) are generally applied to introduce the sample index into the nucleic acid molecule using conventional nucleic acid amplification methods. The amplification may be performed in one or more reaction mixtures (e.g., multiple microwells in an array). The molecular barcode and / or sample index may be introduced simultaneously or in any sequential order. In some embodiments, the molecular barcode and / or sample index are introduced before and / or after performing a sequence capture step. In some embodiments, only the molecular barcode is introduced before probe capture and the sample index is introduced after performing a sequence capture step. In some embodiments, both the molecular barcode and the sample index are introduced before performing the probe-based capture step. In some embodiments, the sample index is introduced after performing the sequence capture step. In some embodiments, the molecular barcode is incorporated into the nucleic acid molecule (e.g., cfDNA molecule) in the sample by ligation (e.g., blunt-end ligation or sticky-end ligation) through an adaptor. In some embodiments, the sample index is incorporated into the nucleic acid molecule (e.g., cfDNA molecule) in the sample by overlap extension polymerase chain reaction (PCR). Typically, the sequence capture protocol involves introducing a single-stranded nucleic acid molecule complementary to a target nucleic acid sequence (e.g., a coding sequence of a genomic region), where mutations in such region are associated with a cancer type.

[0074] In some embodiments, the tag can be located at one or both ends of the sample nucleic acid molecule.In some embodiments, the tag is a predetermined, random, or semi-random sequence oligonucleotide.In some embodiments, the tag is less than about 500, less than 200, less than 100, less than 50, less than 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 nucleotide in length.The tag can be randomly or non-randomly attached to the sample nucleic acid.

[0075] In some embodiments, each sample is uniquely tagged with a sample index or a combination of sample indexes. In some embodiments, each nucleic acid molecule of a sample or sub-sample is uniquely tagged with a molecular barcode or a combination of molecular barcodes. In other embodiments, multiple molecular barcodes may be used that are not necessarily unique to each other (e.g., non-unique molecular barcodes). In these embodiments, the molecular barcodes are generally attached (e.g., by ligation) such that the molecular barcode-sequence combination is attached to each individual molecule to generate a unique sequence that can be tracked individually. Detection of the non-uniquely tagged molecular barcodes in combination with endogenous sequence information (e.g., the beginning (start) and / or end (stop) points corresponding to the sequence of the original nucleic acid molecule in the sample, the subsequence of the sequence read at one or both ends, the length of the sequence read, and / or the length of the original nucleic acid molecule in the sample) typically allows for a unique identity to be assigned to a particular molecule. The length of the individual sequence reads, or the number of base pairs, is also used as needed to assign a unique identity to a given molecule. As described herein, fragments derived from a single strand of nucleic acid may be assigned unique identification information, thereby allowing subsequent identification of fragments derived from the parental and / or complementary strand.

[0076] In some embodiments, molecular barcodes are introduced to molecules in a sample at an expected ratio of a set of identifiers (e.g., a combination of unique or non-unique molecular barcodes). In one exemplary format, about 2 to about 1,000,000 different molecular barcodes, about 5 to about 150 different molecular barcodes, or about 20 to about 50 different molecular barcodes are used that are ligated to both ends of a target molecule. Alternatively, about 25 to about 1,000,000 different molecular barcodes may be used. For example, 20 to 50 molecular barcodes and 20 to 50 molecular barcodes may be used such that both ends of a target molecule are tagged with one of the 20 to 50 different molecular barcodes. Such a number of identifiers is typically sufficient to ensure a high probability (e.g., at least 94%, 99.5%, 99.99%, or 99.999%) that different molecules with the same start and end points will be tagged with different combinations of identifiers. In some embodiments, about 80%, about 90%, about 95%, or about 99% of the molecules have the same combination of molecular barcodes.

[0077] In some embodiments, the assignment of unique or non-unique molecular barcodes in the reactions is performed using methods and systems described, for example, in U.S. Patent Application Nos. 20010053519, 20030152490, and 20110160078, as well as U.S. Patent Nos. 6,582,908, 7,537,898, 9,598,731, and 9,902,992, each of which is incorporated herein by reference in its entirety. Alternatively, in some embodiments, different nucleic acid molecules of a sample can be distinguished using only intrinsic sequence information (e.g., start and / or end positions, subsequences at one or both ends of the sequence, and / or length). C. Nucleic Acid Amplification

[0078] Sample nucleic acids flanked by adaptors are typically amplified by PCR and other amplification methods using nucleic acid primers that bind to primer binding sites in the adaptors flanking the DNA molecules to be amplified. In some embodiments, the amplification method involves cycles of extension, denaturation, and annealing by temperature cycling, or may be isothermal, for example, as in the case of transcriptional amplification. Other exemplary amplification methods that are optionally utilized include ligase chain reaction, strand displacement amplification, nucleic acid sequence-based amplification, and self-sustaining sequence-based replication, among other approaches.

[0079] To introduce a molecular barcode and / or a sample index into a nucleic acid molecule using conventional nucleic acid amplification methods, one or more iterations of an amplification cycle are generally applied. The amplification is typically performed in one or more reaction mixtures. The molecular barcode and the sample index are optionally introduced simultaneously or in any sequential order. In other embodiments, the molecular barcode and the sample index are introduced before and / or after performing a sequence capture step. In some embodiments, only the molecular barcode is introduced before probe capture, and the sample index is introduced after performing a sequence capture step. In certain embodiments, both the molecular barcode and the sample index are introduced before performing a probe-based capture step. In some embodiments, the sample index is introduced after performing a sequence capture step. Typically, a sequence capture protocol involves introducing a single-stranded nucleic acid molecule complementary to a target nucleic acid sequence (e.g., a coding sequence of a genomic region), where mutations in such a region are associated with a cancer type. Typically, the amplification reaction produces a plurality of non-uniquely or uniquely tagged nucleic acid amplicons comprising a molecular barcode and a sample index, ranging in size from about 200 nucleotides (nt) to about 700 nt, 250 nt to about 350 nt, or about 320 nt to about 550 nt. In some embodiments, the size of the amplicon is about 300 nt. In some embodiments, the size of the amplicon is about 500 nt. D. Nucleic Acid Enrichment

[0080] In some embodiments, sequences are enriched prior to sequencing the nucleic acid. Enrichment is performed for a specific target region or non-specifically ("target sequence") as required. In some embodiments, target regions of interest may be enriched with nucleic acid capture probes ("baits") selected for a panel of one or more bait sets using a differential tiling and capture scheme. In a differential tiling and capture scheme, typically, different relative concentrations of bait sets are used to differentially tile (e.g., at different "resolutions") across the genomic regions associated with the baits, and a set of constraints (e.g., sequencer constraints such as sequencing load, utilization of each bait, etc.) are applied to capture the target nucleic acid at the desired stage of downstream sequencing. These target genomic regions of interest include natural or synthetic nucleotide sequences of nucleic acid constructs as required. In some embodiments, biotin-labeled beads with probes for one or more regions of interest can be used to capture the target sequences, followed by amplification of these regions as required to enrich for the regions of interest.

[0081] Sequence capture typically involves the use of oligonucleotide probes that hybridize to target nucleic acid sequences. In certain embodiments, the probe set strategy involves tiling probes across the region of interest. Such probes can be, for example, about 60 to about 120 nucleotides in length. The depth of the set can be about 2x, 3x, 4x, 5x, 6x, 8x, 9x, 10x, 15x, 20x, 50x, or greater than 50x. In general, the effectiveness of sequence capture depends in part on the length of the sequence in the target molecule that is complementary (or nearly complementary) to the sequence of the probe. E. Nucleic Acid Sequencing

[0082] The sample nucleic acid, pre-amplified or unamplified (optionally flanked by adapters), is generally subjected to sequencing.Sequencing methods or optionally commercially available sequencing formats can be used, for example, Sanger sequencing, high throughput sequencing, pyrosequencing, sequencing by synthesis, single molecule sequencing, nanopore-based sequencing, semiconductor sequencing, sequencing by ligation, sequencing by hybridization, RNA-Seq (Illumina), Digital Gene Expression (Helicos), next generation sequencing (NGS), single molecule sequencing by synthesis (SMSS) (Helicos), massively parallel sequencing, Clonal Single Molecule Array (Solexa), shotgun sequencing, Ion Torrent, Oxford Nanopore, Roche Genia, Maxim-Gilbert sequencing, primer walking, sequencing using PacBio, SOLiD, Ion Torrent or Nanopore platform. Sequencing reactions can be performed in a variety of sample processing units, which may include multi-lane, multi-channel, multi-well, or other means of processing multiple sample sets substantially simultaneously. Sample processing units may also include multiple sample chambers to allow multiple runs to be processed simultaneously.

[0083] The sequencing reaction may be performed on one or more nucleic acid fragment types or regions known to contain markers for cancer or other diseases. The sequencing reaction may also be performed on any nucleic acid fragment present in the sample. The sequencing reaction may be performed on at least about 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 99%, 99.9% or 100% of the genome. In other cases, the sequencing reaction may be performed on less than about 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 99%, 99.9% or 100% of the genome.

[0084] Simultaneous sequencing reaction can be carried out using multiplex sequencing technology.In some embodiments, cell-free polynucleotide is sequenced at least about 1000 times, 2000 times, 3000 times, 4000 times, 5000 times, 6000 times, 7000 times, 8000 times, 9000 times, 10000 times, 50000 times, or 100,000 times sequencing reaction.In other embodiments, cell-free polynucleotide is sequenced less than about 1000 times, less than 2000 times, less than 3000 times, less than 4000 times, less than 5000 times, less than 6000 times, less than 7000 times, less than 8000 times, less than 9000 times, less than 10000 times, less than 50000 times, or less than 100,000 times sequencing reaction.Sequencing reaction is typically carried out consecutively or simultaneously. Subsequent data analysis is generally carried out on all or part of sequencing reactions.In some embodiments, data analysis is carried out on at least about 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 50000 or 100,000 sequencing reactions.In other embodiments, data analysis can be carried out on less than about 1000, less than 2000, less than 3000, less than 4000, less than 5000, less than 6000, less than 7000, less than 8000, less than 9000, less than 10000, less than 50000 or less than 100,000 sequencing reactions. An exemplary read depth is about 1,000 to about 50,000 reads per locus (base position).

[0085] In some embodiments, the nucleic acid population is prepared for sequencing by enzymatically forming blunt ends on double-stranded nucleic acids with single-stranded overhangs at one or both ends. In these embodiments, the nucleic acid population is typically treated with an enzyme having 5'-3' DNA polymerase activity and 3'-5' exonuclease activity in the presence of nucleotides (e.g., A, C, G, and T or U), which may be present in a readily incorporated form, e.g., multiple nucleoside triphosphate (dNTP) form. Exemplary enzymes or optional catalytic fragments include Klenow large fragment and T4 polymerase. In the 5' overhang, the enzyme typically extends the recessed 3' end on the opposite strand until it is the same length as the 5' end, generating a blunt end. In the 3' overhang, the enzyme typically digests from the 3' end to the 5' end of the opposite strand, or optionally beyond. If the digestion proceeds beyond the 5' end of the opposite strand, the gap can be filled by an enzyme with the same polymerase activity as used for the 5' overhang. The formation of blunt ends in double-stranded nucleic acids facilitates, for example, adapter binding and subsequent amplification.

[0086] In some embodiments, the population of nucleic acids is further processed, such as converting single stranded nucleic acids to double stranded and / or converting RNA to DNA. These forms of nucleic acid are also optionally ligated to adaptors and amplified.

[0087] The nucleic acid that has been subjected to the above blunt end formation process, with or without prior amplification, and optionally other nucleic acids in the sample, can be sequenced to generate sequenced nucleic acid. Sequenced nucleic acid can refer to either the sequence (i.e., sequence information) of a nucleic acid, or the nucleic acid that has been sequenced. Sequencing can be performed to generate sequence data of individual nucleic acid molecules in a sample directly or indirectly from the consensus sequence of the amplification products of individual nucleic acid molecules in a sample.

[0088] In some embodiments, double-stranded nucleic acids with single-stranded overhangs in a sample are bound at both ends to adapters containing molecular barcodes after blunt end formation, and the nucleic acid sequence and the molecular barcode introduced by the adapter are determined by sequencing. Blunt-ended DNA molecules are optionally ligated to the blunt ends of at least partially double-stranded adapters (e.g., Y-shaped or bell-shaped adapters). Alternatively, the blunt ends of the sample nucleic acid and adapters may be overhanging with complementary nucleotides (e.g., for sticky end ligation) to facilitate ligation.

[0089] The nucleic acid sample is typically contacted with a sufficient number of adapters to make it unlikely that any two copies of the same nucleic acid will receive the same combination of adapter barcodes (i.e., molecular barcodes) from adapters attached to both ends. Using adapters in this manner allows for the identification of a family of nucleic acid sequences that have the same start and end points on the reference nucleic acid and are linked to the same combination of molecular barcodes. Such a family represents the sequence of the amplification product of the nucleic acid in the sample before amplification. The sequences of the family members, modified by blunt end formation and adapter attachment, can be compiled to derive a consensus nucleotide or a complete consensus sequence of the nucleic acid molecules in the original sample. In other words, a nucleotide occupying a particular position of the nucleic acid in the sample is determined to be the consensus of the nucleotides occupying the corresponding position in the family member sequences. A family can include sequences of one or both strands of a double-stranded nucleic acid. When a family member includes sequences of both strands from a double-stranded nucleic acid, the sequence of one strand is converted to its complementary sequence in order to compile the entire sequence to derive a consensus nucleotide or consensus sequence. Some families contain only a single member sequence, in which case this sequence may be interpreted as the sequence of the nucleic acid in the sample before amplification, or alternatively, families with only a single member sequence may be excluded from subsequent analysis.

[0090] Nucleotide variations in the sequenced nucleic acid may be determined by comparing the sequenced nucleic acid with a reference sequence. The reference sequence is often a known sequence, such as a known whole genome sequence or partial genome sequence from a subject (e.g., a whole genome sequence of a human subject). The reference sequence may be, for example, hG19 or hG38. The sequenced nucleic acid may represent a sequence determined directly for the nucleic acid in the sample, or a consensus of sequences of amplification products of such nucleic acid, as described above. Comparison may be performed at one or more designated positions of interest in the reference sequence. A subset of sequenced nucleic acids may be identified, including positions that correspond to the designated positions of the reference sequence when each sequence is maximally aligned. Within such a subset, it may be determined which, if any, sequenced nucleic acids contain nucleotide variations at the designated positions, as well as, if necessary, which, if any, sequenced nucleic acids contain reference nucleotides (i.e., the same as in the reference sequence). If the number of sequenced nucleic acids in the subset containing the nucleotide variant exceeds a selected threshold, the mutated nucleotide can be called at the designated position. The threshold can be a simple number, e.g., at least 1, 2, 3, 4, 5, 6, 7, 9, or 10 sequenced nucleic acids in the subset containing the nucleotide variant, among other possibilities, or a ratio, e.g., at least 0.5, 1, 2, 3, 4, 5, 10, 15, or 20 sequenced nucleic acids in the subset containing the nucleotide variant. The comparison can be repeated for any designated position of interest in the reference sequence. Optionally, the comparison can be performed for designated positions that occupy at least about 20, 100, 200, or 300 consecutive positions on the reference sequence, e.g., about 20 to 500, or about 50 to 300 consecutive positions.

[0091] Further details regarding nucleic acid sequencing, including the formats and applications described herein, can be found in, for example, Levy et al., Annual Review of Genomics and Human Genetics, 17:95-115 (2016); Liu et al., J. of Biomedicine and Biotechnology, Volume 2012, Article ID 251364:1-11 (2012); Voelkerding et al., Clinical Chem., 55:641-658 (2009); MacLean et al., Nature Rev. Microbiol., 7:287-296 (2009); Astier et al., J Am Chem. Soc., 128(5):1705-10(2006), U.S. Patent No. 6,210,891, U.S. Patent No. 6,258,568, U.S. Patent No. 6,833,246, U.S. Patent No. 7,115,400, U.S. Patent No. 6,969,488, U.S. Patent No. 5,912,148, U.S. Patent No. 6,130,073, U.S. Patent No. 7,169,560, U.S. Patent No. 7,282,337, U.S. Patent No. 7,482,120 Nos. 7,501,245, 6,818,395, 6,911,345, 7,501,245, 7,329,492, 7,170,050, 7,302,146, 7,313,308, and 7,476,503, each of which is incorporated by reference in its entirety. F. Analysis

[0092] Sequencing according to embodiments of the present disclosure generates multiple reads. The reads of the present invention generally comprise sequences of nucleotide data less than about 150 bases in length, or less than about 90 bases in length. In certain embodiments, the reads are about 80 to about 90 bases in length, such as about 85 bases in length. In some embodiments, the methods of the present invention are applied to very short reads, i.e., reads less than about 50 or about 30 bases in length. The sequence read data may include sequence data as well as meta information. The sequence read data may be stored in any suitable file format, including, for example, a VCF file, a FASTA file, or a FASTQ file.

[0093] FASTA was originally a computer program for searching sequence databases, and the name FASTA has come to refer to the standard file format. See Pearson and Lipman, 1988, Improved tools for biological sequence comparison, PNAS 85:2444-2448. A FASTA format sequence begins with a line of description, followed by lines of sequence data. The description line is separated from the sequence data by a greater than (">") symbol on the first line. The word following the ">" symbol is the sequence identifier, and the remainder of the line is the description (both are optional). There should be no space between the ">" symbol and the first character of the identifier. It is recommended that all lines of text be less than 80 characters long. The sequence ends when another line begins with ">", indicating the beginning of another sequence.

[0094] The FASTQ format is a text-based format for storing both biological sequences (usually nucleotide sequences) and their corresponding quality scores. The FASTQ format is similar to the FASTA format, but includes a quality score following the sequence data. Both the sequence and the quality score are encoded with a single ASCII character for brevity. The FASTQ format is the de facto standard for storing the output of high-throughput sequencing instruments such as the Illumina Genome Analyzer, as described, for example, in Cock et al. ("The Sanger FASTQ file format for sequences with quality scores, and the Solexa / Illumina FASTQ variants," Nucleic acids Res 38(6):1767-1771, 2009), which is incorporated herein by reference in its entirety.

[0095] For FASTA and FASTQ files, the meta information includes a description line, but not a sequence data line. In some embodiments, for FASTQ files, the meta information includes a quality score. For FASTA and FASTQ files, the sequence data begins after the description line and is typically represented using some subset of the IUPAC ambiguity code, with "-" added as needed. In a preferred embodiment, the sequence data will use the letters A, T, C, G, and N, and will include "-" or U as needed (e.g., to indicate a gap or uracil).

[0096] In some embodiments, at least one master sequence read file and output file are stored as plain text files (e.g., using encoding such as ASCII; ISO / IEC646; EBCDIC; UTF-8; or UTF-16). The computer system provided herein may include a text editor program capable of opening plain text files. A text editor program may refer to a computer program that can display the contents of a text file (e.g., a plain text file) on a computer screen and allow a human to edit the text (e.g., using a monitor, keyboard, and mouse). Exemplary text editors include, but are not limited to, Microsoft Word, emacs, pico, vi, BBEdit, and TextWrangler. Preferably, the text editor program can display the plain text file on a computer screen and present the meta information and sequence reads in a human-readable format (e.g., using alphanumeric characters as humans would write, rather than binary coded).

[0097] Although the methods have been discussed with reference to FASTA or FASTQ files, the methods and systems of the present invention may be used to compress any suitable sequence file format, including, for example, files in the Variant Call Format (VCF) format. A typical VCF file will include a header section and a data section. The header includes any number of meta-information lines, each of which begins with the characters "##", and is separated by tabs from field definition lines, each of which begins with a single "#" character. The field definition lines specify eight required columns, and the body section includes lines of data that make up the columns defined in the field definition lines. The VCF format is described by Danecek et al. ("The variant call format and VCFtools", Bioinformatics 27(15):2156-2158, 2011), which is incorporated herein by reference in its entirety. The header section may be treated as meta-information that is written to the compressed file, and the data section may be treated as the lines, each of which will be stored in the master file only if it is unique.

[0098] Certain embodiments of the present invention provide assembly of sequence reads. In assembly by alignment, for example, reads are aligned with each other or aligned to a reference. By aligning each read and then aligning to a reference genome, all reads are positioned in relation to each other to create an assembly. In addition, alignment or mapping of sequence reads to a reference sequence can also be used to identify variant sequences within sequence reads. Identification of variant sequences can be used in combination with the methods and systems described herein to further aid in the diagnosis or prognosis of a disease or condition, or to guide treatment decisions.

[0099] In some embodiments, any or all of these steps are automated. Alternatively, the method of the present invention is incorporated, in whole or in part, into one or more dedicated programs, for example, each written in a compiled language such as C++ as needed, then compiled and delivered as a binary. The method of the present invention may be implemented, in whole or in part, as a module within an existing sequence analysis platform, or functionally executed within an existing sequence analysis platform. In certain embodiments, the method of the present invention includes multiple steps that are all performed automatically in response to a single initiating cue (e.g., a triggering event or combination of events caused by a human action, another computer program, or a machine). Thus, the present invention provides a method in which any of the steps, or any combination of the steps, occurs automatically in response to a cue. "Automatically" generally means without human input, influence, or interaction (i.e., only in response to an initial or previous cueing human action).

[0100] The system also includes various output formats, including accurate and sensitive interpretation of the nucleic acid of interest. The output of the search can be provided in the format of a computer file. In certain embodiments, the output is a FASTA file, a FASTQ file, or a VCF file. The output can be processed to generate a text file or an XML file that includes sequence data, such as the sequence of the nucleic acid aligned to the sequence of the reference genome. In other embodiments, the processing results in an output that includes a coordinate or string that describes one or more mutations in the nucleic acid of interest compared to the reference genome. Sequence alignments may include the Simple UnGapped Alignment Report (SUGAR), the Verbose Useful Labeled Gapped Alignment Report (VULGAR), and the Compact Idiosyncratic Gapped Alignment Report (CIGAR) (Ning et al., Genome Research 11(10):1725-9, 2001, which are incorporated by reference in their entireties. These strings are implemented, for example, in the Exonerate sequence alignment software by the European Bioinformatics Institute (Hinxton, UK).

[0101] In some embodiments, a sequence alignment (e.g., a sequence alignment map (SAM) or binary alignment map (BAM) file) is generated that includes a CIGAR column (the SAM format is described, for example, by Li et al., "The Sequence Alignment / Map format and SAMtools," Bioinformatics, 25(16):2078-9, 2009, which is incorporated herein by reference in its entirety). In some embodiments, CIGAR displays or includes gapped alignments, one per line. CIGAR is a condensed pairwise alignment format reported as the CIGAR column. The CIGAR column is useful for displaying long (e.g., genomic) pairwise alignments. The CIGAR column is used in the SAM format to display the alignment of reads to a reference genome sequence.

[0102] The CIGAR string follows the established motif. Each letter is preceded by a number to indicate the number of bases in the event. The letters used may include M, I, D, N, and S (M=match; I=insertion; D=deletion; N=gap; S=substitution). The CIGAR string describes the sequence of matches / mismatches and deletions (or gaps). For example, the CIGAR string 2MD3M2D2M would mean an alignment containing 2 matches, 1 deletion (the number 1 is omitted to save space), 3 matches, 2 deletions, and 2 matches.

[0103] In some embodiments, the results of the systems and methods disclosed herein are used as input to generate a report. The report may be in paper or electronic format. For example, information about the allelic imbalance status of the sample determined by the methods or systems disclosed herein may be displayed in such a report. Alternatively, or in addition, information about the presence or absence of contamination in the sample as determined by the methods or systems disclosed herein may be displayed in such a report. The methods or systems disclosed herein may further include a step of transmitting such a report to a third party (e.g., the subject from whom the sample originated, or a medical professional).

[0104] Various steps of the methods disclosed herein or performed by the systems disclosed herein may be performed at the same or different times, in the same or different geographic locations (e.g., countries), and / or by the same or different persons.

[0105] The method may also be used to determine or monitor the effectiveness of a treatment by the relative amounts of therapeutic nucleic acid constructs at different time points.

[0106] FIG. 3 illustrates a computer system 301 that is programmed or otherwise configured to perform the methods provided herein.

[0107] The computer system 301 may be programmed or otherwise configured to execute the architecture for training neural networks using biological sequences, storage, and molecular phenotypes. The computer system 301 can control various aspects of the present disclosure, such as, for example, (a) sequencing a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying, for at least a portion of the plurality of aligned sequence reads, germline variants present in the sample at a mutant allele fraction (MAF), where each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants identified in (c) that are between a plurality of discrete ranges of MAF values; and (e) detecting the allelic imbalance in the sample based on a predetermined criterion by filtering the set of germline variants identified in (c) based on at least the quantitative measure of (d). The computer system 301 may be a user's electronic device or may be a computer system located remotely from the electronic device, which may be a mobile electronic device.

[0108] The computer system 301 includes a central processing unit (CPU, or, as used herein, "processor" and "computer processor") 305, which may be a single-core or multi-core processor, or multiple processors for parallel processing. The computer system 301 also includes memory or storage 310 (e.g., random access memory, read-only memory, flash memory), an electronic storage unit 315 (e.g., hard disk), a communication interface 320 (e.g., network adapter) for communicating with one or more other systems, and peripheral devices 325 (e.g., cache, other memory, data storage, and / or electronic display adapter, etc.). The memory 310, the storage unit 315, the interface 320, and the peripheral devices 325 are in communication with the CPU 305 through a communication bus (solid line) (e.g., motherboard, etc.). The storage unit 315 may be a data storage unit (or data repository) for storing data. The computer system 301 may be operatively connected to a computer network ("network") 330 with the aid of the communication interface 320. Network 330 may be the Internet, an Internet and / or an extranet, or an intranet and / or an extranet in communication with the Internet. Network 330 is, in some cases, a telecommunications and / or data network. Network 330 may include one or more computer servers, which may enable distributed computing such as cloud computing. Network 330 may, in some cases, with the aid of computer system 301, implement a peer-to-peer (P2P) network, which may enable devices connected to computer system 301 to operate as clients or servers.

[0109] The CPU 305 may execute a sequence of machine-readable instructions, which may be embodied in a program or software. The instructions may be stored in a storage location, such as the memory 310. The instructions may be directed to the CPU 305, which may then be programmed or otherwise configured to execute the methods of the present disclosure. Examples of operations performed by the CPU 305 include fetch, decode, execute, and writeback.

[0110] The CPU 305 may be part of a circuit (e.g., an integrated circuit). One or more other components of the system 301 may also be included in the circuit. In some cases, the circuit is an application specific integrated circuit (ASIC).

[0111] The storage unit 315 can store files (e.g., drivers, libraries, and saved programs). The storage unit 315 can store user data (e.g., user preferences and user programs). The computer system 301 may, in some cases, include one or more additional data storage units external to the computer system 301 (e.g., located on a remote server in communication with the computer system 301 through an intranet or the Internet).

[0112] Computer system 301 can communicate with one or more remote computer systems through network 330. For example, computer system 301 can communicate with a user's remote computer system. Examples of remote computer systems include a personal computer (e.g., a portable PC), a slate or tablet PC (e.g., Apple® iPad®, Samsung® Galaxy Tab), a phone, a smart phone (e.g., Apple® iPhone®, Android®-enabled device, Blackberry®), or a PDA (personal digital assistant). A user can access computer system 301 through network 330.

[0113] The methods described herein may be performed by machine (e.g., a computer processor) executable code stored in an electronic storage area of ​​the computer system 301 (e.g., in memory 310 or electronic storage unit 315). The machine executable or machine readable code may be provided in the form of software. In use, the code may be executed by the processor 305. In some cases, the code is read from the storage unit 315 and stored in memory 310 for ready access by the processor 305. In some circumstances, the electronic storage unit 315 may be eliminated and the machine executable instructions may be stored in memory 310.

[0114] The code may be pre-compiled and configured for use on a machine having a processor adapted to execute the code, or it may be compiled during run-time. The code may be provided in a programming language, which may be selected to enable the code to be executed in a pre-compiled or run-time compiled manner.

[0115] Aspects of the systems and methods provided herein, such as the computer system 301, may be embodied in programming. Various aspects of this technology may be considered to be "products" or "articles of manufacture," typically in the form of machine (or processor) executable code and / or associated data carried or embodied in some type of machine-readable medium. The machine executable code may be stored in an electronic storage unit, such as a memory (e.g., read-only memory, random access memory, flash memory), or a hard disk. A "storage" type medium may include any or all tangible memory (e.g., various semiconductor memories, tape drives, disk drives, etc.) of a computer, processor, etc., or their associated modules, which may provide non-transitory storage at any point in time for software programming. All or a portion of the software may be communicated at times over the Internet or various other long distance communication networks. Such communication may, for example, enable the software to be loaded from one computer or processor to another, for example, from a management server or host computer to the computer platform of an application server. Thus, another type of medium that may carry software elements includes light waves, radio waves, or electromagnetic waves (e.g., as used to physically interface between local devices through wired and optical landline networks, and by various wireless links). The physical elements that carry such waves (e.g., wired or wireless links, optical links, etc.) may also be considered media that carry software. As used herein, unless specifically limited to non-transitory tangible "storage" media, terms such as computer or machine "readable medium" refer to any medium that participates in providing instructions to a processor for execution.

[0116] Thus, the machine-readable medium (e.g., computer-executable code) can take many forms, including, but not limited to, tangible storage media, carrier wave media, or physical transmission media. Non-volatile storage media include, for example, optical or magnetic disks (e.g., any of the storage devices in any computer, etc.), which can be used to build, for example, the databases shown in the figures. Volatile storage media include dynamic memories, such as the main memory of such a computer platform. Tangible communication media include coaxial cables; copper wire and fiber optics, including the wires that make up a bus in a computer system. Carrier wave transmission media can be in the form of electric or electromagnetic signals, or acoustic or light waves, such as those generated during radio frequency (RF) and infrared (IR) data communications. Thus, common forms of computer readable media include, for example, a floppy disk, a flexible disk, a hard disk, a magnetic tape, any other magnetic medium; a CD-ROM, a DVD or a DVD-ROM, any other optical medium; punch cards, paper tape, any other physical storage medium having a pattern of holes; a RAM, a ROM, a PROM and an EPROM, a FLASH-EPROM, any other memory chip or cartridge; a carrier wave carrying data or instructions; a cable or link carrying such a carrier wave; or any other medium from which a computer can read program code and / or data. Many of these forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to a processor for execution.

[0117] Computer system 301 may include or be in communication with electronic display 335 that includes user interface (UI) 340. Examples of UIs include, but are not limited to, graphical user interfaces (GUIs) and web-based user interfaces.

[0118] The disclosed method and system may be implemented as one or more algorithms. The algorithm may be implemented as software when executed by the central processing unit 305. The algorithm may, for example, (a) align at least a portion of a plurality of sequence reads from a sequencer to a reference sequence to generate a plurality of aligned sequence reads; (b) for at least a portion of the plurality of aligned sequence reads, identify a set of germline variants in the sample by identifying germline variants present in the sample at a mutant allele fraction (MAF) or minor allele frequency, where each germline variant in the set of germline variants has a corresponding MAF or minor allele frequency value; (c) determine a quantitative measure of the set of germline variants identified in (b), whose MAF or minor allele frequency values ​​are between a plurality of separate ranges; and (d) detect the allelic imbalance in the sample based on a predetermined criterion by filtering the set of germline variants identified in (b) based on at least the quantitative measure of (c).

[0119] Although the above description has been described with reference to specific embodiments, these specific embodiments are exemplary only and not limiting, and the concepts illustrated in the examples may be applied to other examples and implementations. EXAMPLES

[0120] Example 1: Discrimination between samples with allelic imbalance and samples with contamination When samples are assayed using conventional cell-free DNA analysis methods, any sample with more than two germline variants present at MAFs within the range of somatic MAFs (which may be less than about 15%) requires manual review to assess whether the sample is in a "possibly contaminated" state. Such an approach flags different samples that contain multiple such germline variants, such as (1) samples with assay-level contamination, (2) samples with a second genome (e.g., from a graft, transfusion, or fetus), and (3) samples that show allelic imbalance as a result of loss of heterozygosity (LoH). Furthermore, when samples are analyzed by conventional cfDNA assay methods, samples in such cases cannot be identified. For example, both samples with a second genome and samples that show allelic imbalance as a result of LoH would be erroneously considered to have assay-level contamination, thereby requiring repeat sample assay for confirmation purposes. Therefore, this approach can increase assay turnaround time and costs by overcalling contaminating samples and thus requiring re-assaying of samples that in fact have allelic imbalance rather than contamination.

[0121] In the case of samples without copy number variation or alteration, somatic variants may be measured directly from the tumor source. However, if copy number variation or alteration exists in the sample, if such variation includes a germline variant that causes LoH, MAF measurement may be distorted (e.g., MAF measurement may shift and be off by 50%), thereby inducing false positive contamination evaluation and re-assay analysis of the sample. Such allelic imbalance may be seen in patients with CNV resulting from LoH (which is related to copy number) or CN-LoH (copy-neutral LoH) (e.g., resulting from gene exchange between two chromosome arms such that chromosomal information is kept constant). For example, detection of such LoH (which indicates that a gene loses its allele (e.g., loses gene function)) may have important implications for treatment selection, monitoring, and evaluation.

[0122] Using the disclosed method and system, a sample containing cell-free DNA molecules is assayed, and the results are analyzed using a decision tree to distinguish between samples with allele imbalance and samples with contamination. Figure 2 shows an example of a workflow 200 for detecting the presence or absence of allele imbalance or contamination in a cell-free DNA sample. The workflow 200 can include determining (as in operation 202) a quantitative measurement of the germline variant for the cell-free DNA molecules of the sample that is between a plurality of separate ranges of MAF values. Next, the workflow 200 can include determining (as in operation 204) a value of max_CNV (maximum CNV level of all genes measured across the samples), min_CNV (minimum CNV level of all genes measured across the samples), or frac_diploid (diploid gene fraction) at the sample level. Next, the workflow 200 may include determining whether a first criterion is met (as in operation 206), e.g., whether the measurements of germline variants and the values ​​of max_CNV, min_CNV, or frac_diploid meet a particular criterion. If the determination in operation 206 is "yes" (i.e., the first criterion is positive), the workflow proceeds to operation 208; alternatively, if the determination in operation 206 is "no" (i.e., the first criterion is negative), the workflow proceeds to operation 212. Next, the workflow 200 may include determining whether a second criterion is met (as in operation 208), e.g., whether the allelic imbalance candidate (e.g., the cfDNA sample being analyzed to detect the presence or absence of allelic imbalance or contamination) has a germline variant that meets a low MAF criterion. If the determination at operation 208 is "yes" (i.e., the second criterion is positive), the workflow proceeds to operation 210; alternatively, if the determination at operation 208 is "no" (i.e., the second criterion is negative), the workflow proceeds to operation 212. The workflow 200 may then include, for example, generating an output or indication that the sample has allelic imbalance (as in operation 210).Alternatively, workflow 200 may include generating an output or indication that the sample has contamination (e.g., assay-level contamination or contamination with a second genome) (as in operation 212).

[0123] In some embodiments, all the criteria in the decision tree are applied. The first criterion in the decision tree is applied to identify samples that may be contaminated. The second criterion in the decision tree is applied to evaluate the number of germline variants within any of a number of separate ranges (e.g., windows) of MAF values, including about 3% to about 40% MAF and about 60% to about 97% MAF. If the number is large and supported by copy number, such samples may have allelic imbalance. The third criterion in the decision tree is applied to detect extreme cases where a very large number of copy number changes may result in germline variants with MAFs less than about 3%.

[0124] A first set of over 20,000 clinical samples is processed using a 73-gene cell-free DNA (cfDNA) next-generation sequencing (NGS) panel (Guardant Health, Redwood City, CA). From this first set, a training set of 224 samples is selected that have been manually re-assayed to identify allelically imbalanced and contaminated samples. For example, if manual re-assaying results in a given sample no longer displaying a possible contaminant mark, the first run can be identified as indeed likely contaminated. In addition, some patients are contacted to confirm the second genomic status (e.g., transplant, transfusion, or fetal). The contamination status of each of the 224-sample training set is manually reviewed. From this first set, a test set of 2,300 samples is selected, of which 37 samples were originally flagged as possible contaminants.

[0125] In some embodiments, the cell-free DNA assay produces a plurality of genetic variants, including germline variants and somatic variants. Among these plurality of genetic variants, the germline or somatic status of a given genetic variant may be determined (e.g., identified) using a beta binomial distribution model that estimates the mean and variance of MAF values ​​for common germline SNPs located in the vicinity of the candidate variant under consideration. Further details regarding the beta binomial distribution, which is optionally adapted for use in the implementation model of the methods and related aspects disclosed herein, are also described, for example, in International Patent Application No. PCT / US2018 / 052087, filed September 20, 2018, which is incorporated herein by reference in its entirety.

[0126] To identify potentially contaminated samples, a first criterion is applied to assess whether a given sample has more than two common germline single nucleotide polymorphisms (SNPs) with a variant allele fraction (MAF) of less than 15%. If this first criterion is met, a second criterion is applied to assess whether the sample (a) has more than 21 germline variants within any of a plurality of separate ranges (e.g., windows) of MAF values, including about 3% to about 40% MAF and about 60% to about 97% MAF, and (b) the genes within these separate ranges in the sample have a maximum CNV level greater than 0.22, a minimum CNV level less than -0.14, or a diploid gene fraction (e.g., diploid fraction) less than 0.7. The aforementioned thresholds may be determined using a training dataset of a large number of samples (e.g., about 50 samples, about 100 samples, about 150 samples, about 200 samples, about 250 samples) where the contamination / allelic imbalance status of these samples is known and / or these ranges provide maximum accuracy.

[0127] The second criterion may include a quantitative measurement indicative of copy number (e.g., resulting from allelic imbalance or loss of heterozygosity). A quantitative measurement indicative of copy number may include a total of measurements of genomic disruptions (e.g., estimates of the total of copy number changes) (e.g., may be expressed as CNV or diploid fraction); a quantitative measurement obtained by binning by chromosome or chromosome arm; or a quantitative measurement obtained by observing disruptions throughout the genome, measuring the relative amount of distortion in each disruption, and predicting from such measurements the likelihood that another gene on the same chromosome may be altered to the same extent (e.g., as a result of copy-neutral LoH (CN-LoH)). The second criterion evaluates whether there is evidence that copy number alterations may move germline variants into a wider MAF window (e.g., about 3% to about 40% or about 60% to about 97%).

[0128] If this second criterion is met, a third criterion is used to assess whether the sample has either (a) no germline variant with a MAF less than about 3%, or (b) a germline variant with a MAF less than about 3% and an absolute copy number average greater than about 10 (e.g., copy number average greater than about 10 or less than about -10) for the same germline variant. The third criterion assesses whether the extreme case occurs where a very large number of copy number changes can result in a germline variant with a MAF less than about 3%. If the third criterion is met, the sample is identified as having allelic imbalance (e.g., an allelic imbalance sample). If the third criterion is not met, the sample is identified as having contamination (e.g., a truly contaminated sample).

[0129] The performance of the method for detecting contaminated samples (e.g., samples without allelic imbalance) is shown below for a training dataset of 224 samples (selected from a larger set of at least 20,000 different samples) (Table 1), and a testing dataset of at least 2,300 different samples (Table 2).

[0130] Table 1 [Table 1]

[0131] Table 2 [Table 2]

[0132] By applying the methods disclosed herein to distinguish between samples with allelic imbalance and samples with contamination, the overcall rate of the cell-free DNA assay is reduced by 20% while maintaining perfect sensitivity of 100% in detecting samples with true contamination.

[0133] If the liquid biopsy assay changes (e.g., in sequencing depth and panel of common SNPs), the methods and systems of the present disclosure may be retrained as necessary to obtain a reasonable set of thresholds (e.g., for application of one or more criteria of a decision tree to distinguish between samples with allelic imbalance and samples with contamination). Example 2: Detection of allele-specific loss of heterozygosity (LoH) in cell-free DNA (cfDNA)

[0134] Loss of heterozygosity (LoH) is a common feature in tumor biology and can occur frequently due to defects in Homologous Recombination Repair (HRR), resulting in uniparental deletions that manifest as LoH. In the absence of a driving force, allele loss is equally likely and therefore the proportion of retention and loss of a given allele in a population will be equal, although allele-specific loss (or retention) can occur.

[0135] A set of over 70,000 whole blood samples was obtained from patients with advanced solid tumors and assayed using a 73-gene cell-free DNA (cfDNA) next-generation sequencing (NGS) panel (Guardant Health, Redwood City, CA). By implementing the methods disclosed herein, the resulting ctDNA data (including observed allele frequencies and copy number variations) was analyzed using a database of tumor-associated variants to identify allele-specific losses.

[0136] Database analysis revealed that LoH often manifests as allelic imbalance, where the observed mutant allele fraction (MAF) of the retained allele is greater than 50% of the observed allele frequency and the observed mutant allele fraction (MAF) of the lost allele is less than 50% in an individual sample. This imbalance occurs because allele frequency is a relative measure, so the loss of one allele leads to a proportional increase in the amount of remaining alleles. Population analysis revealed that most allele losses are indiscriminate, but certain alleles are more likely to be retained or lost.

[0137] As an example, in a set of more than 90,000 whole blood samples analyzed, 56 variants of the BRCA1 gene were observed in one or more individual samples of this set, but for each variant, the MAF measured for the given variant in all individual samples that have the given variant is less than 50%, suggesting the possibility of allele-specific loss.For example, the BRCA1 P209L variant was observed in 9 individual samples of this set of more than 90,000 whole blood samples, and the MAF measured for the BRCA1 P209L variant in each of these 9 individual samples is less than 50%.Detection of allele-specific loss from ctDNA data provides insight into the underlying tumor biology and the selective pressure that leads to tumor evolution during the course of treatment.

[0138] Although preferred aspects of the invention have been shown and described herein, it will be apparent to those skilled in the art that the embodiments are presented by way of example only. It is not intended that the invention be limited by the specific examples shown herein. Although the invention has been described with reference to the foregoing specification, the description and illustration of the embodiments herein are not intended to be construed in a limiting sense. Numerous variations, changes, and substitutions will readily occur to those skilled in the art without departing from the invention. Moreover, it will be understood that all aspects of the invention are not limited to the specific depictions, configurations, or relative proportions shown herein, which depend upon a variety of conditions and variables. It should be understood that various alternatives to the embodiments of the invention described herein may be employed in the practice of the invention. Thus, it is intended that any such alternatives, modifications, variations, or equivalents are also encompassed by the present invention. It is intended that the following claims define the scope of the invention, and that methods and structures within the scope of the claims and their equivalents are encompassed within the scope of the claims.

Claims

1. 1. A method for detecting the presence or absence of allelic imbalance in a sample from a subject, comprising: (a) sequencing a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample to generate a plurality of sequence reads; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying, for at least a portion of the plurality of aligned sequence reads, germline variants present in the sample at a variant allele fraction (MAF), wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d) to detect the presence or absence of the allelic imbalance in the sample based on predefined criteria. Including, the quantitative measure of (d) comprises a number of sets of germline variants that fall between a plurality of discrete ranges of MAF values; (e) the detecting comprises detecting one or more quantitative measures indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, wherein the quantitative measures are selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average; The method of claim 1, wherein the predetermined criteria is one or more criteria selected from the group consisting of a quantitative measure of (d) greater than a predetermined germline variant threshold, a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and an average copy number for the same germline variant, and the same germline variant has a MAF of less than 3%.

2. 10. The method of claim 1, further comprising detecting the presence or absence of contamination or a second genome in the sample if the absence of the allelic imbalance is detected in the sample.

3. 3. The method of claim 1 or 2, wherein the set of germline variants comprises at least 1,000 different germline variants.

4. 4. The method of any one of claims 1 to 3, wherein the set of germline variants comprises germline variants selected from the group consisting of single nucleotide variants (SNVs), insertions or deletions (indels), and fusions.

5. The method according to any one of claims 1 to 4, wherein the sample is a body fluid sample selected from the group consisting of blood, plasma, serum, urine, saliva, mucosal secretions, sputum, stool, and tears.

6. The method of any one of claims 1 to 5, wherein the subject has a disease or disorder.

7. The method of claim 6, wherein the disease is cancer.

8. 8. The method of any one of claims 1 to 7, further comprising amplifying the cell-free DNA molecules prior to sequencing.

9. 10. The method of claim 8, further comprising selectively enriching said cell-free DNA molecules, or said amplified cell-free DNA molecules, for a set of loci prior to sequencing.

10. 10. The method of any one of claims 1-9, further comprising attaching one or more adapters comprising molecular barcodes to the cell-free DNA molecules prior to sequencing.

11. 11. The method of claim 10, wherein the one or more adaptors are randomly attached to both ends of the cell-free DNA molecule.

12. 11. The method of claim 10, wherein the cell-free DNA molecules are uniquely barcoded with a molecular barcode.

13. 11. The method of claim 10, wherein the cell-free DNA molecules are non-uniquely barcoded with a molecular barcode.

14. 11. The method of claim 10, wherein each molecular barcode comprises a predetermined or semi-random oligonucleotide sequence that, in combination with the diversity of molecules sequenced from a selected region, allows for the identification of a unique cell-free DNA molecule.

15. 15. The method of claim 14, wherein the selected region comprises germline variants found in COSMIC, The Cancer Genome Atlas (TCGA), or Exome Aggregation Consortium (ExAC).

16. 16. The method of any one of claims 1 to 15, wherein the plurality of discrete ranges of MAF values ​​comprises a first range of 3% to 40% and a second range of 60% to 97%.

17. 17. The method of claim 16, wherein the quantitative measurements of (d) include a set of a plurality of said germline variants that fall between a plurality of discrete ranges of MAF values.

18. 18. The method of claim 17, wherein the predetermined germline variant threshold is 21.

19. 19. The method of any one of claims 1 to 18, wherein the one or more quantitative measures indicative of the CNV or the diploid genes comprise two or more quantitative measures selected from the group consisting of: maximum CNV level across the samples, minimum CNV level across the samples, diploid gene fraction, and copy number average.

20. 20. The method of claim 19, wherein the one or more quantitative measurements indicative of the CNV or the diploid genes comprise three or more quantitative measurements selected from the group consisting of: maximum CNV level across the samples, minimum CNV level across the samples, diploid gene fraction, and copy number average.

21. 21. The method of claim 19 or 20, wherein the predetermined criteria comprises one or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than 3%.

22. 22. The method of any one of claims 1 to 21, wherein the predetermined criteria comprises two or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than 3%.

23. 23. The method of claim 22, wherein the predetermined criteria comprises three or more criteria selected from the group consisting of: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than 3%.

24. 24. The method of claim 23, wherein the predetermined criteria include the following: a maximum CNV level across the samples is greater than a predetermined maximum CNV threshold; a minimum CNV level across the samples is less than a predetermined minimum CNV threshold; a diploid gene fraction is less than a predetermined diploid fraction threshold; and an absolute value of the copy number average at the same germline variant is greater than a predetermined copy number average threshold and a MAF of the same germline variant is less than 3%.

25. 25. The method of any one of claims 1 to 24, wherein the predetermined criteria comprises one or more thresholds selected from the group consisting of a maximum CNV threshold of 0.22, a minimum CNV threshold of -0.14, a diploid fraction threshold of 0.7, and a copy number average threshold of 10.

26. 26. The method of claim 25, wherein the predetermined criteria comprises two or more thresholds selected from the group consisting of a maximum CNV threshold of 0.22, a minimum CNV threshold of -0.14, a diploid fraction threshold of 0.7, and a copy number average threshold of 10.

27. 27. The method of claim 26, wherein the predetermined criteria comprises three or more thresholds selected from the group consisting of a maximum CNV threshold of 0.22, a minimum CNV threshold of -0.14, a diploid fraction threshold of 0.7, and a copy number average threshold of 10.

28. 28. The method of claim 27, wherein the predetermined criteria include the following thresholds: a maximum CNV threshold of 0.22, a minimum CNV threshold of -0.14, a diploid fraction threshold of 0.7, and a copy number average threshold of 10.

29. 3. The method of claim 2, further comprising detecting the presence of said contaminant or said second genome in said sample with a positive predictive value (PPV) of at least 60%.

30. 3. The method of claim 2, further comprising detecting the absence of said contamination or said second genome in said sample with a negative predictive value (NPV) of at least 90%.

31. 3. The method of claim 2, further comprising detecting the presence of said contamination or said second genome in said sample with a sensitivity of at least 90%.

32. 3. The method of claim 2, further comprising detecting the presence of said contamination or said second genome in said sample with a sensitivity of at least 99%.

33. 3. The method of claim 2, further comprising detecting the absence of said contamination or said second genome in said sample with a specificity of at least 35%.

34. The germline variant is (i) determining total and mutant allele counts for nucleic acid variants from the cell-free (cf) DNA molecules; (ii) identifying associated variables of said nucleic acid variants from said cfDNA molecules; (iii) determining a quantitative value for said associated variable of said nucleic acid variant; (iv) generating a statistical model for the predicted number of germline variant alleles at the genomic locus of the nucleic acid variant; (v) generating a probability value for the nucleic acid variant based at least in part on the statistical model for the predicted germline variant allele count, the quantitative value for the associated variable of the nucleic acid variant, and at least one of the total allele count and the variant allele count for the nucleic acid variant; and (vi) classifying the nucleic acid variant as (1) of somatic origin if the p-value for the nucleic acid variant is less than a predetermined threshold, or (2) of germline origin if the p-value for the nucleic acid variant is equal to or greater than a predetermined threshold. The method of any one of claims 1 to 33, further comprising identifying by:

35. The method of any one of claims 1 to 34, further comprising detecting an allele-specific loss in the sample based on at least one of the set of germline variants identified as present in a given MAF in (c).

36. 36. The method of claim 35, wherein the allele-specific loss in the sample is detected based on the at least one of the set of germline variants being present in the sample from the subject at a MAF below 50%.

37. 37. The method of claim 36, wherein the allele-specific loss in the sample is detected based on the at least one of the set of germline variants being present in the sample from the subject, and in each of one or more samples from an additional one or more subjects, at a MAF below 50%.

38. 38. The method of any one of claims 35 to 37, wherein said at least one of said sets of germline variants is found in COSMIC, TCGA (The Cancer Genome Atlas), or ExAC (Exome Aggregation Consortium).

39. 39. The method of claim 38, wherein said at least one of said set of germline variants is a BRCA1 gene variant.

40. 40. The method of claim 39, wherein the BRCA1 gene variant is BRCA1 P209L.

41. The method of any one of claims 1 to 40, wherein at least part of the method is performed by a computer system.

42. 1. A system, which when executed by at least one electronic processor, comprises at least (a) obtaining a plurality of sequence reads corresponding to a plurality of cell-free deoxyribonucleic acid (DNA) molecules from a sample of a subject; (b) aligning at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying a set of germline variants in the sample by identifying, for at least a portion of the plurality of aligned sequence reads, germline variants present in the sample at a variant allele fraction (MAF), wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d) to detect the presence or absence of allelic imbalance in the sample based on predefined criteria. A controller including, or capable of accessing, a computer-readable medium including non-transitory computer-executable instructions for implementing the method of the present invention; the quantitative measure of (d) comprises a number of sets of germline variants that fall between a plurality of discrete ranges of MAF values; (e) the detecting comprises detecting one or more quantitative measures indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, wherein the quantitative measures are selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average; The system, wherein the predetermined criteria is one or more criteria selected from the group consisting of a quantitative measure of (d) greater than a predetermined germline variant threshold, a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and an average copy number for the same germline variant, and the same germline variant has a MAF of less than 3%.

43. 43. The system of Claim 42, further comprising a nucleic acid sequencer operably connected to the controller, wherein the nucleic acid sequencer is configured to process the plurality of cell-free DNA molecules from the sample to generate the plurality of sequence reads.

44. 44. The system of claim 42 or 43, further comprising: generating a report, wherein the non-transitory computer-executable instructions, when executed by at least one electronic processor,

45. 45. The system of claim 44, wherein the non-transitory computer-executable instructions, when executed by at least one electronic processor, further perform: communicating the report to a third party.

46. 1. A method for detecting the presence or absence of allelic imbalance in a sample from a subject, comprising: (a) accessing, by a computer system, a plurality of sequencing reads generated from a plurality of cell-free deoxyribonucleic acid (DNA) molecules from the sample; (b) aligning, by the computer system, at least a portion of the plurality of sequence reads to a reference sequence to generate a plurality of aligned sequence reads; (c) identifying, by the computer system, germline variants present in the sample at a variant allele fraction (MAF) for at least a portion of the plurality of aligned sequence reads, thereby identifying a set of germline variants in the sample, wherein each germline variant in the set of germline variants has a corresponding MAF value; (d) determining, by the computer system, a quantitative measure of the set of germline variants that fall between a plurality of discrete ranges of MAF values ​​identified in (c); and (e) detecting, by the computer system, the presence or absence of the allelic imbalance in the sample based on predetermined criteria by filtering the set of germline variants identified in (c) based on at least the quantitative measurements of (d). Including, the quantitative measure of (d) comprises a number of sets of germline variants that fall between a plurality of discrete ranges of MAF values; (e) the detecting comprises detecting one or more quantitative measures indicative of copy number variation (CNV) or diploid genes from the plurality of aligned sequence reads, wherein the quantitative measures are selected from the group consisting of a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and a copy number average; The method of claim 1, wherein the predetermined criteria is one or more criteria selected from the group consisting of a quantitative measure of (d) greater than a predetermined germline variant threshold, a maximum CNV level across the samples, a minimum CNV level across the samples, a diploid gene fraction, and an average copy number for the same germline variant, and the same germline variant has a MAF of less than 3%.

47. The method of any one of claims 1 to 41 and 46, further comprising generating a report.

48. 48. The method of claim 47, further comprising communicating the report to a third party.

Images