Computational deconvolution of DNA methylation data from mixed DNA samples.

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

A Fast, Affordable, and Minimally Invasive Diagnostic Test for Cancer of Unknown Primary Using DNA Methylation Profiling.

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

Digital PCR-based evaluation of nucleic acid extraction kit performance for the co-purification of cell-free DNA and RNA.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

miR-99b-5p, miR-380-3p, and miR-485-3p are novel chemosensitizing miRNAs in high-risk neuroblastoma.

Neuroblastoma is a deadly childhood cancer arising in the developing sympathetic nervous system. High-risk patients are currently treated with intensive chemotherapy, which is curative in only 50% of children and leaves some surviving patients with life-long side effects. microRNAs (miRNAs) are critical regulators of neural crest development and are deregulated during neuroblastoma tumorigenesis, making miRNA-based drugs an attractive therapeutic avenue. A functional screen of >1,200 miRNA mimics was conducted in neuroblastoma cell lines to discover miRNAs that sensitized cells to low doses (30% inhibitory concentration [IC30]) of doxorubicin and vincristine chemotherapy used in the treatment of the disease. Three miRNAs, miR-99b-5p, miR-380-3p, and miR-485-3p, had potent chemosensitizing activity with doxorubicin in multiple models of high-risk neuroblastoma. These miRNAs underwent genomic loss in a subset of neuroblastoma patients, and low expression predicted poor survival outcome. In vitro functional assays revealed each of these miRNAs enhanced the anti-proliferative and pro-apoptotic effects of doxorubicin. We used RNA sequencing (RNA-seq) to show that miR-99b-5p represses neuroblastoma dependency genes LIN28B and PHOX2B both in vitro and in patient-derived xenograft (PDX) tumors. Luciferase reporter assays demonstrate that PHOX2B is a direct target of miR-99b-5p. We anticipate that restoring the function of the tumor-suppressive miRNAs discovered here may be a valuable therapeutic strategy for the treatment of neuroblastoma patients.

Mitotic Dysregulation at Tumor Initiation Creates a Therapeutic Vulnerability to Combination Anti-Mitotic and Pro-Apoptotic Agents for MYCN-Driven Neuroblastoma.

MYCN amplification occurs in approximately 20-30% of neuroblastoma patients and correlates with poor prognosis. The TH-MYCN transgenic mouse model mimics the development of human high-risk neuroblastoma and provides strong evidence for the oncogenic function of MYCN. In this study, we identified mitotic dysregulation as a hallmark of tumor initiation in the pre-cancerous ganglia from TH-MYCN mice that persists through tumor progression. Single-cell quantitative-PCR of coeliac ganglia from 10-day-old TH-MYCN mice revealed overexpression of mitotic genes in a subpopulation of premalignant neuroblasts at a level similar to single cells derived from established tumors. Prophylactic treatment using antimitotic agents barasertib and vincristine significantly delayed the onset of tumor formation, reduced pre-malignant neuroblast hyperplasia, and prolonged survival in TH-MYCN mice. Analysis of human neuroblastoma tumor cohorts showed a strong correlation between dysregulated mitosis and features of MYCN amplification, such as MYC(N) transcriptional activity, poor overall survival, and other clinical predictors of aggressive disease. To explore the therapeutic potential of targeting mitotic dysregulation, we showed that genetic and chemical inhibition of mitosis led to selective cell death in neuroblastoma cell lines with MYCN over-expression. Moreover, combination therapy with antimitotic compounds and BCL2 inhibitors exploited mitotic stress induced by antimitotics and was synergistically toxic to neuroblastoma cell lines. These results collectively suggest that mitotic dysregulation is a key component of tumorigenesis in early neuroblasts, which can be inhibited by the combination of antimitotic compounds and pro-apoptotic compounds in MYCN-driven neuroblastoma.

Distinct Notch1 and BCL11B requirements mediate human γδ/αß T cell development.

