Fragmentation landscape of cell-free DNA revealed by deconvolutional analysis of end motifs.

Cell-free DNA (cfDNA) fragmentation is nonrandom, at least partially mediated by various DNA nucleases, forming characteristic cfDNA end motifs. However, there is a paucity of tools for deciphering the relative contributions of cfDNA cleavage patterns related to underlying fragmentation factors. In this study, through non-negative matrix factorization algorithm, we used 256 5' 4-mer end motifs to identify distinct types of cfDNA cleavage patterns, referred to as "founder" end-motif profiles (F-profiles). F-profiles were associated with different DNA nucleases based on whether such patterns were disrupted in nuclease-knockout mouse models. Contributions of individual F-profiles in a cfDNA sample could be determined by deconvolutional analysis. We analyzed 93 murine cfDNA samples of different nuclease-deficient mice and identified six types of F-profiles. F-profiles I, II, and III were linked to deoxyribonuclease 1 like 3 (DNASE1L3), deoxyribonuclease 1 (DNASE1), and DNA fragmentation factor subunit beta (DFFB), respectively. We revealed that 42.9% of plasma cfDNA molecules were attributed to DNASE1L3-mediated fragmentation, whereas 43.4% of urinary cfDNA molecules involved DNASE1-mediated fragmentation. We further demonstrated that the relative contributions of F-profiles were useful to inform pathological states, such as autoimmune disorders and cancer. Among the six F-profiles, the use of F-profile I could inform the human patients with systemic lupus erythematosus. F-profile VI could be used to detect individuals with hepatocellular carcinoma, with an area under the receiver operating characteristic curve of 0.97. F-profile VI was more prominent in patients with nasopharyngeal carcinoma undergoing chemoradiotherapy. We proposed that this profile might be related to oxidative stress.

Epigenetic analysis of cell-free DNA by fragmentomic profiling.

Cell-free DNA (cfDNA) fragmentation patterns contain important molecular information linked to tissues of origin. We explored the possibility of using fragmentation patterns to predict cytosine-phosphate-guanine (CpG) methylation of cfDNA, obviating the use of bisulfite treatment and associated risks of DNA degradation. This study investigated the cfDNA cleavage profile surrounding a CpG (i.e., within an 11-nucleotide [nt] window) to analyze cfDNA methylation. The cfDNA cleavage proportion across positions within the window appeared nonrandom and exhibited correlation with methylation status. The mean cleavage proportion was ∼twofold higher at the cytosine of methylated CpGs than unmethylated ones in healthy controls. In contrast, the mean cleavage proportion rapidly decreased at the 1-nt position immediately preceding methylated CpGs. Such differential cleavages resulted in a characteristic change in relative presentations of CGN and NCG motifs at 5' ends, where N represented any nucleotide. CGN/NCG motif ratios were correlated with methylation levels at tissue-specific methylated CpGs (e.g., placenta or liver) (Pearson's absolute r > 0.86). cfDNA cleavage profiles were thus informative for cfDNA methylation and tissue-of-origin analyses. Using CG-containing end motifs, we achieved an area under a receiver operating characteristic curve (AUC) of 0.98 in differentiating patients with and without hepatocellular carcinoma and enhanced the positive predictive value of nasopharyngeal carcinoma screening (from 19.6 to 26.8%). Furthermore, we elucidated the feasibility of using cfDNA cleavage patterns to deduce CpG methylation at single CpG resolution using a deep learning algorithm and achieved an AUC of 0.93. FRAGmentomics-based Methylation Analysis (FRAGMA) presents many possibilities for noninvasive prenatal, cancer, and organ transplantation assessment.

Universal Targeted Haplotyping by Droplet Digital PCR Sequencing and Its Applications in Noninvasive Prenatal Testing and Pharmacogenetics Analysis.

