Cell-free DNA (cfDNA) in human plasma is a class of biomarkers with many current and potential future diagnostic applications. Recent studies have shown that cfDNA molecules are not randomly fragmented and possess information related to their tissues of origin. Pathologies causing death of cells from particular tissues result in perturbations in the relative distribution of DNA from the affected tissues. Such tissue-of-origin analysis is particularly useful in the development of liquid biopsies for cancer. It is therefore of value to accurately determine the relative contributions of the tissues to the plasma DNA pool in a simultaneous manner. In this work, we report that in open chromatin regions, cfDNA molecules show characteristic fragmentation patterns reflected by sequencing coverage imbalance and differentially phased fragment end signals. The latter refers to differences in the read densities of sequences corresponding to the orientation of the upstream and downstream ends of cfDNA molecules in relation to the reference genome. Such cfDNA fragmentation patterns preferentially occur in tissue-specific open chromatin regions where the corresponding tissues contributed DNA into the plasma. Quantitative analyses of such signals allow measurement of the relative contributions of various tissues toward the plasma DNA pool. These findings were validated by plasma DNA sequencing data obtained from pregnant women, organ transplantation recipients, and cancer patients. Orientation-aware plasma DNA fragmentation analysis therefore has potential diagnostic applications in noninvasive prenatal testing, organ transplantation monitoring, and cancer liquid biopsy.

SignificanceCell-free DNA fragmentation is a nonrandom process. We showed that cell-free DNA fragments with ends at certain genomic coordinates had higher likelihoods of being derived from hepatocellular carcinoma. Other coordinates were associated with cell-free DNA molecules originating from the liver. Quantitative assessment of cell-free DNA molecules bearing these respective groups of end signatures correlated with the amounts of tumor-derived or liver-derived DNA in plasma. There were millions of tumor-associated plasma DNA end coordinates across the genome. Due to their high prevalence, they were more readily detectable than somatic mutations as a cancer signature in plasma. Hence, detection of tumor-associated plasma DNA ends may offer a cost-effective means of capturing evidence for the presence of cancer through liquid biopsy assessment.

Urinary cell-free (cf) DNA holds great potential as a completely noninvasive form of liquid biopsy. Knowledge of the composition of cfDNA by tissue of origin is useful for guiding its clinical uses. We conducted a global survey of urinary cfDNA composition using genomewide bisulfite sequencing. While previous studies focused on detecting cfDNA from a single source at a time, genomewide tissue specific methylation signatures allow us to simultaneously deduce the proportional contribution from each contributing tissue. The proportional contributions derived from methylation deconvolution are highly correlated with those calculated using allograft-derived donor-specific genetic markers in the urine of hematopoetic stem cell and renal transplant recipients. We found a large variation of proportional contributions from different tissues. We then assessed if cfDNA undergoes time-dependent fragmentation in urine by conducting in vitro incubation experiments. In vitro incubation at 37°C showed that urinary cfDNA concentration decreased under first order kinetics with a half-life of 2.6 to 5.1h. This is reflected in parallel by a decrease in the proportion of long fragments and increase in amplitude of 10bp periodicity seen in the cfDNA size profile. This global survey of urinary cfDNA has deepened our understanding of the composition, degradation and variation of cfDNA in the urinary tract and has laid a foundation for the use of genomewide urinary cfDNA sequencing as a molecular diagnostics tool.