γδ and αß T cells have unique roles in immunity and both originate in the thymus from T-lineage committed precursors through distinct but unclear mechanisms. Here, we show that Notch1 activation is more stringently required for human γδ development compared to αß-lineage differentiation and performed paired mRNA and miRNA profiling across 11 discrete developmental stages of human T cell development in an effort to identify the potential Notch1 downstream mechanism. Our data suggest that the miR-17-92 cluster is a Notch1 target in immature thymocytes and that miR-17 can restrict BCL11B expression in these Notch-dependent T cell precursors. We show that enforced miR-17 expression promotes human γδ T cell development and, consistently, that BCL11B is absolutely required for αß but less for γδ T cell development. This study suggests that human γδ T cell development is mediated by a stage-specific Notch-driven negative feedback loop through which miR-17 temporally restricts BCL11B expression and provides functional insights into the developmental role of the disease-associated genes BCL11B and the miR-17-92 cluster in a human context.

Recurrent chromosomal imbalances provide selective advantage to human embryonic stem cells under enhanced replicative stress conditions.

Human embryonic stem cells (hESCs) and embryonal tumors share a number of common features, including a compromised G1/S checkpoint. Consequently, these rapidly dividing hESCs and cancer cells undergo elevated levels of replicative stress, inducing genomic instability that drives chromosomal imbalances. In this context, it is of interest that long-term in vitro cultured hESCs exhibit a remarkable high incidence of segmental DNA copy number gains, some of which are also highly recurrent in certain malignancies such as 17q gain (17q+). The selective advantage of DNA copy number changes in these cells has been attributed to several underlying processes including enhanced proliferation. We hypothesized that these recurrent chromosomal imbalances become rapidly embedded in the cultured hESCs through a replicative stress driven Darwinian selection process. To this end, we compared the effect of hydroxyurea-induced replicative stress vs normal growth conditions in an equally mixed cell population of isogenic euploid and 17q + hESCs. We could show that 17q + hESCs rapidly overtook normal hESCs. Our data suggest that recurrent chromosomal segmental gains provide a proliferative advantage to hESCs under increased replicative stress, a process that may also explain the highly recurrent nature of certain imbalances in cancer.

Correction to: The pitfalls and promise of liquid biopsies for diagnosing and treating solid tumors in children: as review.

In the original version of this article, a reader pointed out that there was a mistake in the phrasing in a paragraph. This could potentially be harmful to children. The authors agree to change the wording. "vitreous fluid" will be changed to "aqueous humor".

The pitfalls and promise of liquid biopsies for diagnosing and treating solid tumors in children: a review.

Cell-free DNA profiling using patient blood is emerging as a non-invasive complementary technique for cancer genomic characterization. Since these liquid biopsies will soon be integrated into clinical trial protocols for pediatric cancer treatment, clinicians should be informed about potential applications and advantages but also weaknesses and potential pitfalls. Small retrospective studies comparing genetic alterations detected in liquid biopsies with tumor biopsies for pediatric solid tumor types are encouraging. Molecular detection of tumor markers in cell-free DNA could be used for earlier therapy response monitoring and residual disease detection as well as enabling detection of pathognomonic and therapeutically relevant genomic alterations.Conclusion: Existing analyses of liquid biopsies from children with solid tumors increasingly suggest a potential relevance for molecular diagnostics, prognostic assessment, and therapeutic decision-making. Gaps remain in the types of tumors studied and value of detection methods applied. Here we review the current stand of liquid biopsy studies for pediatric solid tumors with a dedicated focus on cell-free DNA analysis. There is legitimate hope that integrating fully validated liquid biopsy-based innovations into the standard of care will advance patient monitoring and personalized treatment of children battling solid cancers.What is Known:• Liquid biopsies are finding their way into routine oncological screening, diagnosis, and disease monitoring in adult cancer types fast.• The most widely adopted source for liquid biopsies is blood although other easily accessible body fluids, such as saliva, pleural effusions, urine, or cerebrospinal fluid (CSF) can also serve as sources for liquid biopsiesWhat is New:• Retrospective proof-of-concept studies in small cohorts illustrate that liquid biopsies in pediatric solid tumors yield tremendous potential to be used in diagnostics, for therapy response monitoring and in residual disease detection.• Liquid biopsy diagnostics could tackle some long-standing issues in the pediatric oncology field; they can enable accurate genetic diagnostics in previously unbiopsied tumor types like renal tumors or brain stem tumors leading to better treatment strategies.

