Multimodal analysis of circulating tumor nucleic acid molecules

Abstract

In an embodiment, a method for detecting the presence of ctDNA from cancer cells of a subject includes the steps of: (a) obtaining a sample of cell-free DNA from the subject; (b) subjecting the sample to library preparation to allow subsequent sequencing of the cell-free methylated DNA; (c) optionally adding a first amount of filler DNA to the sample, wherein at least a portion of the filler DNA is methylated, and then optionally further denaturing the sample; (d) capturing the cell-free methylated DNA using a binding agent selective for methylated polynucleotides; (e) sequencing the captured cell-free methylated DNA; (f) comparing the sequences of the captured cell-free methylated DNA with control cell-free methylated DNA sequences from healthy individuals and individuals with cancer; and (g) identifying the presence of DNA from cancer cells if there is a statistically significant similarity between one or more sequences of the captured cell-free methylated DNA and the cell-free methylated DNA sequence from the individual with cancer, wherein in at least one of the capturing, comparing, or identifying steps, the cell-free methylated DNA of interest is restricted to a subpopulation according to a fragment length metric.

Description

[Technical field]

[0001] cross reference This application claims the benefit of U.S. Provisional Patent Application No. 63 / 041,151, filed June 19, 2020, which is incorporated by reference in its entirety. [Background technology]

[0002] Circulating tumor DNA (ctDNA) is increasingly demonstrating its potential as a non-invasive tumor-specific biomarker for routine clinical use. ctDNA is primarily derived from tumor cells undergoing cell death and is released into the circulation of various body fluids, including blood. In most cancer patients, the majority of blood-derived cell-free DNA originates from peripheral blood leukocytes (PBLs). Thus, detection and quantification of ctDNA requires identification of tumor-derived genetic and epigenetic alterations. Furthermore, the fraction of ctDNA observed can range from less than 0.1% to 90% of the total cell-free DNA at diagnosis, depending on several factors including the primary site of the tumor and disease burden. ctDNA has provided non-invasive access to the molecular landscape and disease burden of tumors. Methods to detect ctDNA with high sensitivity, especially in subjects with low abundance of ctDNA, are needed.

[0003] Incorporation by Reference All publications, patents, and patent applications mentioned in this specification are herein incorporated by reference to the same extent as if each individual publication, patent, or patent application was specifically and individually indicated to be incorporated by reference. Summary of the Invention [Problem to be solved by the invention]

[0004] In an embodiment, there is provided a method for detecting the presence of ctDNA from a cancer cell of a subject, comprising: (a) obtaining a sample of cell-free DNA from a subject; (b) subjecting the sample to library preparation to allow subsequent sequencing of cell-free methylated DNA; (c) a first amount of filler DNA may be added to the sample, at least a portion of the filler DNA may be methylated, and then the sample may be further denatured; (d) capturing the cell-free methylated DNA with a binding agent selective for methylated polynucleotides; (e) sequencing the captured cell-free methylated DNA; (f) comparing the sequence of the captured cell-free methylated DNA with control cell-free methylated DNA sequences from healthy individuals and individuals with cancer; (g) identifying the presence of DNA from cancer cells if there is a statistically significant similarity between one or more sequences of the captured cell-free methylated DNA and a cell-free methylated DNA sequence from an individual with cancer; wherein in at least one of the capturing, comparing or identifying steps, the cell-free methylated DNA of the subject is restricted to a subpopulation according to a fragment length metric.

[0005] In one embodiment, the present disclosure provides a method for determining whether a subject has a disease or is at risk of having a disease. The method includes: sequencing a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from a subject to generate at least one profile selected from the group consisting of (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile; and processing the at least one profile to determine whether the subject has the disease or is at risk of having a disease with at least 80% sensitivity or at least about 90% specificity, wherein the cell-free nucleic acid sample comprises less than 30 nanograms (ng) / milliliter (ml) of a plurality of nucleic acid molecules.

[0006] In some embodiments, the cell-free nucleic acid sample comprises less than 10 ng / ml of said plurality of nucleic acid molecules. In some embodiments, the cell-free nucleic acid sample comprises less than 5 ng / ml of said plurality of nucleic acid molecules. In some embodiments, the cell-free nucleic acid sample comprises less than 1 ng / ml of said plurality of nucleic acid molecules. In some embodiments, the subjecting to (a) produces at least two profiles selected from the group consisting of (i), (ii) and (iii). In some embodiments, the at least two profiles comprise said methylation profile and said fragment length profile.

[0007] In some embodiments, the at least two profiles comprise the mutation profile and the fragment length profile. In some embodiments, the at least two profiles comprise the methylation profile and the mutation profile. In some embodiments, the subjecting step of (a) produces the methylation profile, the mutation profile and the fragment length profile.

[0008] In another aspect, the disclosure provides a method of processing a cell-free nucleic acid sample of a subject to determine whether the subject has or is at risk of having a disease, the method comprising obtaining the cell-free nucleic acid sample comprising a plurality of nucleic acid molecules, subjecting the plurality of nucleic acid molecules or derivatives thereof to sequencing to generate a plurality of sequencing reads, computationally processing the plurality of sequencing reads to identify (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile for the plurality of nucleic acid molecules, and using at least the methylation profile, the mutation profile, and the fragment length profile to determine whether the subject has or is at risk of having the disease.

[0009] In some embodiments, the disease comprises cancer. In some embodiments, the cancer is adrenal cancer, anal cancer, bile duct cancer, bladder cancer, bone cancer, brain / CNS tumors, breast cancer, Castleman's disease, cervical cancer, colon / rectal cancer, endometrial cancer, esophageal cancer, Ewing's family of tumors, eye cancer, gallbladder cancer, gastrointestinal carcinoid tumors, gastrointestinal stromal tumors (gist), gestational trophoblastic disease, Hodgkin's disease, Kaposi's sarcoma, kidney cancer, laryngeal and hypopharyngeal cancer, leukemia (acute lymphocytic, acute myeloid, chronic lymphocytic, chronic myeloid, chronic myelomonocytic), liver cancer, lung cancer (non-small cell, small cell, lung carcinoid tumors), lymphoma, lymphoma of the skin, malignant mesothelioma, The cancer is selected from the group consisting of cancers selected from the group consisting of multiple myeloma, myelodysplastic syndromes, nasal and paranasal sinus cancer, nasopharyngeal cancer, neuroblastoma, non-Hodgkin's lymphoma, oral and oropharyngeal cancer, osteosarcoma, ovarian cancer, penile cancer, pituitary cancer, prostate cancer, retinoblastoma, rhabdomyosarcoma, salivary gland cancer, sarcoma-adult soft tissue cancer, skin cancer (basal and squamous cell, melanoma, Merkel cell), small intestine cancer, gastric cancer, testicular cancer, thymic cancer, thyroid cancer, uterine sarcoma, vaginal cancer, vulvar cancer, Waldenstrom's macroglobulinemia, Wilms' tumor, squamous cell carcinoma, and head and neck squamous cell carcinoma. In some embodiments, the cancer is squamous cell carcinoma. In some embodiments, the cancer is head and neck squamous cell carcinoma.

[0010] In some embodiments, the plurality of cell-free nucleic acid molecules comprises circulating tumor nucleic acid molecules. In some embodiments, the circulating tumor nucleic acid comprises circulating tumor DNA. In some embodiments, the circulating tumor nucleic acid comprises circulating tumor RNA. In some embodiments, the methylation profile comprises a plurality of differentially methylated regions (DMRs). In some embodiments, the plurality of DMRs is from ctDNA. In some embodiments, the plurality of DMRs from peripheral blood leukocytes are removed from the methylation profile. In some embodiments, the plurality of DMRs comprises at least about 56 genomic regions having a hypomethylation level compared to corresponding genomic regions from normal healthy subjects. In some embodiments, the plurality of DMRs comprises at least about 941 genomic regions having a hypermethylation level compared to corresponding genomic regions from normal healthy subjects. In some embodiments, the DMR comprises a size of at least about 300 bp. In some embodiments, the DMR comprises a size of at least about 100 bp to at least about 200 bp. In some embodiments, the DMR comprises a size of at least about 100 bp to at least about 150 bp. In some embodiments, the DMR comprises at least 8 CpG genomic islands. In some embodiments, a normal healthy subject comprises the same set of risk factors as said subject.

[0011] In some embodiments, the mutation profile includes missense, nonsense, deletion, insertion, duplication, inversion, frameshift, or repeat expansion mutants. In some embodiments, any mutations present in a genomic DNA sample obtained from a plurality of peripheral blood leukocytes, the plurality of peripheral blood leukocytes being obtained from the subject, are removed from the mutation profile. In some embodiments, any mutations derived from clonal hematopoiesis are removed from the mutation profile. In some embodiments, the mutation profile does not include mutations in genes DNMT3A, TET2, or ASXL1. In some embodiments, the mutation profile does not include canonical cancer driver genes. In some embodiments, the mutation profile includes non-canonical cancer driver genes, the non-canonical genes being GRIN3A or MYC.

[0012] In some embodiments, the fragment length profile comprises selecting cell-free nucleic acid molecules based on a fragment length range of at least about 80 bp to 170 bp. In some embodiments, the fragment length profile comprises selecting cell-free nucleic acid molecules based on a fragment length range of at least about 100 bp to 150 bp. In some embodiments, circulating tumor nucleic acid molecules are enriched.

[0013] In some embodiments, the method further comprises mixing the cell-free nucleic acid sample with filler DNA molecules to produce a DNA mixture. In some embodiments, the filler DNA molecules comprise a length of about 50 bp to 800 bp. In some embodiments, the filler DNA molecules comprise a length of about 100 bp to 600 bp. In some embodiments, the filler DNA molecules comprise at least about 5% methylated filler DNA molecules. In some embodiments, the filler DNA molecules comprise at least about 20% methylated filler DNA. In some embodiments, the filler DNA molecules comprise at least about 30% methylated filler DNA. In some embodiments, the filler DNA molecules comprise at least about 50% methylated filler DNA.

[0014] In some embodiments, the method further comprises incubating the DNA mixture with a binder configured to bind methylated nucleotides to generate an enriched sample. In some embodiments, the binder comprises a protein comprising a methyl-CpG binding domain. In some embodiments, the protein is an MBD2 protein. In some embodiments, the binder comprises an antibody. In some embodiments, the antibody is a 5-MeC antibody. In some embodiments, the antibody is a 5-hydroxymethylcytosine antibody. In some embodiments, the sequencing does not comprise bisulfite sequencing. In some embodiments, the cell-free nucleic acid sample comprises a blood sample. In some embodiments, the blood sample comprises a plasma sample. In some embodiments, the method further comprises detecting the origin of the cancer tissue.

[0015] In some embodiments, the method further comprises generating a report comprising a prognosis of survival of the subject. In some embodiments, the method further comprises administering a treatment to the subject. In some embodiments, following treatment of the disease, the method further comprises providing a second report indicating whether the treatment is effective.

[0016] In another aspect, the disclosure provides a method for determining whether a subject has or is at risk of having a condition, the method comprising: assaying cell-free nucleic acid molecules obtained from at least a portion of a sample from the subject; detecting a methylation level of at least a portion of said cell-free nucleic acid molecule that falls within a differentially methylated region (DMR) listed in Table 5; and comparing, using at least one computer processor, the methylation level detected in (b) with the methylation level of a corresponding portion(s) of said cell-free nucleic acid molecule that falls within a DMR listed in Table 5.

[0017] In some embodiments, the cell-free nucleic acid molecule comprises ctDNA. In some embodiments, the method comprises performing a sequencing analysis, wherein the sequencing analysis comprises cell-free Methylated DNA ImmunoPrecipitation (cfMeDIP) sequencing. In some embodiments, the detecting step comprises measuring the methylation level of at least a portion of the nucleic acid molecule contained in 6 or more, 10 or more, 15 or more, 20 or more, 30 or more, 40 or more, 50 or more, 60 or more, 70 or more, 80 or more, 90 or more, or 100 or more DMRs listed in Table 5.

[0018] In another aspect, the disclosure provides a method of determining whether a subject has a higher survival rate after undergoing treatment for a disease, comprising: assaying cell-free nucleic acid molecules from at least a portion of a sample from the subject; detecting a methylation level of at least a portion of the cell-free nucleic acid molecule that falls within a differentially methylated region (DMR) listed in Table 6; and processing, using at least one computer processor, the methylation level detected in (b) into a methylation level of a corresponding portion(s) of the cell-free nucleic acid molecule that falls within the DMR listed in Table 6.

[0019] In some embodiments, the cell-free nucleic acid molecule comprises ctDNA. In some embodiments, the detecting step comprises providing a complex methylation score (CMS). In some embodiments, the CMS comprises the sum of the beta values ​​of the DMRs listed in Table 6. In some embodiments, a higher CMS indicates a lower survival rate of the subject. In some embodiments, the CMS is independent of the amount of ctDNA present. In some embodiments, the disease is squamous cell carcinoma. In some embodiments, the cancer is head and neck squamous cell carcinoma.

[0020] In another aspect, the disclosure provides a system for determining whether a subject has or is at risk of having a disease, the system comprising one or more computer processors individually or collectively programmed to perform a process comprising subjecting a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from the subject to sequencing to generate at least one of: (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile, and processing the at least one profile to determine whether the subject has or is at risk of having the disease with at least 80% sensitivity or at least about 90% specificity, wherein the cell-free nucleic acid sample comprises less than 30 ng / ml of the plurality of nucleic acid molecules.

[0021] In another aspect, the disclosure provides a system for processing a cell-free nucleic acid sample of a subject to determine whether the subject has or is at risk of having a disease, the system comprising one or more computer processors individually or collectively programmed to perform a process comprising obtaining the cell-free nucleic acid sample comprising a plurality of nucleic acid molecules; subjecting the plurality of nucleic acid molecules or derivatives thereof to sequencing to generate a plurality of sequencing reads; computer processing the plurality of sequencing reads to identify (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile for the plurality of nucleic acid molecules; and using at least the methylation profile, the mutation profile, and the fragment length profile to determine whether the subject has or is at risk of having a disease.

[0022] These and other features of preferred embodiments of the present invention will become more apparent in the following detailed description, which refers to the accompanying drawings. [Brief description of the drawings]

[0023] [Figure 1] Utilization of PBL filtering for detection of ctDNA by CAPP-Seq. A) Mutant allele fraction of candidate SNVs identified in matched patient plasma and / or PBL. Pearson correlation was performed for SNVs strictly found in both matched patient plasma and PBL. Candidate SNVs found only in patient plasma are shown inside dashed red boxes. B) Occurrence of candidate SNVs identified in both matched patient plasma and PBL. The top histogram shows the number of SNVs per patient, and the right histogram shows the number of patients with a particular gene mutated. C) Average MAF of candidate SNVs across cfDNA (red circles) and PBL (blue circles) of HNSCC patients before and after removal of PBL-associated SNVs. Patients with SNVs not present after PBL filtering show false positive detection of ctDNA. E) Occurrence of selected PBL-filtered SNVs identified in 20 / 32 HNSCC patients. The top and right histograms show histograms as previously described in (B). F) Average mutant allele percentage of PBL-filtered SNVs across all HNSCC patients. For each SNV per patient, the percentage of mutant alleles was calculated by the fraction of reads containing the SNV of interest compared to reads containing native sequences overlapping the SNV base pair position. [Diagram 2]Utilization of PBL filtering for detection of ctDNA by CAPP-Seq. B) Mutant allele fraction of candidate SNVs identified in matched patient plasma and / or PBL. Pearson correlation was performed for SNVs strictly found in both matched patient plasma and PBL. Candidate SNVs found only in patient plasma are shown inside dashed red boxes. C) Occurrence of candidate SNVs identified in both matched patient plasma and PBL. The top histogram shows the number of SNVs per patient, and the right histogram shows the number of patients with a particular gene mutated. D) Average MAF of candidate SNVs across cfDNA (red circles) and PBL (blue circles) of HNSCC patients before and after removal of PBL-associated SNVs. Patients with SNVs not present after PBL filtering show false positive detection of ctDNA. E) Occurrence of selected PBL-filtered SNVs identified in 20 / 32 HNSCC patients. The top and right histograms show histograms as previously described in (B). F) Average mutant allele percentage of PBL-filtered SNVs across all HNSCC patients. For each SNV per patient, the percentage of mutant alleles was calculated by the fraction of reads containing the SNV of interest compared to reads containing native sequences overlapping the SNV base pair position. [Diagram 3]Identification of informative regions for detection of ctDNA by cfMeDIP-seq. B) Pearson correlation of 300 bp non-overlapping windows with 8 or more CpGs from patients and healthy donors against FaDu genomic DNA (gDNA) [1 × 1 × 52 comparisons], mismatched PBL gDNA [1 × 51 × 52 comparisons], and matched PBL gDNA [1 × 1 × 52 comparisons] MeDIP-seq profiles (n=52). C) Performance of in silico PBL depletion on healthy donor (right) and HNSCC (left) PBL MeDIP-seq profiles. Absolute methylation scores were calculated by MeDEStrand (Methods) from MeDIP-seq counts. The 300 bp non-overlapping windows before PBL depletion (blue) correspond to all windows from chromosome 1 to chromosome 22 with 8 or more CpGs (n=702,488). The 300 bp non-overlapping window after PBL depletion (red) contains an additional filter where the median absolute methylation across PBLs of healthy donors is <0.1 (n=99,997). D) Workflow of ctDNA detection by differential methylation analysis of cfMeDIP-seq profiles of HNSCC and healthy donors. To identify HNSCC-associated cfDNA methylation, cfMeDIP-seq profiles from HNSCC patients with detectable SNVs by CAPP-Seq (i.e., CAPP-Seq positive, n=20) were compared to healthy donors (n=20) within the PBL depletion window. Hyper- and hypomethylated regions are indicated as regions with higher or lower methylation in the HNSCC cohort compared to healthy donors at FDR<10%. E) Permutation analysis of hypermethylated regions annotated by CpG sites (n=10,000 total permutations). Significant enrichment / depletion is shown as z-score observed with p-value less than 0.05. F, Permutation analysis of hypermethylated regions within tumor-specific methylated cytosines from TCGA (n=1000 permutations total). Significant enrichment / depletion is shown as z-score observed with p-value less than 0.05. [Figure 4]Concordance of ctDNA detection and abundance between CAPP-Seq and cfMeDIP-seq profiles. A) Median fragment lengths of detected SNVs across HNSCC patients by CAPP-seq. For each patient, the median fragment length of each SNV and the matched reference allele was measured. The distribution of median fragment lengths for each mutation or matched reference allele is shown for each patient. The ends of the box and centerline define the upper and lower quartiles and median, respectively. In cases with a single SNV, the colored line indicates the median fragment length containing the SNV or matched reference allele, respectively. B) Distribution of fragment lengths within HNSCC hypermethylated regions by cfMeDIP-seq. Fragment lengths from healthy donors were pooled prior to analysis, after which each box represents an individual HNSCC cfMeDIP-seq profile. The ends of the box and centerline define the upper and lower quartiles and median, respectively. Individual HNSCC samples are ordered based on the average increase in methylation (RPKM) within the hypermethylated regions. The dashed blue line defines the median fragment length across all healthy donors. C) Ratio of enrichment of excess DMR regions by 100-150 bp fragments compared to enrichment of excess DMR regions by 100-220 bp fragments. Ratios were converted to percent increase / decrease for ease of interpretation. D) Ratio of enrichment of excess DMR regions by 100-150 bp fragments compared to enrichment of excess DMR regions by 100-220 bp fragments. + signs indicate HNSCC patients with detectable ctDNA by CAPP-Seq (CAPP-Seq positive). E) Supervised hierarchical classification of cfMeDIP-seq profiles restricted to 100bp-150bp by log-transformed RPKM values ​​across hypermethylated regions in HNSCC. RPKM values ​​for each cfMeDIP-seq profile were log2-transformed before Euclidean transformation and clustered using Ward's method. Methylation clusters were defined with a threshold of k=4. F) Relationship between mean variant allele frequency and mean RPKM from identified SNVs and hypermethylated regions by CAPP-seq and cfMeDIP-seq (restricted to 100–150 bp), respectively.Dots indicate individual samples from HNSCC or healthy donor plasma. The solid red line and shaded grey area indicate the fitted linear regression model and associated 95% confidence intervals, respectively. G) AUROC analysis based on methylation values ​​(restricted to 100-150 bp) within hypermethylated regions of HNSCC comparing HNSCC to cfMeDIP-seq profiles of healthy donors. Detection of ctDNA was defined when the mean of methylation exceeded the maximum value across healthy donors. H) Kaplan-Meier curve analysis of overall survival of patients within methylation cluster 1+2+3 compared to methylation cluster 4. I+J) Comparison of median fragment lengths from CAPP-Seq and cfMeDIP-seq profiles (I) and median fragment lengths from CAPP-Seq to the ratio of 100-150:151-220 bp from cfMeDIP-seq profiles (J). Dots defined individual HNSCC samples within methylation clusters 1 and 2. The solid red line and shaded grey area indicate the fitted linear regression model and 95% confidence intervals, respectively. [Diagram 5]Prognostic utility of specific methylated regions within ctDNA detected by cfMeDIP-seq. A) Relationship between mean mutant allele fraction and mean RPKM from identified mutated and hypermethylated regions by CAPP-seq and cfMeDIP-seq (restricted to 100-150bp), respectively. Dots indicate individual samples from HNSCC or healthy control plasma. Solid red line: fitted linear regression model. Grey boundaries: 95% confidence interval. B) Kaplan-Meier analysis showing overall survival of patients with detectable ctDNA by both CAPP-Seq and cfMeDIP-seq (mean methylation above healthy controls within hyper-DMRs). C) Identification of prognostic regions based on disease-specific survival by multivariate Cox proportional hazards regression analysis across HNSCC primary tumors presented by TCGA (n=520). Regions were defined as 300bp windows as previously described. 450K data of human methylation were obtained from TCGA and β values ​​from probe IDs overlapping with each region were averaged. Candidate regions for prognostic analysis were selected based on elevated methylation across primary tumors (n=520) compared to solid adjacent normal tissues (n=50) (Wilcoxon test, adjusted p-value <0.05, log2FC>1). G-H) Spearman correlation from methylation of specific 300 bp regions (boxes) to RNA expression of specific transcripts. Regions with absolute R-values ​​≥0.3 (indicated with grey dashed lines) were labeled as significantly associated. Methylated regions (n=520) that were prognostic for disease-specific survival of HNSCC patients presented by TCGA are outlined in red. Prognostic regions further associated with RNA expression are shown in solid red. Examples of prognostic methylated regions associated with RNA expression; (G) OSR1, (H) LINC01391. E) Kaplan-Meier curves of overall survival for HNSCC-TCGA patients based on total methylation across five regions affecting expression of ZNF323 / ZSCAN1, LINC01391, GATA-AS1, OSR1, and STK3 / MST2, respectively. Patients were stratified based on being either below (Blw med. blue) or above (Abv med. red) the median total methylation of the five regions previously identified in (D) across all primary tumors.F) Kaplan-Meier curves of overall survival as described in (E) for the HNSCC plasma cohort with detectable ctDNA by CAPP-Seq. To calculate total methylation across the five genes with prognostic association, RPKM values ​​were appropriately scaled across all previously identified over-DMR regions prior to survival analysis. [Figure 6] Clinical utility of ctDNA detection by cfMeDIP-seq for longitudinal monitoring. A) ctDNA dynamics are typically observed across patients throughout treatment. Complete clearance was defined as a change from detected ctDNA at diagnosis to a decrease in ctDNA abundance below the detection threshold (i.e., 0.2%) at the first available on-treatment / post-treatment time point. Partial clearance was defined as a change from detected ctDNA at diagnosis to a decrease in ctDNA abundance above the detection threshold (>90%) at the first available on-treatment / post-treatment time point. No clearance was defined as an increase in ctDNA abundance in on-treatment / post-treatment samples compared to at diagnosis. lastFU=sample collection at last follow-up, RT=radiotherapy. B) Change in ctDNA abundance at diagnosis to the first available on-treatment / post-treatment time point across HNSCC patients (n=30). Red lines display patients with no clearance dynamics, grey lines display patients with clearance / partial clearance dynamics. C, Kaplan-Meier curves of recurrence-free survival. Patients were stratified based on clearance kinetics (i.e., no clearance vs. clearance / partial clearance). [Figure 7]Comparison of cfMeDIP-seq analyses performed on all or ctDNA-enriched fragments. ctDNA-enriched fragments are defined as fragments ranging from 100 to 150 bp in length. A) Mutant allele frequency of mutations identified by CAPP-Seq vs mean RPKM of previously identified excess DMRs in HNSCC in cfMeDIP-seq profiles including all fragments (left) or ctDNA-enriched fragments (right). B) Area under the curve analysis (AUROC) for ctDNA detection in cfMeDIP-seq profiles of HNSCC versus healthy donors (CAPP-Seq positive only: red, CAPP-Seq positive and negative: blue). Results of cross-validation analysis using CAPP-Seq positive patients are also shown (replicates = 50). Analyses are shown for cfMeDIP-seq profiles using all fragments (left) or ctDNA-enriched fragments (right). C) Kaplan-Meier analysis of recurrence-free survival based on longitudinal cfMeDIP-seq profiling using all fragments (left) or ctDNA-enriched fragments. Patients were classified as positive for post-treatment ctDNA if they demonstrated methylation abundance within previously identified excess DMRs greater than 0.2%. [Figure 8] 1 illustrates a computer system that can be programmed or otherwise configured to perform the methods provided herein. [Figure 9] Characterization of cell-free DNA samples isolated from HNSCC and healthy donors. A) Schematic defining blood isolation time points. B) cfDNA yield (normalized to per mL of plasma) across time points for HNSCC patients as well as healthy donors (i.e., "normals"). [Figure 10] Analysis of the number of SNVs per HNSCC patient covered by the CAPP-Seq selector, assessed among all 364 patients in the TCGA cohort of HNSC (blue diamonds) or using leave-one-out cross-validation (LOOCV; red squares). [Figure 11] Incidence of all PBL-filtered SNVs identified in 20 / 32 HNSCC patients (related to Figure 2E ). [Figure 12] Relevant figures for identifying informative regions (related to Figure 3B and Figure 3C). A) Median RPKM of genome-wide (chromosomes 1-22) 300 bp non-overlapping bins based on n or more CpGs. B) Differential methylation analysis between HNSCC and PBLs from healthy donors within the PBL depletion window as described in Figure 2B and Methods. Hypomethylated regions (i.e., regions with elevated methylation in PBLs from healthy donors) are shown in blue. [Figure 13] Figure 2D: Related results of differential methylation analysis between HNSCC and healthy donor cfDNA samples within the PBL depletion window. A) DMRs were defined based on the original 300bp non-overlapping windows used for the initial analysis. DMRs directly adjacent to each other were binned to their respective widths (i.e., two 300bp windows were defined, each independently, to have a length of 600bp). B) Permutation analysis of the CpG features defined in Figure 2E based on hypomethylated regions. [Figure 14] Supervised hierarchical clustering of TCGA primary tumors based on the identification of cancer-specific differentially methylated cytosines. Cancer_type (columns) refers to the classification of each primary tumor or PBL sample, and cancer_DMCs (rows) refers to the cancer-specific differentially methylated cytosines identified for each cancer type (PBLs are excluded). [Figure 15] Figure associated with Figure 4. A) Median fragment length of CAPP-Seq identified SNVS per patient compared to mean mutant allele fraction. B) Median fragment length within excess cfMeDIP-seq DMRS per patient compared to mean RPKM of excess DMRs. [Figure 16]Related figures of CAPP-Seq and cfMeDIP-seq concordance analysis (Figure 4E). A) Area under the curve values ​​obtained from cross-validation analysis (n=50) of differentially methylated region calls between CAPP-Seq positive HNSCC cfDNA samples and healthy donors. B) Kaplan-mir analysis of overall survival of HNSCC patients based on detection of ctDNA by CAPP-Seq. C) and D) Mean RPKM and mean mutant allele fraction of HNSCC patient samples stratified based on methylation clusters (Figure 4D). [Figure 17] Identification of regions of potential clinical utility (related to Figure 6). A) Genomic tracks of genes currently used in commercial liquid biopsy tests that overlap with plasma-derived excess DMRs from HNSCC primary tumors in TCGA as well as our HNSCC cohort. The dark blue bar below with arrows indicates the direction of transcription of the specific gene. The red bar indicates the position of a 300 bp window that overlaps with plasma-derived excess DMRs of our HNSCC cohort as well as primary tumors from TCGA. B-D) Spearman correlation from methylation of specific 300 bp regions (boxes) to RNA expression of specific transcripts. Regions with absolute R values ​​≥ 0.3 (indicated by grey dashed lines) were labeled as significant associations. Methylated regions (n=520) that were prognostic for disease-specific survival of HNSCC patients presented by TCGA are shown with a red outline. Prognostic regions that were further associated with RNA expression are shown in solid red. Figures were generated for all five genes with prognostic methylated regions associated with RNA expression. (B) GATA2-AS1, (C) ZNF323, (D) STK3. [Figure 18] Zoom in on Figure 6A showing changes in ctDNA abundance by cfMeDIP-seq across treatment for all HNSCC patients (n=32). DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0024] In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be understood that the present invention may be practiced without these specific details.

