The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons.

BACKGROUND: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. RESULTS: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. CONCLUSIONS: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Gene Expression Differences between Ductal Carcinoma in Situ with and without Progression to Invasive Breast Cancer.

To understand the molecular alterations driving the progression of ductal carcinoma in situ (DCIS), we compared patients with pure DCIS and patients with DCIS and synchronous invasive breast cancer (IBC). Twelve patients with extensive pure DCIS were included as a representation of indolent lesions with limited invasive capacity. These cases were matched with 12 patients with a limited DCIS component and IBC, representing lesions with a high invasive potential. Matching included age and surrogate DCIS subtypes. Gene expression profiling was performed on DCIS cells to identify transcriptional differences between these two groups. The identified genes were validated by immunohistochemistry. Nine genes showed significantly different expression. Most of these genes were highly expressed in DCIS samples with IBC, including PLAU (P = 0.002), COL1A1 (P = 0.006), KRT81 (P = 0.009), S100A7 (P = 0.015), SCGB1D2 (P = 0.023), KRT18 (P = 0.029), and NOTCH3 (P = 0.044), whereas EGFR and CXCL14 showed a higher expression in cases with pure DCIS (P = 0.015 and P = 0.028, respectively). This difference was only significant for SCGB1D2 (P = 0.009). Hierarchical clustering revealed distinct clustering of patients with and without invasion. Patients with pure DCIS have a different gene expression pattern as compared to patients with DCIS and synchronous IBC. These genes may pinpoint to driver pathway(s) that play an important role in DCIS progression.

T lymphocytes facilitate brain metastasis of breast cancer by inducing Guanylate-Binding Protein 1 expression.

The discovery of genes and molecular pathways involved in the formation of brain metastasis would direct the development of therapeutic strategies to prevent this deadly complication of cancer. By comparing gene expression profiles of Estrogen Receptor negative (ER-) primary breast tumors between patients who developed metastasis to brain and to organs other than brain, we found that T lymphocytes promote the formation of brain metastases. To functionally test the ability of T cells to promote brain metastasis, we used an in vitro blood-brain barrier (BBB) model. By co-culturing T lymphocytes with breast cancer cells, we confirmed that T cells increase the ability of breast cancer cells to cross the BBB. Proteomics analysis of the tumor cells revealed Guanylate-Binding Protein 1 (GBP1) as a key T lymphocyte-induced protein that enables breast cancer cells to cross the BBB. The GBP1 gene appeared to be up-regulated in breast cancer of patients who developed brain metastasis. Silencing of GBP1 reduced the ability of breast cancer cells to cross the in vitro BBB model. In addition, the findings were confirmed in vivo in an immunocompetent syngeneic mouse model. Co-culturing of ErbB2 tumor cells with activated T cells induced a significant increase in Gbp1 expression by the cancer cells. Intracardial inoculation of the co-cultured tumor cells resulted in preferential seeding to brain. Moreover, intracerebral outgrowth of the tumor cells was demonstrated. The findings point to a role of T cells in the formation of brain metastases in ER- breast cancers, and provide potential targets for intervention to prevent the development of cerebral metastases.

Integrative analysis of genomics and proteomics data on clinical breast cancer tissue specimens extracted with acid guanidinium thiocyanate-phenol-chloroform.

Acid guanidinium thiocyanate, phenol, and chloroform extraction (AGPC) is a commonly used procedure to extract RNA from fresh frozen tissues and cell lines. In addition, DNA and proteins can be recovered, which makes AGPC an attractive source for integrative analysis on tissues of which little material is available, such as clinical specimens. Despite this potential, AGPC has only scarcely been used for proteomic analysis, mainly due to difficulties in extracting proteins. We have used a quantitative mass spectrometry method to show that proteins can readily be recovered from AGPC extracted tissues with high recovery and repeatability, which allows this method to be used for global proteomic profiling. Protein expression data for a selected number of clinically relevant markers, of which transcript and protein levels are known to be correlated, were in agreement with genomic and transcriptomic data obtained from the same AGPC lysate. Furthermore, global proteomic profiling successfully discriminated breast tumor tissues according to their clinical subtype. Lastly, a reference gene set of differentially expressed transcripts was strongly enriched in the differentially abundant proteins in our cohort. AGPC lysates are therefore well suited for comparative protein and integrative analyses.