BACKGROUND: The analysis of haplotypes of variants is important for pharmacogenomics analysis and noninvasive prenatal testing for monogenic diseases. However, there is a lack of robust methods for targeted haplotyping. METHODS: We developed digital PCR haplotype sequencing (dHapSeq) for targeted haplotyping of variants, which is a method that compartmentalizes long DNA molecules into droplets. Within one droplet, 2 target regions are PCR amplified from one template molecule, and their amplicons are fused together. The fused products are then sequenced to determine the phase relationship of the single nucleotide polymorphism (SNP) alleles. The entire haplotype of 10s of SNPs can be deduced after the phase relationship of individual SNPs are determined in a pairwise manner. We applied dHapSeq to noninvasive prenatal testing in 4 families at risk for thalassemia and utilized it to detect NUDT15 diplotypes for predicting drug tolerance in pediatric acute lymphoblastic leukemia (72 cases and 506 controls). RESULTS: For SNPs within 40 kb, phase relation can be determined with 100% accuracy. In 7 trio families, the haplotyping results for 97 SNPs spanning 185 kb determined by dHapSeq were concordant with the results deduced from the genotypes of both parents and the fetus. In 4 thalassemia families, a 19.3-kb Southeast Asian deletion was successfully phased with 97 downstream SNPs, enabling noninvasive determination of fetal inheritance using relative haplotype dosage analysis. In the NUDT15 analysis, the variant status and phase of the variants were successfully determined in all cases and controls. CONCLUSIONS: The dHapSeq represents a robust and scalable haplotyping approach with numerous clinical and research applications.

Droplet digital PCR is a cost-effective method for analyzing long cell-free DNA in maternal plasma: Application in preeclampsia.

OBJECTIVE: Long cell-free DNA (cfDNA) can be found in the plasma of pregnant women and cancer patients. We investigated if droplet digital PCR (ddPCR) can analyze such molecules for diagnostic purposes using preeclampsia as a model. METHOD: Plasma samples from ten preeclamptic and sixteen normal pregnancies were analyzed. Two ddPCR assays targeting a single-copy gene, VCP, and one ddPCR assay targeting LINE-1 repetitive regions were used to measure the percentages of long cfDNA >533, 1001, and 170 bp, respectively. The LINE-1 assay was developed as guided by in silico PCR analyses to better differentiate preeclamptic and normal pregnancies. RESULTS: Preeclamptic patients had a significantly lower median percentage of long cfDNA than healthy pregnant controls, as determined by the LINE-1 170 bp assay (28.9% vs. 35.1%, p < 0.0001) and the VCP 533 bp assay (6.6% vs. 8.7%, p = 0.014). The LINE-1 assay provided a better differentiation than the VCP 533 bp assay (area under ROC curves, 0.94 vs. 0.79). CONCLUSION: ddPCR is a cost-effective approach for unlocking diagnostic information carried by long cfDNA in plasma and may have applications for the detection of preeclampsia. Further longitudinal studies with larger cohorts are required to assess the clinical utility of this test.

Sequencing-based counting and size profiling of plasma Epstein-Barr virus DNA enhance population screening of nasopharyngeal carcinoma.

Circulating tumor-derived DNA testing for cancer screening has recently been demonstrated in a prospective study on identification of nasopharyngeal carcinoma (NPC) among 20,174 asymptomatic individuals. Plasma EBV DNA, a marker for NPC, was detected using real-time PCR. While plasma EBV DNA was persistently detectable in 97.1% of the NPCs identified, ∼5% of the general population had transiently detectable plasma EBV DNA. We hypothesized that EBV DNA in plasma of subjects with or without NPC may have different molecular characteristics. We performed target-capture sequencing of plasma EBV DNA and identified differences in the abundance and size profiles of EBV DNA molecules within plasma of NPC and non-NPC subjects. NPC patients had significantly higher amounts of plasma EBV DNA, which showed longer fragment lengths. Cutoff values were established from an exploratory dataset and tested in a validation sample set. Adopting an algorithm that required a sample to concurrently pass cutoffs for EBV DNA counting and size measurements, NPCs were detected at a positive predictive value (PPV) of 19.6%. This represented superior performance compared with the PPV of 11.0% in the prospective screening study, which required participants with an initially detectable plasma EBV DNA result to be retested within 4 weeks. The observed differences in the molecular nature of EBV DNA molecules in plasma of subjects with or without NPC were successfully translated into a sequencing-based test that had a high PPV for NPC screening and achievable through single time-point testing.

