Background DNA methylation detection in liquid biopsies provides a highly promising and much needed means for real-time monitoring of disease load in advanced cancer patient care. Compared to the often-used somatic mutations, tissue- and cancer-type specific epigenetic marks affect a larger part of the cancer genome and generally have a high penetrance throughout the tumour. Here, we describe the successful application of the recently described MeD-seq assay for genome-wide DNA methylation profiling on cell-free DNA (cfDNA). The compatibility of the MeD-seq assay with different types of blood collection tubes, cfDNA input amounts, cfDNA isolation methods, and vacuum concentration of samples was evaluated using plasma from both metastatic cancer patients and healthy blood donors (HBDs). To investigate the potential value of cfDNA methylation profiling for tumour load monitoring, we profiled paired samples from 8 patients with resectable colorectal liver metastases (CRLM) before and after surgery. Results The MeD-seq assay worked on plasma-derived cfDNA from both EDTA and CellSave blood collection tubes when at least 10 ng of cfDNA was used. From the 3 evaluated cfDNA isolation methods, both the manual QIAamp Circulating Nucleic Acid Kit (Qiagen) and the semi-automated Maxwell® RSC ccfDNA Plasma Kit (Promega) were compatible with MeD-seq analysis, whereas the QiaSymphony DSP Circulating DNA Kit (Qiagen) yielded significantly fewer reads when compared to the QIAamp kit (p < 0.001). Vacuum concentration of samples before MeD-seq analysis was possible with samples in AVE buffer (QIAamp) or water, but yielded inconsistent results for samples in EDTA-containing Maxwell buffer. Principal component analysis showed that pre-surgical samples from CRLM patients were very distinct from HBDs, whereas post-surgical samples were more similar. Several described methylation markers for colorectal cancer monitoring in liquid biopsies showed differential methylation between pre-surgical CRLM samples and HBDs in our data, supporting the validity of our approach. Results for MSC, ITGA4, GRIA4, and EYA4 were validated by quantitative methylation specific PCR. Conclusions The MeD-seq assay provides a promising new method for cfDNA methylation profiling. Potential future applications of the assay include marker discovery specifically for liquid biopsy analysis as well as direct use as a disease load monitoring tool in advanced cancer patients.

Background: Biomarkers predicting treatment response may be used to stratify pancreatic ductal adenocarcinoma (PDAC) patients for therapy. The aim of this study was to identify circulating tumor DNA (ctDNA) mutations that associate with tumor progression during FOLFIRINOX chemotherapy, and overall survival (OS). Methods: Circulating cell-free DNA was analyzed with a 57 gene next-generation sequencing panel using plasma samples of 48 PDAC patients of all disease stages. Patients received FOLFIRINOX as initial treatment. Chemotherapy response was determined on CT scans as disease control ( n = 30) or progressive disease ( n = 18) within eight cycles of FOLFIRINOX, based on RECIST 1.1 criteria. Results: Detection of a TP53 ctDNA mutation before start of FOLFIRINOX [odds ratio (OR) 10.51, 95% confidence interval (CI) 1.40–79.14] and the presence of a homozygous TP53 Pro72Arg germline variant (OR 6.98, 95% CI 1.31–37.30) were predictors of early tumor progression during FOLFIRINOX in multivariable analysis. Five patients presented with the combination of a TP53 ctDNA mutation before start of FOLFIRINOX and the homozygous Pro72Arg variant. All five patients showed progression during FOLFIRINOX. The combination of the TP53 mutation and TP53 germline variant was associated with shorter survival (median OS 4.4 months, 95% CI 2.6–6.2 months) compared with patients without any TP53 alterations (median OS 13.0 months, 95% CI 8.6–17.4 months). Conclusion: The combination of a TP53 ctDNA mutation before start of FOLFIRINOX and a homozygous TP53 Pro72Arg variant is a promising biomarker, associated with early tumor progression during FOLFIRINOX and poor OS. The results of this exploratory study need to be validated in an independent cohort.

The BCL11A gene encodes a transcriptional repressor with essential functions in multiple tissues during human development. Haploinsufficiency for BCL11A causes Dias-Logan syndrome (OMIM 617101), an intellectual developmental disorder with hereditary persistence of fetal hemoglobin (HPFH). Due to the severe phenotype, disease-causing variants in BCL11A occur de novo. We describe a patient with a de novo heterozygous variant, c.1453G&gt;T, in the BCL11A gene, resulting in truncation of the BCL11A-XL protein (p.Glu485X). The truncated protein lacks the 3 C-terminal DNA-binding zinc fingers and the nuclear localization signal, rendering it inactive. The patient displayed high fetal hemoglobin (HbF) levels (12.1-18.7% of total hemoglobin), in contrast to the parents who had HbF levels of 0.3%. We used cultures of patient-derived erythroid progenitors to determine changes in gene expression and chromatin accessibility. In addition, we investigated DNA methylation of the promoters of the γ-globin genes HBG1 and HBG2. HUDEP1 and HUDEP2 cells were used as models for fetal and adult human erythropoiesis, respectively. Similar to HUDEP1 cells, the patient’s cells displayed Assay for Transposase-Accessible Chromatin (ATAC) peaks at the HBG1/2 promoters and significant expression of HBG1/2 genes. In contrast, HBG1/2 promoter methylation and genome-wide gene expression profiling were consistent with normal adult erythropoiesis. We conclude that HPFH is the major erythroid phenotype of constitutive BCL11A haploinsufficiency. Given the essential functions of BCL11A in other hematopoietic lineages and the neuronal system, erythroid-specific targeting of the BCL11A gene has been proposed for reactivation of γ-globin expression in β-hemoglobinopathy patients. Our data strongly support this approach.

Transcriptomics can be combined with TRA and TRB clonotype analysis at the single cell level. The aim of this study was to validate this approach on the ICELL8 Single-Cell system and to evaluate its usefulness to analyse clinical paucicellular samples. For this purpose, we carefully selected T cell lines with defined TRA/TRB clonotypes as well as clinical samples enriched for CD3 + T cells that possess a complex TCR repertoire. Low cell numbers of the different samples were dispensed in a chip on the ICELL8 Single-Cell System. Two sequencing libraries were generated from each single cell cDNA preparation, one for the TRA/TRB repertoire and one for the 5 ′ ends of transcripts, and subsequently sequenced. Transcriptome analysis revealed that the cell lines on average express 2,268 unique genes/cell and T cells of clinical samples 770 unique genes/cell. The expected combined TRA/TRB clonotype was determined for on average 71% of the cells of the cell lines. In the clinical samples the TRA/TRB repertoire was more complex than those of the cell lines. Furthermore, the TRB clonotype distribution of the clinical samples was positively correlated to frequencies of TCRVβ families in CD3 + T cells obtained by a flow cytometry-based approach (Spearman's Rho correlation coefficient 0.81, P = 6.49 * 10 −7). Combined analyses showed that transcriptome-based cell type-specific clusters in clinical samples corresponded to clinical features such as CMV status. In conclusion, we showed that the ICELL8 Single-Cell System enabled combined interrogation of both TRA/TRB repertoire and transcriptome of paucicellular clinical samples. This opens the way to study the response of single T cells within heterogeneous samples for both their transcriptome and TRA/TRB clonotypes in disease or upon treatment.