DC-SCRIPT is a novel regulator of the tumor suppressor gene CDKN2B and induces cell cycle arrest in ERα-positive breast cancer cells.

Breast cancer is one of the most common causes of cancer-related deaths in women. The estrogen receptor (ERα) is well known for having growth promoting effects in breast cancer. Recently, we have identified DC-SCRIPT (ZNF366) as a co-suppressor of ERα and as a strong and independent prognostic marker in ESR1 (ERα gene)-positive breast cancer patients. In this study, we further investigated the molecular mechanism on how DC-SCRIPT inhibits breast cancer cell growth. DC-SCRIPT mRNA levels from 190 primary ESR1-positive breast tumors were related to global gene expression, followed by gene ontology and pathway analysis. The effect of DC-SCRIPT on breast cancer cell growth and cell cycle arrest was investigated using novel DC-SCRIPT-inducible MCF7 breast cancer cell lines. Genome-wide expression profiling of DC-SCRIPT-expressing MCF7 cells was performed to investigate the effect of DC-SCRIPT on cell cycle-related gene expression. Findings were validated by real-time PCR in a cohort of 1,132 ESR1-positive breast cancer patients. In the primary ESR1-positive breast tumors, DC-SCRIPT expression negatively correlated with several cell cycle gene ontologies and pathways. DC-SCRIPT expression strongly reduced breast cancer cell growth in vitro, breast tumor growth in vivo, and induced cell cycle arrest. In addition, in the presence of DC-SCRIPT, multiple cell cycles related genes were differentially expressed including the tumor suppressor gene CDKN2B. Moreover, in 1,132 primary ESR1-positive breast tumors, DC-SCRIPT expression also correlated with CDKN2B expression. Collectively, these data show that DC-SCRIPT acts as a novel regulator of CDKN2B and induces cell cycle arrest in ESR1-positive breast cancer cells.

KRAS and BRAF mutation status in circulating colorectal tumor cells and their correlation with primary and metastatic tumor tissue.

Although anti-EGFR therapy has established efficacy in metastatic colorectal cancer, only 10-20% of unselected patients respond. This is partly due to KRAS and BRAF mutations, which are currently assessed in the primary tumor. To improve patient selection, assessing mutation status in circulating tumor cells (CTCs), which possibly better represent metastases than the primary tumor, could be advantageous. We investigated the feasibility of KRAS and BRAF mutation detection in colorectal CTCs by comparing three sensitive methods and compared mutation status in matching primary tumor, liver metastasis and CTCs. CTCs were isolated from blood drawn from 49 patients before liver resection using CellSearch™. DNA and RNA was isolated from primary tumors, metastases and CTCs. Mutations were assessed by co-amplification at lower denaturation temperature-PCR (Transgenomic™), real-time PCR (EntroGen™) and nested Allele-Specific Blocker (ASB-)PCR and confirmed by Sanger sequencing. In 43 of the 49 patients, tissue RNA and DNA was of sufficient quantity and quality. In these 43 patients, discordance between primary and metastatic tumor was 23% for KRAS and 7% for BRAF mutations. RNA and DNA from CTCs was available from 42 of the 43 patients, in which ASB-PCR was able to detect the most mutations. Inconclusive results in patients with low CTC counts limited the interpretation of discrepancies between tissue and CTCs. Determination of KRAS and BRAF mutations in CTCs is challenging but feasible. Of the tested methods, nested ASB-PCR, enabling detection of KRAS and BRAF mutations in patients with as little as two CTCs, seems to be superior.

Tumor-agnostic ctDNA levels by mFAST-SeqS in first-line HR-positive, HER2 negative metastatic breast cancer patients as a biomarker for survival.

