GCparagon: evaluating and correcting GC biases in cell-free DNA at the fragment level.

Analyses of cell-free DNA (cfDNA) are increasingly being employed for various diagnostic and research applications. Many technologies aim to increase resolution, e.g. for detecting early-stage cancer or minimal residual disease. However, these efforts may be confounded by inherent base composition biases of cfDNA, specifically the over - and underrepresentation of guanine (G) and cytosine (C) sequences. Currently, there is no universally applicable tool to correct these effects on sequencing read-level data. Here, we present GCparagon, a two-stage algorithm for computing and correcting GC biases in cfDNA samples. In the initial step, length and GC base count parameters are determined. Here, our algorithm minimizes the inclusion of known problematic genomic regions, such as low-mappability regions, in its calculations. In the second step, GCparagon computes weights counterbalancing the distortion of cfDNA attributes (correction matrix). These fragment weights are added to a binary alignment map (BAM) file as alignment tags for individual reads. The GC correction matrix or the tagged BAM file can be used for downstream analyses. Parallel computing allows for a GC bias estimation below 1 min. We demonstrate that GCparagon vastly improves the analysis of regulatory regions, which frequently show specific GC composition patterns and will contribute to standardized cfDNA applications.

STAT3/LKB1 controls metastatic prostate cancer by regulating mTORC1/CREB pathway.

Prostate cancer (PCa) is a common and fatal type of cancer in men. Metastatic PCa (mPCa) is a major factor contributing to its lethality, although the mechanisms remain poorly understood. PTEN is one of the most frequently deleted genes in mPCa. Here we show a frequent genomic co-deletion of PTEN and STAT3 in liquid biopsies of patients with mPCa. Loss of Stat3 in a Pten-null mouse prostate model leads to a reduction of LKB1/pAMPK with simultaneous activation of mTOR/CREB, resulting in metastatic disease. However, constitutive activation of Stat3 led to high LKB1/pAMPK levels and suppressed mTORC1/CREB pathway, preventing mPCa development. Metformin, one of the most widely prescribed therapeutics against type 2 diabetes, inhibits mTORC1 in liver and requires LKB1 to mediate glucose homeostasis. We find that metformin treatment of STAT3/AR-expressing PCa xenografts resulted in significantly reduced tumor growth accompanied by diminished mTORC1/CREB, AR and PSA levels. PCa xenografts with deletion of STAT3/AR nearly completely abrogated mTORC1/CREB inhibition mediated by metformin. Moreover, metformin treatment of PCa patients with high Gleason grade and type 2 diabetes resulted in undetectable mTORC1 levels and upregulated STAT3 expression. Furthermore, PCa patients with high CREB expression have worse clinical outcomes and a significantly increased risk of PCa relapse and metastatic recurrence. In summary, we have shown that STAT3 controls mPCa via LKB1/pAMPK/mTORC1/CREB signaling, which we have identified as a promising novel downstream target for the treatment of lethal mPCa.

How to detect cancer early using cell-free DNA.

Cell-free circulating DNA (cfDNA) analyses may enable multi-cancer early detection (MCED), but the most suitable method is unknown. In this issue of Cancer Cell, Jamshidi et al. compare ten machine-learning classifiers using various cfDNA features and conclude that cfDNA methylation patterns are most promising for MCED testing.

A clinician's handbook for using ctDNA throughout the patient journey.