[0025] The present disclosure provides methods, systems, and kits for multimodal analysis of ctDNA in determining the likelihood that a subject has cancer with high sensitivity and / or high specificity. Additionally, the present disclosure provides methods, systems, and kits for detecting minimal residual disease (MRD) after treatment of cancer and assessing whether treatment of such cancer is therapeutically effective.

[0026] Identifying specific molecular features from ctDNA before treatment can inform prognosis and / or predict response to treatment, whereas detection of ctDNA after treatment can facilitate identification of MRD and aid in identifying patients at high risk of recurrence and / or death. To achieve robust sensitivity, most clinical studies utilize ctDNA detection methods that interrogate a small number of regions, matched tumor profiling, and / or cases of high ctDNA abundance. However, for cancers with low levels of ctDNA or lacking common / known abnormalities across patients, additional strategies may be utilized to achieve a similar degree of sensitivity. Genome-wide profiling techniques can help improve sensitivity by covering a very large number of regions. However, the amount of cell-free DNA and sequencing depth required to achieve detection of less than 1% of the fragments have been cost-prohibitive.

[0027] Two personalized genome-wide profiling techniques capable of sensitive ctDNA detection have been described. The first, CAncer Personalized Profiling by deep Sequencing (CAPP-Seq), utilizes broad hybrid capture probes targeting over 100 genes to identify low-frequency allelic variants. The second, cell-free Methylated DNA ImmunoPrecipitation sequencing (cfMeDIP-seq), enriches for methylated cfDNA fragments by using anti-5-methylcytosine (anti-5mC) antibodies. Identification of mutations or hypermethylation events by each of these methods has its own advantages. Mutations can distinguish ctDNA from healthy sources of cell-free DNA due to their irreversible nature, provided that appropriate error suppression tools are used and any contribution of mutations from clonal hematopoiesis is considered. DNA hypermethylation events potentially affect a higher number of recurrent genomic regions in cancer, contributing to their ability to inform tumor origin by cell-free DNA analysis. Moreover, hypermethylation events in the vicinity of cancer driver genes may affect their expression, thereby potentially reflecting cancer behavior and providing prognostic value.To date, no studies have utilized the combination of both mutation-based and methylation-based methods to improve tumor-naive detection and characterization of ctDNA in localized cancers.

[0028] Utilizing fluid-based biomarkers for prognostication, risk stratification, and disease monitoring may improve patient outcomes by guiding treatment decisions without the need for invasive tumor sampling. Circulating tumor (ct)DNA has shown promise, especially as a liquid biopsy tool, but paired tumor profiling is frequently required in patients with low disease burden, such as those with localized nonmetastatic cancer. We hypothesized that multimodal analysis of genetic and epigenetic features from plasma cell-free DNA may enable widespread application of tumor-naive ctDNA profiling. Mutation- and methylation-based profiling identified ctDNA in 65% of patients with localized head and neck cancer. Results from both approaches were quantitative and strongly correlated, and their combined analysis revealed common features of tumor-derived DNA fragments. Furthermore, the ctDNA methylome revealed dynamic patterns of tumor histology, putative prognostic biomarkers, and treatment response. These findings will aid future non-invasive biomarker discovery efforts and inform the clinical implementation of ctDNA for localized cancers.

[0029] Certain methods for capturing cell-free methylated DNA are described in Applicant's WO 2017 / 190215 and WO 2019 / 010564, both of which are incorporated by reference.

[0030] Specifically, we utilize both CAPP-Seq and cfMeDIP-seq to perform tumor-naive ctDNA detection within a cohort of patients with localized head and neck squamous cell carcinoma (HNSCC). HNSCC is a clinically heterogeneous disease that frequently recurs after curative treatment and can greatly benefit from ctDNA detection to better inform treatment decisions and disease management. 33We demonstrate that utilizing both methods in parallel, as well as matched PBL profiling, can achieve reliable tumor-naive ctDNA detection. Furthermore, we show that the combined analysis reveals common molecular signatures of tumor-derived DNA fragments. Finally, we show that the ctDNA methylome reveals dynamic patterns of tumor histology, putative prognostic biomarkers, and treatment response, imparting a blueprint for future biomarker research in other disease contexts.

[0031] In an embodiment, there is provided a method for detecting the presence of ctDNA from a cancer cell of a subject, comprising: (a) obtaining a sample of cell-free DNA from a subject; (b) subjecting the sample to library preparation to allow subsequent sequencing of cell-free methylated DNA; (c) a first amount of filler DNA may be added to the sample, at least a portion of the filler DNA may be methylated, and then the sample may be further denatured; (d) capturing the cell-free methylated DNA with a binding agent selective for methylated polynucleotides; (e) sequencing the captured cell-free methylated DNA; (f) comparing the sequence of the captured cell-free methylated DNA with control cell-free methylated DNA sequences from healthy individuals and individuals with cancer; (g) identifying the presence of DNA from cancer cells if there is a statistically significant similarity between one or more sequences of the captured cell-free methylated DNA and a cell-free methylated DNA sequence from an individual with cancer; wherein in at least one of the capturing, comparing or identifying steps, the cell-free methylated DNA of the subject is restricted to a subpopulation according to a fragment length metric.

[0032] Various sequencing techniques are known to those skilled in the art, such as polymerase chain reaction (PCR) followed by Sanger sequencing. Next generation sequencing (NGS) techniques, also known as high-throughput sequencing, are also available, including various sequencing techniques, including Illumina (Solexa) sequencing, Roche 454 sequencing, Ion torrent:Proton / PGM sequencing, SOLiD sequencing, and long-read sequencing (Oxford Nanopore and Pactbio). NGS allows for DNA and RNA sequencing much faster and cheaper than previously used Sanger sequencing. In some embodiments, the sequencing is optimized for short-read sequencing.

[0033] As used herein, the term "subject" refers to any member of the animal kingdom.Thus, the methods and described herein are applicable to both human and veterinary diseases and animal models.Preferred subject is a "patient", i.e., a living human being who is being investigated to determine whether a disease or condition requires treatment or medical care, or who is receiving medical care for a disease or condition (e.g., cancer).

[0034] As used herein, the term "genome" generally refers to genomic information from a subject, which may be, for example, at least a portion or the entirety of the subject's genetic information. A genome may be encoded in either DNA or RNA. A genome may include coding regions (e.g., protein-encoding) as well as non-coding regions. A genome may include the sequences of all chromosomes together in an organism. For example, the human genome usually has a total of 46 chromosomes. All these sequences together may constitute the human genome.

[0035] As used herein, the term "nucleic acid" refers to a polynucleotide comprising two or more nucleotides, either deoxyribonucleotides (dNTPs) or ribonucleotides (rNTPs), or analogs thereof, in a polymeric form of nucleotides of any length. Non-limiting examples of nucleic acids include deoxyribonucleic acid (DNA), ribonucleic acid (RNA), coding or non-coding regions of a gene or gene fragment, loci defined from linkage analysis, exons, introns, messenger RNA (mRNA), transfer RNA, ribosomal RNA, short interfering RNA (siRNA), short hairpin RNA (shRNA), microRNA (miRNA), ribozymes, cDNA, recombinant nucleic acids, branched nucleic acids, plasmids, vectors, isolated DNA of any sequence, isolated RNA of any sequence, nucleic acid probes and primers. Nucleic acids may contain one or more modified nucleotides, such as methylated nucleotides and nucleotide analogs. If present, modifications to the nucleotide structure may be made before or after construction of the nucleic acid. The sequence of nucleotides of a nucleic acid may be interrupted by non-nucleotide components. The nucleic acid may be further modified after polymerization, for example by conjugation or binding with a reporter factor. A "mutant" nucleic acid is a polynucleotide that has the same nucleotide sequence as its original nucleic acid, except that at least one nucleotide is modified, for example deleted, inserted or substituted, respectively. A mutant may have a nucleotide sequence that is at least about 80%, 90%, 95%, or 99% identical to the nucleotide sequence of the original nucleic acid.

[0036] Cell-free methylated DNA is DNA that circulates freely in the bloodstream and is methylated at various regions of the DNA. A sample, for example a plasma sample, may be taken and cell-free methylated DNA analyzed. Studies have revealed that much of the circulating nucleic acid in blood originates from necrotic or apoptotic cells, and a significant increase in the level of nucleic acid from apoptosis is observed in diseases such as cancer. In particular, in the case of cancers where circulating DNA bears hallmarks of disease including mutations in cancer genes, microsatellite alterations, and in the case of certain cancers, viral genome sequences, DNA or RNA in plasma are increasingly being investigated as potential biomarkers of disease. For example, quantitative assays of low levels of circulating tumor DNA in total circulating DNA may serve as a better marker for detecting colorectal cancer recurrence compared to carcinoembryonic antigen, the standard biomarker used clinically. Circulating cfDNA may include circulating tumor DNA (ctDNA).

[0037] As used herein, "library preparation" includes end-repair, A-tailing, adapter ligation, or any other preparation performed on cell-free DNA to allow for subsequent sequencing of the DNA.

[0038] As used herein, "filler DNA" may be non-coding DNA or may consist of amplicons.

[0039] In some embodiments, the fragment length metric is the fragment length. In some preferred embodiments, the subject's cell-free methylated DNA is limited to fragments having a length of <170bp, <165bp, <160bp, <155bp, <150bp, <145bp, <140bp, <135bp, <130bp, <125bp, <120bp, <115bp, <110bp, <105bp, or <100bp. In other preferred embodiments, the subject's cell-free methylated DNA is limited to fragments having a length of about 100 to about 150bp, 110 to 140bp, or 120 to 130bp.

[0040] In some embodiments, the fragment length metric is a fragment length distribution of the subject's cell-free methylated DNA. In some preferred embodiments, the subject's cell-free methylated DNA is restricted to fragments within the bottom 50, 45, 40, 35, 30, 25, 20, 15, or 10 percentile based on length.

[0041] In some embodiments, the subject's cell-free methylated DNA is further restricted to fragments of differentially methylated regions (DMRs).

[0042] In some embodiments, the restriction of cell-free methylated DNA of interest is during the capture step.

[0043] In some embodiments, the limitation of the cell-free methylated DNA of the subject is during the comparison step.

[0044] In some embodiments, restriction of the subject's cell-free methylated DNA is during the identification step.

[0045] In some embodiments, the comparison step is based on fitting using a statistical classifier. A statistical classifier using DNA methylation data can be used to assign samples to a particular disease state, such as a cancer type or subtype. For the purpose of cancer type or subtype classification, the classifier consists of one or more DNA methylation variables (i.e., features) in a statistical model, the output of which has one or more thresholds for distinguishing between different disease states. The particular feature(s) and threshold(s) used in the statistical classifier can be derived from prior knowledge of the cancer type or subtype, from prior knowledge of the features most likely to be informative, from machine learning, or from a combination of two or more of these approaches.

[0046] In some embodiments, the classifier is derived by machine learning. Preferably, the classifier is an elastic net classifier, a lasso, a support vector machine, a random forest, or a neural network.

[0047] The genomic space analyzed can be genome-wide or preferably limited to regulatory regions (i.e., FANTOM5 enhancers, CpG islands, CpG shores, and CpG shelves).

[0048] Preferably, the percentage of spike-in methylated DNA recovered is included as a covariate to control for variation in pull-down efficiency.

[0049] For classifiers that can distinguish multiple cancer types (or subtypes) from one another, the classifier preferably consists of differentially methylated regions derived from pairwise comparisons of each type (or subtype) of interest.

[0050] In some embodiments, control cell-free methylated DNA sequences from healthy individuals and individuals with cancer are included in a database of differentially methylated regions (DMRs) between healthy and cancer individuals.

[0051] In some embodiments, the control cell-free methylated DNA sequences from healthy individuals and cancer individuals are limited to control cell-free methylated DNA sequences that are differentially methylated between healthy individuals and cancer individuals in DNA derived from cell-free DNA obtained from bodily fluids, such as serum, cerebrospinal fluid, urine stool, sputum, pleural fluid, peritoneal fluid, tears, sweat, Pap smear fluid, endoscopic brushing fluid, and the like, preferably from plasma.