This prospective cohort study reports aneuploidy score by mFast-SeqS as a strong prognostic marker in MBC patients. mFAST-SeqS is an affordable and easily implementable method for the assessment of total ctDNA levels and, as such, provides an alternative prognostic tool. One mixed cohort (cohort A, n = 45) starting any type of treatment in any line of therapy and one larger cohort (cohort B, n = 129) consisting of patients starting aromatase inhibitors (AI) as first-line therapy were used. mFAST-SeqS was performed using plasma of blood in which CTCs (CellSearch) were enumerated. The resulting aneuploidy score was correlated with categorized CTC count and associated with outcome. The aneuploidy score was significantly correlated with CTC count, but discordance was observed in 31.6% when applying cut-offs of 5. In both cohorts, aneuploidy score was a significant prognostic marker for both PFS and OS. In the Cox regression models, the HR for aneuploidy score for PFS was 2.52 (95% CI: 1.56-4.07), and the HR for OS was 2.37 (95% CI: 1.36-4.14). Results presented here warrant further investigations into the clinical utility of this marker in MBC patients.

mFast-SeqS-based aneuploidy score in circulating cell-free DNA is a prognostic biomarker in prostate cancer.

Multiple prognostic biomarkers, including circulating tumour cell (CTC) counts, exist in metastatic castration-resistant prostate cancer (mCRPC) patients, but none of them have been implemented into daily clinical care. The modified fast aneuploidy screening test-sequencing system (mFast-SeqS), which yields a genome-wide aneuploidy score, is able to reflect the fraction of cell-free tumour DNA (ctDNA) within cell-free DNA (cfDNA) and may be a promising biomarker in mCRPC. In this study, we investigated the prognostic value of dichotomized aneuploidy scores (< 5 vs. ≥ 5) as well as CTC counts (< 5 vs. ≥ 5) in 131 mCRPC patients prior to treatment with cabazitaxel. We validated our findings in an independent cohort of 50 similarly treated mCRPC patients. We observed that, similar to the dichotomized CTC count [HR: 2.92; 95% confidence interval (CI);1.84-4.62], dichotomized aneuploidy scores (HR: 3.24; CI: 2.12-4.94) significantly correlated with overall survival in mCRPC patients. We conclude that a dichotomized aneuploidy score from cfDNA is a prognostic marker for survival in mCRPC patients within our discovery cohort and in an independent mCRPC validation cohort. Therefore, this easy and robust minimally-invasive assay can be readily implemented as a prognostic marker in mCRPC. A dichotomized aneuploidy score might also be used as a stratification factor in clinical studies to account for tumour load.

Machine learning-based somatic variant calling in cell-free DNA of metastatic breast cancer patients using large NGS panels.

Next generation sequencing of cell-free DNA (cfDNA) is a promising method for treatment monitoring and therapy selection in metastatic breast cancer (MBC). However, distinguishing tumor-specific variants from sequencing artefacts and germline variation with low false discovery rate is challenging when using large targeted sequencing panels covering many tumor suppressor genes. To address this, we built a machine learning model to remove false positive variant calls and augmented it with additional filters to ensure selection of tumor-derived variants. We used cfDNA of 70 MBC patients profiled with both the small targeted Oncomine breast panel (Thermofisher) and the much larger Qiaseq Human Breast Cancer Panel (Qiagen). The model was trained on the panels' common regions using Oncomine hotspot mutations as ground truth. Applied to Qiaseq data, it achieved 35% sensitivity and 36% precision, outperforming basic filtering. For 20 patients we used germline DNA to filter for somatic variants and obtained 245 variants in total, while our model found seven variants, of which six were also detected using the germline strategy. In ten tumor-free individuals, our method detected in total one (potentially germline) variant, in contrast to 521 variants detected without our model. These results indicate that our model largely detects somatic variants.

Liquid Biopsies for Early Response Evaluation of Radium-223 in Metastatic Prostate Cancer.

PURPOSE: Reliable biomarkers for response monitoring during radium-223 treatment in patients with metastatic castration-resistant prostate cancer (mCRPC) are lacking. Circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA), obtained from liquid biopsies, are shown to have prognostic value in mCRPC. The aim of this study was to determine the value of CTCs and ctDNA for response evaluation of radium-223. METHODS: In this prospective multicenter study, longitudinal blood draws and imaging were performed in 28 patients with mCRPC and predominantly bone disease, who were treated with radium-223. CTCs were counted (CELLSEARCH CTC test), while fraction of ctDNA was estimated by measuring aneuploidy of cell-free DNA (cfDNA; modified Fast Aneuploidy Screening Test-Sequencing System). CTC counts and aneuploidy score (AS) were categorized as low (<5) and high (≥5). Primary and secondary clinical end points were failure-free survival (FFS), and overall survival (OS) and development of extraosseous metastases, respectively. Additionally, CTC count and AS were related to alkaline phosphatase (ALP) and total tumor volume in bone (TTVbone) on positron emission tomography-computed tomography with 68gallium prostate-specific membrane antigen. RESULTS: FFS was longer in patients with a low CTC count or AS either at baseline or after 12 weeks, whereas for OS, only a significant association with CTC count was observed. Liquid biopsy results correlated well with ALP and TTVbone at baseline, but not with change in both parameters after three cycles of radium-223. AS and CTC count were significantly correlated. CONCLUSION: CTC count and AS of cfDNA at baseline and during treatment predict clinical response to radium-223 in patients with mCRPC, warranting future evaluation of their value in treatment guidance.