BACKGROUND: The promise of precision cancer medicine presently centers around the genomic sequence of a patient's tumor being translated into timely, actionable information to inform clinical care. The analysis of cell-free DNA from liquid biopsy, which contains circulating tumor DNA (ctDNA) in patients with cancer, has proven to be amenable to various settings in oncology. However, open questions surrounding the clinical validity and utility of plasma-based analyses have hindered widespread clinical adoption. MAIN BODY: Owing to the rapid evolution of the field, studies supporting the use of ctDNA as a biomarker throughout a patient's journey with cancer have accumulated in the last few years, warranting a review of the latest status for clinicians who may employ ctDNA in their precision oncology programs. In this work, we take a step back from the intricate coverage of detection approaches described extensively elsewhere and cover basic concepts around the practical implementation of next generation sequencing (NGS)-guided liquid biopsy. We compare relevant targeted and untargeted approaches to plasma DNA analysis, describe the latest evidence for clinical validity and utility, and highlight the value of genome-wide ctDNA analysis, particularly as it relates to early detection strategies and discovery applications harnessing the non-coding genome. CONCLUSIONS: The maturation of liquid biopsy for clinical application will require interdisciplinary efforts to address current challenges. However, patients and clinicians alike may greatly benefit in the future from its incorporation into routine oncology care.

Biallelic BUB1 mutations cause microcephaly, developmental delay, and variable effects on cohesion and chromosome segregation.

Budding uninhibited by benzimidazoles (BUB1) contributes to multiple mitotic processes. Here, we describe the first two patients with biallelic BUB1 germline mutations, who both display microcephaly, intellectual disability, and several patient-specific features. The identified mutations cause variable degrees of reduced total protein level and kinase activity, leading to distinct mitotic defects. Both patients' cells show prolonged mitosis duration, chromosome segregation errors, and an overall functional spindle assembly checkpoint. However, while BUB1 levels mostly affect BUBR1 kinetochore recruitment, impaired kinase activity prohibits centromeric recruitment of Aurora B, SGO1, and TOP2A, correlating with anaphase bridges, aneuploidy, and defective sister chromatid cohesion. We do not observe accelerated cohesion fatigue. We hypothesize that unresolved DNA catenanes increase cohesion strength, with concomitant increase in anaphase bridges. In conclusion, BUB1 mutations cause a neurodevelopmental disorder, with clinical and cellular phenotypes that partially resemble previously described syndromes, including autosomal recessive primary microcephaly, mosaic variegated aneuploidy, and cohesinopathies.

Biallelic truncating variants in ATP9A cause a novel neurodevelopmental disorder involving postnatal microcephaly and failure to thrive.

BACKGROUND: Genes implicated in the Golgi and endosomal trafficking machinery are crucial for brain development, and mutations in them are particularly associated with postnatal microcephaly (POM). METHODS: Exome sequencing was performed in three affected individuals from two unrelated consanguineous families presenting with delayed neurodevelopment, intellectual disability of variable degree, POM and failure to thrive. Patient-derived fibroblasts were tested for functional effects of the variants. RESULTS: We detected homozygous truncating variants in ATP9A. While the variant in family A is predicted to result in an early premature termination codon, the variant in family B affects a canonical splice site. Both variants lead to a substantial reduction of ATP9A mRNA expression. It has been shown previously that ATP9A localises to early and recycling endosomes, whereas its depletion leads to altered gene expression of components from this compartment. Consistent with previous findings, we also observed overexpression of ARPC3 and SNX3, genes strongly interacting with ATP9A. CONCLUSION: In aggregate, our findings show that pathogenic variants in ATP9A cause a novel autosomal recessive neurodevelopmental disorder with POM. While the physiological function of endogenous ATP9A is still largely elusive, our results underline a crucial role of this gene in endosomal transport in brain tissue.

Novel subtype of mucopolysaccharidosis caused by arylsulfatase K (ARSK) deficiency.

