Short single-stranded DNAs with putative non-canonical structures comprise a new class of plasma cell-free DNA

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

A review on the impact of single-stranded library preparation on plasma cell-free diversity for cancer detection

In clinical oncology, cell-free DNA (cfDNA) has shown immense potential in its ability to noninvasively detect cancer at various stages and monitor the progression of therapy. Despite the rapid improvements in cfDNA liquid biopsy approaches, achieving the required sensitivity to detect rare tumor-derived cfDNA still remains a challenge. For next-generation sequencing, the perceived presentation of cfDNA is strongly linked to the extraction and library preparation protocols. Conventional double-stranded DNA library preparation (dsDNA-LP) focuses on assessing ~167bp double-stranded mononucleosomal (mncfDNA) and its other oligonucleosomal cell-free DNA counterparts in plasma. However, dsDNA-LP methods fail to include short, single-stranded, or nicked DNA in the final library preparation, biasing the representation of the actual cfDNA populations in plasma. The emergence of single-stranded library preparation (ssDNA-LP) strategies over the past decade has now allowed these other populations of cfDNA to be studied from plasma. With the use of ssDNA-LP, single-stranded, nicked, and ultrashort cfDNA can be comprehensively assessed for its molecular characteristics and clinical potential. In this review, we overview the current literature on applications of ssDNA-LP on plasma cfDNA from a potential cancer liquid biopsy perspective. To this end, we discuss the molecular principles of single-stranded DNA adapter ligation, how library preparation contributes to the understanding of native cfDNA characteristics, and the potential for ssDNA-LP to improve the sensitivity of circulating tumor DNA detection. Additionally, we review the current literature on the newly reported species of plasma ultrashort single-stranded cell-free DNA plasma, which appear biologically distinct from mncfDNA. We conclude with a discussion of future perspectives of ssDNA-LP for liquid biopsy endeavors.

What fraction of cellular DNA turnover becomes cfDNA?

Cell-free DNA (cfDNA) tests use small amounts of DNA in the bloodstream as biomarkers. While it is thought that cfDNA is largely released by dying cells, the proportion of dying cells' DNA that reaches the bloodstream is unknown. Here, we integrate estimates of cellular turnover rates to calculate the expected amount of cfDNA. By comparing this to the actual amount of cell type-specific cfDNA, we estimate the proportion of DNA reaching plasma as cfDNA. We demonstrate that <10% of the DNA from dying cells is detectable in plasma, and the ratios of measured to expected cfDNA levels vary a thousand-fold among cell types, often reaching well below 0.1%. The analysis suggests that local clearance, presumably via phagocytosis, takes up most of the dying cells' DNA. Insights into the underlying mechanism may help to understand the physiological significance of cfDNA and improve the sensitivity of liquid biopsies.

Ultra-short cell-free DNA fragments enhance cancer early detection in a multi-analyte blood test combining mutation, protein and fragmentomics

OBJECTIVES

Cancer morbidity and mortality can be reduced if the cancer is detected early. Cell-free DNA (cfDNA) fragmentomics emerged as a novel epigenetic biomarker for early cancer detection, however, it is still at its infancy and requires technical improvement. We sought to apply a single-strand DNA sequencing technology, for measuring genetic and fragmentomic features of cfDNA and evaluate the performance in detecting multiple cancers.

METHODS

Blood samples of 364 patients from six cancer types (colorectal, esophageal, gastric, liver, lung, and ovarian cancers) and 675 healthy individuals were included in this study. Circulating tumor DNA mutations, cfDNA fragmentomic features and a set of protein biomarkers were assayed. Sensitivity and specificity were reported by cancer types and stages.

RESULTS

Circular Ligation Amplification and sequencing (CLAmp-seq), a single-strand DNA sequencing technology, yielded a population of ultra-short fragments (<100 bp) than double-strand DNA preparation protocols and reveals a more significant size difference between cancer and healthy cfDNA fragments (25.84 bp vs. 16.05 bp). Analysis of the subnucleosomal peaks in ultra-short cfDNA fragments indicates that these peaks are regulatory element "footprints" and correlates with gene expression and cancer stages. At 98 % specificity, a prediction model using ctDNA mutations alone showed an overall sensitivity of 46 %; sensitivity reaches 60 % when protein is added, sensitivity further increases to 66 % when fragmentomics is also integrated. More improvements observed for samples representing earlier cancer stages than later ones.