Computational deconvolution of transcriptomics data from mixed cell populations.

Summary: Gene expression analyses of bulk tissues often ignore cell type composition as an important confounding factor, resulting in a loss of signal from lowly abundant cell types. In this review, we highlight the importance and value of computational deconvolution methods to infer the abundance of different cell types and/or cell type-specific expression profiles in heterogeneous samples without performing physical cell sorting. We also explain the various deconvolution scenarios, the mathematical approaches used to solve them and the effect of data processing and different confounding factors on the accuracy of the deconvolution results. Contact: katleen.depreter@ugent.be. Supplementary information: Supplementary data are available at Bioinformatics online.

Correction: Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

The original version of this Article contained an error in the spelling of the author Anja K. Mayer, which was incorrectly given as Anja Kathrin Mayer. This has now been corrected in both the PDF and HTML versions of the Article.

Promoter-associated proteins of EPAS1 identified by enChIP-MS - A putative role of HDX as a negative regulator.

Presence of perivascular neuroblastoma cells with high expression of hypoxia inducible factor (HIF)-2α correlates with distant metastasis and aggressive disease. Regulation of HIFs are traditionally considered to occur post-translationally, but we have recently shown that HIF-2α is unconventionally regulated also at the transcriptional level in neuroblastoma cells. Regulatory factors binding directly to EPAS1 (encoding HIF-2α) to promote transcription are yet to be defined. Here, we employ the novel CRISPR/Cas9-based engineered DNA-binding molecule-mediated chromatin immunoprecipitation (enChIP) - mass spectrometry (MS) methodology to, in an unbiased fashion, identify proteins that associate with the EPAS1 promoter under normoxic and hypoxic conditions. Our enChIP analysis resulted in 27 proteins binding to the EPAS1 promoter in neuroblastoma cells. In agreement with a general hypoxia-driven downregulation of gene transcription, the majority (24 out of 27) of proteins dissociate from the promoter at hypoxia. Among them were several nucleosome-associated proteins suggesting a general opening of chromatin as one explanation to induced EPAS1 transcription at hypoxia. Of particular interest from the list of released factors at hypoxia was the highly divergent homeobox (HDX) transcription factor, that we show inversely correlates with HIF-2α in neuroblastoma cells. We propose a putative model where HDX negatively regulates EPAS1 expression through a release-of-inhibition mechanism.

Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

PurposePart of the hidden genetic variation in heterogeneous genetic conditions such as inherited retinal diseases (IRDs) can be explained by copy-number variations (CNVs). Here, we explored the genomic landscape of IRD genes listed in RetNet to identify and prioritize those genes susceptible to CNV formation.MethodsRetNet genes underwent an assessment of genomic features and of CNV occurrence in the Database of Genomic Variants and literature. CNVs identified in an IRD cohort were characterized using targeted locus amplification (TLA) on extracted genomic DNA.ResultsExhaustive literature mining revealed 1,345 reported CNVs in 81 different IRD genes. Correlation analysis between rankings of genomic features and CNV occurrence demonstrated the strongest correlation between gene size and CNV occurrence of IRD genes. Moreover, we identified and delineated 30 new CNVs in IRD cases, 13 of which are novel and three of which affect noncoding, putative cis-regulatory regions. Finally, the breakpoints of six complex CNVs were determined using TLA in a hypothesis-neutral manner.ConclusionWe propose a ranking of CNV-prone IRD genes and demonstrate the efficacy of TLA for the characterization of CNVs on extracted DNA. Finally, this IRD-oriented CNV study can serve as a paradigm for other genetically heterogeneous Mendelian diseases with hidden genetic variation.

