Methods, compositions and systems for calibrating epigenetic compartment assays

Abstract

In one aspect, a method for assessing epigenetic state-based partitions of nucleic acid molecules in a polynucleotide sample is provided, comprising: (a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample to create a spike-in sample; (b) partitioning the nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; (c) enriching a subset of molecules from the plurality of partitioned sets to generate enriched molecules, wherein the enriched molecules comprise a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the polynucleotide sample; (d) sequencing the enriched molecules to generate sequencing reads; and (e) analyzing the sequencing reads to generate one or more epigenetic partition scores for the epigenetic control nucleic acid molecules.

Description

[Technical field]

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS This application claims priority to U.S. Provisional Patent Application No. 62 / 753,826, filed October 31, 2018, which is incorporated by reference in its entirety. [Background technology]

[0002] background Cancer is a leading cause of disease worldwide. Tens of millions of people worldwide are diagnosed with cancer each year, and more than half will ultimately die from the disease. In many countries, cancer ranks second only to cardiovascular disease as the leading cause of death. For many cancers, early detection is associated with improved outcomes.

[0003] Cancer can be caused by the accumulation of genetic mutations in an individual's normal cells, at least some of which result in inappropriate regulation of cell division. Such mutations generally include copy number variations (CNVs), single nucleotide variants (SNVs), gene fusions, insertions and / or deletions (indels), epigenetic mutations including 5-cytosine methylation (5-methylcytosine) and DNA association with chromatin and transcription factors.

[0004] Cancer is often detected by tumor biopsy followed by analysis of cells, markers or DNA extracted from the cells. However, it has been proposed more recently that cancer can also be detected from cell-free nucleic acids in bodily fluids such as blood or urine. Such tests have the advantage of being non-invasive and can be performed without identifying suspected cancer cells by biopsy. However, such liquid biopsy tests are complicated by the fact that the amount of nucleic acid in bodily fluids is very low and that the form in which the nucleic acid is present is heterogeneous (e.g., RNA and DNA, single-stranded and double-stranded, as well as various states of post-replicative modifications and association with proteins such as histones). Summary of the Invention [Means for solving the problem]

[0005] Abstract In one aspect, the disclosure provides a method for assessing epigenetic state-based partitions of nucleic acid molecules in a polynucleotide sample, the method comprising: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample, thereby creating a spike-in sample; b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; c) enriching at least a subset of molecules from the plurality of partitioned sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the polynucleotide sample; d) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; e) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the epigenetic control nucleic acid molecules; and f) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs.

[0006] In another aspect, the disclosure provides a method for assessing epigenetic state-based compartmentalization of nucleic acid molecules in a sample of polynucleotides, comprising the steps of: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the sample of polynucleotides, thereby creating a spike-in sample; b) compartmentalizing at least a subset of nucleic acid molecules of the spike-in sample into a plurality of compartmentalized sets; and c) enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, wherein the set of enriched molecules is a set of epigenetic control nucleic acid molecules and a set of nucleic acids from the sample of polynucleotides. The method includes the steps of: (a) determining whether a set of epigenetic control nucleic acid molecules from a polynucleotide sample comprises a set of endogenous control molecules; (b) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; (c) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of epigenetic control nucleic acid molecules and the endogenous control molecules; and (d) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0007] In another aspect, the disclosure provides a method for assessing epigenetic state-based partitions of nucleic acid molecules in a polynucleotide sample, the method comprising: a) partitioning nucleic acid molecules from at least a subset of the polynucleotide sample into a plurality of partitioned sets; b) enriching at least a subset of molecules from the plurality of partitioned sets to generate a set of enriched molecules, the set of enriched molecules comprising a group of nucleic acid molecules from the polynucleotide sample, the group of nucleic acid molecules from the polynucleotide sample comprising a set of endogenous control molecules; c) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; d) analyzing a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the set of endogenous control molecules; and e) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs.

[0008] In another aspect, the disclosure provides a method for evaluating partitions based on the epigenetic state of nucleic acid molecules in a sample of polynucleotides, the method comprising: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the sample of polynucleotides, thereby generating a spike-in sample; b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; c) sequencing at least a subset of the partition molecules to generate a set of sequencing reads; d) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the epigenetic control nucleic acid molecules; and e) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs. In some embodiments, the method further comprises, prior to the sequencing step, enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the sample of polynucleotides.

[0009] In another aspect, the disclosure provides a method for assessing partitions based on the epigenetic state of nucleic acid molecules in a sample of polynucleotides, the method comprising: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the sample of polynucleotides, thereby generating a spike-in sample; b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; c) sequencing at least a subset of the partition molecules to generate a set of sequencing reads; d) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the set of epigenetic control nucleic acid molecules and endogenous control molecules; and e) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs. In some embodiments, the method further comprises, prior to sequencing, enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules derived from the sample of polynucleotides, wherein the group of nucleic acid molecules derived from the sample of polynucleotides comprises a set of endogenous control molecules.

[0010] In another aspect, the present disclosure provides a method for evaluating epigenetic state-based partitions of nucleic acid molecules in a polynucleotide sample, comprising: a) partitioning molecules from at least one subset of the polynucleotide sample into a plurality of partition sets; b) sequencing at least one subset of the set of enriched molecules to generate a set of sequencing reads; d) analyzing the subset of the set of sequencing reads to generate one or more epigenetic partition scores for a set of endogenous control molecules; and e) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs. In some embodiments, the method further comprises enriching at least one subset of molecules from the plurality of partition sets prior to sequencing to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of nucleic acid molecules from the polynucleotide sample, and wherein the group of nucleic acid molecules from the polynucleotide sample comprises a set of endogenous control molecules.

[0011] In another aspect, the disclosure provides a method for assessing epigenetic state-based partitioning of nucleic acid molecules in a polynucleotide sample, the method comprising: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample, thereby creating a spike-in sample; b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; c) enriching at least a subset of molecules from the plurality of partitioned sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the polynucleotide sample; and d) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads. In some embodiments, the method further comprises: e) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the epigenetic control nucleic acid molecule; and f) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0012] In another aspect, the disclosure provides a method for assessing epigenetic state-based partitioning of nucleic acid molecules in a polynucleotide sample, the method comprising: a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample, thereby creating a spike-in sample; b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; c) enriching at least a subset of molecules from the plurality of partitioned sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules derived from the polynucleotide sample, wherein the group of nucleic acid molecules derived from the polynucleotide sample comprises a set of endogenous control molecules; and d) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads. In some embodiments, the method further includes e) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of epigenetic control nucleic acid molecules and endogenous control molecules; and f) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0013] In another aspect, the present disclosure provides a method for evaluating epigenetic status-based partitions of nucleic acid molecules in a polynucleotide sample, comprising: a) partitioning nucleic acid molecules from at least one subset of the polynucleotide sample into a plurality of partition sets; b) enriching at least one subset of molecules from the plurality of partition sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of nucleic acid molecules from the polynucleotide sample, and the group of nucleic acid molecules from the polynucleotide sample comprises a set of endogenous control molecules; and c) sequencing at least one subset of the set of enriched molecules to generate a set of sequencing reads. In some embodiments, the method further comprises: d) analyzing a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the set of endogenous control molecules; and e) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs.

[0014] In some embodiments, the analyzing step comprises estimating the number / fraction of epigenetic control nucleic acid molecules and / or endogenous control molecules that are in a given epigenetic state in at least one of the compartmentalized sets.

[0015] In some embodiments, the method further comprises tagging the nucleic acid molecules in a compartmentalization set of the plurality of compartmentalization sets with a set of tags to create a population of tagged nucleic acid molecules, the tagged nucleic acid molecules comprising one or more tags. In some embodiments, the set of tags (molecular barcodes) used in a first compartmentalization set of the plurality of compartmentalization sets is different from the set of tags (molecular barcodes) used in a second compartmentalization set of the plurality of compartmentalization sets. In some embodiments, the set of tags is attached to the nucleic acid molecules by ligating an adaptor to the nucleic acid molecules, where the adaptor comprises one or more tags (molecular barcodes). The tag (molecular barcode) sequence used may be correlated with the compartmentalization set, for example, the tag (molecular barcode) used in one compartmentalization set is not used in another compartmentalization set.

[0016] In some embodiments, the method further comprises g) classifying the partition method as (i) a success if each of the one or more epigenetic partition scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules falls within a corresponding epigenetic partition cutoff range, or (ii) a failure if at least one of the one or more epigenetic partition scores of the set of epigenetic control molecules and / or endogenous control molecules is outside a corresponding epigenetic partition cutoff range.

[0017] In some embodiments, the set of epigenetic control nucleic acid molecules comprises two or more subsets of epigenetic control nucleic acid molecules, where one subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules that comprise an epigenetic modification region.

[0018] In some embodiments, the sequencing of the plurality of enriched molecules is performed by a nucleic acid sequencer. In some embodiments, the nucleic acid sequencer is a next generation sequencer.

[0019] In another aspect, the disclosure provides a set of epigenetic control nucleic acid molecules comprising two or more subsets of epigenetic control nucleic acid molecules, wherein a subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules that comprise an epigenetic modification region.

[0020] In another aspect, the disclosure provides a population of nucleic acids comprising: (i) a set of epigenetic control nucleic acid molecules comprising two or more subsets of epigenetic control nucleic acid molecules, wherein a subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules that comprise an epigenetically modified region; and (ii) a set of nucleic acid molecules in a sample of polynucleotides from a subject.

[0021] In some embodiments, the epigenetic control nucleic acid molecule further comprises an identifier region. In some embodiments, the identifier region is present on one or both sides of the epigenetic modification region of the epigenetic control nucleic acid molecule.

[0022] In some embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule in at least one subset comprises at least one nucleotide with an epigenetic modification. In some embodiments, the subsets comprise epigenetic control nucleic acid molecules with the same number of nucleotides with epigenetic modifications. In some embodiments, the number of nucleotides with epigenetic modifications in the first subset is different from the number of nucleotides with epigenetic modifications in the second subset. In some embodiments, the nucleotides with epigenetic modifications comprise methylated nucleotides. In some embodiments, the methylated nucleotides comprise 5-methylcytosine. In some embodiments, the methylated nucleotides comprise 5-hydroxymethylcytosine.

[0023] In some embodiments, the identifier region of the epigenetic control nucleic acid molecule comprises a molecular barcode. In some embodiments, the identifier region further comprises at least one epigenetic state barcode. In some embodiments, the identifier region comprises one or more primer binding sites.

[0024] In some embodiments, the epigenetically modified regions of the epigenetic control nucleic acid molecules in the two or more subsets comprise identical nucleic acid sequences.

[0025] In some embodiments, the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in the first subset comprise nucleic acid sequences that are distinguishable from the nucleic acid sequences of the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in the second subset.

[0026] In some embodiments, the epigenetic modification is DNA methylation.

[0027] In some embodiments, each subset of epigenetic control nucleic acid molecules is present in equimolar concentrations. In some embodiments, each subset of epigenetic control nucleic acid molecules is present in non-equimolar concentrations.

[0028] In some embodiments, the number of methylated nucleotides in the epigenetic control nucleic acid molecules in at least one of the subsets is 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, at least 12, at least 15, at least 20, at least 25, at least 30, at least 40, or at least 50.

[0029] In some embodiments, the epigenetic control nucleic acid molecule comprises a sequence corresponding to lambda phage DNA, a region of the human genome, or a combination of both.

[0030] In some embodiments, the epigenetic state is a methylation level of a nucleic acid molecule. In some embodiments, the plurality of compartmentalized sets comprises nucleic acid molecules of the spike-in sample that are compartmentalized based on the methylation level of the nucleic acid molecule.

[0031] In some embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule comprises a length of about 160 bp.

[0032] In some embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule comprises a nucleic acid sequence that corresponds to a non-human genome.

[0033] In some embodiments, the polynucleotide sample is selected from the group consisting of DNA sample, RNA sample, polynucleotide sample, cell-free DNA sample and cell-free RNA sample.In some embodiments, the polynucleotide sample is selected from the group consisting of DNA sample, RNA sample, polynucleotide sample, cell-free DNA sample and cell-free RNA sample.In some embodiments, the cell-free DNA is between 1ng and 500ng.

[0034] In some embodiments, the epigenetic control nucleic acid molecule is between 1 femtomole and 200 femtomoles.

[0035] In some embodiments, the compartmentalization involves compartmentalizing the nucleic acid molecules based on differential binding affinity of the nucleic acid molecules to binding agents that preferentially bind to nucleic acid molecules that include nucleotides with epigenetic modifications.

[0036] In another aspect, the disclosure provides a system for evaluating a method of partitioning nucleic acid molecules in a sample of polynucleotides based on epigenetic state, comprising a communication interface that receives, over a communication network, a set of sequencing reads of a spike-in sample generated by a nucleic acid sequencer, the set of sequencing reads comprising: (i) at least a first population of sequencing reads generated from polynucleotides originating from the sample, wherein the sequencing reads from the first population comprise tag sequences and sequences derived from the polynucleotides originating from the sample; and (ii) at least a second population of sequencing reads generated from epigenetic control nucleic acid molecules, wherein the sequencing reads generated from the second population comprise an epigenetic modification region and optionally an identifier region. The present invention provides a system comprising: a communications interface including a group; and a computer in communication with the communications interface, the computer including one or more computer processors, and a computer-readable medium including machine-executable code that, when executed by the one or more computer processors, implements a method comprising: receiving, via a communications network, a set of sequencing reads from a first population and a second population of sequencing reads from a nucleic acid sequencer; analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for an epigenetic control nucleic acid molecule and / or an endogenous control molecule; and comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0037] In another aspect, the disclosure provides a system including a controller that includes or can access a computer-readable medium including non-transitory computer-executable instructions that, when executed by at least one electronic processor, perform at least the steps of: (a) obtaining a set of sequencing reads of a spike-in sample generated by a nucleic acid sequencer, the spike-in sample including sample polynucleotides and epigenetic control nucleic acid molecules, the set of sequencing reads including (i) a first population of sequencing reads generated from the sample polynucleotides and (ii) a second population of sequencing reads generated from the epigenetic control nucleic acid molecules; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the epigenetic control nucleic acid molecules and / or endogenous control molecules; and (c) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0038] In another aspect, the disclosure provides a system including a controller that includes or can access a computer-readable medium including non-transitory computer-executable instructions that, when executed by at least one electronic processor, perform at least the steps of: (a) obtaining a set of sequencing reads of a sample generated by a nucleic acid sequencer, wherein the set of sequencing reads includes sequencing reads generated from polynucleotides of the sample; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for an endogenous control molecule; and (c) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0039] In some embodiments, the system further comprises the step of g) generating a compartment method result state based on the comparison of the epigenetic compartment scores. In some embodiments, the compartment method result state is classified as (i) success if one or more epigenetic compartment scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules fall within the range of the corresponding epigenetic compartment cutoff, or (ii) failure if at least one of the one or more epigenetic compartment scores of the epigenetic control molecules and / or endogenous control molecules is outside the range of the corresponding epigenetic compartment cutoff.

[0040] In some embodiments, the epigenetic compartment score comprises a fraction or percentage of the number of hypermethylated epigenetic control nucleic acid molecules and / or hypermethylated control molecules in a compartmentalization set. In some embodiments, the epigenetic compartment score comprises a fraction or percentage of the number of hypomethylated epigenetic control nucleic acid molecules and / or hypomethylated control molecules in a compartmentalization set. In some embodiments, the compartmentalization set is a hypermethylated compartmentalization set. In some embodiments, the compartmentalization set is a hypomethylated compartmentalization set. In some embodiments, the epigenetic compartment score is a 0 CG score. In some embodiments, the epigenetic compartment score is a hypo score. In some embodiments, the epigenetic compartment score is a methyl-half. In some embodiments, the epigenetic compartment score is a methyl-5.

[0041] In some embodiments, the epigenetic compartment cutoff for 0 CG score is 0.01%, 0.02%, 0.05%, 0.1%, 0.2%, 0.3%, 0.4%, 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1%, 2%, 5%, at least 5% or at least 10%. In some embodiments, the epigenetic compartment cutoff for hypo score is 0.1%, 0.5%, 1%, 2%, 3%, 4%, 5%, 7% or at least 10%. In some embodiments, the epigenetic compartment cutoff for methyl-half is 5, 10, 15, 20, 25, 30, 35 or 40 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 is 5, 10, 20, 30, 40 or 50 mCG.

[0042] Additional aspects and advantages of the present disclosure will become readily apparent to those skilled in the art from the following detailed description, in which only illustrative embodiments of the present disclosure are shown and described. As will be understood, the present disclosure is capable of other and different embodiments, and its several details are capable of modification in various obvious respects, all without departing from the present disclosure. Accordingly, the drawings and description are to be regarded as illustrative in nature, and not restrictive.

[0043] In some embodiments, the results of the system and / or method disclosed herein are used as input to generate a report. The report can be in paper or electronic format. For example, such report can present information about the partition of the nucleic acid molecule determined by the method or system disclosed herein and / or information derived from the partition of the nucleic acid molecule. The method or system disclosed herein can further include transmitting the report to the subject from whom the sample is derived or to a third party such as a health care practitioner.

[0044] The various steps of the methods disclosed herein or steps performed by the systems disclosed herein may occur simultaneously or at different times and / or in the same geographic location, e.g., country, or in different geographic locations, e.g., countries. The various steps of the methods disclosed herein may be performed by the same person or by different people.

[0045] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate certain embodiments and, together with the specification, serve to explain certain principles of the methods, computer-readable media, and systems disclosed herein. The description presented herein is better understood when read in conjunction with the accompanying drawings, which are included by way of example and not by way of limitation. It will be understood that like reference numerals identify like components throughout the figures, unless otherwise indicated by the context. It will also be understood that some or all of the drawings may be schematic for illustrative purposes and do not necessarily depict the actual relative size or location of the elements shown. [Brief description of the drawings]

[0046] [Figure 1] 1A and 1B are schematic diagrams of fully methylated (FIG. 1A) and half methylated (FIG. 1B) CpG pairs in double-stranded DNA.

[0047] [Diagram 2] FIG. 2 is a flow chart depicting a method for assessing compartments of a sample of polynucleotides according to certain embodiments of the present disclosure.

[0048] [Diagram 3] FIG. 3 is a flow chart depicting a method for assessing compartments of a sample of polynucleotides according to certain embodiments of the present disclosure.

[0049] [Figure 4]FIG. 4 is a flow chart depicting a method for assessing compartments of a sample of polynucleotides according to certain embodiments of the present disclosure.

[0050] [Diagram 5] FIG. 5 is a schematic diagram of an epigenetic control nucleic acid molecule suitable for use in some embodiments of the present disclosure.

[0051] [Figure 6] FIG. 6 is a schematic diagram of an epigenetic control nucleic acid molecule suitable for use in some embodiments of the present disclosure.

[0052] [Figure 7] FIG. 7 is a schematic diagram of an epigenetic control nucleic acid molecule suitable for use in some embodiments of the present disclosure.

[0053] [Figure 8] FIG. 8 is a schematic diagram of an example system suitable for use in some embodiments of the present disclosure.

[0054] [Figure 9] Figures 9A, 9B and 9C are graphical representations of epigenetic compartment scores of epigenetic control nucleic acid molecules belonging to subsets 1, 2, 3, 4, 5 and 6 within the high compartmentalization set (Figure 9A), the medium compartmentalization set (Figure 9B) and the low compartmentalization set (Figure 9C).

[0055] [Figure 10] 10A and 10B are graphical representations of the fraction of hypermethylated control molecules of Sample 1 within the high compartmentalization set (FIG. 10A) and the low compartmentalization set (FIG. 10B).

[0056] [Figure 11] 11A and 11B are graphical representations of the fraction of hypermethylated control molecules of Sample 2 within the high compartmentalization set (FIG. 11A) and the low compartmentalization set (FIG. 11B).

[0057] definition To facilitate understanding of this disclosure, certain terms are first defined below. Additional definitions for the following terms and other terms may be set forth throughout this specification. In the event that a definition of a term below conflicts with a definition in an application or patent incorporated by reference, the definition set forth in this application should be used to understand the meaning of the term.

[0058] As used in this specification and the appended claims, the singular forms "a," "an," and "the" include plural referents unless the context clearly dictates otherwise. Thus, for example, reference to "a method" includes one or more methods and / or steps of the type described herein and / or that will become apparent to those skilled in the art upon reading this disclosure, and so forth.

[0059] It is also understood that the terminology used herein is only for describing specific embodiments and is not limiting.Furthermore, unless otherwise defined, all scientific and technical terms used herein have the same meaning as commonly understood by those skilled in the art to which this disclosure pertains.In describing and claiming the method, computer readable medium, and system, the following terminology and its grammatical variants are used in accordance with the following definitions.

[0060] About: As used herein, "about" or "approximately" refers to a value or element similar to the stated reference value or element when applied to one or more values ​​or elements of interest. In certain embodiments, the term "about" or "approximately" refers to a range of values ​​or elements that are within 25%, 20%, 19%, 18%, 17%, 16%, 15%, 14%, 13%, 12%, 11%, 10%, 9%, 8%, 7%, 6%, 5%, 4%, 3%, 2%, 1%, or less in either direction (greater or less) of the stated reference value or element, unless otherwise specified or otherwise clear from the context (except where such number exceeds 100% of the possible values ​​or elements).

[0061] Adapter: As used herein, "adapter" refers to a short nucleic acid (e.g., less than about 500 nucleotides, less than about 100 nucleotides, or less than about 50 nucleotides in length), which is generally at least partially double-stranded and attached to either or both ends of a given sample nucleic acid molecule. The adapter may include a nucleic acid primer binding site that allows amplification of the nucleic acid molecule flanked by the adapter at both ends, and / or a sequencing primer binding site, including a primer binding site for sequencing applications, such as various next generation sequencing (NGS) applications. The adapter may also include a binding site for a capture probe, such as an oligonucleotide attached to a flow cell support or the like. The adapter may also include a nucleic acid tag, as described herein. The nucleic acid tag is generally positioned relative to the amplification primer and sequencing primer binding sites such that the nucleic acid tag is included in the amplicon and sequence reads of a given nucleic acid molecule. Adapters of the same or different sequences can be ligated to each end of the nucleic acid molecule. In some embodiments, adapters of the same sequence, except for different nucleic acid tags, are ligated to each end of the nucleic acid molecule. In some embodiments, the adaptor is a Y-shaped adaptor with one end blunt-ended or tailed as described herein for joining with a nucleic acid molecule, where the nucleic acid molecule is also blunt-ended or tailed with one or more complementary nucleotides. In yet other exemplary embodiments, the adaptor is a bell-shaped adaptor that includes a blunt or tailed end for joining with a nucleic acid molecule to be analyzed. Other examples of adaptors include T-tailed adaptors and C-tailed adaptors.

[0062] Amplify: As used herein, "amplify" or "amplification," in reference to nucleic acids, refers to the production of multiple copies of a polynucleotide or a portion of a polynucleotide, generally starting from a small amount of polynucleotide (e.g., a single polynucleotide molecule), where an amplification product or amplicon is generally detectable. Amplification of polynucleotides encompasses a variety of chemical and enzymatic processes. Amplification includes, but is not limited to, polymerase chain reaction (PCR).

[0063] Barcode: As used herein, "barcode" or "molecular barcode," in reference to nucleic acids, refers to a nucleic acid molecule that contains a sequence that can serve as a molecular identifier. For example, individual "barcode" sequences are typically added to DNA fragments during next generation sequencing (NGS) library preparation, thus allowing sequencing reads to be identified and sorted prior to final data analysis.

[0064] Cancer type: As used herein, "cancer type" refers to a type or subtype of cancer defined, for example, by histopathological examination. A cancer type may be classified by any conventional criteria, for example, based on occurrence in a given tissue (e.g., blood cancer, central nervous system (CNS), brain cancer, lung cancer (small cell and non-small cell), skin cancer, nasal cancer, pharyngeal cancer, liver cancer, bone cancer, lymphoma, pancreatic cancer, bowel cancer, rectal cancer, thyroid cancer, bladder cancer, kidney cancer, oral cancer, stomach cancer, breast cancer, prostate cancer, ovarian cancer, lung cancer, intestinal cancer, etc.). Cancers may be defined based on the type of cancer, such as soft tissue cancer, neuroendocrine cancer, gastroesophageal cancer, head and neck cancer, gynecological cancer, colorectal cancer, urothelial cancer, solid state cancer, heterogeneous cancer, homogeneous cancer), unknown primary, etc., and / or based on being of the same cell lineage (e.g., carcinoma, sarcoma, lymphoma, cholangiocarcinoma, leukemia, mesothelioma, melanoma, or glioblastoma), and / or based on cancers exhibiting cancer markers such as, but not limited to, Her2, CA15-3, CA19-9, CA-125, CEA, AFP, PSA, HCG, hormone receptors, and NMP-22. Cancers may also be classified by stage (e.g., stage 1, 2, 3, or 4) and whether primary or secondary.