[0052] sample The sample may be any biological sample isolated from a subject. For example, the sample may include, but is not limited to, bodily fluids, whole blood, platelets, serum, plasma, stool, red blood cells, white blood cells, endothelial cells, tissue biopsies, synovial fluid, lymphatic fluid, ascites fluid, interstitial or extracellular fluid, fluid in the space between cells, such as gingival fluid, bone marrow, cerebrospinal fluid, saliva, mucus, sputum, semen, sweat, urine, fluid from nose brushing, fluid from pep smear, or any other bodily fluid. The bodily fluid may include saliva, blood, or serum. The sample may also be a tumor sample, which may be obtained from a subject by various approaches, including, but not limited to, venipuncture, excretion, ejaculation, massage, biopsy, needle aspiration, lavage, scraping, surgical incision, or intervention or other approaches. The sample may be a cell-free sample (e.g., substantially free of cells). The DNA sample may be denatured, for example, using sufficient heat.

[0053] In some embodiments, the present disclosure provides a system, method, or kit that includes or uses one or more biological samples.As used herein, one or more samples can include any material that contains or is suspected to contain nucleic acid.A sample can include a biological sample obtained from a subject.In some embodiments, the biological sample is a liquid sample.

[0054] In some embodiments, the sample comprises less than about 100 ng, 90 ng, 80 ng, 75 ng, 70 ng, 60 ng, 50 ng, 40 ng, 30 ng, 20 ng, 10 ng, 5 ng, 1 ng, or any amount in between the number of cell-free nucleic acid molecules. Further, in some embodiments, the sample comprises less than about 1 pg, less than about 5 pg, less than about 10 pg, less than about 20 pg, less than about 30 pg, less than about 40 pg, less than about 50 pg, less than about 100 pg, less than about 200 pg, less than about 500 pg, less than about 1 ng, less than about 5 ng, less than about 10 ng, less than about 20 ng, less than about 30 ng, less than about 40 ng, less than about 50 ng, less than about 100 ng, less than about 200 ng, less than about 500 ng, less than about 1000 ng, or any amount in between the number of cell-free nucleic acid molecules.

[0055] In some embodiments, the disclosure includes methods and systems for filling a sample with an amount of filler DNA to generate a mixture sample, the mixture sample comprising at least about 50 ng, 55 ng, 60 ng, 65 ng, 70 ng, 75 ng, 80 ng, 85 ng, 90 ng, 95 ng, 100 ng, 120 ng, 140 ng, 160 ng, 180 ng, 200 ng, or any amount between the total amount of nucleic acid mixture. In some embodiments, the filler DNA comprises at least about 5%, 10%, 15%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or 100% methylated filler DNA, with the remainder being unmethylated filler DNA, preferably 5%-50%, 10%-40%, or 15%-30% methylated filler DNA. In some embodiments, the mixture sample comprises filler DNA in an amount of 20 ng to 100 ng, preferably 30 ng to 100 ng, more preferably 50 ng to 100 ng. In some embodiments, the cell-free DNA from the sample and the first amount of filler DNA together comprise at least 50 ng of total DNA, preferably at least 100 ng of total DNA.

[0056] In some embodiments, the filler DNA is 50 bp to 800 bp in length, preferably 100 bp to 600 bp in length, more preferably 200 bp to 600 bp in length. In some embodiments, the filler DNA is double-stranded. The filler DNA is double-stranded. For example, the filler DNA can be junk DNA. Also, the filler DNA can be endogenous or exogenous DNA. For example, the filler DNA is non-human DNA, and in a preferred embodiment, is λ DNA. As used herein, "λ DNA" refers to enterobacteria phage λ DNA. In some embodiments, the filler DNA is not aligned with human DNA.

[0057] In some embodiments, samples may be taken before and / or after treatment of a subject with a disease or disorder. Samples may be taken from a subject during treatment or treatment regimen. Multiple samples may be taken from a subject to monitor the effectiveness of treatment over time. Samples may be taken from subjects known or suspected to have a disease or disorder where clinical testing does not provide a definitive positive or negative diagnosis. Samples may be taken from subjects suspected of having a disease or disorder. Samples may be taken from subjects experiencing unexplained symptoms such as fatigue, nausea, weight loss, aches and pains, weakness or bleeding. Samples may be taken from subjects with explained symptoms. Samples may be taken from subjects at risk of developing a disease or disorder due to factors such as family history, age, hypertension or prehypertension, diabetes or prediabetes, overweight or obesity, environmental exposures, lifestyle risk factors (e.g., smoking, alcohol intake, or drug use), or the presence of other risk factors.

[0058] In some embodiments, a sample can be taken at a first time point and sequenced, and then another sample can be taken at a later time point and sequenced.Such methods can be used for long-term monitoring purposes, for example, to track the onset or progression of disease.In some embodiments, disease progression can be tracked before, after, or during the course of treatment to determine the effectiveness of treatment.For example, the methods described herein can be performed on a subject before and after medical treatment to measure disease progression or regression in response to medical treatment.

[0059] After obtaining a sample from a subject, the sample can be processed to generate a dataset indicative of a disease or disorder of the subject. For example, the presence or absence or quantitative evaluation of cell-free nucleic acid molecules (e.g., ctDNA molecules) of the sample at a panel of cancer-associated genomic loci or microbe-associated loci can indicate cancer in the subject. Processing the sample obtained from the subject can include (i) subjecting the sample to conditions sufficient to isolate, enrich or extract a plurality of cell-free nucleic acid molecules, and (ii) assaying the plurality of cell-free nucleic acid molecules to generate a dataset (e.g., nucleic acid sequence). In some embodiments, the plurality of cell-free nucleic acid molecules are extracted from the sample and subjected to sequencing to generate a plurality of sequencing reads.

[0060] In some embodiments, the cell-free nucleic acid molecule may include cell-free ribonucleic acid (cfRNA) or cell-free deoxyribonucleic acid (cfDNA). The cell-free nucleic acid molecule (e.g., cfRNA or cfDNA) may be extracted from the sample by various methods. The cell-free nucleic acid molecule may be enriched by a plurality of probes configured to enrich nucleic acid (e.g., RNA or DNA) molecules corresponding to a panel of cancer-associated genomic loci. The probes may have sequence complementarity with nucleic acid sequences from one or more of the panel of cancer-associated genomic loci. A panel of cancer-associated genomic loci may include at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, at least 11, at least 12, at least 13, at least 14, at least 15, at least 16, at least 17, at least 18, at least 19, at least 20, at least about 25, at least about 30, at least about 35, at least about 40, at least about 45, 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 100, or more, different cancer-associated genomic loci. The probes may be nucleic acid molecules (e.g., RNA or DNA) having sequence complementarity to the nucleic acid sequences (e.g., RNA or DNA) of one or more genomic loci (e.g., cancer-associated genomic loci). These nucleic acid molecules may be primers or enrichment sequences. Assaying a sample using probes selective for one or more genomic loci (e.g., cancer-associated genomic loci or microbe-associated loci) may involve the use of array hybridization, polymerase chain reaction (PCR) or nucleic acid sequencing (e.g., RNA sequencing or DNA sequencing).

[0061] Sequencing of Nucleic Acid Molecules The present disclosure provides methods and techniques for determining the sequence of nucleotide bases of one or more polynucleotides. A polynucleotide can be, for example, a nucleic acid molecule, such as deoxyribonucleic acid (DNA) or ribonucleic acid (RNA), including variants or derivatives thereof (e.g., single-stranded DNA). Sequencing can be performed by a variety of currently available systems, such as, but not limited to, sequencing systems by Illumina®, Pacific Biosciences (PacBio®), Oxford Nanopore®, or Life Technologies (Ion Torrent®). Additionally, any sequencing method that indicates fragment length, such as paired-end sequencing, can be utilized. Alternatively or additionally, sequencing can be performed using nucleic acid amplification, polymerase chain reaction (PCR) (e.g., digital PCR, quantitative PCR, or real-time PCR), or isothermal amplification. Such systems can provide a plurality of raw genetic data corresponding to the genetic information of a subject (e.g., human) as generated by the system from a sample provided by the subject. In some examples, such systems provide sequencing reads (also "reads" herein). A read may include a series of nucleic acid bases that correspond to the sequence of a sequenced nucleic acid molecule. In some situations, the systems and methods provided herein may be used in conjunction with proteomic information.

[0062] In some embodiments, the sequencing reads are obtained by next-generation sequencing or next-generation sequencing. In some embodiments, the sequencing method comprises CAncer Personalized Profiling by deep Sequencing (CAPP-Seq), a next-generation sequencing-based method used to quantify cancer circulating DNA (ctDNA). This method can be generalized to any cancer type known to have recurrent mutations and can detect one molecule of mutant DNA in 10,000 molecules of healthy DNA. In some embodiments, the sequencing method comprises cfMeDIP sequencing as described by Shen et al., sensitive tumor detection and classification using plasma cell-free DNA methylomes, (2018) Nature, incorporated herein in its entirety. In some embodiments, the sequencing comprises bisulfite sequencing.

[0063] In some embodiments, sequencing includes modification of the nucleic acid molecule or fragment thereof, for example, by ligating a barcode, a unique molecular identifier (UMI), or another tag to the nucleic acid molecule or fragment thereof. Ligating a barcode, UMI, or tag to one end of the nucleic acid molecule or fragment thereof can facilitate analysis of the nucleic acid molecule or fragment thereof after sequencing. In some embodiments, the barcode is a unique barcode (e.g., a UMI). In some embodiments, the barcode is not unique, and the sequence of the barcode can be used in conjunction with endogenous sequence information, such as the start and stop sequences of the target nucleic acid (e.g., the target nucleic acid is adjacent to the barcode, and the sequence of the barcode is associated with the sequence of the start and end of the target nucleic acid to generate a uniquely tagged molecule). The barcode, UMI, or tag can be a known sequence used to associate a polynucleotide or fragment thereof with an input or target nucleic acid molecule or fragment thereof. The barcode, UMI, or tag can include natural or non-natural (e.g., modified) nucleotides (e.g., as described herein). The sequence of the barcode may be included within the adapter sequence such that the sequence of the barcode may be included within the sequencing read. The sequence of the barcode may include at least 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, or more nucleotides in length. In some cases, the sequence of the barcode may be of sufficient length and sufficiently different from the sequence of another barcode to allow identification of a sample based on the sequence of the associated barcode. The sequence of the barcode, or a combination of sequences of the barcode, may be used to tag and subsequently identify an "original" nucleic acid molecule or a fragment thereof (e.g., a nucleic acid molecule or a fragment thereof present in a sample obtained from a subject). In some cases, the sequence of the barcode, or a combination of sequences of the barcode, may be used in conjunction with endogenous sequence information to identify the original nucleic acid molecule or a fragment thereof. For example, the sequence of the barcode, or a combination of sequences of the barcode, may be used in conjunction with endogenous sequences adjacent to the barcode, UMI, or tag (e.g., the start and end of the endogenous sequence).

[0064] Processing the nucleic acid molecule or a fragment thereof may include performing amplification of the nucleic acid. For example, any type of nucleic acid amplification reaction can be used to amplify the target nucleic acid molecule or a fragment thereof to generate an amplification product. Non-limiting examples of nucleic acid amplification methods include reverse transcription, primer extension, polymerase chain reaction (PCR), ligase chain reaction, asymmetric amplification, rolling circle amplification, and multiple displacement amplification (MDA). Examples of PCR include, but are not limited to, quantitative PCR, real-time PCR, digital PCR, emulsion PCR, hot start PCR, multiplex PCR, asymmetric PCR, nested PCR, and assembly PCR. Amplification of nucleic acids may include one or more reagents, such as one or more primers, probes, polymerases, buffers, enzymes, and deoxyribonucleotides. Amplification of nucleic acids may be isothermal or may include thermal cycling. and / or have the length of the endogenous sequence.

[0065] Methylation Profile The present disclosure provides a method, system and kit for generating a methylation profile of a subject having or suspected of having a disease / condition, and the methylation profile can be used to determine whether the subject has or is at risk of having a disease / condition. Before using cfMeDIP-seq, the sample disclosed herein is subjected to library preparation. Briefly, after end repair and A-tailing, the sample is ligated to a nucleic acid adaptor and digested with an enzyme. As described above in the sample section, the prepared library can be combined with a filler nucleic acid (e.g., filler λDNA) to minimize the effect of low abundance ctDNA in the prepared library and create a mixed sample. In some embodiments, when the disease / condition is a localized (non-metastatic) cancer, the amount of ctDNA is low and cannot be easily and accurately measured and quantified. The mixed sample is made at least about 50ng, 80ng, 100ng, 120ng, 150ng or 200ng and subjected to further enrichment.

[0066] The methods, systems, and kits described herein are useful in treating, but are not limited to, adrenal gland cancer, anal cancer, bile duct cancer, bladder cancer, bone cancer, brain / CNS tumors, breast cancer, Castleman's disease, cervical cancer, colon / rectal cancer, endometrial cancer, esophageal cancer, Ewing's family of tumors, eye cancer, gallbladder cancer, gastrointestinal carcinoid tumors, gastrointestinal stromal tumors (gist), gestational trophoblastic disease, Hodgkin's disease, Kaposi's sarcoma, kidney cancer, laryngeal and hypopharyngeal cancer, leukemia (acute lymphocytic, acute myeloid, chronic lymphocytic, chronic myeloid, chronic myelomonocytic), liver cancer, lung cancer (non-small cell, small cell, lung carcinoma), and other cancers. The present invention is applicable to a variety of cancers, including lymphoma, lymphoma of the skin, malignant mesothelioma, multiple myeloma, myelodysplastic syndromes, nasal and paranasal sinus cancer, nasopharyngeal cancer, neuroblastoma, non-Hodgkin's lymphoma, oral and oropharyngeal cancer, osteosarcoma, ovarian cancer, penile cancer, pituitary cancer, prostate cancer, retinoblastoma, rhabdomyosarcoma, salivary gland cancer, sarcoma - adult soft tissue cancer, skin cancer (basal and squamous cell, melanoma, Merkel cell), small intestine cancer, gastric cancer, testicular cancer, thymic cancer, thyroid cancer, uterine sarcoma, vaginal cancer, vulvar cancer, Waldenstrom's macroglobulinemia, and Wilms' tumor. In an embodiment, the cancer is head and neck squamous cell carcinoma.

[0067] Binders can be used to enrich mixed samples. In some embodiments, the binder is a protein that contains a methyl-CpG binding domain. One such exemplary protein is the MBD2 protein. As used herein, "methyl-CpG binding domain (MBD)" refers to a specific domain of proteins and enzymes that is approximately 70 residues long and binds to DNA that contains one or more symmetrically methylated CpGs. The MBDs of MeCP2, MBD1, MBD2, MBD4 and BAZ2 mediate binding to DNA, and in the case of MeCP2, MBD1 and MBD2, preferentially to methylated CpGs. The human proteins MECP2, MBD1, MBD2, MBD3 and MBD4 comprise a family of nuclear proteins related by the presence of each of the methyl-CpG binding domains (MBDs). Each of these proteins, except MBD3, can specifically bind to methylated DNA.

[0068] In other embodiments, the binding agent is an antibody, and capturing the cell-free methylated DNA includes immunoprecipitating the cell-free methylated DNA using the antibody. As used herein, "immunoprecipitation" refers to a technique in which an antigen (such as a polypeptide and a nucleotide) is precipitated from a solution using an antibody that specifically binds to that particular antigen. This process can be used to isolate and concentrate a particular protein or DNA from a sample, and requires that the antibody be bound to a solid substrate at some point in the procedure. Solid substrates include, for example, beads, such as magnetic beads. Other types of beads and solid substrates may also be used.

[0069] One exemplary antibody is 5-MeC antibody. For immunoprecipitation procedure, in some embodiments, at least 0.05 μg of antibody is added to the sample. Meanwhile, in more preferred embodiments, at least 0.16 μg of antibody is added to the sample. To confirm the immunoprecipitation reaction, in some embodiments, the method described herein further comprises adding a second amount of control DNA to the sample.

[0070] The enriched sample is further amplified, purified, and sequenced to generate a plurality of sequence reads. The plurality of sequence reads is analyzed to identify a plurality of differentially methylated regions (DMRs). In some embodiments, the plurality of DMRs comprises DMRs derived from cell-free nucleic acid molecules derived from peripheral blood leukocytes (PBLs). In some embodiments, the plurality of DMRs comprises at least about 750,000 non-overlapping nucleic acid fragment windows of about 300 bp. These fragments comprise 8 or more CpG islands. In some embodiments, the DMRs are identified by comparing sequence reads generated from a sample obtained from a patient having a disease / condition with sequence reads generated from a sample obtained from a healthy control. In some embodiments, the healthy control comprises the same set of risk factors for developing the disease / condition. In some embodiments, the plurality of DMRs comprises at least about 997 DMRs, about 941 hypermethylations in HNSCC and 56 hypomethylations in HNSCC (Table 5). Using the same disclosed approach herein, hypermethylated DMRs can be detected for different cancers (e.g., lung cancer, pancreatic cancer, colorectal cancer) and hypomethylated DMRs can be detected for different cancers. [Table 1] TIFF0007665659000002.tif230151TIFF0007665659000003.tif230149TIFF0007665659000004.tif229149TIFF000 7665659000005.tif230150TIFF0007665659000006.tif230149TIFF0007665659000007.tif231149TIFF00076656590 00008.tif231150TIFF0007665659000009.tif228150TIFF0007665659000010.tif231150TIFF0007665659000011.t if231150TIFF0007665659000012.tif231150TIFF0007665659000013.tif232151TIFF0007665659000014.tif230150 TIFF0007665659000015.tif231151TIFF0007665659000016.tif230151TIFF0007665659000017.tif230150TIFF000 7665659000018.tif230150TIFF0007665659000019.tif230150TIFF0007665659000020.tif230151TIFF00076656590 00021.tif230151TIFF0007665659000022.tif230150TIFF0007665659000023.tif230149TIFF0007665659000024.t if230153TIFF0007665659000025.tif230151TIFF0007665659000026.tif231149TIFF0007665659000027.tif123150

[0071] Genomic mutation profile The present disclosure provides methods, systems and kits for generating a mutation profile of a subject having or suspected of having a disease / condition, and the methylation profile can be used to determine whether the subject has or is at risk of having a disease / condition. The samples disclosed herein are subjected to library preparation and next-generation deep sequencing (e.g., CAPP-Seq). Multiple sequencing reads are generated and analyzed. In some embodiments, the deep sequencing can be configured to maximize the identification of genomic mutations associated with the disease / condition. For example, and not intended to be limiting, for head and neck squamous cell carcinoma (HNSCC), a panel of standard HNSCC driver genes can be included in the selector for CAPP-seq. Additionally, for lung cancer, a panel of lung cancer driver genes can be included in the selector for CAPP-seq. Additionally, for pancreatic cancer, a panel of pancreatic cancer driver genes can be included in the selector for CAPP-seq. In some embodiments, the sensitivity of ctDNA detection can be increased by including genes with no known driver effect in a particular cancer type in the selector for CAPP-seq.

[0072] In some embodiments, a relative measure of ctDNA abundance is calculated from the mean mutant allele fraction (MAF). In some embodiments, the mean MAF of mutations identifies subjects whose mutation profile includes at least about 0.01% to at least about 10%. The ctDNA fraction of the samples disclosed herein is at least about 0.01%, 0.02%, 0.03%, 0.04%, 0.05%, 0.06%, 0.07%, 0.08%, 0.09%, 0.1%, 0.15%, 0.2%, 0.5%, 1%, 1.5%, 2%, 2.5%, 3%, 3.5%, 4%, 4.5%, 5%, 5.5%, 6%, 6.5%, 7%, 7.5%, 8%, 8.5%, 9%, 9.5%, 10%, or any percentage therebetween.

[0073] In some embodiments, the generated mutation profile of the subject does not include mutation variants derived from cell-free nucleic acid molecules derived from PBL. In some embodiments, the mutation profile includes genetic polymorphisms such as missense mutations, nonsense mutations, deletion mutations, insertion mutations, duplication mutations, inversion mutations, frameshift mutations, or repeat expansion mutations. In some embodiments, the mutation profile may include mutation variants derived from a fraction of cell-free nucleic acid molecules of a particular size range.

[0074] Fragment length profile In some embodiments, the length of the ctDNA fragment is shorter than that of the cell-free nucleic acid molecule derived from healthy subjects.In some embodiments, the length of the ctDNA that comprises at least one mutation is shorter than that of the cell-free nucleic acid molecule that comprises corresponding reference allele.In some embodiments, the length of the ctDNA fragment that comprises at least one DMR is shorter than that of the fragment of the cell-free nucleic acid molecule that comprises corresponding genomic region.

