Cell-free DNA end characteristics

Abstract

To measure quantities (e.g., relative frequencies) of sequence end motifs of cell-free DNA fragments in a biological sample of an organism.SOLUTION: Different tissue types exhibit different patterns for the relative frequencies of the sequence end motifs. The present disclosure describes techniques for measuring quantities (e.g., relative frequencies) of sequence end motifs of cell-free DNA fragments in a biological sample of an organism for measuring a property of the sample (e.g., fractional concentration of clinically relevant DNA) and / or determining a condition of the organism based on such measurement.SELECTED DRAWING: Figure 21

Description

[Background technology]

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS This application is a PCT application of and claims the benefit of U.S. Provisional Patent Application No. 62 / 782,316, entitled "CELL-FREE DNA END CHARACTERISTICS," filed December 19, 2018, which is incorporated by reference in its entirety for all purposes.

[0002] Plasma DNA is believed to be composed of cell-free DNA released from multiple tissues in the body, including but not limited to hematopoietic tissues, brain, liver, lung, colon, pancreas, etc. (Sun et al, Proc Natl Acad Sci USA. 2015; 112: E5503-12; Lehmann-Werman et al, Proc Natl Acad Sci USA. 2016; 113: E1826-34; Moss et al, Nat Commun. 2018; 9: 5068). It has been demonstrated that plasma DNA molecules (a type of cell-free DNA molecule) are generated through non-random processes, e.g., their size profile shows a 10 bp periodicity occurring with a major peak of 166 bp and a minor peak (Lo et al, Sci Transl Med. 2010; 2: 61ra91; Jiang et al, Proc Natl Acad Sci USA. 2015; 112: E1317-25).

[0003] Very recently, it was reported that a subset of human genomic locations (e.g., locations on a reference genome) are preferentially cleaved, thereby generating plasma DNA fragments with terminal locations that have a relationship to the tissue of origin (Chan et al, Proc Natl Acad Sci USA. 2016;113:E8159-8168; Jiang et al, Proc Natl Acad Sci USA. 2018;doi:10.1073 / pnas.1814616115). Chandrananda et al. (BMC Med Genomics. 2015;8:29) used the de novo discovery software DREME (Bailey, Bioinformatics. 2011;27:1653-9) to mine cell-free DNA data for motifs associated with nuclease cleavage, regardless of tissue type. Summary of the Invention

[0004] The present disclosure describes a technique for measuring the amount (e.g., relative frequency) of sequence end motifs of cell-free DNA fragments in biological samples of an organism, to measure the characteristics of the sample (e.g., the fractional concentration of clinically relevant DNA) and / or to determine the state of the organism based on such measurements. Different tissue types show different patterns of the relative frequency of sequence end motifs. The present disclosure provides various uses for measuring the relative frequency of sequence end motifs of cell-free DNA, for example, in a mixture of cell-free DNA from various tissues. The DNA from one of such tissues can be referred to as clinically relevant DNA.

[0005] Various examples may quantify the amount of sequence motifs (terminal motifs) that represent the terminal sequences of DNA fragments. For example, embodiments may determine the relative frequency of a set of sequence motifs for the terminal sequences of DNA fragments. In various implementations, a preferred set of terminal motifs and / or a pattern of terminal motifs may be determined using genotypic (e.g., tissue-specific alleles) or phenotypic approaches (e.g., using samples with the same condition). The relative frequency of having a preferred set or a particular pattern may be used to classify properties of new samples (e.g., fractional concentration of clinically relevant DNA) or to measure the state of an organism (e.g., fetal gestational age or level of pathology). Thus, embodiments may provide measurements to inform physiological changes, including cancer, autoimmune disease, transplantation, and pregnancy.

[0006] As a further example, sequence end motifs can be used for physical and / or in silico enrichment of biological samples for clinically relevant cell-free DNA fragments. Enrichment can use sequence end motifs that are preferred for clinically relevant tissues such as fetuses, tumors, or transplants. Physical enrichment can use one or more probe molecules that detect a particular set of sequence end motifs, such that the biological sample is enriched for clinically relevant DNA fragments. For in silico enrichment, a group of sequence reads of cell-free DNA fragments that have one of a set of preferred end sequences for clinically relevant DNA can be identified. Particular sequence reads can be stored based on the likelihood that they correspond to clinically relevant DNA, the likelihood describing the sequence reads that contain the preferred sequence end motif. The stored sequence reads can be analyzed to determine the characteristics of the clinically relevant DNA biological sample.

[0007] These and other embodiments of the disclosure are described in detail below. For example, other embodiments relate to systems, devices, and computer-readable media associated with the methods described herein.

[0008] A better understanding of the nature and advantages of the embodiments of the present disclosure can be obtained with reference to the following detailed description and the accompanying drawings. [Brief description of the drawings]

[0009] [Figure 1] 1 shows examples of terminal motifs according to embodiments of the present disclosure.

[0010] [Diagram 2] FIG. 1 shows a schematic diagram of a genotypic difference-based approach for analyzing differential terminal motif patterns between fetal and maternal DNA molecules according to an embodiment of the present disclosure.

[0011] [Diagram 3] 1 shows a bar graph of terminal motif frequency between fetal and maternal DNA molecules according to an embodiment of the present disclosure.

[0012] [Figure 4] 1 shows the top 10 terminal motifs from FIG. 3 for fetal and shared (i.e., fetal plus maternal) sequences according to an embodiment of the present disclosure.

[0013] [Figure 5A] 1 shows a box plot of the entropy between fetal and maternal DNA molecules in pregnant women according to an embodiment of the present invention. [Figure 5B] 1 shows a box plot of the entropy between fetal and maternal DNA molecules in pregnant women according to an embodiment of the present invention.

[0014] [Figure 6A] 1 shows a hierarchical clustering analysis for fetal and maternal DNA molecules according to an embodiment of the present disclosure. [Figure 6B] 1 shows a hierarchical clustering analysis for fetal and maternal DNA molecules according to an embodiment of the present disclosure.

[0015] [Figure 7]7A and 7B show entropy distributions using all motifs for pregnant women across different trimesters, according to an embodiment of the present disclosure, and 7C and 7D show entropy distributions using 10 motifs for pregnant women across different trimesters, according to an embodiment of the present disclosure.

[0016] [Figure 8A] Figure 1 shows the entropy for all fragments across different gestational ages. It was shown that the entropy of plasma DNA fragments in subjects in the third trimester was lower than those in the first and second trimesters (p-value=0.06). [Figure 8B] Figure 1 shows the entropy for Y chromosome derived fragments across different gestational ages. It was shown that the entropy of Y chromosome derived fragments in subjects in the third trimester was lower than those in the first and second trimesters (p-value=0.01).

[0017] [Figure 9] 1 shows the distribution of the top 10 ranked terminal motifs among fetal and maternal DNA molecules across different trimesters, according to an embodiment of the present disclosure. [Figure 10] 1 shows the distribution of the top 10 ranked terminal motifs among fetal and maternal DNA molecules across different trimesters, according to an embodiment of the present disclosure.

[0018] [Figure 11] 1 shows the combined frequency of the top 10 ranked motifs among fetuses and shared molecules across different trimesters according to an embodiment of the present disclosure.

[0019] [Figure 12] FIG. 1 shows a schematic diagram of a genotypic difference-based approach for analyzing differential terminal motif patterns between mutant and shared molecules in plasma DNA of cancer patients according to an embodiment of the present disclosure.

[0020] [Figure 13]1 shows the profile of cancer-associated variants and shared molecular plasma DNA terminal motifs in hepatocellular carcinoma according to embodiments of the present disclosure.

[0021] [Figure 14] 1 shows the radial landscape of cancer-associated variants and shared molecular plasma DNA terminal motifs in hepatocellular carcinoma according to an embodiment of the present disclosure.

[0022] [Figure 15A] 1 shows the top 10 terminal motifs in rank order difference of terminal motif frequency between variant and shared sequences in plasma DNA of HCC patients according to an embodiment of the present disclosure.

[0023] [Figure 15B] 1 shows the combined frequency of eight terminal motifs for HCC patients and pregnant women according to an embodiment of the present disclosure.

[0024] [Figure 16A] 1 shows entropy values ​​for shared and variant fragments for terminal motifs of different sets for HCC cases according to an embodiment of the present disclosure. [Figure 16B] 1 shows entropy values ​​for shared and variant fragments for terminal motifs of different sets for HCC cases according to an embodiment of the present disclosure.

[0025] [Figure 17] A plot of motif diversity score (entropy) versus measured circulating tumor DNA fraction, according to an embodiment of the present disclosure.

[0026] [Figure 18A] 1 shows an entropy analysis using donor-specific fragments according to an embodiment of the present disclosure. [Figure 18B] Hierarchical clustering analysis using donor-specific fragments is shown.

[0027] [Figure 19]1 is a flow chart showing a method for estimating the fractional concentration of clinically relevant DNA in a biological sample of a subject, according to an embodiment of the present disclosure.

[0028] [Figure 20] 1 is a flow chart illustrating a method for determining the gestational age of a fetus by analyzing a biological sample from a female subject carrying the fetus, according to an embodiment of the present disclosure.

[0029] [Figure 21] FIG. 1 shows a schematic diagram of a phenotypic approach for plasma DNA terminal motif analysis according to an embodiment of the present disclosure.

[0030] [Figure 22] 1 shows an example of a frequency profile of 4-mer terminal motifs between HCC and HBV subjects using total plasma DNA molecules, according to an embodiment of the present disclosure.

[0031] [Figure 23A] 1 shows box plots of the combined frequencies of the top 10 plasma DNA 4-mer terminal motifs for various subjects with different levels of cancer, where the levels are: Control: healthy control subjects, HBV: chronic hepatitis B carriers, Cirr: subjects with cirrhosis, eHCC: early-stage HCC, iHCC: immediate-stage HCC, and aHCC: advanced-stage HCC, according to an embodiment of the present disclosure. [Figure 23B] 1 shows a receiver operating characteristic (ROC) curve of the combined frequency of the top 10 plasma DNA 4-mer end motifs between HCC and non-cancer subjects according to an embodiment of the present disclosure.

[0032] [Figure 24A] 1 shows a box plot of the frequency of CCA motifs across different groups, according to an embodiment of the present disclosure. [Figure 24B] 1 shows the ROC curve between non-HCC and HCC groups using the most frequent 3mer motifs (CCAs) present in non-HCC subjects according to an embodiment of the present disclosure.

[0033] [Figure 25A] 1 shows a box plot of entropy values ​​across different groups using 256 4-mer terminal motifs, according to an embodiment of the present disclosure. [Figure 25B] 1 shows a box plot of entropy values ​​across different groups using ten 4-mer terminal motifs, according to an embodiment of the present disclosure.

[0034] [Figure 26A] According to an embodiment of the present disclosure, the entropy values ​​using 3mer motifs across various groups are shown in box plots. The entropy of HCC subjects using 3mer motifs (64 motifs in total) was found to be significantly higher (p-value<0.0001) than that of non-HCC subjects. [Figure 26B] 1 shows the ROC curve using the entropy of 64 3mer motifs between non-HCC and HCC groups according to an embodiment of the present disclosure. The AUC was found to be 0.872.

[0035] [Figure 27A] 1 shows a box plot of motif diversity (entropy) scores using 4mers across different groups according to an embodiment of the present disclosure. [Figure 27B] 1 shows a box plot of motif diversity (entropy) scores using 4mers across different groups according to an embodiment of the present disclosure.

[0036] [Figure 28] 13 shows receiver operating curves for various techniques for distinguishing healthy controls from cancer, according to embodiments of the present disclosure.

[0037] [Figure 29] 1 shows receiver operating curves for MDS analysis using various k-mers according to an embodiment of the present disclosure.

[0038] [Diagram 30]1 shows the performance of MDS-based cancer detection for various tumor DNA fractions according to embodiments of the present disclosure.

[0039] [Diagram 31] 13 shows receiver operating curves for MDS, SVM, and logistic regression analyses according to an embodiment of the present disclosure.

[0040] [Diagram 32] 1 shows a hierarchical clustering analysis of the top 10 ranked terminal motifs across different groups with different levels of cancer, including: Control: healthy control subjects, HBV: chronic hepatitis B carriers, Cirr: subjects with cirrhosis, eHCC: early stage HCC, iHCC: immediate stage HCC, and aHCC: advanced stage HCC, according to an embodiment of the present disclosure.

[0041] [Diagram 33] 33A-33C show hierarchical clustering analysis using all plasma DNA molecules across different groups with different levels of cancer, according to an embodiment of the present disclosure.

[0042] [Diagram 34] 1 shows a hierarchical clustering analysis based on 3mer motifs using all plasma DNA molecules across different groups with different levels of cancer, according to an embodiment of the present disclosure.

[0043] [Figure 35A] 1 shows an entropy analysis using all plasma DNA molecules between healthy control subjects and SLE patients according to an embodiment of the present disclosure. [Figure 35B] 1 shows a hierarchical clustering analysis using all plasma DNA molecules between healthy control subjects and SLE patients according to an embodiment of the present disclosure.

[0044] [Diagram 36]1 shows an entropy analysis using plasma DNA molecules with 10 selected terminal motifs between healthy control subjects and SLE patients according to an embodiment of the present disclosure.

[0045] [Figure 37] 1 shows the ROC curves of a combined analysis including terminal motifs and copy number or methylation according to an embodiment of the present disclosure.

[0046] [Figure 38A] 1 shows an entropy analysis based on 4mers jointly constructed from the ends of sequenced plasma DNA fragments and their adjacent genomic sequences in HCC and non-HCC subjects according to an embodiment of the present disclosure. [Figure 38B] 1 shows a 4mer-based clustering analysis jointly constructed from the ends of sequenced plasma DNA fragments and their adjacent genomic sequences in HCC and non-HCC subjects according to an embodiment of the present disclosure.

[0047] [Figure 39] 2 shows a ROC comparison for the techniques 140 and 160 of FIG. 1 used to define terminal motifs in plasma DNA, according to an embodiment of the present disclosure.

[0048] [Diagram 40] 1 shows an accuracy comparison showing that tissue-specific open chromatin regions improve the discrimination power of plasma DNA end motifs, according to embodiments of the present disclosure.

[0049] [Diagram 41] 1 shows size-band-based plasma DNA terminal motif analysis according to an embodiment of the present disclosure.

[0050] [Diagram 42] 1 is a flow chart illustrating a method for classifying a level of pathology in a biological sample of a subject, according to an embodiment of the present disclosure.

[0051] [Diagram 43]1 is a flow chart illustrating a method of enriching a biological sample for clinically relevant DNA according to an embodiment of the present disclosure.

[0052] [Diagram 44] 37 is a flow chart illustrating a method 3700 of enriching a biological sample for clinically relevant DNA according to an embodiment of the present disclosure.

[0053] [Diagram 45] 1 shows an exemplary plot illustrating the increase in fetal DNA fraction using CCCA terminal motifs, according to embodiments of the present disclosure.

[0054] [Diagram 46] 1 illustrates a measurement system according to an embodiment of the present invention.

[0055] [Figure 47] FIG. 1 illustrates a block diagram of an exemplary computer system that can be used with systems and methods according to embodiments of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0056] term A "tissue" corresponds to a group of cells that group together as a functional unit. Two or more types of cells can be found in a single tissue. Different types of tissues can consist of different types of cells (e.g., liver cells, alveolar cells, or blood cells), but also tissues from different organisms (maternal vs. fetal) or healthy vs. tumor cells. A "reference tissue" can correspond to the tissue used to determine tissue-specific methylation levels. Multiple samples of the same tissue type from different individuals can be used to determine tissue-specific methylation levels of that tissue type.

[0057] A "biological sample" refers to any sample taken from a subject (e.g., a human (or other animal), such as a pregnant woman, a person with or suspected of having cancer, an organ transplant recipient, or a subject suspected of having a disease process involving an organ (e.g., the heart in myocardial infarction, the brain in stroke, or the hematopoietic system in anemia)) and containing one or more nucleic acid molecules of interest. A biological sample can be a body fluid, such as blood, plasma, serum, urine, vaginal fluid, fluid from edema (e.g., of the testes), vaginal washing fluid, pleural fluid, ascites, cerebrospinal fluid, saliva, sweat, tears, sputum, bronchoalveolar lavage fluid, nipple discharge, aspirates from various parts of the body (e.g., thyroid, mammary gland), intraocular fluid (e.g., aqueous humor), etc. Stool samples can also be used. In various embodiments, the majority of DNA in a biological sample enriched for cell-free DNA (e.g., a plasma sample obtained via a centrifugation protocol) may be cell-free, e.g., more than 50%, more than 60%, more than 70%, more than 80%, more than 90%, more than 95%, or more than 99% of the DNA may be cell-free. The centrifugation protocol may include, for example, obtaining a fluid portion at 3,000g x 10 minutes and recentrifuging at 30,000g for another 10 minutes to remove residual cells. As part of the analysis of a biological sample, at least 1,000 cell-free DNA molecules may be analyzed. As another example, at least 10,000 or 50,000 or 100,000 or 500,000 or 1,000,000 or 5,000,000 or more cell-free DNA molecules may be analyzed.

[0058] "Clinically relevant DNA" may refer to the DNA of a particular tissue source that is to be measured, for example, to determine the fractional concentration of such DNA or to classify the phenotype of a sample (e.g., plasma). An example of clinically relevant DNA is fetal DNA in maternal plasma, or tumor DNA in a patient's plasma, or other samples that contain cell-free DNA. Another example includes measuring the amount of graft-related DNA in the plasma, serum, or urine of a transplant patient. Further examples include measuring the fractional concentration of hematopoietic and non-hematopoietic DNA in a subject's plasma, or the fractional concentration of liver DNA fragments (or other tissues) in a sample, or the fractional concentration of brain DNA fragments in cerebrospinal fluid.