[0065] Cell-free nucleic acid: As used herein, "cell-free nucleic acid" refers to nucleic acid that is not contained or otherwise bound within a cell, or in some embodiments, nucleic acid that remains in a sample after removal of intact cells. Cell-free nucleic acid may include, for example, all unencapsulated nucleic acid sourced from a subject's body fluid (e.g., blood, plasma, serum, urine, cerebrospinal fluid (CSF), etc.). Cell-free nucleic acid includes DNA (cfDNA), RNA (cfRNA), and hybrids thereof, including genomic DNA, mitochondrial DNA, circulating DNA, siRNA, miRNA, circulating RNA (cRNA), tRNA, rRNA, small nucleolar RNA (snoRNA), piwi-interacting RNA (piRNA), long non-coding RNA strands (long ncRNA), and / or fragments of any of these. Cell-free nucleic acid may be double-stranded, single-stranded, or hybrids thereof. Cell-free nucleic acid may be released into body fluids through secretion or cell death processes, such as cell necrosis, apoptosis, etc. Some cell-free nucleic acids are released from cancer cells into body fluids, for example, circulating tumor DNA (ctDNA). Other cell-free nucleic acids are released from healthy cells. ctDNA can be unencapsulated tumor-derived fragmented DNA. Cell-free nucleic acids can have one or more epigenetic modifications, for example, cell-free nucleic acids can be acetylated, 5-methylated, and / or hydroxymethylated.

[0066] Intracellular Nucleic Acid: As used herein, "intracellular nucleic acid" refers to nucleic acid that is located within the cell or cells from which it is derived, at least at the time the sample is obtained or harvested from a subject, even if it is later removed (e.g., by cell lysis) as part of a given analytical process.

[0067] Coverage: As used herein, the terms "coverage," "total molecule count," or "total allele count" are used interchangeably. These terms refer to the total number of DNA molecules at a particular genomic location in a given sample.

[0068] CpG pair: As used herein, the term "CpG pair" refers to a dinucleotide CpG (cytosine-phosphate-guanine (i.e., a cytosine followed by a guanine in the 5'→3' direction of a nucleic acid sequence)) dinucleotide on the sense strand of a double-stranded DNA molecule and its complementary CpG on the antisense strand (illustrated in FIG. 1).

[0069] Deoxyribonucleic acid or ribonucleic acid: As used herein, "deoxyribonucleic acid" or "DNA" refers to natural or modified nucleotides that have a hydrogen group at the 2' position of the sugar moiety. DNA generally comprises a chain of nucleotides composed of the four types of nucleotide bases; adenine (A), thymine (T), cytosine (C), and guanine (G). As used herein, "ribonucleic acid" or "RNA" refers to natural or modified nucleotides that have a hydroxyl group at the 2' position of the sugar moiety. RNA generally comprises a chain of nucleotides composed of the four types of nucleotide bases; A, uracil (U), G, and C. As used herein, the term "nucleotide" refers to natural or modified nucleotides. Certain pairs of nucleotides specifically bind to each other in a complementary manner (referred to as complementary base pairing). In DNA, adenine (A) pairs with thymine (T) and cytosine (C) pairs with guanine (G). In RNA, adenine (A) pairs with uracil (U) and cytosine (C) pairs with guanine (G). When a first nucleic acid strand binds to a second nucleic acid strand that is composed of nucleotides complementary to those of the first strand, the two strands combine to form a duplex. As used herein, "nucleic acid sequencing data," "nucleic acid sequencing information," "sequence information," "nucleic acid sequence," "nucleotide sequence," "genomic sequence," "gene sequence," or "fragment sequence" or "nucleic acid sequencing read" refers to any information or data that indicates the order and identity of nucleotide bases (e.g., adenine, guanine, cytosine, and thymine or uracil) within a molecule of nucleic acid, such as DNA or RNA (e.g., a whole genome, whole transcriptome, exome, oligonucleotide, polynucleotide, or fragment).It should be understood that the present teachings contemplate sequence information obtained using any of a variety of available techniques, platforms or technologies, including, but not limited to, capillary electrophoresis, microarrays, ligation-based systems, polymerase-based systems, hybridization-based systems, direct or indirect nucleotide discrimination systems, pyrosequencing, ion- or pH-based detection systems, and electronic signature-based systems.

[0070] Endogenous control molecule: As used herein, "endogenous control molecule" refers to a nucleic acid molecule in a sample of polynucleotides that corresponds to at least one human genomic region with an invariant epigenetic state. In some embodiments, the endogenous control molecule can be consistently highly methylated or hypomethylated across tissues, subjects, and cancers. In some embodiments, an endogenous control molecule that corresponds to a human genomic region that has a region that is consistently highly methylated can be referred to as a "hypomethylated control molecule." In some embodiments, an endogenous control molecule that corresponds to a human genomic region that has a region that is consistently hypomethylated can be referred to as a "hypomethylated control molecule."

[0071] Epigenetic control nucleic acid molecule: As used herein, "epigenetic control nucleic acid molecule" refers to a set of nucleic acid molecules that are added to a sample of polynucleotides to evaluate a section of the sample based on epigenetic modification. For example, the epigenetic modification can be DNA methylation, and the epigenetic control nucleic acid molecule can have different / distinguishable levels of methylation. In some embodiments, the epigenetic control nucleic acid molecule comprises an epigenetic modification region and optionally an identifier region. In some embodiments, the epigenetic control nucleic acid molecule comprises an epigenetic modification region and an identifier region. The epigenetic control nucleic acid molecule can be a synthetic oligonucleotide. In some embodiments, the epigenetic control nucleic acid molecule can have a nucleic acid sequence that does not occur in nature. In some embodiments, the epigenetic control nucleic acid molecule can have a nucleic acid sequence that occurs in nature. In some embodiments, the epigenetic control nucleic acid molecule can have a nucleic acid sequence that corresponds to a non-human genome. As non-limiting examples, these molecules can have either (i) sequences corresponding to regions of lambda phage DNA or the human genome, (ii) sequences that do not occur in nature, and / or (iii) a combination of (i) and (ii).

[0072] Epigenetic modification: As used herein, "epigenetic modification" refers to the modification of the base of a nucleotide or nucleotides in a nucleic acid molecule. The modification can be a chemical modification of the base of the nucleotide. In some cases, the modification can be a methylation of the base of the nucleotide. For example, the modification can be a methylation of cytosine resulting in 5-methylcytosine.

[0073] Epigenetic modification region: As used herein, "epigenetic modification region" refers to a region of an epigenetic control nucleic acid molecule that represents the level / degree of epigenetic modification of the epigenetic control nucleic acid molecule. In some embodiments, the epigenetic modification region may include nucleotides with epigenetic modification. In some embodiments, the epigenetic modification is DNA methylation. In these embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule may have methylated nucleotides. The number of methylated nucleotides in the epigenetic modification region may vary among the epigenetic control nucleic acid molecules. In some embodiments, the epigenetic control nucleic acid molecule may have 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, at least 10, at least 15, at least 20, at least 30, at least 40 or at least 50 methylated nucleotides in the epigenetic modification region. Epigenetic control nucleic acid molecules can be grouped into subsets based on the number of nucleotides with epigenetic modification in epigenetic modification region. The length of epigenetic modification region in different subsets can be the same, for example, approximately 160bp. The length of epigenetic modification region between subsets can be different. For example, epigenetic control nucleic acid molecules can be grouped into three subsets (subsets A, B and C) based on the number of methylated nucleotides in epigenetic modification region. Subsets A, B and C can have epigenetic control nucleic acid molecules with 5, 10 and 15 methylated nucleotides in epigenetic modification region, respectively, and the length of epigenetic modification region in subsets A, B and C can be the same (for example, 160bp) or can be different--100bp, 150bp and 200bp for subsets A, B and C, respectively.

[0074] Epigenetic compartment score: As used herein, "epigenetic compartment score" refers to a numerical value that represents the compartment of nucleic acid molecules that belong to a particular epigenetic state in a given compartmentalization set. In some embodiments, the epigenetic compartment score of the nucleic acid molecules that belong to an epigenetic state is determined for each compartmentalization set. For example, the epigenetic compartment score of the epigenetic control nucleic acid molecules and / or endogenous control molecules that belong to a particular epigenetic state can be determined. The epigenetic compartment score can be a measure of the number (or a statistically estimated number) of nucleic acid molecules that belong to a particular epigenetic state. The epigenetic compartment score can be by fraction or percentage. The epigenetic compartment score can be a measure of the ratio of the number of epigenetic control nucleic acid molecules and / or endogenous control molecules that belong to a particular epigenetic state that are compartmentalized in at least one compartmentalization set to the number of epigenetic control nucleic acid molecules and / or endogenous control molecules that belong to that epigenetic state that are present in the other remaining compartmentalization set(s). In some embodiments, the epigenetic compartment score can be a fraction or percentage of the number of epigenetic control nucleic acid molecules and / or endogenous control molecules that belong to a particular epigenetic state that is partitioned in at least one partitioning set relative to the total number of epigenetic control nucleic acid molecules and / or endogenous control molecules that belong to that epigenetic state in all partitioning sets. In some embodiments, the epigenetic compartment score is determined for each epigenetic state of the epigenetic control nucleic acid molecules and / or endogenous control molecules in each of the partitioning sets. In some embodiments, the epigenetic compartment score is determined for the epigenetic control nucleic acid molecules and / or endogenous control molecules that have one or more particular epigenetic states in one or more partitioning sets. In some embodiments, the epigenetic compartment score is determined for the epigenetic control nucleic acid molecules and / or endogenous control molecules that have a particular epigenetic state in a particular partitioning set.

[0075] In some embodiments, the epigenetic compartment score may be directed to the effectiveness with which molecules with no CG ("zero" CG) are partitioned into the highly compartmentalized set. This score may be referred to as a 0 CG score. In some embodiments, the 0 CG score may be expressed as a fraction or percentage of molecules with no CG in the highly compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure for the fraction of epigenetic control nucleic acid molecules and / or the fraction of hypermethylated control molecules in the highly compartmentalized set that have at least one of the following: (i) 1 methyl CG (the epigenetic compartment score can be referred to as 1 CG score), (ii) two methyl CGs (the epigenetic compartment score can be referred to as the 2 CG score); (iii) three methyl CGs (the epigenetic compartment score can be referred to as the 3 CG score); (iv) four methyl CGs (the epigenetic compartment score can be referred to as the 4 CG score), and (v) 5 methyl CGs (the epigenetic compartment score can be referred to as the 5 CG score).

[0076] In some embodiments, the epigenetic compartment score may be directed to the effectiveness of the low methylation control molecules or low methylation epigenetic control nucleic acid molecules being partitioned into the high methylation compartment set. This score may be referred to as the hypo score. In some embodiments, the hypo score may be expressed as a fraction or percentage of the low methylation control molecules or low methylation epigenetic control nucleic acid molecules in the high methylation compartment set. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required for the high methylation control molecules and / or high methylation epigenetic control nucleic acid molecules in the low methylation compartment set to be less than a specified value, for example, 5%. Using the high methylation control molecules and / or high methylation epigenetic control nucleic acid molecules in the low methylation compartment set being 5%. That is, in the example where the epigenetic compartment score is a measure of the number of methylated CGs required for the high methylation control molecules and / or high methylation epigenetic control nucleic acid molecules in the low methylation compartment set to be less than 5%, this score may be conveniently referred to as methyl-5. In some embodiments, the epigenetic compartment score can be a measure of the number of methylated CGs required for the highly methylated control molecules and / or highly methylated epigenetic control nucleic acid molecules in the highly methylated compartment set to be at least a specified value, for example, 50%. Using 50% of the highly methylated control molecules and / or highly methylated epigenetic control nucleic acid molecules in the highly methylated compartment set. That is, in the example where the epigenetic compartment score is a measure of the number of methylated CGs required for at least 50% of the highly methylated control molecules and / or highly methylated epigenetic control nucleic acid molecules in the highly methylated compartment set, this score can be referred to as methyl-half. A wide range of different values ​​from 0% to 100% (rather than just 5% and 50%) can be used in different embodiments, and corresponding different convenient names can be used to refer to the specified values.

[0077] For example, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides. The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalization sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. For each subset, an epigenetic compartment score is determined for each of the compartmentalization sets (P1, P2 and P3). That is, the epigenetic control nucleic acid molecules belonging to subset A have three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. Similarly, each of subsets B and C has three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. An epigenetic compartment score can also be determined for endogenous control molecules.

[0078] In another embodiment, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides (i.e., each subset has a different epigenetic state). The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalized sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. In this embodiment, the epigenetic score is determined only for the subset A molecules in the P1 compartmentalized set. The epigenetic score can be a measure of the fraction or percentage of the subset A molecules in the P1 compartmentalized set relative to the total number of subset A molecules (in the P1, P2 and P3 compartmentalized sets).

[0079] Epigenetic compartment cutoff: As used herein, "epigenetic compartment cutoff" refers to a predetermined cutoff value or cutoff used to assess the compartment of a nucleic acid molecule belonging to a particular epigenetic state in a particular compartmentalization set. In some embodiments, the epigenetic compartment cutoff is determined by analyzing an in-house sample data set. Each compartmentalization set may have an epigenetic compartment cutoff for a nucleic acid molecule belonging to an epigenetic state. If one or more epigenetic compartment scores (used to assess the compartment) of epigenetic control nucleic acid molecules belonging to one or more epigenetic states fall within the range of their corresponding epigenetic compartment cutoff, the compartmentalization method is successful. Otherwise, the compartmentalization method fails. The epigenetic compartment cutoff varies depending on the epigenetic state of the nucleic acid molecule and the compartmentalization set, i.e., each epigenetic state has its own epigenetic compartment cutoff, and any of the compartmentalization sets has a separate epigenetic compartment cutoff for its epigenetic state. The cutoff can be by percentage or fraction, and the cutoff can be a cutoff range rather than a specific cutoff value.For example, the epigenetic compartment cutoff for the epigenetic control nucleic acid molecule belonging to a specific epigenetic state can be between 70% and 79%, between 10% and 15%, and less than 5% for the compartmentalization sets P1, P2 and P3, respectively.If the epigenetic compartment score of the epigenetic control nucleic acid molecule belonging to that epigenetic state falls within the range of the corresponding epigenetic compartment cutoff, the compartmentalization method is successful.

[0080] Epigenetic state: As used herein, "epigenetic state" refers to the level / degree of epigenetic modification of a nucleic acid molecule. For example, when the epigenetic modification is DNA methylation (or hydroxymethylation), the epigenetic state can refer to the presence or absence of methylation at DNA bases (e.g., cytosine) or the degree of methylation within a nucleic acid sequence (e.g., highly methylated, low methylated, moderately methylated or unmethylated nucleic acid molecule). The epigenetic state can also refer to the number of nucleotides that have epigenetic modification. For example, when the epigenetic modification is DNA methylation, the epigenetic state can refer to the number of methylated nucleotides of a nucleic acid molecule.

[0081] Epigenetic state barcode: As used herein, "epigenetic state barcode" refers to a nucleic acid sequence used to identify the epigenetic state of an epigenetic control nucleic acid molecule. Identification can be achieved by having a predetermined correlation between a particular epigenetic state barcode or barcodes and the epigenetic state of the epigenetic control nucleic acid molecule. This can refer to the number of nucleotides with epigenetic modification in the epigenetic modification region of the epigenetic control nucleic acid molecule. In some embodiments, the identifier region of the epigenetic control nucleic acid molecule comprises at least one epigenetic state barcode. For example, if the epigenetic modification is DNA methylation and a subset of the epigenetic control nucleic acid molecules has five methylated nucleotides, then all of the epigenetic control nucleic acid molecules in that subset will have the same epigenetic state barcode. In some embodiments, the epigenetic state barcode can be used to identify the level / degree of epigenetic modification of the epigenetic modification region of the epigenetic control nucleic acid molecule. The epigenetic control nucleic acid molecules can be grouped into subsets based on the number of cytosines or CpG nucleotides in the epigenetic modification region. In some embodiments, the level of methylation can vary within each subset (e.g., highly methylated, moderately methylated, low methylated, or unmethylated), and each level of methylation can have a separate epigenetic state barcode. For example, within subset A, all of the low methylated epigenetic control nucleic acid molecules will have an epigenetic state barcode, e.g., ESB1, and all of the highly methylated epigenetic control nucleic acid molecules will have another epigenetic state barcode, e.g., ESB3. In this example, the epigenetic state barcode is used to identify the level / degree of methylation.

[0082] Human genomic regions with invariant epigenetic status: As used herein, "human genomic regions with invariant epigenetic status" refers to regions in the human genome that have a particular epigenetic status, and the epigenetic status of the region does not vary / change often, but always remains the same or remains consistent even in different subjects and / or different disease types / stages. For example, human genomic regions with invariant epigenetic status may be mostly methylated or mostly unmethylated.

[0083] Identifier region: As used herein, "identifier region" refers to a region of an epigenetic control nucleic acid molecule that is used to distinguish the epigenetic control nucleic acid molecule from other epigenetic control nucleic acid molecules. The identifier region may have a molecular barcode and / or an epigenetic state barcode. The identifier region may be present on one or both sides of the epigenetic modification region. The molecular barcode serves as an identifier of the epigenetic control nucleic acid molecule, while the epigenetic state barcode serves as an identifier of the epigenetic state of the epigenetic control nucleic acid molecule. The identifier region may have an additional region (primer binding site) that facilitates binding of one or more primers.

[0084] Mutant allele count: As used herein, the term "mutant allele count" refers to the number of DNA molecules that carry a mutant allele at a particular genomic locus.

[0085] Mutant allele fraction: As used herein, "mutant allele fraction", "mutation dose" or "MAF" refers to the fraction of nucleic acid molecules that have an allelic alteration or mutation at a given genomic location / locus in a given sample. MAF is generally expressed as a fraction or percentage. For example, the MAF of a somatic mutant may be less than 0.15.

[0086] Mutation: As used herein, "mutation" refers to a mutation from a known reference sequence, including, for example, single nucleotide variants (SNVs) and mutations such as insertions or deletions (indels). Mutations can be germline mutations or somatic mutations. In some embodiments, the reference sequence for comparison is the wild-type genomic sequence of the species of the subject that provides the test sample, typically the human genome.

[0087] Mutation caller: As used herein, "mutation caller" means an algorithm (typically embedded in software or otherwise computer-implemented) used to identify mutations in test sample data (e.g., sequence information obtained from a subject).

[0088] Neoplasm: As used herein, the terms "neoplasm" and "tumor" are used interchangeably. These terms refer to an abnormal growth of cells in a subject. A neoplasm or tumor can be benign, potentially malignant, or malignant. A malignant tumor is referred to as a cancer or cancerous tumor.

[0089] Next generation sequencing: As used herein, "Next generation sequencing" or "NGS" refers to a sequencing technology that has increased throughput compared to traditional Sanger-based and capillary electrophoresis-based techniques, and has the ability to generate hundreds of thousands of relatively small sequence reads at once.Some examples of next generation sequencing techniques include, but are not limited to, sequencing by synthesis, sequencing by ligation, and sequencing by hybridization.In some embodiments, next generation sequencing includes the use of an instrument that can sequence single molecules.

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

[0091] Compartment: As used herein, "compartment" and "epigenetic compartment" are used interchangeably. This term refers to separating or separating nucleic acid molecules based on the characteristics of the nucleic acid molecules (e.g., the level / degree of epigenetic modification). Compartment can be a physical compartment of molecules. Compartment can involve separating nucleic acid molecules into groups or sets based on the level of epigenetic modification (i.e., epigenetic state). For example, nucleic acid molecules can be partitioned based on the level of methylation of the nucleic acid molecules. In some embodiments, the methods and systems used for compartmentation can be found in PCT Patent Application No. PCT / US2017 / 068329, which is incorporated by reference in its entirety.

[0092] Compartmentalized set: As used herein, "compartmentalized set" refers to a set of nucleic acid molecules that are partitioned into sets / groups based on the differential binding affinity of the nucleic acid molecules to a binding agent. The binding agent preferentially binds to nucleic acid molecules that contain nucleotides with epigenetic modifications. For example, if the epigenetic modification is methylation, the binding agent can be a methyl-binding domain (MBD) protein. In some embodiments, the compartmentalized set can include nucleic acid molecules that belong to a particular level / degree of epigenetic modification (i.e., epigenetic state). For example, the nucleic acid molecules can be partitioned into three sets: one set of highly methylated nucleic acid molecules (or highly methylated nucleic acid molecules), which can be referred to as the high methylation compartmentalized set or high compartmentalized set, another set of low methylated nucleic acid molecules (or low methylated nucleic acid molecules), which can be referred to as the low methylation compartmentalized set or low compartmentalized set, and a third set of moderately methylated nucleic acid molecules, which can be referred to as the moderately methylated compartmentalized set or medium compartmentalized set. In another example, nucleic acid molecules can be partitioned based on the number of nucleotides that have an epigenetic modification, where one compartmentalized set can have nucleic acid molecules that have 9 methylated nucleotides and another compartmentalized set can have unmethylated nucleic acid molecules (zero methylated nucleotides).

[0093] Polynucleotide: As used herein, "polynucleotide," "nucleic acid," "nucleic acid molecule," or "oligonucleotide" refers to a linear polymer of nucleosides (including deoxyribonucleosides, ribonucleosides, or analogs thereof) joined by internucleoside linkages. Generally, a polynucleotide contains at least three nucleosides. Oligonucleotides often range in size from a small number, e.g., 3-4 monomeric units, to hundreds of monomeric units. Whenever a polynucleotide is represented by a string such as "ATGCCTG," it will be understood that the nucleotides are in 5'→3' order from left to right, unless otherwise noted, and that in the case of DNA, "A" denotes deoxyadenosine, "C" denotes deoxycytidine, "G" denotes deoxyguanosine, and "T" denotes deoxythymidine. The letters A, C, G, and T may be used to denote the base itself, a nucleoside, or a nucleoside containing a base, as is standard in the art.

[0094] Reference sequence: As used herein, "reference sequence" refers to a known sequence that is used for comparison with experimentally determined sequences. For example, the known sequence can be the entire genome, chromosome, or any segment thereof. The reference generally comprises at least about 20, at least about 50, at least about 100, at least about 200, at least about 250, at least about 300, at least about 350, at least about 400, at least about 450, at least about 500, at least about 1000, or more than 1000 nucleotides. The reference sequence can be aligned with a single continuous sequence of genome or chromosome, or can comprise non-contiguous segments that are aligned with different regions of genome or chromosome. Examples of reference sequences include human genomes, such as hG19 and hG38.

[0095] Sample: As used herein, "sample" means anything that can be analyzed by the methods and / or systems disclosed herein.

[0096] Sequencing: As used herein, "sequencing" refers to any of a number of techniques used to determine the sequence (e.g., identity and order of monomeric units) of a biological molecule, e.g., a nucleic acid such as DNA or RNA. Examples of sequencing methods include, but are not limited to, targeted sequencing, single molecule real-time sequencing, exon or exome sequencing, intron sequencing, electron microscopy-based sequencing, panel sequencing, transistor-mediated sequencing, direct sequencing, random shotgun sequencing, Sanger dideoxytermination sequencing, whole genome sequencing, sequencing by hybridization, pyrosequencing, capillary electrophoresis, gel electrophoresis, double stranded sequencing, cycle sequencing, single base extension sequencing, and the like. These include long sequencing, solid-phase sequencing, high-throughput sequencing, massively parallel processing signature sequencing, emulsion PCR, co-amplification-PCR at low denaturation temperature (COLD-PCR), multiplex PCR, reversible dye terminator sequencing, paired-end sequencing, short-term sequencing, exonuclease sequencing, ligation sequencing, short-read sequencing, single molecule sequencing, synthesis sequencing, real-time sequencing, reverse terminator sequencing, nanopore sequencing, 454 sequencing, Solexa Genome Analyzer sequencing, SOLiD™ sequencing, MS-PET sequencing, and combinations thereof.In some embodiments, sequencing can be carried out by a genetic analyzer, such as the genetic analyzer commercially available from Ilumina, Inc., Pacific Biosciences, Inc., or Applied Biosystems / Thermo Fisher Scientific, among many others.

[0097] Sequence information: As used herein, "sequence information," in reference to a nucleic acid polymer, means the order and identity of the monomeric units (eg, nucleotides) within that polymer.

[0098] Somatic mutation: As used herein, the terms "somatic mutation" or "somatic mutation" are used interchangeably. These terms refer to mutations in the genome that occur after conception. Somatic mutations can occur in any cell of the body other than germ cells, and therefore are not passed on to progeny.

[0099] Spike-in sample: As used herein, a "spike-in sample" is a sample in which an epigenetic control nucleic acid molecule is added to a sample of polynucleotides derived from a subject.

[0100] Subject: As used herein, "subject" refers to an animal, such as a mammalian species (e.g., human) or avian (e.g., bird) species, or other organism, such as a plant. More specifically, the subject may be a vertebrate, e.g., a mammal, such as a mouse, a primate, a monkey, or a human. Animals include farm animals (e.g., beef cattle, dairy cattle, poultry, horses, pigs, etc.), sport animals, and companion animals (e.g., pet or support animals). A subject may be a healthy individual, an individual having or suspected of having a disease or predisposition to a disease, or an individual in need of treatment or suspected of needing treatment. The terms "individual" or "patient" are intended to be interchangeable with "subject."