Zipper plot: visualizing transcriptional activity of genomic regions.

BACKGROUND: Reconstructing transcript models from RNA-sequencing (RNA-seq) data and establishing these as independent transcriptional units can be a challenging task. Current state-of-the-art tools for long non-coding RNA (lncRNA) annotation are mainly based on evolutionary constraints, which may result in false negatives due to the overall limited conservation of lncRNAs. RESULTS: To tackle this problem we have developed the Zipper plot, a novel visualization and analysis method that enables users to simultaneously interrogate thousands of human putative transcription start sites (TSSs) in relation to various features that are indicative for transcriptional activity. These include publicly available CAGE-sequencing, ChIP-sequencing and DNase-sequencing datasets. Our method only requires three tab-separated fields (chromosome, genomic coordinate of the TSS and strand) as input and generates a report that includes a detailed summary table, a Zipper plot and several statistics derived from this plot. CONCLUSION: Using the Zipper plot, we found evidence of transcription for a set of well-characterized lncRNAs and observed that fewer mono-exonic lncRNAs have CAGE peaks overlapping with their TSSs compared to multi-exonic lncRNAs. Using publicly available RNA-seq data, we found more than one hundred cases where junction reads connected protein-coding gene exons with a downstream mono-exonic lncRNA, revealing the need for a careful evaluation of lncRNA 5'-boundaries. Our method is implemented using the statistical programming language R and is freely available as a webtool.

The immunophenotypic fingerprint of patients with primary antibody deficiencies is partially present in their asymptomatic first-degree relatives.

The etiology of primary antibody deficiencies is largely unknown. Beside rare monogenic forms, the majority of cases seem to have a more complex genetic basis. Whereas common variable immunodeficiency has been investigated in depth, there are only a few reports on milder primary antibody deficiencies such as idiopathic primary hypogammaglobulinemia and IgG subclass deficiency. We performed flow cytometric immunophenotyping in 33 patients with common variable immunodeficiency, 23 with idiopathic primary hypogammaglobulinemia and 21 with IgG subclass deficiency, as well as in 47 asymptomatic first-degree family members of patients and 101 unrelated healthy controls. All three groups of patients showed decreased memory B- and naïve T-cell subsets and decreased B-cell activating factor receptor expression. In contrast, circulating follicular helper T-cell frequency and expression of inducible T-cell co-stimulator and chemokine receptors were only significantly altered in patients with common variable immunodeficiency. Asymptomatic first-degree family members of patients demonstrated similar, albeit intermediate, alterations in naïve and memory B- and T-cell subsets. About 13% of asymptomatic relatives had an abnormal peripheral B-cell composition. Furthermore, asymptomatic relatives showed decreased levels of CD4+ recent thymic emigrants and increased central memory T cells. Serum IgG and IgM levels were also significantly lower in asymptomatic relatives than in healthy controls. We conclude that, in our cohort, the immunophenotypic landscape of primary antibody deficiencies comprises a spectrum, in which some alterations are shared between all primary antibody deficiencies whereas others are only associated with common variable immunodeficiency. Importantly, asymptomatic first-degree family members of patients were found to have an intermediate phenotype for peripheral B- and T-cell subsets.

Asthma inflammatory phenotypes show differential microRNA expression in sputum.