[0075] In some embodiments, sequencing does not utilize bisulfite sequencing, as this would cause degradation of ctDNA fragments and prevent preservation of the ctDNA length distribution. In some embodiments, the fragment lengths of the ctDNA are at least 60-500 bp, 80-300 bp, 90-250 bp, 80-170 bp, or 100-150 bp. In some embodiments, the present disclosure provides enrichment of cell-free nucleic acid samples based on the selection of cell-free molecules of a particular size. In some embodiments, the multimodal analysis includes utilizing the mutation profiles and fragment length profiles described herein by selectively including multiple nucleic acid molecules in the mutation profile based on their fragment lengths. In some embodiments, the multimodal analysis includes utilizing the methylation profiles and fragment length profiles described herein by selectively including multiple nucleic acid molecules in the methylation profile based on their fragment lengths. In some embodiments, multimodal analysis involves utilizing a mutation profile, a methylation profile, and a fragment length profile together by selectively including multiple nucleic acid molecules in a mutation profile based on their fragment lengths, and by selectively including multiple nucleic acid molecules in a methylation profile based on their fragment lengths, respectively.

[0076] Methods and systems for detecting cancer, determining the tissue of origin of a tumor, and providing a prognosis - Patents.com The present disclosure provides methods and systems for determining whether a subject has or is at risk of having a disease, the methods and systems comprising: subjecting a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from the subject to sequencing to generate at least one of (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile; and processing the at least one profile to determine whether the subject has or is at risk of the disease with at least 80% sensitivity or at least about 90% specificity, wherein the cell-free nucleic acid sample comprises less than 30 ng / ml of the plurality of nucleic acid molecules. In some embodiments, the sensitivity is at least about 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between the numbers. In some embodiments, the specificity is at least about 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.1%, 99.2%, 99.3%, 99.4%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between the numbers.

[0077] In some embodiments, the methods and systems include subjecting a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from the subject to sequencing to generate at least two of (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile. The methods result in a sensitivity of at least about 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage in between. In some embodiments, the sensitivity when using two profiles is increased by at least about 0.5%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, or a percentage between any numbers compared to the sensitivity when using one profile. In some embodiments, the sensitivity when using three profiles is increased by at least about 0.5%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, or a percentage between any numbers compared to the sensitivity when using two profiles.

[0078] Further, the method results in a specificity of at least about 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.1%, 99.2%, 99.3%, 99.4%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between numbers. In some embodiments, the specificity when using two profiles is increased by at least about 0.5%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, or any percentage between numbers compared to the specificity when using one profile. In some embodiments, the specificity when three profiles are used is increased by at least about 0.5%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, or a percentage between any numbers compared to the specificity when two profiles are used.

[0079] The present disclosure provides methods and systems for processing a cell-free nucleic acid sample of a subject to determine whether the subject has or is at risk of having a disease, the methods and systems comprising obtaining the cell-free nucleic acid sample comprising a plurality of nucleic acid molecules; subjecting the plurality of nucleic acid molecules or derivatives thereof to sequencing to generate a plurality of sequencing reads; computer processing the plurality of sequencing reads to identify (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile for the plurality of nucleic acid molecules; and using at least the methylation profile, the mutation profile, and the fragment length profile to determine whether the subject has or is at risk of having the disease. In some embodiments, the method provides a sensitivity of at least about 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between the numbers. The method provides a specificity of at least about 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.1%, 99.2%, 99.3%, 99.4%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between the numbers.

[0080] The present disclosure provides methods and systems for determining the tissue origin of a tumor, comprising identifying a plurality of differentially methylated regions (DMRs), the plurality of DMRs being specific to a particular cancer (e.g., breast cancer, colon cancer, prostate cancer, HSNCC) and derived from a fraction of cell-free nucleic acid molecules. In some embodiments, the fraction of cell-free nucleic acid molecules is derived from ctDNA. In some embodiments, the method results in a sensitivity of at least about 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between numbers. The method results in a specificity of at least about 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.1%, 99.2%, 99.3%, 99.4%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or any percentage between numbers.

[0081] The present disclosure describes a method and system for obtaining a prognosis for a subject after undergoing treatment for a disease / condition. For example, the treatment includes surgical removal of tumor, chemotherapy, radiation therapy, or immunotherapy (e.g., TCR, CAR, etc.) designed for a particular type of cancer. In some embodiments, the method or system includes sequencing a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from the subject to generate at least one profile of (i) methylation profile, (ii) mutation profile, and (iii) fragment length profile, and monitor or detect minimal residual disease (MRD) based at least on at least one profile.

[0082] The present disclosure provides methods and systems for determining whether a subject has a disease / condition by assaying cell-free nucleic acid molecules from at least a portion of a sample from the subject, detecting a methylation level of at least a portion of the cell-free nucleic acid molecule that falls within a differentially methylated region (DMR) listed in Table 5, and comparing, using at least one computer processor, the methylation level detected in (b) with the methylation level of a corresponding portion(s) of the cell-free nucleic acid molecule that falls within a DMR listed in Table 5. In some embodiments, the methylation levels of at least about 6 or more, 10 or more, 15 or more, 20 or more, 30 or more, 40 or more, 50 or more, 60 or more, 70 or more, 80 or more, 90 or more, or 100 or more, 200 or more, 300 or more, 400 or more, 500 or more, 600 or more, or 700 or more DMRs listed in Table 5 are measured and compared to the methylation levels of the corresponding DMRs in healthy subjects as discussed herein.

[0083] Once a subject is accurately diagnosed and undergoes treatment to treat cancer, such as surgical removal, chemotherapy, or radiation therapy, it is important to monitor the effectiveness of the treatment and predict the patient's survival rate. In addition, it is important to detect minimal residual disease of cancer cells. The present disclosure provides a method and system for determining whether the subject has a high survival rate after undergoing treatment for a disease, the method and system comprising: assaying cell-free nucleic acid molecules from at least a portion of a sample from the subject; detecting a methylation level of at least a portion of the cell-free nucleic acid molecules contained in a differentially methylated region (DMR) listed in Table 6; and using at least one computer processor, comparing the methylation level detected in (b) with the methylation level of a corresponding portion(s) of the cell-free nucleic acid molecules contained in the DMRs listed in Table 6. In some embodiments, the DMRs listed in Table 6 represent regions associated with genes ZSCAN31, LINC01391, GATA2-AS1, STK3, and OSR1. [Table 2]

[0084] In some embodiments, the method further comprises adding a second amount of control DNA to the sample to confirm the immunoprecipitation reaction.

[0085] As used herein, a "control" can include both a positive control and a negative control, or at least a positive control.

[0086] In some embodiments, the method further comprises adding a second amount of control DNA to the sample to confirm capture of cell-free methylated DNA.

[0087] In some embodiments, identifying the presence of DNA from a cancer cell further comprises identifying the tissue of origin of the cancer cell.

[0088] In some instances, tumor tissue sampling may be difficult or involve significant risks, in which case it may be desirable to diagnose and / or subtype cancers without the need for tumor tissue sampling. For example, lung tumor tissue sampling may require invasive procedures such as mediastinoscopy, thoracotomy, or percutaneous needle biopsy. These procedures may result in the need for hospitalization, chest tubes, mechanical ventilation, antibiotics, or other medical interventions. Some individuals may not undergo the invasive procedures required for tumor tissue sampling due to medical comorbidities or due to priorities. In some instances, the actual procedure of tumor tissue procurement may depend on the suspected cancer subtype. In other instances, cancer subtypes may evolve over time within the same individual, and serial evaluations by invasive tumor tissue sampling procedures are often impractical and not well tolerated by patients. Thus, non-invasive cancer subtyping by blood tests may have many advantageous applications in the practice of clinical oncology.

[0089] Thus, in some embodiments, identifying the tissue of origin of cancer cells further comprises identifying the subtype of cancer.Preferably, the subtype of cancer distinguishes cancer based on stage (e.g., early lung cancer treated with surgery vs. late lung cancer treated with chemotherapy), histology (e.g., small cell carcinoma vs. adenocarcinoma vs. squamous cell carcinoma in lung cancer), gene expression pattern or transcription factor activity (e.g., ER status in breast cancer), copy number abnormality (e.g., HER2 status in breast cancer), specific rearrangement (e.g., FLT3 in AML), specific gene point mutation status (e.g., IDH gene point mutation), and DNA methylation pattern (e.g., MGMT gene promoter methylation in brain cancer).

[0090] In some embodiments, the comparison in step (f) is performed genome-wide.

[0091] In other embodiments, the comparison in step (f) is genome-wide and limited to specific regulatory regions, such as, but not limited to, FANTOM5 enhancers, CpG islands, CpG shores, CpG shelves, or any combination of the foregoing.

[0092] In some embodiments, the methods herein are for use in detecting cancer.

[0093] In some embodiments, the methods herein are for use in monitoring the treatment of cancer.

[0094] Data analysis system and method The methods and systems disclosed herein may include an algorithm or its use. One or more algorithms may be used to classify one or more samples from one or more subjects. One or more algorithms may be applied to data from one or more samples. The data may include biomarker expression data. The methods disclosed herein may include assigning a classification to one or more samples from one or more subjects. Assigning a classification to a sample may include applying an algorithm to a methylation profile, a mutation profile, and a fragment length profile. In some cases, at least one profile is input into a data analysis system that includes a trained algorithm for classifying a sample obtained from a subject as having a disease or minor damage.

[0095] The data analysis system may be a trained algorithm. The algorithm may include a linear classifier. In some cases, the linear classifier includes one or more of a linear discriminant analysis, a Fisher linear discriminant, a naive Bayes classifier, a logistic regression, a perceptron, a support vector machine, or a combination thereof. The linear classifier may be a support vector machine (SVM) algorithm. The algorithm may include a bidirectional classifier. The bidirectional classifier may include one or more decision trees, random forests, Bayesian networks, support vector machines, neural networks, or logistic regression algorithms.

[0096] The algorithm may include one or more of Linear Discriminant Analysis (LDA), Basic Perceptron, Elastic Net, Logistic Regression, (Kernel) Support Vector Machine (SVM), Diagonal Linear Discriminant Analysis (DLDA), Golbin Classifier, Parzen Base, (Kernel) Fisher Discriminant Classifier, k-Nearest Neighbors, Iterative RELIEF, Classification Tree, Maximum Likelihood Classifier, Random Forest, Nearest Centroid, Predictive Analysis of Microarrays (PAM), k-Means Clustering, Fuzzy C-Means Clustering, Gaussian Mixture Model, Gradient Response (GR), Gradient Boosting Method (GBM), Elliptic Net Logistic Regression, Logistic Regression, or combinations thereof. The algorithm may include a Diagonal Linear Discriminant Analysis (DLDA) algorithm. The algorithm may include a Nearest Centroid algorithm. The algorithm may include a Random Forest algorithm. In some embodiments, the performance of logistic regression, random forest, and gradient boosting method (GBM) is superior to that of linear discriminant analysis (LDA), neural network, and support vector machine (SVM) for distinguishing pre-eclampsia and non-pre-eclampsia.

[0097] kit The present disclosure provides a kit for identifying or monitoring a disease or disorder (e.g., cancer) in a subject. The kit may include a probe for identifying a quantitative measure (e.g., indicating the presence, absence, or relative amount) of a sequence in each of a panel of cancer-associated genomic loci of a sample of a subject. The quantitative measure (e.g., indicating the presence, absence, or relative amount) of a sequence in each of a panel of cancer-associated genomic loci of a sample of a subject may indicate a disease or disorder (e.g., cancer) in the subject. The probe may be selective for a sequence in a panel of cancer-associated genomic loci of a sample of a subject (e.g., DMRs listed in Tables 3, 5, and 6). The kit may include instructions for processing the sample using the probe to generate a data set indicating a quantitative measure (e.g., indicating the presence, absence, or relative amount) of a sequence in each of a panel of cancer-associated genomic loci of a sample of a subject.

[0098] The probes in the kit may be selective for sequences in a panel of cancer-associated genomic loci of the sample. The probes in the kit may be configured to selectively enrich for nucleic acid (e.g., RNA or DNA) molecules corresponding to the panel of cancer-associated genomic loci. The probes in the kit may be nucleic acid primers. The probes in the kit may have sequence complementarity with nucleic acid sequences from one or more of the panel of cancer-associated genomic loci or genomic regions. The panel of cancer-associated or microbe-associated genomic loci or genomic regions may include a panel of at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, at least 11, at least 12, at least 13, at least 14, at least 15, at least 16, at least 17, at least 18, at least 19, at least 20, or more different cancer-associated genomic loci or genomic regions.

[0099] The kit instructions may include instructions for assaying the sample using probes selective for sequences in a panel of cancer-associated genomic loci of the cell-free biological sample. The probes may be nucleic acid molecules (e.g., RNA or DNA) having sequence complementarity with nucleic acid sequences (e.g., RNA or DNA) from one or more of the panels of cancer-associated genomic loci. The nucleic acid molecules may be primers or enrichment sequences. The instructions for assaying the cell-free biological sample may include instructions for performing array hybridization, polymerase chain reaction (PCR), or nucleic acid sequencing (e.g., DNA sequencing or RNA sequencing) to process the sample and generate a data set indicative of a quantitative measure (e.g., indicative of presence, absence, or relative amount) of sequences in each of the panels of cancer-associated genomic loci of the sample. The quantitative measure (e.g., indicative of presence, absence, or relative amount) of sequences in each of the panels of cancer-associated genomic loci of the sample may be indicative of a disease or disorder (e.g., cancer).

[0100] The kit instructions may include instructions for measuring and interpreting assay readouts, which may be quantified at one or more of the panel of cancer-associated genomic loci to generate a data set that indicates a quantitative measure of the sequence at each of the panel of cancer-associated genomic loci of the sample (e.g., indicates the presence, absence, or relative amount). For example, quantification of array hybridization or polymerase chain reaction (PCR) corresponding to the panel of cancer-associated genomic loci may generate a data set that indicates a quantitative measure of the sequence at each of the panel of cancer-associated genomic loci of the sample (e.g., indicates the presence, absence, or relative amount). The assay readouts may include quantitative PCR (qPCR) values, digital PCR (dPCR) values, digital droplet PCR (ddPCR) values, fluorescence values, etc., or normalized values ​​thereof.

[0101] Computer Systems In some embodiments, certain steps are performed by a computer processor. The system and method can be implemented in various embodiments. An appropriately configured computing device, as well as associated communication networks, devices, software and firmware, can provide a platform for enabling one or more of the above-described embodiments. By way of example, FIG. 8 illustrates a general-purpose computing device 100 that can include a central processing unit ("CPU") 102 connected to a storage unit 104 and a random access memory 106. The CPU 102 can process an operating system 101, application programs 103, and data 123. The operating system 101, application programs 103, and data 123 can be stored in the storage unit 104 and loaded into the memory 106 as needed. The computing device 100 can further include a graphics processing unit (GPU) 122 operably connected to the CPU 102 and the memory 106, for offloading intensive image processing calculations from the CPU 102 and performing these calculations in parallel with the CPU 102. An operator 107 can interact with the computing device 100 using various input / output devices, such as a video display 108 connected by a video interface 105, and a keyboard 115, a mouse 112, and a disk drive or solid state drive 114 connected by an I / O interface 109. The mouse 112 can be configured to control cursor movement on the video display 108 and to operate various graphical user interface (GUI) controls displayed on the video display 108 with mouse buttons. The disk drive or solid state drive 114 can be configured to accept a computer readable medium 116. The computing device 100 can form part of a network via a network interface 111, enabling the computing device 100 to communicate with other appropriately configured data processing systems (not shown).One or more different types of sensors 135 may be used to receive input from a variety of sources.

[0102] The system and method may be implemented in virtually any manner of computing device, including a desktop computer, laptop computer, tablet computer, or wireless handheld. The system and method may also be implemented as a computer readable / usable medium that includes computer program code for enabling one or more computing devices to perform each of the various process steps in the method according to the invention. When multiple computing devices perform the entire operation, the computing devices are networked to distribute the various steps of the operation. It is understood that the term computer readable medium or computer usable medium includes one or more of any type of physical embodiment of the program code. In particular, the computer readable / usable medium may include program code embedded in one or more data storage units of a computing device, such as one or more portable storage products (e.g., optical disks, magnetic disks, tapes, etc.), memory associated with a computer and / or storage system.

[0103] As used herein, a "processor" may be any type of processor, for example, any type of general-purpose microprocessor or microcontroller (e.g., Intel™ x86, PowerPC™, ARM™ processor, etc.), a digital signal processing (DSP) processor, an integrated circuit, a field programmable gate array (FPGA), or any combination thereof.

[0104] As used herein, "memory" may include any suitable combination of any type of computer memory, either internally or externally located, such as, for example, random access memory (RAM), read only memory (ROM), compact disc read only memory (CDROM), electro-optical memory, magneto-optical memory, erasable programmable read only memory (EPROM), and electrically erasable programmable read only memory (EEPROM). Portions of memory 102 may be organized using a conventional file system and controlled and managed by an operating system that manages the overall operation of the device.

[0105] As used herein, a "computer-readable storage medium" (also referred to as a machine-readable medium, a processor-readable medium, or a computer usable medium having computer-readable program code embodied therein) is a medium capable of storing data in a format readable by a computer or machine. The machine-readable medium may be any suitable tangible, non-transitory medium, including magnetic, optical, or electrical storage media, including diskettes, compact disk read-only memories (CD-ROMs), memory devices (volatile or non-volatile), or similar storage mechanisms. The computer-readable storage medium may include various sets of instructions, code sequences, configuration information, or other data that, when executed, cause a processor to perform steps of methods according to embodiments of the present disclosure. Those skilled in the art will appreciate that other instructions and operations necessary to carry out the described embodiments may also be stored in the computer-readable storage medium. The instructions stored in the computer-readable storage medium may be executed by a processor or other suitable processing device and interface with circuitry to perform the described tasks.

[0106] As used herein, a "data structure" is a particular way of organizing computer data so that it can be used efficiently. A data structure may implement one or more particular abstract data types (ADTs), which specify the operations that may be performed on the data structure and the computational complexity of those operations. In comparison, a data structure is a concrete implementation of the specification provided by the ADTs.

[0107] The advantages of the present invention are further illustrated by the following examples. The examples and their specific details described herein are presented for illustrative purposes only and should not be construed as limiting the scope of the claims of the present invention.

[0108] [Example] material and method Obtaining peripheral blood leukocytes (PBLs) and plasma from HNSCC and healthy donors Patients diagnosed with HNSCC between 2014 and 2016 were identified from the prospective Anthology of Clinical Outcomes (Wong K. et al. 2010). All studies were approved by the Research Ethics Board of the University Health Network. HNSCC patient samples were obtained from the HNC Translational Research program at the Princess Margaret Cancer Centre based on the following criteria: 1) presentation of localized disease at diagnosis, 2) blood collection at diagnosis and at least one post-treatment time point, and 3) a minimum follow-up time of 2 years after diagnosis. All patients received treatment with curative intent, consisting of surgery with or without adjuvant radiotherapy. Healthy donors matched for age, sex, and current smoking status were identified from a prospective lung cancer screening program. 5–10 mL of blood was collected in ethylene-diamine-tetraacetic acid (EDTA) tubes. For HNSCC patients, blood was collected at diagnosis (baseline, BL) as well as 3 months after primary surgery (3M). When applicable, additional blood was drawn prior to adjuvant radiotherapy (PreRT), intermediate adjuvant radiotherapy (MidRT), and / or 12 months (12M) after primary surgery. Plasma was isolated from blood within 1 hour of collection and stored at -80°C until further processing. Peripheral blood leukocytes were also isolated from the same blood draws on HNSCC patients or healthy donors at the time of diagnosis.

[0109] cell culture The HPV-negative HNSCC cell line FaDu was kindly provided by Dr. Bradly Wouters (Princess Margaret Cancer Center) and was cultured in DMEM (Gibco) supplemented with 10% fetal bovine serum and 1% penicillin / streptomycin. FaDu cell cultures were incubated at 37°C in a humidified atmosphere containing 5% CO2. The identity of FaDu cells was confirmed by STR profiling. Cells were subjected to mycoplasma testing (e-MycoTMVALiD Mycoplasma PCR Detection Kit, Intron Bio) before use.

[0110] Isolation of cell-free DNA (cfDNA) and PBL genomic DNA (gDNA) cfDNA was isolated from total plasma using the QIAamp Circulating Nucleic Acid Kit (Qiagen) according to the manufacturer's instructions. Genomic DNA was isolated from PBLs, sheared to 150-200 base pairs using a Covaris M220 focused sonicator, and size-selected by AMPure XP magnetic beads (Beckman Coulter) to remove fragments >300 base pairs. Isolated cfDNA and sheared PBL genomic DNA were quantified by Qubit prior to library generation (Figure 9A and Figure 9B).

