Unsupervised detection of fragment length signatures of circulating tumor DNA using non-negative matrix factorization
Abstract
Sequencing of cell-free DNA (cfDNA) is currently being used to detect cancer by searching both for mutational and non-mutational alterations. Recent work has shown that the length distribution of cfDNA fragments from a cancer patient can inform tumor load and type. Here, we propose non-negative matrix factorization (NMF) of fragment length distributions as a novel and completely unsupervised method for studying fragment length patterns in cfDNA. Using shallow whole-genome sequencing (sWGS) of cfDNA from a cohort of patients with metastatic castration-resistant prostate cancer (mCRPC), we demonstrate how NMF accurately infers the true tumor fragment length distribution as an NMF component - and that the sample weights of this component correlate with ctDNA levels (r=0.75). We further demonstrate how using several NMF components enables accurate cancer detection on data from various early stage cancers (AUC=0.96). Finally, we show that NMF, when applied across genomic regions, can be used to discover fragment length signatures associated with open chromatin.
Data availability
Danish law requires ethical approval of any specific research aim and imposes restrictions on sharing of personal data. This means that the prostate cancer data used in this article cannot be uploaded to international databases. External researchers (academic or commercial) interested in analysing the prostate dataset (including any derivatives of it) will need to contact the Data Access Committee via email to kdso@clin.au.dk. The Data Access Committee is formed of co-authors Karina Dalsgaard Sørensen and Michael Borre, and Ole Halfdan Larsen (Department Head Consultant, Department of Clinical Medicine, Aarhus University). Due to Danish Law, for the authors to be allowed to share the data (pseudonymised) it will require prior approval from The Danish National Committee on Health Research Ethics (or similar) for the specific new research goal. The author (based in Denmark) has to submit the application for ethical approval, with the external researcher(s) as named collaborator(s)). In addition to ethical approval, a Collaboration Agreement and a Data Processing Agreement is required, both of which must be approved by the legal office of the institution of the author (data owner) and the legal office of the institution of the external researcher (data processor). Raw fragment length distributions along with ctDNA% estimates are available in Supplementary File 1.
Article and author information
Author details
Funding
The Independent Research Fund Denmark (Sapere Aude Research Leader)
- Søren Besenbacher
The Danish Cancer Society
- Karina Dalsgaard Sørensen
The Central Denmark Region Health Fund
- Karina Dalsgaard Sørensen
Aarhus Universitet (Graduate School of Health)
- Maibritt Nørgaard
Direktør Emil C. Hertz og Hustru Inger Hertz Fond
- Karina Dalsgaard Sørensen
KV Fonden
- Karina Dalsgaard Sørensen
Raimond og Dagmar Ringgård-Bohns Fond
- Karina Dalsgaard Sørensen
Beckett Fonden
- Karina Dalsgaard Sørensen
Snedkermester Sophus Jacobsen og Hustru Astrid Jacobsens Fond
- Karina Dalsgaard Sørensen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Human subjects: The prostate study was approved by The National Committee on Health Research Ethics (#1901101) and notified to The Danish Data Protection Agency (#1-16-02-366-15). All patients provided written informed consent.
Copyright
© 2022, Renaud et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,544
- views
-
- 497
- downloads
-
- 13
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
Expression of NPRL2/TUSC4, a tumor-suppressor gene, is reduced in many cancers including NSCLC. Restoration of NPRL2 induces DNA damage, apoptosis, and cell-cycle arrest. We investigated NPRL2 antitumor immune responses in aPD1R/KRAS/STK11mt NSCLC in humanized-mice. Humanized-mice were generated by transplanting fresh human cord blood-derived CD34 stem cells into sub-lethally irradiated NSG mice. Lung-metastases were developed from KRAS/STK11mt/aPD1R A549 cells and treated with NPRL2 w/wo pembrolizumab. NPRL2-treatment reduced lung metastases significantly, whereas pembrolizumab was ineffective. Antitumor effect was greater in humanized than non-humanized-mice. NPRL2 + pembrolizumab was not synergistic in KRAS/STK11mt/aPD1R tumors but was synergistic in KRASwt/aPD1S H1299. NPRL2 also showed a significant antitumor effect on KRASmt/aPD1R LLC2 syngeneic-tumors. The antitumor effect was correlated with increased infiltration of human cytotoxic-T, HLA-DR+DC, CD11c+DC, and downregulation of myeloid and regulatory-T cells in TME. Antitumor effect was abolished upon in-vivo depletion of CD8-T, macrophages, and CD4-T cells whereas remained unaffected upon NK-cell depletion. A distinctive protein-expression profile was found after NPRL2 treatment. IFNγ, CD8b, and TBX21 associated with T-cell functions were significantly increased, whereas FOXP3, TGFB1/B2, and IL-10RA were strongly inhibited by NPRL2. A list of T-cell co-inhibitory molecules was also downregulated. Restoration of NPRL2 exhibited significantly slower tumor growth in humanized-mice, which was associated with increased presence of human cytotoxic-T, and DC and decreased percentage of Treg, MDSC, and TAM in TME. NPRL2-stable cells showed a substantial increase in colony-formation inhibition and heightened sensitivity to carboplatin. Stable-expression of NPRL2 resulted in the downregulation of MAPK and AKT-mTOR signaling. Taken-together, NPRL2 gene-therapy induces antitumor activity on KRAS/STK11mt/aPD1R tumors through DC-mediated antigen-presentation and cytotoxic immune-cell activation.
-
- Biochemistry and Chemical Biology
- Cancer Biology
Cancer cells display high levels of oncogene-induced replication stress (RS) and rely on DNA damage checkpoint for viability. This feature is exploited by cancer therapies to either increase RS to unbearable levels or inhibit checkpoint kinases involved in the DNA damage response. Thus far, treatments that combine these two strategies have shown promise but also have severe adverse effects. To identify novel, better-tolerated anticancer combinations, we screened a collection of plant extracts and found two natural compounds from the plant, Psoralea corylifolia, that synergistically inhibit cancer cell proliferation. Bakuchiol inhibited DNA replication and activated the checkpoint kinase CHK1 by targeting DNA polymerases. Isobavachalcone interfered with DNA double-strand break repair by inhibiting the checkpoint kinase CHK2 and DNA end resection. The combination of bakuchiol and isobavachalcone synergistically inhibited cancer cell proliferation in vitro. Importantly, it also prevented tumor development in xenografted NOD/SCID mice. The synergistic effect of inhibiting DNA replication and CHK2 signaling identifies a vulnerability of cancer cells that might be exploited by using clinically approved inhibitors in novel combination therapies.