Correlation of breast cancer susceptibility loci with patient characteristics, metastasis-free survival, and mRNA expression of the nearest genes.

To understand the biology of low-risk breast cancer alleles, and to investigate whether these loci also contribute to disease progression that was once established, we examined the association of SNPs tagging the low-risk breast cancer loci in or near FGFR2, LSP1, MAP3K1, H19, TOX3, POU5F1P1, MYC, and 2q35, with clinical, pathological characteristics, prognosis, and mRNA expression of the nearest genes. Tumor DNA samples of 2,480 breast cancer patients were available. Out of this cohort, 1,290 patients with lymph-node negative disease who did not receive adjuvant systemic therapy, the SNP status was associated with metastasis-free survival (MFS). In 1,401 patients, the mRNA expression levels of FGFR2, LSP1, MAP3K1, H19, TOX3, POU5F1P1, and MYC were determined and correlated with SNP genotypes. The SNP rs2981582 in FGFR2 was significantly associated with positive ER and PgR status (P < 0.001 and P = 0.003, respectively). No other significant associations with patient or tumor characteristics were observed. Only rs2107425 near H19 was significantly associated with shorter MFS in uni- and multi-variate analysis (HR: 1.53, CI: 1.12-2.08, P = 0.006 and HR: 1.59, CI: 1.16-2.20, P = 0.004, respectively), with the more aggressive minor allele displaying a recessive trait. The minor allele of SNP rs3803662 located near the TOX3 gene was associated with lower mRNA expression of this gene. In conclusion, except for the association of rs13283662 with TOX3 gene expression indicating a tumor suppressor role of TOX3, our findings suggest that breast cancer low-risk loci generally do not affect expression of the nearest gene in breast tumor tissue. Also the prognosis of patients is largely not affected by low-risk breast cancer loci except for the SNP near H19. How, this SNP affects prognosis warrants further study as it does not operate through altering H19 mRNA expression.

Gene expression profiling assigns CHEK2 1100delC breast cancers to the luminal intrinsic subtypes.

CHEK2 1100delC is a moderate-risk cancer susceptibility allele that confers a high breast cancer risk in a polygenic setting. Gene expression profiling of CHEK2 1100delC breast cancers may reveal clues to the nature of the polygenic CHEK2 model and its genes involved. Here, we report global gene expression profiles of a cohort of 155 familial breast cancers, including 26 CHEK2 1100delC mutant tumors. In line with previous work, all CHEK2 1100delC mutant tumors clustered among the hormone receptor-positive breast cancers. In the hormone receptor-positive subset, a 40-gene CHEK2 signature was subsequently defined that significantly associated with CHEK2 1100delC breast cancers. The identification of a CHEK2 gene signature implies an unexpected biological homogeneity among the CHEK2 1100delC breast cancers. In addition, all 26 CHEK2 1100delC tumors classified as luminal intrinsic subtype breast cancers, with 8 luminal A and 18 luminal B tumors. This biological make-up of CHEK2 1100delC breast cancers suggests that a relatively limited number of additional susceptibility alleles are involved in the polygenic CHEK2 model. Identification of these as-yet-unknown susceptibility alleles should be aided by clues from the 40-gene CHEK2 signature.

Genome-wide aneuploidy detected by mFast-SeqS in circulating cell-free DNA is associated with poor response to pembrolizumab in patients with advanced urothelial cancer.