Topologic Analysis of Plasma Mitochondrial DNA Reveals the Coexistence of Both Linear and Circular Molecules.

BACKGROUND: Cellular mitochondrial DNA (mtDNA) is organized as circular, covalently closed and double-stranded DNA. Studies have demonstrated the presence of short mtDNA fragments in plasma. It is not known whether circular mtDNA might concurrently exist with linear mtDNA in plasma. METHODS: We elucidated the topology of plasma mtDNA using restriction enzyme BfaI cleavage signatures on mtDNA fragment ends to differentiate linear and circular mtDNA. mtDNA fragments with both ends carrying BfaI cleavage signatures were defined as circular-derived mtDNA, whereas those with no cleavage signature or with 1 cleavage signature were defined as linear-derived mtDNA. An independent assay using exonuclease V to remove linear DNA followed by restriction enzyme MspI digestion was used for confirming the conclusions based on BfaI cleavage analysis. We analyzed the presence of BfaI cleavage signatures on plasma DNA ends in nonhematopoietically and hematopoietically derived DNA molecules by sequencing plasma DNA of patients with liver transplantation and bone marrow transplantation. RESULTS: Both linear and circular mtDNA coexisted in plasma. In patients with liver transplantation, donor-derived (i.e., liver) mtDNA molecules were mainly linear (median fraction, 91%; range, 75%-97%), whereas recipient-derived (i.e., hematopoietic) mtDNA molecules were mainly circular (median fraction, 88%; range, 77%-93%). The proportion of linear mtDNA was well correlated with liver DNA contribution in the plasma DNA pool (r = 0.83; P value = 0.0008). Consistent data were obtained from a bone marrow transplantation recipient in whom the donor-derived (i.e., hematopoietic) mtDNA molecules were predominantly circular. CONCLUSIONS: Linear and circular mtDNA molecules coexist in plasma and may have different tissue origins.

Liver- and Colon-Specific DNA Methylation Markers in Plasma for Investigation of Colorectal Cancers with or without Liver Metastases.

BACKGROUND: Measurement of DNA derived from different tissues in the circulating DNA pool can provide important information regarding the presence of many pathological conditions. However, existing methods involving genome-wide bisulfite sequencing are relatively expensive and may present challenges for large-scale analysis. METHODS: Through identifying differentially methylated regions in the liver and colon compared with other tissues, we identified 2 markers and developed corresponding droplet digital PCR assays. Plasma concentrations of liver-derived and colon-derived DNA were measured for 13 liver transplant recipients, 40 liver cancer patients, and 62 colorectal cancer (CRC) patients (27 with and 35 without liver metastases). RESULTS: In liver transplant recipients, the fractional concentration of liver-derived DNA measured using the liver-specific methylation marker and donor-specific alleles showed good correlation (Pearson R = 0.99). In liver cancer patients, the concentration of liver-derived DNA correlated positively with the maximal dimension of the tumor (Spearman R = 0.74). In CRC patients with and without liver metastasis, the plasma concentrations of colon-derived DNA (median, 138 copies/mL and 4 copies/mL, respectively) were increased compared with the 30 healthy controls (26 had undetectable concentrations). The absolute concentration of liver-derived DNA provided a better differentiation between CRC patients with and without liver metastasis compared with the fractional concentration (area under ROC curve, 0.85 vs 0.75). CONCLUSIONS: Quantitative analysis of plasma DNA with tissue-specific methylation patterns using droplet digital PCR is applicable for the investigation of cancers and assessing organ transplantation. This approach is useful for differentiating patients with and without metastases to other organs.