CONCLUSIONS

These results suggest synergistic properties of protein, genetic and fragmentomics features in the identification of early-stage cancers.

Distinct Features of Plasma Ultrashort Single-Stranded Cell-Free DNA as Biomarkers for Lung Cancer Detection

BACKGROUND

Using broad range cell-free DNA sequencing (BRcfDNA-Seq), a nontargeted next-generation sequencing (NGS) methodology, we previously identified a novel class of approximately 50 nt ultrashort single-stranded cell-free DNA (uscfDNA) in plasma that is distinctly different from 167 bp mononucleosomal cell-free DNA (mncfDNA). We hypothesize that uscfDNA possesses characteristics that are useful for disease detection.

METHODS

Using BRcfDNA-Seq, we examined both cfDNA populations in the plasma of 18 noncancer controls and 14 patients with late-stage nonsmall cell lung carcinoma (NSCLC). In comparison to mncfDNA, we assessed whether functional element (FE) peaks, fragmentomics, end-motifs, and G-Quadruplex (G-Quad) signatures could be useful features of uscfDNA for NSCLC determination.

RESULTS

In noncancer participants, compared to mncfDNA, uscfDNA fragments showed a 45.2-fold increased tendency to form FE peaks (enriched in promoter, intronic, and exonic regions), demonstrated a distinct end-motif-frequency profile, and presented with a 4.9-fold increase in G-Quad signatures. Within NSCLC participants, only the uscfDNA population had discoverable FE peak candidates. Additionally, uscfDNA showcased different end-motif-frequency candidates distinct from mncfDNA. Although both cfDNA populations showed increased fragmentation in NSCLC, the G-Quad signatures were more discriminatory in uscfDNA. Compilation of cfDNA features using principal component analysis revealed that the first 5 principal components of both cfDNA subtypes had a cumulative explained variance of >80%.

CONCLUSIONS

These observations indicate that the distinct biological processes of uscfDNA and that FE peaks, fragmentomics, end-motifs, and G-Quad signatures are uscfDNA features with promising biomarker potential. These findings further justify its exploration as a distinct class of biomarker to augment pre-existing liquid biopsy approaches.

Multidimensional fragmentomic profiling of cell-free DNA released from patient-derived organoids

BACKGROUND

Fragmentomics, the investigation of fragmentation patterns of cell-free DNA (cfDNA), has emerged as a promising strategy for the early detection of multiple cancers in the field of liquid biopsy. However, the clinical application of this approach has been hindered by a limited understanding of cfDNA biology. Furthermore, the prevalence of hematopoietic cell-derived cfDNA in plasma complicates the in vivo investigation of tissue-specific cfDNA other than that of hematopoietic origin. While conventional two-dimensional cell lines have contributed to research on cfDNA biology, their limited representation of in vivo tissue contexts underscores the need for more robust models. In this study, we propose three-dimensional organoids as a novel in vitro model for studying cfDNA biology, focusing on multifaceted fragmentomic analyses.

RESULTS

We established nine patient-derived organoid lines from normal lung airway, normal gastric, and gastric cancer tissues. We then extracted cfDNA from the culture medium of these organoids in both proliferative and apoptotic states. Using whole-genome sequencing data from cfDNA, we analyzed various fragmentomic features, including fragment size, footprints, end motifs, and repeat types at the end. The distribution of cfDNA fragment sizes in organoids, especially in apoptosis samples, was similar to that found in plasma, implying occupancy by mononucleosomes. The footprints determined by sequencing depth exhibited distinct patterns depending on fragment sizes, reflecting occupancy by a variety of DNA-binding proteins. Notably, we discovered that short fragments (< 118 bp) were exclusively enriched in the proliferative state and exhibited distinct fragmentomic profiles, characterized by 3 bp palindromic end motifs and specific repeats.

CONCLUSIONS

In conclusion, our results highlight the utility of in vitro organoid models as a valuable tool for studying cfDNA biology and its associated fragmentation patterns. This, in turn, will pave the way for further enhancements in noninvasive cancer detection methodologies based on fragmentomics.

Multi-faceted attributes of salivary cell-free DNA as liquid biopsy biomarkers for gastric cancer detection

BACKGROUND