[0059] A "sequence read" refers to a chain of nucleotides sequenced from any part or all of a nucleic acid molecule. For example, a sequence read can be a short chain of nucleotides (e.g., about 20-150 nucleotides) sequenced from a nucleic acid fragment, a short chain of nucleotides at one or both ends of a nucleic acid fragment, or the sequencing of an entire nucleic acid fragment present in a biological sample. Sequence reads can be obtained, for example, using sequencing techniques, or in a variety of ways using probes, for example, with hybridization arrays or capture probes, or with amplification techniques such as polymerase chain reaction (PCR) or linear amplification, using single primers or isothermal amplification. At least 1,000 sequence reads can be analyzed as part of the analysis of a biological sample. As other examples, at least 10,000 or 50,000 or 100,000 or 500,000 or 1,000,000 or 5,000,000 or more sequence reads can be analyzed.

[0060] A sequence read may include "terminal sequences" associated with the ends of a fragment. An terminal sequence may correspond to the outermost N bases of a fragment, e.g., 2-30 bases at the end of a fragment. If a sequence read corresponds to an entire fragment, the sequence read may include two terminal sequences. If paired terminal sequencing provides two sequence reads corresponding to the ends of a fragment, each sequence read may include one terminal sequence.

[0061] A "sequence motif" may refer to a short repeating pattern of bases in a DNA fragment (e.g., cell-free DNA fragment). A sequence motif may occur at the end of a fragment and thus may be part of or include a terminal sequence. A "terminal motif" may refer to a sequence motif for a terminal sequence that occurs preferentially at the end of a DNA fragment, potentially for a particular type of tissue. A terminal motif may also occur just before or just after the end of a fragment, thereby still corresponding to a terminal sequence.

[0062] The term "allele" refers to alternative DNA sequences at the same physical genomic locus that may or may not result in different phenotypic characteristics. In any given diploid organism (except for the sex chromosomes of male subjects) with two copies of each chromosome, the genotype of each gene includes pairs of alleles present at that locus that are the same in homozygotes and different in heterozygotes. A population or species of organisms typically contains multiple alleles at each locus in various individuals. Genomic loci at which more than one allele is found in a population are called polymorphic sites. Allelic diversity at a locus can be measured as the number of alleles present (i.e., the degree of polymorphism) or the proportion of heterozygotes in a population (i.e., the heterozygosity rate). As used herein, the term "polymorphism" refers to any inter-individual variation in the human genome, regardless of its frequency. Examples of such variation include, but are not limited to, single nucleotide polymorphisms, simple tandem repeat polymorphisms, insertion-deletion polymorphisms, mutations (which may cause disease), and copy number variation. As used herein, the term "haplotype" refers to a combination of alleles at multiple loci that are transmitted together on the same chromosome or chromosomal region. A haplotype can refer to just one pair of loci, or a chromosomal region, or an entire chromosome or chromosomal arm.

[0063] The term "fractional fetal DNA concentration" is used interchangeably with the terms "fetal DNA percentage" and "fetal DNA fraction" and refers to the percentage of fetal DNA molecules present in a biological sample derived from a fetus (e.g., maternal plasma or serum sample) (Lo et al, Am J Hum Genet. 1998; 62: 768-775, Lun et al, Clin Chem. 2008; 54: 1664-1672). Similarly, tumor fraction or tumor DNA fraction can refer to the fractional concentration of tumor DNA in a biological sample.

[0064] "Relative frequency" can refer to a proportion (e.g., a percentage, a fraction, or a concentration). In particular, the relative frequency of a particular terminal motif (e.g., CCGA) can provide the proportion of cell-free DNA fragments that are associated with the terminal motif CCGA, for example, by having a terminal sequence of CCGA.

[0065] "Aggregate value" can refer to, for example, a collective property of the relative frequencies of a set of terminal motifs. Examples include the mean, median, sum of the relative frequencies, the variation between the relative frequencies (e.g., entropy, standard deviation (SD), coefficient of variation (CV), interquartile range (IQR), or a particular percentile cutoff (e.g., 95th or 99th percentile) among the various relative frequencies), or the difference (e.g., distance) of the relative frequencies from a reference pattern, which can be implemented by clustering.

[0066] A "calibration sample" may correspond to a biological sample in which the fractional concentration of clinically relevant DNA (e.g., tissue-specific DNA fraction) is known or determined via a calibration method, for example, using tissue-specific alleles, such as in transplants, where alleles present in the donor's genome but absent in the recipient's genome may be used as markers for the transplanted organ. As another example, a calibration sample may correspond to a sample in which terminal motifs may be determined. A calibration sample may be used for both purposes.

[0067] A "calibration data point" includes a "calibration value" and a measured or known fractional concentration of clinically relevant DNA (e.g., DNA of a particular tissue type). The calibration value may be determined from a relative frequency (e.g., a tabulated value) determined for a calibration sample in which the fractional concentration of clinically relevant DNA is known. The calibration data points may be defined in various ways, for example, as discrete points or as a calibration function (also called a calibration curve or calibration surface). The calibration function may be derived from an additional mathematical transformation of the calibration data points.

[0068] A "site" (also called a "genomic site") corresponds to a single site, which may be a single base position, or a group of correlated base positions, e.g., a CpG site, or a larger group of correlated base positions. A "locus" may correspond to a region that includes multiple sites. A locus may contain only one site that would make the locus equivalent to the site in that context.

[0069] A "methylation index" for each genomic site (e.g., a CpG site) may refer to the proportion of DNA fragments (e.g., as determined from sequence reads or probes) that exhibit methylation at that site over the total number of reads that cover that site. A "read" may correspond to information obtained from a DNA fragment (e.g., the methylation state of the site). A read may be obtained using a reagent (e.g., a primer or a probe) that preferentially hybridizes to DNA fragments of a particular methylation state. Typically, such a reagent is applied after treatment with a process that differentially modifies or differentially recognizes DNA molecules depending on their methylation state, such as bisulfite conversion, or a methylation-sensitive restriction enzyme, or a methylation-binding protein, or an anti-methylcytosine antibody, or a single molecule sequencing technique that recognizes, for example, methylcytosine and hydroxymethylcytosine.

[0070] The "methylation density" of a region may refer to the number of reads at sites in the region that show methylation divided by the total number of reads that cover sites in this region. The sites may have specific characteristics, for example, CpG sites. Thus, the "CpG methylation density" of a region refers to the number of reads that show CpG methylation divided by the total number of reads that cover CpG sites in this region (e.g., a particular CpG site, a CpG site in a CpG island, or a larger region). For example, the methylation density of each 100 kb bin in the human genome may be determined from the total number of unconverted cytosines (corresponding to methylated cytosines) after bisulfite treatment of the CpG sites as a percentage of all CpG sites covered by sequence reads mapped to the 100 kb region. This analysis may also be performed for other bin sizes, such as 500 bp, 5 kb, 10 kb, 50 kb, or 1 Mb. The region may be the whole genome, or a chromosome, or a portion of a chromosome (e.g., a chromosome arm). The methylation index of a CpG site is the same as the methylation density of the region if the region contains only that CpG site. "Percentage of methylated cytosines" may refer to the number of cytosine sites "C" that are shown to be methylated (e.g., unconverted after bisulfite conversion) relative to the total number of cytosine residues analyzed in the region, i.e., including cytosines outside the CpG context. Methylation index, methylation density, and percentage of methylated cytosines are examples of "methylation levels." Apart from bisulfite conversion, other processes known to those skilled in the art can be used to examine the methylation status of DNA molecules, including, but not limited to, methylation-state sensitive enzymes (e.g., methylation-sensitive restriction enzymes), methylation-binding proteins, single molecule sequencing using methylation-state sensitive platforms (e.g., nanopore sequencing (Schreiber et al, Proc Natl Acad Sci USA. 2013; 110: 18910-18915) and Pacific Biosciences single molecule real-time analysis (Flusberg et al, Nat Methods. 2010; 7: 461-465)).The methylation metric of a DNA molecule may correspond to the percentage of sites (e.g., CpG sites) that are methylated. The methylation metric may be specified as an absolute number or a percentage, which may be referred to as the methylation density of the molecule.

[0071] The term "sequencing depth" refers to the number of times a locus is covered by sequence reads aligned to that locus. A locus can be as small as a nucleotide, or as large as a chromosome arm, or as large as an entire genome. Sequencing depth can be expressed as 50x, 100x, etc., where "x" refers to the number of times a locus is covered by sequence reads. Sequencing depth can also be applied to multiple loci or an entire genome, where x can refer to the average number of times a locus or haploid genome or entire genome is sequenced, respectively. Ultra-deep sequencing can refer to a sequencing depth of at least 100x.

[0072] "Separation value" refers to the difference or ratio that encompasses two values, such as two fractional contributions or two methylation levels. Separation value can be a simple difference or ratio. As an example, the direct ratio of x / y is a separation value, as is x / (x+y). Separation value can include other factors, such as multiplicative factors. As another example, the difference or ratio of a function of values ​​can be used, such as the difference or ratio of the natural logarithms (ln) of two values. Separation value can include difference and ratio.

[0073] "Separation values" and "aggregation values" (e.g., relative frequencies) are two examples of parameters (also called metrics) that provide a measure of the variation of a sample between different classifications (states) and thus can be used to determine various classifications. Aggregation values ​​can be separation values, for example, when the difference is taken between a set of relative frequencies of a sample and a reference set of relative frequencies, as is done in clustering.

[0074] The term "classification" as used herein refers to any number(s) or other feature(s) associated with a particular characteristic of a sample. For example, a "+" symbol (or the word "positive") can mean that the sample is classified as having a deletion or amplification. Classification can be binary (e.g., positive or negative) or can have more levels of classification (e.g., a scale of 1 to 10 or 0 to 1).

[0075] The terms "cutoff" and "threshold" refer to a predetermined number used in an operation. For example, a cutoff size may refer to a size above which fragments are excluded. A threshold may be a value above or below which a particular classification is applied. Either of these terms may be used in any of these contexts. A cutoff or threshold may be a "reference value" or may be derived from a reference value that represents a particular classification or that distinguishes between two or more classifications. Such a reference value may be determined in a variety of ways, as will be understood by those skilled in the art. For example, a metric may be determined for two different cohorts of subjects with different known classifications, and a reference value may be selected as a representative of one classification (e.g., the average) or as a value between two clusters of the metric (e.g., selected to obtain a desired sensitivity and specificity). As another example, a reference value may be determined based on statistical simulation of samples.

[0076] The term "level of cancer" may refer to whether cancer is present (i.e., presence or absence), the stage of cancer, the size of the tumor, whether there is metastasis, the total tumor burden in the body, the response of the cancer to treatment, and / or other measures of the severity of the cancer (e.g., recurrence of cancer). The level of cancer may be a number or other indicia, such as a symbol, an alphabetic letter, and a color. The level may be zero. The level of cancer may also include a pre-malignant or pre-cancerous condition (state). The level of cancer may be used in various ways. For example, screening may check whether cancer is present in a person who was not previously known to have cancer. Evaluation may examine a person who has been diagnosed with cancer to monitor the progression of the cancer over time, study the effectiveness of a therapy, or determine a prognosis. In one embodiment, the prognosis may be expressed as the likelihood that the patient will die from the cancer, or the likelihood that the cancer will progress after a certain duration or time, or the likelihood or extent that the cancer will metastasize. Detection may mean "screening" or checking whether a person with suggestive characteristics of cancer (eg, symptoms or other positive test) has cancer.

