Neuroblastoma, a rare embryonic tumor arising from neural crest development, is responsible for 15% of pediatric cancer-related deaths. Recently, several single-cell transcriptome studies were performed on neuroblastoma patient samples to investigate the cell of origin and tumor heterogeneity. However, these individual studies involved a small number of tumors and cells, limiting the conclusions that could be drawn. To overcome this limitation, we integrated seven single-cell or single-nucleus datasets into a harmonized cell atlas covering 362,991 cells across 61 patients. We use this atlas to decipher the transcriptional landscape of neuroblastoma at single-cell resolution, revealing associations between transcriptomic profiles and clinical outcomes within the tumor compartment. In addition, we characterize the complex immune-cell landscape and uncover considerable heterogeneity among tumor-associated macrophages. Finally, we showcase the utility of our atlas as a resource by expanding it with additional data and using it as a reference for data-driven cell-type annotation.

Esophageal adenocarcinoma (EAC) is an aggressive cancer characterized by a high risk of relapse post-surgery. Current follow-up methods (serum carcinoembryonic antigen detection and PET-CT) lack sensitivity and reliability, necessitating a novel approach. Analyzing cell-free DNA (cfDNA) from blood plasma emerges as a promising avenue. This study aims to evaluate the cost-effective and genome-wide cell-free reduced representation bisulfite sequencing (cfRRBS) method combined with computational deconvolution for effective disease monitoring in EAC patients. cfDNA methylation profiling with cfRRBS was performed on 162 blood plasma samples from 33 EAC cancer patients and 28 blood plasma samples from 20 healthy donors. The estimated tumor fraction for EAC patients at the time of diagnosis was significantly different from the healthy donor plasma samples (one-sided Wilcoxon rank-sum test: p-value = 0.032). Tumor fractions above 15% and focal gains/amplifications in MYC (chr8), KRAS (chr12), EGFR (chr7) and NOTCH2 (chr1) were observed in four samples of distinct patients at the time metastatic disease was detected. This study showed feasibility to estimate tumor fractions in blood plasma of EAC patients based on cfDNA methylation using cfRRBS and computational deconvolution. Nevertheless, in this study only cancer patients with evidence of metastatic disease show high tumor fractions and copy number alterations.

Abstract Pediatric central nervous system tumors remain challenging to diagnose. Imaging approaches do not provide sufficient detail to discriminate between different tumor types, while the histopathological examination of tumor tissue shows high inter-observer variability. Recent studies have demonstrated the accurate classification of central nervous system tumors based on the DNA methylation profile of a tumor biopsy. However, a brain biopsy holds significant risk of bleeding and damaging the surrounding tissues. Liquid biopsy approaches analyzing circulating tumor DNA show high potential as an alternative and less invasive tool to study the DNA methylation pattern of tumors. Here, we explore the potential of classifying pediatric brain tumors based on methylation profiling of the circulating cell-free DNA (cfDNA) in cerebrospinal fluid (CSF). For this proof-of-concept study, we collected cerebrospinal fluid samples from 19 pediatric brain cancer patients via a ventricular drain placed for reasons of increased intracranial pressure. Analyses on the cfDNA showed high variability of cfDNA quantities across patients ranging from levels below the limit of quantification to 40 ng cfDNA per milliliter of CSF. Classification based on methylation profiling of cfDNA from CSF was correct for 7 out of 20 samples in our cohort. Accurate results were mostly observed in samples of high quality, more specifically those with limited high molecular weight DNA contamination. Interestingly, we show that centrifugation of the CSF prior to processing increases the fraction of fragmented cfDNA to high molecular weight DNA. In addition, classification was mostly correct for samples with high tumoral cfDNA fraction as estimated by computational deconvolution (> 40%). In summary, analysis of cfDNA in the CSF shows potential as a tool for diagnosing pediatric nervous system tumors especially in patients with high levels of tumoral cfDNA in the CSF. Further optimization of the collection procedure, experimental workflow and bioinformatic approach is required to also allow classification for patients with low tumoral fractions in the CSF.