Second-line treatment with immune checkpoint inhibition in patients with metastatic urothelial cancer (mUC) has a low success rate (~ 20%). Circulating tumour-derived DNA (ctDNA) levels may guide patient stratification, provided that an affordable and robust assay is available. Here, we investigate whether the modified fast aneuploidy screening test-sequencing system (mFast-SeqS) may provide such an assay. To this end, mFast-SeqS was performed on cell-free DNA (cfDNA) from 74 patients with mUC prior to treatment with pembrolizumab. Results were associated with corresponding tissue-based profiles, plasma-based variant allele frequencies (VAFs) and clinical response. We found that plasma-derived mFast-SeqS-based aneuploidy scores significantly correlated with those observed in the corresponding tumour tissue as well as with the ctDNA level in the plasma. In multivariate logistic regression analysis, a high aneuploidy score was independently associated with lack of clinical benefit from treatment with pembrolizumab. In conclusion, mFast-SeqS provides a patient-friendly, high-throughput and affordable method to estimate ctDNA level. Following independent validation, this test could be used to stratify mUC patients for response prior to the initiation of treatment with pembrolizumab.

MicroRNA-30c expression level is an independent predictor of clinical benefit of endocrine therapy in advanced estrogen receptor positive breast cancer.

MicroRNAs (miRNAs) are small RNA molecules that modulate gene expression and which have been implicated in cancer. We evaluated whether five candidate predictive miRNAs, derived from a pilot study in which 249 miRNAs were assayed, were associated with clinical benefit of tamoxifen therapy in advanced breast cancer. These five miRNAs were measured in an independent series of 246 estrogen receptor (ER)-positive primary breast tumors of patients who received tamoxifen for advanced disease by quantitative Real Time PCR. Univariate analysis showed that higher expression levels of hsa-miR-30a-3p, hsa-miR-30c, and hsa-miR-182 were significantly associated with benefit of tamoxifen treatment and with longer PFS (all P-values <0.01). In multivariate analysis, corrected for the traditional predictive factors, only hsa-miRNA-30c was an independent predictor (P-value <0.01). Finally, in an attempt to understand the biology connected to this miRNA, Global testing pathway analysis showed an association of hsa-miRNA-30c expression with HER and RAC1 signaling pathways. We identified hsa-miRNA-30c as an independent predictor for clinical benefit of tamoxifen therapy in patients with advanced breast cancer. Assessment of tumor levels and connected pathways could be helpful to improve treatment strategies.

Four miRNAs associated with aggressiveness of lymph node-negative, estrogen receptor-positive human breast cancer.

In this study, we quantified 249 mature micro-RNA (miRNA) transcripts in estrogen receptor-positive (ER(+)) primary breast tumors of patients with lymph node-negative (LNN) disease to identify miRNAs associated with metastatic capability. In addition, the prognostic value of the candidate miRNAs was determined in ER(-)/LNN breast cancer. Unsupervised analysis in a prescreening set of 38 patients identified three subgroups predominantly driven by three miRNA signatures: an ER-driven luminal B-associated miRNA signature, a stromal miRNA signature, and an overexpressed miRNA cluster located on chromosome 19q23, but these intrinsic miRNA signatures were not associated with tumor aggressiveness. Supervised analysis in the initial subset and subsequent analysis in additional tumors significantly linked four miRNAs (miR-7, miR-128a, miR-210, and miR-516-3p) to ER(+)/LNN breast cancer aggressiveness (n = 147) and one miRNA (miR-210) to metastatic capability in ER(-)/LNN breast cancer (n = 114) and in the clinically important triple-negative subgroup (n = 69) (all P < 0.05). Bioinformatic analysis coupled miR-210 to hypoxia/VEGF signaling, miR-7 and miR-516-3p to cell cycle progression and chromosomal instability, and miR-128a to cytokine signaling. In conclusion, our work connects four miRNAs to breast cancer progression and to several distinct biological processes involved therein.

A panel of DNA methylation markers for the classification of consensus molecular subtypes 2 and 3 in patients with colorectal cancer.