[0101] For example, the subject may be an individual who has been diagnosed with cancer, an individual who is going to undergo cancer treatment, and / or an individual who has undergone at least one cancer treatment.The subject may be an individual who is in remission from cancer.As another example, the subject may be an individual who has been diagnosed with autoimmune disease.As another example, the subject may be a female individual who is pregnant or planning to become pregnant and has been diagnosed with or suspected of having disease, such as cancer, autoimmune disease. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0102] Detailed Description I. Overview Genomic / epigenetic compartment-based methods can allow simultaneous signal detection for multiple analytes in one assay.However, the detected signals of compartment-based analytes can be poorly resolved and are subject to variable assay conditions that change the sensitivity and specificity of the signals.It is desirable to increase the sensitivity of liquid biopsy assays while reducing the loss of circulating nucleic acid (original material) or data during the process.It is also desirable to provide the ability to compare results across different experiments by controlling assay variability by using one or more controls as described herein.

[0103] The present disclosure provides a method and composition for calibrating epigenetic compartment assay.The present invention includes using a set of epigenetic control nucleic acid molecules with completely resolved genomic / epigenetic features (e.g., distinct number of methylated cytosines in synthetic oligonucleotide duplexes) as control or reference to increase the signal sensitivity and specificity of the sample being analyzed.These molecules can also be used to evaluate the compartment of nucleic acid molecules in a sample based on epigenetic modification and to determine the epigenetic state of nucleic acid molecule(s) in a sample.

[0104] Nucleic acid molecules, such as cell-free polynucleotides, can differ based on epigenetic features, such as methylation. Nucleic acids can have different nucleotide sequences, e.g., specific genes or loci. Features can differ in terms of degree. For example, DNA molecules can differ in the extent of their epigenetic modifications. The degree of modification can refer to the number of modification events to which a molecule is subjected, such as the number of methylation groups (degree of methylation) or other epigenetic changes. For example, methylated DNA can be in a state of hypomethylation or hypermethylation. A form can be characterized by a combination of features, e.g., single-stranded-unmethylated or double-stranded-methylated. Sorting molecules based on one or a combination of features can be useful for multidimensional analysis of single molecules. These methods can be adapted to multiple forms and / or modifications of nucleic acids in a sample, thus obtaining sequence information for multiple forms. Also, the method preserves the identity of the original multiple forms or modification states throughout processing and analysis, thus allowing analysis of nucleic acid base sequences to be combined with epigenetic analysis. Some methods involve separating, tagging, and then pooling different forms or modification states, thereby reducing the number of processing steps required to analyze the multiple forms present in a sample.Analyzing the multiple forms of nucleic acid in a sample provides more information, partly because more molecules are analyzed (which may be important when the total amount of nucleic acid available is very small), but also because different forms or modification states provide different information (e.g., mutations may only exist in RNA), and different types of information (e.g., genetic and epigenetic) can be correlated with each other, thereby generating higher accuracy, certainty, or discovering new correlations with medical conditions.

[0105] A feature of a nucleic acid molecule can be a modification, which can include various chemical modifications (i.e., epigenetic modifications). Non-limiting examples of chemical modifications can include covalent DNA modifications, including, but not limited to, DNA methylation. In some embodiments, DNA methylation includes the addition of a methyl group to a cytosine at a CpG site (a cytosine-phosphate-guanine site (i.e., a cytosine followed by a guanine in the 5'→3' direction of a nucleic acid sequence)). In some embodiments, DNA methylation includes, for example, the addition of a methyl group to a cytosine at a CpG site (a cytosine-phosphate-guanine site (i.e., a cytosine followed by a guanine in the 5'→3' direction of a nucleic acid sequence). 6 In some embodiments, DNA methylation is 5-methylation (modification of the fifth carbon of the six-membered carbon ring of cytosine). In some embodiments, 5-methylation is addition of a methyl group to the 5C position of cytosine to create 5-methylcytosine (m5c). In some embodiments, methylation includes derivatives of m5c. Derivatives of m5c include, but are not limited to, 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC), and 5-carboxylcytosine (5-caC). In some embodiments, DNA methylation is 3C methylation (modification of the third carbon of the six-membered carbon ring of cytosine). In some embodiments, 3C methylation is addition of a methyl group to the 3C position of cytosine to create 3-methylcytosine (3mC). Methylation can also occur at non-CpG sites, for example, methylation can occur at CpA, CpT, or CpC sites. DNA methylation can change the activity of methylated DNA regions. For example, if DNA in a promoter region is methylated, gene transcription can be silenced. DNA methylation is crucial for normal development, and abnormal methylation can disrupt epigenetic regulation. Disruption, e.g., repression, of epigenetic regulation can result in diseases such as cancer. Promoter methylation in DNA can indicate cancer.

[0106] A CpG pair is a dinucleotide CpG (cytosine-phosphate-guanine, i.e., a cytosine followed by a guanine in the 5'→3' direction of a nucleic acid sequence) on the sense strand of a double-stranded DNA molecule and its complementary CpG on the antisense strand. A CpG pair is either fully methylated or half methylated. FIG. 1 is a schematic diagram of fully methylated and half methylated CpG pairs in double-stranded DNA. FIG. 1A shows a fully methylated CpG pair 103, where the cytosine nucleotides of both CpG pairs, strands 101 and 102, are methylated (M-methylcytosine; G-guanine). FIG. 1B shows a half methylated CpG pair 104, where the cytosine nucleotides of the CpG pair on one strand 101 are methylated but the cytosine nucleotides of the CpG pair on the complementary strand 102 are unmethylated (C-unmethylated cytosine; G-guanine).

[0107] CpG dinucleotides are under-abundant in the normal human genome, and the majority of CpG dinucleotide sequences are transcriptionally inactive (e.g., DNA heterochromatic regions within pericentromeric portions of chromosomes and within repeat elements) and are methylated, but many CpG islands, especially around transcription start sites (TSSs), are protected from such methylation.

[0108] Cancer can be manifested by epigenetic alterations such as methylation. Examples of methylation changes in cancer include localized increases in DNA methylation in CpG islands at the TSS of genes involved in normal growth control, DNA repair, cell cycle regulation, and / or cell differentiation. This hypermethylated state may be associated with abnormal loss of transcriptional capacity of the involved genes, and occurs at least as frequently as point mutations and deletions as a cause of altered gene expression. DNA methylation profiling can be used to detect regions of the genome with different degrees of methylation ("differentially methylated regions" or "DMRs") that are altered during development or perturbed by disease, such as cancer or any cancer-related disease.

[0109] Methylation profiling can involve determining methylation patterns across different regions of a genome. For example, molecules can be partitioned based on the degree of methylation (e.g., the relative number of methylated nucleotides per molecule) and sequenced, and then the sequences of the molecules in the different partitions can be mapped to a reference genome. This can indicate regions of the genome that are more or less highly methylated compared to other regions. In this way, genomic regions can differ in their degree of methylation as opposed to individual molecules. In addition to methylation, other epigenetic modifications can be profiled similarly.

[0110] Nucleic acid molecules in a sample can be sorted or partitioned based on one or more features. Partitioning nucleic acid molecules in a sample can increase rare signals. For example, by partitioning a sample into hypermethylated and hypomethylated nucleic acid molecules, genetic mutations that are present in hypermethylated DNA but less present (or absent) in hypomethylated DNA can be more easily detected. By analyzing multiple fractions of a sample, multidimensional analysis of single molecules can be performed, thus achieving higher sensitivity. Partitioning can include physically partitioning nucleic acid molecules into subsets or groups based on the presence or absence of a genomic feature. Sorting can include physically partitioning nucleic acid molecules into groups based on the degree to which a genomic feature, such as an epigenetic modification, is present. A sample can be sorted or partitioned into one or more group groups based on features indicative of differential gene expression or disease state. Samples can be separated based on a feature, or a combination thereof, that results in a difference in signal between normal and disease states during analysis of nucleic acids, such as cell-free DNA ("cfDNA"), non-cfDNA, tumor DNA, circulating tumor DNA ("ctDNA"), and cell-free nucleic acid ("cfNA").

[0111] The present disclosure provides methods, compositions and systems for assessing or evaluating compartments of nucleic acid molecules and determining the epigenetic state (e.g., methylation state) and the number of epigenetically modified nucleotides (e.g., number of methylated nucleotides) in the nucleic acid molecules. The methods may include partitioning the nucleic acid molecules into different compartmentalized sets based on one or more epigenetic modifications, and then sequencing (single or together) to analyze the nucleic acid molecules in each compartment. In some embodiments, the compartments of nucleic acids are enriched for specific target genomic regions. In some embodiments, the compartments of nucleic acid molecules are amplified before and / or after enrichment. In some embodiments, the compartmentalized sets can be differentially tagged with molecular barcodes and recombined in a mixture of differentially tagged compartmentalized sets before enrichment. The methods can be used for various applications, such as prognosing, diagnosing and / or monitoring disease. In some embodiments, the disease is cancer.

[0112] The partitioning method of nucleic acid molecules can be evaluated by using epigenetic control nucleic acid molecules. Epigenetic control nucleic acid molecules are synthetic nucleic acid molecules that can have epigenetically modified nucleotides. In some embodiments, epigenetic control nucleic acid molecules can include nucleic acid molecules with different epigenetic states. Epigenetic state can refer to the level / degree of epigenetic modification of a nucleic acid molecule. For example, when the epigenetic modification is DNA methylation, the epigenetic state can refer to highly methylated, low methylated or moderately methylated nucleic acid molecules. Epigenetic state can also refer to the number of nucleotides that have epigenetic modification. For example, when the epigenetic modification is DNA methylation, the epigenetic state can refer to the number of methylated nucleotides of a nucleic acid molecule. Epigenetic modification can be any modification of the base of a nucleotide(s) that does not involve a change in the sequence and / or base pairing specificity of a nucleic acid molecule. The modification can be a chemical modification of the base of a nucleotide. In some cases, the modification can be a methylation of the base of a nucleotide. For example, the modification can be the methylation of cytosine resulting in 5-methylcytosine.

[0113] In some embodiments, the epigenetic control nucleic acid molecule is a synthetic molecule, and the sequence of the epigenetic control nucleic acid molecule and the position and number of epigenetically modified nucleotides in the epigenetic control nucleic acid molecule are already known prior to analysis. Thus, by adding the epigenetic control nucleic acid molecule to a sample of polynucleotides and tracking the epigenetic control nucleic acid molecule in the compartmentalized set, the effect of the compartment of the epigenetic control nucleic acid molecule can be analyzed.