BACKGROUND: Asthma is classified according to severity and inflammatory phenotype and is likely to be distinguished by specific microRNA (miRNA) expression profiles. OBJECTIVE: We sought to associate miRNA expression in sputum supernatants with the inflammatory cell profile and disease severity in asthmatic patients. METHODS: We investigated miRNA expression in sputum supernatants of 10 healthy subjects, 17 patients with mild-to-moderate asthma, and 9 patients with severe asthma using high-throughput, stem-loop, reverse transcriptase quantitative real-time PCR miRNA expression profiling (screening cohort, n = 36). Differentially expressed miRNAs were validated in an independent cohort (n = 60; 10 healthy subjects and 50 asthmatic patients). Cellular miRNA origin was examined by using in situ hybridization and reverse transcriptase quantitative real-time PCR. The functional role of miRNAs was assessed by using in silico analysis and in vitro transfecting miRNA mimics in human bronchial epithelial cells. RESULTS: In 2 independent cohorts expression of miR-629-3p, miR-223-3p, and miR-142-3p was significantly upregulated in sputum of patients with severe asthma compared with that in healthy control subjects and was highest in patients with neutrophilic asthma. Expression of the 3 miRNAs was associated with sputum neutrophilia, and miR-223-3p and miR-142-3p expression was associated also with airway obstruction (FEV1/forced vital capacity). Expression of miR-629-3p was localized in the bronchial epithelium, whereas miR-223-3p and miR-142-3p were expressed in neutrophils, monocytes, and macrophages. Transfecting human bronchial epithelial cells with miR-629-3p mimic induced epithelial IL-8 mRNA and protein expression. IL-1ß and IL-8 protein levels were significantly increased in sputum of patients with severe asthma and were positively associated with sputum neutrophilia. CONCLUSIONS: Expression of miR-223-3p, miR-142-3p, and miR-629-3p is increased in sputum of patients with severe asthma and is linked to neutrophilic airway inflammation, suggesting that these miRNAs contribute to this asthma inflammatory phenotype.

Individual patient risk stratification of high-risk neuroblastomas using a two-gene score suited for clinical use.

Several gene expression-based prognostic signatures have been described in neuroblastoma, but none have successfully been applied in the clinic. Here we have developed a clinically applicable prognostic gene signature, both with regards to number of genes and analysis platform. Importantly, it does not require comparison between patients and is applicable amongst high-risk patients. The signature is based on a two-gene score (R-score) with prognostic power in high-stage tumours (stage 4 and/or MYCN-amplified diagnosed after 18 months of age). QPCR-based and array-based analyses of matched cDNAs confirmed cross platform (array-qPCR) transferability. We also defined a fixed cut-off value identifying prognostically differing subsets of high-risk patients on an individual patient basis. This gene expression signature independently contributes to the current neuroblastoma classification system, and if prospectively validated could provide further stratification of high-risk patients, and potential upfront identification of a group of patients that are in need of new/additional treatment regimens.

MYCN and HDAC2 cooperate to repress miR-183 signaling in neuroblastoma.

MYCN is a master regulator controlling many processes necessary for tumor cell survival. Here, we unravel a microRNA network that causes tumor suppressive effects in MYCN-amplified neuroblastoma cells. In profiling studies, histone deacetylase (HDAC) inhibitor treatment most strongly induced miR-183. Enforced miR-183 expression triggered apoptosis, and inhibited anchorage-independent colony formation in vitro and xenograft growth in mice. Furthermore, the mechanism of miR-183 induction was found to contribute to the cell death phenotype induced by HDAC inhibitors. Experiments to identify the HDAC(s) involved in miR-183 transcriptional regulation showed that HDAC2 depletion induced miR-183. HDAC2 overexpression reduced miR-183 levels and counteracted the induction caused by HDAC2 depletion or HDAC inhibitor treatment. MYCN was found to recruit HDAC2 in the same complexes to the miR-183 promoter, and HDAC2 depletion enhanced promoter-associated histone H4 pan-acetylation, suggesting epigenetic changes preceded transcriptional activation. These data reveal miR-183 tumor suppressive properties in neuroblastoma that are jointly repressed by MYCN and HDAC2, and suggest a novel way to bypass MYCN function.