[0077] "Level of pathology" may refer to the amount, degree, and severity of pathology associated with an organism, and the level may be as described above for cancer. Another example of pathology is the rejection of transplanted organs. Other examples of pathology may include autoimmune attacks (e.g., lupus nephritis or multiple sclerosis, which damage the kidney), inflammatory diseases (e.g., hepatitis), fibrotic processes (e.g., cirrhosis), fatty infiltration (e.g., fatty liver disease), degenerative processes (e.g., Alzheimer's disease), and ischemic tissue damage (e.g., myocardial infarction or stroke). The healthy state of a subject may be considered as a classification without pathology.

[0078] The term "about" or "approximately" may mean within an acceptable error range of a particular value as determined by one of ordinary skill in the art, which depends in part on the limitations of the method of measuring or determining the value, i.e., the measurement system. For example, "about" may mean within 1 or more than 1 standard deviation, as is customary in the art. Alternatively, "about" may mean a range of up to 20%, up to 10%, up to 5%, or up to 1% of a given value. Alternatively, particularly with respect to biological systems or processes, the term "about" or "approximately" may mean within an order of magnitude, within 5-fold, and more preferably within 2-fold of a value. When a particular value is described in the present application and claims, unless otherwise indicated, the term "about" within an acceptable error range of the particular value should be assumed. The term "about" may have the meaning commonly understood by one of ordinary skill in the art. The term "about" may refer to ±10%. The term "about" may refer to ±5%. Detailed Description

[0079] The present disclosure describes a technique for measuring the amount (e.g., relative frequency) of terminal motifs of cell-free DNA fragments in biological samples of organisms to measure the characteristics of the sample and / or to determine the state of the organism based on such measurement. Different tissue types show different patterns of relative frequency of sequence motifs. The present disclosure provides various uses for measuring the relative frequency of cell-free DNA terminal motifs, for example, in a mixture of cell-free DNA from various tissues. The DNA derived from one of such tissues can be called clinically relevant DNA.

[0080] Clinically relevant DNA of a particular tissue (e.g., fetus, tumor, or transplanted organ) exhibits a particular pattern of relative frequency, which can be measured as an aggregate value. Other DNA in the sample may exhibit a different pattern, thereby allowing the amount of clinically relevant DNA in the sample to be measured. Thus, in one example, the fractional concentration (e.g., percentage) of clinically relevant DNA can be determined based on the relative frequency of the terminal motifs. The fractional concentration can be a number, a numerical range, or other classification, e.g., high, medium, or low, or whether the fractional concentration exceeds a threshold. In various implementations, the aggregate value can be the sum of the relative frequencies of a set of terminal motifs, the variance of the relative frequencies of all or a set of terminal motifs (e.g., entropy, also called motif diversity score), or the difference (e.g., total distance) from a reference pattern, e.g., an array (vector) of the relative frequencies of a calibration sample(s) with known fractional concentrations. Such an array can be considered as a reference set of relative frequencies. Such differences can be used in classifiers such as hierarchical clustering, support vector machines, logistic regression, etc. By way of example, clinically relevant DNA can be the DNA of a fetus, a tumor, a transplanted organ, or other tissue (eg, hematopoietic or liver).

[0081] In another example, the level of pathology can be determined using the relative frequency of motifs.Organisms with different phenotypes can show different patterns of motif relative frequency of cell-free DNA fragments.The tabulated value of the relative frequency of terminal motifs can be compared with a reference value to classify phenotype.In various implementations, the tabulated value can be the sum of relative frequency, the variance of relative frequency, or the difference of relative frequency from a reference set.Examples of pathology include cancer and autoimmune diseases such as SLE.

[0082] In another example, the motif relative frequency can be used to determine the gestational age of the fetus. As the gestational age of the fetus increases in the maternal sample, the aggregate value of the relative frequency of the terminal motifs changes. Such aggregate value can be determined as described above and elsewhere.

[0083] Considering that cell-free DNA fragments from a particular tissue have a set of preferred specific terminal motifs, the preferred terminal motifs can be used to enrich the sample for DNA from the particular tissue (clinically relevant DNA). Such enrichment can be performed via physical manipulation to enrich the physical sample. Some embodiments can capture and / or amplify cell-free DNA fragments with terminal sequences that match the set of preferred terminal motifs, for example, using primers or adapters. Other examples are described herein.

[0084] In some embodiments, enrichment may be performed in silico. For example, the system may receive sequence reads and filter the reads based on terminal motifs to obtain a subset of sequence reads that are enriched for corresponding DNA fragments from clinically relevant DNA. If a DNA fragment has a terminal sequence that includes a preferred terminal motif, it may be identified as more likely to originate from the tissue of interest. As described herein, the likelihood may be further determined based on the methylation and size of the DNA fragment.

[0085] The use of such terminal motifs may avoid the need for a reference genome, which may be necessary when using terminal positions (Chan et al, Proc Natl Acad Sci USA. 2016; 113: E8159-8168; Jiang et al, Proc Natl Acad Sci USA. 2018; doi: 10.1073 / pnas.1814616115). Furthermore, because the number of terminal motifs may be less than the number of preferred terminal positions in the reference genome, more statistics may be collected for each terminal motif, improving accuracy.

[0086] For example, Chandrananda et al. found high similarity between maternal and fetal fragments with respect to position-specific nucleotide patterns in terms of mononucleotide frequency in a 51 bp (20 bp upstream / downstream) region surrounding the fragment start site (Chandrananda et al, BMC Med Genomics. 2015;8:29), which is surprising because it means that the use of those methods based on mononucleotide frequency around the termini cannot provide information on the tissue of origin of the cell-free DNA fragments.

[0087] I. Cell-free DNA end motifs Terminal motifs refer to the terminal sequences of cell-free DNA fragments, e.g., sequences of K bases at either end of the fragment. Terminal sequences can be k-mers with various numbers of bases, e.g., 1, 2, 3, 4, 5, 6, 7, etc. Terminal motifs (or "sequence motifs") refer to the sequence itself, as opposed to a specific location in the reference genome. Thus, the same terminal motif can occur at multiple locations throughout the reference genome. Terminal motifs can be determined using the reference genome, for example, to identify the base immediately before the start position or immediately after the end position. Such bases correspond to the ends of cell-free DNA fragments, for example, as they are identified based on the terminal sequences of the fragments.

[0088] FIG. 1 shows examples of terminal motifs according to embodiments of the present disclosure. FIG. 1 shows two methods of defining 4mer terminal motifs to be analyzed. In technique 140, 4mer terminal motifs are constructed directly from the first 4 bp sequence of each end of the plasma DNA molecule. For example, the first 4 nucleotides or the last 4 nucleotides of the sequenced fragment can be used. In technique 160, 4mer terminal motifs are constructed jointly by utilizing a 2mer sequence from the sequenced end of the fragment and another 2mer sequence from the genomic region adjacent to the end of the fragment. In other embodiments, other types of motifs can be used, for example, 1mer, 2mer, 3mer, 5mer, 6mer, and 7mer terminal motifs.

[0089] As shown in Figure 1, cell-free DNA fragments 110 are obtained using a purification process of the blood sample, such as by centrifugation. In addition to plasma DNA fragments, other types of cell-free DNA molecules from, for example, serum, urine, saliva, and other such cell-free samples mentioned herein may be used. In one embodiment, the DNA fragments may be blunt-ended.

[0090] At block 120, the DNA fragment is subjected to paired-end sequencing. In some embodiments, paired-end sequencing may generate two sequence reads from the two ends of the DNA fragment, e.g., 30-120 bases per sequence read. These two sequence reads may form a pair of reads of the DNA fragment (molecule), each sequence read including terminal sequence for each end of the DNA fragment. In other embodiments, the entire DNA fragment may be sequenced, thereby providing a single sequence read including terminal sequence for both ends of the DNA fragment.

[0091] At block 130, the sequence reads may be aligned to a reference genome. This alignment is to illustrate different ways to define sequence motifs and may not be used in some embodiments. The alignment procedure may be performed using various software packages such as BLAST, FASTA, Bowtie, BWA, BFAST, SHRiMP, SSAHA2, NovoAlign, and SOAP.

[0092] The technique 140 shows sequence reads of a sequenced fragment 141 with alignment to the genome 145. Taking the 5' end as the start, a first terminal motif 142 (CCCA) is at the start of the sequenced fragment 141. A second terminal motif 144 (TCGA) is at the tail of the sequenced fragment 141. Such terminal motifs may arise in one embodiment when an enzyme recognizes CCCA and then makes a cut just before the first C. In that case, CCCA is preferentially at the end of the plasma DNA fragment. For TCGA, an enzyme may recognize it and make a cut after the A.

[0093] The technique 160 shows sequence reads of a sequenced fragment 161 with alignment to the genome 165. Taking the 5' end as the start, the first terminal motif 162 (CGCC) has a first portion (CG) that occurs just before the start of the sequenced fragment 161, and a second portion (CC) that is part of the terminal sequence of the start of the sequenced fragment 161. The second terminal motif 164 (CCGA) has a first portion (GA) that occurs just after the tail of the sequenced fragment 161, and a second portion (CC) that is part of the terminal sequence of the tail of the sequenced fragment 161. Such a terminal motif may arise in one embodiment when an enzyme recognizes CGCC and then cleaves between G and C. In that case, CC will preferentially be present at the end of the plasma DNA fragment that is immediately preceded by a CG, thereby providing a terminal motif of CGCC. For the second terminal motif 164 (CCGA), the enzyme may cleave between C and G. In that case, CC will preferentially be present at the end of the plasma DNA fragment. For technique 160, the number of bases from the adjacent genomic regions and the sequenced plasma DNA fragments can be varied and are not necessarily limited to a fixed ratio, e.g., instead of 2:2, the ratio can be 2:3, 3:2, 4:4, 2:4, etc.

[0094] The more nucleotides contained in the cell-free DNA end signature, the more specific the motif, since the probability of having six bases ordered in the correct configuration in a genome is lower than the probability of having two bases ordered in the correct configuration in a genome. Thus, the selection of the length of the end motif may be governed by the sensitivity and / or specificity required for the intended application.

[0095] Since the terminal sequences are used to align sequence reads to a reference genome, any sequence motifs determined from the terminal sequence, or immediately before / after, are still determined from the terminal sequence. Thus, technique 160 creates an association of the terminal sequence to other bases, and the reference is used as a mechanism to create that association. The difference between techniques 140 and 160 is which two terminal motifs a particular DNA fragment is assigned to, which affects the particular values ​​for relative frequency. However, as long as a consistent technique is used for the training data as is used in production, the overall results (e.g., fractional concentration of clinically relevant DNA, classification of level of pathology, etc.) are not affected by how the DNA fragments are assigned to terminal motifs.

[0096] The counted number of DNA fragments having terminal sequences corresponding to a particular terminal motif can be counted (e.g., stored in an array in memory) to determine the relative frequency. As described in more detail below, the relative frequency of terminal motifs for cell-free DNA fragments can be analyzed. Differences in the relative frequency of terminal motifs have been detected for different types of tissues and different phenotypes, such as different levels of pathology. The difference can be quantified by the amount of DNA fragments having a particular terminal motif or the overall pattern across a set of terminal motifs (e.g., all possible combinations of k-mers corresponding to the length used), such as variance (e.g., entropy, also called motif diversity score).

[0097] II. Genotype Difference-Based Approach Different tissue types have been identified as having different terminal motifs. Described herein are methods for using terminal motifs to determine the fractional concentration of clinically relevant DNA, such as fetal DNA, tumor DNA, DNA from transplanted organs, or DNA from a particular organ.

[0098] To identify terminal motifs that are preferential for a particular type of clinically relevant DNA, the genotypic differences can be used to identify DNA fragments as originating from clinically relevant tissues. Once a DNA fragment is detected as originating from a clinically relevant tissue, the terminal motifs of the DNA fragment can be determined. Analysis of the relative frequency of the terminal motifs reveals that the relative frequency of the terminal motifs varies with different tissues. As described below, quantification of the relative frequency differences can be used in combination with a calibration sample(s) in which the fractional concentration of clinically relevant DNA is known (e.g., measured by another technique such as tissue-specific alleles) to determine the classification of the fractional concentration of clinically relevant DNA in the biological sample.

[0099] Although measurement of the fractional concentration of clinically relevant DNA in a calibration sample may be required, the resulting calibration value (e.g., as part of a calibration function) may be used to determine the fractional concentration of a new sample without identifying the alleles that are unique to the clinically relevant DNA. In this way, fractional concentrations may be determined in a more robust manner.

[0100] A. Pregnancy Genotypic differences between the maternal and fetal genomes can be used to distinguish fetal and maternal DNA molecules, for example by exploiting informative single nucleotide polymorphism (SNP) sites where the mother is homozygous (AA) and the fetus is heterozygous (AB).

[0101] FIG. 2 shows a schematic diagram of a genotypic difference-based approach for analyzing differential terminal motif patterns between fetal and maternal DNA molecules according to an embodiment of the present disclosure. As shown in FIG. 2, fetal-specific molecules 205 carrying fetal-specific alleles (B) can be determined. On the other hand, shared molecules 207 carrying shared alleles (A) can be determined, which represent DNA molecules mainly derived from the mother, since fetal DNA molecules are generally the minority in the maternal plasma DNA pool. Therefore, the characteristics of any molecules derived from shared molecules reflect the characteristics of the maternal background DNA molecules (i.e., DNA molecules derived from the hematopoietic system). In addition to alleles, other fetal-specific markers (e.g., epigenetic markers) can be used.

[0102] Using the technique 140 of FIG. 1, 4-mer terminal motifs were analyzed. 256 terminal motifs were analyzed. The percentage of each 4-mer motif was calculated and the frequency was compared across the 256 motifs using a bar graph shown as bar graph 220. Such a bar graph provides the relative frequency (%) at which each 4-mer occurs as a terminal motif. For ease of illustration, only a few 4-mers are shown. The relative frequency (sometimes simply referred to as "frequency") may be determined by (number of DNA fragments with terminal motif) / total number of DNA fragments analyzed, with a factor of two in the denominator to count both ends. Such a percentage may be considered a relative frequency because it relates to the ratio of one amount (e.g., count) for the first terminal motif to the amount of one or more other motifs (potentially including the first terminal motif). As can be seen, the terminal motif 222 has significant differences in relative frequency among DNA fragments of various tissue types. Such differences may be used for various purposes, such as enriching a sample for fetal DNA or determining fetal DNA concentration.

[0103] The relative frequency values ​​shown in the bar graph 220 may be stored in an array with 256 values. A counter may exist for each terminal motif in the set of terminal motifs, and the counter for a particular terminal motif is incremented each time a new DNA fragment has the terminal motif that corresponds to that counter. The set of motifs may be selected in various ways, for example, all terminal motifs, or a smaller set, such as those that occur most frequently in the reference sample or those that show the greatest separation in the reference sample.

[0104] Various quantification techniques can be used to provide a measure of the relative frequency of a sample, and such quantification techniques can be used to classify the amount of cell-free DNA from clinically relevant DNA. One example of quantification technique includes the sum of the relative frequency of a set of terminal motifs, also referred to herein as composite frequency. As an example, such a set can be the terminal motifs that are identified as occurring most frequently in a particular tissue type, or having the greatest separation between two tissue types. Weighted sums can also be used. Weights can be pre-determined or variable, for example, the weight of a given frequency can depend on the frequency itself. Entropy is such an example.

[0105] In another embodiment, an entropy-based analysis 230 can be used to capture the differences in the profile of terminal motifs between fetal and maternal DNA molecules. Entropy is an example of variance / diversity. To analyze the distribution of motif frequencies (e.g., a total of 256 motifs), one definition of entropy uses the following equation:

number

[0106] In this example, if the 256 motifs are equally present in frequency, the entropy reaches a maximum value (i.e., 5.55). In contrast, if the 256 motifs have a skewed distribution in frequency, the entropy decreases. For example, if a particular motif accounts for 99% and other motifs make up the remaining 1%, in this formulation the entropy decreases to 0.11, although other formulations, such as using no log or only log, may be used. Thus, a decrease in the entropy of motif frequency implies an increase in skewness in the frequency distribution over the terminal motifs. An increasing entropy of motif frequency suggests that the frequency over the motifs shifts towards equal probability of those motifs. Thus, the entropy of motif frequency measures how evenly present the abundance of terminal motifs is in plasma DNA. The greater the degree of uniformity in motif frequency, the higher the entropy value is expected. In other words, a decrease in the entropy of motif frequency implies an increase in skewness over the terminal motifs in terms of their frequency.

[0107] In various other examples, the standard deviation (SD), coefficient of variation (CV), interquartile range (IQR) or a particular percentile cutoff (e.g., 95th or 99th percentile) between the frequencies of various motifs can be used to assess the landscape changes of terminal motif patterns between fetal and maternal DNA molecules. Such various examples provide a measure of dispersion / diversity in the relative frequencies for a set of terminal motifs. Considering the definition of entropy in FIG. 2, if only one terminal motif has a non-zero count, the entropy has a minimum value. If other terminal motifs appear in some DNA fragments, the entropy will increase. In the absence of selection (e.g., random distribution for all terminal motifs in one hypothetical scenario where all have the same frequency), the entropy will be at a maximum value. In this way, the entropy quantifies the overall selectivity of the terminal sequences of cell-free DNA fragments for terminal motifs.

[0108] Plot 235 shows the entropy values ​​of the shared sequences (mainly maternal) and fetal sequences. The shared sequences contain less fetal DNA (approximately 5% if the original sample contains 10% fetal DNA) than the fetal sequences that have nearly 100% fetal DNA, within the margin of error for the genotyping measurement. Given this separation, the higher the concentration of fetal DNA in the sample, the greater the difference in the entropy values. This relationship between fetal DNA concentration and entropy can be used to determine fetal DNA concentration, for example, as measured using one or more calibration values. For example, the concentration of clinically relevant DNA may be measured on a calibration sample using another technique (resulting in a calibration value), which may not be generally applicable, such as using Y chromosome DNA for male fetuses, or using previously identified mutations for tumor tissue. Given the entropy measurements for the calibration sample, a comparison of the two entropy values ​​(one for the test sample and one for the calibration sample) can provide a fractional concentration for the test sample using the concentration measured in the calibration sample. Further details of such use of calibration values ​​and calibration functions are provided below.

[0109] In yet another embodiment, a clustering-based analysis 240 may be employed. The vertical axis corresponds to 4-mer motifs and the horizontal axis corresponds to different samples with different classifications of fetal DNA concentration, for example. The color corresponds to the relative frequency of a particular 4-mer motif for a particular sample, for example a red calibration sample 242 having a higher concentration than a green calibration sample 244, which has a lower value.

[0110] Clustering-based analysis may take advantage of the assumption that the similarity of the frequency profiles of the 256 4-mer terminal motifs is relatively high either within fetal or maternal DNA molecules (i.e., within-group molecular characteristics) compared to the similarity between fetal and maternal DNA molecules (i.e., between-group molecular characteristics). Thus, calibration samples of individuals characterized with terminal motifs derived from shared sequences (e.g., higher concentration of shared sequences) were expected to differ from calibration samples of individuals characterized with terminal motifs derived from fetal-specific sequences (e.g., lower concentration of shared sequences and therefore higher fetal). Each individual corresponded to a vector (i.e., a 256-dimensional vector) containing the 256 terminal motifs and their corresponding frequencies. Examples of clustering techniques include, but are not limited to, hierarchical clustering, centroid-based clustering, distribution-based clustering, and density-based clustering. The different clusters may correspond to different amounts of fetal DNA in the sample because they have different patterns of relative frequency due to differences in the frequency of terminal motifs between maternal and fetal DNA fragments.

[0111] To evaluate the differences in terminal motifs between fetal and maternal DNA molecules, maternal buffy coat and fetal samples were genotyped using a microarray platform (Human Omni2.5, Illumina), respectively, and matched plasma DNA samples were sequenced. Peripheral blood samples from 10 pregnant women from each of the first (12–14 weeks), second (20–23 weeks), and third (38–40 weeks) trimesters were obtained, and plasma and maternal buffy coat samples from each situation were collected. A median of 195,331 informative SNPs (range: 146,428–202,800) for which the mother was homozygous and the fetus was heterozygous were obtained. Plasma DNA molecules carrying fetal-specific alleles were identified as fetal-specific DNA molecules. Plasma DNA molecules carrying shared alleles were identified and are believed to be primarily maternal-derived DNA molecules. The median fetal DNA fraction among these samples was 17.1% (range: 7.0%-46.8%). A median of 103 million mapped paired-end reads (range: 52-186 million) were obtained for each situation. The terminal motifs of each plasma DNA molecule were determined by bioinformatically interrogating the 4-mer sequences closest to the fragment ends. The results of the analysis of this sample set are provided below.

[0112] 1. Differences in the relative frequency of rank orders We considered that the top terminal motifs in the ranked difference in motif frequency between fetal and maternal DNA molecules would be useful for the detection or enrichment of fetal and maternal DNA molecules. Therefore, we ranked the terminal motifs in terms of the difference in frequency between fetal and maternal DNA molecules in one pregnant woman at a sequencing depth of 270x. Fetal and shared sequences were identified according to informative SNPs using the same method as above.

[0113] FIG. 3 shows a bar graph of terminal motif frequency between fetal and maternal DNA molecules according to an embodiment of the present disclosure. Data was obtained from one pregnant woman at a sequencing depth of 270x. The vertical axis corresponds to the frequency percentage of a given 4mer motif determined from the number of DNA fragments (determined from sequence reads), which is the percentage of the given 4mer motif divided by the total number of terminal sequences of the analyzed DNA fragments (e.g., twice the number of DNA fragments). The horizontal axis corresponds to 256 4mers. The 4mers are arranged in order of decreasing frequency of shared sequences, and FIG. 3 is divided into two parts with different scales used for the vertical axis. Differences in the frequency of terminal motifs were observed between fetal DNA molecules (those with fetal-specific alleles) and maternal DNA molecules (those with shared alleles).

[0114] Figure 4 shows the top 10 terminal motifs from Figure 3 for fetal and shared (i.e. fetal plus maternal) sequences according to an embodiment of the present disclosure. The vertical axis is shifted to start at 1% frequency. The top 10 terminal motifs are CCCA, CCAG, CCTG, CCAA, CCCT, CCTT, CCAT, CAAA, CCTC, and CCAC. As can be seen, some terminal motifs have a larger difference between the shared and fetal specific sequences than others. Thus, the terminal motif with the largest difference may be used to distinguish between maternal and fetal DNA, as opposed to simply the most frequent terminal motif.

[0115] 2. Use of Entropy Next, for the various samples, the entropy of DNA molecules with shared alleles, and the entropy of DNA molecules with fetal-specific alleles, the former identified as maternal and the latter identified as fetal, are analyzed. For each sample, two data points are obtained: the entropy of the fetal DNA molecule and the entropy of the shared DNA molecule (labeled "maternal").

[0116] Figure 5A shows that the entropy of the terminal motifs in fetal DNA molecules is lower than that in maternal DNA molecules (p-value < 0.0001), suggesting a higher skewness in the distribution of terminal motifs originating from maternal DNA molecules. The entropy in Figure 5A is determined using all 256 motifs, with 4mers used in these examples, for a given sample and for a given pool of fetal or maternal DNA molecules.

[0117] Similar to plot 235 in FIG. 2, the difference in entropy for the two tissue types indicates that entropy can be used to determine the fractional concentration of fetal DNA in a mixture of cell-free DNA fragments (e.g., plasma or serum). As described above, the pool identified as fetal DNA has a higher percentage of fetal DNA (e.g., nearly 100%) than the maternal pool. The entropy values ​​determined for the types of pools are different. Thus, there is a relationship between entropy and fetal DNA concentration. This relationship can be determined as a calibration function based on the measurements of fetal DNA concentration of calibration samples (calibration values) and the corresponding entropy values ​​(examples of relative frequency), and the calibration values ​​and relative frequencies can form calibration data points. Calibration samples with different fetal DNA concentrations have different entropy values. A calibration function can be fitted to the calibration data points such that the newly measured relative frequencies (e.g., entropy) can be input into the calibration function to provide an output of fetal DNA concentration.

[0118] Figure 5B shows the entropy when using the relative frequency of the 10 motifs of Figure 4. As shown, the relationship changes with fetal sequences having higher entropy for this given set of 10 terminal motifs. The fractional concentration of fetal DNA can still be determined, but a different calibration function would be used. Thus, the set of motifs used for calibration needs to be the same as that used later, i.e., when measuring fractional concentration based on entropy, or other summary value of the relative frequency of the set.

[0119] 3. Clustering Furthermore, we performed a hierarchical clustering analysis on pregnant women. Each pregnant woman was characterized by a 256-dimensional vector containing the terminal motif frequencies of all 4-mers. Indeed, individuals characterized by fetal-specific sequences and terminal motifs derived from maternal DNA molecules could be clustered into two groups.

[0120] Figures 6A and 6B show hierarchical clustering analysis of first trimester fetal and maternal DNA molecules according to an embodiment of the present disclosure. Figure 6A shows a hierarchical clustering analysis based on 256 4mer end motif frequencies. The vertical axis corresponds to the 4mer motifs and the horizontal axis corresponds to the various portions of the different samples (i.e., 620 fetal-specific (yellow) and 610 shared (blue) sequences). The colors correspond to the relative frequency of a particular 4mer motif in a particular portion of the sample.

[0121] The different parts (fetal-specific and shared) have different fetal DNA concentrations and therefore different classifications for the concentration of fetal DNA. If such clustering is performed using calibration samples, the fetal DNA concentration can be measured, for example, as described in the entropy section above. Each calibration sample has a corresponding vector of length equal to the number of motifs used (e.g., 256 for all 4-mers or potentially only a subset of 4-mers, as other k-mers can be used, but may have the greatest difference between fetal and shared sequences).

[0122] FIG. 6B shows a zoomed-in visualization for the hierarchical clustering analysis based on the 256 4-mer terminal motif frequencies. Each row represents one type of terminal motif (i.e., different terminal motifs). Each column represents pregnant subjects. The gradient colors indicate the frequency of the terminal motifs, with red representing the highest frequency and green representing the lowest frequency. As can be seen, the two parts (fetal and shared) representing samples with different fetal DNA concentrations are neatly clustered into two separate clusters, indicating a good accuracy with which samples with different levels of fetal DNA concentration can be discriminated.

[0123] 4. Samples from various trimesters In addition to being able to distinguish between samples with different fraction concentrations, some embodiments can distinguish between different samples from pregnant subjects at different gestational ages (e.g., which trimester, or whether it is just the third trimester).

[0124] 7A and 7B show the entropy distribution using all motifs of pregnant women across different trimesters according to an embodiment of the present disclosure. Interestingly, the entropy value of the number of terminal motifs determined using fetal-specific fragments seems to be related to gestational age (p-value: 0.024, first trimester data vs. pooled data from second and third trimesters), but that from shared fragments (mainly maternal DNA) does not seem to be associated with gestational age (p-value: 1, first trimester data vs. pooled data from second and third trimesters). Later in pregnancy, fetal DNA concentration is generally higher. Therefore, there may be a correlation between concentration and gestational age.

[0125] For fetal-specific fragments, the entropy was decreased in the second and third trimesters compared to the first trimester. Thus, fetal fragments may convey gestational age. And since the shared fragments have essentially constant entropy (e.g., because they are primarily maternal fragments and / or because of maternal physiology-related changes in terminal motifs that counteract such fetal signals), the change in entropy for all fragments reflects gestational age due to changes in fetal fragments. Although such a relationship of entropy between various trimesters shows less change due to the presence of maternal fragments, the relationship still exists. However, if fetal-specific alleles can be identified (e.g., by identifying alleles occurring in similar percentages in male fetuses or expected fetal DNA concentrations, or using paternal genotype information), a more pronounced relationship exists (e.g., as shown in Figure 7B).

[0126] 7C and 7D show entropy distributions using 10 motifs for pregnant women across different trimesters according to an embodiment of the present disclosure. The 10 motifs were selected by ranking determined from the shared fragments. These figures show that even though the specific selection of motifs may decrease the association (as opposed to the increase in FIG. 7B), the entropy of the different trimesters for the fetal-specific fragments still changes.

[0127] Figure 8A shows the entropy of all fragments across various gestational ages according to an embodiment of the present disclosure. The entropy is determined using all 256 4-mer terminal motifs. It was shown that the entropy of plasma DNA fragments in subjects in the third trimester is lower (p-value=0.06) compared to those in the first and second trimesters. And the average of the second trimester is lower than the first trimester. Thus, when all of the fetal fragments are included (as opposed to the shared fragments in Figure 7A), the entropy provides a gestational age.

[0128] Figure 8B shows the entropy of Y-chromosome-derived fragments across different gestational ages. It is shown that the entropy of Y-chromosome-derived fragments in subjects in the third trimester is lower than that in the first and second trimesters (p-value=0.01). These samples that have been filtered for fetal molecules (using fetal-specific sequences from Y-chromosome) show a greater separation between the third and second trimesters.

[0129] 9 and 10 show the distribution of the top 10 ranked terminal motifs between fetal and maternal DNA molecules across different trimesters according to an embodiment of the present disclosure. The top 10 terminal motifs of the ranked difference in motif frequency between fetal and maternal DNA molecules were mined from one single deep sequenced pregnancy case. These top 10 terminal motifs were then used to analyze each sample.

[0130] The percentage of fetal and shared DNA molecules carrying these terminal motifs of interest was calculated in an independent cohort of 10 pregnant women from each of the first (12-14 weeks), second (20-23 weeks), and third (38-40 weeks) trimesters. There were numerous terminal motifs that were found to be higher in fetal DNA molecules compared to shared molecules, suggesting that those terminal motifs have a specific relationship to the tissue of origin. For example, the median CAAA% was found to be consistently higher in fetal DNA molecules than in shared molecules (mainly maternal) across the first (1.26% vs. 1.11%), second (1.24% vs. 1.11%), and third (1.24% vs. 1.15%) trimesters. Thus, the terminal motif CAAA may be identified as a marker indicating an increased likelihood that a particular DNA fragment with a terminal sequence of CAAA is of fetal origin.

[0131] Certain terminal motifs show a more pronounced relationship with gestational age. For example, fetal DNA molecules with the terminal motif CCCA show a continuous (monotonically) increase with gestational age, as do CCAG, CCTG, CCAA, CCCT, and CCAC. However, CCTT does not show a continuous increase, as the median value declines in the second trimester and then increases in the third trimester.

[0132] In another embodiment, the top 10 ranked terminal motifs can be combined to look at the differences between fetal and maternal DNA molecules across different trimesters.

[0133] FIG. 11 shows the combined frequency of the top 10 ranked motifs between fetal and shared molecules across different trimesters according to an embodiment of the present disclosure. As shown in FIG. 11, the difference in combined frequency of the top 10 ranked motifs between fetal and maternal DNA molecules was found to be relatively larger in both the second trimester (p-value: 0.013) and third trimester (P-value: 0.0019) compared to the first trimester (p-value: 0.92). The frequency for fetal molecules increases continuously from the first trimester to the second trimester to the third trimester, but this continuous relationship is not shown for shared molecules. This indicates that different physiological conditions (e.g., gestational age) affect the terminal motifs derived from different tissues of origin.

[0134] B. Oncology Genotypic tools devised in the pregnancy context may also be applied in the oncology context.

[0135] Figure 12 shows a schematic diagram of a genotypic difference-based approach to analyze differential terminal motif patterns between mutant and shared molecules in plasma DNA of cancer patients according to an embodiment of the present disclosure. As shown in Figure 12, tumor-specific molecules 1205 carrying tumor-specific alleles (B) can be determined. On the other hand, shared molecules 1207 carrying shared alleles (A) can be determined, which would mainly represent healthy-derived DNA molecules since tumor DNA molecules are generally in the minority in the plasma DNA pool.

[0136] As an example, mutant sequences (i.e., plasma DNA carrying cancer-associated mutations) and shared sequences (mainly DNA of hematopoietic origin) may be identified. Cancer-associated mutations may be defined as variants present in tumor tissue (hepatocellular carcinoma, HCC) but absent in normal cells (e.g., buffy coat). For example, in an HCC patient, assuming that the genotype of tumor tissue is "AG" at a particular genomic locus and the genotype of buffy coat cells is "AA", the "G" present specifically in tumor tissue is considered to be a cancer-associated mutation and the "A" is considered to be a shared wild-type allele. In various implementations, mutant sequences may be obtained by sequencing tissue biopsies from the tumor or by analyzing acellular samples such as plasma or serum, for example, as described in U.S. Patent Publication No. 2014 / 0100121.

[0137] The frequency profile of terminal motifs between variant and shared sequences was determined in HCC patients whose plasma DNA was sequenced at a depth of 220x. Bar graph 1220 provides the relative frequency (%) at which each 4-mer occurs as a terminal motif for variant and shared sequences. Such relative frequencies can be determined as in bar graph 220 of FIG. 2 above. As can be seen, terminal motifs 1222 have significant differences in relative frequency between DNA fragments of different tissue types. Such differences can be used for various purposes, for example, to enrich samples for tumor DNA or to measure tumor DNA concentration.

[0138] In another embodiment, an entropy-based analysis 1230 can be used to capture differences in the profile of terminal motifs between tumor and shared DNA molecules, similar to FIG. 2. Plot 1235 shows entropy values ​​for the shared and tumor sequences. The difference in entropy or other dispersion metric can provide the tumor fraction concentration, for example, using a calibration function.

[0139] In yet another embodiment, a clustering-based analysis 1240 can be performed similar to the fetal analysis of Figure 2. A classification for the amount of tumor sequences in a sample can be determined based on a new sample belonging to a reference cluster with a known classification of the tumor fraction.

[0140] 1. Differences in the relative frequency of rank orders Figure 13 shows the plasma DNA terminal motif profile of cancer-associated variants and shared molecules in hepatocellular carcinoma according to an embodiment of the present disclosure. There are many terminal motifs that have been observed to vary between variant and shared sequences, including but not limited to CCCA, CCAG, CCAA, CCTG, CCTT, CCCT, CAAA, CCAT, TAAA, AAAA motifs. Figure 13 shows similar information to Figure 3, but for clinically relevant DNA, it is tumor DNA instead of fetal DNA.

[0141] 14 shows a radial view of plasma DNA terminal motifs of cancer-associated variants and shared molecules in hepatocellular carcinoma according to an embodiment of the present disclosure. The various terminal motifs are listed on the periphery, and the frequency of the terminal motifs is shown at various radial lengths. The terminal motifs are sorted by frequency of wild-type (wt) alleles in non-tumor (e.g., healthy) cells. Frequency value 1410 corresponds to the wt allele, and frequency value 1420 corresponds to the mutant (mut) allele. This radial view shows a significant difference in the relative frequency of terminal motifs in the mutant sequences compared to the wild-type (shared) sequences.

[0142] Figure 15A shows the top 10 terminal motifs in the rank order difference of terminal motif frequency between mutant and shared sequences in plasma DNA of HCC patients according to an embodiment of the present disclosure. The top terminal motifs are determined for the shared sequences in the reference sample. As shown, the top terminal motifs are CCCA, CCAG, CCAA, CCTG, CCTT, CCCT, CAAA, CCAT, TAAA, and AAAA. The difference in relative frequency varies between terminal motifs. For example, the motif (CCCA) that shows the largest difference between mutant and shared sequences is found to be 1.9% and 1.6%, respectively, suggesting that such motifs are reduced by 15% in mutant sequences compared to shared sequences (mainly wild-type sequences from blood cells).

[0143] Figure 15B shows the composite frequency of eight terminal motifs for HCC patients and pregnant women according to an embodiment of the present disclosure. The composite frequency is an exemplary summary value, for example, as the sum of the relative frequencies of a set of terminal motifs. As can be seen, there is a separation in the composite frequency for two classes of sequences in each of these two scenarios, wild type (WT) and mutant, and between maternal and fetal sequences. The separation of the composite frequency for wild type (WT) and mutant is greater than the separation for maternal and fetal sequences.

[0144] This composite frequency behaves similarly to the entropy plot for fetal analysis. Thus, Figure 15B shows another example of a relative frequency summary that can be used to determine the fractional concentration of clinically relevant DNA. And the wt vs. mutant relationship in Figure 15B shows that the fractional concentration of other clinically relevant DNA (e.g., tumor DNA) can also be determined.

[0145] 2. Use of entropy Figures 16A and 16B show entropy values ​​for shared and variant fragments of different sets of terminal motifs for HCC cases according to an embodiment of the present disclosure. As with fetal sequences, the relationship between the entropy for two sequences varies depending on the set of terminal motifs used. Figure 16A uses all 256 terminal motifs for 4mers. The entropy is higher for variant fragments because the frequency distribution for variant fragments is more uniform (e.g., flatter). Also, the entropy of shared fragments is lower because the frequency distribution is more skewed.

[0146] Figure 16B uses the top 10 terminal motifs of 4mers occurring in HCC subjects for shared fragments. The entropy relationship is reversed for the top 10 motifs. Both Figures 16A and 16B show that the calibration analysis for determining fetal DNA concentration can also be used to determine tumor DNA concentration.

[0147] As mentioned above, a high entropy value indicates high diversity in the terminal motifs. The motif diversity score (MDS) can be used to estimate the fractional concentration of clinically relevant DNA (e.g., fetal, transplant, tumor, etc.) in biological samples of circulating cell-free DNA.

[0148] 17 is a plot of motif diversity score versus measured circulating tumor DNA fraction according to an embodiment of the present disclosure. For each of a plurality of calibration samples, a calibration data point 1705 was measured. The calibration data point includes the motif diversity score for the sample and the fractional concentration of clinically relevant DNA, in this case the tumor DNA fraction. The tumor DNA fraction was estimated based on the software package ichorCNA, which measures tumor DNA fraction in plasma DNA by utilizing cancer-associated copy number aberrations (Adalsteinsson et al. 2017).

[0149] A given sample can be a healthy control sample with no tumor DNA, or a sample from a patient with a tumor in which the tumor DNA fraction is non-zero, i.e., tumor DNA and other (e.g., healthy) DNA are present. The MDS values ​​of plasma DNA from HCC patients were found to be positively correlated with the tumor DNA fraction (Spearman's ρ: 0.597, p-value: 0.0002). This is illustrated by the calibration function 1710 (a linear function in this example).

[0150] The calibration function 1710 can be used to determine the tumor DNA fraction in the new test sample for which the motif diversity score was measured. The calibration function 1710 can be determined by a functional fit to the calibration data points 1705, for example using regression.

[0151] In some examples, the calculated value X of MDS for a new test sample can be used as an input to a function F(X), where F is a calibration function (curve). The output of F(X) is the fractional concentration. An error range can be provided that can be different for each X value, thereby providing a range of values ​​as the output of F(X). In other examples, the fractional concentration corresponding to a measured value of 0.95 for MDS in a new sample can be determined as the average concentration calculated from the calibration data points at MDS 0.95. As another example, the calibration data points 1705 can be used to provide a range of fractional DNA concentrations for a particular calibration value, which can be used to determine whether the fractional concentration exceeds a threshold amount.

[0152] C. Transplant Genotypic techniques can also be applied to monitor transplants, e.g., liver transplants. SNP sites where the recipient is homozygous and the donor is heterozygous will allow for the determination of donor-specific DNA molecules and mainly hematopoietic DNA in the plasma in transplant patients.

[0153] FIG. 18A shows an entropy analysis using donor-specific fragments according to an embodiment of the present disclosure. FIG. 18B shows a hierarchical clustering analysis using donor-specific fragments. As shown in FIG. 18A and FIG. 18B, in the context of liver transplantation, liver-specific DNA molecules were observed to have different characteristics from shared sequences (mainly blood-derived DNA). The entropy of plasma DNA terminal motifs was generally found to be lower in donor-specific DNA molecules (liver DNA) compared to shared sequences (FIG. 18A). Individuals characterized by terminal motifs derived from liver-specific DNA molecules were clustered together, and individuals characterized by terminal motifs derived from shared DNA molecules were clustered into another group.

[0154] D. Classification of fraction concentrations As described above, the relative frequency of a set of one or more terminal motifs can be used to determine a classification of the fractional concentration of clinically relevant DNA.

[0155] 19 is a flow chart illustrating a method 1900 for estimating the fractional concentration of clinically relevant DNA in a biological sample of a subject according to an embodiment of the present disclosure. The biological sample may contain clinically relevant DNA and other DNA that is acellular. In other examples, the biological sample may not contain clinically relevant DNA and the estimated fractional concentration may show zero or a low percentage of clinically relevant DNA. Method 1900 and any other method aspects described herein may be implemented by a computer system.

[0156] At block 1910, a plurality of cell-free DNA fragments from the biological sample are analyzed to obtain sequence reads. The sequence reads may include terminal sequences corresponding to the ends of the plurality of cell-free DNA fragments. By way of example, the sequence reads may be obtained using sequencing or probe-based techniques, any of which may include, for example, enrichment via amplification or capture probes.

[0157] Sequencing can be carried out in various ways, for example, by using massively parallel sequencing or next generation sequencing, by using single molecule sequencing, and / or by using double-stranded or single-stranded DNA sequencing library preparation protocols.Those skilled in the art will understand the various sequencing techniques that can be used.As part of sequencing, it is possible that some of the sequence reads can correspond to cellular nucleic acid.

[0158] Sequencing can be targeted sequencing, for example as described herein.For example, biological sample can be enriched for DNA fragments from specific regions.Enrichment can include using capture probes that bind to part or the whole genome, for example as defined by reference genome.

[0159] A statistically significant number of cell-free DNA molecules can be analyzed to provide accurate determination of fractional concentration.In some embodiments, at least 1,000 cell-free DNA molecules are analyzed.In other embodiments, at least 10,000 or 50,000 or 100,000 or 500,000 or 1,000,000 or 5,000,000 or more cell-free DNA molecules can be analyzed.

[0160] At block 1920, for each of the plurality of cell-free DNA fragments, a sequence motif is determined for each of one or more terminal sequences of the cell-free DNA fragment. The sequence motif may include an N base position (e.g., 1, 2, 3, 4, 5, 6, etc.). By way of example, the sequence motif may be determined by analyzing sequence reads at ends corresponding to the ends of the DNA fragments, correlating signals with particular motifs (e.g., if probes are used), and / or aligning the sequence reads to a reference genome, e.g., as described in FIG. 1.

[0161] For example, after sequencing by a sequencing device, sequence reads can be received by a computer system that can be communicatively coupled to the sequencing device that performs the sequencing, for example, via wired or wireless communication or a removable storage device. In some implementations, one or more sequence reads including both ends of the nucleic acid fragment can be received. The location of the DNA molecule can be determined by mapping (aligning) one or more sequence reads of the DNA molecule to a respective portion, e.g., a specific region, of the human genome. In other embodiments, a specific probe (e.g., after PCR or other amplification) can indicate the location or a specific terminal motif, such as via a specific fluorescent color. The identification can be that the cell-free DNA molecule corresponds to one of a set of sequence motifs.

[0162] At block 1930, the relative frequency of a set of one or more sequence motifs corresponding to the terminal sequences of the plurality of cell-free DNA fragments is determined. The relative frequency of the sequence motifs can provide the proportion of the plurality of cell-free DNA fragments having terminal sequences corresponding to the sequence motifs. The set of one or more sequence motifs can be identified using a reference set of one or more reference samples. Genotypic differences can be determined so that differences between the terminal motifs of clinically relevant DNA and other DNA (e.g., healthy DNA, maternal DNA, or DNA of a subject who received a transplanted organ) can be identified, but without knowing the fractional concentration of clinically relevant DNA for the reference sample. A particular terminal motif can be selected based on the difference (e.g., selecting the terminal motif with the largest absolute or percentage difference). Examples of relative frequencies are described throughout this disclosure.

[0163] In some implementations, the sequence motif comprises N base positions, and the set of one or more sequence motifs comprises all combinations of N bases. In some examples, N can be an integer equal to or greater than 2 or 3. The set of one or more sequence motifs can be the top M (e.g., 10) most frequently occurring sequence motifs in one or more calibration samples or other reference samples not used to calibrate the fractional concentrations.

[0164] At block 1940, a summary value of the relative frequencies of the set of one or more sequence motifs is determined. Exemplary summary values ​​are described throughout the disclosure, including, for example, an entropy value (motif diversity score), a sum of relative frequencies, and a multidimensional data point corresponding to a vector of counts for the set of motifs (e.g., a vector of 256 counts for 245 possible 4-mer motifs, or 64 counts for 64 possible 3-mer motifs). When the set of one or more sequence motifs includes multiple sequence motifs, the summary value may include a sum of the relative frequencies of the set.

[0165] As an example, when the set of one or more sequence motifs includes a plurality of sequence motifs, the aggregate value may include the sum of the relative frequencies of the set. As another example, the aggregate value may correspond to the variance of the relative frequencies. For example, the aggregate value may include an entropy term. The entropy term may include a sum of terms, each of which may include a relative frequency multiplied by the logarithm of the relative frequency. As another example, the aggregate value may include a final output or intermediate output of a machine learning model, such as a clustering model.

[0166] At block 1950, a classification of the fractional concentration of clinically relevant DNA in the biological sample is determined by comparing the tabulated value to one or more calibration values. The one or more calibration values ​​may be determined from one or more calibration samples in which the fractional concentration of clinically relevant DNA is known (e.g., measured). The comparison may be to a plurality of calibration values. The comparison may occur by inputting the tabulated value into a calibration function that is fitted to calibration data that provides a change in the tabulated value relative to a change in the fractional concentration of clinically relevant DNA in the sample. As another example, the one or more calibration values ​​may correspond to one or more tabulated values ​​of the relative frequency of a set of one or more sequence motifs measured using cell-free DNA fragments in one or more calibration samples.

[0167] The calibration value may be calculated as an aggregate value for each calibration sample. A calibration data point may be determined for each sample, the calibration data point including the calibration value and the fractional concentration measured for the sample. These calibration data points may be used in the method 1900 or may be used to determine the final calibration data points (e.g., as defined via a function fit). For example, a linear function may be fitted to the calibration value as a function of the fractional concentration. The linear function may define the calibration data points used in the method 1900. The new aggregate value of the new sample may be used as an input to the function as part of the comparison to provide the output fractional concentration. Thus, the one or more calibration values ​​may be a plurality of calibration values ​​of a calibration function determined using the fractional concentrations of clinically relevant DNA of a plurality of calibration samples.

[0168] As another example, the new tabulated value can be compared to the average tabulated value of samples having the same classification (e.g., within the same range) of fractional concentrations, and if the new tabulated value is closer to this average than a calibration value to the average of another classification, the new sample can be determined to be of the same concentration as the closest calibration value. Such techniques can be used when performing clustering. For example, the calibration value can be a representative value for a cluster corresponding to a particular classification of fractional concentrations.

[0169] Determining the calibration data points may include measuring fractional concentrations, for example, as follows: For each calibration sample of one or more calibration samples, the fractional concentration of clinically relevant DNA may be measured in the calibration sample. A summary value of the relative frequency of a set of one or more sequence motifs may be determined by analyzing cell-free DNA fragments from the calibration sample as part of obtaining the calibration data points, thereby determining one or more summary values. Each calibration data point may specify the measured fractional concentration of clinically relevant DNA in the calibration sample and the summary value determined for the calibration sample. The one or more calibration values ​​may be one or more summary values ​​or may be determined using one or more summary values ​​(e.g., when using a calibration function). Measuring fractional concentrations may be performed by various methods as described herein, for example, by using alleles specific to clinically relevant DNA.

[0170] In various embodiments, measuring the fractional concentration of clinically relevant DNA can be performed using tissue-specific alleles or epigenetic markers, or using DNA fragment sizes, for example, as described in U.S. Patent Publication No. 2013 / 0237431, which is incorporated by reference in its entirety. Tissue-specific epigenetic markers can include DNA sequences that are indicative of tissue-specific DNA methylation patterns in the sample.

[0171] In various embodiments, the clinically relevant DNA may be selected from the group consisting of fetal DNA, tumor DNA, DNA from transplanted organs, and specific tissue types (e.g., from specific organs). The clinically relevant DNA may be of a specific tissue type, for example, the specific tissue type is liver or hematopoietic. If the subject is a pregnant woman, the clinically relevant DNA may be placental tissue corresponding to fetal DNA. As another example, the clinically relevant DNA may be tumor DNA derived from an organ with cancer.

[0172] Generally, the one or more calibration values ​​determined from the one or more calibration samples are preferably generated using an assay similar to that used for the biological (test) sample whose fractional concentration is being measured. For example, a sequencing library can be generated in the same way. Two examples of processing techniques are GeneRead (www.qiagen.com / us / shop / sequencing / generead-size-selection-kit / #orderinginformation) and SPRI (Solid Phase Reversible Immobilization, AMPure Beads, www.beckman.hk / reagents_depr / genomic_depr / cleanup-and-size-selection / pcr). GeneRead can remove short DNA, mainly tumor fragments, which can affect the relative frequency of terminal motifs in fetal and transplant cases as well as wild type and mutant fragments.

[0173] E. Determination of Gestational Age As described above in Figures 7A, 7B, and 8-10, fetal-specific fragment motifs can be used to infer gestational age.

[0174] 20 is a flow chart illustrating a method 2000 for determining the gestational age of a fetus by analyzing a biological sample from a female subject carrying a fetus, according to an embodiment of the present disclosure. The biological sample includes cell-free DNA molecules from the female subject and the fetus.

[0175] At block 2010, a plurality of cell-free DNA fragments from a biological sample are analyzed to obtain sequence reads. The sequence reads may include terminal sequences corresponding to ends of the plurality of cell-free DNA fragments. Block 2010 may be performed in a manner similar to block 1910 of FIG. 19.

[0176] Before, after, or as part of the analysis, the plurality of cell-free DNA fragments may be identified as originating from the fetus, for example, as described above for Figures 2 and 5A. This may filter for DNA fragments that are fetal or most likely fetal. As an example, the plurality of cell-free DNA fragments may be identified using fetal-specific alleles or fetal-specific epigenetic markers. As another example, for each of the sequence reads, the likelihood that the sequence read corresponds to a fetus may be determined based on the end sequence of the sequence read containing a sequence motif of a set of one or more sequence motifs. Other criteria may also be used, for example, as described in Section II.E. The likelihood may be compared to a threshold, and the sequence read may be identified as originating from the fetus if the likelihood exceeds the threshold. Further details of methods for enriching a sample for clinically relevant DNA may be found in Section IV.

[0177] At block 2020, for each of the plurality of cell-free DNA fragments, a sequence motif is determined for each of the one or more end sequences of the cell-free DNA fragment. Block 2020 may be performed in a manner similar to block 2020 of FIG.

[0178] At block 2030, a relative frequency of a set of one or more sequence motifs corresponding to end sequences of the plurality of cell-free DNA fragments is determined. The relative frequency of the sequence motifs can provide a proportion of the plurality of cell-free DNA fragments having end sequences corresponding to the sequence motifs. Block 2030 can be performed in a manner similar to block 1930 of FIG.

[0179] A summary of the relative frequencies of the set of one or more sequence motifs is determined at block 2040. Block 2040 may be performed in a manner similar to block 1940 of FIG.

[0180] At block 2050, one or more calibration data points are obtained. Each calibration data point may specify a gestational age (e.g., a trimester as described in the figure above) that corresponds to the aggregate value. As described above, the one or more calibration data points may be determined from a plurality of calibration samples having known gestational ages and including cell-free DNA molecules. In some implementations, the one or more calibration data points may be a plurality of calibration data points that form a calibration function that approximates a measured aggregate value determined from cell-free DNA molecules in a plurality of calibration samples having known gestational ages.

[0181] At block 2060, the aggregate value is compared to the calibration value of at least one calibration data point. For example, the new aggregate value of the new sample can be compared to the average for the third trimester as determined in FIG. 8A. As another example, the calibration value of the at least one calibration data point can correspond to an aggregate value measured using cell-free DNA molecules in at least one of a plurality of calibration samples. The comparison of the aggregate value can be for a plurality of calibration values, e.g., each corresponding to one of a plurality of calibration samples. The comparison can occur by inputting the aggregate value into a function fit (calibration function) to the calibration data that provides the change in the aggregate value with respect to gestational age. The comparison can be performed in a manner similar to that described for method 1900, e.g., with respect to block 1950.

[0182] At block 2070, the gestational age of the fetus is estimated based on the comparison. If the new aggregate value is closest to the third trimester average (or other calibration value used), then the new sample may be determined to be third trimester. As another example, the new aggregate value may be compared to a calibration function (e.g., a linear function) that is fitted to the data in FIG. 8A or other similar figure. The function may output the gestational age, for example, as the Y value of the linear function. Other examples provided herein for using a calibration function may also be used in the context of determining gestational age.

[0183] III. Phenotypic Approaches Using genotype-based analysis of pregnant subjects, cancer subjects, as well as liver transplants, the presence of plasma DNA terminal motifs was associated with the tissue of origin. We reasoned that in cancer patients, tumor DNA is released into the blood circulation, altering the original normal presentation of plasma DNA terminal motifs. However, we do not exclude the possibility that other aspects of cancer pathobiology, such as the tumor microenvironment (infiltrating T cells, B cells, neutrophils, etc.), may generate different terminal motifs and affect the terminal motif landscape. Thus, analysis of plasma DNA terminal motifs between cancer and non-cancer control subjects would reveal the power of classifying HCC from control subjects.

[0184] Figure 21 shows a schematic diagram of a phenotypic approach for plasma DNA terminal motif analysis according to an embodiment of the present disclosure. Figure 21 is similar to Figures 2 and 12, for example, relative frequencies can be plotted, variance values ​​(e.g., entropy) can be determined, and clustering can be performed.

[0185] In Figure 21, the terminal motifs (e.g., 4mers) predicted from plasma DNA molecules are used and compared between cancer and control subjects, which eliminates the limitations of genotypic markers and allows for broad application in many clinical scenarios, such as autoimmune diseases (e.g., systemic lupus erythematosus, SLE) and transplantation. Using a phenotypic approach using all sequenced plasma DNA fragments, entropy and clustering analysis can be performed in a very similar analytical procedure to that performed in the genotypic difference-based approach. In this context, entropy and clustering analysis are compared between control and diseased subjects.

[0186] The diseased molecules 2105 are from one or more subjects determined to have the disease. The control molecules 2107 are from one or more subjects not having the disease. The relative frequency of the set of terminal motifs is determined for the two pools of molecules. Bar graph 1220 provides the relative frequency (%) at which each 4-mer occurs as a terminal motif in the control and diseased sequences. Such relative frequencies may be determined as described above for bar graph 220 in FIG. 2. As can be seen, the terminal motifs 2122 have significant differences in relative frequency among DNA fragments of various tissue types. Such differences may be used for various purposes, for example, to classify new samples as diseased or non-diseased, or some other level of disease.

[0187] To capture the differences in the profile of terminal motifs between tumor and shared DNA molecules, an entropy-based analysis 2130 can be used, similar to Fig. 2. Plot 2135 shows the entropy values ​​for control and diseased subjects. Differences in entropy or other dispersion metrics can provide a classification of the level of pathology associated with the disease.

[0188] In yet another embodiment, a clustering-based analysis 2140 may be performed, similar to the fetal analysis of Figure 2 and the tumor analysis of Figure 12. A classification of the level of pathology may be determined based on new samples belonging to reference clusters whose classification is known.

[0189] Thus, in one example of a relative frequency summary, each individual may be characterized by a vector containing 256 frequencies for 4-mer terminal motifs (i.e., a 256-dimensional vector). In other examples, standard deviations (SD), coefficients of variation (CV), interquartile ranges (IQR), or specific percentile cutoffs (e.g., 95th or 99th percentiles) between various motif frequencies may be used to assess changes in the profile of terminal motif patterns between disease and control groups. Other examples of summary values ​​are provided in other sections and are applicable here.

[0190] A. Oncology In some embodiments, the disease (pathology) may be cancer. Accordingly, some embodiments may classify the level of cancer.

[0191] 1. Differences in the relative frequency of rank orders FIG. 22 shows an example of a frequency profile of 4-mer terminal motifs between hepatocellular carcinoma (HCC) and hepatitis B virus (HBV) subjects using all plasma DNA molecules, according to an embodiment of the present disclosure. FIG. 22 compares the frequency of 256 terminal motifs in HCC patients with one HBV subject. As a similar plot, the vertical axis is the motif frequency and the horizontal axis corresponds to each terminal motif. In FIG. 22, the motifs are ranked in ascending order based on the average motif frequency in non-HCC subjects. The bottom plot continues from the top plot, but is scaled differently for ease of illustration.

[0192] There were many terminal motifs that showed abnormalities in HCC patients. For example, the top 10 ranked terminal motifs that showed increased frequency in HCC patients (TGGG, TAAA, AAAA, GAAA, GGAG, TAGA, GCAG, TGGT, GCTG, and GAGA) changed an average of 1.22-fold with a range of 1.12-1.35-fold changes, and the top 10 terminal motifs that showed decreased frequency in HCC patients (CCCA, CCAG, CCAA, CCCT, CCTG, CCAC, CCAT, CCCC, CCTC, and CCTT) changed an average of 1.23-fold with a range of 1.16-1.29-fold changes. A set of such top motifs that show an increase in its frequency (or a decrease as a separate set) in the HCC group compared to the non-cancer group can be used to classify new subjects for cancer. As another example, the ranking process can select all motifs that show an increase in HCC and rank them in descending order according to AUC between HCC and non-HCC subjects. Then, the top 10 motifs are selected based on the AUC values.

[0193] To test the diagnostic potential using plasma DNA terminal motifs, we sequenced 20 healthy control subjects (Controls), 22 chronic hepatitis B carriers (HBV), 12 cirrhotic subjects (Cirr), 24 early stage HCC (eHCC), 11 immediate stage HCC (iHCC), and 7 advanced stage HCC (aHCC) with a median of 215 million paired reads (range: 97-1,681 million).

[0194] FIG. 23A shows a box plot of the combined frequency of the top 10 plasma DNA 4-mer terminal motifs for various subjects with different levels of cancer according to an embodiment of the present disclosure. The top 10 plasma DNA 4-mer terminal motifs were selected based on the data in FIG. 22, i.e., frequency in HBV subjects. The combined frequency is the sum of the frequencies of the 10 terminal motifs for a given subject. The combined frequency of the top 10 terminal motifs was found to be significantly decreased in HCC patients compared to non-cancer subjects (p-value<0.0001). Importantly, using this terminal motif analysis, 58.3% of eHCC patients could be identified with a specificity of 95%. Furthermore, different stages of cancer could be detected. For example, advanced HCC has a significantly lower value than eHCC and iHCC.

[0195] Figure 23B shows the receiver operating characteristic (ROC) curve of the combined frequency of the top 10 plasma DNA 4-mer terminal motifs between HCC and non-cancer subjects according to an embodiment of the present disclosure. The area under the curve (AUC) of the ROC curve was found to be 0.91, indicating that plasma DNA terminal motifs indeed have clinical potential to distinguish HCC from non-cancer subjects. In another embodiment, the combined frequency of the seven terminal motifs with the greatest separation between HCC and non-HCC subjects provides an AUC of 0.92.

[0196] Figure 24A shows the box plot of the frequency of CCA motifs across different groups according to an embodiment of the present disclosure. The most frequent 3mer motif (CCA) in the non-HCC group was shown to be significantly lower in the HCC group (p-value < 0.0001). Figure 24B shows the ROC curve between non-HCC and HCC groups using the most frequent 3mer motif (CCA) present in non-HCC subjects according to an embodiment of the present disclosure. The AUC was found to be 0.915. The most frequent 4mer (CCCA) also provides a similar AUC of 0.91.

[0197] 2. Use of Entropy (Motif Diversity Score) FIG. 25A shows box plots of entropy values ​​across different groups using 256 4-mer terminal motifs according to an embodiment of the present disclosure. All 256 motifs of 4-mers were used. As shown in FIG. 25A, entropy values ​​were significantly increased (p-value<0.0001) in HCC patients (mean: 5.242, range: 5.164-5.29) compared to non-HCC subjects (mean: 5.203, range: 5.124-5.253). Importantly, using this terminal motif analysis, 41.7% of eHCC patients could be identified with 95% specificity. Entropy was generally increased in HCC, IHCC, and advanced stage HCC groups compared to the non-HCC group. Furthermore, different stages of cancer could be detected. For example, advanced HCC has significantly higher values ​​than eHCC and iHCC.

[0198] Figure 25B shows the boxplot of entropy values ​​across various groups using 10 4mer terminal motifs according to an embodiment of the present disclosure. Here, HCC subjects have reduced entropy compared to non-HCC subjects. Thus, the set of terminal motifs used can change the relationship from increasing to decreasing. For example, using the top 10 motifs reduces the entropy of the HCC group. In any case, there is diagnostic power between HCC and non-HCC groups, as well as advanced HCC, compared to early stages of HCC.

[0199] Figure 26A shows the boxplot of entropy values ​​using 3mer motifs across various groups according to an embodiment of the present disclosure. The entropy of HCC subjects using 3mer motifs (64 motifs in total) was found to be significantly higher (p-value<0.0001) than that of non-HCC subjects. Figure 26B shows the ROC curve using the entropy of 64 3mer motifs between non-HCC and HCC groups according to an embodiment of the present disclosure. The AUC was found to be 0.872.

[0200] As noted above, a higher entropy value indicates a higher diversity of terminal motifs. As a further illustration of the ability of embodiments to use motif diversity scores to discriminate between various cancer types and control (e.g., healthy) samples, data from a published study was used.

[0201] Figures 27A and 27B show box plots of motif diversity scores using 4mers across various groups according to an embodiment of the present disclosure. All 256 4mers were used to determine the motif diversity scores. When MDS analysis was performed using the sequencing results of plasma DNA downloaded from a published study (Song et al. 2017), an increase in plasma DNA end diversity could be generally observed among various cancer types, which may reflect the fact that different tumor cells from different anatomical sites shed their DNA into the blood circulation (Bettegowda et al. 2014). The cancers analyzed were hepatocellular carcinoma (HCC), lung cancer (LC), breast cancer (BC), gastric cancer (GC), glioblastoma multiforme (GBM), pancreatic cancer (PC), and colorectal cancer (CRC).

[0202] To further test the generalizability of the MDS changes across various cancer types, an independent cohort was further sequenced using 40 plasma DNA samples from other cancer types, including patients with colorectal cancer (n=10), lung cancer (n=10), nasopharyngeal cancer (n=10), and head and neck squamous cell carcinoma (n=10), with a median of 42 million paired end-terminal reads (range: 19-65 million). As shown in Figure 27B, the MDS values ​​of the cancer patient group (median: 0.943, range: 0.939-0.949) were significantly higher than the non-cancer control group (median: 0.941, range: 0.933-0.946, p-value < 0.0001, Wilcoxon sum rank test).

[0203] 28 shows receiver operating curves for various techniques for distinguishing healthy controls from cancer according to embodiments of the present disclosure. There were a total of 129 samples including healthy controls (n=38), HBV carriers (n=17), patients with hepatocellular carcinoma (n=34), colorectal cancer (n=10), lung cancer (n=10), nasopharyngeal carcinoma (n=10), and head and neck squamous cell carcinoma (n=10). Interestingly, the MDS-based method 2801 (AUC=0.85) appeared to have the best performance compared to other fragment metrics including fragment size 2803 (AUC=0.74, p-value=0.0040, DeLong test) (Yu et al. 2017b), end-favorable fragments 2804 (AUC=0.52, p-value<0.0001) (Jiang et al. 2018) and orientation-aware trait cell-free fragment signal, OCF, 2802 (AUC=0.68, p-value=0.0013) (Sun et al. 2019). The combined analysis 2805 identified a subject as having cancer if any one of the techniques classified the subject as having cancer.

[0204] The accuracy of the MDS analysis for distinguishing cancer and non-cancer remains relatively good for motifs of various lengths. The analysis was performed using MDS for 1mer to 5mer.

[0205] FIG. 29 shows receiver operating curves for MDS analysis using various k-mers according to an embodiment of the present disclosure. MDS values ​​estimated from 1-5mer motifs also have the power to distinguish patients with and without cancer. 1mer analysis 2901 provides 0.81 AUC. 2mer analysis 2902 provides 0.85 AUC. 3mer analysis 2903 provides 0.85 AUC. 4mer analysis 2904 provides 0.85 AUC. 5mer analysis 2905 provides 0.81 AUC.

[0206] We also investigated the effect of tumor DNA fraction on the performance of MDS-based cancer detection according to computer simulations.

[0207] Figure 30 shows the performance of MDS-based cancer detection for various tumor DNA fractions according to an embodiment of the present disclosure. As shown in Figure 30, the performance of cancer detection gradually improved as the tumor DNA fraction in plasma DNA increased. For example, the area under the curve (AUC) of ROC was only 0.52 for patients with a tumor DNA fraction of 0.1%, but for patients with a tumor DNA fraction of 3%, the AUC increased to 0.9, already approaching a maximum at a tumor fraction of 5%, and further increasing at higher concentrations.

[0208] 3. Machine Learning (SVM, Regression, and Clustering) To further investigate whether a classifier can be constructed to detect cancer patients using plasma DNA terminal motifs, a classifier was constructed to distinguish patients with cancer (n=55) and patients without cancer (n=74) using 256 plasma DNA terminal motifs, using support vector machines (SVM) and logistic regression that takes into account the magnitude and direction of each terminal motif, respectively. The SVM analysis identifies a hyperplane that best distinguishes cancer patients and non-cancer patients at 256 dimensional locations, and the training data points are the frequencies of each of the 256 motifs in 4mers. The logistic regression determines the coefficient by which each of the 256 frequencies is multiplied, and also determines the cutoff of the resulting output of the logistic function, which may be a weighted sum of the multiplied frequencies, or may receive the weighted sum as input. Such a logistic function may be a sigmoid function or other activation function, as is well known to those skilled in the art.

[0209] To minimize the problem of overfitting, a leave-one-out procedure was adopted to evaluate the performance using receiver operating characteristic (ROC) curve analysis. The leave-one-out procedure was carried out according to the following procedure. Of a sample size of N, one sample was left out as a test sample, and the remaining samples (N-1) were used to train a classifier based on SVM and logistic regression using 256 plasma DNA terminal motifs. The trained classifier was then used to determine whether the left-out sample was classified as having been taken from a subject with or without cancer. To test the trained classifier from the remaining samples, one sample was systematically left out as a test sample. Thus, prediction results for each sample could be obtained, and accuracy was calculated from the prediction results.

[0210] 31 shows receiver operating curves for MDS, SVM, and logistic regression analyses according to an embodiment of the present disclosure. A slight increase in AUC using a classifier with 256 terminal motifs (AUC=0.89 for both SVM and logistic regression) was observed compared to the MDS-based analysis (AUC=0.85).

[0211] As another machine learning technique, we used clustering based on the frequency of terminal motifs.

[0212] FIG. 32 shows a hierarchical clustering analysis of the top 10 ranked terminal motifs across different groups with different levels of cancer according to an embodiment of the present disclosure. As shown, HCC subjects (eHCC: early stage HCC 3205, iHCC: immediate stage HCC 3230, and aHCC: advanced stage HCC 3225) generally cluster together, and non-HCC (healthy control subjects, HBV: chronic hepatitis B carriers) generally cluster together. For example, the cluster on the right is early stage HCC 3205 (yellow). In the center left, there are mainly controls 3210, HBV 3215, and cirrhosis 3220. The distinct clustering patterns between HCC and non-HCC groups suggested that the terminal motifs reflect disease-related preferences in plasma DNA terminal motifs, suggesting potential diagnostic power for plasma DNA terminal motifs. In addition to connectivity-based hierarchical clustering as a statistical technique, other clustering techniques can be used, such as centroid-based clustering, distribution-based clustering, and density-based clustering.

[0213] 33A-33C show hierarchical clustering analysis using all plasma DNA molecules across different groups with different levels of cancer according to an embodiment of the present disclosure. FIG. 33A shows hierarchical clustering analysis based on 256 4-mer terminal motif frequencies. FIG. 33B shows a co-in visualization for hierarchical clustering analysis based on 256 4-mer terminal motif frequencies. Each row represents one type of terminal motif. Each column represents an individual plasma DNA sample. The gradient colors indicate the frequency of the terminal motifs, with red representing the highest frequency and green representing the lowest frequency. FIG. 33C shows a principal component analysis (PCA) of HCC and non-HCC subjects using terminal motifs. The principal components are linear combinations of the 256 motifs that provide the greatest variance, e.g., a weighted sum of the frequencies is obtained.

[0214] Since HCC and non-HCC subjects seem to form two different clusters, the terminal motifs derived from all plasma DNA molecules are important indicators for distinguishing HCC from non-HCC subjects. Figures 33A and 33B show that HCC subjects 3305 (red) tend to be clustered into one group, and non-HCC subjects 3310 (blue) tend to be clustered into another group. In Figure 33C, PCA analysis also showed that HCC and non-HCC subjects tend to be clustered into two different groups. PC1 and PC2 correspond to different linear combinations (such as weighted averages) of relative frequencies, which may represent a certain histogram pattern of relative frequencies. Figure 33C shows that linear combinations (or other transformations) can be performed before performing clustering or using cutoff values ​​or cutoff planes. Thus, the transformed relative frequencies can be used to determine aggregate values.

[0215] Figure 34 shows a hierarchical clustering analysis based on 3mer motifs using all plasma DNA molecules across different groups with different levels of cancer according to an embodiment of the present disclosure. For ease of illustration, only the top part of the heatmap is shown. As shown, HCC subjects (eHCC: early stage HCC 3405, iHCC: immediate stage HCC 3430, and aHCC: advanced stage HCC 3425) generally cluster together, and non-HCC (healthy control subjects 3410, HBV 3415: chronic hepatitis B carriers, and cirrhosis 3420) generally cluster together.

[0216] Based on these findings, machine learning (e.g., deep learning) models can be used to train cancer classifiers using 256-dimensional vectors containing plasma DNA terminal motifs, including, but not limited to, support vector machines (SVMs), decision trees, naive Bayes classification, logistic regression, clustering algorithms, PCA, singular value decomposition (SVD), t-distributed stochastic neighbor embedding (tSNE), artificial neural networks, and ensemble methods that classify new data points by configuring a set of classifiers and weighting their predictions. Once a cancer classifier is trained based on the "256-dimensional vector-based matrix" containing a set of cancer and non-cancer patients, it can predict the probability of cancer for new patients.

[0217] In such uses of machine learning algorithms, the aggregate value may correspond to a probability or distance (e.g., when using SVM) that may be compared to a reference value. In other embodiments, the aggregate value may correspond to an early output in a model (e.g., an early layer of a neural network) that is compared to a cutoff between two classifications or to a representative value of a given classification.

[0218] B. Immune disease monitoring Figure 35A shows an entropy analysis using all plasma DNA molecules between healthy control subjects and SLE patients according to an embodiment of the present disclosure. Figure 35B shows a hierarchical clustering analysis using all plasma DNA molecules between healthy control subjects and SLE patients according to an embodiment of the present disclosure.

[0219] Global abnormality analysis of plasma DNA terminal motifs, including entropy (Figure 35A, p-value: 0.00014) and clustering analysis (Figure 35B), showed that SLE patients could be distinguished from healthy control subjects. For example, the entropy increases for subjects with SLE (Figure 35A). And two clusters are generally formed on the left (SLE 3510) and right (control / normal 3505). Thus, autoimmune diseases change plasma DNA fragmentation patterns, thereby indicating the discriminatory power of plasma DNA terminal motifs between SLE and control subjects.

[0220] FIG. 36 shows an entropy analysis using plasma DNA molecules with 10 selected terminal motifs between healthy control subjects and SLE patients according to an embodiment of the present disclosure. The top 10 motifs with the highest relative frequency for control subjects were used. As with other phenotypes, the set of motifs may affect whether the SLE entropy is high or low. Considering that the 10 motifs were selected as the ones with the highest values ​​for controls, the entropy will be high because the values ​​are similar to each other (i.e., because of ranking). Also, the SLE entropy will be lower with more variation, for example, because they are not ranked for SLE subjects. If the top 10 motifs were selected using SLE samples, the inverse relationship may exist. Thus, the level of autoimmune disease (e.g., SLE) may be determined using the aggregate value of the relative frequency.

[0221] C. Synergy analysis of terminal motifs and conventional metrics We tested whether combined analysis of plasma DNA terminal motifs and other metrics (copy number aberrations (CNAs), hypomethylation, and hypermethylation) would improve the performance of noninvasive cancer detection. For example, decision tree-based classification could be used for the combined analysis.

[0222] FIG. 37 shows the ROC curves of combined analysis including terminal motifs and copy number or methylation for HCC and non-HCC subjects according to an embodiment of the present disclosure. The terminal motif analysis uses the motif diversity score determined using all 356 motifs of 4mers. The combined analysis identifies cancer when either analysis results in a classification of cancer. The combined analysis of terminal motifs and methylation analysis (AUC: 0.94) or terminal motifs and CNA analysis (AUC: 0.93) outperformed the analysis using only terminal motifs (AUC: 0.86). The methylation analysis used the number of 1 Mb bins with hypomethylation (defined as methylation density z-score <-3) above the number of normal controls as the cutoff number of abnormal bins to distinguish cancer from non-cancer. The CNA analysis used the number of 1 Mb bins with z-score > 3 or < -3 as the cutoff number of abnormal bins to distinguish cancer from non-cancer. Further details on methylation analysis can be found in US Patent Publication 2014 / 0080715 and on CNA analysis can be found in US Patent Publication US2013 / 0040824.

[0223] An example of decision tree-based classification is described. For example, a random forest algorithm can be used to estimate cutoffs for each metric, such as CNA, hypomethylation, hypermethylation, size (e.g., as described in U.S. Patent Publication 2013 / 0237431), terminal motifs, and fragmentation patterns (e.g., as described in U.S. Patent Publications 2017 / 0024513 and 2019 / 0341127 and U.S. Patent Application 16 / 519,912). Each metric has a specific cutoff. Take one metric (hypomethylation) as an example, a case can be classified as cancer or non-cancer depending on whether the metric is below or above the cutoff. A metric represents a node in the decision tree. For example, after a sample has traveled through all the nodes in the entire tree, a majority of votes (e.g., the number of nodes indicating cancer is greater than the number of nodes indicating non-cancer) can provide the final classification.

[0224] D. Examples of alternative methods for defining terminal motifs in plasma DNA To demonstrate the feasibility of using an alternative method to define terminal motifs in plasma DNA, technique 160 in Figure 1 was employed to analyze HCC and non-HCC subjects, including 20 healthy control subjects (Control), 22 chronic hepatitis B carriers (HBV), 12 cirrhotic subjects (Cirr), 24 early stage HCC (eHCC), 11 immediate stage HCC (iHCC), and 7 advanced stage HCC (aHCC).

[0225] FIG. 38A shows an entropy analysis based on 4mers jointly constructed from the ends of sequenced plasma DNA fragments and their adjacent genomic sequences in HCC and non-HCC subjects according to an embodiment of the present disclosure. The entropy was determined using all 256 terminal motifs. Similar to the analysis using the technique 140 of FIG. 1 to define motifs, the entropy of HCC subjects is different from non-cancer subjects. Also, advanced HCC is significantly different from eHCC and iHCC. FIG. 38B shows a clustering analysis based on 4mers jointly constructed from the ends of sequenced plasma DNA fragments and their adjacent genomic sequences in HCC subjects 3810 and non-HCC subjects 3805 according to an embodiment of the present disclosure.

[0226] Figure 39 shows a ROC comparison for techniques 140 and 160 of Figure 1 used to define terminal motifs of plasma DNA according to an embodiment of the present disclosure. The same subjects as in Figure 38A were used and classified using entropy analysis using 4mers. Method (i) corresponds to technique 140 and method (ii) corresponds to technique 160. A slightly inferior performance (AUC: 0.815 vs. 0.856) was observed using technique 160 of Figure 1 compared to technique 140 in Figure 1.

[0227] E. Filtering to improve discrimination Certain criteria may be used to filter out certain DNA fragments (other than terminal motifs), e.g., to provide high accuracy of sensitivity and specificity. As an example, terminal motif analysis may be limited to DNA fragments originating from open chromatin regions of a particular tissue, e.g., as determined by reads that align completely or partially within one of multiple open chromatin regions. For example, any read that has at least one nucleotide that overlaps with an open chromatin region may be defined as a read within an open chromatin region. A typical open chromatin region is about 300 bp according to DNase I hypersensitive sites. The size of the open chromatin region may vary depending on the technology used to define the open chromatin region, e.g., ATAC-seq (Assay for Transposase Accessible Chromatin Sequencing) vs. DNaseI-Seq.

[0228] As another example, DNA fragments of a particular size can be selected to perform terminal motif analysis on, which, as shown below, increases the resolution of the tabulation of relative frequencies of terminal motifs, thereby improving precision.

[0229] Further examples may use the methylation characteristics of DNA fragments. Fetal and tumor DNA are generally hypomethylated. An embodiment may determine the methylation metric (e.g., density) of a DNA fragment (e.g., as a percentage or absolute number of site(s) that are methylated on the DNA fragment). Also, DNA fragments may be selected for use in terminal motif analysis based on the measured methylation density. For example, DNA fragments may be used only if the methylation density is above a threshold value.

[0230] Whether a DNA fragment contains sequence variation (e.g., base substitutions, insertions, or deletions) compared to a reference genome can also be used for filtering.

[0231] Various filtering criteria may be used in combination. For example, each criterion may need to be met, or at least a certain number of criteria may need to be met. In another implementation, the probability that a fragment corresponds to clinically relevant DNA (e.g., fetal, tumor, or transplant) may be determined, and a threshold imposed probability must be met before the DNA fragment is used in terminal motif analysis. As a further example, the contribution of a DNA fragment to a frequency counter of a particular terminal motif may be weighted based on probability (e.g., instead of adding one, add the probability of having a value less than 1). Thus, DNA fragments with a particular terminal motif will be weighted higher and / or have a higher probability. Such enrichment is further described below.

[0232] 1. Terminal motifs across tissue-specific chromatin regions Because different tissues have preferred fragmentation patterns during apoptosis (Chan et al., Proc Natl Acad Sci USA. 2016;113:E8159-8168; Jiang et al., Proc Natl Acad Sci USA. 2018;doi:10.1073 / pnas.1814616115), we further reasoned that the selection of specific genomic regions for plasma DNA terminal motif analysis would further improve the discriminatory power in classifying diseased and control subjects. Taking the detection of HCC patients as an example, open chromatin regions from blood and liver were used.

[0233] Figure 40 shows a comparison of accuracy showing that tissue-specific open chromatin regions improve the discrimination of plasma DNA terminal motifs in HCC and non-cancer patients according to an embodiment of the present disclosure. Analysis was performed on the entropy of all 256 motifs using 4mers and the combined frequency of the top 10 motifs. For liver open chromatin results, sequence reads were retained (i.e., not filtered out) if the reads had at least one nucleotide overlapping with one of the liver open chromatin regions.

[0234] The power of terminal motifs derived from plasma DNA molecules that overlap with open chromatin regions in the liver yields the best performance with an AUC of 0.918 using the combined frequency of the top 10 ranked motifs. In contrast, the discriminatory power of terminal motifs derived from plasma DNA molecules of all 256 motifs without any selection had the lowest AUC of 0.855.

[0235] Thus, if a particular tissue is being screened for cancer, DNA fragments from the open chromatin of that particular tissue (or at least if the terminal sequences are in the open chromatin regions) can be used to perform the analysis, whereas DNA fragments that are not in these identified regions are not used. The liver was used here because the cancer was HCC. The location of the DNA fragments can be determined by aligning the sequence reads to a reference genome, which may identify open chromatin regions from the literature or a database.

[0236] 2. Size-band-based Terminal Motif Analysis The frequency of certain terminal motifs has been shown to vary depending on the size range (size band) being analyzed; for example, the percentage of CCCAs shows this behavior. This means that size band-based terminal motif analysis may affect performance in using plasma DNA terminal motifs to distinguish cancer patients from non-cancer subjects. To account for this possibility, a series of size ranges, including but not limited to 50-80 bp, 81-110 bp, 111-140 bp, 141-170 bp, 171-200 bp, and 201-230 bp, are tested to investigate how the size band being analyzed affects the overall diagnostic performance.

[0237] FIG. 41 illustrates a size band based plasma DNA terminal motif analysis according to an embodiment of the present disclosure. Classification using motif diversity scores (entropy) is determined using 256 motifs of 4-mers. Although various ranges are listed in FIG. 41, other ranges may be used. 50-80 analysis 4101 provides 0.826 AUC. 81-110 analysis 4102 provides 0.537 AUC. 111-140 analysis 4103 provides 0.551 AUC. 141-170 analysis 4104 provides 0.716 AUC. 171-200 analysis 4105 provides 0.769 AUC. 201-230 analysis 4106 provides 0.756 AUC.

[0238] Such size ranges may be used for techniques that enrich for clinically relevant DNA. For example, selecting DNA molecules between 50 and 80 bases would enrich the sample for tumor DNA. Instead of a single size range, multiple disjoint size ranges may be used. Such enrichment may be the reason for better AUC in the 50-80 base vs. 81-110 base size ranges.

[0239] Terminal motifs derived from plasma DNA molecules within the 50-80 bp range appeared to give the highest discriminatory power for detecting HCC from non-HCC subjects (AUC: 0.83). Thus, embodiments may filter DNA fragments to select fragments in a particular size range, and then use the selected DNA fragments (reads) to determine relative frequency and subsequent manipulations. By way of example, size filtering may be performed via physical separation or by using sequence reads to determine size (e.g., length if the entire fragment is sequenced, or by aligning paired ends to a reference). Examples of physical enrichment of short DNA include band excision in gel electrophoresis, collection of eluates at specific retention times in capillary electrophoresis, post liquid chromatography, or by microfluidics.

[0240] F. Classification of the Level of Pathology 42 is a flow chart illustrating a method 4200 for classifying a level of pathology in a biological sample of a subject, the biological sample comprising cell-free DNA, according to an embodiment of the present disclosure. Aspects of method 4200 may be performed in a manner similar to method 1900 of FIG. 19 and method 2000 of FIG. 20.

[0241] At block 4210, a plurality of cell-free DNA fragments from the biological sample are analyzed to obtain sequence reads. The sequence reads include terminal sequences corresponding to ends of the plurality of cell-free DNA fragments. Block 4210 may be performed in a manner similar to block 1910 of FIG. 19.

[0242] For each of the plurality of cell-free DNA fragments, a sequence motif is determined for each of the one or more end sequences of the cell-free DNA fragments at block 4220. Block 4220 may be performed in a manner similar to block 1920 of FIG.

[0243] At block 4230, a relative frequency of a set of one or more sequence motifs corresponding to end sequences of the plurality of cell-free DNA fragments is determined. The relative frequency of the sequence motifs may provide a percentage of the plurality of cell-free DNA fragments having end sequences corresponding to the sequence motifs. Block 4230 may be performed in a manner similar to block 1930 of FIG. 19. For example, the set of one or more sequence motifs may include N base positions. The set of one or more sequence motifs may include all combinations of N bases. N may be an integer greater than or equal to 3, and other integers.

[0244] As another example, the set of one or more sequence motifs can be the top M sequence motifs with the greatest difference between two types of DNA determined in one or more reference samples, for example, all of the motifs showing the greatest positive difference (e.g., top 10 or other number) or all of the motifs with the greatest negative difference. M can be an integer equal to or greater than 1. For methods 1900 and 2000, the two types of DNA can be clinically relevant DNA and other DNA. For method 4200, the two types of DNA can be from two reference samples with different classifications for the level of pathology. As a further example, the set of one or more sequence motifs can be the top M most frequent sequence motifs occurring in one or more reference samples, for example, as shown in FIG. 22, the reference sample is a non-cancer sample such as an HBV sample.

[0245] At block 4240, an aggregate value of the relative frequencies of the set of one or more sequence motifs is determined. Block 4240 may be performed in a manner similar to block 1940 of Figure 19. Examples of aggregate values ​​are described throughout this disclosure and include entropy, composite frequency, difference (e.g., distance) from a reference pattern of relative frequency that may be implemented in clustering or using SVM, a value (e.g., probability) determined from the difference or output in a machine learning model (e.g., an intermediate or final layer in a neural network) compared to a cutoff between two classifications or compared to a representative value of a given classification.

[0246] When the set of one or more sequence motifs includes a plurality of sequence motifs, the summary value can include the sum of the relative frequencies of the set. The sum can be a weighted sum. For example, the summary value can include an entropy term that includes the sum of the terms that include the weighted sum. Each term can include the relative frequency multiplied by the logarithm of the relative frequency. The summary value can correspond to the variance of the relative frequencies.

[0247] In another example, the aggregate value includes a final or intermediate output of a machine learning model. In various implementations, the machine learning model uses clustering, a support vector machine, or logistic regression.

[0248] At block 4250, a classification of a level of pathology may be determined for the subject based on the comparison of the aggregate value with the reference value. As an example, the pathology may be cancer or an autoimmune disorder. As an example, the level may be no cancer, early stage, intermediate stage, or advanced stage. The classification may then select one of the levels. Thus, the classification may be determined from multiple levels of cancer, including multiple stages of cancer. As an example, the cancer may be hepatocellular carcinoma, lung cancer, breast cancer, gastric cancer, glioblastoma multiforme, pancreatic cancer, colorectal cancer, nasopharyngeal carcinoma, and head and neck squamous cell carcinoma. As an example, the autoimmune disorder may be systemic lupus erythematosus.

[0249] In a further example, the level of pathology corresponds to a fractional concentration of clinically relevant DNA associated with the pathology. For example, the level of pathology can be cancer and the clinically relevant DNA can be tumor DNA. The reference value can be a calibration value determined from a calibration sample as described for method 1900.

[0250] In some embodiments, cell-free DNA is filtered to identify a number of cell-free DNA fragments. Examples of filtering are described in the above section. For example, filtering can be based on methylation (density or whether a particular site is methylated), size, or the region from which the DNA fragment originates. Cell-free DNA can be filtered for DNA fragments from open chromatin regions of a particular tissue.

[0251] IV. Concentration Selection of DNA fragments from a particular tissue that exhibit a particular set of terminal motifs can be used to enrich a sample for DNA from that particular tissue. Thus, embodiments can enrich a sample for clinically relevant DNA. For example, only DNA fragments with particular terminal sequences can be sequenced, amplified, and / or captured using the assay. As another example, filtering of sequence reads can be performed in a manner similar to that described, for example, in section III.E.

[0252] A. Physical concentration Physical enrichment can be carried out in various ways, for example, through targeted sequencing or PCR, such as by using specific primers or adapters. If a specific terminal motif of the terminal sequence is detected, an adapter can be added to the end of the fragment. Then, when sequencing is carried out, only the DNA fragments with adapters are sequenced (or at least mainly sequenced), thereby providing targeted sequencing.

[0253] As another example, primers that hybridize to a set of specific terminal motifs can be used. These primers can then be used to perform sequencing or amplification. Capture probes that correspond to specific terminal motifs can also be used to capture DNA molecules with those terminal motifs for further analysis. Some embodiments can link short oligonucleotides to the ends of plasma DNA molecules. Then, the probes can be designed to only recognize sequences that are partly terminal motifs and partly linked oligonucleotides.

[0254] Some embodiments may use CRISPR-based diagnostic techniques, for example using guide RNAs to identify sites corresponding to preferred terminal motifs in clinically relevant DNA, and then using nucleases to cleave the DNA fragments, as can be done using Cas-9 or Cas-12. For example, an adapter may be used to recognize the terminal motif, and CRISPR / Cas9 or Cas-12 may be used to cleave the terminal motif / adapter hybrid and create a universally recognizable end for further enrichment of molecules at the desired end.

[0255] 43 is a flow chart illustrating a method 4300 of enriching a biological sample for clinically relevant DNA according to an embodiment of the present disclosure. The biological sample includes clinically relevant DNA molecules and other DNA molecules that are acellular. The method 4300 may perform the enrichment using a specific assay.

[0256] At block 4310, a plurality of cell-free DNA fragments are received from a biological sample. Clinically relevant DNA fragments (e.g., fetal or tumor) have end sequences that include sequence motifs that occur at a higher relative frequency than other DNA (e.g., maternal DNA, healthy DNA, or blood cells). As an example, data from Figures 3 and 13 may be used. Thus, sequence motifs may be used to enrich for clinically relevant DNA.

[0257] At block 4320, the plurality of cell-free DNA fragments are subjected to one or more probe molecules that detect sequence motifs in the terminal sequences of the plurality of cell-free DNA fragments. Such use of probe molecules can result in obtaining detected DNA fragments. In one example, the one or more probe molecules can include one or more enzymes that interrogate the plurality of cell-free DNA fragments and add new sequences that are used to amplify the detected DNA fragments. In another example, the one or more probe molecules can be attached to a surface to detect sequence motifs in the terminal sequences by hybridization.

[0258] At block 4330, the detected DNA fragments are used to enrich the biological sample for clinically relevant DNA fragments. As an example, enriching the biological sample for clinically relevant DNA fragments using the detected DNA fragments may include amplifying the detected DNA fragments. As another example, the detected DNA fragments may be captured and the undetected DNA fragments may be discarded.

[0259] B. In silico enrichment In silico enrichment may use various criteria to select or discard specific DNA fragments. Such criteria include terminal motifs, open chromatin regions, size, sequence diversity, methylation, and other epigenetic features. Epigenetic features include all modifications of the genome that do not involve changes in the DNA sequence. Criteria may specify cutoffs that require certain properties, such as, for example, a certain size range, a methylation metric above or below a certain amount, a combination of methylation states of two or more CpG sites (e.g., a methylation haplotype (Guo et al, Nat Genet. 2017;49:635-42)), or have a combined probability above a threshold. Such enrichment may also include weighting DNA fragments based on such probability.

[0260] For example, enriched samples can be used to classify pathology (as described above), as well as to identify tumor or fetal mutations, or for tag counting to detect amplification / deletion of chromosomes or chromosomal regions.For example, if a particular terminal motif or set of terminal motifs is associated with liver cancer (i.e., higher relative frequency than non-cancer or other cancers), an embodiment for performing cancer screening can weight such DNA fragments higher than DNA fragments that do not have this preferred one or this preferred set of terminal motifs.

[0261] 44 is a flow chart illustrating a method 4400 of enriching a biological sample for clinically relevant DNA according to an embodiment of the present disclosure. The biological sample includes clinically relevant DNA molecules and other DNA molecules that are acellular. The method 4400 may perform the enrichment using certain criteria of sequence reads.

[0262] At block 4410, a plurality of cell-free DNA fragments from the biological sample are analyzed to obtain sequence reads. The sequence reads include end sequences corresponding to ends of the plurality of cell-free DNA fragments. Block 4410 may be performed in a manner similar to block 1910 of FIG. 19.

[0263] For each of the plurality of cell-free DNA fragments, a sequence motif is determined for each of the one or more end sequences of the cell-free DNA fragments at block 4420. Block 4420 may be performed in a manner similar to block 1920 of FIG.

[0264] At block 4430, a set of one or more sequence motifs that occur in the clinically relevant DNA at a higher relative frequency than in other DNA is identified. The set of sequence motif(s) may be identified by genotypic or phenotypic techniques described herein. Calibration or reference samples may be used to rank and select sequence motifs selective for clinically relevant DNA.

[0265] At block 4440, a group of sequence reads having a set of one or more sequence motifs in the end sequences is identified. This may be considered the first stage of filtering.

[0266] At block 4450, sequence reads having a likelihood corresponding to clinically relevant DNA that exceeds a threshold may be stored. The likelihood may be determined using the set of terminal motif(s). For example, for each sequence read of a group of sequence reads, the likelihood of the sequence read corresponding to clinically relevant DNA may be determined based on the terminal sequence of the sequence read containing one or more sequence motifs of the set of sequence motifs. The likelihood may be compared to a threshold. As an example, the threshold may be determined empirically. For example, for a sample in which the concentration of clinically relevant DNA may be measured for a group of sequence reads, various thresholds may be tested. An optimal threshold may maximize the concentration while maintaining a certain percentage of the total number of sequence reads. The threshold may be determined by one or more given percentiles (5, 10, 90, or 95) of the concentration of one or more terminal motifs present in healthy controls or a control group without disease but exposed to similar etiological risk factors. The threshold may be a regression or probability score.

[0267] If the likelihood exceeds a threshold, the sequence read may be stored in memory (e.g., a file, table, or other data structure), thereby obtaining a stored sequence read. Sequence reads with a likelihood below a threshold may be discarded or not stored in the memory location of the retained reads, or a field in the database may include a flag indicating a low threshold of reads so that later analysis may exclude such reads. By way of example, the likelihood may be determined using various techniques, such as odds ratios, z-scores, or probability distributions.

[0268] At block 4460, the stored sequence reads may be analyzed to determine a characteristic of the clinically relevant DNA biological sample as described in other flowcharts, e.g., as described herein. Methods 1900, 2000, and 4200 are such examples. For example, the characteristic of the clinically relevant DNA biological sample may be a fractional concentration of clinically relevant DNA. As another example, the characteristic may be a level of pathology of the subject from which the biological sample was obtained, the level of pathology being associated with the clinically relevant DNA. As another example, the characteristic may be a gestational age of the fetus of the pregnant woman from whom the biological sample was obtained.

[0269] Other criteria can be used to determine likelihood.The size of a plurality of cell-free DNA fragments can be measured using sequence reads.The likelihood that a particular sequence read corresponds to clinically relevant DNA can be further based on the size of the cell-free DNA fragment that corresponds to a particular sequence read.

[0270] Methylation may also be used. Thus, embodiments may measure one or more methylation states at one or more sites of cell-free DNA fragments that correspond to a particular sequence read. The likelihood that a particular sequence read corresponds to clinically relevant DNA may further be based on one or more methylation states. As a further example, whether a read is within an identified set of open chromatin regions may be used as a filter.

[0271] FIG. 45 shows an exemplary plot showing the increase in fetal DNA fraction using CCCA terminal motifs according to an embodiment of the present disclosure. The vertical axis is the fetal DNA fraction for the samples tested. The two sets of data are for (1) all fragments overlapping with the informative SNP (i.e., fragments with fetal-specific alleles) and (2) fragments with CCCA terminal motifs and overlapping with the informative SNP. Thus, the data on the left provides the actual fetal DNA fraction in the whole sample, and the data on the right provides the data of the in silico enriched sample. In this example, it can be determined that the likelihood is above a threshold if the terminal motif is CCCA. More motifs can be used in a similar manner, e.g., as a group indicating the likelihood is above a threshold.

[0272] The median relative increase in fetal DNA fraction was 3.2% (IQR: 1.3-6.4%). The relative increase in fetal DNA fraction was (ba) / a * 100, where a is the original fetal DNA fraction calculated by all fragments overlapping with informative SNPs for which the mother is homozygous and the fetus is heterozygous, and b is the fetal DNA fraction calculated by fragments tagged by CCCA motifs that are abundant in fetal DNA molecules.

[0273] For any of the methods described herein, the sequence motif for each of one or more end sequences of the cell-free DNA fragments can be performed using a reference genome (e.g., via technique 160 of FIG. 1 ). Such a technique involves aligning one or more sequence reads corresponding to the cell-free DNA fragments to the reference genome, identifying one or more bases in the reference genome adjacent to the end sequence, and using the end sequence and the one or more bases to determine the sequence motif.

[0274] V. Example System FIG. 46 illustrates a measurement system 4600, according to an embodiment of the invention. The system shown includes a sample 4605, such as cell-free DNA molecules in a sample holder 4610, which may contact an assay 4608 to provide a signal of a physical property 4615. An example of a sample holder may be a flow cell, which includes probes and / or primers of an assay, or a tube through which droplets (along with droplets containing an assay) travel. The physical property 4615 (e.g., fluorescence intensity, voltage, or current) from the sample is detected by a detector 4620. The detector 4620 may measure at intervals (e.g., periodic intervals) to obtain data points that constitute a data signal. In one embodiment, an analog-to-digital converter converts the analog signal from the detector to digital form multiple times. The sample holder 4610 and the detector 4620 may form an assay device, e.g., a sequencing device that performs sequencing according to embodiments described herein. The data signal 4625 is transmitted from the detector 4620 to a logic system 4630. The data signal 4625 may be stored in a local memory 4635 , an external memory 4640 , or a storage device 4645 .

[0275] Logic system 4630 may be or include a computer system, ASIC, microprocessor, etc. It may also include or be coupled to a display (e.g., monitor, LED display, etc.), and user input devices (e.g., mouse, keyboard, buttons, etc.). Logic system 4630 and other components may be part of a stand-alone or networked computer system, or may be directly attached to or incorporated into a device (e.g., a sequencing device) including detector 4620 and / or sample holder 4610. Logic system 4630 may also include software executing in processor 4650. Logic system 4630 may include a computer readable medium storing instructions for controlling system 4600 to perform any of the methods described herein. For example, logic system 4630 may provide commands to a system including sample holder 4610 such that sequencing or other physical operations are performed. Such physical operations may be performed in a particular order, such as, for example, reagents are added and removed in a particular order. Such physical manipulations can be performed by a robotic system, including, for example, a robotic arm, such that the system can be used to obtain samples and perform assays.

[0276] Any of the computer systems referred to herein may utilize any suitable number of subsystems. An example of such subsystems is shown in FIG. 47 for computer system 10. In some embodiments, a computer system includes a single computer device, and the subsystems may be components of the computer device. In other embodiments, a computer system may include multiple computer devices, each of which is a subsystem and has internal components. Computer systems may include desktop and laptop computers, tablets, mobile phones, and other mobile devices.

[0277] The subsystems shown in FIG. 47 are interconnected via a system bus 75. Additional subsystems are shown, such as a printer 74, a keyboard 78, a storage device(s) 79, a monitor 76 (e.g., a display screen such as LED) connected to a display adapter 82, and others. Peripherals and input / output (I / O) devices coupled to the I / O controller 71 may be connected to the computer system by any number of means known in the art, such as an input / output (I / O) port 77 (e.g., USB, FireWire). For example, the I / O port 77 or an external interface 81 (e.g., Ethernet, Wi-Fi, etc.) may be used to connect the computer system 10 to a wide area network such as the Internet, a mouse input device, or a scanner. The interconnection via the system bus 75 allows the central processor 73 to communicate with each subsystem and control the execution of instructions from the system memory 72 or storage device(s) 79 (e.g., a fixed disk such as a hard drive or optical disk), and the exchange of information between the subsystems. The system memory 72 and / or storage device(s) 79 may embody a computer-readable medium. Another subsystem is a data collection device 85, such as a camera, microphone, and accelerometer, and the like. Any of the data referred to herein may be output from one component to another or to a user.

[0278] A computer system may include multiple identical components or subsystems connected together, for example, by an external interface 81, by an internal interface, or through a storage device that may be connected or removed from one component to another. In some embodiments, computer systems, subsystems, or devices may communicate over a network. In such an example, one computer may be considered a client and another computer a server, each of which may be part of the same computer system. The client and server may each include multiple systems, subsystems, or components.

[0279] Aspects of the embodiments may be implemented using hardware circuitry in the form of control logic (e.g., application specific integrated circuits or field programmable gate arrays) and / or using computer software with general-purpose programmable processors in a modular or integrated manner. As used herein, a processor may include a single-core processor, a multi-core processor on the same integrated chip, or multiple processing units on a single circuit board or networked, as well as dedicated hardware. Based on the present disclosure and teachings provided herein, those skilled in the art will recognize and understand other ways and / or methods for implementing embodiments of the present invention using hardware and combinations of hardware and software.

[0280] Any of the software components or functions described in this application may be implemented as software code executed by a processing device using any suitable computer language, such as, for example, Java, C, C++, C#, Objective-C, Swift, or a scripting language, such as, for example, Perl or Python, using conventional or object-oriented techniques. The software code may be stored as a series of instructions or commands on a computer-readable medium for storage and / or transmission. Suitable non-transitory computer-readable media may include random access memory (RAM), read-only memory (ROM), magnetic media (such as a hard drive or floppy disk), or optical media (such as a compact disk (CD) or DVD (digital versatile disk)), or Blu-ray disks and flash memories, and the like. The computer-readable medium may be any combination of such storage or transmission devices.

[0281] Such a program may also be coded and transmitted using a carrier signal adapted for transmission over wired, optical, and / or wireless networks following various protocols, including the Internet. Thus, a computer-readable medium may be created using a data signal coded with such a program. A computer-readable medium coded with a program code may be packaged with a compatible device or provided separately from other devices (e.g., via Internet download). Any such computer-readable medium may be present on or within a single computer product (e.g., a hard drive, a CD, or an entire computer system), or may be present on or within different computer products in a system or network. A computer system may include a monitor, printer, or other suitable display for providing any of the results described herein to a user.

[0282] Any of the methods described herein may be implemented in whole or in part using a computer system including one or more processors that can be configured to perform the steps. Thus, embodiments may be directed to a computer system configured to perform the steps of any of the methods described herein, with potentially different components performing each step or each group of steps. Although presented as numbered steps, steps of the methods herein may be performed simultaneously or at different times, or in different orders. In addition, portions of these steps may be combined with portions of other steps from other methods. Also, all or portions of a step may be optional. In addition, any of the steps of any of the methods may be implemented using modules, units, circuits, or other means of systems for performing these steps.

[0283] The specific details of particular embodiments may be combined in any suitable manner without departing from the spirit and scope of the embodiments of the invention, however, other embodiments of the invention may be directed to specific embodiments relating to each individual aspect, or specific combinations of these individual aspects.

[0284] The above description of exemplary embodiments of the present disclosure has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosure to the precise form described, and many modifications and variations are possible in light of the above teachings.

[0285] The descriptions "a," "an," or "the" are intended to mean "one or more," unless specifically indicated to the contrary. The use of "or" is intended to mean "including, or," rather than "except, or," unless specifically indicated to the contrary. A reference to a "first" element does not necessarily require that a second element is provided. Further, a reference to a "first" or "second" element does not limit the referenced element to a particular location unless expressly stated. The term "based on" is intended to mean "based at least in part on."

[0286] All patents, patent applications, publications, and specifications mentioned herein are incorporated by reference in their entirety for all purposes. None is admitted to be prior art.

Claims

1. 1. A method of classifying a level of pathology in a biological sample of a subject, the biological sample comprising cell-free DNA, the method comprising: analyzing a plurality of cell-free DNA fragments from the biological sample to obtain sequence reads, the sequence reads including terminal sequences corresponding to terminals of the plurality of cell-free DNA fragments; determining, for each of the plurality of cell-free DNA fragments, a sequence motif for each of one or more terminal sequences of the cell-free DNA fragments; determining a relative frequency of a set of one or more sequence motifs that correspond to the end sequences of the plurality of cell-free DNA fragments, wherein the relative frequency of a sequence motif provides a proportion of the plurality of cell-free DNA fragments having end sequences that correspond to the sequence motif; determining a summary of the relative frequencies of the set of one or more sequence motifs; and determining a classification of a level of pathology for the subject based on a comparison of the aggregated value to a reference value.

2. 10. The method of claim 1, further comprising filtering the cell-free DNA to identify the plurality of cell-free DNA fragments.

3. The method of claim 2 , wherein the filtering is based on the size or region from which the DNA fragments originate.

4. 4. The method of claim 3, wherein the cell-free DNA is filtered for DNA fragments derived from open chromatin regions of a particular tissue.

5. The method of claim 1 , wherein the pathology is cancer.

6. 6. The method of claim 5, wherein the cancer is hepatocellular carcinoma, lung cancer, breast cancer, gastric cancer, glioblastoma multiforme, pancreatic cancer, colorectal cancer, nasopharyngeal carcinoma, and head and neck squamous cell carcinoma.

7. 6. The method of claim 5, wherein the classification is determined from multiple levels of cancer, including multiple stages of cancer.

8. The method of claim 1 , wherein the pathology is an autoimmune disorder.

9. The method of claim 8, wherein the autoimmune disorder is systemic lupus erythematosus.

10. The method of claim 1 , wherein the level of pathology corresponds to a fractional concentration of clinically relevant DNA associated with the pathology.

11. 1. A method of estimating the fractional concentration of clinically relevant DNA in a biological sample of a subject, said biological sample comprising said clinically relevant DNA and other DNA that is acellular, said method comprising: analyzing a plurality of cell-free DNA fragments from the biological sample to obtain sequence reads, the sequence reads including terminal sequences corresponding to terminals of the plurality of cell-free DNA fragments; determining, for each of the plurality of cell-free DNA fragments, a sequence motif for each of one or more terminal sequences of the cell-free DNA fragments; determining a relative frequency of a set of one or more sequence motifs that correspond to the end sequences of the plurality of cell-free DNA fragments, wherein the relative frequency of a sequence motif provides a proportion of the plurality of cell-free DNA fragments having end sequences that correspond to the sequence motif; determining a summary of the relative frequencies of the set of one or more sequence motifs; and determining a classification of the fractional concentration of clinically relevant DNA in the biological sample by comparing the aggregated value to one or more calibration values ​​determined from one or more calibration samples having known fractional concentrations of clinically relevant DNA.

12. 12. The method of claim 11, wherein the clinically relevant DNA is selected from the group consisting of fetal DNA, tumor DNA, DNA from a transplanted organ, and a specific tissue type.

13. The method of claim 11 , wherein the clinically relevant DNA is of a specific tissue type.

14. The method of claim 13, wherein the specific tissue type is hepatic or hematopoietic.

15. 12. The method of claim 11, wherein the subject is a pregnant female and the clinically relevant DNA is placental tissue.

16. The method of claim 11 , wherein the clinically relevant DNA is tumor DNA derived from an organ having cancer.

17. 12. The method of claim 11, wherein the one or more calibration values ​​are a plurality of calibration values ​​of a calibration function determined using fractional concentrations of clinically relevant DNA of a plurality of calibration samples.

18. 12. The method of claim 11, wherein the one or more calibration values ​​correspond to one or more summary values ​​of the relative frequencies of the set of one or more sequence motifs measured using cell-free DNA fragments in the one or more calibration samples.

19. For each calibration sample of the one or more calibration samples, determining the fractional concentration of clinically relevant DNA in the calibration sample; 12. The method of claim 11, further comprising determining the aggregate value of the relative frequencies of the set of one or more sequence motifs by analyzing cell-free DNA fragments from the calibration sample as part of obtaining calibration data points, thereby determining one or more aggregate values, wherein each calibration data point specifies the measured fractional concentration of clinically relevant DNA in the calibration sample and the aggregate value determined for the calibration sample, and wherein the one or more calibration values ​​are or are determined using the one or more aggregate values.

20. 20. The method of claim 19, wherein the measuring of the fractional concentration of clinically relevant DNA in the calibration sample is performed using alleles specific for the clinically relevant DNA.

21. 1. A method for determining the gestational age of a fetus by analyzing a biological sample from a female subject carrying a fetus, the biological sample comprising cell-free DNA molecules from the female subject and the fetus, the method comprising: analyzing a plurality of cell-free DNA fragments from the biological sample to obtain sequence reads, the sequence reads including terminal sequences corresponding to terminals of the plurality of cell-free DNA fragments; determining, for each of the plurality of cell-free DNA fragments, a sequence motif for each of one or more terminal sequences of the cell-free DNA fragments; determining a relative frequency of a set of one or more sequence motifs that correspond to the end sequences of the plurality of cell-free DNA fragments, wherein the relative frequency of a sequence motif provides a proportion of the plurality of cell-free DNA fragments having end sequences that correspond to the sequence motif; determining a summary of the relative frequencies of the set of one or more sequence motifs; obtaining one or more calibration data points, each calibration data point designating a gestational age corresponding to an aggregate value, the one or more calibration data points having known gestational ages and determined from a plurality of calibration samples comprising cell-free DNA molecules; comparing the aggregated value to a calibration value for at least one calibration data point; and estimating a gestational age of the fetus based on the comparison.

22. 22. The method of claim 21, wherein the one or more calibration data points are a plurality of calibration data points that form a calibration function that approximates a measured aggregate value determined from the cell-free DNA molecules in the plurality of calibration samples of known gestational age.

23. 22. The method of claim 21, wherein the aggregate value is compared to a plurality of calibration values, each corresponding to one of the plurality of calibration samples.

24. 22. The method of claim 21 , wherein the calibration value of the at least one calibration data point corresponds to the aggregate value measured using the cell-free DNA molecules in at least one of the plurality of calibration samples.

25. 22. The method of claim 21, further comprising identifying the plurality of cell-free DNA fragments as derived from the fetus.

26. 26. The method of claim 25, wherein the plurality of cell-free DNA fragments is identified using fetal-specific alleles or fetal-specific epigenetic markers.

27. the plurality of cell-free DNA fragments For each of the sequence reads, determining a likelihood that the sequence read corresponds to the fetus based on an end sequence of the sequence read that includes a sequence motif of the set of one or more sequence motifs; comparing the likelihood to a threshold; and identifying the sequence read as originating from the fetus when the likelihood exceeds the threshold.

28. The method of any one of claims 1 to 27, wherein the set of one or more sequence motifs comprises N base positions, and the set of one or more sequence motifs comprises all combinations of N bases, where N is an integer greater than or equal to 3.

29. 28. The method of any one of claims 1 to 27, wherein the set of one or more sequence motifs are the top M sequence motifs having the largest difference between two types of DNA determined in one or more reference samples, where M is an integer equal to or greater than 1.

30. 30. The method of claim 29, wherein said two types of DNA are said clinically relevant DNA and said other DNA.

31. 30. The method of claim 29, wherein the two types of DNA are derived from two reference samples having different classifications for the level of the pathology.

32. 28. The method of any one of claims 1 to 27, wherein the set of one or more sequence motifs are the top M most frequent sequence motifs occurring in one or more reference samples, where M is an integer equal to or greater than 1.

33. 33. The method of any one of claims 28 to 32, wherein the set of one or more sequence motifs comprises a plurality of sequence motifs, and the summary value comprises the sum of the relative frequencies of the set.

34. 34. The method of claim 33, wherein the sum is a weighted sum.

35. 35. The method of claim 34, wherein the aggregate value comprises an entropy term, the entropy term comprising a sum of terms comprising the weighted sum, each term comprising a relative frequency multiplied by the logarithm of the relative frequency.

36. A method according to any preceding claim, wherein the aggregate value corresponds to the variance of the relative frequencies.

37. The method of any one of claims 1 to 35, wherein the aggregated value comprises a final or intermediate output of a machine learning model.

38. 38. The method of claim 37, wherein the machine learning model uses clustering, a support vector machine, or logistic regression.

39. 1. A method of enriching a biological sample for clinically relevant DNA, said biological sample comprising said clinically relevant DNA and other DNA that is acellular, said method comprising: analyzing a plurality of cell-free DNA fragments from the biological sample to obtain sequence reads, the sequence reads including terminal sequences corresponding to terminals of the plurality of cell-free DNA fragments; determining, for each of the plurality of cell-free DNA fragments, a sequence motif for each of one or more terminal sequences of the cell-free DNA fragments; identifying a set of one or more sequence motifs that are present in said clinically relevant DNA at a higher relative frequency than in said other DNA; identifying a group of said sequence reads having said set of one or more sequence motifs in end sequences; For each sequence read in the group of sequence reads, determining a likelihood that the sequence read corresponds to the clinically relevant DNA based on an end sequence of the sequence read that contains a sequence motif of the set of one or more sequence motifs; comparing the likelihood to a threshold; storing the sequence read when the likelihood exceeds the threshold, thereby obtaining a stored sequence read; and analyzing the stored sequence reads to determine a characteristic of the clinically relevant DNA of the biological sample.

40. 40. The method of claim 39, wherein the characteristic of the clinically relevant DNA of the biological sample is (1) a fractional concentration of the clinically relevant DNA, (2) a level of pathology of the subject from which the biological sample was obtained, the level of pathology associated with the clinically relevant DNA, or (3) a gestational age of a fetus of a pregnant female from which the biological sample was obtained.

41. 40. The method of claim 39, further comprising measuring a size of the plurality of cell-free DNA fragments using the sequence reads, and wherein determining the likelihood that a particular sequence read corresponds to the clinically relevant DNA fragment is further based on a size of the cell-free DNA fragment that corresponds to the particular sequence read.

42. 40. The method of claim 39, further comprising measuring one or more methylation states at one or more sites in the cell-free DNA fragments corresponding to a particular sequence read, wherein determining the likelihood that the particular sequence read corresponds to the clinically relevant DNA is further based on the one or more methylation states.

43. Determining the sequence motif for each of one or more end sequences of the cell-free DNA fragments, aligning one or more sequence reads corresponding to the cell-free DNA fragments to a reference genome; identifying one or more bases in the reference genome adjacent to the end sequence; and determining the sequence motif using the terminal sequences and the one or more bases.

44. 1. A method of enriching a biological sample for clinically relevant DNA, said biological sample comprising said clinically relevant DNA and other DNA that is acellular, said method comprising: receiving a plurality of cell-free DNA fragments from the biological sample, wherein the clinically relevant DNA fragments have terminal sequences that include sequence motifs that occur at a higher relative frequency than the other DNA fragments; subjecting the plurality of cell-free DNA fragments to one or more probe molecules that detect the sequence motif in the terminal sequences of the plurality of cell-free DNA fragments, thereby obtaining detected DNA fragments; and using said detected DNA fragments to enrich said biological sample for said clinically relevant DNA fragments.

45. using the detected DNA fragments to enrich the biological sample for the clinically relevant DNA fragments; 45. The method of claim 44, comprising amplifying the detected DNA fragments.

46. 46. ​​The method of claim 45, wherein the one or more probe molecules comprise one or more enzymes that interrogate the plurality of cell-free DNA fragments and add new sequences that are used to amplify the detected DNA fragments.

47. using the detected DNA fragments to enrich the biological sample for the clinically relevant DNA fragments; capturing the detected DNA fragments; and discarding undetected DNA fragments.

48. 48. The method of claim 47, wherein one or more probe molecules are attached to a surface and detect the sequence motif in the terminal sequence by hybridization.

49. A computer product comprising a computer readable medium storing a number of instructions for controlling a computer system to implement the method of any one of claims 1 to 48.

50. 50. A computer product according to claim 49; and one or more processors for executing instructions stored on the computer-readable medium.

51. A system comprising means for carrying out the method according to any one of claims 1 to 48.

52. A system comprising one or more processors configured to carry out the method of any one of claims 1 to 48.

53. A system comprising modules for carrying out respectively the steps of the method according to any one of claims 1 to 48.

Images