[0114] Thus, in one aspect, the disclosure provides a method for assessing epigenetic state-based partitions of nucleic acid molecules in a polynucleotide sample, the method comprising: (a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample, thereby creating a spike-in sample; (b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of partitioned sets; (c) enriching a subset of molecules from the plurality of partitioned sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the polynucleotide sample; (d) sequencing the set of enriched molecules to generate a set of sequencing reads; (e) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic partition scores for the epigenetic control nucleic acid molecules; and (f) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs. In these embodiments, the partitioning of the sample nucleic acid molecules and the epigenetic control nucleic acid molecules is necessarily performed simultaneously. In some embodiments, the analyzing step comprises estimating the number / fraction of epigenetic control nucleic acid molecules that are in a given epigenetic state in at least one of the partitioned sets.

[0115] 2 illustrates an example embodiment of a method 200 for assessing epigenetic state-based compartments of nucleic acid molecules in a sample of polynucleotides. In 201, epigenetic control nucleic acid molecules whose compartments are to be assessed are added to the sample to generate a spike-in sample.

[0116] In some embodiments, the epigenetic control nucleic acid molecule may include one or more subsets of nucleic acid molecules that have different levels of epigenetic state (i.e., different numbers of epigenetically modified nucleotides). In some embodiments, the epigenetic control nucleic acid molecule may include nucleic acid molecules that differ in sequence and / or length. In other embodiments, the epigenetic control nucleic acid molecule may include nucleic acid molecules that are identical in sequence or identical in length.

[0117] At 202, at least a subset of nucleic acid molecules of the spike-in sample, including both epigenetic control nucleic acid molecules and nucleic acid molecules from the sample of polynucleotides, are partitioned or sorted into a plurality of partitioned sets based on the epigenetic state of the molecules. The partitions can be based on the presence or absence of epigenetic modifications and / or based on the degree of epigenetic modifications. Examples of epigenetic modifications include, but are not limited to, the presence or absence of methylation, the level of methylation, and the type of methylation (5' cytosine). In some embodiments, the epigenetic modification can be DNA methylation. In these embodiments, the molecules of the spike-in sample are partitioned based on different levels of methylation (number of differently methylated nucleotides). In some embodiments, the spike-in sample can be partitioned into two or more partitioned sets (e.g., at least three, four, five, six, or seven partitioned sets). In some embodiments, the partitions are based on the differential binding affinity of the nucleic acid molecules to the binding agent. Examples of binding agents include, but are not limited to, methyl-binding domains (MBDs) and methyl-binding proteins (MBPs). Examples of MBPs contemplated herein include, but are not limited to, the following: (a) MeCP2 is a protein that binds preferentially to 5-methyl-cytosine over unmodified cytosine; (b) RPL26, PRP8 and the DNA mismatch repair protein MHS6 bind 5-hydroxymethyl-cytosine preferentially over unmodified cytosine; (c) FOXK1, FOXK2, FOXP1, FOXP4 and FOXI3 preferentially bind 5-formyl-cytosine over unmodified cytosine (Iurlaro et al., Genome Biol. 14, R119 (2013)); and (d) an antibody specific for one or more methylated nucleotide bases.

[0118] For some affinity agents and modifications, binding to the agent occurs essentially in an all-or-none manner depending on whether the nucleic acid has the modification, but there may be various degrees of separation. In such an embodiment, nucleic acids with an overabundance of the modification bind to the agent to a greater extent than nucleic acids with an underabundance of the modification. Alternatively, nucleic acids with modifications may bind in an all-or-none manner. However, different levels of modifications can then be sequentially eluted from the binding agent.

[0119] For example, in some embodiments, the partitions can be binary or based on the degree / level of modification. For example, a methyl-binding domain protein (e.g., MethylMiner Methylated DNA Enrichment Kit (ThermoFisher Scientific)) can be used to partition all methylated fragments from unmethylated fragments. Additional partitioning can then involve eluting fragments with different levels of methylation by adjusting the salt concentration of the solution with the methyl-binding domain and bound fragments. As the salt concentration increases, fragments with higher methylation levels elute.

[0120] In some embodiments, the compartmentalization involves compartmentalizing the nucleic acid molecules based on differential binding affinity of the nucleic acid molecules to binding agents that preferentially bind to nucleic acid molecules that include nucleotides with epigenetic modifications.

[0121] In some embodiments, the compartmentalized sets represent nucleic acids with different degrees of modification (over- or under-representation of the modification). Over- and under-representation can be defined by the number of modifications that a nucleic acid has compared to the median number of modifications per strand in the population. For example, if the median number of 5-methylcytosine nucleotides in nucleic acid molecules in a sample is 2, then nucleic acid molecules that contain more than two 5-methylcytosine residues are over-represented in this modification, and nucleic acids with one or zero 5-methylcytosine residues are under-represented. The effect of affinity separation is to partition for nucleic acids in the binding phase, where the modification is over-represented, and for nucleic acids in the non-binding phase (i.e., in solution), where the modification is under-represented. Nucleic acids in the binding phase can be eluted before further processing.

[0122] Using the MethylMiner Methylated DNA Enrichment Kit (ThermoFisher Scientific), different levels of methylation can be partitioned using sequential elution. For example, the low methylated compartment (no methylation) can be separated from the methylated compartment by contacting the nucleic acid population with MBD from the kit attached to magnetic beads. The beads are used to separate the methylated nucleic acids from the unmethylated nucleic acids. Then, one or more elution steps are performed sequentially to elute the nucleic acids with different levels of methylation. For example, the first set of methylated nucleic acids can be eluted at a salt concentration of about 150 mM or about 160 mM or higher, for example, at least 150 mM, 200 mM, 300 mM, 400 mM, 500 mM, 600 mM, 700 mM, 800 mM, 900 mM, 1000 mM, or 2000 mM. After elution of such methylated nucleic acids, magnetic separation is once again used to separate the more highly methylated nucleic acids from the less highly methylated nucleic acids. The elution and magnetic separation steps themselves can be repeated to create various compartments, such as hypomethylated compartments (representing no methylation), methylated compartments (representing low levels of methylation), and hypermethylated compartments (representing high levels of methylation).

[0123] In some methods, the nucleic acids bound to the agent used for affinity separation are subjected to a washing step, in which nucleic acids that are weakly bound to the affinity agent are washed away. Such nucleic acids can be enriched for nucleic acids that have modifications that are somewhat close to the average or median (i.e., halfway between the nucleic acids that remain bound to the solid phase and those that are not bound to the solid phase when the sample is first contacted with the agent). The result of affinity separation is at least two, and sometimes three or more compartments of nucleic acids with different degrees of modification.

[0124] The compartments of nucleic acid molecules can be analyzed by sequencing the compartmented nucleic acid molecules, or by digital droplet PCR (ddPCR), or by quantitative PCR (qPCR).Before analyzing the compartments, the nucleic acid molecules in the compartmentalized set can be enriched, so that the signal from the nucleic acid molecule of interest can be increased, and thus the sensitivity can be improved.In 203, at least a subset of the nucleic acid molecules in the multiple compartmentalized sets is enriched, and thus the epigenetic control nucleic acid molecules and the nucleic acid molecules from the sample of polynucleotides that belong to the region of interest are enriched.

[0125] In some embodiments, each of the multiple compartmentalization sets is differentially tagged prior to enrichment. The tagged compartmentalization sets are then pooled together for collective sample preparation and / or sequencing. Differential tagging of the compartmentalization sets helps to keep track of the nucleic acid molecules that belong to a particular compartmentalization set. The tags are usually provided as components of the adapter. Nucleic acid molecules in different compartmentalization sets receive different tags that can distinguish members of one compartmentalization set from members of another compartmentalization set. The tags linked to nucleic acid molecules of the same compartment set can be the same or different from each other. However, if different from each other, the tags can have a portion of their sequence in common so that the molecules to which they are attached are identified as belonging to a particular compartmentalization set. For example, if the molecules of the spike-in sample are partitioned into two compartmentalization sets, P1 and P2, the molecules in P1 can be tagged with A1, A2, A3, etc., and the molecules in P2 can be tagged with B1, B2, B3, etc. Such a tagging system allows for differentiation between compartmentalization sets and between molecules within a compartmentalization set.

[0126] At least one subset of the enriched molecules is sequenced in 204. The sequence information obtained includes the sequence of the nucleic acid molecule and the tag attached to the nucleic acid molecule. From the sequence of the tag attached to the nucleic acid molecule, the tag can be correlated with a compartmentalized set of nucleic acid molecules. The sequence information is used to identify epigenetic control nucleic acid molecules and their corresponding compartmentalized sets. This information is used to analyze the compartments of the epigenetic control nucleic acid molecules. In 205, one or more epigenetic compartment scores of the epigenetic control nucleic acid molecules belonging to one or more epigenetic states in one or more compartmentalized sets are determined. In some embodiments, the sensitivity and / or specificity of the compartmentalization method can be assessed by the epigenetic compartment score. The epigenetic compartment score is a score that represents the compartment of the nucleic acid molecule belonging to a particular epigenetic state. The epigenetic compartment score of the nucleic acid molecule belonging to the epigenetic state is determined for each compartmentalized set. For example, the epigenetic compartment score of the epigenetic control nucleic acid molecule belonging to a particular epigenetic state can be determined. The epigenetic compartment score can be a measure of the number (or a statistically estimated number) of nucleic acid molecules belonging to a particular epigenetic state. The epigenetic compartment score can be a fraction or percentage. The epigenetic compartment score can be a measure of the ratio of the number of epigenetic control nucleic acid molecules belonging to a particular epigenetic state that are partitioned in at least one partitioning set to the number of epigenetic control nucleic acid molecules belonging to that epigenetic state that are present in the other remaining partitioning set(s). In some embodiments, the epigenetic compartment score can be a fraction or percentage of the number of epigenetic control nucleic acid molecules belonging to a particular epigenetic state that are partitioned in at least one partitioning set to the total number of epigenetic control nucleic acid molecules belonging to that epigenetic state in all partitioning sets. In some embodiments, the epigenetic compartment score is determined for each epigenetic state of epigenetic control nucleic acid molecules in each of the partitioning sets.In some embodiments, the epigenetic compartment score is determined for epigenetic control nucleic acid molecules having one or more particular epigenetic states in one or more compartmentalized sets. In some embodiments, the epigenetic compartment score is determined for epigenetic control nucleic acid molecules having a particular epigenetic state in a particular compartmentalized set.

[0127] In some embodiments, the epigenetic compartment score may address the effectiveness with which molecules with no CG ("zero" CG) are partitioned into the highly compartmentalized set. This score may be referred to as a 0 CG score. In some embodiments, the 0 CG score may be expressed as a fraction or percentage of molecules with no CG in the highly compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the fraction of epigenetic control nucleic acid molecules and / or the fraction of hypermethylated control molecules that have at least one of the following: In the hypermethylated compartmentalization set (i.e., the highly methylated compartmentalization set), (i) 1 methyl CG (the epigenetic compartment score can be referred to as 1 CG score), (ii) two methyl CGs (the epigenetic compartment score can be referred to as the 2 CG score); (iii) three methyl CGs (the epigenetic compartment score can be referred to as the 3 CG score); (iv) 4 methyl CGs (the epigenetic compartment score can be referred to as the 4 CG score); and (v) 5 methyl CGs (the epigenetic compartment score can be referred to as the 5 CG score).

[0128] In some embodiments, the epigenetic compartment score may refer to the effectiveness of the hypomethylated (i.e., low methylated) epigenetic control nucleic acid molecules being partitioned into the hypermethylated compartmentalized set. This score may be referred to as the hypo score. In some embodiments, the hypo score may be expressed as a fraction or percentage of the hypomethylated epigenetic control nucleic acid molecules in the hypermethylated compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to achieve less than 5% hypermethylated epigenetic control nucleic acid molecules in the hypomethylated compartmentalized set. This score may be referred to as methyl-5. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to achieve at least 50% hypermethylated epigenetic control nucleic acid molecules in the hypermethylated compartmentalized set. This score may be referred to as methyl-half.

[0129] For example, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides. The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalization sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. For each subset, an epigenetic compartment score is determined for each of the compartmentalization sets (P1, P2 and P3). That is, the epigenetic control nucleic acid molecules belonging to subset A have three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. Similarly, each of subsets B and C has three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. An epigenetic compartment score can also be determined for endogenous control molecules.

[0130] In another embodiment, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides (i.e., each subset has a different epigenetic state). The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalized sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. In this embodiment, the epigenetic score is determined only for the subset A molecules in the P1 compartmentalized set. The epigenetic score can be a fraction or percentage measure of the subset A molecules in the P1 compartmentalized set relative to the total number of subset A molecules (in the P1, P2 and P3 compartmentalized sets).

[0131] The epigenetic compartment score can be any value or range between 0 and 1 (in fractions) or between 0 and 100% (in percentages). In some embodiments, the epigenetic compartment score can be by the number of methylated CGs (e.g., for methyl-half and methyl-5).

[0132] In 206, the epigenetic compartment score of the epigenetic control nucleic acid molecule is compared to an epigenetic compartment cutoff (predetermined cutoff) to evaluate the partitioning method. The epigenetic compartment cutoff is a predetermined cutoff value or cutoff range used to evaluate the partition of the nucleic acid molecule belonging to a specific epigenetic state, and each partitioning set has an epigenetic compartment cutoff for the nucleic acid molecule belonging to the epigenetic state. The epigenetic compartment cutoff varies depending on the epigenetic state of the nucleic acid molecule and the partitioning set, i.e., each epigenetic state has its own epigenetic compartment cutoff, and any of the partitioning sets has a separate epigenetic compartment cutoff for its epigenetic state. The cutoff can be by percentage or fraction, and the cutoff can be a cutoff range rather than a specific cutoff value. For example, the epigenetic compartment cutoff for epigenetic control nucleic acid molecules belonging to a particular epigenetic state can be between 70% and 79%, between 10% and 15%, and less than 5% for the compartmentalization sets P1, P2, and P3, respectively. If the epigenetic compartment score of the epigenetic control nucleic acid molecule belonging to that epigenetic state falls within the range of the corresponding epigenetic compartment cutoff, the compartmentalization method is successful. In some embodiments, the epigenetic compartment cutoff for a 0 CG score can be 0.01%, 0.02%, 0.05%, 0.1%, 0.2%, 0.3%, 0.4%, 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1%, 2%, 5%, at least 5% or at least 10%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score can be 0.01%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.02%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.03%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.04%.In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.05%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.2%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.3%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.4%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.6%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.7%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.8%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.9%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 1%.

[0133] In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%, 0.5%, 1%, 2%, 3%, 4%, 5%, 7%, or at least 10%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 2%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 3%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 4%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 5%.

[0134] In some embodiments, the cutoff may be by the number of methylated CGs (e.g., for methyl-5 and methyl-half). In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5, 10, 20, 30, 40 or 50 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 40 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 50 mCG.

[0135] In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5, 10, 15, 20, 25, 30, 35 or 40 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 15 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 25 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 35 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 40 mCG.

[0136] In some embodiments, the partitioning method can be successfully classified if one or more epigenetic partition scores of epigenetic control nucleic acid molecules belonging to one or more epigenetic states within one or more partitioning sets fall within the range of the corresponding epigenetic partition cutoffs. Otherwise, if the epigenetic partition scores are outside the cutoff for all partitioning sets, the partitioning method can be classified as a failure. For example, there are two subsets of epigenetic control nucleic acid molecules, subset A and B, and each subset has a different degree of epigenetic modification (i.e., each subset has a different epigenetic state). These epigenetic control nucleic acid molecules can be partitioned into two partitioning sets, P1 and P2. For the molecules belonging to subset A, one or two epigenetic partition scores (e.g., S1 and S2) are determined for each of the partitioning sets P1 and P2 based on their partitions. Similarly, for the molecules belonging to subset B, one epigenetic partition score for P1 and one for P2, two epigenetic partition scores (e.g., S3 and S4) are determined. Each subset of molecules having a particular epigenetic state has a predetermined epigenetic partition cutoff for each of the partitioning sets. In this example, the epigenetic control nucleic acid molecules of subset A have two epigenetic partition cutoffs, C1 and C2 (for the two partitioning sets P1 and P2), and similarly, the epigenetic control nucleic acid molecules of subset B have two epigenetic partition cutoffs, C3 and C4. The epigenetic partition scores of both subsets are compared with their corresponding epigenetic partition cutoffs. In this example, the partitioning method is considered successful only if all four epigenetic partition scores fall within the range of their corresponding epigenetic partition cutoffs, i.e., in this example, if S1 < C1 and S2 < C2 and S3 < C3 and S4 < C4. Otherwise, if the epigenetic partition scores are outside the cutoff for all partitioning sets, the partitioning method can be classified as a failure.

[0137] In another embodiment, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides (i.e., each subset has a different epigenetic state). The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalized sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. In this embodiment, the epigenetic score is determined only for the subset A molecules in the P1 compartmentalized set. This epigenetic score can be a measure of the fraction or percentage of the subset A molecules in the P1 compartmentalized set relative to the total number of subset A molecules (in the P1, P2 and P3 compartmentalized sets). If this epigenetic compartment score falls within its corresponding epigenetic compartment cutoff range, the compartmentalization method is classified as successful. Otherwise, the compartmentalization method is classified as unsuccessful.

[0138] In another aspect, the disclosure provides a method for assessing epigenetic state-based compartmentalization of nucleic acid molecules in a polynucleotide sample, comprising the steps of: (a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the polynucleotide sample, thereby creating a spike-in sample; (b) compartmentalizing at least a subset of nucleic acid molecules of the spike-in sample into a plurality of compartmentalized sets; and (c) enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, wherein the set of enriched molecules is a set of epigenetic control nucleic acid molecules and a set of nucleic acids from the polynucleotide sample. The method includes the steps of: (d) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; (e) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of epigenetic control nucleic acid molecules and endogenous control molecules; and (f) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs. In these embodiments, the partitioning of the nucleic acid molecules of the sample and the epigenetic control nucleic acid molecules is necessarily performed simultaneously. In some embodiments, the analyzing step includes estimating the number / fraction of epigenetic control nucleic acid molecules and / or endogenous control molecules in a given epigenetic state in at least one of the partitioned sets.

[0139] FIG. 3 illustrates an example embodiment of a method 300 for evaluating a partition based on the epigenetic state of a nucleic acid molecule in a sample of polynucleotides. In this embodiment, both the partition of the epigenetic control nucleic acid molecule and the endogenous control molecule are analyzed to evaluate the partition method. There are regions in the human genome that have a specific epigenetic state, and the epigenetic state of the region often does not vary / change and always remains the same / remains consistent even in different subjects and / or different disease types / stages. The nucleic acid molecule in the sample of polynucleotides that corresponds to such a human genome region with an invariant epigenetic state is referred to as an endogenous control molecule. In 301, the epigenetic control nucleic acid molecule is added to the sample of polynucleotides whose partition is to be evaluated to generate a spike-in sample.

[0140] In some embodiments, the epigenetic control nucleic acid molecule may include one or more subsets of nucleic acid molecules that have different levels of epigenetic state (i.e., different numbers of epigenetically modified nucleotides). In some embodiments, the epigenetic control nucleic acid molecule may include nucleic acid molecules that differ in sequence and / or length. In other embodiments, the epigenetic control nucleic acid molecule may include nucleic acid molecules that are identical in sequence or identical in length.

[0141] At 302, at least a subset of nucleic acid molecules of the spike-in sample, including both epigenetic control nucleic acid molecules and nucleic acid molecules from the sample of polynucleotides, are partitioned or sorted into a plurality of partition sets based on the epigenetic state of the molecules. The partitions can be based on the presence or absence of epigenetic modifications and / or based on the degree of epigenetic modifications. Examples of epigenetic modifications include, but are not limited to, the presence or absence of methylation, the level of methylation, and the type of methylation (5' cytosine). In some embodiments, the epigenetic modification can be DNA methylation. In these embodiments, the molecules of the spike-in sample are partitioned based on different levels of methylation (number of differently methylated nucleotides). In some embodiments, the spike-in sample can be partitioned into two or more partition sets (e.g., at least 3, 4, 5, 6, or 7 partition sets). In some embodiments, the partitions are based on the differential binding affinity of the nucleic acid molecules to the binding agent.

[0142] The compartments of nucleic acid molecules can be analyzed by sequencing the compartmented nucleic acid molecules, by digital droplet PCR (ddPCR), or by quantitative PCR (qPCR). Before analyzing the compartments, the nucleic acid molecules in the compartmentalized set can be enriched, thereby increasing the signal from the nucleic acid molecule of interest and thus improving the sensitivity. In 303, at least a subset of the nucleic acid molecules in the multiple compartmentalized sets is enriched, thus enriching the epigenetic control nucleic acid molecules belonging to the region of interest, the endogenous control molecules (from the sample of polynucleotides) and the nucleic acid molecules from other samples of polynucleotides.

[0143] In some embodiments, each of the multiple compartmentalization sets is differentially tagged before enrichment. The tagged compartmentalization sets are then pooled together for collective sample preparation and / or sequencing. Differential tagging of the compartmentalization sets helps to keep track of the nucleic acid molecules that belong to a particular compartmentalization set. The tags are usually provided as components of the adapter. The nucleic acid molecules in different compartmentalization sets receive different tags that can distinguish the members of one compartmentalization set from the members of another compartment cassette. The tags linked to the nucleic acid molecules of the same compartment set can be the same or different from each other. However, when different from each other, the tags can have a part of their sequence in common so that the molecules to which they are attached can be identified as belonging to a particular compartmentalization set.

[0144] At 304, at least a subset of the enriched molecules is sequenced. The sequence information obtained includes the sequence of the nucleic acid molecule and the tag attached to the nucleic acid molecule. From the sequence of the tag attached to the nucleic acid molecule, the tag can be correlated with a compartmentalized set of nucleic acid molecules. The sequence information is used to identify epigenetic control nucleic acid molecules and endogenous control molecules and their corresponding compartmentalized sets. This information is used to analyze the compartments of the epigenetic control nucleic acid molecules and endogenous control molecules. At 305, one or more epigenetic compartment scores of the epigenetic control nucleic acid molecules and endogenous control molecules that belong to one or more epigenetic states in one or more compartmentalized sets are determined. In some embodiments, the sensitivity and / or specificity of the compartmentalization method can be assessed by the epigenetic compartment score. The epigenetic compartment score is a score that represents the compartment of the nucleic acid molecule that belongs to a particular epigenetic state. The epigenetic compartment score of the nucleic acid molecule that belongs to the epigenetic state is determined for each compartmentalized set. For example, the epigenetic compartment score of epigenetic control nucleic acid molecules and endogenous control molecules belonging to a particular epigenetic state can be determined. The epigenetic compartment score can be a measure of the number (or a statistically estimated number) of nucleic acid molecules belonging to a particular epigenetic state. The epigenetic compartment score can be a fraction or percentage. The epigenetic compartment score can be (i) for epigenetic control nucleic acid molecules, a measure of the ratio of the number of epigenetic control nucleic acid molecules belonging to a particular epigenetic state that are partitioned into at least one partitioning set to the number of epigenetic control nucleic acid molecules belonging to that epigenetic state that are present in the other remaining partitioning set(s), and (ii) for endogenous control molecules, a measure of the ratio of the number of endogenous control molecules belonging to a particular epigenetic state that are partitioned into at least one partitioning set to the number of endogenous control molecules belonging to that epigenetic state that are present in the other remaining partitioning set(s).In some embodiments, the epigenetic compartment score can be (i) for epigenetic control nucleic acid molecules, the number of epigenetic control nucleic acid molecules that belong to a particular epigenetic state that is partitioned in at least one partitioning set, relative to the total number of epigenetic control nucleic acid molecules that belong to that epigenetic state in all partitioning sets, or (ii) for endogenous control molecules, the number of endogenous control molecules that belong to a particular epigenetic state that is partitioned in at least one partitioning set, relative to the total number of endogenous control molecules that belong to that epigenetic state in all partitioning sets, relative to the fraction or percentage. In some embodiments, the epigenetic compartment score is determined for each epigenetic state of the epigenetic control nucleic acid molecules and endogenous control molecules in each of the partitioning sets. In some embodiments, the epigenetic compartment score is determined for the epigenetic control nucleic acid molecules and endogenous control molecules that have one or more particular epigenetic states in one or more partitioning sets. In some embodiments, an epigenetic compartment score is determined for epigenetic control nucleic acid molecules and endogenous control molecules having a particular epigenetic state within a particular compartmentalized set.

[0145] In some embodiments, the epigenetic compartment score may address the effectiveness with which molecules with no CG ("zero" CG) are partitioned into the highly compartmentalized set. This score may be referred to as a 0 CG score. In some embodiments, the 0 CG score may be expressed as a fraction or percentage of molecules with no CG in the highly compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the fraction of epigenetic control nucleic acid molecules and / or the fraction of hypermethylated control molecules that have at least one of the following: In the hypermethylated compartmentalization set (i.e., the highly methylated compartmentalization set), (vi) 1 methyl CG (the epigenetic compartment score can be referred to as 1 CG score); (vii) two methyl CGs (the epigenetic compartment score can be referred to as the 2 CG score); (viii) three methyl CGs (the epigenetic compartment score can be referred to as the 3 CG score); (ix) 4 methyl CGs (the epigenetic compartment score can be referred to as the 4 CG score); and (x) 5 methyl CGs (the epigenetic compartment score can be referred to as the 5 CG score).

[0146] In some embodiments, the epigenetic compartment score may be directed to the effectiveness of the hypomethylated control molecules or hypomethylated epigenetic control nucleic acid molecules being partitioned into the hypermethylated compartmentalized set. This score may be referred to as the hypo score. In some embodiments, the hypo score may be expressed as a fraction or percentage of the hypomethylated control molecules or hypomethylated epigenetic control nucleic acid molecules in the hypermethylated compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to achieve less than 5% of the hypermethylated control molecules and / or hypermethylated epigenetic control nucleic acid molecules in the hypomethylated compartmentalized set. This score may be referred to as the methyl-5. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to achieve at least 50% of the hypermethylated control molecules and / or hypermethylated epigenetic control nucleic acid molecules in the hypermethylated compartmentalized set. This score may be referred to as the methyl-half.

[0147] For example, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides. The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalization sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. For each subset, an epigenetic compartment score is determined for each of the compartmentalization sets (P1, P2 and P3). That is, the epigenetic control nucleic acid molecules belonging to subset A have three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. Similarly, each of subsets B and C has three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. An epigenetic compartment score can also be determined for endogenous control molecules.

[0148] In another embodiment, three subsets of epigenetic control nucleic acid molecules (subsets A, B and C) are used, each subset having a different number of methylated nucleotides (i.e., each subset has a different epigenetic state). The epigenetic control nucleic acid molecules in these three subsets can be partitioned into three compartmentalized sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. In this embodiment, the epigenetic score is determined only for the subset A molecules in the P1 compartmentalized set. The epigenetic score can be a fraction or percentage measure of the subset A molecules in the P1 compartmentalized set relative to the total number of subset A molecules (in the P1, P2 and P3 compartmentalized sets).

[0149] The epigenetic compartment score can be any value or range between 0 and 1 (by fraction) or between 0 and 100% (by percentage). In some embodiments, the epigenetic compartment score can be by the number of methylated CGs (e.g., for methyl-half and methyl-5).

[0150] In 306, the epigenetic compartment scores of the epigenetic control nucleic acid molecule and the endogenous control molecule are compared to their corresponding epigenetic compartment cutoffs (predetermined cutoffs) to evaluate the partitioning method. The epigenetic compartment cutoff is a predetermined cutoff value or cutoff range used to evaluate the compartments of the nucleic acid molecules belonging to a particular epigenetic state, and each partitioning set has an epigenetic compartment cutoff for the nucleic acid molecules belonging to the epigenetic state. The epigenetic compartment cutoff varies depending on the epigenetic state of the nucleic acid molecule and the partitioning set, i.e., each epigenetic state has its own epigenetic compartment cutoff, and any of the partitioning sets has a separate epigenetic compartment cutoff for its epigenetic state. The cutoffs can be by percentage or fraction, and the cutoffs can be a cutoff range rather than a specific cutoff value. For example, the epigenetic compartment cutoffs for epigenetic control nucleic acid molecules belonging to a particular epigenetic state can be between 70% and 79%, between 10% and 15%, and less than 5% for the compartmentalization sets P1, P2, and P3, respectively. If the epigenetic compartment scores of the epigenetic control nucleic acid molecules belonging to that epigenetic state fall within the range of the corresponding epigenetic compartment cutoffs, the compartmentalization method is successful.

[0151] In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.01%, 0.02%, 0.05%, 0.1%, 0.2%, 0.3%, 0.4%, 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1%, 2%, 5%, at least 5% or at least 10%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.01%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.02%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.03%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.04%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.05%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.2%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.3%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.4%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.6%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.7%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.8%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.9%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 1%.

[0152] In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%, 0.5%, 1%, 2%, 3%, 4%, 5%, 7%, or at least 10%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 2%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 3%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 4%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 5%.

[0153] In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5, 10, 20, 30, 40 or 50 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 40 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 50 mCG.

[0154] In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5, 10, 15, 20, 25, 30, 35 or 40 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 15 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 25 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 35 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 40 mCG.

[0155] In some embodiments, if the epigenetic compartment scores of one or more epigenetic control nucleic acid molecules and endogenous control molecules belonging to one or more epigenetic states in one or more compartmentalization sets fall within the range of the corresponding epigenetic compartment cutoffs, the compartmentalization method can be classified as successful. Otherwise, the compartmentalization method can be classified as unsuccessful.

[0156] In another aspect, the present disclosure provides a method for evaluating the epigenetic status-based partition of nucleic acid molecules in a polynucleotide sample, comprising: (a) partitioning at least one subset of nucleic acid molecules of the polynucleotide sample into a plurality of partitioning sets; (c) enriching at least one subset of molecules from the plurality of partitioning sets to generate a set of enriched molecules, wherein the set of enriched molecules comprises a group of nucleic acid molecules from the polynucleotide sample, and the group of nucleic acid molecules from the cell-free polynucleotide sample comprises a set of endogenous control molecules; (d) sequencing at least one subset of the set of enriched molecules to generate a set of sequencing reads; (e) analyzing at least one subset of the set of sequencing reads to generate one or more epigenetic partition scores for the set of endogenous control molecules; and (f) comparing the one or more epigenetic partition scores to one or more epigenetic partition cutoffs. In these embodiments, the partitioning of the nucleic acid molecules of the sample and the epigenetic control nucleic acid molecules is always performed simultaneously. In some embodiments, the analyzing step comprises estimating the number / fraction of endogenous control molecules that are in a given epigenetic state in at least one of the compartmentalized sets.

[0157] FIG. 4 illustrates an example embodiment of a method 400 for evaluating partitions based on epigenetic states of nucleic acid molecules in a polynucleotide sample. In this embodiment, partitions of endogenous control molecules in a polynucleotide sample are analyzed to evaluate the partition method. There are regions in the human genome that have a specific epigenetic state, and the epigenetic state of the region often does not vary / change, but always remains the same / remains consistent even in different subjects and / or different disease types / stages. Nucleic acid molecules in the polynucleotide sample that correspond to such human genome regions with an invariant epigenetic state are referred to as endogenous control molecules. In 401, a polynucleotide sample from a subject is taken into consideration. In 402, the nucleic acid molecules of at least a subset of the polynucleotide sample are partitioned or separated into a plurality of partitioned sets based on the epigenetic state of the molecules. The partitions can be based on the presence or absence of epigenetic modifications and / or based on the degree of epigenetic modifications. Examples of epigenetic modifications include, but are not limited to, the presence or absence of methylation, the level of methylation, and the type of methylation (5' cytosine). In some embodiments, the epigenetic modification can be DNA methylation. In these embodiments, the molecules of the spike-in sample are partitioned based on different levels of methylation (number of differently methylated nucleotides). In some embodiments, the spike-in sample can be partitioned into two or more partition sets (e.g., at least 3, 4, 5, 6, or 7 partition sets). In some embodiments, the partitions are based on the differential binding affinity of the nucleic acid molecules to the binding agent.

[0158] The compartments of nucleic acid molecules can be analyzed by sequencing the compartmented nucleic acid molecules or by digital droplet PCR (ddPCR). Before analyzing the compartments, the nucleic acid molecules in the compartmentalized set can be enriched, thereby increasing the signal from the nucleic acid molecule of interest and thus improving the sensitivity. In 403, at least one subset of the nucleic acid molecules in the multiple compartmentalized sets is enriched, thus enriching the endogenous control molecules (from the polynucleotide sample) and the nucleic acid molecules from other polynucleotide samples that belong to the region of interest.

[0159] In some embodiments, each of the multiple compartmentalization sets is differentially tagged before enrichment. The tagged compartmentalization sets are then pooled together for collective sample preparation and / or sequencing. Differential tagging of the compartmentalization sets helps to keep track of the nucleic acid molecules that belong to a particular compartmentalization set. The tags are usually provided as components of the adapter. The nucleic acid molecules in different compartmentalization sets receive different tags that can distinguish the members of one compartmentalization set from the members of another compartment cassette. The tags linked to the nucleic acid molecules of the same compartment set can be the same or different from each other. However, when different from each other, the tags can have a part of their sequence in common so that the molecules to which they are attached can be identified as belonging to a particular compartmentalization set.

[0160] At least one subset of the enriched molecules is sequenced at 404. The sequence information obtained includes the sequence of the nucleic acid molecule and the tag attached to the nucleic acid molecule. From the sequence of the tag attached to the nucleic acid molecule, the tag can be correlated with the compartmentalized set of the nucleic acid molecule. The sequence information is used to identify the endogenous control molecules and their corresponding compartmentalized set. This information is used to analyze the compartments of the endogenous control molecules. At 405, one or more epigenetic compartment scores of the endogenous control molecules belonging to one or more compartmentalized sets are determined. In some embodiments, the sensitivity and / or specificity of the compartmentalization method can be assessed by the epigenetic compartment score. The epigenetic compartment score is a score that represents the compartment of the nucleic acid molecule that belongs to a particular epigenetic state. In some embodiments, the epigenetic compartment score of the nucleic acid molecule that belongs to the epigenetic state is determined for each compartmentalized set. For example, the epigenetic compartment score of the endogenous control molecule that belongs to a particular epigenetic state can be determined. The epigenetic compartment score can be a measure of the number (or a statistically estimated number) of nucleic acid molecules that belong to a particular epigenetic state. The epigenetic compartment score can be a fraction or percentage. The epigenetic compartment score can be a measure of the ratio of the number of endogenous control molecules that belong to a particular epigenetic state that are partitioned in at least one partitioning set to the number of endogenous control molecules that belong to that epigenetic state that are present in the other remaining partitioning set(s). In some embodiments, the epigenetic compartment score can be a fraction or percentage of the number of endogenous control molecules that belong to a particular epigenetic state that are partitioned in at least one partitioning set to the total number of endogenous control molecules that belong to that epigenetic state in all partitioning sets. In some embodiments, the epigenetic compartment score is determined for each epigenetic state of the endogenous control molecules in each of the partitioning sets. In some embodiments, the epigenetic compartment score is determined for endogenous control molecules that have one or more particular epigenetic states in one or more partitioning sets.In some embodiments, the epigenetic compartment score is determined with respect to an endogenous control molecule having a particular epigenetic state within a particular compartmentalization set.

[0161] In some embodiments, the epigenetic compartment score may address the effectiveness with which molecules with no CG ("zero" CG) are partitioned into the highly compartmentalized set. This score may be referred to as a 0 CG score. In some embodiments, the 0 CG score may be expressed as a fraction or percentage of molecules with no CG in the highly compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the fraction of highly methylated control molecules that have at least one of the following: In the hypermethylated compartmentalization set (i.e., the highly methylated compartmentalization set), (xi) 1 methyl CG (the epigenetic compartment score can be referred to as 1 CG score); (xii) 2 methyl CGs (the epigenetic compartment score can be referred to as the 2 CG score); (xiii) three methyl CGs (the epigenetic compartment score can be referred to as the 3 CG score); (xiv) 4 methyl CGs (the epigenetic compartment score can be referred to as the 4 CG score); and (xv) 5 methyl CGs (the epigenetic compartment score can be referred to as the 5 CG score).

[0162] In some embodiments, the epigenetic compartment score may be directed to the effectiveness of the hypomethylated control molecules being partitioned into the hypermethylated compartmentalized set. This score may be referred to as the hypo score. In some embodiments, the hypo score may be expressed as a fraction or percentage of the hypomethylated control molecules in the hypermethylated compartmentalized set. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to have less than 5% of the hypermethylated control molecules in the hypomethylated compartmentalized set. This score may be referred to as the methyl-5. In some embodiments, the epigenetic compartment score may be a measure of the number of methylated CGs required to have at least 50% of the hypermethylated control molecules in the hypermethylated compartmentalized set. This score may be referred to as the methyl-half.

[0163] For example, two subsets of endogenous control molecules (subsets A and B) are analyzed, each subset differing in the level / degree of methylation (i.e., each subset has a different epigenetic state). The endogenous control molecules in these two subsets can be partitioned into three compartmentalization sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. For each subset, an epigenetic compartment score is determined for each of the compartmentalization sets (P1, P2 and P3). That is, the epigenetic control nucleic acid molecules belonging to subset A have three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3. Similarly, subset B has three epigenetic compartment scores, i.e., one epigenetic compartment score for each of the three compartmentalization sets, P1, P2 and P3.

[0164] In another embodiment, three subsets of endogenous control molecules (subsets A, B and C) are analyzed, each with a different level / degree of methylation (i.e., each subset has a different epigenetic state). The endogenous control molecules in these three subsets can be partitioned into three compartmentalized sets, P1, P2 and P3, based on their binding affinity to methyl-binding proteins. In this embodiment, an epigenetic score is determined only for endogenous molecules of subset A in the P1 compartmentalized set. This epigenetic score can be a fraction or percentage measure of the endogenous control molecules of subset A in the P1 compartmentalized set relative to the total number of endogenous control molecules of subset A (in the P1, P2 and P3 compartmentalized sets).

[0165] The epigenetic compartment score can be any value or range between 0 and 1 (by fraction) or between 0 and 100% (by percentage). In some embodiments, the epigenetic compartment score can be by the number of methylated CGs (e.g., for methyl-half and methyl-5).

[0166] In 406, the epigenetic compartment scores of the endogenous control molecules are compared to their corresponding epigenetic compartment cutoffs (predetermined cutoffs) to evaluate the partitioning method. The epigenetic compartment cutoffs are predetermined cutoff values ​​or cutoff ranges used to evaluate the partitions of the nucleic acid molecules belonging to a particular epigenetic state, and each partitioning set has an epigenetic compartment cutoff for the nucleic acid molecules belonging to the epigenetic state. The epigenetic compartment cutoffs vary depending on the epigenetic state of the nucleic acid molecules and the partitioning sets, i.e., each epigenetic state has its own epigenetic compartment cutoff, and any of the partitioning sets has a separate epigenetic compartment cutoff for its epigenetic state. The cutoffs can be by percentage or fraction, and the cutoffs can be a cutoff range rather than a specific cutoff value. For example, the epigenetic compartment cutoffs for endogenous control molecules belonging to a particular epigenetic state may be between 70% and 79%, between 10% and 15%, and less than 5% for compartmentalization sets P1, P2, and P3, respectively. If the epigenetic compartment scores of endogenous control molecules belonging to that epigenetic state fall within the range of the corresponding epigenetic compartment cutoffs, the compartmentalization method is successful.

[0167] In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.01%, 0.02%, 0.05%, 0.1%, 0.2%, 0.3%, 0.4%, 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1%, 2%, 5%, at least 5% or at least 10%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.01%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.02%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.03%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.04%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.05%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.2%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.3%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.4%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.6%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.7%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.8%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 0.9%. In some embodiments, the epigenetic compartment cutoff for a 0 CG score may be 1%.

[0168] In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%, 0.5%, 1%, 2%, 3%, 4%, 5%, 7%, or at least 10%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 0.5%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 1%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 2%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 3%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 4%. In some embodiments, the epigenetic compartment cutoff for the hypo score may be 5%.

[0169] In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5, 10, 20, 30, 40 or 50 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 40 mCG. In some embodiments, the epigenetic compartment cutoff for methyl-5 may be 50 mCG.

[0170] In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5, 10, 15, 20, 25, 30, 35 or 40 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 5 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 10 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 15 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 20 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 25 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 30 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 35 mCG. In some embodiments, the epigenetic compartment cutoff for the methyl-half score may be 40 mCG.

[0171] In some embodiments, if one or more epigenetic compartment scores of endogenous control molecules belonging to one or more epigenetic states in one or more compartmentalization sets fall within the range of the corresponding epigenetic compartment cutoffs, the compartmentalization method can be classified as successful. Otherwise, the compartmentalization method can be classified as unsuccessful.

[0172] In another aspect, the disclosure provides a method for determining an epigenetic state of a nucleic acid molecule(s) in a sample of polynucleotides, the method comprising: (a) adding a set of epigenetic control nucleic acid molecules to nucleic acid molecules in the sample of polynucleotides, thereby creating a spike-in sample; (b) partitioning at least a subset of nucleic acid molecules of the spike-in sample into a plurality of compartmentalized sets; and (c) enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, wherein the enriched set of molecules is a set of epigenetic control nucleic acid molecules. The present invention provides a method for estimating a probability for an epigenetic state of a sample nucleic acid molecule using a plurality of compartmentalization profiles of epigenetic control nucleic acid molecules in different epigenetic states in a plurality of compartmentalization sets. In these embodiments, the compartmentalization of the sample nucleic acid molecules and the epigenetic control nucleic acid molecules is necessarily performed simultaneously.

[0173] In some embodiments, the analyzing step includes determining the number or fraction of epigenetic control nucleic acid molecules per epigenetic state in the multiple compartmentalized sets. A compartment profile may refer to a representation of the fraction / number of epigenetic control nucleic acid molecules in each epigenetic state in two or more compartmentalized sets. In some embodiments, the compartment profile further includes information on the number of nucleotides with epigenetic modifications in the epigenetic control nucleic acid molecule, the position of the nucleotides with epigenetic modifications in the epigenetic control nucleic acid molecule and / or the sequence composition of the epigenetic control nucleic acid molecule. This compartment profile can be used to estimate the probability for the epigenetic state of the nucleic acid molecule in the sample. In some embodiments, if the epigenetic modification is methylation, the compartment profile can be used to estimate the probability for the methylation state (i.e., the level / degree of methylation or the number of methylated nucleotides) of the nucleic acid molecule of the sample.

[0174] In another aspect, the disclosure provides a method for determining an epigenetic state of a nucleic acid molecule(s) in a polynucleotide sample, the method comprising: (a) partitioning nucleic acid molecules from at least a subset of the sample into a plurality of compartmentalized sets; (b) enriching at least a subset of molecules from the plurality of compartmentalized sets to generate a set of enriched molecules, the set of enriched molecules comprising a group of nucleic acid molecules from the polynucleotide sample, the group of nucleic acid molecules from the polynucleotide sample comprising a set of endogenous control molecules; (c) sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; (e) analyzing at least a subset of the set of sequence reads to generate a plurality of compartment profiles of endogenous control molecules in epigenetic states in different plurality of compartmentalized sets; and (f) estimating a probability for the epigenetic state of the nucleic acid molecule using the plurality of compartmentalized profiles of the endogenous control molecules.

[0175] In some embodiments, the analyzing step includes determining the number of endogenous control molecules per epigenetic state in the multiple compartmentalized sets. The compartment profile may refer to a representation of the fraction / number of endogenous control molecules in each epigenetic state in two or more compartmentalized sets. In some embodiments, the compartment profile further includes information regarding the number of nucleotides with epigenetic modifications in the epigenetic control nucleic acid molecule, the position of the nucleotides with epigenetic modifications in the epigenetic control nucleic acid molecule and / or the sequence composition of the epigenetic control nucleic acid molecule. In some embodiments, the number of methylated CpGs in the endogenous control molecule is determined based on previous experimental data and / or from the literature. This compartment profile can be used to estimate the probability for the epigenetic state of the nucleic acid molecule in the sample. In some embodiments, if the epigenetic modification is methylation, the compartment profile can be used to estimate the probability for the methylation state (i.e., the level / degree of methylation or the number of methylated nucleotides) of the nucleic acid molecule of the sample.

[0176] In some embodiments, endogenous control molecules (e.g., high methylation control molecules and low methylation control molecules) can be used to estimate the methylation status of the nucleic acid molecules of the sample. When there are three compartmentalization sets, P1, P2 and P3, compartment profiles of the high methylation control molecules can be generated for P1, P2 and P3 based on the fraction of the high methylation control molecules in each of the three compartmentalization sets and the number of methylated CpGs present in the high methylation control molecules. Similarly, for the low methylation control molecules, compartment profiles of the low methylation control molecules can be generated for P1, P2 and P3 based on the fraction of the low methylation control molecules in each of the three compartmentalization sets and the number of unmethylated CpGs present in the low methylation control molecules. In some embodiments, when endogenous control molecules are used, the number of methylated CpGs in the endogenous control molecules is determined based on previous experimental data and / or from literature. These six compartment profiles can be used to estimate the probability regarding the level / degree of methylation or the number of methylated nucleotides present in a particular region in the nucleic acid molecules of the sample.

[0177] In some embodiments, epigenetic control nucleic acid molecules (e.g., highly methylated epigenetic control nucleic acid molecules and low methylated epigenetic control nucleic acid molecules) can be used to estimate the methylation state of the nucleic acid molecules of the sample. When there are three compartmentalization sets, P1, P2 and P3, compartment profiles of highly methylated epigenetic control nucleic acid molecules can be generated for P1, P2 and P3 based on the fraction of highly methylated epigenetic control nucleic acid molecules in each of the three compartmentalization sets and the number of methylated CpGs present in the highly methylated epigenetic control nucleic acid molecules. Similarly, for low methylated epigenetic control nucleic acid molecules, compartment profiles of low methylated epigenetic control nucleic acid molecules can be generated for P1, P2 and P3 based on the fraction of low methylated epigenetic control nucleic acid molecules in each of the three compartmentalization sets and the number of unmethylated CpGs present in the low methylated epigenetic control nucleic acid molecules. These six compartment profiles can be used to estimate the probability as to the level / degree of methylation or the number of methylated nucleotides present in a particular region within a sample nucleic acid molecule. II. Epigenetic Control Nucleic Acid Molecules

[0178] Epigenetic control nucleic acid molecules are used as control or reference molecules to evaluate the compartments of nucleic acid molecules in a sample based on epigenetic modifications. These epigenetic control nucleic acid molecules can also be used to determine the epigenetic state of the nucleic acid molecule(s) in the sample. For example, the epigenetic modification can be DNA methylation, and the epigenetic control nucleic acid molecules can have different / distinguishable levels of methylation. The epigenetic control nucleic acid molecules can be synthetic oligonucleotides. In some embodiments, the epigenetic control nucleic acid molecules can have nucleic acid sequences that do not occur in nature. In some embodiments, the epigenetic control nucleic acid molecules can have nucleic acid sequences that occur in nature. In some embodiments, the epigenetic control nucleic acid molecules can have nucleic acid sequences that correspond to non-human genomes. For example, these molecules can have either (i) sequences that correspond to regions of lambda phage DNA or the human genome, (ii) sequences that do not occur in nature, and / or (iii) a combination of (i) and (ii). Additionally, epigenetic control nucleic acid molecules can be grouped into subsets, each of which can have a particular number of nucleotides representing a degree of epigenetic modification that is different from the number of nucleotides representing a degree of epigenetic modification in any other set.

[0179] In another aspect, the present disclosure provides a set of epigenetic control nucleic acid molecules comprising one or more subsets of epigenetic control nucleic acid molecules, each subset comprising a plurality of epigenetic control nucleic acid molecules, each epigenetic control nucleic acid molecule comprising an epigenetic modification region. The epigenetic modification region is a region of the epigenetic control nucleic acid molecule that represents the epigenetic state of the epigenetic control nucleic acid molecule. The epigenetic state is the level / degree of epigenetic modification of a nucleic acid molecule. For example, when the epigenetic modification is DNA methylation, the epigenetic state can refer to a highly methylated, low methylated or moderately methylated nucleic acid molecule. The epigenetic state can also refer to the number of nucleotides that have epigenetic modification. For example, when the epigenetic modification is DNA methylation, the epigenetic state can refer to the number of methylated nucleotides of a nucleic acid molecule.

[0180] In some embodiments, the epigenetic control nucleic acid molecule comprises at least one of the following: (i) an epigenetic modification region and (ii) an identifier region. In some embodiments, the epigenetic modification region comprises nucleotides having an epigenetic modification. In some embodiments, the epigenetic modification is DNA methylation. In these embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule may have methylated nucleotides. The number of methylated nucleotides in the epigenetic modification region may vary between epigenetic control nucleic acid molecules. In some embodiments, the epigenetic control nucleic acid molecule may have 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, at least 10, at least 15, at least 20, at least 30, at least 40 or at least 50 methylated nucleotides in the epigenetic modification region. The epigenetic control nucleic acid molecules can be grouped into subsets based on the epigenetic state (i.e., the number of nucleotides with epigenetic modification / level of epigenetic modification) in the epigenetic modification region. The length of the epigenetic modification region in different subsets can be the same, for example, approximately 160bp. The length of the epigenetic modification region between subsets can be different. For example, the epigenetic control nucleic acid molecules can be grouped into three subsets (subsets A, B and C) based on the number of methylated nucleotides in the epigenetic modification region. Subsets A, B and C can have epigenetic control nucleic acid molecules with 5, 10 and 15 methylated nucleotides in the epigenetic modification region, respectively, and the length of the epigenetic modification region in subsets A, B and C can be the same (for example, 160bp) or can be different - 100bp, 150bp and 200bp for subsets A, B and C, respectively.

[0181] In certain embodiments, epigenetic control nucleic acid molecules can be grouped into subsets, each subset representing a degree of epigenetic modification, and the number of polynucleotides in each subset is different from the number of nucleotides in every other set. In some embodiments, the number of methylated nucleotides in a subset is 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, at least 12, at least 15, at least 20, at least 25, at least 30, at least 40, or at least 50. In some embodiments, the epigenetic modification region of the epigenetic control nucleic acid molecule in at least one subset comprises at least one nucleotide with epigenetic modification. In some embodiments, the subset of at least one epigenetic control nucleic acid molecule can comprise nucleotides that do not have any epigenetic modification (i.e., nucleotides that are not epigenetically modified) in the epigenetic modification region of the epigenetic control nucleic acid molecule.

[0182] In some embodiments, the epigenetic modification regions of all epigenetic control nucleic acid molecules in the subset contain the same number of nucleotides with epigenetic modification. In some embodiments, the number of nucleotides with epigenetic modification in the first subset is different from the number of nucleotides with epigenetic modification in the second subset. In some embodiments, the epigenetic modification regions of the epigenetic control nucleic acid molecules in one or more subsets contain the same nucleic acid sequence. In some embodiments, the epigenetic modification regions of the epigenetic control nucleic acid molecules in the first subset contain nucleic acid sequences that are distinguishable from the nucleic acid sequences of the epigenetic modification regions of the epigenetic control nucleic acid molecules in the second subset.

[0183] In some embodiments, the epigenetic modification regions of the epigenetic control nucleic acid molecules in one or more subsets are the same length and have the same sequence composition, but the number of nucleotides with epigenetic modification can be different in each of one or more subsets.In some embodiments, the epigenetic modification regions of the epigenetic control nucleic acid molecules in one or more subsets are the same length and have the same number of nucleotides with epigenetic modification, but the position of the nucleotide with epigenetic modification can be different in each of one or more subsets.In some embodiments, the epigenetic modification regions of the epigenetic control nucleic acid molecules in one or more subsets are the same length, have the same number of nucleotides with epigenetic modification, and the position of the nucleotide with epigenetic modification can be the same, but the nucleotides adjacent to both sides of the nucleotide with epigenetic modification can be different in each of one or more subsets.

[0184] In some embodiments, each subset of epigenetic control nucleic acid molecules is present in equimolar concentrations. In some embodiments, each subset of epigenetic control nucleic acid molecules is present in non-equimolar concentrations. In some embodiments, the epigenetic modification is DNA methylation. In some embodiments, the nucleotide having the epigenetic modification comprises a methylated nucleotide. In some embodiments, the methylated nucleotide comprises 5-methylcytosine. In some embodiments, the methylated nucleotide comprises 5-hydroxymethylcytosine. In some embodiments, the methylated nucleotide comprises N 6 -Contains methyladenine.

[0185] In some embodiments, the epigenetic control nucleic acid molecule further comprises an identifier region. The identifier region is a region of the epigenetic control nucleic acid molecule that is used to distinguish the epigenetic control nucleic acid molecule from other epigenetic control nucleic acid molecules. The identifier region may have a molecular barcode and / or an epigenetic state barcode. The identifier region may be present on one or both sides of the epigenetic modification region. The molecular barcode serves as an identifier of the epigenetic control nucleic acid molecule, while the epigenetic state barcode serves as an identifier of the epigenetic state of the epigenetic control nucleic acid molecule. The epigenetic state barcode is a type of barcode (nucleic acid sequence) that is used to identify the epigenetic state of the epigenetic control nucleic acid molecule. In some embodiments, the epigenetic state barcode is one that can identify (by a predetermined correlation) the number of nucleotides with epigenetic modification in the epigenetic modification region of the epigenetic control nucleic acid molecule. In some embodiments, the epigenetic state barcode can identify the level of epigenetic modification in the epigenetic modification region of the epigenetic control nucleic acid molecule. In some embodiments, the identifier region of the epigenetic control nucleic acid molecule comprises an epigenetic state barcode. For example, if the epigenetic modification is DNA methylation and a subset of the epigenetic control nucleic acid molecules has 5 methylated nucleotides, then all of the epigenetic control nucleic acid molecules in the subset will have the same epigenetic state barcode. In some embodiments, the epigenetic state barcode can be used to identify the level / degree of epigenetic modification in the epigenetic modification region of the epigenetic control nucleic acid molecule. The epigenetic control nucleic acid molecules can be grouped into subsets based on the number of cytosine or CpG nucleotides in the epigenetic modification region.In some embodiments, within each subset, the level of methylation may vary (e.g., highly methylated, moderately methylated, and low methylated), and each level of methylation may have a separate epigenetic state barcode. For example, within subset A, the low methylated epigenetic control nucleic acid molecules would all have an epigenetic state barcode, e.g., ESB1, and the highly methylated epigenetic control nucleic acid molecules would all have another epigenetic state barcode, e.g., ESB3. In this example, the epigenetic state barcodes are used to identify the level / degree of methylation. The molecular barcodes within the identifier region may be unique barcodes (each molecule has a unique barcode) or non-unique barcodes. The molecular barcodes may be any length between 2 nucleotides and 50 nucleotides. In some embodiments, the molecular barcodes may be at least 2 nucleotides, at least 3 nucleotides, at least 4 nucleotides, at least 5 nucleotides, at least 6 nucleotides, at least 7 nucleotides, at least 8 nucleotides, at least 9 nucleotides, or at least 10 nucleotides. In some embodiments, the epigenetic state barcode can be at least 2 nucleotides, at least 3 nucleotides, at least 4 nucleotides, at least 5 nucleotides, at least 6 nucleotides, at least 7 nucleotides, or at least 8 nucleotides.

[0186] Figure 5 is a schematic diagram of epigenetic control nucleic acid molecules suitable for use in some embodiments of the present disclosure. The epigenetic control nucleic acid molecules described herein have a similar length to the epigenetic control nucleic acid molecules of the sample being assayed, and all of the subsets have the same sequence composition to reduce the effect of any sequence-specific compartment. In Figure 5, as an example, the epigenetic control nucleic acid molecules are grouped into four subsets, subsets 1, 2, 3 and 4. The epigenetic control nucleic acid molecules in Figure 5 are double-stranded DNA molecules. For illustrative purposes, only one representative of the epigenetic control nucleic acid molecules in each subset is shown in the figure. In this embodiment, the sequence of the epigenetic modification region of the epigenetic control nucleic acid molecule is the same in all of the subsets. The epigenetic modification region of the epigenetic control nucleic acid molecule in all four subsets has five CpG pairs. The '---' region in the double-stranded DNA sequence represents any other sequence apart from the CpG pairs, M represents 5-methylcytosine, C represents cytosine, and G represents guanine. In FIG. 5, the epigenetic state (level of methylation) of the epigenetic control nucleic acid molecules in a subset is different from the epigenetic state of the other subsets. Subset 1 has zero methylated CpG pairs, subset 2 has one fully methylated CpG pair, subset 3 has three fully methylated CpG pairs, and subset 4 has five fully methylated CpG pairs. In this embodiment, the identifier regions are present on both sides of the epigenetic modification region. The identifier regions on both sides have an epigenetic state barcode (ESB), while the molecular barcode (MB) is only on one side. Molecular barcodes are used as identifiers for epigenetic control nucleic acid molecules, with each epigenetic control nucleic acid molecule having a unique molecular barcode (i.e., molecule 1 has MB1, molecule 2 has MB2, molecule 3 has MB3, etc.). Epigenetic state barcodes can be used as identifiers of the epigenetic state of epigenetic control nucleic acid molecules, where the epigenetic state barcodes are used to identify the number of fully methylated CpG pairs in the epigenetic control nucleic acid molecule.All epigenetic control nucleic acid molecules in subset 1 have zero methylated CpG pairs, and therefore all epigenetic control nucleic acid molecules in subset 1 have the same epigenetic state barcode-ESB1. Similarly, all epigenetic control nucleic acid molecules in subsets 2, 3 and 4 have one, three and five fully methylated CpG pairs, respectively. Thus, all epigenetic control nucleic acid molecules in subsets 2, 3 and 4 have epigenetic state barcodes that are ESB2, ESB3 and ESB4, respectively. In this example, the same epigenetic state barcode is present on both sides of the epigenetic modification region.

[0187] In some embodiments, the molecular barcodes may be on one or both sides of the epigenetic modification region. In some embodiments, the epigenetic state barcodes may be on one or both sides of the epigenetic modification region. In some embodiments, the epigenetic state barcodes on both sides of the epigenetic modification region may be the same or different and / or may be attached randomly.

[0188] In some embodiments, the identifier region may have an additional region (primer binding site) that facilitates binding of one or more primers. In some embodiments, the identifier region primer binding site in one subset is different from the primer binding site in the other subset. In some embodiments, if the epigenetic control nucleic acid molecules have different epigenetic states within a subset, the primer binding site may be different for each epigenetic state within the molecule, i.e., each unique epigenetic state has a unique primer binding site. In some embodiments, these primer binding sites are used to analyze the sections of the epigenetic control nucleic acid molecule. In some embodiments, instead of analyzing the sections of the epigenetic control nucleic acid molecule by sequencing, the sections of the epigenetic control nucleic acid molecule may be analyzed by digital droplet PCR (ddPCR) using primers that bind to these primer binding states.

[0189] In some embodiments, the epigenetic control nucleic acid molecules can be grouped into subsets such that the epigenetic control nucleic acid molecules within each subset have a sequence, but the epigenetic state of the epigenetic control nucleic acid molecules within each subset can vary.

[0190] FIG. 6 is a schematic diagram of an epigenetic control nucleic acid molecule that may be suitable for use in certain embodiments of the present disclosure. The epigenetic control nucleic acid molecule described herein can also take into account the sequence composition and number of CpG pairs between sections of a nucleic acid molecule / the effect of fully methylated CpG pairs. In FIG. 6, as an example, the epigenetic control nucleic acid molecules are grouped into three subsets, subsets 1, 2 and 3. The epigenetic control nucleic acid molecules in FIG. 6 are double-stranded DNA molecules. For the sake of illustration, only one representative of the epigenetic control nucleic acid molecule for every epigenetic state in each subset is shown in the figure. In this embodiment, the epigenetic modification regions of the epigenetic control nucleic acid molecules in subsets 1, 2 and 3 are different in length. The epigenetic modification regions of the epigenetic control nucleic acid molecules in subsets 1, 2 and 3 have one, three and five CpG pairs, respectively. The '---' region in the double-stranded DNA sequence represents any other sequence apart from the CpG pairs, M represents 5-methylcytosine, C represents cytosine, and G represents guanine. In FIG. 6, within each subset, the epigenetic control nucleic acid molecules are in different epigenetic states, for example, low methylation, moderate methylation, and high methylation. The epigenetic control nucleic acid molecules of subset 1 are in two different epigenetic states, namely low methylation (zero methylated CpG pairs) and high methylation (one fully methylated CpG pair). The epigenetic control nucleic acid molecules of subset 2 are in three different epigenetic states, namely low methylation (zero methylated CpG pairs), moderate methylation (one fully methylated CpG pair), and high methylation (three fully methylated CpG pairs). The epigenetic control nucleic acid molecules of subset 3 are in three different epigenetic states, namely, a low methylated state (one fully methylated CpG pair), a moderately methylated state (three fully methylated CpG pairs) and a highly methylated state (five fully methylated CpG pairs), where the identifier regions are present on either side of the epigenetic modification region.The identifier regions on both sides have an epigenetic state barcode (ESB) and a molecular barcode (MB). The molecular barcode is used as an identifier of the epigenetic control nucleic acid molecule, with each epigenetic control nucleic acid molecule having a unique molecular barcode (i.e., molecule 1 has MB1, molecule 2 has MB2, molecule 3 has MB3, etc.). The epigenetic state barcode is used as an identifier of the epigenetic state of the epigenetic control nucleic acid molecule. Here, the epigenetic state barcode is used to identify the degree / level of methylation of the epigenetic control nucleic acid molecule, i.e., whether it is in a low methylated state, a moderately methylated state, or a highly methylated state. All low methylated epigenetic control nucleic acid molecules in subsets 1, 2, and 3 have the same epigenetic state barcode-ESB1. Subsets 2 and 3 have epigenetic control nucleic acid molecules that are moderately methylated, and all of these molecules have the same epigenetic state barcode-ESB2 (subset 1 has no moderately methylated states, and therefore no epigenetic control nucleic acid molecules have the ESB2 epigenetic state barcode). Thus, from the sequence of the epigenetic control nucleic acid molecule and the sequence of the epigenetic state barcode, the epigenetic state of the epigenetic control nucleic acid molecule and the subset to which the epigenetic nucleic acid molecule belongs can be identified.

[0191] Furthermore, the identifier region may have a primer binding site. Different primer binding sites can be used to distinguish different epigenetic states within each subset and between subsets. For example, the low methylated epigenetic control nucleic acid molecule in subset 1 may have primer binding sites-Pr1 and Pr2 on both sides of the epigenetic modification region. The highly methylated epigenetic control nucleic acid molecule in subset 1 may have primer binding sites-Pr3 and Pr4 on both sides of the epigenetic modification region. Similarly, in subset 2, the low, moderate and highly methylated epigenetic control nucleic acid molecules may have primer binding sites Pr5&Pr6, P7&Pr8 and Pr9&Pr19 on both sides of the epigenetic modification region, respectively. Similarly, in subset 3, the low, moderate and highly methylated epigenetic control nucleic acid molecules may have primer binding sites Pr11&Pr12, P13&Pr14 and Pr15&Pr16 on both sides of the epigenetic modification region, respectively. Also, from separate primer sets used for different epigenetic state molecules in different subsets, a measure of the number of epigenetic control nucleic acid molecules belonging to a particular epigenetic state in a particular subset can be estimated by either ddPCR or quantitative PCR (qPCR).In this embodiment, from the epigenetic state barcode sequence and the sequence of the epigenetic modification region, the number of CpG pairs in the epigenetic modification region and the number of fully methylated CpG pairs in the epigenetic modification region can be determined.

[0192] FIG. 7 is a schematic diagram of an epigenetic control nucleic acid molecule suitable for use in some embodiments of the present disclosure. The epigenetic control nucleic acid molecule described herein can take into account the position-specific effect of fully methylated CpG pairs between sections of a nucleic acid molecule. In FIG. 7, the epigenetic control nucleic acid molecules are grouped into five subsets. The sequence length and sequence composition of the epigenetic modification region of the epigenetic control nucleic acid molecule are the same in all subsets. Each subset has two fully methylated CpG pairs, but the positions of the two fully methylated CpG pairs vary from subset to subset (i.e., the distance between the two fully methylated CpG pairs varies from subset to subset). In subset 1, the two fully methylated CpG pairs are far apart, while in subset 4, the two fully methylated CpG pairs are very close to each other. Here, the identifier regions are present on both sides of the epigenetic modification region. The identifier regions on both sides have an epigenetic state barcode (ESB) and a molecular barcode (MB). Molecular barcodes are used as identifiers for individual epigenetic control nucleic acid molecules, and each epigenetic control nucleic acid molecule has a unique molecular barcode, i.e., molecule 1 has MB1, molecule 2 has MB2, molecule 3 has MB3, etc. These subsets have different binding affinities based on the influence of the position of the fully methylated CpG pair. Here, the epigenetic state barcodes can be used to identify the position of the fully methylated CpG pair. All of the epigenetic control nucleic acid molecules in subset 1 have two fully methylated CpG pairs at the same position, and therefore, the epigenetic control nucleic acid molecules in subset 1 have the same epigenetic state barcode-ESB1. Similarly, all of the epigenetic control nucleic acid molecules in subsets 2, 3, 4, and 5 have epigenetic state barcodes, which are ESB2, ESB3, and ESB4, respectively. In this example, the same epigenetic state barcodes are present on both sides of the epigenetic modification region.

[0193] In another aspect, the disclosure provides a population of nucleic acids comprising a set of epigenetic control nucleic acid molecules, the set comprising one or more subsets of epigenetic control nucleic acid molecules, each subset comprising a plurality of epigenetic control nucleic acid molecules, each epigenetic control nucleic acid molecule comprising an epigenetic modified region, and a set of nucleic acid molecules in a sample of polynucleotides from a subject.

[0194] In some embodiments, the epigenetic control nucleic acid molecule can have either (i) a sequence corresponding to a region of lambda phage DNA or the human genome, (ii) a sequence that does not occur in nature, and / or (iii) a combination of (i) and (ii). In some embodiments, the epigenetic control nucleic acid molecule can include a sequence that does not occur in nature.

[0195] In some embodiments, the polynucleotide sample is a DNA sample, an RNA sample, a cell-free polynucleotide sample, a cell-free DNA sample or a cell-free RNA sample. In some embodiments, the polynucleotide sample is a cell-free DNA sample.

[0196] In some embodiments, the cell free DNA is at least 1 ng, at least 5 ng, at least 10 ng, at least 15 ng, at least 20 ng, at least 30 ng, at least 50 ng, at least 75 ng, at least 100 ng, at least 150 ng, at least 200 ng, at least 250 ng, at least 300 ng, at least 350 ng, at least 400 ng, at least 450 ng, or at least 500 ng.

[0197] In some embodiments, the amount of the epigenetic control nucleic acid molecule is at least 1 femtomole, at least 2 femtomole, at least 5 femtomole, at least 10 femtomole, at least 15 femtomole, at least 20 femtomole, at least 50 femtomole, at least 75 femtomole, at least 100 femtomole, at least 125 femtomole, at least 150 femtomole, or at least 200 femtomole. III. General Features of the Method A. Sample

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

[0199] In some embodiments, the sample volume of bodily fluid obtained from the subject depends on the desired depth of reads to the region to be sequenced. Example volumes are about 0.4-40 milliliters (mL), about 5-20 mL, about 10-20 mL. For example, the volume can be about 0.5 mL, about 1 mL, about 5 mL, about 10 mL, about 20 mL, about 30 mL, about 40 mL, or more milliliters. The volume of sampled plasma is generally between about 5 mL and about 20 mL.

[0200] Samples can contain varying amounts of nucleic acid. Generally, the amount of nucleic acid in a given sample is equated to a number of genome equivalents. For example, a sample of about 30 nanograms (ng) of DNA contains about 10,000 (10 4 ), and in the case of cfDNA, about 200 billion (2×10 11 Similarly, a sample of about 100 ng of DNA may contain about 30,000 haploid human genome equivalents, or in the case of cfDNA, about 600 billion individual molecules.

[0201] In some embodiments, the sample comprises nucleic acid from different sources, for example, nucleic acid from cells and nucleic acid from acellular sources (such as blood samples). Generally, the sample comprises nucleic acid carrying mutations. For example, the sample optionally comprises DNA carrying germline mutations and / or somatic mutations. Generally, the sample comprises DNA carrying cancer-associated mutations (e.g., cancer-associated somatic mutations).

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

[0203] The cell-free nucleic acids generally have a size distribution between about 100 and about 500 nucleotides in length, with molecules between about 110 and about 230 nucleotides in length accounting for about 90% of the molecules in a sample, with a mode of about 168 nucleotides in length (in samples from human subjects) and a second minor peak ranging from about 240 to about 440 nucleotides in length. In some embodiments, the cell-free nucleic acids are from about 160 to about 180 nucleotides in length, or from about 320 to about 360 nucleotides in length, or from about 440 to about 480 nucleotides in length.

[0204] In some embodiments, the cell-free nucleic acid is isolated from the body fluid through a partitioning step that separates the cell-free nucleic acid found in solution from intact cells and other insoluble components in the body fluid. In some embodiments, the partitioning includes techniques such as centrifugation or filtration. Alternatively, the cells in the body fluid can be lysed, and the cell-free and intracellular nucleic acids can be processed together. Generally, after the addition of a buffer and a washing step, the cell-free nucleic acid can be precipitated, for example, with alcohol. In some embodiments, an additional cleaning step, such as a silica-based column, is used to remove contaminants or salts. Non-specific bulk carrier nucleic acid is added throughout the reaction, for example, as needed, to optimize aspects of the exemplary procedure, such as yield. After such processing, the sample generally contains various forms of nucleic acid, including double-stranded DNA, single-stranded DNA and / or single-stranded RNA. If necessary, the single-stranded DNA and / or single-stranded RNA is converted to double-stranded form and thus included in subsequent processing and analysis steps. B. Tagging

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

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

[0207] In some embodiments, each sample is uniquely tagged with a sample index or a combination of sample indexes. In some embodiments, the nucleic acid molecules of each sample or subsample are uniquely tagged with a molecular barcode or a combination of molecular barcodes. In other embodiments, multiple molecular barcodes can be used such that the molecular barcodes are not necessarily unique to each other within the plurality (e.g., non-unique molecular barcodes). In these embodiments, the molecular barcodes are generally attached (e.g., by ligation) to individual molecules such that the combination of the molecular barcode and the sequence to which it can be attached creates a unique sequence that can be tracked individually. Detection of the combination of the non-uniquely tagged molecular barcodes and the intrinsic sequence information (e.g., the first (start) and / or last (end) genomic location / location corresponding to the sequence of the original nucleic acid molecule in the sample, the partial sequence of the sequence read at one or both ends, the length of the sequence read, and / or the length of the original nucleic acid molecule in the sample) generally allows a unique identity to be assigned to a particular molecule. In some embodiments, detection of a combination of a non-uniquely tagged molecular barcode and endogenous sequence information (e.g., the first (start) and / or last (end) region of the alignment of the sequence read with the reference sequence, a subsequence of the sequence read at one or both ends, the length of the sequence read, and / or the length of the original nucleic acid molecule in the sample) generally allows a unique identity to be assigned to a particular molecule. In some embodiments, the first region includes a genomic start position of the sequencing read, which determines the 5' end of the sequencing read where alignment to the reference sequence begins, and the last region includes a genomic end position of the sequencing read, which determines the 3' end of the sequencing read where alignment to the reference sequence ends. In some embodiments, the first region includes the first 1 base position, the first 2 base positions, the first 5 base positions, the first 10 base positions, the first 15 base positions, the first 20 base positions, the first 25 base positions, the first 30 base positions, or at least the first 30 base positions of the 5' end of the sequencing read that is aligned with the reference sequence.In some embodiments, the final region includes the last 1 base position, the last 2 base positions, the last 5 base positions, the last 10 base positions, the last 15 base positions, the last 20 base positions, the last 25 base positions, the last 30 base positions or at least the last 30 base positions of the 3' end of the sequencing read that is aligned to the reference sequence.

[0208] Optionally, the length or number of base pairs of each sequence read is used to assign a unique identity to a given molecule. As described herein, fragments derived from a single strand of a nucleic acid having been assigned a unique identity may allow subsequent identification of fragments derived from the parental and / or complementary strands.

[0209] In some embodiments, molecular barcodes are introduced with a set of identifiers (e.g., combinations of unique or non-unique molecular barcodes) and predicted proportions of molecules in a sample. One example format uses about 2 to about 1,000,000 different molecular barcode sequences, or about 5 to about 150 different molecular barcode sequences, or about 20 to about 50 different molecular barcode sequences ligated to both ends of a target molecule. Alternatively, about 25 to about 1,000,000 different molecular barcode sequences can be used. For example, 20-50 x 20-50 molecular barcode sequences (i.e., one of 20-50 different molecular barcode sequences can be attached to each end of a target molecule) can be used. Such a number of identifiers is generally sufficient to provide a high probability (e.g., at least 94%, 99.5%, 99.99%, or 99.999%) that different molecules with the same start and end points will have different combinations of identifiers. In some embodiments, about 80%, about 90%, about 95%, or about 99% of the molecules have the same combination of molecular barcodes.

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

[0211] An epigenetic state barcode (ESB) is a type of tag that is attached to an epigenetic modification region of an epigenetic control nucleic acid molecule. An ESB can be used as an identifier of the epigenetic state of an epigenetic control nucleic acid molecule. An ESB can refer to the number of nucleotides with epigenetic modification in an epigenetic modification region of an epigenetic control nucleic acid molecule. In some embodiments, an identifier region of an epigenetic control nucleic acid molecule comprises at least one epigenetic state barcode. In some embodiments, an ESB is a part of an identifier region of an epigenetic control nucleic acid molecule. For example, if the epigenetic modification is DNA methylation and a subset of epigenetic control nucleic acid molecules has five methylated nucleotides, all of the epigenetic control nucleic acid molecules in the subset will have the same epigenetic state barcode. In some embodiments, an epigenetic state barcode can be used to identify the level / degree of epigenetic modification of an epigenetic modification region of an epigenetic control nucleic acid molecule. The epigenetic control nucleic acid molecules can be grouped into subsets based on the number of cytosines or CpG nucleotides in the epigenetic modification region. In some embodiments, the level of methylation can vary within each subset (e.g., highly methylated, moderately methylated, and low methylated), and each level of methylation can have a separate epigenetic state barcode. For example, within subset A, all low methylated epigenetic control nucleic acid molecules will have an epigenetic state barcode, e.g., ESB1, and all highly methylated epigenetic control nucleic acid molecules will have another epigenetic state barcode, e.g., ESB3. In this example, epigenetic state barcodes are used to identify the level / degree of methylation.

[0212] In some embodiments, assignment of unique or non-unique molecular barcodes in reactions is performed using methods and systems described, for example, in U.S. Patent Application Nos. 20010053519, 20030152490, and 20110160078, as well as U.S. Patent Nos. 6,582,908, 7,537,898, 9,598,731, and 9,902,992, each of which is incorporated by reference herein in its entirety. C. Amplification

[0213] Sample nucleic acid can be amplified by PCR and other amplification methods using nucleic acid primers that are flanked by adaptors and that bind to primer binding sites in the adaptors that flank the DNA molecules to be amplified.In some embodiments, the amplification method involves cycles of extension, denaturation, and annealing, which can be by thermocycling or isothermal, for example, as in transcription amplification methods.Other examples of amplification methods that can be used as needed include ligase chain reaction, strand displacement amplification, nucleic acid sequence-based amplification, and self-sustaining sequence-based replication.

[0214] Generally, the amplification reaction produces a plurality of non-uniquely or uniquely tagged nucleic acid amplicons having molecular barcodes and sample indexes, the sizes of which range from about 150 nucleotides (nt) to about 700 nt, 250 nt to about 350 nt, or about 320 nt to about 550 nt. In some embodiments, the size of the amplicon is about 180 nt. In some embodiments, the size of the amplicon is about 200 nt. D. enrichment

[0215] In some embodiments, sequences are enriched before sequencing the nucleic acids. Enrichment is performed for specific target regions or non-specifically ("target sequences") as required. In some embodiments, targeted regions of interest can be enriched using differential tiling and capture schemes with nucleic acid capture probes ("baits") selected for one or more bait set panels. Differential tiling and capture schemes generally use bait sets at various relative concentrations to differentially tile (e.g., at different "resolutions") across genomic regions associated with the baits, subject to a set of constraints (e.g., sequencer constraints such as sequencing load, availability of each bait, etc.) to capture the targeted nucleic acids at a desired level for downstream sequencing. These targeted genomic regions of interest include natural or synthetic nucleotide sequences of nucleic acid constructs as required. In some embodiments, biotin-labeled beads with probes for one or more regions of interest can be used to capture target sequences, and these regions are then optionally amplified to enrich the regions of interest.

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

[0217] The sample nucleic acid, optionally sandwiched between adapters, is generally subjected to sequencing, with or without pre-amplification.Optionally, the sequencing method or commercially available format used includes, for example, Sanger sequencing, high throughput sequencing, pyrosequencing, sequencing by synthesis, single molecule sequencing, nanopore-based sequencing, semiconductor sequencing, sequencing by ligation, sequencing by hybridization, RNA-Seq (Illumina), Digital Gene Expression (Helicos), next generation sequencing (NGS), single molecule synthesis sequencing (SMSS) (Helicos), massively parallel sequencing, Clonal Single Molecule Array (Solexa), shotgun sequencing, Ion Torrent, Oxford Nanopore, Roche Genia, Maxim-Gilbert sequencing, primer walking, sequencing using PacBio, SOLiD, Ion Torrent or Nanopore platform. Sequencing reactions can be performed on a variety of sample processing devices that may include multiple lanes, multiple channels, multiple wells, or other means of processing multiple sample sets substantially simultaneously. Sample processing devices may also include multiple sample chambers to allow multiple operations to be processed simultaneously.

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

[0219] Simultaneous sequencing reaction can be carried out using multiplex sequencing technique.In some embodiments, the sequencing of cell-free polynucleotide is carried out at least about 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 50000 or 100,000 sequencing reaction.In other embodiments, the sequencing of cell-free polynucleotide is carried out at less than about 1000, less than about 2000, less than about 3000, less than about 4000, less than about 5000, less than about 6000, less than about 7000, less than about 8000, less than about 9000, less than about 10000, less than about 50000 or less than about 100,000 sequencing reaction.Sequencing reaction is generally carried out sequentially or simultaneously. Subsequent data analysis is generally carried out on all or part of the sequencing reaction.In some embodiments, data analysis is carried out on at least about 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 50000, or 100,000 sequencing reactions.In other embodiments, data analysis can be carried out on less than about 1000, less than about 2000, less than about 3000, less than about 4000, less than about 5000, less than about 6000, less than about 7000, less than about 8000, less than about 9000, less than about 10000, less than about 50000, or less than about 100,000 sequencing reactions.Examples of read depth are about 1000 to about 50000 reads per locus (e.g., base position). F. Analysis

[0220] Sequencing generates a plurality of sequencing reads or reads. Sequencing reads or reads can include sequences of nucleotide data that are less than about 150 bases long, or less than about 90 bases long. In some embodiments, reads are between about 80 bases and about 90 bases long, for example, about 85 bases long. In some embodiments, the method of the present disclosure is applied to very short reads, for example, less than about 50 bases long or less than about 30 bases long. Sequencing read data can include sequence data as well as meta information. Sequence read data can be stored in any suitable file format, including, for example, VCF file, FASTA file, or FASTQ file.

[0221] FASTA may refer to a computer program for searching sequence databases, and the name FASTA may also refer to a standard file format. For example, FASTA is described, for example, in Pearson & Lipman, 1988, Improved tools for biological sequence comparison, PNAS 85:2444-2448, which is incorporated herein by reference in its entirety. A sequence in FASTA format begins with a line of description, followed by lines of sequence data. The description line is separated from the sequence data by a greater than (">") symbol in the first column. The word following the ">" symbol is the identifier for the sequence, and the rest of the line is the description (both can be omitted). There should be no space between the ">" and the first character of the identifier. It is recommended that all lines of the string be less than 80 characters. The sequence ends when another line begins with ">", which indicates the beginning of another sequence.

[0222] FASTQ format is a text-based format for storing both biological sequences (usually nucleotide sequences) and their corresponding quality scores. FASTQ format is similar to FASTA format, but with quality scores next to sequence data. For simplicity, both sequence characters and quality scores are coded with a single ASCII character. FASTQ format is the de facto standard for storing the output of high-throughput sequencing instruments such as Illumina Genome Analyzer, as described, for example, in Cock et al. ("The Sanger FASTQ file format for sequences with quality scores, and the Solexa / Illumina FASTQ variants," Nucleic Acids Res38 (6): 1767-1771, 2009), the entire contents of which are incorporated herein by reference.

[0223] For FASTA and FASTQ files, the meta information includes a description line and does not include a line of sequence data. In some embodiments, for FASTQ files, the meta information includes a quality score. For FASTA and FASTQ files, the sequence data begins after the description line and is generally indicated using some subset of the IUPAC ambiguity code, with "-" as needed. In some embodiments, the sequence data can use A, T, C, G, and N characters, and optionally includes "-" or U as needed (e.g., to represent a gap or uracil).

[0224] In some embodiments, at least one master sequence read file and output file are saved as plain text files (e.g., using encoding such as ASCII; ISO / IEC 646; EBCDIC; UTF-8; or UTF-16). The computer system presented in this disclosure may include a text editor program capable of opening plain text files. A text editor program may refer to a computer program that can present the contents of a text file (e.g., a plain text file) on a computer screen and allow a person to edit the text (e.g., using a monitor, keyboard, and mouse). Examples of text editors include, without limitation, Microsoft word, emacs, pico, vi, BBEdit, and TextWrangler. A text editor program is one that can display plain text files on a computer screen and present meta information and sequence reads in a human-readable format (e.g., using alphanumeric characters that can be printed or written by a person, rather than being binary coded).

[0225] Although the methods are discussed in relation to FASTA or FASTQ files, the methods and systems of the present disclosure can be used to compress any suitable sequence file format, including, for example, Variant Call Format (VCF) format files. A typical VCF file can include a header section and a data section. The header contains any number of meta information lines, each beginning with a "##" character, separated by a TAB into field definition lines beginning with a single "#" character. The field definition lines name eight required columns, and the body section contains lines of data that aggregate the columns defined by the field definition lines. The VCF format is described, for example, in Danecek et al. ("The variant call format and VCF tools," Bioinformatics 27(15):2156-2158, 2011), which is incorporated herein by reference in its entirety. The header section can be treated as meta information for writing to the compressed file, and the data section can be treated as lines, each of which is saved in the master file only if it is unique.

[0226] Some embodiments provide for the assembly of sequencing reads. In the assembly by alignment, for example, the sequencing reads are aligned with each other or aligned with a reference sequence. By sequentially aligning each read with a reference genome, all of the reads are positioned relative to each other to create an assembly. Furthermore, the alignment or mapping of sequencing reads with a reference sequence can also be used to identify variant sequences in sequencing reads. The identification of variant sequences can be used in combination with the methods and systems described herein to further assist in the diagnosis or prognosis of disease or condition, or to guide treatment decisions.

[0227] In some embodiments, any or all of the steps are automated. Alternatively, the disclosed method can be fully or partially embodied in one or more dedicated programs, each written in a compiled language such as C++, as appropriate, and then compiled and distributed as a binary. The disclosed method can be fully or partially implemented as a module in an existing sequence analysis platform or by invoking a function in an existing sequence analysis platform. In some embodiments, the disclosed method includes multiple steps, all of which are automatically invoked in response to a single initiation cue (e.g., one or a combination of triggering events provided by a human action, another computer program, or a machine). Thus, the disclosure provides a method in which any or any combination of steps or steps can be performed automatically in response to a cue. "Automatically" generally means without human input, influence, or interaction (e.g., only in response to the original or pre-cue human action).

[0228] The disclosed method may also include various forms of output, including accurate and sensitive interpretation of the nucleic acid sample of interest. The output of the search may be provided in the format of a computer file. In some embodiments, the output is a FASTA file, a FASTQ file, or a VCF file. The output may be processed to create a text file or an XML file that contains sequence data, such as the sequence of the nucleic acid aligned to the sequence of the reference genome. In other embodiments, the processing results in an output that contains coordinates or strings that describe one or more mutations in the nucleic acid of interest compared to the reference genome. The alignment string may include Simple UnGapped Alignment Report (SUGAR), Verbose Useful Labeled Gapped Alignment Report (VULGAR), and Compact Idiosyncratic Gapped Alignment Report (CIGAR) (e.g., as described in Ning et al., Genome Research 11 (10): 1725-9, 2001, which is incorporated herein by reference in its entirety). These strings can be implemented, for example, in the Exonerate sequence alignment software of the European Bioinformatics Institute (Hinxton, UK).

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

[0230] CIGAR strings may follow established motifs. Each letter is preceded by a number, which indicates the base count of the event. Letters used may include M, I, D, N, and S (M=match; I=insertion; D=deletion; N=gap; S=substitution). CIGAR strings define the sequence of matches / mismatches and deletions (or gaps). For example, the CIGAR string 2MD3M2D2M may indicate that the alignment contains 2 matches, 1 deletion (the number 1 is omitted to save some space), 3 matches, 2 deletions, and 2 matches.

[0231] In some embodiments, a nucleic acid population is prepared for sequencing by enzymatically forming blunt ends on double-stranded nucleic acids with single-stranded overhangs at one or both ends. In these embodiments, the population is typically treated with an enzyme having 5'-3' DNA polymerase activity and 3'-5' exonuclease activity in the presence of nucleotides (e.g., A, C, G, and T or U). Examples of enzymes or catalytic fragments that can be used as needed include Klenow large fragment and T4 polymerase. In the 5' overhang, the enzyme typically extends the recessed 3' end on the opposite strand until it lines up with the 5' end, creating a blunt end. In the 3' overhang, the enzyme typically digests from the 3' end to the 5' end of the opposite strand, and sometimes beyond. If this digestion proceeds beyond the 5' end of the opposite strand, the gap can be filled by an enzyme with the same polymerase activity used for the 5' overhang. The formation of blunt ends in double-stranded nucleic acids facilitates, for example, the attachment of adaptors and subsequent amplification.

[0232] In some embodiments, the population of nucleic acids is subjected to additional processing, such as converting single-stranded nucleic acids to double-stranded nucleic acids and / or converting RNA to DNA (e.g., complementary DNA or cDNA). These forms of nucleic acids are also optionally ligated to adapters and amplified.

[0233] Nucleic acid can be subjected to the above-mentioned blunt end forming process with or without prior amplification, and if necessary, other nucleic acids in the sample can be sequenced to generate sequenced nucleic acid.Sequenced nucleic acid can refer to either the sequence (e.g., sequence information) of nucleic acid or the nucleic acid whose sequence is determined.Sequencing can be carried out so that the sequence data of each nucleic acid molecule in the sample is obtained directly or indirectly from the consensus sequence of the amplification product of each nucleic acid molecule in the sample.

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

[0235] In general, the nucleic acid sample is contacted with a sufficient number of adapters so that the probability that any two copies of the same nucleic acid receive the same combination of adapter barcodes from adapters linked to both ends is low (e.g., less than about 1% or less than about 0.1%). This use of adapters may allow for the identification of a family of nucleic acid sequences that have the same start and end points in the reference nucleic acid and are linked to the same combination of barcodes. Such a family may represent the sequence of the amplification product of the nucleic acid in the sample before amplification. As modified by blunt end formation and adapter attachment, the sequences of the family members can be compiled to derive a consensus nucleotide(s) or a complete consensus sequence for the nucleic acid molecules in the original sample. In other words, a nucleotide occupying a particular position of a nucleic acid in a sample can be determined to be the consensus of the nucleotides occupying its corresponding position in the family member sequences. A family may include sequences of one or both strands of a double-stranded nucleic acid. When a family member comprises the sequence of both strands of double-stranded nucleic acid, the sequence of one strand can be converted to its complement for the purpose of compiling the sequence and deriving consensus nucleotide(s) or sequence.Some families only comprise a single member sequence.In this case, this sequence can be considered as the sequence of the nucleic acid in the sample before amplification.Alternatively, families with only a single member sequence can be excluded from further analysis.

[0236] Nucleotide variations (e.g., SNVs or indels) in the sequenced nucleic acid can be determined by comparing the sequenced nucleic acid to a reference sequence. The reference sequence is often a known sequence, such as a known full genome sequence or partial genome sequence from a subject (e.g., a full genome sequence from a human subject). The reference sequence can be, for example, hG19 or hG38. The sequenced nucleic acid represents a sequence determined directly for the nucleic acid in the sample, as described above, or a consensus of sequences of an amplification product of such nucleic acid. The comparison can be performed at one or more designated positions of the reference sequence. A subset of sequenced nucleic acids can be identified that includes a position corresponding to a designated position of the reference sequence when the respective sequences are maximally aligned. Within such a subset, it can be determined which, if any, sequenced nucleic acids include a nucleotide variation at a designated position, and optionally which, if any, include a reference nucleotide (e.g., the same as in the reference sequence). If the number of sequenced nucleic acids including a nucleotide variant in the subset exceeds a selected threshold, a variant nucleotide can be called at the designated position. The threshold value can be, among other possibilities, a simple number of sequenced nucleic acids that contain a nucleotide variant in the subset, such as at least 1, 2, 3, 4, 5, 6, 7, 8, 9, or 10, or a ratio of sequenced nucleic acids that contain a nucleotide variant in the subset, such as at least 0.5, 1, 2, 3, 4, 5, 10, 15, or 20. Comparison can be repeated for any designated position of interest in the reference sequence. Sometimes, comparison can be performed for designated positions that occupy at least about 20, 100, 200, or 300 consecutive positions of the reference sequence, such as about 20 to 500, or about 50 to 300 consecutive positions.

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

[0238] The method of the present disclosure can be implemented using or with the aid of a computer system.For example, such a method can include: (a) adding a set of epigenetic control nucleic acid molecules to the nucleic acid molecules in a polynucleotide sample, thereby creating a spike-in sample; (b) partitioning the nucleic acid molecules of the spike-in sample into a plurality of partition sets; (c) enriching a subset of molecules from the plurality of partition sets to generate a plurality of enriched molecules, wherein the plurality of enriched molecules comprises a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from the polynucleotide sample; (d) sequencing the plurality of enriched molecules to generate a plurality of sequencing reads; (e) analyzing the plurality of sequencing reads to generate epigenetic compartment scores of the plurality of epigenetic control nucleic acid molecules; and (f) using a computer processor to perform the plurality of epigenetic compartment scores and a plurality of epigenetic compartment cutoffs. In this embodiment, the system includes components for adding, compartmentalizing, enriching, and sequencing epigenetic control nucleic acid molecules.

[0239] In another embodiment, a system for evaluating a method of partitioning nucleic acid molecules in a sample of polynucleotides based on epigenetic state includes a communications interface that receives, via a communications network, a set of sequencing reads of a spike-in sample generated by a nucleic acid sequencer, the set of sequencing reads including: (i) at least a first population of sequencing reads generated from polynucleotides originating from the sample, where the sequencing reads from the first population comprise tag sequences and sequences derived from the polynucleotides originating from the sample; and (ii) at least a second population of sequencing reads derived from epigenetic control nucleic acid molecules, where the sequencing reads generated from the second population comprise an epigenetic modification region and optionally an identifier region. A system comprising: a communications interface; and a computer in communication with the communications interface, the computer including one or more computer processors; and a computer-readable medium including machine-executable code that, when executed by the one or more computer processors, implements a method comprising: (i) receiving, via a communications network, a set of sequencing reads from a first population and a second population of sequencing reads from a nucleic acid sequencer; (ii) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for an epigenetic control nucleic acid molecule and / or an endogenous control molecule; and (iii) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0240] In another embodiment, a system includes a controller that includes or can access a computer-readable medium including non-transitory computer-executable instructions that, when executed by at least one electronic processor, perform at least the steps of: (a) obtaining a set of sequencing reads of a spike-in sample generated by a nucleic acid sequencer, the spike-in sample comprising sample polynucleotides and epigenetic control nucleic acid molecules, the set of sequencing reads comprising (i) a first population of sequencing reads generated from the sample polynucleotides and (ii) a second population of sequencing reads derived from the epigenetic control nucleic acid molecules; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the epigenetic control nucleic acid molecules and / or endogenous control molecules; and (c) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0241] In another embodiment, a system includes a controller that includes or can access a computer-readable medium including non-transitory computer-executable instructions that, when executed by at least one electronic processor, perform at least the steps of: (a) obtaining a set of sequencing reads of a sample generated by a nucleic acid sequencer, wherein the set of sequencing reads includes sequencing reads generated from polynucleotides of the sample; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for an endogenous control molecule; and (c) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0242] In some embodiments, the system further comprises the step of g) generating a compartment method result state based on the comparison of the epigenetic compartment scores. In some embodiments, the compartment method result state is classified as (i) success if one or more epigenetic compartment scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules fall within a corresponding epigenetic compartment cutoff range, or (ii) failure if at least one of the one or more epigenetic compartment scores of the epigenetic control molecules and / or endogenous control molecules is outside a corresponding epigenetic compartment cutoff range.

[0243] 8 illustrates a computer system 801 programmed or otherwise configured to implement the methods of the present disclosure. The computer system 801 is capable of coordinating various aspects of sample preparation, sequencing, and / or analysis. In some embodiments, the computer system 801 is configured to perform sample preparation and sample analysis, including nucleic acid sequencing.

[0244] The computer system 801 includes a central processing unit (CPU, also referred to herein as "processor" and "computer processor") 805, which may be a single-core processor or a multi-core processor, or multiple processors for parallel processing. The computer system 801 also includes memory or memory locations 810 (e.g., random access memory, read-only memory, flash memory), an electronic storage unit 815 (e.g., hard disk), a communication interface 820 (e.g., network adapter) for communicating with one or more other systems, and peripherals 825 such as cache, other memory, data storage, and / or electronic display adapters. The memory 810, storage unit 815, interface 820, and peripherals 825 communicate with the CPU 805 through a bus (solid lines), such as a communication network or a motherboard. The storage unit 815 may be a data storage unit (or data repository) for data storage. The computer system 801 may be operatively coupled to a computer network 430 using the communication interface 820. The computer network 830 may be the Internet, an Internet and / or an extranet, or an intranet and / or an extranet in communication with the Internet. The computer network 830 may in some cases be a telecommunications and / or data network. The computer network 830 may include one or more computer servers, thereby enabling distributed computing such as cloud computing. The computer network 830 may in some cases implement a peer-to-peer network with the computer system 801, thereby enabling devices coupled to the computer system 801 to operate as clients or servers.

[0245] The CPU 805 is capable of executing a sequence of machine-readable instructions, which may be embodied in a program or software. The instructions may be stored in a memory location, such as the memory 810. Examples of operations performed by the CPU 405 may include fetch, decode, execute, and writeback.

[0246] The storage unit 815 may store files such as drivers, libraries, and saved programs. The storage unit 815 may store programs and recorded sessions created by users, as well as output(s) associated with the programs. The storage unit 815 may store user data, such as user preferences and user programs. The computer system 801 may, in some cases, include one or more additional data storage units external to the computer system 801, such as, for example, those located on a remote server that communicates with the computer system 801 over an intranet or the Internet. Data may be migrated from one location to another, for example, using a communications network or physical data migration (e.g., using a hard drive, thumb drive, or other data storage mechanism).

[0247] The computer system 801 may communicate with one or more remote computer systems through the network 830. For illustrative purposes, the computer system 801 may communicate with a remote computer system of a user (e.g., an operator). Examples of remote computer systems include a personal computer (e.g., a portable PC), a slate or tablet PC (e.g., Apple® iPad®, Samsung® Galaxy Tab), a telephone, a smartphone (e.g., Apple® iPhone®, Android-enabled devices, Blackberry®), or a personal digital assistant. A user may access the computer system 801 through the network 830.

[0248] The methods described herein may be implemented by machine-executable (e.g., computer processor) code stored in an electronic storage location of computer system 801, such as memory 810 or electronic storage unit 815. The machine-executable or machine-readable code may be provided in the form of software. In use, the code may be executed by processor 805. In some cases, the code may be retrieved from storage unit 815 and stored in memory 810 for ready access by processor 805. In some situations, electronic storage unit 815 may be omitted and machine-executable instructions may be stored in memory 810.

[0249] In one aspect, the disclosure provides a non-transitory computer-readable medium comprising computer-executable instructions that, when executed by at least one electronic processor, performs a method comprising: (a) obtaining a set of sequencing reads generated by a nucleic acid sequencer; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for an epigenetic control nucleic acid molecule; and (f) comparing the one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs.

[0250] The code may be precompiled and configured for use on a machine having a processor adapted to execute the code, or it may be compiled during run time. The code may be provided in a programming language that can be selected to allow the code to be executed in a precompiled or compiled-at-once fashion.

[0251] Aspects of the systems and methods described herein, such as the computer system 801, can be embodied in programming. Various aspects of the technology can generally be considered as "products" or "articles of manufacture" in the form of machine (or processor) executable code and / or associated data carried or embodied in some type of machine-readable medium. The machine executable code can be stored in an electronic storage unit, such memory (e.g., read-only memory, random access memory, flash memory) or hard disk. A "storage" type medium can include any or all tangible memory of a computer, processor, etc., or modules associated therewith, such as various semiconductor memories, tape drives, disk drives, etc., that can provide non-transitory storage at any time for software programming.

[0252] All or a portion of the software may be communicated from time to time over the Internet or various other telecommunications networks. Such communication may, for example, enable loading of the software from one computer or processor to another, for example, from a management server or host computer to a computer platform of an application server. Thus, other types of media that may bear software elements include light waves, radio waves, and electromagnetic waves, such as those used across physical interfaces between local devices through wired and optical landline networks, and by various air links. Physical elements that carry such waves, such as wired or wireless links, optical links, etc., may also be considered media bearing software. As used herein, except when limited to non-transitory, terms such as tangible "storage" media, computer or machine "readable media" refer to any medium that participates in providing instructions to a processor for execution.

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

[0254] The computer system 801 may include or be in communication with an electronic display that includes, for example, a user interface (UI) for providing one or more results of a sample analysis. Examples of UIs include, without limitation, graphical user interfaces (GUIs) and web-based user interfaces.

[0255] Additional details regarding computer systems and networks, databases, and computer program products can be found, for example, in Peterson, Computer Networks: A Systems Approach, Morgan Kaufmann, 5th Ed. (2011); Kurose, Computer Networking: A Top-Down Approach, Pearson, 7 th Ed. (2016), Elmasri, Fundamentals of Database Systems, Addison Wesley, 6th Ed. (2010), Coronel, Database Systems: Design, Implementation, & Management, Cengage Learning, 11 th Ed. (2014), Tucker, Programming Languages,McGraw-Hill Science / Engineering / Math, 2nd Ed. (2006), and Rhoton, Cloud Computing Architected: Solution Design Handbook, Recursive Press (2011). V. Application A. Cancer and other diseases

[0256] In some embodiments, the methods and systems disclosed herein can be used to identify customized or targeted therapies for treating a given disease or condition in a patient based on the classification of whether the nucleic acid variant is of somatic or germline origin. Generally, the disease under consideration is a type of cancer. Non-limiting examples of such cancers include biliary tract cancer, bladder cancer, transitional cell carcinoma, urothelial carcinoma, brain cancer, glioma, astrocytoma, breast cancer, metaplastic carcinoma, cervical cancer, cervical squamous cell carcinoma, rectal cancer, colorectal cancer, colon cancer, hereditary nonpolyposis colorectal cancer, colorectal adenocarcinoma, gastrointestinal stromal tumor (GIST), endometrial cancer, endometrial stromal sarcoma, esophageal cancer, esophageal squamous cell carcinoma, esophageal adenocarcinoma, intraocular melanoma, uveal melanoma, gallbladder cancer, gallbladder adenocarcinoma, renal cell carcinoma, renal clear cell carcinoma, transitional cell carcinoma, urothelial carcinoma, Wilms' tumor, leukemia, acute lymphocytic leukemia (ALL), acute myeloid leukemia (AML), chronic lymphocytic leukemia (CLL), chronic myelogenous leukemia (CML), chronic myelomonocytic ... ML), liver cancer, hepatoma, hepatoma, hepatocellular carcinoma, cholangiocarcinoma, hepatoblastoma, lung cancer, non-small cell lung cancer (NSCLC), mesothelioma, B cell lymphoma, non-Hodgkin's lymphoma, diffuse large B cell lymphoma, mantle cell lymphoma, T cell lymphoma, non-Hodgkin's lymphoma, precursor T lymphoblastic lymphoma / leukemia, peripheral T cell lymphoma, multiple myeloma, nasopharyngeal carcinoma (NPC), neuroblastoma, oropharyngeal carcinoma, oral squamous cell carcinoma, osteosarcoma, ovarian cancer, pancreatic cancer, pancreatic ductal adenocarcinoma, pseudopapillary neoplasm, acinar cell carcinoma, prostate cancer, prostate adenocarcinoma, skin cancer, melanoma, malignant melanoma, cutaneous melanoma, small intestine carcinoma, gastric cancer, gastrointestinal stromal tumor (GIST), uterine cancer, or uterine sarcoma.

[0257] Non-limiting examples of other genetically-based diseases, disorders, or conditions that may be optionally evaluated using the methods and systems disclosed herein include achondroplasia, alpha 1 antitrypsin deficiency, antiphospholipid syndrome, autism, autosomal dominant polycystic kidney disease, Charcot-Marie-Tooth disease (CMT), cri-a-cat syndrome, Crohn's disease, cystic fibrosis, Dercum's disease, Down's syndrome, Duane's syndrome, Duchenne muscular dystrophy, factor V Leiden thrombophilia, familial hypercholesterolemia, familial Mediterranean fever, fragile These include syndrome X, Gaucher disease, hemochromatosis, hemophilia, holoprosencephaly, Huntington's disease, Klinefelter syndrome, Marfan syndrome, myotonic dystrophy, neurofibromatosis, Noonan syndrome, osteogenesis imperfecta, Parkinson's disease, phenylketonuria, Poland anomaly, porphyria, progeria, retinitis pigmentosa, severe combined immunodeficiency (SCID), sickle cell disease, spinal muscular atrophy, Tay-Sachs, thalassemia, trimethylaminuria, Turner syndrome, palatocardiofacial syndrome, WAGR syndrome, and Wilson's disease. B. Treatment and Related Procedures

[0258] In certain embodiments, the methods disclosed herein relate to identifying and administering customized treatments to patients, taking into account the status of somatic or germline origin of nucleic acid variants. In some embodiments, essentially any cancer treatment (e.g., surgery, radiation, chemotherapy, and / or the like) can be included as part of these methods. In general, customized treatments include at least one immunotherapy (or immunotherapeutic agent). Immunotherapy generally refers to a method of enhancing immune response against a given cancer type. In certain embodiments, immunotherapy refers to a method of enhancing T cell response against tumors or cancer.

[0259] In certain embodiments, the somatic or germline origin status of the nucleic acid variant from the sample from the subject can be compared with the comparator database obtained from a reference population to identify customized or targeted treatment for the subject.Generally, the reference population comprises patients with the same type of cancer or disease as the test subject, and / or patients who are undergoing or have undergone the same treatment as the test subject.If the nucleic acid variant and the comparator result meet certain classification criteria (e.g., substantially or nearly match), customized or targeted treatment (or multiple treatments) can be identified.

[0260] In certain embodiments, the customized treatment described herein is generally administered parenterally (e.g., intravenously or subcutaneously). Pharmaceutical compositions containing immunotherapeutic agents are generally administered intravenously. Certain therapeutic agents are administered orally. However, the customized treatment (e.g., immunotherapeutic agent, etc.) can also be administered by any method known in the art, including, for example, buccal, sublingual, rectal, vaginal, urethral, ​​topical, intraocular, intranasal, and / or intraauricular, and administration can include tablets, capsules, granules, aqueous suspensions, gels, sprays, suppositories, salves, ointments, etc. EXAMPLES

[0261] Example 1 Assessment of cell-free DNA sample compartments A cell-free DNA sample from a patient is analyzed in this example. A spike-in sample is created by combining a cell-free DNA sample with a set of epigenetic control nucleic acid molecules. In this example, the epigenetic control nucleic acid molecules are double-stranded DNA molecules, and the set of epigenetic control nucleic acid molecules is a pool of six different subsets (subset 1 to subset 6) of epigenetic control nucleic acid molecules. Subset 1, subset 2, subset 3, subset 4, subset 5, and subset 6 include epigenetic control nucleic acid molecules with 0, 1, 3, 5, 7, and 9 methylated cytosines (5-methylcytosines) in the epigenetic modification region. The epigenetic control nucleic acid molecules have a molecular barcode at one end of the epigenetic modification region, and epigenetic state barcodes are present at both ends of the epigenetic modification region. The molecular barcodes used in this example are unique molecular barcodes, i.e., each epigenetic control nucleic acid molecule has a distinct molecular barcode.

[0262] The spike-in sample is then combined with a methyl-binding domain (MBD) buffer and magnetic beads conjugated with MBD protein and incubated overnight. During this incubation, the MBD protein binds to the methylated DNA (if present in the cell-free DNA sample) and the methylated epigenetic control nucleic acid molecule. The unmethylated or less methylated DNA is washed off the beads with a buffer containing increasing concentrations of salt. Finally, a buffer with a high salt concentration is used to wash the heavily methylated DNA off the MBD protein. These washes result in three compartments of DNA with increasing methylation (three compartmentalization sets, low, medium, and high). The compartmentalized DNA present in the compartmentalization sets includes DNA from the cell-free DNA sample and the epigenetic control nucleic acid molecule. The compartmentalized DNA in the three compartmentalization sets is clarified to remove salt and concentrated in preparation for the enzymatic step of library preparation.

[0263] After enrichment of the DNA in the compartmentalized sets, the terminal overhangs of the compartmentalized DNA are extended and an adenosine residue is added to the 3' end of the fragments. The 5' end of each fragment is phosphorylated. These modifications render the compartmentalized DNA ligable. DNA ligase and adapters are added to ligate each compartmentalized DNA molecule to an adapter at each end. These adapters contain non-unique barcodes, and each compartmentalized set is ligated to an adapter with a non-unique barcode that is distinguishable from the barcodes in the adapters used for the other compartmentalized sets. After ligation, the three compartmentalized sets are pooled together and amplified by PCR.

[0264] After PCR, the amplified DNA is again clarified, concentrated, and then enriched. Once enriched, the amplified DNA is combined with a salt buffer and a biotinylated RNA probe targeting a specific region of interest, and the epigenetic control nucleic acid molecule and this mixture are incubated overnight. The biotinylated RNA probe is captured by streptavidin magnetic beads and separated from the uncaptured amplified DNA by a series of salt washes, thereby enriching the sample. After enrichment, a sample index is incorporated into the enriched molecules via PCR amplification. After PCR amplification, the amplified molecules from different samples (within a batch) are pooled together and sequenced using an Illumina NovaSeq sequencer.

[0265] The sequence reads generated by the sequencer are then analyzed using bioinformatics tools / algorithms to generate epigenetic compartment scores for epigenetic control nucleic acid molecules belonging to each subset present in each of the three compartmentalization sets. Figure 9A shows a graphical plot of epigenetic compartment scores for epigenetic control nucleic acid molecules belonging to each of the six subsets described in this example (subset 1, subset 2, subset 3, subset 4, subset 5 and subset 6) in the high compartmentalization set. Figure 9B shows a graphical plot of epigenetic compartment scores for epigenetic control nucleic acid molecules belonging to each of the six subsets in the medium compartmentalization set. Figure 9C shows a graphical plot of epigenetic compartment scores for epigenetic control nucleic acid molecules belonging to each of the six subsets in the low compartmentalization set. As shown in Figure 9, the epigenetic compartment scores for subset 1 in the high compartmentalization set, the medium compartmentalization set and the low compartmentalization set are about 0.1%, 0.3% and 99.6%, respectively. The predetermined epigenetic compartment cutoffs for the subset 1 in the high compartmentalization set, the medium compartmentalization set and the low compartmentalization set are ≦0.3%, ≦0.5% and ≧97%, respectively. Herein, the epigenetic compartment score is expressed as a percentage. The epigenetic compartment score of the subset 1 in each of the compartmentalization sets is compared to the corresponding epigenetic compartment cutoff for the subset 1. That is, the epigenetic compartment score of the subset 1 in the high compartmentalization set (0.1%) is compared to the epigenetic compartment cutoff for the subset 1 in the high compartmentalization set (≦0.3%). The epigenetic compartment score of the subset 1 in the high compartmentalization set (0.1%) falls within the range of the epigenetic compartment cutoff for the subset 1 in the high compartmentalization set (≦0.3%). Similarly, the epigenetic compartment score of the subset 1 in the medium compartmentalization set and the low compartmentalization set is compared to the respective epigenetic compartment cutoff for the subset 1 in the medium compartmentalization set and the low compartmentalization set.Similarly, compare the epigenetic compartment scores of subset 2, subset 3, subset 4, subset 5 and subset 6 with the respective epigenetic compartment cutoffs in all three compartmentalization sets.Therefore, have a total of 18 epigenetic compartment scores (6×3=18 for 6 subsets in three compartmentalization sets), and compare each of these epigenetic compartment scores with the corresponding epigenetic compartment cutoffs.It is found that all 18 epigenetic compartment scores fall within the range of the respective epigenetic compartment cutoffs.Therefore, the compartmentalization method performed on the cell-free DNA sample analyzed in this example is classified as successful. Example 2 Assessment of cell-free DNA sample compartments

[0266] A set of cell-free DNA samples from a set of patients is analyzed in this example. In this example, no epigenetic control nucleic acid molecule is used. Instead, endogenous control molecules in cell-free DNA samples are used to assess the compartments of cell-free DNA samples. Each patient's cell-free DNA sample is combined with a methyl-binding domain (MBD) buffer and magnetic beads conjugated with MBD protein, and incubated overnight. During this incubation, MBD protein binds to methylated DNA (if present in cell-free DNA sample) and methylated epigenetic control nucleic acid molecule. Unmethylated or low methylated DNA is washed off the beads with a buffer containing increasing concentrations of salt. Finally, a high salt concentration buffer is used to wash heavily methylated DNA off the MBD protein. These washes result in three compartments of DNA with increasing methylation (three compartmentalization sets, low, medium, high). The compartmentalized DNA present in the compartmentalization set includes DNA from cell-free DNA sample and epigenetic control nucleic acid molecule. The compartmentalized DNA in the three compartmentalized sets is clarified to remove salts and concentrated in preparation for the enzymatic step of library preparation.

[0267] After enrichment of the DNA in the compartmentalized sets, the terminal overhangs of the compartmentalized DNA are extended and an adenosine residue is added to the 3' end of the fragments. The 5' end of each fragment is phosphorylated. These modifications render the compartmentalized DNA ligable. DNA ligase and adapters are added to ligate each compartmentalized DNA molecule to an adapter at each end. These adapters contain non-unique barcodes, and each compartmentalized set is ligated to an adapter with a non-unique barcode that is distinguishable from the barcodes in the adapters used for the other compartmentalized sets. After ligation, the three compartmentalized sets are pooled together and amplified by PCR.

[0268] After PCR, the amplified DNA is again clarified, concentrated, and then enriched. Once enriched, the amplified DNA is combined with a salt buffer and a biotinylated RNA probe targeting a specific region of interest, and the epigenetic control nucleic acid molecule and this mixture are incubated overnight. The biotinylated RNA probe is captured by streptavidin magnetic beads and separated from the uncaptured amplified DNA by a series of salt washes, thereby adding the sample. After enrichment, a sample index is incorporated into the enriched molecules via PCR amplification. After PCR amplification, the amplified molecules from different samples (within a batch) are pooled together and sequenced using an Illumina NovaSeq sequencer.

[0269] The sequence reads generated by the sequencer are then analyzed using bioinformatics tools / algorithms to generate epigenetic compartment scores for one or more endogenous control molecules. In this example, methyl-half and methyl-5 are used as epigenetic compartment scores. Figure 10A shows a graphical plot of the fraction of highly methylated control molecules for sample 1 in the high compartmentalization set, where sample 1 has a methyl-half score of 11. Figure 10B shows a graphical plot of the fraction of highly methylated molecules for sample 1 in the low compartmentalization set, where sample 1 has a methyl-5 score of 13. Figure 11A shows a graphical plot of the fraction of highly methylated control molecules for sample 2 in the high compartmentalization set, where sample 2 has a methyl-half score of 13. Figure 11B shows a graphical plot of the fraction of highly methylated molecules for sample 2 in the low compartmentalization set, where sample 2's methyl-5 score cannot be determined (as shown in Figure 11B). In this example, the epigenetic compartment cutoffs for methyl-half and methyl-5 are 15 methylated CGs and 20 methylated CGs, respectively. The methyl-half and methyl-5 scores of sample 1 fall within the range of the corresponding epigenetic compartment cutoffs. However, for sample 2, the methyl-half score falls within the range of its corresponding epigenetic compartment cutoffs, but the methyl-5 score does not fall within the range of its corresponding epigenetic compartment cutoffs. Therefore, the compartment method of sample 1 is classified as successful, and the compartment method of sample 2 is classified as unsuccessful.

[0270] Although preferred embodiments of the present invention have been shown and described herein, it will be apparent to those skilled in the art that such embodiments are provided by way of example only. The present invention is not limited by the specific examples provided herein. Although the present invention has been described in conjunction with the above specification, the description and illustration of the embodiments herein are not intended to be construed in a limiting sense. Numerous modifications, changes, and substitutions will readily occur to those skilled in the art without departing from the present invention. Furthermore, it should be understood that all aspects of the present invention are not limited to the specific depictions, configurations, or relative proportions set forth herein, which depend upon a variety of conditions and variables. It should be understood that various alternatives to the embodiments of the present disclosure described herein can be used in the practice of the present invention. It is therefore intended that the present disclosure should embrace any such alternatives, modifications, variations, or equivalents. It is intended that the following claims define the scope of the present invention, and that methods and structures falling within the scope of these claims and their equivalents are covered thereby.

[0271] Although the foregoing disclosure has been described in some detail by way of illustration and example for purposes of clarity and understanding, it will be apparent to those skilled in the art upon reading the disclosure that various changes in form and detail may be made without departing from the true scope of the disclosure and may be practiced within the scope of the appended claims. For example, the features, steps, elements, or other aspects of the methods, systems, computer readable media, and / or components may all be used in various combinations.

[0272] All patents, patent applications, websites, other publications or documents, accession numbers, etc. cited herein are incorporated by reference in their entirety for all purposes as if each item was specifically and individually indicated to be incorporated by reference. Where different versions of sequences are associated with accession numbers at different times, the version associated with the accession number at the effective filing date of this application is intended. By effective filing date, we mean the filing date of the priority application that references the accession number prior to the actual filing date or, where applicable, the filing date of the priority application. Similarly, where different versions of publications, websites, etc. are published at different times, the version published most recently prior to the effective filing date of this application is intended unless otherwise specified. The present invention provides, for example, the following items. (Item 1) A set of epigenetic control nucleic acid molecules comprising two or more subsets of epigenetic control nucleic acid molecules, a subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules comprising an epigenetic modification region; A set of epigenetic control nucleic acid molecules. (Item 2) The set of epigenetic control nucleic acid molecules according to item 1, wherein the epigenetic control nucleic acid molecules further comprise an identifier region. (Item 3) 3. The set of epigenetic control nucleic acid molecules according to item 2, wherein the identifier region is present on one or both sides of the epigenetic modification region of the epigenetic control nucleic acid molecule. (Item 4) A set of epigenetic control nucleic acid molecules described in any of the preceding items, wherein the epigenetic modification region of the epigenetic control nucleic acid molecules in at least one subset contains at least one nucleotide having an epigenetic modification. (Item 5) 5. The set of epigenetic control nucleic acid molecules according to item 4, wherein the subset comprises epigenetic control nucleic acid molecules having the same number of nucleotides with epigenetic modifications. (Item 6) 5. The set of epigenetic control nucleic acid molecules according to item 4, wherein the number of nucleotides having an epigenetic modification in the first subset is different from the number of nucleotides having an epigenetic modification in the second subset. (Item 7) 3. The set of epigenetic control nucleic acid molecules according to item 2, wherein the identifier regions of the epigenetic control nucleic acid molecules comprise molecular barcodes. (Item 8) A set of epigenetic control nucleic acid molecules according to any one of the preceding items, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in the two or more subsets comprise identical nucleic acid sequences. (Item 9) 3. The set of epigenetic control nucleic acid molecules according to item 2, wherein the identifier region further comprises at least one epigenetic state barcode. (Item 10) 3. The set of epigenetic control nucleic acid molecules according to item 2, wherein the identifier region comprises one or more primer binding sites. (Item 11) 3. The set of epigenetic control nucleic acid molecules according to item 1 or 2, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a first subset comprise nucleic acid sequences that are distinguishable from the nucleic acid sequences of the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a second subset. (Item 12) The set of epigenetic control nucleic acid molecules according to any one of the preceding items, wherein the epigenetic modification is DNA methylation. (Item 13) The set of epigenetic control nucleic acid molecules according to any one of the preceding items, wherein the nucleotides having epigenetic modifications include methylated nucleotides. (Item 14) 14. The set of epigenetic control nucleic acid molecules according to item 13, wherein the methylated nucleotides include 5-methylcytosine. (Item 15) 14. The set of epigenetic control nucleic acid molecules according to item 13, wherein the methylated nucleotides include 5-hydroxymethylcytosine. (Item 16) 2. The set of epigenetic control nucleic acid molecules according to item 1, wherein each subset of epigenetic control nucleic acid molecules is present in equimolar concentrations. (Item 17) 2. The set of epigenetic control nucleic acid molecules according to item 1, wherein each subset of epigenetic control nucleic acid molecules is present in a non-equimolar concentration. (Item 18) 14. The set of epigenetic control nucleic acid molecules according to item 13, wherein the number of methylated nucleotides in the epigenetic control nucleic acid molecules in at least one of the subsets is 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, at least 12, at least 15, at least 20, at least 25, at least 30, at least 40 or at least 50. (Item 19) 2. The set of epigenetic control nucleic acid molecules according to item 1, wherein the epigenetic control nucleic acid molecules comprise sequences corresponding to lambda phage DNA, human genomic regions, or a combination of both. (Item 20) (i) a set of epigenetic control nucleic acid molecules comprising two or more subsets of epigenetic control nucleic acid molecules, a set of epigenetic control nucleic acid molecules, wherein a subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules that include an epigenetic modification region; (ii) a set of nucleic acid molecules in a sample of polynucleotides derived from a subject; A population of nucleic acids comprising: (Item 21) 21. The population of nucleic acids of item 20, wherein the epigenetic control nucleic acid molecules further comprise an identifier region. (Item 22) 22. The population of nucleic acids of item 21, wherein the identifier region is present on one or both sides of the epigenetic modification region of the epigenetic control nucleic acid molecule. (Item 23) A population of nucleic acids according to any of the preceding items, wherein the epigenetic modification region of the epigenetic control nucleic acid molecule comprises at least one nucleotide having an epigenetic modification. (Item 24) 24. The population of nucleic acids of claim 23, wherein the subset comprises epigenetic control nucleic acid molecules having the same number of nucleotides with epigenetic modifications. (Item 25) 24. The population of nucleic acids of item 23, wherein the number of nucleotides having an epigenetic modification in the first subset is different from the number of nucleotides having an epigenetic modification in the second subset. (Item 26) 22. The population of nucleic acids of claim 21, wherein the identifier regions of the epigenetic control nucleic acid molecules comprise molecular barcodes. (Item 27) The population of nucleic acids of any one of the preceding paragraphs, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in the two or more subsets comprise identical nucleic acid sequences. (Item 28) 22. The population of nucleic acids of item 21, wherein the identifier region further comprises at least one epigenetic state barcode. (Item 29) 22. The population of nucleic acids of item 21, wherein the identifier region comprises one or more primer binding sites. (Item 30) 22. The population of nucleic acids of item 20 or 21, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a first subset comprise nucleic acid sequences that are distinguishable from the nucleic acid sequences of the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a second subset. (Item 31) The population of nucleic acids according to any one of the preceding items, wherein the epigenetic modification is DNA methylation. (Item 32) The population of nucleic acids according to any one of the preceding items, wherein the nucleotides having an epigenetic modification include methylated nucleotides. (Item 33) 33. The population of nucleic acids according to claim 32, wherein the methylated nucleotides comprise 5-methylcytosine. (Item 34) 33. The population of nucleic acids according to claim 32, wherein the methylated nucleotides comprise 5-hydroxymethylcytosine. (Item 35) 21. The population of nucleic acids of item 20, wherein each subset of epigenetic control nucleic acid molecules is present in equimolar concentrations. (Item 36) 21. The population of nucleic acids of item 20, wherein each subset of epigenetic control nucleic acid molecules is present in a non-equimolar concentration. (Item 37) 33. The population of nucleic acids of claim 32, wherein the number of methylated nucleotides in the epigenetic control nucleic acid molecules in at least one of the subsets is 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, at least 12, at least 15, at least 20, at least 25, at least 30, at least 40, or at least 50. (Item 38) 21. The population of nucleic acids of item 20, wherein the epigenetic control nucleic acid molecules comprise sequences corresponding to lambda phage DNA, human genomic regions, or a combination of both. (Item 39) 1. A method for assessing epigenetic state-based compartments of nucleic acid molecules in a sample of polynucleotides, comprising: a. adding a set of epigenetic control nucleic acid molecules to the nucleic acid molecules in the sample of polynucleotides, thereby creating a spike-in sample; b. partitioning at least a subset of nucleic acid molecules of said spike-in sample into a plurality of compartmentalized sets; c. enriching at least a subset of molecules from said plurality of compartmentalized sets to generate a set of enriched molecules, said enriched set of molecules comprising a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules from said sample of polynucleotides; d. sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; e. analyzing at least a subset of said set of sequencing reads to generate one or more epigenetic compartment scores for said epigenetic control nucleic acid molecules; f. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; The method includes: (Item 40) 1. A method for assessing epigenetic state-based compartments of nucleic acid molecules in a sample of polynucleotides, comprising: a. adding a set of epigenetic control nucleic acid molecules to the nucleic acid molecules in the sample of polynucleotides, thereby creating a spike-in sample; b. partitioning at least a subset of nucleic acid molecules of said spike-in sample into a plurality of compartmentalized sets; c. enriching at least a subset of molecules from said plurality of compartmentalized sets to generate a set of enriched molecules, said enriched set of molecules comprising a group of epigenetic control nucleic acid molecules and a group of nucleic acid molecules derived from said sample of polynucleotides, said group of nucleic acid molecules derived from said sample of polynucleotides comprising a set of endogenous control molecules; d. sequencing at least a subset of the set of enriched molecules to generate a set of sequencing reads; e. analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of epigenetic control nucleic acid molecules and the endogenous control molecules; f. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; The method includes: (Item 41) 1. A method for assessing epigenetic state-based compartments of nucleic acid molecules in a sample of polynucleotides, comprising: a. partitioning nucleic acid molecules from at least a subset of said polynucleotide samples into a plurality of compartmentalized sets; b. enriching at least a subset of molecules from said plurality of compartmentalized sets to generate a set of enriched molecules, said enriched set of molecules comprising a group of nucleic acid molecules from said sample of polynucleotides, said group of nucleic acid molecules from said sample of polynucleotides comprising a set of endogenous control molecules; c. sequencing at least a subset of the enriched set of molecules to generate a plurality of sequencing reads; d. analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of endogenous control molecules; e. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; The method includes: (Item 42) 42. The method of claim 39, 40, or 41, further comprising tagging nucleic acid molecules in a compartmentalization set of the plurality of compartmentalization sets with a set of tags to create a population of tagged nucleic acid molecules, wherein the tagged nucleic acid molecules comprise one or more tags. (Item 43) Item 43. The method of item 42, wherein the set of tags used for a first partitioning set of the plurality of partitioning sets is different from the set of tags used for a second partitioning set of the plurality of partitioning sets. (Item 44) 44. The method of claim 43, wherein the set of tags is attached to the nucleic acid molecule by ligating an adaptor to the nucleic acid molecule, the adaptor comprising one or more tags. (Item 45) g) (i) classifying as success if the one or more epigenetic compartment scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules fall within the range of the corresponding epigenetic compartment cutoff, or (ii) classifying as failure if at least one of the one or more epigenetic compartment scores of the set of epigenetic control molecules and / or endogenous control molecules is outside the range of the corresponding epigenetic compartment cutoff. (Item 46) the set of epigenetic control nucleic acid molecules comprises two or more subsets of epigenetic control nucleic acid molecules; 41. The method of claim 39 or 40, wherein a subset of the two or more subsets of epigenetic control nucleic acid molecules comprises a plurality of epigenetic control nucleic acid molecules that contain epigenetic modification regions. (Item 47) 47. The method of claim 46, wherein the epigenetic control nucleic acid molecule further comprises an identifier region. (Item 48) 48. The method of claim 47, wherein the identifier region is present on one or both sides of the epigenetic modification region of the epigenetic control nucleic acid molecule. (Item 49) The method according to any of the preceding items, wherein the epigenetic modification region of the epigenetic control nucleic acid molecule in at least one subset comprises at least one nucleotide having an epigenetic modification. (Item 50) 50. The method of claim 49, wherein the subset comprises epigenetic control nucleic acid molecules having the same number of nucleotides with epigenetic modifications. (Item 51) 50. The method of claim 49, wherein the number of nucleotides having an epigenetic modification in the first subset is different from the number of nucleotides having an epigenetic modification in the second subset. (Item 52) 48. The method of claim 47, wherein the identifier region of the epigenetic control nucleic acid molecule comprises a molecular barcode. (Item 53) The method of any one of the preceding paragraphs, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in the two or more subsets comprise identical nucleic acid sequences. (Item 54) 48. The method of claim 47, wherein the identifier region further comprises at least one epigenetic state barcode. (Item 55) 48. The method of claim 47, wherein the identifier region comprises one or more primer binding sites. (Item 56) 47. The method of claim 46, wherein the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a first subset comprise nucleic acid sequences that are distinguishable from the nucleic acid sequences of the epigenetic modification regions of the plurality of epigenetic control nucleic acid molecules in a second subset. (Item 57) The method according to any one of the preceding items, wherein the epigenetic modification is DNA methylation. (Item 58) The method according to any one of the preceding claims, wherein the nucleotide having an epigenetic modification comprises a methylated nucleotide. (Item 59) 59. The method of claim 58, wherein the methylated nucleotide comprises 5-methylcytosine. (Item 60) 59. The method of claim 58, wherein the methylated nucleotide comprises 5-hydroxymethylcytosine. (Item 61) 47. The method of claim 46, wherein each subset of epigenetic control nucleic acid molecules is present in equimolar concentrations. (Item 62) 47. The method of claim 46, wherein each subset of epigenetic control nucleic acid molecules is present in a non-equimolar concentration. (Item 63) 59. The method of claim 58, wherein the number of methylated nucleotides in the epigenetic control nucleic acid molecules in at least one of the subsets is 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, at least 12, at least 15, at least 20, at least 25, at least 30, at least 40, or at least 50. (Item 64) 47. The method of claim 46, wherein the epigenetic control nucleic acid molecule comprises a sequence corresponding to a genomic region of lambda phage DNA, a human genomic region, or a combination of both. (Item 65) The method according to any of the preceding claims, wherein the epigenetic state is a methylation level of the nucleic acid molecule. (Item 66) 42. The method of claim 39, 40 or 41, wherein the plurality of compartmentalized sets comprises nucleic acid molecules of the spike-in sample partitioned based on the methylation levels of the nucleic acid molecules. (Item 67) 47. The method of claim 46, wherein the epigenetic modification region of the epigenetic control nucleic acid molecule is about 160 bp in length. (Item 68) 42. The method of claim 39, 40 or 41, wherein sequencing of the plurality of enriched molecules is performed by a nucleic acid sequencer. (Item 69) Item 69. The method of item 68, wherein the nucleic acid sequencer is a next-generation sequencer. (Item 70) 47. The method of claim 46, wherein the epigenetic modification region of the epigenetic control nucleic acid molecule comprises a nucleic acid sequence corresponding to a non-human genome. (Item 71) 20. The method of claim 19, wherein the sample of polynucleotides is selected from the group consisting of a sample of DNA, a sample of RNA, a sample of polynucleotides, a sample of cell-free DNA, and a sample of cell-free RNA. (Item 72) 2. The method of any of the preceding claims, wherein the sample of polynucleotides is cell-free DNA. (Item 73) 73. The method of claim 72, wherein the number of methylated nucleotides in the epigenetic control nucleic acid molecules in at least one of the subsets is 0, 2, 4, 6, 8, 10, 12, 14, at least 16, at least 20, at least 30, at least 40, or at least 50. (Item 74) 73. The method of claim 72, wherein the cell-free DNA is between 1 ng and 500 ng. (Item 75) The method of any of the preceding items, wherein the epigenetic control nucleic acid molecule is between 1 femtomole and 200 femtomoles. (Item 76) 20. The method of any preceding claim, wherein the compartmentalization comprises compartmentalizing the nucleic acid molecules based on differential binding affinity of the nucleic acid molecules to a binding agent that preferentially binds to nucleic acid molecules that include nucleotides having epigenetic modifications. (Item 77) 1. A system for evaluating a method for partitioning nucleic acid molecules in a sample of polynucleotides based on their epigenetic state, comprising: a communications interface that receives, via a communications network, a set of sequencing reads for the spike-in sample generated by the nucleic acid sequencer, a communications interface, the set of sequencing reads comprising: (i) at least a first population of sequencing reads generated from polynucleotides originating from the sample, where the sequencing reads from the first population comprise tag sequences and sequences derived from the polynucleotides originating from the sample; and (ii) at least a second population of sequencing reads generated from epigenetic control nucleic acid molecules, where the sequencing reads generated from the second population comprise an epigenetic modification region and optionally an identifier region; a computer in communication with the communications interface, the computer comprising: one or more computer processors; and when executed by the one or more computer processors, (i) receiving, via the communications network, a set of sequencing reads from the first and second populations of sequencing reads from the nucleic acid sequencer; (ii) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the epigenetic control nucleic acid molecules and / or endogenous control molecules; (iii) comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; A computer readable medium including machine executable code for implementing a method including Including, the system. (Item 78) When executed by at least one electronic processor, the method includes at least (a) obtaining a set of sequencing reads of a spike-in sample generated by a nucleic acid sequencer, the spike-in sample comprising polynucleotides of a sample and epigenetic control nucleic acid molecules, the set of sequencing reads comprising (i) a first population of sequencing reads generated from the polynucleotides of the sample and (ii) a second population of sequencing reads generated from the epigenetic control nucleic acid molecules; (b) analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the epigenetic control nucleic acid molecules and / or endogenous control molecules; (c) comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; A system including a controller that includes or has access to a computer-readable medium that includes non-transitory computer-executable instructions that implement the method of the present invention. (Item 79) When executed by at least one electronic processor, the method includes at least a. obtaining a set of sequencing reads of a sample generated by a nucleic acid sequencer, the set of sequencing reads including sequencing reads generated from polynucleotides of the sample; b. analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the endogenous control molecule; c. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; A system including a controller that includes or has access to a computer-readable medium that includes non-transitory computer-executable instructions that implement the method of the present invention. (Item 80) g) generating a result state of the compartment method based on the comparison of the epigenetic compartment scores. (Item 81) 5. The system of claim 4, wherein the result status of the partitioning method is classified as (i) success if the one or more epigenetic partition scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules fall within the range of the corresponding epigenetic partition cutoff, or (ii) failure if at least one of the one or more epigenetic partition scores of the epigenetic control molecules and / or endogenous control molecules is outside the range of the corresponding epigenetic partition cutoff. (Item 82) The system of any one of the preceding claims, wherein the epigenetic compartment score comprises a fraction or percentage of the number of hypermethylated epigenetic control nucleic acid molecules and / or hypermethylated control molecules within the compartmentalized set. (Item 83) 81. The system of any one of items 78 to 80, wherein the epigenetic compartment score comprises a fraction or percentage of the number of hypomethylated epigenetic control nucleic acid molecules and / or hypomethylated control molecules within a compartmentalized set. (Item 84) 84. The system of claim 82 or 83, wherein the compartmentalization set is a hypermethylated compartmentalization set. (Item 85) 84. The system of claim 82 or 83, wherein the compartmentalization set is a hypomethylated compartmentalization set. (Item 86) 81. The system of any one of items 78 to 80, wherein the epigenetic compartment score is a 0 CG score. (Item 87) 81. The system of any one of items 78 to 80, wherein the epigenetic compartment score is a hypo score. (Item 88) 81. The system of any one of items 78 to 80, wherein the epigenetic compartment score is a methyl-half. (Item 89) 81. The system of any one of items 78 to 80, wherein the epigenetic compartment score is methyl-5. (Item 90) 87. The system of item 86, wherein the epigenetic compartment cutoff for the 0 CG score is 0.01%, 0.02%, 0.05%, 0.1%, 0.2%, 0.3%, 0.4%, 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1%, 2%, 5%, at least 5% or at least 10%. (Item 91) 88. The system of claim 87, wherein the epigenetic compartment cutoff for the hypo score is 0.1%, 0.5%, 1%, 2%, 3%, 4%, 5%, 7% or at least 10%. (Item 92) 89. The system of item 88, wherein the epigenetic compartment cutoff for the methyl-half is 5, 10, 15, 20, 25, 30, 35, or 40 mCG. (Item 93) 90. The system of claim 89, wherein the epigenetic compartment cutoff for methyl-5 is 5, 10, 20, 30, 40 or 50 mCG. (Item 94) 94. The method or system of any one of items 39 to 93, further comprising the step of generating a report optionally comprising information about the compartment of the nucleic acid molecule and / or information derived from the compartment of the nucleic acid molecule. (Item 95) Item 95. The method or system of item 94, further comprising the step of transmitting the report to the subject from whom the sample was derived or to a third party, such as a health care practitioner.

Claims

1. 1. A method for evaluating a compartment of nucleic acid molecules in a sample of polynucleotides based on epigenetic state, said epigenetic state being a level or degree of epigenetic modification, said method comprising: a. partitioning nucleic acid molecules from at least a subset of said polynucleotide samples into a plurality of partition sets; b. enriching at least a subset of molecules from the plurality of compartmentalized sets that comprise a particular target genomic region to generate a set of enriched molecules, the set of enriched molecules comprising a group of nucleic acid molecules from the sample of polynucleotides, the group of nucleic acid molecules from the sample of polynucleotides comprising a set of endogenous control molecules corresponding to at least one human genomic region whose epigenetic state is unchanged; c. sequencing at least a subset of the enriched set of molecules to generate a plurality of sequencing reads; d. analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the set of endogenous control molecules; e. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; f. (i) classifying the method as successful if the one or more epigenetic section scores of the set of epigenetic control nucleic acid molecules and / or the endogenous control molecules fall within a corresponding epigenetic section cutoff range, or (ii) classifying the method as unsuccessful if at least one of the one or more epigenetic section scores of the set of epigenetic control molecules and / or the endogenous control molecules is outside a corresponding epigenetic section cutoff range. A method comprising:

2. adding a set of epigenetic control nucleic acid molecules to the nucleic acid molecules in the sample of polynucleotides, thereby creating a spike-in sample; Step a. is performed on the spike-in sample; In step b., the enriched set of molecules further comprises a group of epigenetic control nucleic acid molecules; 2. The method of claim 1, wherein step d. further generates one or more epigenetic compartment scores for said epigenetic control nucleic acid molecules.

3. i) the set of epigenetic control nucleic acid molecules comprises two or more subsets of epigenetic control nucleic acid molecules, one subset of the two or more subsets of epigenetic control nucleic acid molecules comprising a plurality of epigenetic control nucleic acid molecules comprising an epigenetic modification region; ii) the epigenetic control nucleic acid molecule further comprises an identifier region; iii) the epigenetically modified region of the epigenetic control nucleic acid molecule in at least one subset comprises at least one nucleotide having an epigenetic modification; and / or iv) the plurality of compartmentalized sets comprises nucleic acid molecules of the spike-in sample compartmentalized based on the methylation levels of the nucleic acid molecules; The method of claim 2.