[0111] Sequencing library preparation 5–10 ng or 10–20 ng of DNA were used as input for cfMeDIP-seq or CAPP-seq, respectively. Input DNA was prepared for library generation using the KAPA HyperPrep Kit (KAPA Biosystems) with some modifications. Library adapters were utilized that incorporate random 2 bp sequences flanking both strands of the input DNA upon ligation followed by a constant 1 bp T sequence 5'. To minimize adapter dimerization during ligation, library adapters were added at a molar ratio of 100:1 adapter:DNA (~0.07 uM per 10 ng cfDNA) and incubated overnight at 4 °C for 17 h. After post-ligation cleanup, input DNA was eluted in 40 μL of elution buffer (EB, 10 mM Tris-HCl, pH 8.0–8.5) prior to library generation.

[0112] Generation of CAPP-seq libraries CAPP-seq library creation was performed as described in Newman et al. 2014 with some modifications. Libraries were PCR amplified for 10 cycles and up to 12 indexed amplified libraries were pooled together at 500–1000 ng. After addition of COT DNA and blocking oligos, the pooled libraries were SpeedVaced to evaporate all liquid and resuspended in 13 µL of resuspension mix (8.5 µL of 2X hybridization buffer, 3.4 µL of hybridization component A, 1.1 µL of nuclease-free water). Prior to hybridization, 4 µL of hybridization probe (i.e., HNSCC selector) was added to the resuspension mix for a total of 17 µL. After hybridization and PCR amplification / cleanup, libraries were eluted in 30 µL of IDTE, pH 8.0 (1x TE solution). The multiplexed libraries were sequenced on an Illumina NextSeq / NovaSeq / HiSeq4000 in 2x75 / 100 / 125 paired runs, respectively. The design of the HNSCC selector incorporated frequently recurrent genomic alterations in HNSCC from the COSMIC database as well as the E6 and E7 regions of the HPV-16 genome (Figure 11).

[0113] Alignment and quality control of CAPP-seq libraries The first two base pairs at each 5' end of unaligned paired reads, corresponding to the incorporated random molecular barcodes, were extracted and matched to generate 4-bp molecular identifiers (UMIs). A third T base pair spacer was also removed before alignment. Paired reads were aligned to the human genome (genome assembly GRCh37 / hg19) by BWA-mem, sorted and indexed by SAMtools (v 1.3.1), and recalibrated for base quality scores using the Genome Analysis ToolKit (GATK) BaseRecalibrator (v 3.8) according to best practices (see ). Overlapping sequences obtained from BAM files were collapsed based on their UMIs and labeled as singletons, single-stranded consensus sequences (SSCSs) or overlapping consensus sequences (DCSs) by ConsensusCruncher. 44 Quality control of each library was assessed by different metrics obtained from FastQC (Babraham Bioinformatics) as well as by different scripts to obtain capture efficiency (CollectHsMetrics, Picard 2.10.9), depth of coverage (DepthOfCoverage, GATK 3.8) and base pair position error rate (ides-bgreport.pl, Newman et al., 2016).

[0114] Detection of somatic nucleotide variants (SNVs) and quantification of ctDNA Removal of potential sequencing errors was performed by integrated digital error suppression (iDES) as described in Newman et al. 2016. Background polishing was performed by utilizing 20 healthy donor cfDNA samples as a training cohort (Figure 12). To prevent outliers from affecting downstream analysis, candidate SNVs in the bottom 15 or top 85 percentile of sequencing depth (<1500x, >5000x) across HNSCC cfDNA or PBL gDNA samples, as well as genes with average sequencing depth <500x, were excluded from the analysis. To account for clonal hematopoiesis, non-germline mutations were defined as having a mutant allele fraction of <10% in plasma. Candidate SNVs in HNSCC cfDNA samples were identified based on ≥3 criteria, read supporting duplex support and complete absence in matched PBL gDNA samples. The mutant allele fraction (MAF) of an identified SNV was calculated by dividing the number of reads corresponding to the alternative allele by the sum of the reads corresponding to the alternative and reference alleles. For each HNSCC cfDNA sample with an identifiable SNV, the average MAF across the SNV was calculated and used as a measure of ctDNA abundance. For cfDNA samples with only one identifiable SNV, the calculated MAF was used. Many of the detectable cancer-derived mutations may not be homozygous or clonal within the tumor, and for these reasons, the average MAF may be an underestimate of the true abundance of ctDNA in cell-free DNA.

[0115] Generation of cfMeDIP-seq libraries The cfMeDIP-seq protocol was performed as described in Shen et al. 2019 with modifications to the library preparation steps as described in "Sequencing Library Preparation". Multiplexed libraries were sequenced on an Illumina NextSeq / NovaSeq / HiSeq4000 in 2 × 75 / 100 / 125 paired runs, respectively. For generalizability, the cfMeDIP-seq library is described as any MeDIP-seq preparation method utilizing 5-10 ng of input DNA, regardless of source (i.e., cfDNA, gDNA).

[0116] Alignment and quality control of cfMeDIP-seq libraries For alignment and quality control of CAPP-seq libraries, unaligned paired reads were processed, aligned, sorted and indexed as described previously. Duplicate sequences from BAM files were collapsed by SAMtools. Quality control of each library was assessed by various metrics obtained from FastQC (Babraham Bioinformatics) and the R package MEDIPS (reference), including CpG coverage (MEDIPS.seqCoverage) and enrichment (MEDIPS.CpGenrich).

[0117] Selection of informative regions in cfMeDIP-seq profiles Fragments generated from paired reads of cfMeDIP-seq libraries were counted in non-overlapping 300 base pair windows by MEDIPS (MEDIPS.createSet), scaled by Reads Per Kilobase per Million (RPKM), and exported as WIG format (MEDIPS.exportWIG). WIG files from each sample were imported by R and collated as a matrix. The analysis was limited to cfDNA and PBL samples from 20 healthy donors to allow application in a non-disease context. Informative regions were based on the criteria of CpG density and correlation of RPKM values ​​between cfDNA and matched PBL. A sliding window based on CpG density (≥ n CpGs) was used to select a minimum threshold of ≥ 8 CpGs.

[0118] Calculation of absolute methylation from cfMeDIP-seq libraries Fragments from paired reads of cfMeDIP-seq libraries were counted as previously described in Selection of Informative Regions in cfMeDIP-seq Profiles and scaled to absolute methylation levels by the MeDEStrand R package. To calculate absolute methylation from counts, a logistic regression model was used to estimate the bias of DNA pulldown based on the density of CpGs (i.e., CpG density bias) (MeDEStrand.calibrationCurve). Based on the estimated CpG density bias, methylation within each window was corrected for fragments from positive and negative DNA strands. Windows with corrected fragments were log-transformed and scaled to values ​​between 0 and 1 to represent absolute methylation (MeDEStrand.binMethyl). Absolute methylation levels from each cfMeDIP-seq sample were exported as WIG-like files (i.e., WIG file format without track lines).

[0119] Design and performance evaluation of in silico PBL depletion To enrich for windows within the disease context, methylation from PBLs was removed by a process called "in silico PBL depletion". The analysis was restricted to PBL samples obtained from a cohort of 20 healthy donor samples to allow application in a non-cancer specific context. Our strategy for in silico PBL depletion was implemented as follows: 1. For each informative window described in the selection of informative regions in the cfMeDIP-seq profile, calculate the median absolute methylation across healthy donor PBL samples.

[0120] 2. Define the PBL depletion window based on the criterion of absolute median methylation <0.1.

[0121] 3. Restrict analysis of cfDNA samples inside the PBL depletion window.

[0122] The performance of the in silico PBL depletion strategy was evaluated by comparing the absolute methylation distribution of pre- and post-depletion PBL samples from a cohort of healthy donors used as a training set with an HNSCC cohort used as a validation set.

[0123] Differential methylation analysis To allow for robust detection of HNSCC-associated differentially methylated regions (DMRs), the analysis was restricted to HNSCC patients with detectable SNVs in plasma by CAPP-seq (n=20 / 32). Differential methylation analysis was restricted to informative regions after in silico PBL depletion. A collated matrix of binned fragment counts from HNSCC and healthy donor cfDNA samples, generated as previously described in the selection of informative regions in cfMeDIP-seq profiles, was utilized for identification of DMRs by the DESeq2 R package. Prefiltering was performed by removing regions with <10 counts across all cfDNA samples. A single factor defined as the condition (HNSCC vs. healthy donor) was used for contrast during differential methylation analysis. Briefly, differential methylation analysis was performed by scaling samples based on a size factor and variance estimate, followed by fitting a negative binomial general linear model. For each window, P values ​​were calculated between HNSCC and healthy donor status by Wald test. P values ​​within regions beyond the default Cook's distance cutoff were omitted from the calculation of adjusted P values ​​(Benjamini-Hochberg). Significantly hyper- or hypomethylated regions (hyper-DMR / hypo-DMR) in HNSCC cfDNA samples were defined as windows with adjusted P values ​​<0.1.

[0124] Enrichment of CpG signatures within cfDNA hypermethylated regions in HNSCC CpG features such as island, shore, shelf, and open sea (interCGI) are defined in the AnnotationHub R package (reference) (hg19_cpgs annotation). Using an in-house R package utilizing the “annotatr” and “GenomicRanges” R packages, we labeled the ID coordinates of each hypermethylated window (i.e., “chr.start.end”) within the PBL-depleted regions with overlapping CpG features (Figure 13).

[0125] To determine the enrichment probability of observed feature overlap versus the null distribution, 1000 random samplings were performed. For each sampling, an equal number of bins were selected based on the number of hypermethylation windows while maintaining the same distribution of CpGs. For each CpG feature across samplings, the observed overlaps were used to generate respective null distributions, which were subsequently transformed to a z-score scale. The overlaps of observed hypermethylated regions for each CpG feature were also z-score transformed, and informative statistics were derived from the null distributions. The estimated P-value of the observed overlap from the hypermethylation window was calculated as the number of random samplings with overlaps above / below the observed overlap of the null distribution.

[0126] Enrichment of cfDNA hypermethylated regions in HNSCC with cancer-specific hypermethylated cytosines from the Tumor Cancer Genome Atlas (TCGA) File information from publicly available hm450k profiles of all primary tumors from breast (BRCA), colorectal (COAD), head and neck (HNSC), prostate (PRAD), pancreatic (PAAD), lung adeno (LUAD), and lung squamous (LUSC) were downloaded from TCGA. Due to the majority of our HNSCC cohort presenting with oral tumors, files from the HNSC group were restricted to patients with primary sites in the "floor of the mouth" (n=55). An equal number of hm450k files were randomly selected from each of the remaining cancer types as well as from a separate database of healthy PBLs (GEO series GSE67393). A manifest of the downloaded files is provided in (Figure 14).

[0127] To generate "tumor-specific" hypermethylated cytosines, differential methylation analysis with limma was performed for each cancer type and compared individually with other cancer types as well as PBL (i.e., contrast agent). For a given contrast, a linear model is fitted for each probe-targeted cytosine incorporating the residual variance and sample beta value, and then a P value for the observed difference between contrasts is calculated by empirical Bayes smoothing. Hypermethylated cytosines with elevated methylation in a given cancer type for an individual comparison were defined by a log fold change of ≥ 0.25 and an adjusted P value (Benjamini-Hochberg) < 0.01. Hypermethylated cytosines specific to an individual cancer type were named "tumor-specific". For LUSC, LUAD, and PAAD cases, no or few tumor-specific hypermethylated cytosines were identified (0, 15, 18) and therefore were omitted from further analysis. For comparison with cfMeDIP-seq libraries, base pair positions from tumor-specific hypermethylated cytosines were overlaid with information windows after in silico PBL depletion, as described in Design of In-silico PBL Depletion and Evaluation of Performance.

[0128] Enrichment of overlap for HNSCC cfDNA hypermethylated regions with tumor-specific regions from TCGA was assessed by 10,000 random samplings using the same methodology described in Enrichment of CpG features within cfDNA hypermethylated regions in HNSCC.

[0129] Sensitivity and specificity of ctDNA detection by cfMeDIP-seq For cfMeDIP-seq libraries from a cohort of 32 HNSCC and 20 healthy donor cfDNA samples, ctDNA detection was defined based on the observation of the mean RPKM across HNSCC cfDNA hypermethylated regions of individual HNSCC cfDNA samples greater than the mean maximum RPKM across healthy donor cfDNA samples. The sensitivity and specificity of ctDNA detection based on this definition was evaluated by receiver operating characteristic (ROC) curve analysis. To minimize any confounding results due to the potential lack of ctDNA release in a subset of patients, ROC curve analysis was also performed only on 20 of the 32 HNSCC cfDNA samples with detectable ctDNA by CAPP-seq. Cross-validation was performed to evaluate the accuracy of ctDNA detection by DMR analysis. Briefly, CAPP-Seq positive patients and healthy donors were randomly assigned to a training set (60%, n=24) and a validation set (40%, n=16) while maintaining similar ctDNA abundance (determined by CAPP-Seq) between both sets. Excess DMRs were identified by differential methylation analysis between HNSCC and healthy donor samples in the training set. The sensitivity of ctDNA detection within these excess DMRs was evaluated as previously described in the validation set to obtain AUROC values ​​(Figure 2C). Random sampling was performed a total of 50 times.

[0130] Fragment length analysis of ctDNA detected by CAPP-seq and cfMeDIP-seq For each HNSCC cfDNA CAPP-seq library, we measured the median fragment lengths from all supporting paired reads of the specified SNVs (i.e., singletons, SCSs, DCSs) and paired reads containing the reference allele. When median fragment lengths were reported for patients with >1 SNV, we calculated the median across median fragment lengths from each SNV. For each HNSCC cfDNA, cfMeDIP-seq library, we calculated the median fragment length from all fragments mapping to previously determined HNSCC cfDNA hypermethylated regions. Due to the relative absence of methylation within the cohort of 20 healthy donors, fragment lengths for each healthy donor cfMeDIP-seq library were cross-checked prior to any calculations. In both types of libraries, fragment length analysis was limited to cfDNA within the first peak (i.e., <220 base pairs).

[0131] Enrichment of fragments (100-150 bp or 100-220 bp) within the excess DMR was calculated as follows: A null distribution of expected counts was generated from random 300 bp bins within the previously designed PBL depletion window, with identical counts and CpG density distributions, from a total of 30 samplings. Observed counts for each sample were determined based on read counts across the excess DMR. Enrichment was calculated for each sample based on the average observed count divided by the average expected count.

[0132] Supervised Hierarchical Clustering Prior to clustering, a pseudocount of 0.1 was added to all RPKM values ​​of cfMeDIP-seq libraries to allow for log2 transformation. Values ​​were scaled by Euclidean transformation and clustered by Ward's method. An arbitrary number of three distinct clusters were chosen (k=3), named methylation clusters 1–3, and used for subsequent analysis.

[0133] Metrics of ctDNA detection and quantification for clinical outcomes in HNSCC patients The potential clinical utility of ctDNA detection was evaluated by three metrics: 1) detection of SNVs by CAPP-seq, 2) detection of increased mean RPKM in hypermethylated regions by cfMeDIP-seq. For comparative analysis, patients were stratified based on the following criteria: 1) presence or absence of SNVs, 2) methylation cluster 1 vs. methylation cluster 2+3. Patient characteristics are listed in Table 1. [Table 3] TIFF0007665659000030.tif24235 [Table 4]

[0134] Cross-validation of ctDNA-derived methylation by cfMeDIP-seq analysis Receiver operating characteristic (ROC) curve analysis was performed to evaluate the robustness of cfMeDIP-seq for identifying ctDNA-derived methylation. To minimize confounding results due to low / absent ctDNA, the analysis was restricted to HNSCC patients with detectable ctDNA by CAPP-seq. The cfMeDIP-seq profiles of patients and healthy controls were split into a training set (HNSCC: n=12 / 20; healthy controls: n=12 / 20) and a test set (HNSCC: n=8 / 20; healthy controls: n=8 / 20). The training and test sets were balanced for ctDNA abundance as determined by CAPP-Seq analysis. A total of 50 splits were performed with ROC curve analysis performed at each iteration.

[0135] Identification of prognostic regions in HNSCC by TCGA analysis All available HNSCC cases from TCGA with matched legacy hm450k and RNA expression data were selected (n=520). Survival data were obtained from Jianfang et al.. For hm450k data, methylation was summarized in 300bp regions as previously described by calculating the average beta value among probe IDs for a given region. To identify regions hypermethylated in HNSCC primary tumors compared to adjacent normal tissues, independent Wilcoxon tests were performed for each region. Regions with an adjusted p-value <0.05 (Holms method) as well as a log fold change ≥1 in primary tumors compared to adjacent normal tissues were selected for subsequent analysis. To identify hypermethylated regions associated with prognosis, multivariate Cox regression was performed by selecting regions with p-value <0.05, taking into account age, sex, and clinical stage. Survival analysis was limited to a maximum follow-up period of 5 years after diagnosis, reflecting what was observed within the HNSCC cfDNA cohort. To further identify prognostic regions associated with altered gene expression, Spearman's correlation was calculated for each region's hm450k primary tumor profile against the matched RNA expression profiles of transcripts within a 2Kb window. Regions with absolute Rho values ​​>0.3 and false discovery rates <0.05 were selected, resulting in the final identification of five prognostic regions associated with expression of ZNF323 / ZSCAN31, LINC01395, GATA2-AS1, OSR1, and STK3 / MST2. For TCGA patient profiles, a composite methylation score (CMS) was obtained by calculating the sum of the β values ​​across all five prognostic regions. For cfMeDIP-seq profiles, the RPKM values ​​across all 943 excess DMRs were scaled to a sum of 1, and the CMS was obtained by calculating the sum of these scaled RPKM values ​​across all five prognostic regions.

[0136] Longitudinal monitoring of post-treatment plasma samples by cfMeDIP-seq cfMeDIP-seq libraries were successfully generated for 30 / 32 patients (Figure 17A-D). For the remaining two patients, insufficient material was isolated from plasma and / or did not meet quality metrics. ctDNA quantification of post-treatment cfMeDIP-seq libraries was performed as described previously and mean values ​​of RPKM across hypermethylated regions identified by differential methylation analysis were calculated. To facilitate interpretation, both pre- and post-treatment cfMeDIP-seq libraries were converted to percent DNA values ​​based on linear regression against the mean MAF calculated by matched CAPP-Seq profiles. To achieve high confidence detection of residual disease, a minimum ctDNA fraction of 0.2% was required in post-treatment samples, which corresponded to the maximum mean value of RPKM observed across all healthy controls.

[0137] Results and Discussion Multimodal profiling of cell-free DNA in localized HNSCC To investigate the ability of multimodal profiling to characterize ctDNA in the setting of localized cancer, we recruited 32 HNSCC patients into a prospective observational study in which peripheral blood samples were collected at consecutive time points (Figure 9A, Table 1). All patients were treated with surgery, and a subset received adjuvant radiotherapy (n=14) or chemoradiotherapy (n=11). With a median follow-up of 43.2 months, 9 / 32 patients (28%) had recurrence (Insurance Statistics 2-year recurrence-free survival: 88%).

[0138] As most patients exhibited a heavy smoking history, which is well documented to alter the genomic / epigenomic landscape of somatic tissues and contribute to precancerous lesions, we also analyzed blood samples from 20 risk-matched healthy donors previously enrolled in a lung cancer screening program. 34~37Plasma-derived cell-free DNA and PBL-derived genomic DNA (gDNA) were co-isolated from blood and subjected to quantification and analysis (Supplementary Figure 1A). In contrast to other studies that demonstrated significantly elevated levels of total plasma cell-free DNA in metastatic disease compared to healthy controls, no significant differences were observed between our HNSCC cohort and healthy donors. 38~41 (Supplementary Figure 1B).

[0139] Multimodal profiling of cell-free DNA and gDNA of PBLs from patients and healthy controls was performed (Figure 1). By subjecting the same samples to both mutational and methylome profiling, we were able to evaluate their contribution to tumor-naive detection and characterization of ctDNA. Mutations and methylation were profiled independently using CAncer Personalized Profiling by deep Sequencing (CAPP-Seq) and cell-free Methylated DNA ImmunoPrecipitation and high-throughput sequencing (cfMeDIP-seq), respectively. Furthermore, paired-end sequencing was utilized for both methodologies to obtain the length of the sequenced cell-free DNA fragments.

[0140] Tumor-Naive Detection of Mutation-Based ctDNA from Pre-Treatment Plasma We first evaluated approaches to improve the reliability of mutation-based ctDNA detection without confirmation within matched tumor samples. Recent studies have shown that genes frequently targeted for ctDNA detection, such as TP53, may harbor mutations derived from clonally expanded PBLs. Furthermore, because ctDNA contains both genetic and epigenetic features of tumors, we reasoned that orthogonal analysis of both features in patient cell-free DNA may increase the reliability of ctDNA detection. Therefore, to achieve tumor-naive detection of low abundance ctDNA with high confidence, we independently profiled mutations and methylation by CAPP-Seq and cfMeDIP-seq, respectively, for both cfDNA and matched PBLs.