Consensus molecular subtypes (CMSs) can guide precision treatment of colorectal cancer (CRC). We aim to identify methylation markers to distinguish between CMS2 and CMS3 in patients with CRC, for which an easy test is currently lacking. To this aim, fresh-frozen tumor tissue of 239 patients with stage I-III CRC was analyzed. Methylation profiles were obtained using the Infinium HumanMethylation450 BeadChip. We performed adaptive group-regularized logistic ridge regression with post hoc group-weighted elastic net marker selection to build prediction models for classification of CMS2 and CMS3. The Cancer Genome Atlas (TCGA) data were used for validation. Group regularization of the probes was done based on their location either relative to a CpG island or relative to a gene present in the CMS classifier, resulting in two different prediction models and subsequently different marker panels. For both panels, even when using only five markers, accuracies were > 90% in our cohort and in the TCGA validation set. Our methylation marker panel accurately distinguishes between CMS2 and CMS3. This enables development of a targeted assay to provide a robust and clinically relevant classification tool for CRC patients.

Circular RNA in Chemonaive Lymph Node Negative Colon Cancer Patients.

Circular RNAs (circRNAs) appear important in tumor progression of colon cancer (CC). We identified an extensive catalog of circRNAs in 181 chemonaive stage I/II colon tumors, who underwent curative surgery between 2007 and 2014. We identified circRNAs from RNAseq data, investigated common biology related to circRNA expression, and studied the association between circRNAs and relapse status, tumor stage, consensus molecular subtypes (CMS), tumor localization and microsatellite instability (MSI). We identified 2606 unique circRNAs. 277 circRNAs (derived from 260 genes) were repeatedly occurring in at least 20 patients of which 153 showed a poor or even negative (R < 0.3) correlation with the expression level of their linear gene. The circular junctions for circSATB2, circFGD6, circKMT2C and circPLEKHM3 were validated by Sanger sequencing. Multiple correspondence analysis showed that circRNAs were often co-expressed and that high diversity in circRNAs was associated with favorable disease-free survival (DFS), which was confirmed by Cox regression analysis (Hazard Ratio (HR) 0.60, 95% CI 0.38-0.97, p = 0.036). Considering individual circRNAs, absence of circMGA was significantly associated with relapse, whereas circSATB2, circNAB1, and circCEP192 were associated with both MSI and CMS. This study represents a showcase of the potential clinical utility of circRNAs for prognostic stratification in patients with stage I-II colon cancer and demonstrated that high diversity in circRNAs is associated with favorable DFS.

Detection of Aneuploidy in Cerebrospinal Fluid from Patients with Breast Cancer Can Improve Diagnosis of Leptomeningeal Metastases.

PURPOSE: Detection of leptomeningeal metastasis is hampered by limited sensitivities of currently used techniques: MRI and cytology of cerebrospinal fluid (CSF). Detection of cell-free tumor DNA in CSF has been proposed as a tumor-specific candidate to detect leptomeningeal metastasis at an earlier stage. The aim of this study was to investigate mutation and aneuploidy status in CSF-derived cell-free DNA (cfDNA) of patients with breast cancer with a clinical suspicion of leptomeningeal metastasis. EXPERIMENTAL DESIGN: cfDNA was isolated from stored remnant CSF and analyzed by targeted next-generation sequencing (NGS; n = 30) and the modified fast aneuploidy screening test-sequencing system (mFAST-SeqS; n = 121). The latter method employs selective amplification of long interspaced nuclear elements sequences that are present throughout the genome and allow for fast and cheap detection of aneuploidy. We compared these results with the gold standard to diagnose leptomeningeal metastasis: cytology. RESULTS: Leptomeningeal metastasis was cytology proven in 13 of 121 patients. Low DNA yields resulted in insufficient molecular coverage of NGS for the majority of samples (success rate, 8/30). The mFAST-SeqS method, successful in 112 of 121 (93%) samples, detected genome-wide aneuploidy in 24 patients. Ten of these patients had cytology-proven leptomeningeal metastasis; 8 additional patients were either concurrently diagnosed with central nervous system metastases by radiological means or developed these soon after the lumbar puncture. The remaining six cases were suspected of leptomeningeal metastasis, but could not be confirmed by cytology or imaging. Aneuploidy was associated with development of leptomeningeal metastasis and significantly worse overall survival. CONCLUSIONS: Aneuploidy in CSF-derived cfDNA may provide a promising biomarker to improve timely detection of leptomeningeal metastasis.

High-throughput and affordable genome-wide methylation profiling of circulating cell-free DNA by methylated DNA sequencing (MeD-seq) of LpnPI digested fragments.

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.