Recent advances in circulating cell-free DNA (cfDNA) analysis from biofluids have opened new avenues for liquid biopsy (LB). However, current cfDNA LB assays are limited by the availability of existing information on established genotypes associated with tumor tissues. Certain cancers present with a limited list of established mutated cfDNA biomarkers, and thus, nonmutated cfDNA characteristics along with alternative biofluids are needed to broaden the available cfDNA targets for cancer detection. Saliva is an intriguing and accessible biofluid that has yet to be fully explored for its clinical utility for cancer detection.

METHODS

In this report, we employed a low-coverage single stranded (ss) library NGS pipeline "Broad-Range cell-free DNA-Seq" (BRcfDNA-Seq) using saliva to comprehensively investigate the characteristics of salivary cfDNA (ScfDNA). The identification of cfDNA features has been made possible by applying novel cfDNA processing techniques that permit the incorporation of ultrashort, ss, and jagged DNA fragments. As a proof of concept using 10 gastric cancer (GC) and 10 noncancer samples, we examined whether ScfDNA characteristics, including fragmentomics, end motif profiles, microbial contribution, and human chromosomal mapping, could differentiate between these two groups.

RESULTS

Individual and integrative analysis of these ScfDNA features demonstrated significant differences between the two cohorts, suggesting that disease state may affect the ScfDNA population by altering nuclear cleavage or the profile of contributory organism cfDNA to total ScfDNA. We report that principal component analysis integration of several aspects of salivary cell-free DNA fragmentomic profiles, genomic element profiles, end-motif sequence patterns, and distinct oral microbiome populations can differentiate the two populations with a p value of < 0.0001 (PC1).

CONCLUSION

These novel features of ScfDNA characteristics could be clinically useful for improving saliva-based LB detection and the eventual monitoring of local or systemic diseases.

A highly efficient scheme for library preparation from single-stranded DNA

Although methods for sequencing library preparation from double-stranded DNA are well established, those from single-stranded DNA (ssDNA) have not been well studied. Further, the existing methods have limitations in efficiency and yield. Therefore, we developed a highly efficient procedure for sequencing library preparation from ssDNA. In this method, the first adaptor tagging of ssDNA is performed using terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation, which we reported recently. After complementary strand synthesis using the adaptor-tagged ssDNA, second adaptor tagging via Vaccinia virus topoisomerase I (VTopoI or TOPO)-based adaptor ligation is performed. With additional steps for degradation, repression, and removal of the adaptor dimer, the proposed TACS-TOPO scheme realizes adaptor dimer-free sequencing library preparation from ssDNA samples of 24 pg. The TACS-TOPO scheme was successfully applied to cell-free DNA analysis with amplification-free library preparation from 50 µL of human serum. A modified TACS-TOPO scheme was also applied to DNA extracted from ancient human bones, bringing two to eight times more library yields than those using a conventional library preparation protocol. The procedures for preparing VTopoI and its complex with a double-stranded oligonucleotide adaptor are also described. Overall, the proposed TACS-TOPO scheme can facilitate practical and sensitive sequencing analysis of ssDNA.

Cross-platform comparisons for targeted bisulfite sequencing of MGISEQ-2000 and NovaSeq6000

BACKGROUND

An accurate and reproducible next-generation sequencing platform is essential to identify malignancy-related abnormal DNA methylation changes and translate them into clinical applications including cancer detection, prognosis, and surveillance. However, high-quality DNA methylation sequencing has been challenging because poor sequence diversity of the bisulfite-converted libraries severely impairs sequencing quality and yield. In this study, we tested MGISEQ-2000 Sequencer's capability of DNA methylation sequencing with a published non-invasive pancreatic cancer detection assay, using NovaSeq6000 as the benchmark.

RESULTS

We sequenced a series of synthetic cell-free DNA (cfDNA) samples with different tumor fractions and found MGISEQ-2000 yielded data with similar quality as NovaSeq6000. The methylation levels measured by MGISEQ-2000 demonstrated high consistency with NovaSeq6000. Moreover, MGISEQ-2000 showed a comparable analytic sensitivity with NovaSeq6000, suggesting its potential for clinical detection. As to evaluate the clinical performance of MGISEQ-2000, we sequenced 24 clinical samples and predicted the pathology of the samples with a clinical diagnosis model, PDACatch classifier. The clinical model performance of MGISEQ-2000's data was highly consistent with that of NovaSeq6000's data, with the area under the curve of 1. We also tested the model's robustness with MGISEQ-2000's data when reducing the sequencing depth. The results showed that MGISEQ-2000's data showed matching robustness of the PDACatch classifier with NovaSeq6000's data.