DNA of Erythroid Origin Is Present in Human Plasma and Informs the Types of Anemia.

BACKGROUND: There is much interest in the tissue of origin of circulating DNA in plasma. Data generated using DNA methylation markers have suggested that hematopoietic cells of white cell lineages are important contributors to the circulating DNA pool. However, it is not known whether cells of the erythroid lineage would also release DNA into the plasma. METHODS: Using high-resolution methylation profiles of erythroblasts and other tissue types, 3 genomic loci were found to be hypomethylated in erythroblasts but hypermethylated in other cell types. We developed digital PCR assays for measuring erythroid DNA using the differentially methylated region for each locus. RESULTS: Based on the methylation marker in the ferrochelatase gene, erythroid DNA represented a median of 30.1% of the plasma DNA of healthy subjects. In subjects with anemia of different etiologies, quantitative analysis of circulating erythroid DNA could reflect the erythropoietic activity in the bone marrow. For patients with reduced erythropoietic activity, as exemplified by aplastic anemia, the percentage of circulating erythroid DNA was decreased. For patients with increased but ineffective erythropoiesis, as exemplified by ß-thalassemia major, the percentage was increased. In addition, the plasma concentration of erythroid DNA was found to correlate with treatment response in aplastic anemia and iron deficiency anemia. Plasma DNA analysis using digital PCR assays targeting the other 2 differentially methylated regions showed similar findings. CONCLUSIONS: Erythroid DNA is a hitherto unrecognized major component of the circulating DNA pool and is a noninvasive biomarker for differential diagnosis and monitoring of anemia.

Testing for EGFR Variants in Pleural and Pericardial Effusion Cell-Free DNA in Patients With Non-Small Cell Lung Cancer.

Importance: Molecular testing in non-small cell lung cancer (NSCLC) is commonly limited by inadequate tumor sample. Plasma cell-free DNA (cfDNA) genotyping as a complementary test is specific but only moderately sensitive. Genotyping of cfDNA in pleural and pericardial effusion (PE-cfDNA) can further optimize molecular diagnostic yield and reduce the need for repeated biopsies. Objective: To prospectively validate droplet digital polymerase chain reaction (ddPCR) for detection of sensitizing EGFR variants and acquired Thr790Met variant (T790M) from PE-cfDNA in patients with NSCLC. Design, Setting, and Participants: This prospective diagnostic validation study was conducted between September 6, 2016, and January 21, 2021 at 2 major Hong Kong cancer centers. Patients with advanced NSCLC with both wild-type and variant EGFR status and exudative PE who underwent thoracocentesis or pericardiocentesis were randomly enrolled. Patients were either EGFR-tyrosine kinase inhibitor (TKI) naive (cohort 1) or EGFR-TKI treated but osimertinib naive (cohort 2). Enrolled patients underwent pleural- or pericardial-fluid and blood sampling for ddPCR EGFR testing. EGFR status results with ddPCR testing of PE-cfDNA and blood were compared with EGFR status in matched tumor biopsy or PE cell block samples. Main Outcomes and Measures: Specificity, sensitivity, and concordance of PE-cfDNA for detection of sensitizing EGFR variants and acquired T790M variation. Results: Among 171 patients (54% female) enrolled, there were 104 in cohort 1 and 67 in cohort 2. In cohort 1, 37% (38/102) were EGFR-variant positive; PE-cfDNA showed 97% sensitivity (95% CI, 92%-100%), 97% specificity (95% CI, 93%-100%), and 97% concordance (ĸ = 0.94, P < .001) for the detection of sensitizing EGFR variants. It was more sensitive than plasma in detecting sensitizing EGFR variants (97% vs 74%, P < .001). In cohort 2, 38% (15 of 40) were positive for the EGFR T790M variant; PE-cfDNA showed 87% sensitivity (95% CI, 69%-100%), 60% specificity (95% CI, 41%-79%), and 70% concordance (ĸ = 0.42, P = .004) for acquired T790M. The EGFR T790M variant was detected in 51% of PE-cfDNA vs 25% of PE cell block samples. Conclusions and Relevance: In this diagnostic study, EGFR variants could be accurately detected from PE-cfDNA in patients with NSCLC. More EGFR T790M was detected in PE-cfDNA than in guideline-recommended PE cell block preparations. These results suggest that PE-cfDNA can complement plasma and tumor genotyping for detecting EGFR variants in patients with advanced NSCLC.