[0141] To evaluate the sensitivity of ctDNA detection in HPV-negative HNSCC without prior knowledge from the tumor, we first measured mutation abundance in baseline plasma samples (Figure 2A). CAPP-Seq was performed using a sequencing panel designed to maximize the number of HNSCC-associated mutations (Table 3 and Figure 10). We also employed established error suppression methodologies to remove background base substitution errors. [Table 5] TIFF0007665659000033.tif231155TIFF0007665659000034.tif230154TIFF0007665659000035.tif232156TIFF0007665659000036.tif232155TIFF000 7665659000037.tif232155TIFF0007665659000038.tif230155TIFF0007665 659000039.tif232155TIFF0007665659000040.tif231155TIFF00076656590 00041.tif231155TIFF0007665659000042.tif234156TIFF0007665659000043.tif232155TIFF0007665659000044.tif232157TIFF0007665659000045.t if232154TIFF0007665659000046.tif233155TIFF0007665659000047.tif232156TIFF0007665659000048.tif231155TIFF0007665659000049.tif226154

[0142] We profiled plasma and PBL samples from HNSCC patients and healthy donors at diagnosis by CAPP-Seq, utilizing 10-30 ng of input DNA. To achieve sensitive detection of ctDNA at low abundance, we applied a CAPP-Seq selector optimized to maximize the number of detected mutations in HNSCC (Table 2 and Figure 10). We further improved analytical sensitivity by integrated digital error suppression (iDES), incorporating custom molecular barcodes and removing background base substitution errors identified within healthy donor plasma samples (Methods). [Table 6] TIFF0007665659000051.tif232152TIFF0007665659000052.tif214150

[0143] After selecting candidate somatic single nucleotide variants (SNVs) based on plasma profiling and removal of potential germline mutations, we characterized potential false positives due to clonal hematopoiesis (CH) by comparison with matched PBL profiles. Of 24 patients with identifiable candidate SNVs, 10 showed identical SNVs in matched PBL profiles with highly correlated mutant allele fractions (MAFs) (R=0.94, p=1.392e -07 , Fig. 2B). With the exception of PIK3CA, genes with these SNVs were unique to each patient (Fig. 2C). Because genes commonly affected by CH, such as DNMT3A, TET2, and ASXL1, were not included within the CAPP-Seq selectors, our findings of patient-specific SNVs within matched cfDNA and PBL samples further highlight the advantages of this approach over gene-level filtering. Plasma samples from four patients were strictly positive for SNVs derived from CH (Fig. 2D), suggesting that matched PBL profiling may significantly minimize false-positive detection of low-abundance ctDNA.

[0144] After removing candidate SNVs potentially reflecting CH, ctDNA was detected in the plasma of 20 patients (median [range]: 3 [1-10] SNVs per patient). To assess the validity of these SNVs, we compared our results with whole-exome sequencing data from 279 HNSCC tumors published by The Cancer Genome Atlas (TCGA). 45 We observed similarities in hypermutated genes, including TP53 (65% vs. 72%), PIK3CA (20% vs. 21%), FAT1 (15% vs. 23%) and NOTCH1 (10% vs. 19%) (Figure 2E). Interestingly, two patients showed single SNVs not found within these genes (GRIN3A and MYC, Figure 11), demonstrating the additional utility of gene profiling of unknown / non-driver effects to increase the sensitivity of ctDNA detection.

[0145] When ctDNA abundance was calculated based on the mean MAF of SNVs, ctDNA levels ranged from 0.14% to 4.83% (Figure 2F). This lower limit of detection is similar to that previously described by others utilizing tumor-naive CAPP-Seq analysis, estimated at approximately 0.14%. Including patients with undetectable ctDNA, the median ctDNA abundance across our HNSCC cohort was 0.49%, similar to that observed in localized NSCLC by CAPP-Seq.

[0146] Tumor-Naive Detection of Methylation-Based ctDNA from Baseline Plasma Next, we sought to define ctDNA-associated methylation patterns in HNSCC and healthy control samples. Because CAPP-Seq results showed the impact of false-positive mutations arising from PBLs, we reasoned that reduction of false-positive ctDNA-associated methylation could be achieved by removal of PBL-derived DNA methylation signals. Therefore, we used matched PBL MeDIP-seq profiles from HNSCC and healthy control samples to suppress their contribution to cell-free DNA methylation signals (Figure 3A) and assessed whether analysis of matched PBLs could also enable methylation-based ctDNA detection (Figure 3A). Pre-treatment HNSCC and healthy donor plasma as well as PBLs were profiled by cfMeDIP-seq utilizing 5-10 ng of input DNA. Methylation abundance was defined in non-overlapping 300-bp windows spanning chromosomes 1–22 ( n = 9,603,454 windows) with read counts normalized to reads per kilobase per million (RPKM) ( Methods ) as previously described.

[0147] Since the anti-5mC antibody utilized for methylation pull-down preferentially binds to DNA fragments with increasing CpG density, including CpG islands, we first characterized this interaction to identify regions likely to be highly represented within the cfMeDIP-seq data. We also applied MeDIP-seq to the HNSCC cell line FaDu to evaluate preferential binding of cancer-derived methylated DNA fragments. Comparing DNA fragment pull-down abundance (median RPKM) across windows with various numbers of CpGs, we observed increased enrichment up to 8 CpGs or more for both PBL and FaDu (Figures 12A and 12B). FaDu showed greater enrichment compared to PBL at 8 CpGs or more per 300 bp window. This result is consistent with the established phenomenon of CpG island hypermethylation in cancer cells, including FaDu. Based on these observations, we determined that windows with 8 or more CpGs (n=702,488) were most likely to be informative for ctDNA detection and were therefore utilized for all subsequent analyses.

[0148] For patients with localized cancer, the majority of plasma cell-free DNA is derived from PBLs. Therefore, we sought to utilize PBL MeDIP-seq profiles to bioinformatically suppress this contribution to cell-free DNA signals. We compared the RPKM values ​​for each window in cfMeDIP-seq profiles made from HNSCC and healthy donor cfDNA with MeDIP-seq profiles made from FaDu (1x1 comparison), unpaired PBLs (1x51 comparison) or paired PBLs (1x1 comparison). With PBLs being the main contributor to plasma cell-free DNA, genome-wide methylation profiles were highly correlated between plasma cell-free DNA and paired or unpaired PBLs (mode R=0.92 and R=0.91, respectively). The strength of these correlations likely reflects the contribution of PBLs to plasma cfDNA outside of their known size. In contrast, the correlation between plasma cell-free DNA and FaDu was weaker (modal R = 0.78) (Figure 3B).

[0149] To select a threshold for reduced methylation across PBLs while accounting for preferential pulldown, we scaled and normalized PBL cfMeDIP-seq profiles to absolute methylation levels (0-1) based on logistic regression modeling with the MeDEStrand R package (Methods). We selected 99,997 windows that exhibited a median absolute methylation of less than 0.1 across healthy donor PBLs. When these windows were applied to the excluded HNSCC PBLs, we observed an absolute methylation distribution similar to that of the utilized healthy donor PBLs (Figure 3B), demonstrating the generalizability of this approach. Similarly, none of these windows individually showed significantly higher methylation across HNSCC PBLs compared to healthy donor PBLs (Figure 3C and Figure 12B), limiting any sources of HNSCC-specific PBL methylation that may confound ctDNA detection. In other words, these results confirm that the major source of cfDNA methylation in both control and focally confined HPV-negative HNSCC plasma originated from PBLs and that bioinformatic removal of PBL-derived methylation could limit the signal that confounded ctDNA quantification.

[0150] Tumor-Naive Detection of Methylation-Based ctDNA Pretreatment To identify common ctDNA-derived hypermethylated regions within the HNSCC cohort, we performed differential methylation analysis comparing HNSCC patients and CAPP-Seq-detectable ctDNA (n=20) with healthy donors. Utilizing 99,994 300-bp windows depleted for methylation in PBL, we identified ctDNA-derived differentially methylated regions (DMRs) by comparing 20 HNSCC patients with CAPP-Seq-detectable ctDNA and 20 healthy controls. Overall, we identified 997 differentially methylated regions (DMRs) (hypermethylated: 941, hypomethylated: 56) across HNSCC samples (Figure 3C). Approximately half of the hypermethylated regions (hyperDMRs) were found to be directly adjacent to each other, with blocks of hypermethylation spanning lengths of up to 1800 base pairs (Figure 13A). These data suggest the presence of CpG islands within the identified over-represented DMRs. Conversely, adjacent hypomethylated regions (hypoDMRs) were not observed. Among the over-represented DMRs of 300 bp, 47.5% were present in contiguous blocks of hypermethylated signals extending up to 1800 bp in length (Figure 13A), and typically exhibited CpG islands ranging from 300 to 3000 bp in length. Indeed, CpG islands were significantly enriched for over-represented DMRs (Figure 3E). In contrast, CpG islands were significantly depleted for under-represented DMRs (Figure 13B).

[0151] To determine whether these over-represented DMRs are indeed enriched for CpG islands, we next assessed enrichment of over-represented DMRs for CpG islands, shore, shelf, and open sea by substitution analysis (Methods). As expected, significant enrichment of CpG islands, as well as significant depletion of shore and open sea, was observed within over-represented DMRs (Figure 3E). In contrast, under-represented DMRs were significantly enriched for open sea and depleted for CpG islands, in line with the hypomethylation of CpG-sparse regions frequently observed across cancers (Supplementary Figure 5B).

[0152] Finally, since methylation of specific regions may distinguish tissues of origin as previously described using cfMeDIP-seq, we also investigated whether the hyper-DMRs contain regions specific to HNSCC or other cancers. To identify tumor-specific methylation regions, we utilized HumanMethylation450K (hm450k) data generated from primary tumors provided by TCGA (Methods). Comparing primary tumors from invasive breast cancer (BRCA), colon adenocarcinoma (COAD), lung squamous cell carcinoma (LUSC), prostate adenocarcinoma (PRAD), HNSCC, pancreatic adenocarcinoma (PAAD), and PBL, we identified sufficient hypermethylated CpGs (≧50) specific to BRCA, COAD, PRAD, and HNSCC (Methods) (Figure 14). As expected, we observed a significant enrichment of plasma-derived DMRs overlapping with HNSC-specific hypermethylated CpGs, as well as a significant depletion of overlap across BRCA, COAD and PRAD-specific hypermethylated CpGs (Figure 3F), suggesting that the excess DMRs contain regions specific to HNSCC origin compared to a variety of other cancer types.

[0153] Mutation-based and methylation-based ctDNA detection are particularly concordant More studies have described ctDNA associated with reduced fragment lengths compared to healthy sources of plasma cell-free DNA, providing additional metrics for robust tumor-naive detection. As targeted sequencing has previously been shown to detect ctDNA at reduced fragment lengths, we first utilized our CAPP-Seq profiles to determine whether we might observe a similar trend in HNSCC patients. For each SNV identified per patient (Figure 2E), we measured the median fragment length including the SNV allele as well as the overlapping reference allele. For cases where multiple SNVs were identified in a patient's sample, we used the median across all SNVs and their reference alleles. In accordance with previous findings, we observed a consistent reduction in ctDNA fragment size compared to healthy cell-free DNA across patients (median [range] Δ=-17.5 [1-58] bp) (Figure 4A). There was no significant association between the mean MAF of these mutations and fragment length (Figure 15A).

[0154] Unlike bisulfite-based DNA methylation approaches, cfMeDIP-seq does not cause DNA degradation and therefore maintains the original fragment size distribution. This provides a novel opportunity to map DNA methylation and fragment length simultaneously. We assessed the distribution of fragment lengths within the plasma-derived excess DMRs previously identified for each patient. Due to the nature of these regions having low methylation across our healthy donors, DNA fragments across donors were combined for comparison. Similar to the mutation-based analysis, we observed a reduction in fragment lengths from 19 / 20 CAPP-Seq positive patients compared to grouped healthy controls (median [range] Δ=-7 [1-21] bp) (Figure 4B). This indicated a smaller reduction in fragment lengths compared to the mutation-based analysis, likely due to a partial contribution by healthy tissue of cell-free DNA fragments within the excess DMRs. Supporting this idea, samples with the shortest excess DMR fragments showed higher methylated ctDNA abundance (Pearson r = -0.64, p = 0.002) (Figure 15B). When we used the ratio of small fragments (100bp to 150bp) to large fragments (151bp to 220bp) for our excess DMRs, a previously described approach to enrich for ctDNA, we observed a similar trend of enrichment of ctDNA across the majority of CAPP-Seq-positive HNSCC samples (median [range] = 28 [-8 to 63]%) (Figure 4C).

[0155] To assess how the plasma cell-free DNA excess DMRs identified in our HNSCC cohort may differ across individuals within these small fragments (100-150 bp), we first performed hierarchical clustering. Four major clusters emerged utilizing the ConsensusClusterPlus R package, each with different levels of methylation across the excess DMRs (Figure 4E and Figure S16C). Similarly, the three clusters were defined by distinct ctDNA abundance as determined by CAPP-Seq (Figure S16D), suggesting a potential relationship between the mean excess DMR methylation and mutation-based ctDNA abundance.

[0156] Next, we investigated whether fragment lengths were concordant between ctDNA molecules identified by both CAPP-Seq and cfMeDIP-seq, potentially providing an additional layer of validation for our multimodal approach. To minimize the possibility that background DNA fragments could confound the calculated fragment lengths of ctDNA in cfMeDIP-seq profiles, we limited our analysis to patients (n=10 HNSCC patients) with above-median methylation levels across the over-represented DMRs. Surprisingly, ctDNA fragment lengths were highly concordant between paired CAPP-Seq and cfMeDIP-seq profiles for each patient (Figure 4C) (Pearson r=0.86, p=0.0016) (CAPP-Seq: 43 distinct mutations, cfMeDIP-seq: 941 over-represented DMRs).

[0157] To further characterize the relationship between excess DMR methylation levels and mutation-based ctDNA abundance, we compared the mean RPKM across the 941 excess DMRs to the mean MAF determined by CAPP-Seq for each patient. Similar to the trends observed between methylation clusters, we observed a significant positive correlation (Pearson correlation, R = 0.85, p = 5e-10) (Figure 4F). To evaluate the sensitivity of ctDNA detection within these excess DMRs by cfMeDIP-seq, we compared the mean RPKM between our HNSCC cohort and healthy donors. For CAPP-Seq-positive patients (n = 20), ctDNA detection was particularly concordant (AUC = 0.998), with a slight decrease in performance upon inclusion of CAPP-Seq-negative patients (n = 12) (AUC = 0.944) (Figure 4G). Cross-validation (n=50 sampling) across CAPP-Seq positive patients and healthy donors resulted in a median AUC of 0.984 (Figure 16A), demonstrating the robustness of the approach disclosed herein.

[0158] Based on these observations, we assessed whether ctDNA could be enriched within cfMeDIP-seq profiles by limiting the analysis to cell-free DNA fragments of shortened length. We assessed the proportion of cell-free DNA fragments within excess DMRs consisting of small (100-150 bp) fragments, as similar methods have been described to enrich for ctDNA using unmethylation-based approaches. Indeed, this resulted in ctDNA enrichment across the majority (median [range] = 28 [-8 to 63]%) of CAPP-Seq-positive HNSCC samples, but not for any of the healthy controls (Figure 4D). Thus, in silico size selection of cell-free DNA fragments may enrich for ctDNA within cfMeDIP-seq libraries and contribute to tumor-naive multimodal ctDNA analysis.

[0159] It has been previously described that in patients with localized non-metastatic cancer, detection of ctDNA by CAPP-Seq at diagnosis is associated with poor prognosis. Similarly, levels of ctDNA, as assessed by methylation of SHOX2 and SEPT9, are associated with poor prognosis in HNSCC. Therefore, we asked whether detection or quantification of ctDNA by CAPP-Seq and cfMeDIP-seq at diagnosis is associated with clinical outcomes within our HNSCC cohort. Indeed, detection of ctDNA by CAPP-Seq (i.e., CAPP-Seq positive vs. CAPP-Seq negative) (hazard ratio [HR] = 7.6, log-rank p = 0.026; Supplementary Figure 8D), as well as increased methylation within our previously identified excess DMRs (i.e., methylation cluster 1 + 2 + 3 vs. methylation cluster 4) (HR = 4.51, p = 0.038; Figure 4G) correlated with shorter survival. Consistent with this finding, the average RPKM across excess DMRs correlated with cancer stage (Supplementary Figure 8E).

[0160] Next, we compared the median fragment lengths of ctDNA identified by either mutation- or methylation-based profiling. To minimize the possibility that background DNA fragments might confound the calculated fragment lengths of ctDNA in cfMeDIP-seq profiles, we selected patients with high ctDNA abundance as defined by hierarchical clustering (i.e., methylation clusters 1 and 2, Figure 4D, Supplementary Figure 8A-B). With this approach, ctDNA fragment lengths were notably concordant between paired CAPP-Seq and cfMeDIP-seq profiles for each patient (R = 0.83, p = 0.0016) (Figure 4H), despite completely different genomic regions being represented using these two profiling techniques. Moreover, similar to our analysis using fragments of all lengths, we observed the same relationship between small fragment ratios by CAPP-Seq (R = -0.79, p = 0.0038) and ctDNA fragment lengths (Figure 4I).

[0161] These results suggest that the similar reduction in fragment lengths observed from ctDNA detected by CAPP-Seq and cfMeDIP-seq may be a result of intrinsic properties of the tumor rather than due to genomic region, and that utilizing shorter fragment lengths may contribute to more specific identification of ctDNA.