CONCLUSIONS

Taken together, MGISEQ-2000 demonstrated similar data quality, consistency of the methylation levels, comparable analytic sensitivity, and matching clinical performance, supporting its application in future non-invasive early cancer detection investigations by detecting distinct methylation patterns of cfDNAs.

Saliva diagnostics: Salivaomics, saliva exosomics, and saliva liquid biopsy

BACKGROUND

Each day, humans produce approximately 0.5 through 1.5 liters of saliva, a biofluid that is rich in biological omic constituents. Our lack of understanding how omic biomarkers migrate from diseased tissue to the saliva has impeded the clinical translation of saliva testing. Although such biomarkers must be conveyed via the vascular and lymphatic systems to the salivary glands, the molecular mechanisms that underlie this transport remain unclear. Although COVID-19 highlighted the need for rapid and reliable testing for infectious diseases, it represents only one of the many health conditions that potentially can be diagnosed using a saliva sample.

TYPES OF STUDIES REVIEWED

The authors discuss salivaomics, saliva exosomics, and the mechanisms on which saliva diagnostics are based and introduce a novel electrochemical sensing technology that may be exploited for saliva liquid biopsy.

RESULTS

The utility of saliva for screening for lung cancer is under investigation. Saliva testing may be used to stratify patients, monitor treatment response, and detect disease recurrence. The authors also highlight the landscapes of saliva-based SARS-CoV-2 testing and ultrashort cell-free DNA and outline how these fields are likely to evolve in the near future.

PRACTICAL IMPLICATIONS

Breakthroughs in the study of saliva research, therefore, will facilitate clinical deployment of saliva-based testing.

Genome-Wide Profiling of Endogenous Single-Stranded DNA Using the SSiNGLe-P1 Method

Endogenous single-stranded DNA (essDNA) can form in a mammalian genome as the result of a variety of molecular processes and can both play important roles inside the cell as well as have detrimental consequences to genome integrity, much of which remains to be fully understood. Here, we established the SSiNGLe-P1 approach based on limited digestion by P1 endonuclease for high-throughput genome-wide identification of essDNA regions. We applied this method to profile essDNA in both human mitochondrial and nuclear genomes. In the mitochondrial genome, the profiles of essDNA provide new evidence to support the strand-displacement model of mitochondrial DNA replication. In the nuclear genome, essDNA regions were found to be enriched in certain types of functional genomic elements, particularly, the origins of DNA replication, R-loops, and to a lesser degree, in promoters. Furthermore, interestingly, many of the essDNA regions identified by SSiNGLe-P1 have not been annotated and thus could represent yet unknown functional elements.

Multi-Faceted Attributes of Salivary Cell-free DNA as Liquid Biopsy Biomarkers for Gastric Cancer Detection

Background

Recent advances in circulating cell-free DNA (cfDNA) analysis from biofluids have opened new avenues for liquid biopsy (LB). However, current cfDNA LB assays are limited by the availability of existing information on established genotypes associated with tumor tissues. Certain cancers present with a limited list of established mutated cfDNA biomarkers, and thus, nonmutated cfDNA characteristics along with alternative biofluids are needed to broaden the available cfDNA targets for cancer detection. Saliva is an intriguing and accessible biofluid that has yet to be fully explored for its clinical utility for cancer detection.

Methods

In this report, we employed a low-coverage single stranded (ss) library NGS pipeline "Broad-Range cell-free DNA-Seq" (BRcfDNA-Seq) using saliva to comprehensively investigate the characteristics of salivary cfDNA (ScfDNA). The identification of cfDNA features has been made possible by applying novel cfDNA processing techniques that permit the incorporation of ultrashort, ss, and jagged DNA fragments. As a proof of concept using 10 gastric cancer (GC) and 10 noncancer samples, we examined whether ScfDNA characteristics, including fragmentomics, end motif profiles, microbial contribution, and human chromosomal mapping, could differentiate between these two groups.

Results