4. 4. The method of claim 1, wherein the compartmentalization comprises compartmentalizing the nucleic acid molecules based on their differential binding affinity to a binding agent that preferentially binds to nucleic acid molecules that include nucleotides having epigenetic modifications.

5. 5. The method of any one of claims 1-4, further comprising tagging nucleic acid molecules in a compartmentalized set of the plurality of compartmentalized sets with a set of tags to create a population of tagged nucleic acid molecules, wherein the tagged nucleic acid molecules comprise one or more tags.

6. i) the epigenetic modification is DNA methylation; ii) the nucleotide having the epigenetic modification comprises a methylated nucleotide; and / or iii) The method according to any one of claims 1 to 5, wherein the epigenetic state is the methylation level of the nucleic acid molecule.

7. The epigenetic compartment score is i) the number of hypermethylated epigenetic control nucleic acid molecules and / or hypermethylated control molecules in the compartmentalized set, or ii) the number of hypomethylated epigenetic control nucleic acid molecules and / or hypomethylated control molecules in the compartmentalized set and / or The partitioned set comprises: a) a set of hypermethylated compartments, or b) Hypomethylated Compartment Set The method according to any one of claims 1 to 6, wherein

8. The epigenetic compartment score is i) 0 CG score; ii) hypo score; iii) methyl-half, or iv) methyl-5. The method according to any one of claims 1 to 7,