Currently, we cannot provide a conclusive diagnosis for 3% to 5% of people who are confronted with cancer. These patients have cancer of unknown primary (CUP), ie, a metastasized cancer for which the tissue of origin cannot be determined. Studies have shown that the DNA methylation profile is a unique "fingerprint" that can be used to classify tumors. Here we used cell-free reduced representation bisulfite sequencing (cfRRBS), a technique that allows us to identify the methylation profile starting from minimal amounts of highly fragmented DNA, for CUP diagnosis on formalin-fixed paraffin-embedded (FFPE) tissue and liquid biopsies. We collected 80 primary tumor FFPE samples covering 16 tumor entities together with 15 healthy plasma samples to use as a custom cfRRBS reference data set. Entity-specific methylation regions are defined for each entity to build a classifier based on nonnegative least squares deconvolution. This classification framework was tested on 30 FFPE, 19 plasma, and 40 pleural and peritoneal effusion samples of both known metastatic tumors and clinical CUPs for which pathological investigation finally resulted in a cancer diagnosis. Using this framework, 27 of 30 FFPE (all CUPs) and 16 of 19 plasma samples (10/13 CUPs) obtained an accurate diagnosis, with a minimal DNA input of 400 pg. Diagnosis of the 40 pleural and peritoneal effusion samples is possible in 9 of 27 samples with negative/inconclusive cytology (6/13 CUPs), showing that cell-free DNA (cfDNA) methylation profiling could complement routine cytologic analysis. However, a low "cfDNA - high-molecular weight DNA ratio" has a considerable impact on the prediction accuracy. Moreover, the accuracy improves significantly if the predicted tumor percentage is >7%. This proof-of-concept study shows the feasibility of using DNA methylation profiling on FFPE and liquid biopsy samples such as blood, ascites, and pleural effusions in a fast and affordable way. Our novel RRBS-based technique requires minimal DNA input, can be performed in <1 week, and is highly adaptable to specific diagnostic problems as we only use 5 FFPE references per tumor entity. We believe that cfRRBS methylation profiling could be a valuable addition to the pathologist's toolbox in the diagnosis of CUPs.

Transcriptomic profiling of blood immune cells offers a promising alternative to invasive, sampling bias-prone tissue-based biomarker assays for predicting immune checkpoint inhibitor (ICI) therapy response in non-small cell lung cancer (NSCLC) patients. However, the optimal analytical approach to identify systemic correlates of response still needs to be explored. We collected peripheral blood mononuclear cells and whole blood (WB) samples from 33 ICI-treated NSCLC patients before ICI treatment and at the first response evaluation. After bulk polyadenylated RNA-sequencing, we assessed differences in gene expression profiles between non-responders and responders using differential expression analysis, single sample gene set enrichment analysis (ssGSEA), and cell type deconvolution. We evaluated gene expression values, ssGSEA scores, and deconvolved cell type proportions to distinguish non-responders from responders via ROC curve (AUC) analysis, training a logistic regression classification model. Gene expression values and deconvolved proportions yielded the best results with WB samples after treatment (AUC = 0.87 and 0.85, respectively). Overall, ssGSEA scores showed superior classification performance across all sample types and timepoints (AUC > 0.7). In conclusion, transcriptomic analysis through ssGSEA demonstrated the best performance as a non-invasive biomarker for predicting clinical benefit in ICI-treated NSCLC patients, with gene expression and deconvolution on post-treatment WB samples also showing promising results.

It is generally assumed that extracellular long RNA molecules in biofluids are fragmented. Few studies have indirectly hinted at the existence of possibly functional, intact long RNA transcripts. In search for such RNA molecules, we developed a long-read full transcript sequencing workflow for low-input and low-quality samples. We applied our method to human blood plasma, urine, and their isolated macromolecular compartments, in parallel with total RNA sequencing. This approach enabled us to find intact messenger RNA molecules in human biofluids and macromolecular compartments. We showed that the full-length transcriptome of human urine and blood plasma differs, but we also reveal intact messenger RNA molecules shared between biofluids. In addition, we show that these intact molecules are differentially distributed over fractionated macromolecular compartments. This study provides a foundation for future extracellular RNA studies to elucidate the human biofluid full-length transcriptome. ### Competing Interest Statement The authors have declared no competing interest.

Blood plasma, one of the most studied liquid biopsies, contains various molecules that have biomarker potential for cancer detection, including cell-free DNA (cfDNA) and cell-free RNA (cfRNA). As the vast majority of cell-free nucleic acids in circulation are non-cancerous, a laboratory workflow with a high detection sensitivity of tumor-derived nucleic acids is a prerequisite for precision oncology. One way to meet this requirement is by the combined analysis of cfDNA and cfRNA from the same liquid biopsy sample. So far, no study has systematically compared the performance of cfDNA and cfRNA co-purification to increase sensitivity. First, we set up a framework using digital PCR (dPCR) technology to quantify cfDNA and cfRNA from human blood plasma in order to compare cfDNA/cfRNA co-purification kit performance. To that end, we optimized two dPCR duplex assays, designed to quantify both cfDNA and cfRNA with the same assays, by ensuring that primers and probes are located within a highly abundant exon. Next, we applied our optimized workflow to evaluate the co-purification performance of two manual and two semi-automated methods over a range of plasma input volumes (0.06–4 mL). Some kits result in higher nucleic acid concentrations in the eluate, while consuming only half of the plasma volume. The combined nucleic acid quantification systematically results in higher nucleic acid concentrations as compared to a parallel quantification of cfDNA and cfRNA in the eluate. We provide a framework to evaluate the performance of cfDNA/cfRNA co-purification kits and have tested two manual and two semi-automated co-purification kits in function of the available plasma input amount and the intended use of the nucleic acid eluate. We demonstrate that the combined quantification of cfDNA and cfRNA has a benefit compared to separate quantification. We foresee that the results of this study are instrumental for clinical applications to help increase mutation detection sensitivity, allowing improved disease detection and monitoring.