Individual and integrative analysis of these ScfDNA features demonstrated significant differences between the two cohorts, suggesting that disease state may affect the ScfDNA population by altering nuclear cleavage or the profile of contributory organism cfDNA to total ScfDNA. We report that principal component analysis integration of several aspects of salivary cell-free DNA fragmentomic profiles, genomic element profiles, end-motif sequence patterns, and distinct oral microbiome populations can differentiate the two populations with a p value of < 0.0001 (PC1).

Conclusion

These novel features of ScfDNA characteristics could be clinically useful for improving saliva-based LB detection and the eventual monitoring of local or systemic diseases.

Tracing the Origin of Cell-Free DNA Molecules through Tissue-Specific Epigenetic Signatures

All cell and tissue types constantly release DNA fragments into human body fluids by various mechanisms including programmed cell death, accidental cell degradation and active extrusion. Particularly, cell-free DNA (cfDNA) in plasma or serum has been utilized for minimally invasive molecular diagnostics. Disease onset or pathological conditions that lead to increased cell death alter the contribution of different tissues to the total pool of cfDNA. Because cfDNA molecules retain cell-type specific epigenetic features, it is possible to infer tissue-of-origin from epigenetic characteristics. Recent research efforts demonstrated that analysis of, e.g., methylation patterns, nucleosome occupancy, and fragmentomics determined the cell- or tissue-of-origin of individual cfDNA molecules. This novel tissue-of origin-analysis enables to estimate the contributions of different tissues to the total cfDNA pool in body fluids and find tissues with increased cell death (pathologic condition), expanding the portfolio of liquid biopsies towards a wide range of pathologies and early diagnosis. In this review, we summarize the currently available tissue-of-origin approaches and point out the next steps towards clinical implementation.

Plasma contains ultrashort single-stranded DNA in addition to nucleosomal cell-free DNA

Plasma cell-free DNA is being widely explored as a biomarker for clinical screening. Currently, methods are optimized for the extraction and detection of double-stranded mononucleosomal cell-free DNA of ∼160bp in length. We introduce uscfDNA-seq, a single-stranded cell-free DNA next-generation sequencing pipeline, which bypasses previous limitations to reveal a population of ultrashort single-stranded cell-free DNA in human plasma. This species has a modal size of 50nt and is distinctly separate from mononucleosomal cell-free DNA. Treatment with single-stranded and double-stranded specific nucleases suggests that ultrashort cell-free DNA is primarily single-stranded. It is distributed evenly across chromosomes and has a similar distribution profile over functional elements as the genome, albeit with an enrichment over promoters, exons, and introns, which may be suggestive of a terminal state of genome degradation. The examination of this cfDNA species could reveal new features of cell death pathways or it can be used for cell-free DNA biomarker discovery.

Characteristics, origin, and potential for cancer diagnostics of ultrashort plasma cell-free DNA

Current evidence suggests that plasma cell-free DNA (cfDNA) is fragmented around a mode of 166 bp. Data supporting this view has been mainly acquired through the analysis of double-stranded cfDNA. The characteristics and diagnostic potential of single-stranded and damaged double-stranded cfDNA in healthy individuals and cancer patients remain unclear. Here, through a combination of high-affinity magnetic bead-based DNA extraction and single-stranded DNA sequencing library preparation (MB-ssDNA), we report the discovery of a large proportion of cfDNA fragments centered at ∼50 bp. We show that these "ultrashort" cfDNA fragments have a greater relative abundance in plasma of healthy individuals (median = 19.1% of all sequenced cfDNA fragments, n = 28) than in plasma of patients with cancer (median = 14.2%, n = 21, P < 0.0001). The ultrashort cfDNA fragments map to accessible chromatin regions of blood cells, particularly in promoter regions with the potential to adopt G-quadruplex (G4) DNA secondary structures. G4-positive promoter chromatin accessibility is significantly enriched in ultrashort plasma cfDNA fragments from healthy individuals relative to patients with cancers (P < 0.0001), in whom G4-cfDNA enrichment is inversely associated with copy number aberration-inferred tumor fractions. Our findings redraw the landscape of cfDNA fragmentation by identifying and characterizing a novel population of ultrashort plasma cfDNA fragments. Sequencing of MB-ssDNA libraries could facilitate the characterization of gene regulatory regions and DNA secondary structures via liquid biopsy. Our data underline the diagnostic potential of ultrashort cfDNA through classification for cancer patients.