9. 9. The method of any one of claims 1 to 8, wherein the sequencing of the plurality of enriched molecules is performed by a nucleic acid sequencer.

10. The sample of polynucleotides comprises: i) selected from the group consisting of a sample of DNA, a sample of RNA, a sample of polynucleotides, a sample of cell-free DNA, and a sample of cell-free RNA; and / or ii) The method according to any one of claims 1 to 9, wherein the DNA is cell-free DNA.

11. The method of claim 1, further comprising the step of generating a report.

12. The system, when executed by at least one electronic processor, comprises at least a. Obtaining a set of sequencing reads of a sample generated by a nucleic acid sequencer, the set of sequencing reads comprising sequencing reads generated from polynucleotides of the sample, the sample being compartmentalized and enriched prior to sequencing according to steps a and b of claim 1; b. analyzing at least a subset of the set of sequencing reads to generate one or more epigenetic compartment scores for the endogenous control molecule; c. comparing said one or more epigenetic compartment scores to one or more epigenetic compartment cutoffs; d. generating a compartmentalization method result state based on the comparison of the epigenetic compartment scores, the compartmentalization method result state being classified as (i) success if the one or more epigenetic compartment scores of the set of endogenous control molecules fall within a range of a corresponding epigenetic compartment cutoff, or (ii) failure if at least one of the one or more epigenetic compartment scores of the endogenous control molecules is outside the range of the corresponding epigenetic compartment cutoff. A system including a controller that includes or has access to a computer-readable medium that includes non-transitory computer-executable instructions that implement the method of the present invention.

13. The sample is a spike-in sample and is compartmentalized and enriched according to steps a and b of claim 2; Step b further generates one or more epigenetic compartment scores for said epigenetic control nucleic acid molecules; and The system of claim 12, wherein the result status of step d is classified as (i) success if the one or more epigenetic compartment scores of the set of epigenetic control nucleic acid molecules and / or endogenous control molecules fall within a corresponding epigenetic compartment cutoff range, or (ii) failure if at least one of the one or more epigenetic compartment scores of the epigenetic control molecules and / or endogenous control molecules is outside a corresponding epigenetic compartment cutoff range.

14. The epigenetic compartment score is i) 0 CG score; ii) hypo score; iii) methyl-half, or iv) methyl-5. The system according to claim 12 or 13.

Images