BACKGROUND: Mucopolysaccharidoses (MPS) are monogenic metabolic disorders that significantly affect the skeleton. Eleven enzyme defects in the lysosomal degradation of glycosaminoglycans (GAGs) have been assigned to the known MPS subtypes (I-IX). Arylsulfatase K (ARSK) is a recently characterised lysosomal hydrolase involved in GAG degradation that removes the 2-O-sulfate group from 2-sulfoglucuronate. Knockout of Arsk in mice was consistent with mild storage pathology, but no human phenotype has yet been described. METHODS: In this study, we report four affected individuals of two unrelated consanguineous families with homozygous variants c.250C>T, p.(Arg84Cys) and c.560T>A, p.(Leu187Ter) in ARSK, respectively. Functional consequences of the two ARSK variants were assessed by mutation-specific ARSK constructs derived by site-directed mutagenesis, which were ectopically expressed in HT1080 cells. Urinary GAG excretion was analysed by dimethylene blue and electrophoresis, as well as liquid chromatography/mass spectrometry (LC-MS)/MS analysis. RESULTS: The phenotypes of the affected individuals include MPS features, such as short stature, coarse facial features and dysostosis multiplex. Reverse phenotyping in two of the four individuals revealed additional cardiac and ophthalmological abnormalities. Mild elevation of dermatan sulfate was detected in the two subjects investigated by LC-MS/MS. Human HT1080 cells expressing the ARSK-Leu187Ter construct exhibited absent protein levels by western blot, and cells with the ARSK-Arg84Cys construct showed markedly reduced enzyme activity in an ARSK-specific enzymatic assay against 2-O-sulfoglucuronate-containing disaccharides as analysed by C18-reversed-phase chromatography followed by MS. CONCLUSION: Our work provides a detailed clinical and molecular characterisation of a novel subtype of mucopolysaccharidosis, which we suggest to designate subtype X.

Evolutionary conservation in noncoding genomic regions.

Humans may share more genomic commonalities with other species than previously thought. According to current estimates, ~5% of the human genome is functionally constrained, which is a much larger fraction than the ~1.5% occupied by annotated protein-coding genes. Hence, ~3.5% of the human genome comprises likely functional conserved noncoding elements (CNEs) preserved among organisms, whose common ancestors existed throughout hundreds of millions of years of evolution. As whole-genome sequencing emerges as a standard procedure in genetic analyses, interpretation of variations in CNEs, including the elucidation of mechanistic and functional roles, becomes a necessity. Here, we discuss the phenomenon of noncoding conservation via four dimensions (sequence, regulatory conservation, spatiotemporal expression, and structure) and the potential significance of CNEs in phenotype variation and disease.

Profiling of circulating tumor DNA and tumor tissue for treatment selection in patients with advanced and refractory carcinoma: a prospective, two-stage phase II Individualized Cancer Treatment trial.