[0162] Application of multimodal ctDNA detection for prognostic determination To evaluate the potential clinical application of tumor-naive multimodal ctDNA analysis, we compared ctDNA with clinical outcomes in a HNSCC cohort. cfMeDIP-seq profiles based on fragment length information were strongly associated with MAF in matched CAPP-Seq profiles (Pearson r = 0.85, p = 3 × 10-9), suggesting that methylation intensity within the 941 excess DMRs indeed reflects ctDNA abundance (Figure 5C). Importantly, cross-validation analysis confirmed the robustness of these excess DMRs for detecting ctDNA (Figure 16C). Patients with ctDNA detected in baseline plasma by both mutation- and methylation-based methods (n=19) were significantly more likely to have advanced disease (i.e., stages III-IVA) (n=18 / 19) when compared to patients without detectable ctDNA (n=8 / 13) (Fisher's exact test p=0.028) and showed dramatically worse overall survival (hazard ratio [HR]=7.55, 95% confidence interval [CI]=[0.95-59.94], log-rank p=0.025) (Figure 5G). In comparison, stage alone failed to predict patients with worse overall survival (HR=2.59, 95% CI=[0.32-20.46], log-rank p=0.35) (Figure 16D), further demonstrating the potential clinical utility of multimodal ctDNA profiling.

[0163] Due to the known effect of DNA methylation on gene expression and resulting functional activity of cancer drivers, we reasoned that ctDNA methylation patterns at specific loci may have prognostic significance independent of ctDNA abundance. To assess whether our previously identified excess DMRs contain specific regions associated with prognosis independent of ctDNA abundance, we examined DNA methylation, RNA expression and clinical outcome data obtained by TCGA for all available HNSCC patients (n=520) (Figure 5C). First, we calculated the average value of β across all CpGs contained within different 300bp windows from the TCGA hm450k methylation array data. Restricting the analysis to probe hm450k regions (n=764 / 941) overlapping with plasma-derived excess DMRs, we identified 483 hypermethylated regions in primary tumors (n=520) compared with adjacent normal tissues (n=50) (Wilcoxon test, FDR<0.05, log2FC>1). We observed that several of these hypermethylated regions overlap or are located near CpGs within genes profiled by commercially available methylation-based ctDNA diagnostic tests, including SEPT9 and SHOX2, previously evaluated in HNSCC, as well as TWIST1 and ONECUT2 (Figure 17A). These results provide further evidence supporting the potential clinical relevance of our plasma-derived excess DMRs.

[0164] To further explore the potential clinical utility of these hypermethylated regions commonly held by our HNSCC cohort and TCGA HNSC hm450k profiles, we performed univariate Cox proportional hazards regression across all TCGA HNSCC patients with available hm450k profiles and disease-specific survival (DSS) outcomes (n=493 / 520). We identified 33 regions significantly associated with DSS (p<0.05). To further select prognostic regions likely to have a functional role in tumorigenesis, we compared the methylation levels of each region (n=33) with the expression of surrounding gene transcripts within 2kb. We then used the TCGA HNSCC cohort to identify a subset of 483 DMRs associated with (1) prognosis in multivariate Cox regression and (2) the expression of neighboring gene transcripts. Five regions were identified as fulfilling both criteria, and increased methylation in each region led to higher expression of ZNF323 / ZSCAN31, LINC01391, and GATA2-AS1 (Figure 5G, Figures S17A-S17C), and lower expression of STK3 / MST2 and OSR1 (Figure 5H) (Figure 5D). Regions associated with decreased and increased expression as a result of methylation were found to reside within the promoter or first exon / intron and gene body, respectively. We constructed a composite methylation score (CMS) from these five regions (Table 6) and stratified the TCGA HNSCC cohort according to this score (Figure 5E). Higher CMS was significantly associated with poorer survival outcome (HR=1.67, 95% CI=[1.25, 2.21], log-rank p=3.4×10 -4 ).

[0165] Finally, we also assessed whether CMS could provide similar prognostic information when applied to ctDNA. To enrich for ctDNA, we limited the analysis of cfMeDIP-seq libraries to fragments of 100–150 bp length as described above (Figure 4E). To account for the relative contribution of ctDNA methylation levels obtained by the five putative prognostic markers, we normalized the cfMeDIP-seq RPKM values ​​from these regions across the 941 excess DMRs. This resulted in a similar trend where higher CMS was marginally associated with poorer survival (log-rank p=0.1; HR=3.06) (Figure 5F), suggesting that increased methylation of these putative prognostic regions identified from TCGA may also be informative within the cfMeDIP-seq profile. Furthermore, these results highlight how plasma cell-free DNA methylome profiling can be leveraged in combination with existing multi-omic cancer databases for biomarker discovery.

[0166] Monitoring disease after final treatment by cfMeDIP-seq Because cfMeDIP-seq achieved sensitive and quantitative ctDNA detection in HNSCC patients, we reasoned that, similar to CAPP-seq, cfMeDIP-seq may also be capable of monitoring treatment-related changes in ctDNA abundance. To quantify the percent of ctDNA within post-treatment cfMeDIP-seq profiles, we applied a linear transformation of the mean RPKM across previously identified plasma-derived excess DMRs (n=941) and restricted the size of fragments between 100-150 bp to further enrich for ctDNA. A detection threshold of 0.2% ctDNA was calculated based on the maximum of the mean values ​​of RPKM observed across all healthy controls. For HNSCC patients (n=20) with CAPP-Seq positive in one or more available post-treatment samples, cfMeDIP-seq was performed utilizing 10 ng of input cfDNA.

[0167] Measuring the change in ctDNA abundance throughout treatment, we observed various kinetics indicating complete clearance (CC), partial clearance (PC; reduction >90%) or no clearance (NC) (Figure 6A Supplementary Figure 10). Of 18 eligible patients, 5 (28%) showed no clearance (Figure 6B). Patients with no clearance were more likely to experience disease recurrence compared to patients with complete or partial clearance (HR=8.73, 95%CI=[1.5, 50.92], log-rank p=0.0046) (Figure 6C). Interestingly, all patients with higher ctDNA abundance at the time of last sampling compared to the time of diagnosis showed disease recurrence. Moreover, the only patient within this group who did not demonstrate disease recurrence was lost to follow-up but died within 1 year of treatment for unknown reasons. Of the 13 patients with undetectable post-treatment ctDNA by cfMeDIP-seq, nine remained disease-free with a median follow-up of 44.4 months (minimum = 12.2, maximum = 58.7). Of the other four patients, one had persistent disease in regional lymph nodes and the other experienced a relapse 3.5 to 7.7 months (median 7.4 months) after last collection. Of note, these relapses among patients with undetectable post-treatment ctDNA were significantly delayed compared with the four relapses among patients with detectable post-treatment ctDNA (median [range]: 3.0 [1.7 to 5.2] months after last collection). Collectively, these results demonstrate that cell-free DNA methylome profiling of plasma by cfMeDIP-seq can be used to assess response to definitive treatment and to identify patients at high risk for rapid relapse.

[0168] Consideration The widespread implementation of ctDNA in clinical practice could be accelerated by methods that can be applied across patients and in the absence of tumor material. In the described study, we evaluated the ability of a multimodal genome-wide cell-free DNA profiling technique for tumor-naive detection of ctDNA within an exploratory cohort of HNSCC patients with low ctDNA. We show that incorporation of matched PBLs improves ctDNA detection using both mutations (i.e., CAPP-Seq) as well as DNA methylation (i.e., cfMeDIP-seq). Furthermore, by utilizing CAPP-Seq to stratify patients with detectable and undetectable ctDNA, we achieved robust identification of ctDNA-derived methylation patterns. We show for the first time that biophysical properties of plasma cell-free DNA that reflect tumor origin (i.e., shortened fragment length) are conserved across molecular aberrations and detection platforms. Tumor-naive ctDNA detection and quantification will find multiple clinical applications, investigating the prognostic associations of ctDNA abundance and methylation patterns.

[0169] Tumor-naive ctDNA detection currently encounters several limitations due to low ctDNA abundance. Recent studies have profiled paired PBL and / or healthy control plasma to identify mutations derived from clonal hematopoiesis, which are the main cause of false-positive detection of ctDNA. However, the incorporation of orthogonal metrics can further improve accuracy and clinical applicability. Here, we evaluated the ability of multimodal genome-wide cell-free DNA profiling techniques for tumor-naive ctDNA detection within a cohort of HNSCC patients with low ctDNA abundance. We demonstrated a high degree of concordance between ctDNA metrics (abundance and fragment length) detected by mutation-based and methylation-based profiling methods. Furthermore, we showed that tumor-naive multimodal ctDNA profiling may present value by identifying putative prognostic biomarkers independent of ctDNA abundance, as well as by monitoring ctDNA abundance in serial samples.

[0170] Tumor-naive detection of ctDNA has many practical advantages in both research and clinical settings. Recent studies have utilized matched tumor profiling for validation of ctDNA-derived regions identified at low abundance in early disease to improve sensitivity. However, one limitation of these approaches is the number of informative regions that are lost due to tumor sampling heterogeneity, which may be further exacerbated when applied to post-treatment ctDNA derived from previously unsampled subclones. Furthermore, the clinical benefits of these tumor-informed detection methods are limited to cancers that are readily accessible by biopsy, circumventing one of the main strengths of non-invasive liquid biopsies. By utilizing a tumor-naive multimodal profiling strategy, we achieved similar results in early cancers without the drawbacks of tumor-informed methods.

[0171] This is the first study to utilize mutation and methylation profiling for comprehensive detection of ctDNA from a cohort of patients with localized cancer. Extending this multimodal profiling approach to other cancer types and disease settings will be important for the continued development of liquid biopsy. Furthermore, while numerous ctDNA studies in HNSCC have been described utilizing detection methods based on mutation, methylation or HPV profiling, we describe the first application of a genome-wide mutation / methylation profiling method to identify previously known targets (i.e., TP53 mutations or SEPT9 / SHOX2 methylation) in addition to less well-explored / unexplored targets.

[0172] Tumor-naive detection of ctDNA has many practical advantages in both research and clinical settings. Although tumor mutation profiling can identify patient-specific markers for ctDNA detection at low abundance, such personalized approaches rely on high-purity tumor samples from cancer types with sufficient mutational burden. Mutation profiling for personalized assay design can be costly and time-consuming, and rarely accounts for genomic heterogeneity within primary tumors or across metastatic clones. Furthermore, ctDNA detection methods that rely on access to tumor tissue diminish an important advantage of noninvasive liquid biopsy. By integrating independent cell-free DNA characteristics, we have achieved highly sensitive detection of ctDNA in early-stage cancers without the drawbacks of tumor-informed methods.

[0173] In our analysis, we selected patients with detectable ctDNA by CAPP-Seq to identify ctDNA-derived methylation patterns using cfMeDIP-seq. This approach provided further validation of the tumor-derived nature of plasma cell-free DNA in our cohort. ctDNA methylation patterns could quantify ctDNA abundance in a similar manner to ctDNA mutations. Furthermore, methylation patterns revealed tumor origin and identified putative prognostic and dynamic biomarkers. The combination of CAPP-Seq and cfMeDIP-seq enabled detailed molecular characterization of low abundance ctDNA. Mutation-based ctDNA quantification contributed to the discovery of HNSCC-specific excess DMRs in plasma, some of which were confirmed to be prognostic even after adjusting for ctDNA abundance. Thus, simultaneous mutation and methylation profiling can complement each other by revealing quantitative, tissue-specific, and prognostic ctDNA biomarkers. Furthermore, methylome profiling may prove particularly useful in cancer types with few recurrent or clonal mutations.

[0174] Similar to previous studies, we also observed a decrease in ctDNA fragment length while comparing with cell-free DNA from healthy donors using both mutation- and methylation-based approaches. Unlike healthy cell-free DNA, which is consistently on average around 166-167 bp, the length of ctDNA between patients can be highly variable. Factors influencing ctDNA fragment length include position-dependent fragmentation. 49 , metastatic versus non-metastatic disease 73 , as well as the dysregulated dynamics of various intracellular / extracellular DNases involved in healthy cell-free DNA fragmentation. 74Interestingly, we observed high concordance between the fragment lengths of ctDNA identified by CAPP-Seq and cfMeDIP-seq in eligible patients, even though both techniques explore different regions and tumor-derived abnormalities. These compelling data provide further evidence for the relevance and reproducibility of plasma cell-free DNA fragmentation in cancer patients.

[0175] We observed that detectable ctDNA by CAPP-Seq or elevated ctDNA abundance by cfMeDIP-seq was associated with poor prognosis within our HNSCC cohort. These results support the finding that detection of ctDNA by methylation 56 , as well as increased abundance due to copy number aberrations 75 or HPV detection 76 This is in accordance with previous HNSCC ctDNA studies that identified patients at high risk, and there was an incomplete correlation with tumor stage, suggesting that other unmeasured features of tumor biology may contribute to ctDNA abundance.

[0176] To our knowledge, no study has previously identified prognostic regions in cell-free DNA of HNSCC, independent of ctDNA detection / abundance, possibly in part due to limitations of commonly used ctDNA detection methods. We demonstrated that cell-free DNA methylome profiles, in conjunction with TCGA data, can serve as a discovery tool to identify novel prognostic methylation biomarkers in HNSCC. A composite methylation score composed of five DMRs demonstrated consistent prognostic associations across methylation detection platforms (hm450k and cfMeDIP-seq) and biospecimen types (cell-free DNA of tumor tissue and plasma). Although larger cohorts are needed in the future to validate our findings, this study indicates that genome-wide identification of methylated regions by cfMeDIP-seq may enable the discovery of novel prognostic biomarkers.

[0177] The performance of cfMeDIP-seq was evaluated in relation to disease prognosis. By applying a stringent threshold of approx. 0.2% ctDNA or more after treatment as detectable disease, we were able to predict disease recurrence for 4 of 9 patients. For the remaining 5 patients who relapsed (n=4) or had persistent disease (n=1) who failed to have detectable ctDNA after treatment, we typically observed a longer time to recurrence, suggesting that the fraction of ctDNA at these time points may have been below the lower limit of detection for cfMeDIP-seq. In subsequent studies utilizing cfMeDIP-seq for monitoring tumor-naive disease, more frequent collection of plasma after treatment may help address these limitations.

[0178] Because we have demonstrated the potential clinical utility of multimodal profiling in localized disease and HNSCC, these methods will contribute to future biomarker discovery and ultimately clinical utility for patients with various cancer types. This study makes multiple notable contributions. First, it combines cell-free DNA mutation, methylation and fragment length analysis. Furthermore, we systematically profiled plasma samples and paired PBLs from both HNSCC patients and risk-matched healthy controls. These analyses revealed important insights regarding the optimal handling of multimodal profiling for ctDNA detection and characterization. For example, our unique approach of removing contributing methylation signals from leukocytes and enriching for tumor-derived methylation using fragment length signatures will prove useful for future studies.

[0179] In conclusion, we demonstrate that tumor-naive CAPP-Seq profiling of ctDNA allows for high-confidence identification of ctDNA-derived methylation by cfMeDIP-seq. Taking advantage of the strength of epigenetic profiling by cfMeDIP-seq, we further show that these ctDNA-derived methylated regions demonstrate potential as markers of tumor origin, prognosis, and treatment response. Incorporating some of the approaches we described for improved sensitivity of ctDNA detection by cfMeDIP-seq in HNSCC, such as PBL depletion window and limiting analysis to short fragments, may also be applied to a variety of other localized cancers for clinical benefit. The disclosed framework is broadly applicable to other clinical situations where tumor tissue availability is limited.

[0180] Although preferred embodiments of the invention have been described herein, those skilled in the art will recognize that modifications can be made to the invention without departing from the spirit of the invention or the scope of the appended claims. All documents disclosed herein, including those in the following references, are incorporated by reference. In one aspect, the present invention provides the following. [Item 1] 1. A method for detecting the presence of circulating tumor deoxyribonucleic acid (ctDNA) from cancer cells in a subject, comprising: (a) obtaining a cell-free deoxyribonucleic acid (DNA) sample from said subject; (b) subjecting the sample to library preparation to allow subsequent sequencing of the cell-free methylated DNA; (c) capturing the cell-free methylated DNA with a binding agent selective for methylated polynucleotides; (d) sequencing the captured cell-free methylated DNA; (e) computing the sequence of the captured cell-free methylated DNA with control cell-free methylated DNA sequences from healthy individuals and individuals with cancer; and (f) identifying said presence of DNA from a cancer cell if there is a statistically significant similarity between one or more sequences of the captured cell-free methylated DNA and a cell-free methylated DNA sequence from an individual with cancer; In at least one of (d), (f) and (g), the subject's cell-free methylated DNA is restricted to a subpopulation according to fragment length metrics. [Item 2] 2. The method of claim 1, further comprising adding a first amount of filler DNA to the sample, wherein at least a portion of the filler DNA is methylated, and then optionally further denaturing the sample. [Item 3] 2. The method of claim 1, wherein the fragment length metric is the fragment length. [Item 4] 3. The method of claim 2, wherein the subject's cell-free methylated DNA is limited to fragments having a length of <170 base pairs (bp), <165 bp, <160 bp, <155 bp, <150 bp, <145 bp, <140 bp, <135 bp, <130 bp, <125 bp, <120 bp, <115 bp, <110 bp, <105 bp, or <100 bp. [Item 5] 3. The method according to item 2, wherein the subject cell-free methylated DNA is limited to fragments having a length of about 100 to about 150 bp, 110 to 140 bp, or 120 to 130 bp. [Item 6] 2. The method of claim 1, wherein the fragment length metric is the fragment length distribution of the subject's cell-free methylated DNA. [Item 7] 6. The method of claim 5, wherein the subject's cell-free methylated DNA is restricted to fragments within the bottom 50, 45, 40, 35, 30, 25, 20, 15, or 10 percentile based on length. [Item 8] 7. The method according to any one of items 1 to 6, wherein the subject's cell-free methylated DNA is further restricted to fragments within differentially methylated regions (DMRs). [Item 9] 8. The method according to any one of items 1 to 7, wherein the cell-free methylated DNA of the subject is further restricted during the capturing step. [Item 10] 8. The method according to any one of items 1 to 7, wherein the subject's cell-free methylated DNA is further restricted during the comparing step. [Item 11] 8. The method according to any one of items 1 to 7, wherein the limiting is during the identifying. [Item 12] 11. The method according to any one of items 1 to 10, wherein the sample is derived from blood or plasma of the subject. [Item 13] 12. The method of any one of items 1 to 11, wherein (f) comprises using a statistical classifier. [Item 14] Item 13. The method of item 12, wherein the classifier is derived by machine learning. [Item 15] 15. The method according to any one of items 1 to 14, wherein the control cell-free methylated DNA sequences from healthy individuals and cancer individuals are included in a database of differentially methylated regions (DMRs) between healthy individuals and cancer individuals. [Item 16] 16. The method of any one of items 1 to 15, wherein the control cell-free methylated DNA sequences from healthy individuals and cancer individuals are limited to the control cell-free methylated DNA sequences that are differentially methylated between healthy individuals and cancer individuals in DNA derived from cell-free DNA. [Item 17] 17. The method of claim 16, wherein the control cell-free methylated DNA sequence is differentially methylated in DNA derived from plasma between healthy individuals and cancer individuals. [Item 18] 18. The method according to any one of items 1 to 17, wherein the sample has less than 100 ng, 75 ng or 50 ng of cell-free DNA. [Item 19] 19. The method of any one of items 1 to 18, wherein the first amount of filler DNA comprises at least about 5%, 10%, 15%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or 100% methylated filler DNA, with the remainder being unmethylated filler DNA, preferably 5% to 50%, 10% to 40%, or 15% to 30% methylated filler DNA. [Item 20] 19. The method according to any one of items 1 to 18, wherein the first amount of filler DNA is between 20 ng and 100 ng, preferably between 30 ng and 100 ng, more preferably between 50 ng and 100 ng. [Item 21] 21. The method according to any one of items 1 to 20, wherein the cell-free DNA and the first amount of filler DNA obtained from the sample together comprise at least 50 ng of total DNA, preferably at least 100 ng of total DNA. [Item 22] 22. The method according to any one of items 1 to 21, wherein the filler DNA is 50 bp to 800 bp in length, preferably 100 bp to 600 bp in length, more preferably 200 bp to 600 bp in length. [Item 23] 23. The method of any one of items 1 to 22, wherein the filler DNA is double-stranded. [Item 24] 12. The method according to any one of items 1 to 11, wherein the filler DNA is junk DNA. [Item 25] 13. The method of any one of items 1 to 12, wherein the filler DNA is endogenous or exogenous DNA. [Item 26] 26. The method of claim 25, wherein the filler DNA is non-human DNA, preferably lambda DNA. [Item 27] 27. The method of any one of items 1 to 26, wherein the filler DNA does not align with human DNA. [Item 28] 28. The method according to any one of items 1 to 27, wherein the binding agent is a protein comprising a methyl-CpG binding domain. [Item 29] 29. The method according to any one of items 1 to 28, wherein the protein is MBD2 protein. [Item 30] 30. The method according to any one of items 1 to 29, wherein (d) comprises immunoprecipitating the cell-free methylated DNA using an antibody. [Item 31] 31. The method according to item 30, comprising adding at least 0.05 μg, preferably at least 0.16 μg, of the antibody to the sample for immunoprecipitation. [Item 32] 31. The method of claim 30, wherein the antibody is a 5-MeC antibody. [Item 33] 31. The method of claim 30, further comprising adding a second amount of control DNA to the sample after (c) to confirm the immunoprecipitation reaction. [Item 34] 33. The method of any one of items 1 to 32, further comprising adding a second amount of control DNA to the sample after (c) to confirm the capture of cell-free methylated DNA. [Item 35] 35. The method according to any one of items 1 to 34, wherein identifying the presence of DNA derived from a cancer cell further comprises identifying a tissue of origin of the cancer cell. [Item 36] 36. The method of claim 35, wherein identifying the tissue of origin of the cancer cells further comprises identifying a subtype of the cancer. [Item 37] 37. The method of claim 36, wherein the cancer subtype distinguishes the cancer based on stage, histology, gene expression pattern, copy number abnormality, rearrangement, or point mutation status. [Item 38] 38. The method according to any one of items 1 to 37, wherein (f) is performed genome-wide. [Item 39] 38. The method according to any one of items 1 to 37, wherein (f) is restricted to a specific regulatory region from a genome-wide perspective. [Item 40] 40. The method of claim 39, wherein the regulatory region is a FANTOM5 enhancer, a CpG island, a CpG shore, a CpG shelf, or any combination of the foregoing. [Item 41] 41. The method according to any one of items 1 to 40, wherein steps (f) and (g) are carried out by a computer processor. [Item 42] The cancers include adrenal gland cancer, anal cancer, bile duct cancer, bladder cancer, bone cancer, brain / CNS tumors, breast cancer, Castleman's disease, cervical cancer, colon / rectal cancer, endometrial cancer, esophageal cancer, Ewing's family of tumors, eye cancer, gallbladder cancer, gastrointestinal carcinoid tumors, gastrointestinal stromal tumors (gist), gestational trophoblastic disease, Hodgkin's disease, Kaposi's sarcoma, kidney cancer, laryngeal and hypopharyngeal cancer, leukemia (acute lymphocytic, acute myeloid, chronic lymphocytic, chronic myeloid, chronic myelomonocytic), liver cancer, lung cancer (non-small cell, small cell, pulmonary carcinoid tumors), lymphoma, lymphoma of the skin, malignant mesothelioma 42. The method according to any one of items 1 to 41, wherein the cancer is selected from the group consisting of: multiple myeloma, myelodysplastic syndromes, nasal cavity and paranasal sinus cancer, nasopharyngeal cancer, neuroblastoma, non-Hodgkin's lymphoma, oral cavity and oropharyngeal cancer, osteosarcoma, ovarian cancer, penile cancer, pituitary cancer, prostate cancer, retinoblastoma, rhabdomyosarcoma, salivary gland cancer, sarcoma - adult soft tissue cancer, skin cancer (basal and squamous cell, melanoma, Merkel cell), small intestine cancer, gastric cancer, testicular cancer, thymic cancer, thyroid cancer, uterine sarcoma, vaginal cancer, vulvar cancer, Waldenstrom's macroglobulinemia, Wilms' tumor. [Item 43] 42. The method according to any one of items 1 to 41, wherein the cancer is head and neck squamous cell carcinoma. [Item 44] 44. The method according to any one of items 1 to 43, for use in the detection of said cancer. [Item 45] 44. The method according to any one of items 1 to 43, for use in monitoring the treatment of cancer. [Item 46] 1. A method for determining whether a subject has or is at risk for having a disease, comprising: (a) subjecting a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from the subject to sequencing to generate at least one profile selected from the group consisting of (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile; and (b) processing the at least one profile to determine whether the subject has or is at risk for the disease with at least 80% sensitivity or at least about 90% specificity, wherein the cell-free nucleic acid sample comprises less than 30 nanograms (ng) / milliliter (ml) of the plurality of nucleic acid molecules. [Item 47] 47. The method of claim 46, wherein the cell-free nucleic acid sample comprises less than 10 ng / ml of the plurality of nucleic acid molecules. [Item 48] 47. The method of claim 46, wherein the cell-free nucleic acid sample comprises less than 5 ng / ml of the plurality of nucleic acid molecules. [Item 49] 47. The method of claim 46, wherein the cell-free nucleic acid sample comprises less than 1 ng / ml of the plurality of nucleic acid molecules. [Item 50] Item 47. The method of item 46, wherein the step of subjecting to (a) generates at least two profiles selected from the group consisting of (i), (ii) and (iii). [Item 51] 51. The method of claim 50, wherein the at least two profiles comprise the methylation profile and the fragment length profile. [Item 52] 51. The method of claim 50, wherein the at least two profiles comprise the mutation profile and the fragment length profile. [Item 53] 51. The method of claim 50, wherein the at least two profiles comprise the methylation profile and the mutation profile. [Item 54] 47. The method of claim 46, wherein the subjecting step (a) generates the methylation profile, the mutation profile, and the fragment length profile. [Item 55] 1. A method of processing a cell-free nucleic acid sample from a subject to determine whether the subject has or is at risk for having a disease, comprising: (a) obtaining the cell-free nucleic acid sample comprising a plurality of nucleic acid molecules; (b) subjecting the plurality of nucleic acid molecules or derivatives thereof to sequencing to generate a plurality of sequencing reads; (c) computationally processing the plurality of sequencing reads to identify (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile for the plurality of nucleic acid molecules; and (d) using at least said methylation profile, said mutation profile and said fragment length profile to determine whether said subject has or is at risk of having said disease. [Item 56] 56. The method of any one of items 46 to 55, wherein the disease comprises cancer. [Item 57] The cancers include adrenal gland cancer, anal cancer, bile duct cancer, bladder cancer, bone cancer, brain / CNS tumors, breast cancer, Castleman's disease, cervical cancer, colon / rectal cancer, endometrial cancer, esophageal cancer, Ewing's family of tumors, eye cancer, gallbladder cancer, gastrointestinal carcinoid tumors, gastrointestinal stromal tumors (gist), gestational trophoblastic disease, Hodgkin's disease, Kaposi's sarcoma, kidney cancer, laryngeal and hypopharyngeal cancer, leukemia (acute lymphocytic, acute myeloid, chronic lymphocytic, chronic myeloid, chronic myelomonocytic), liver cancer, lung cancer (non-small cell, small cell, pulmonary carcinoid tumors), lymphoma, lymphoma of the skin, malignant mesothelioma, multiple myeloma, myelodysplastic syndromes ... 57. The method of claim 56, wherein the cancer is selected from the group consisting of adult soft tissue cancer, nasal and paranasal sinus cancer, nasopharyngeal cancer, neuroblastoma, non-Hodgkin's lymphoma, oral cavity and oropharyngeal cancer, osteosarcoma, ovarian cancer, penile cancer, pituitary cancer, prostate cancer, retinoblastoma, rhabdomyosarcoma, salivary gland cancer, sarcoma - adult soft tissue cancer, skin cancer (basal and squamous cell, melanoma, Merkel cell), small intestine cancer, gastric cancer, testicular cancer, thymic cancer, thyroid cancer, uterine sarcoma, vaginal cancer, vulvar cancer, Waldenstrom's macroglobulinemia, Wilms' tumor, squamous cell carcinoma, and head and neck squamous cell carcinoma. [Item 58] 58. The method of claim 57, wherein the cancer is squamous cell carcinoma. [Item 59] 59. The method of claim 58, wherein the cancer is head and neck squamous cell carcinoma. [Item 60] 57. The method of any one of items 46 to 56, wherein the plurality of cell-free nucleic acid molecules comprises circulating tumor nucleic acid molecules. [Item 61] 61. The method of item 60, wherein the circulating tumor nucleic acid comprises circulating tumor DNA. [Item 62] 61. The method of item 60, wherein the circulating tumor nucleic acid comprises circulating tumor RNA. [Item 63] 63. The method of any of items 46 to 62, wherein the methylation profile comprises a plurality of differentially methylated regions (DMRs). [Item 64] 64. The method of claim 63, wherein the multiple DMRs are derived from ctDNA. [Item 65] 64. The method of claim 63, wherein multiple DMRs derived from peripheral blood leukocytes are removed from the methylation profile. [Item 66] 64. The method of claim 63, wherein the plurality of DMRs comprises at least about 56 genomic regions having a hypomethylation level compared to corresponding genomic regions from a normal healthy subject. [Item 67] 55. The method of claim 54, wherein the plurality of DMRs comprises at least about 941 genomic regions having hypermethylation levels compared to corresponding genomic regions from normal healthy subjects. [Item 68] 64. The method of item 63, wherein the DMR comprises a size of at least about 300 bp. [Item 69] 70. The method of claim 68, wherein the DMR comprises a size of at least about 100 bp to at least about 200 bp. [Item 70] 70. The method of claim 68, wherein the DMR comprises a size of at least about 100 bp to at least about 150 bp. [Item 71] 64. The method of item 63, wherein the DMR comprises at least 8 CpG genomic islands. [Item 72] 68. The method of any of items 66 or 67, wherein said normal healthy subject comprises the same set of risk factors as said subject. [Item 73] 73. The method of any of items 45 to 72, wherein the mutation profile comprises missense mutants, nonsense mutants, deletion mutants, insertion mutants, duplication mutants, inversion mutants, frameshift mutants, or repeat expansion mutants. [Item 74] 73. The method according to any of items 45 to 72, wherein any mutations present in a genomic DNA sample obtained from a plurality of peripheral blood leukocytes, said plurality of peripheral blood leukocytes being obtained from said subject, are removed from said mutation profile. [Item 75] 73. The method according to any of items 45 to 72, wherein any variants derived from clonal hematopoiesis are removed from the mutation profile. [Item 76] 76. The method of claim 75, wherein the mutation profile does not include mutations of genes DNMT3A, TET2, or ASXL1. [Item 77] 76. The method of claim 75, wherein the mutation profile does not include canonical cancer driver genes. [Item 78] 76. The method of item 75, wherein the mutation profile comprises a non-canonical cancer driver gene, and the non-canonical gene is GRIN3A or MYC. [Item 79] 79. The method of any of items 46 to 78, comprising selecting cell-free nucleic acid molecules based on the fragment length profile being at least about a fragment length range of 80 bp to 170 bp. [Item 80] 79. The method of any of items 46 to 78, comprising selecting cell-free nucleic acid molecules based on the fragment length profile being a fragment length range of at least about 100 bp to 150 bp. [Item 81] 81. The method of any of items 79 or 80, wherein the circulating tumor nucleic acid molecules are enriched. [Item 82] 82. The method of any of items 46 to 81, further comprising mixing the cell-free nucleic acid sample with filler DNA molecules to produce a DNA mixture. [Item 83] 83. The method of item 82, wherein the filler DNA molecule comprises a length of about 50 bp to 800 bp. [Item 84] 83. The method of item 82, wherein the filler DNA molecule comprises a length of about 100 bp to 600 bp. [Item 85] 83. The method of claim 82, wherein the filler DNA molecules comprise at least about 5% methylated filler DNA molecules. [Item 86] 83. The method of claim 82, wherein the filler DNA molecules comprise at least about 20% methylated filler DNA. [Item 87] 83. The method of claim 82, wherein the filler DNA molecules comprise at least about 30% methylated filler DNA. [Item 88] 83. The method of claim 82, wherein the filler DNA molecules comprise at least about 50% methylated filler DNA. [Item 89] 89. The method of any of items 46 to 88, further comprising incubating the DNA mixture with a binding agent configured to bind to methylated nucleotides to produce an enriched sample. [Item 90] 90. The method of claim 89, wherein the binding agent comprises a protein comprising a methyl-CpG binding domain. [Item 91] 90. The method of claim 89, wherein the protein is MBD2 protein. [Item 92] 90. The method of claim 89, wherein the binding agent comprises an antibody. [Item 93] 90. The method of claim 89, wherein the antibody is a 5-MeC antibody. [Item 94] 90. The method of claim 89, wherein the antibody is a 5-hydroxymethylcytosine antibody. [Item 95] 95. The method according to any of items 46 to 94, wherein said sequencing does not comprise bisulfite sequencing. [Item 96] 95. The method according to any of items 46 to 94, wherein the cell-free nucleic acid sample comprises a blood sample. [Item 97] 97. The method of claim 96, wherein the blood sample comprises a plasma sample. [Item 98] 98. The method according to any one of items 46 to 97, further comprising detecting the origin of the cancer tissue. [Item 99] 98. The method of any of items 46 to 97, further comprising generating a report comprising a prognosis of survival of the subject. [Item 100] 98. The method of any of items 46 to 97, further comprising administering a treatment to the subject. [Item 101] 98. The method of any of items 46 to 97, further comprising, following treatment of the disease, providing a second report indicating whether the treatment is effective. [Item 102] 1. A method for determining whether a subject has or is at risk for having a condition, comprising: (a) assaying cell-free nucleic acid molecules obtained from at least a portion of a sample from said subject; (b) detecting the methylation level of at least a portion of said cell-free nucleic acid molecules that are contained in a differentially methylated region (DMR) listed in Table 5; and (c) comparing, using at least one computer processor, the methylation level detected in (b) with the methylation level of a corresponding portion(s) of the cell-free nucleic acid molecule that is contained in a DMR listed in Table 5. [Item 103] 103. The method of claim 102, wherein the cell-free nucleic acid molecule comprises ctDNA. [Item 104] 103. The method of claim 102, comprising performing said sequencing analysis, wherein said sequencing analysis comprises cell-free Methylated DNA ImmunoPrecipitation (cfMeDIP) sequencing. [Item 105] 103. The method of claim 102, wherein the detecting step comprises measuring a methylation level of at least a portion of the nucleic acid molecule contained in 6 or more, 10 or more, 15 or more, 20 or more, 30 or more, 40 or more, 50 or more, 60 or more, 70 or more, 80 or more, 90 or more, or 100 or more DMRs listed in Table 5. [Item 106] 1. A method for determining whether a subject has a higher survival rate after undergoing treatment for a disease, comprising: (a) assaying cell-free nucleic acid molecules obtained from at least a portion of a sample from said subject; (b) detecting the methylation level of at least a portion of said cell-free nucleic acid molecules that are contained in a differentially methylated region (DMR) listed in Table 6; and (c) processing, using at least one computer processor, the methylation level detected in (b) with the methylation level of a corresponding portion(s) of the cell-free nucleic acid molecule that is contained in a DMR listed in Table 6. [Item 107] 107. The method of claim 106, wherein the cell-free nucleic acid molecule comprises ctDNA. [Item 108] 107. The method of claim 106, wherein the detecting step comprises presenting a complex methylation score (CMS). [Item 109] Item 108. The method of item 107, wherein the CMS comprises the sum of the β values ​​of the DMRs listed in Table 6. [Item 110] The method of item 107, wherein a higher CMS indicates a lower survival rate of the subject. [Item 111] The method of claim 107, wherein the CMS is independent of the abundance of ctDNA. [Item 112] 112. The method according to any of items 102 to 111, wherein the disease is squamous cell carcinoma. [Item 113] Item 113. The method of item 112, wherein the cancer is head and neck squamous cell carcinoma. [Item 114] 114. The method of any one of items 102 to 113, further comprising selecting the cell-free nucleic acid molecules based on a fragment length range of at least about 80 bp to 170 bp. [Item 115] 1. A system for determining whether a subject has or is at risk of having a disease, comprising: subjecting a plurality of nucleic acid molecules from a cell-free nucleic acid sample obtained from said subject to sequencing to generate at least one of: (i) a methylation profile; (ii) a mutation profile; and (iii) a fragment length profile; and processing said at least one profile to determine whether said subject has or is at risk for said disease with at least 80% sensitivity or at least about 90% specificity, wherein said cell-free nucleic acid sample comprises less than 30 ng / ml of said plurality of nucleic acid molecules. A system comprising one or more computer processors individually or collectively programmed to carry out a process comprising: [Item 116] 1. A system for processing a cell-free nucleic acid sample from a subject to determine whether the subject has or is at risk for having a disease, comprising: (a) obtaining the cell-free nucleic acid sample comprising a plurality of nucleic acid molecules; (b) subjecting the plurality of nucleic acid molecules or derivatives thereof to sequencing to generate a plurality of sequencing reads; (c) computationally processing the plurality of sequencing reads to identify (i) a methylation profile, (ii) a mutation profile, and (iii) a fragment length profile for the plurality of nucleic acid molecules; and (d) a system comprising one or more computer processors individually or collectively programmed to perform a process comprising using at least the methylation profile, the mutation profile, and the fragment length profile to determine whether the subject has or is at risk of having the disease.

Claims

1. 1. A method of treating a cell-free nucleic acid sample from a subject to test the risk of said subject having a disease, comprising: (a) subjecting a plurality of cell-free nucleic acid molecules from the cell-free nucleic acid sample to sequencing to generate a plurality of sequencing reads; (b) computationally processing the plurality of sequencing reads to identify, for the plurality of cell-free nucleic acid molecules, (i) a methylation profile, (ii) a mutation profile, wherein any mutations present in genomic DNA samples obtained from a plurality of peripheral blood leukocytes obtained from the subject are removed from the mutation profile, and (iii) a fragment length profile, wherein the fragment length profile comprises a length of each cell-free nucleic acid molecule in the plurality of cell-free nucleic acid molecules; and (c) comparing at least said methylation profile, said mutation profile and said fragment length profile with a healthy control to determine whether said subject has or is at risk of having said disease.

2. The method of claim 1 , wherein the disease comprises cancer.

3. The cancers include adrenal gland cancer, anal cancer, bile duct cancer, bladder cancer, bone cancer, brain / CNS tumors, breast cancer, Castleman's disease, cervical cancer, colon / rectal cancer, endometrial cancer, esophageal cancer, Ewing's family of tumors, eye cancer, gallbladder cancer, gastrointestinal carcinoid tumors, gastrointestinal stromal tumors (gist), gestational trophoblastic disease, Hodgkin's disease, Kaposi's sarcoma, kidney cancer, laryngeal and hypopharyngeal cancer, leukemia (acute lymphocytic, acute myeloid, chronic lymphocytic, chronic myeloid, chronic myelomonocytic), liver cancer, lung cancer (non-small cell, small cell, pulmonary carcinoid tumors), lymphoma, lymphoma of the skin, malignant mesothelioma, multiple myeloma, myelodysplastic syndromes ...

3. The method of claim 2, wherein the cancer is selected from the group consisting of adult soft tissue cancer, nasal and paranasal sinus cancer, nasopharyngeal cancer, neuroblastoma, non-Hodgkin's lymphoma, oral cavity and oropharyngeal cancer, osteosarcoma, ovarian cancer, penile cancer, pituitary cancer, prostate cancer, retinoblastoma, rhabdomyosarcoma, salivary gland cancer, sarcoma - adult soft tissue cancer, skin cancer (basal and squamous cell, melanoma, Merkel cell), small intestine cancer, gastric cancer, testicular cancer, thymic cancer, thyroid cancer, uterine sarcoma, vaginal cancer, vulvar cancer, Waldenstrom's macroglobulinemia, Wilms' tumor, squamous cell carcinoma, and squamous cell carcinoma of the head and neck.

4. 4. The method of claim 3, wherein the cancer is squamous cell carcinoma, and the cancer may be head and neck squamous cell carcinoma.

5. The method of any one of claims 1 to 2, wherein the plurality of cell-free nucleic acid molecules comprises circulating tumor nucleic acid molecules.

6. The method of claim 5 , wherein the circulating tumor nucleic acid comprises circulating tumor DNA.

7. The method of claim 5 , wherein the circulating tumor nucleic acid comprises circulating tumor RNA.

8. The method of any one of claims 1 to 7, wherein the methylation profile comprises a plurality of differentially methylated regions (DMRs) derived from the plurality of cell-free nucleic acid molecules.

9. The method of claim 8, wherein the multiple DMRs are derived from ctDNA.

10. 9. The method of claim 8, wherein a plurality of DMRs derived from peripheral blood leukocytes are removed from the methylation profile.

11. 9. The method of claim 8, wherein the plurality of DMRs comprises at least 56 genomic regions that have hypomethylation levels compared to corresponding genomic regions from normal healthy subjects.

12. 2. The method of claim 1, wherein the plurality of DMRs comprises at least 941 genomic regions that have hypermethylation levels compared to corresponding genomic regions from normal healthy subjects.

13. 9. The method of claim 8, wherein one DMR in the plurality of DMRs comprises a size of at least 300 bp, a size of at least 100 bp to at least 200 bp, a size of at least 100 bp to at least 150 bp.

14. The method of claim 8, wherein one DMR in the plurality of DMRs comprises at least 8 CpG genomic islands.

15. The method of claim 11 or 12, wherein the normal, healthy subject comprises the same set of risk factors as the subject.

16. 16. The method of any one of claims 1 to 15, wherein the mutation profile comprises missense variants, nonsense variants, deletion variants, insertion variants, duplication variants, inversion variants, frameshift variants, or repeat expansion variants.

17. The method of any one of claims 1 to 15, wherein any variants derived from clonal hematopoiesis are removed from the mutation profile.

18. 18. The method of claim 17, wherein the mutation profile does not include mutations in genes DNMT3A, TET2, or ASXL1.

19. 18. The method of claim 17, wherein the mutation profile does not include canonical cancer driver genes.

20. 18. The method of claim 17, wherein the mutation profile comprises a non-canonical cancer driver gene, and the non-canonical gene is GRIN3A or MYC.

21. 21. The method of any one of claims 1 to 20, comprising selecting a portion of the plurality of cell-free nucleic acid molecules based on the fragment length profile having a fragment length range of at least 80 bp to 170 bp.

22. 21. The method of any of claims 1 to 20, comprising selecting a portion of the plurality of cell-free nucleic acid molecules based on the fragment length profile having a fragment length range of at least 100 bp to 150 bp.

23. 23. The method of claim 21 or 22, wherein the circulating tumor nucleic acid molecules are enriched.

24. 24. The method of any one of claims 1 to 23, further comprising mixing the cell-free nucleic acid sample with filler DNA molecules to produce a DNA mixture.

25. 25. The method of claim 24, wherein the filler DNA molecule comprises a length of 50 bp to 800 bp.

26. 25. The method of claim 24, wherein the filler DNA molecule comprises a length of 100 bp to 600 bp.

27. 25. The method of claim 24, wherein the filler DNA molecules comprise at least 5%, at least 20%, at least 30%, or at least 50% methylated filler DNA molecules.

28. 28. The method of any one of claims 1 to 27, further comprising incubating the DNA mixture with a binding agent configured to bind to methylated nucleotides to produce an enriched sample.

29. 29. The method of claim 28, wherein the binding agent comprises a protein that includes a methyl-CpG binding domain.

30. 29. The method of claim 28, wherein the protein is an MBD2 protein.

31. 29. The method of claim 28, wherein the binding agent comprises an antibody, which may be a 5-MeC antibody or a 5-hydroxymethylcytosine antibody.

32. 32. The method of any one of claims 1-31, wherein said determining comprises sequencing said plurality of sequencing reads, and wherein said sequencing does not comprise bisulfite sequencing.

33. The method of any one of claims 1 to 31, wherein the cell-free nucleic acid sample comprises a blood sample, and the blood sample may comprise a plasma sample.

34. The method of any one of claims 1 to 33, further comprising detecting the origin of the cancer tissue.

35. The method of any one of claims 1 to 33, further comprising generating a report comprising a prognosis of survival of the subject.

Images