Applications of genetic-epigenetic tissue mapping for plasma DNA in prenatal testing, transplantation and oncology.

We developed genetic-epigenetic tissue mapping (GETMap) to determine the tissue composition of plasma DNA carrying genetic variants not present in the constitutional genome through comparing their methylation profiles with relevant tissues. We validated this approach by showing that, in pregnant women, circulating DNA carrying fetal-specific alleles was entirely placenta-derived. In lung transplant recipients, we showed that, at 72 hr after transplantation, the lung contributed only a median of 17% to the plasma DNA carrying donor-specific alleles, and hematopoietic cells contributed a median of 78%. In hepatocellular cancer patients, the liver was identified as the predominant source of plasma DNA carrying tumor-specific mutations. In a pregnant woman with lymphoma, plasma DNA molecules carrying cancer mutations and fetal-specific alleles were accurately shown to be derived from the lymphocytes and placenta, respectively. Analysis of tissue origin for plasma DNA carrying genetic variants is potentially useful for noninvasive prenatal testing, transplantation monitoring, and cancer screening.

Epigenetic Biomarkers in Cell-Free DNA and Applications in Liquid Biopsy.

Cell-free circulating DNA (cfDNA) in plasma has gained global interest as a diagnostic material for noninvasive prenatal testing and cancer diagnosis, or the so-called "liquid biopsy". Recent studies have discovered a great number of valuable genetic and epigenetic biomarkers for cfDNA-based liquid biopsy. Considering that the genetic biomarkers, e.g., somatic mutations, usually vary from case to case in most cancer patients, epigenetic biomarkers that are generalizable across various samples thus possess certain advantages. In this study, we reviewed the most recent studies and advances on utilizing epigenetic biomarkers for liquid biopsies. We first reviewed more traditional methods of using tissue/cancer-specific DNA methylation biomarkers and digital PCR or sequencing technologies for cancer diagnosis, as well as tumor origin determination. In the second part, we discussed the emerging novel approaches for exploring the biological basis and clinical applications of cfDNA fragmentation patterns. We further provided our comments and points of view on the future directions on epigenetic biomarker development for cfDNA-based liquid biopsies.

DCDB 2.0: a major update of the drug combination database.

Experience in clinical practice and research in systems pharmacology suggested the limitations of the current one-drug-one-target paradigm in new drug discovery. Single-target drugs may not always produce desired physiological effects on the entire biological system, even if they have successfully regulated the activities of their designated targets. On the other hand, multicomponent therapy, in which two or more agents simultaneously interact with multiple targets, has attracted growing attention. Many drug combinations consisting of multiple agents have already entered clinical practice, especially in treating complex and refractory diseases. Drug combination database (DCDB), launched in 2010, is the first available database that collects and organizes information on drug combinations, with an aim to facilitate systems-oriented new drug discovery. Here, we report the second major release of DCDB (Version 2.0), which includes 866 new drug combinations (1363 in total), consisting of 904 distinctive components. These drug combinations are curated from ∼140,000 clinical studies and the food and drug administration (FDA) electronic orange book. In this update, DCDB collects 237 unsuccessful drug combinations, which may provide a contrast for systematic discovery of the patterns in successful drug combinations. Database URL: http://www.cls.zju.edu.cn/dcdb/