BACKGROUND: Molecular profiling (MP) represents an opportunity to match patients to a targeted therapy and when tumor tissue is unavailable, circulating tumor deoxyribonucleic acid (ctDNA) can be harnessed as a non-invasive analyte for this purpose. We evaluated the success of a targeted therapy selected by profiling of ctDNA and tissue in patients with advanced and refractory carcinoma. PATIENTS AND METHODS: A blood draw as well as an optional tissue biopsy were obtained for MP. Whole-genome sequencing and a cancer hotspot panel were performed, and publicly available databases were used to match the molecular profile to targeted treatments. The primary endpoint was the progression-free survival (PFS) ratio (PFS on MP-guided therapy/PFS on the last evidence-based therapy), whereas the success of the targeted therapy was defined as a PFS ratio ⩾1.2. To test the impact of molecular profile-treatment matching strategies, we retrospectively analyzed selected cases via the CureMatch PreciGENE™ decision support algorithm. RESULTS: Interim analysis of 24 patients yielded informative results from 20 patients (83%). A potential tumor-specific drug could be matched in 11 patients (46%) and eight (33%) received a matched treatment. Median PFS in the matched treatment group was 61.5 days [interquartile range (IQR) 49.8-71.0] compared with 81.5 days (IQR 68.5-117.8) for the last evidence-based treatment, resulting in a median PFS ratio of 0.7 (IQR 0.6-0.9). Hence, as no patient experienced a PFS ratio ⩾1.2, the study was terminated. Except for one case, the CureMatch analysis identified either a two-drug or three-drug combination option. CONCLUSIONS: Our study employed a histotype-agnostic approach to harness molecular profiling data from both ctDNA and metastatic tumor tissue. The outcome results indicate that more innovative approaches to study design and matching algorithms are necessary to achieve improved patient outcomes.EU Clinical Trials Registry (https://www.clinicaltrialsregister.eu): EudraCT: 2014-005341-44.

Dynamic Changes of Circulating Tumor DNA Predict Clinical Outcome in Patients With Advanced Non-Small-Cell Lung Cancer Treated With Immune Checkpoint Inhibitors.

PURPOSE: Immune checkpoint inhibitors (ICIs) are increasingly being used in non-small-cell lung cancer (NSCLC), yet biomarkers predicting their benefit are lacking. We evaluated if on-treatment changes of circulating tumor DNA (ctDNA) from ICI start (t0) to after two cycles (t1) assessed with a commercial panel could identify patients with NSCLC who would benefit from ICI. PATIENTS AND METHODS: The molecular ctDNA response was evaluated as a predictor of radiographic tumor response and long-term survival benefit of ICI. To maximize the yield of ctDNA detection, de novo mutation calling was performed. Furthermore, the impact of clonal hematopoiesis (CH)-related variants as a source of biologic noise was investigated. RESULTS: After correction for CH-related variants, which were detected in 75 patients (44.9%), ctDNA was detected in 152 of 167 (91.0%) patients. We observed only a fair agreement of the molecular and radiographic response, which was even more impaired by the inclusion of CH-related variants. After exclusion of those, a ≥ 50% molecular response improved progression-free survival (10 v 2 months; hazard ratio [HR], 0.55; 95% CI, 0.39 to 0.77; P = .0011) and overall survival (18.4 v 5.9 months; HR, 0.44; 95% CI, 0.31 to 0.62; P < .0001) compared with patients not achieving this end point. After adjusting for clinical variables, ctDNA response and STK11/KEAP1 mutations (HR, 2.08; 95% CI, 1.4 to 3.0; P < .001) remained independent predictors for overall survival, irrespective of programmed death ligand-1 expression. A landmark survival analysis at 2 months (n = 129) provided similar results. CONCLUSION: On-treatment changes of ctDNA in plasma reveal predictive information for long-term clinical benefit in ICI-treated patients with NSCLC. A broader NSCLC patient coverage through de novo mutation calling and the use of a variant call set excluding CH-related variants improved the classification of molecular responders, but had no significant impact on survival.

On-treatment measurements of circulating tumor DNA during FOLFOX therapy in patients with colorectal cancer.

We addressed a significant unknown feature of circulating tumor DNA (ctDNA), i.e., how ctDNA levels change during chemotherapy, by serially monitoring ctDNA in patients with colorectal cancer during the 48-h application of FOLFOX. Surprisingly, we did not observe a spike in ctDNA as a sign of a responsive tumor, but instead ctDNA levels initially decreased and remained low in patients with stable disease or partial response. Our observations reveal further insights into cell destruction during chemotherapy with important implications for the management of patients.

Comparison of three commercial decision support platforms for matching of next-generation sequencing results with therapies in patients with cancer.

OBJECTIVE: Precision oncology depends on translating molecular data into therapy recommendations. However, with the growing complexity of next-generation sequencing-based tests, clinical interpretation of somatic genomic mutations has evolved into a formidable task. Here, we compared the performance of three commercial clinical decision support tools, that is, NAVIFY Mutation Profiler (NAVIFY; Roche), QIAGEN Clinical Insight (QCI) Interpret (QIAGEN) and CureMatch Bionov (CureMatch). METHODS: In order to obtain the current status of the respective tumour genome, we analysed cell-free DNA from patients with metastatic breast, colorectal or non-small cell lung cancer. We evaluated somatic copy number alterations and in parallel applied a 77-gene panel (AVENIO ctDNA Expanded Panel). We then assessed the concordance of tier classification approaches between NAVIFY and QCI and compared the strategies to determine actionability among all three platforms. Finally, we quantified the alignment of treatment suggestions across all decision tools. RESULTS: Each platform varied in its mode of variant classification and strategy for identifying druggable targets and clinical trials, which resulted in major discrepancies. Even the frequency of concordant actionable events for tier I-A or tier I-B classifications was only 4.3%, 9.5% and 28.4% when comparing NAVIFY with QCI, NAVIFY with CureMatch and CureMatch with QCI, respectively, and the obtained treatment recommendations differed drastically. CONCLUSIONS: Treatment decisions based on molecular markers appear at present to be arbitrary and dependent on the chosen strategy. As a consequence, tumours with identical molecular profiles would be differently treated, which challenges the promising concepts of genome-informed medicine.

Characterization of circulating breast cancer cells with tumorigenic and metastatic capacity.

Functional studies giving insight into the biology of circulating tumor cells (CTCs) remain scarce due to the low frequency of CTCs and lack of appropriate models. Here, we describe the characterization of a novel CTC-derived breast cancer cell line, designated CTC-ITB-01, established from a patient with metastatic estrogen receptor-positive (ER+ ) breast cancer, resistant to endocrine therapy. CTC-ITB-01 remained ER+ in culture, and copy number alteration (CNA) profiling showed high concordance between CTC-ITB-01 and CTCs originally present in the patient with cancer at the time point of blood draw. RNA-sequencing data indicate that CTC-ITB-01 has a predominantly epithelial expression signature. Primary tumor and metastasis formation in an intraductal PDX mouse model mirrored the clinical progression of ER+ breast cancer. Downstream ER signaling was constitutively active in CTC-ITB-01 independent of ligand availability, and the CDK4/6 inhibitor Palbociclib strongly inhibited CTC-ITB-01 growth. Thus, we established a functional model that opens a new avenue to study CTC biology.

Technical Evaluation of Commercial Mutation Analysis Platforms and Reference Materials for Liquid Biopsy Profiling.

Molecular profiling from liquid biopsy, in particular cell-free DNA (cfDNA), represents an attractive alternative to tissue biopsies for the detection of actionable targets and tumor monitoring. In addition to PCR-based assays, Next Generation Sequencing (NGS)-based cfDNA assays are now commercially available and are being increasingly adopted in clinical practice. However, the validity of these products as well as the clinical utility of cfDNA in the management of patients with cancer has yet to be proven. Within framework of the Innovative Medicines Initiative (IMI) program CANCER-ID we evaluated the use of commercially available reference materials designed for ctDNA testing and cfDNA derived from Diagnostic Leukaphereses (DLA) for inter- and intra-assay as well as intra- and inter-laboratory comparisons. In three experimental setups, a broad range of assays including ddPCR, MassARRAY and various NGS-based assays were tested. We demonstrate that both reference materials with predetermined VAFs and DLA samples are extremely useful for the performance assessment of mutation analysis platforms. Moreover, our data indicate a substantial variability of NGS assays with respect to sensitivity and specificity highlighting the importance of extensive validation of the test performance before offering these tests in clinical routine practice.

Shallow Whole-Genome Sequencing from Plasma Identifies FGFR1 Amplified Breast Cancers and Predicts Overall Survival.

Background: Focal amplification of fibroblast growth factor receptor 1 (FGFR1) defines a subgroup of breast cancers with poor prognosis and high risk of recurrence. We sought to demonstrate the potential of circulating cell-free DNA (cfDNA) analysis to evaluate FGFR1 copy numbers from a cohort of 100 metastatic breast cancer (mBC) patients. Methods: Formalin-fixed paraffin-embedded (FFPE) tissue samples were screened for FGFR1 amplification by FISH, and positive cases were confirmed with a microarray platform (OncoscanTM). Subsequently, cfDNA was evaluated by two approaches, i.e., mFAST-SeqS and shallow whole-genome sequencing (sWGS), to estimate the circulating tumor DNA (ctDNA) allele fraction (AF) and to evaluate the FGFR1 status. Results: Tissue-based analyses identified FGFR1 amplifications in 20/100 tumors. All cases with a ctDNA AF above 3% (n = 12) showed concordance for FGFR1 status between tissue and cfDNA. In one case, we were able to detect a high-level FGFR1 amplification, although the ctDNA AF was below 1%. Furthermore, high levels of ctDNA indicated an association with unfavorable prognosis based on overall survival. Conclusions: Screening for FGFR1 amplification in ctDNA might represent a viable strategy to identify patients eligible for treatment by FGFR inhibition, and mBC ctDNA levels might be used for the evaluation of prognosis in clinical drug trials.

Cell-Free DNA and Apoptosis: How Dead Cells Inform About the Living.

Cell-free DNA (cfDNA) is evolving into a widely used prognostic and predictive biomarker, particularly in oncology. However, its versatile clinical use precedes a profound understanding of the underlying biology of cfDNA release. There is much evidence to suggest that cfDNA is mainly derived from dying (i.e., apoptotic) cells. However, numerous cancer studies have shown that cfDNA is informative about acquired resistance to given therapies, which is present in living, proliferating tumor subclones. To explain this contradiction, we review current insights regarding cfDNA release, in particular the interplay between apoptosis and proliferation. We describe how improved knowledge about cfDNA biology could be used for novel therapeutic strategies and how this may affect patient management.

Publisher Correction: Inference of transcription factor binding from cell-free DNA enables tumor subtype prediction and early detection.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Cell-free DNA analysis reveals POLR1D-mediated resistance to bevacizumab in colorectal cancer.

BACKGROUND: Bevacizumab, a monoclonal antibody against soluble VEGFA, is an approved and commonly administered anti-angiogenic drug in patients with metastasized colorectal cancer (mCRC). The survival benefit of anti-VEGF therapy in mCRC patients is limited to a few months, and acquired resistance mechanisms are largely unknown. Here, we employed whole-genome sequencing of plasma DNA to evaluate the tumor genome of patients undergoing treatment with bevacizumab to determine novel aberrations associated with resistance. METHODS: Using longitudinal plasma analyses, we studied the evolution of tumor genomes in a mCRC cohort (n = 150) and conducted analyses of CRC cases from The Cancer Genome Atlas (TCGA) database (n = 619) to identify associations between genomic aberrations and clinical features. We employed whole-genome sequencing to identify the most frequently occurring focal somatic copy number alterations (SCNAs). Using the TCGA data as a comparative and supporting dataset, we defined the minimally amplified overlapping region and studied the mechanistic consequences of copy number gain of the involved genes in this segment. In addition, we established an in vitro cell model and conducted downstream gene expression and cell viability assays to confirm our findings from the patient dataset. RESULTS: We observed a recurrent focal amplification (8.7% of cases) on chromosome 13q12.2. Analysis of CRC cases from the TCGA database suggested that this amplicon is associated with more advanced stages. We confirmed that this 13q12.2 amplicon frequently emerges later during the clinical course of disease. After defining the minimally amplified region, we observed that the amplification and expression of one gene, POLR1D, impacted cell proliferation and resulted in upregulation of VEGFA, an important regulator of angiogenesis which has been implicated in the resistance to bevacizumab treatment. In fact, in several patients, we observed the emergence of this 13q12.2 amplicon under bevacizumab treatment, which was invariably associated with therapy resistance. CONCLUSIONS: Non-invasive analyses of cell-free DNA from patients undergoing treatment with bevacizumab enabled the tracking of evolving tumor genomes and helped identify a recurrent focal SCNA of clinical relevance. Here, we describe a novel resistance mechanism against a widely applied treatment in patients with mCRC which will impact the clinical management of patients.