Beyond basics: Key mutation selection features for successful tumor-informed ctDNA detection.

Tumor-informed mutation-based approaches are frequently used for detection of circulating tumor DNA (ctDNA). Not all mutations make equally effective ctDNA markers. The objective was to explore if prioritizing mutations using mutational features-such as cancer cell fraction (CCF), multiplicity, and error rate-would improve the success rate of tumor-informed ctDNA analysis. Additionally, we aimed to develop a practical and easily implementable analysis pipeline for identifying and prioritizing candidate mutations from whole-exome sequencing (WES) data. We analyzed WES and ctDNA data from three tumor-informed ctDNA studies, one on bladder cancer (Cohort A) and two on colorectal cancer (Cohorts I and N). The studies included 390 patients. For each patient, a unique set of mutations (median mutations/patient: 6, interquartile 13, range: 1-46, total n = 4023) were used as markers of ctDNA. The tool PureCN was used to assess the CCF and multiplicity of each mutation. High-CCF mutations were detected more frequently than low-CCF mutations (Cohort A: odds ratio [OR] 20.6, 95% confidence interval [CI] 5.72-173, p = 1.73e-12; Cohort I: OR 2.24, 95% CI 1.44-3.52, p = 1.66e-04; and Cohort N: OR 1.78, 95% CI 1.14-2.79, p = 7.86e-03). The detection-likelihood was additionally improved by selecting mutations with multiplicity of two or above (Cohort A: OR 1.55, 95% CI 1. 14-2.11, p = 3.85e-03; Cohort I: OR 1.78, 95% CI 1.23-2.56, p = 1.34e-03; and Cohort N: OR 1.94, 95% CI 1.63-2.31, p = 2.83e-14). Furthermore, selecting the mutations for which the ctDNA detection method had the lowest error rates, additionally improved the detection-likelihood, particularly evident when plasma cell-free DNA tumor fractions were below 0.1% (p = 2.1e-07). Selecting mutational markers with high CCF, high multiplicity, and low error rate significantly improve ctDNA detection likelihood. We provide free access to the analysis pipeline enabling others to perform qualified prioritization of mutations for tumor-informed ctDNA analysis.

DREAMS: deep read-level error model for sequencing data applied to low-frequency variant calling and circulating tumor DNA detection.

Circulating tumor DNA detection using next-generation sequencing (NGS) data of plasma DNA is promising for cancer identification and characterization. However, the tumor signal in the blood is often low and difficult to distinguish from errors. We present DREAMS (Deep Read-level Modelling of Sequencing-errors) for estimating error rates of individual read positions. Using DREAMS, we develop statistical methods for variant calling (DREAMS-vc) and cancer detection (DREAMS-cc). For evaluation, we generate deep targeted NGS data of matching tumor and plasma DNA from 85 colorectal cancer patients. The DREAMS approach performs better than state-of-the-art methods for variant calling and cancer detection.

Error Characterization and Statistical Modeling Improves Circulating Tumor DNA Detection by Droplet Digital PCR.

BACKGROUND: Droplet digital PCR (ddPCR) is a widely used and sensitive application for circulating tumor DNA (ctDNA) detection. As ctDNA is often found in low abundance, methods to separate low-signal readouts from noise are necessary. We aimed to characterize the ddPCR-generated noise and, informed by this, create a sensitive and specific ctDNA caller. METHODS: We built 2 novel complimentary ctDNA calling methods: dynamic limit of blank and concentration and assay-specific tumor load estimator (CASTLE). Both methods are informed by empirically established assay-specific noise profiles. Here, we characterized noise for 70 mutation-detecting ddPCR assays by applying each assay to 95 nonmutated samples. Using these profiles, the performance of the 2 new methods was assessed in a total of 9447 negative/positive reference samples and in 1311 real-life plasma samples from colorectal cancer patients. Lastly, performances were compared to 7 literature-established calling methods. RESULTS: For many assays, noise increased proportionally with the DNA input amount. Assays targeting transition base changes were more error-prone than transversion-targeting assays. Both our calling methods successfully accounted for the additional noise in transition assays and showed consistently high performance regardless of DNA input amount. Calling methods that were not noise-informed performed less well than noise-informed methods. CASTLE was the only calling method providing a statistical estimate of the noise-corrected mutation level and call certainty. CONCLUSIONS: Accurate error modeling is necessary for sensitive and specific ctDNA detection by ddPCR. Accounting for DNA input amounts ensures specific detection regardless of the sample-specific DNA concentration. Our results demonstrate CASTLE as a powerful tool for ctDNA calling using ddPCR.

Publisher Correction: MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

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

MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

Transcriptional characterization and classification has potential to resolve the inter-tumor heterogeneity of colorectal cancer and improve patient management. Yet, robust transcriptional profiling is difficult using formalin-fixed, paraffin-embedded (FFPE) samples, which complicates testing in clinical and archival material. We present MethCORR, an approach that allows uniform molecular characterization and classification of fresh-frozen and FFPE samples. MethCORR identifies genome-wide correlations between RNA expression and DNA methylation in fresh-frozen samples. This information is used to infer gene expression information in FFPE samples from their methylation profiles. MethCORR is here applied to methylation profiles from 877 fresh-frozen/FFPE samples and comparative analysis identifies the same two subtypes in four independent cohorts. Furthermore, subtype-specific prognostic biomarkers that better predicts relapse-free survival (HR = 2.66, 95%CI [1.67-4.22], P value < 0.001 (log-rank test)) than UICC tumor, node, metastasis (TNM) staging and microsatellite instability status are identified and validated using DNA methylation-specific PCR. The MethCORR approach is general, and may be similarly successful for other cancer types.

A genetically inducible porcine model of intestinal cancer.

Transgenic porcine cancer models bring novel possibilities for research. Their physical similarities with humans enable the use of surgical procedures and treatment approaches used for patients, which facilitates clinical translation. Here, we aimed to develop an inducible oncopig model of intestinal cancer. Transgenic (TG) minipigs were generated using somatic cell nuclear transfer by handmade cloning. The pigs encode two TG cassettes: (a) an Flp recombinase-inducible oncogene cassette containing KRAS-G12D, cMYC, SV40LT - which inhibits p53 - and pRB and (b) a 4-hydroxytamoxifen (4-OHT)-inducible Flp recombinase activator cassette controlled by the intestinal epithelium-specific villin promoter. Thirteen viable transgenic minipigs were born. The ability of 4-OHT to activate the oncogene cassette was confirmed in vitro in TG colonic organoids and ex vivo in tissue biopsies obtained by colonoscopy. In order to provide proof of principle that the oncogene cassette could also successfully be activated in vivo, three pigs were perorally treated with 400 mg tamoxifen for 2 × 5 days. After two months, one pig developed a duodenal neuroendocrine carcinoma with a lymph node metastasis. Molecular analysis of the carcinoma and metastasis confirmed activation of the oncogene cassette. No tumor formation was observed in untreated TG pigs or in the remaining two treated pigs. The latter indicates that tamoxifen delivery can probably be improved. In summary, we have generated a novel inducible oncopig model of intestinal cancer, which has the ability to form metastatic disease already two months after induction. The model may be helpful in bridging the gap between basic research and clinical usage. It opens new venues for longitudinal studies of tumor development and evolution, for preclinical assessment of new anticancer regimens, for pharmacology and toxicology assessments, as well as for studies into biological mechanisms of tumor formation and metastasis.

Establishment and characterization of models of chemotherapy resistance in colorectal cancer: Towards a predictive signature of chemoresistance.

Current standard treatments for metastatic colorectal cancer (CRC) are based on combination regimens with one of the two chemotherapeutic drugs, irinotecan or oxaliplatin. However, drug resistance frequently limits the clinical efficacy of these therapies. In order to gain new insights into mechanisms associated with chemoresistance, and departing from three distinct CRC cell models, we generated a panel of human colorectal cancer cell lines with acquired resistance to either oxaliplatin or irinotecan. We characterized the resistant cell line variants with regards to their drug resistance profile and transcriptome, and matched our results with datasets generated from relevant clinical material to derive putative resistance biomarkers. We found that the chemoresistant cell line variants had distinctive irinotecan- or oxaliplatin-specific resistance profiles, with non-reciprocal cross-resistance. Furthermore, we could identify several new, as well as some previously described, drug resistance-associated genes for each resistant cell line variant. Each chemoresistant cell line variant acquired a unique set of changes that may represent distinct functional subtypes of chemotherapy resistance. In addition, and given the potential implications for selection of subsequent treatment, we also performed an exploratory analysis, in relevant patient cohorts, of the predictive value of each of the specific genes identified in our cellular models.

Cellular disposal of miR23b by RAB27-dependent exosome release is linked to acquisition of metastatic properties.

Exosomes are small secreted vesicles that can transfer their content to recipient cells. In cancer, exosome secretion has been implicated in tumor growth and metastatic spread. In this study, we explored the possibility that exosomal pathways might discard tumor-suppressor miRNA that restricts metastatic progression. Secreted miRNA characterized from isogenic bladder carcinoma cell lines with differing metastatic potential were uncoupled from binding to target transcripts or the AGO2-miRISC complex. In metastatic cells, we observed a relative increase in secretion of miRNA with tumor-suppressor functions, including miR23b, miR224, and miR921. Ectopic expression of miR23b inhibited invasion, anoikis, angiogenesis, and pulmonary metastasis. Silencing of the exocytotic RAB family members RAB27A or RAB27B halted miR23b and miR921 secretion and reduced cellular invasion. Clinically, elevated levels of RAB27B expression were linked to poor prognosis in two independent cohorts of patients with bladder cancer. Moreover, highly exocytosed miRNA from metastatic cells, such as miR23b, were reduced in lymph node metastases compared with patient-matched primary tumors and were correlated with increments in miRNA-targeted RNA. Taken together, our results suggested that exosome-mediated secretion of tumor-suppressor miRNA is selected during tumor progression as a mechanism to coordinate activation of a metastatic cascade.

Putative cis-regulatory drivers in colorectal cancer.

The cis-regulatory effects responsible for cancer development have not been as extensively studied as the perturbations of the protein coding genome in tumorigenesis. To better characterize colorectal cancer (CRC) development we conducted an RNA-sequencing experiment of 103 matched tumour and normal colon mucosa samples from Danish CRC patients, 90 of which were germline-genotyped. By investigating allele-specific expression (ASE) we show that the germline genotypes remain important determinants of allelic gene expression in tumours. Using the changes in ASE in matched pairs of samples we discover 71 genes with excess of somatic cis-regulatory effects in CRC, suggesting a cancer driver role. We correlate genotypes and gene expression to identify expression quantitative trait loci (eQTLs) and find 1,693 and 948 eQTLs in normal samples and tumours, respectively. We estimate that 36% of the tumour eQTLs are exclusive to CRC and show that this specificity is partially driven by increased expression of specific transcription factors and changes in methylation patterns. We show that tumour-specific eQTLs are more enriched for low CRC genome-wide association study (GWAS) P values than shared eQTLs, which suggests that some of the GWAS variants are tumour specific regulatory variants. Importantly, tumour-specific eQTL genes also accumulate more somatic mutations when compared to the shared eQTL genes, raising the possibility that they constitute germline-derived cancer regulatory drivers. Collectively the integration of genome and the transcriptome reveals a substantial number of putative somatic and germline cis-regulatory cancer changes that may have a role in tumorigenesis.

Functional screening identifies miRNAs influencing apoptosis and proliferation in colorectal cancer.

MicroRNAs (miRNAs) play a critical role in many biological processes and are aberrantly expressed in human cancers. Particular miRNAs function either as tumor suppressors or oncogenes and appear to have diagnostic and prognostic significance. Although numerous miRNAs are dys-regulated in colorectal cancer (CRC) only a small fraction has been characterized functionally. Using high-throughput functional screening and miRNA profiling of clinical samples the present study aims at identifying miRNAs important for the control of cellular growth and/or apoptosis in CRC. The high-throughput functional screening was carried out in six CRC cell lines transfected with a pre-miR library including 319 synthetic human pre-miRs. Phenotypic alterations were evaluated by immunostaining of cleaved cPARP (apoptosis) or MKI67 (proliferation). Additionally, TaqMan Human MicroRNA Array Set v2.0 was used to profile the expression of 667 miRNAs in 14 normal colon mucosa and 46 microsatellite stable stage II CRC patients. Among the miRNAs that induced growth arrest and apoptosis in the CRC cell lines, and at same time were dys-regulated in the clinical samples, miR-375 was selected for further analysis. Independent in vitro analysis of transient and stable transfected CRC cell lines confirmed that miR-375 reduces cell viability through the induction of apoptotic death. We identified YAP1 as a direct miR-375 target in CRC and show that HELLS and NOLC1 are down-stream targets. Knock-down of YAP1 mimicked the phenotype induced by miR-375 over-expression indicating that miR-375 most likely exerts its pro-apoptotic role through YAP1 and its anti-apoptotic down-stream targets BIRC5 and BCL2L1. Finally, in vivo analysis of mouse xenograft tumors showed that miR-375 expression significantly reduced tumor growth. We conclude that the high-throughput screening successfully identified miRNAs that induce apoptosis and/or inhibit proliferation in CRC cells. Finally, combining the functional screening with profiling of CRC tissue samples we identified clinically relevant miRNAs and miRNA targets in CRC.

High expression of microRNA-625-3p is associated with poor response to first-line oxaliplatin based treatment of metastatic colorectal cancer.

The backbone of current cytotoxic treatment of metastatic colorectal cancer (mCRC) consists of a fluoropyrimidine together with either oxaliplatin (XELOX/FOLFOX) or irinotecan (XELIRI/FOLFIRI). With an overall objective response rate of approximately 50% for either treatment combination, a major unsolved problem is that no predictors of response to these treatments are available. To address this issue, we profiled 742 microRNAs in laser-capture microdissected cancer cells from responding and non-responding patients receiving XELOX/FOLFOX as first-line treatment for mCRC, and identified, among others, high expression of miR-625-3p, miR-181b and miR-27b to be associated with poor clinical response. In a validation cohort of 94 mCRC patients treated first-line with XELOX, high expression of miR-625-3p was confirmed to be associated with poor response (OR = 6.25, 95%CI [1.8; 21.0]). Independent analyses showed that miR-625-3p was not dysregulated between normal and cancer samples, nor was its expression associated with recurrence of stage II or III disease, indicating that miR-625-3p solely is a response marker. Finally, we also found that these miRNAs were up-regulated in oxaliplatin resistant HCT116/oxPt (miR-625-3p, miR-181b and miR-27b) and LoVo/oxPt (miR-181b) colon cancer cell lines as compared with their isogenic parental cells. Altogether, our results suggest an association between miR-625-3p and response to first-line oxaliplatin based chemotherapy of mCRC.

MiRNA-362-3p induces cell cycle arrest through targeting of E2F1, USF2 and PTPN1 and is associated with recurrence of colorectal cancer.

Colorectal cancer (CRC) is one of the leading causes of cancer deaths in Western countries. A significant number of CRC patients undergoing curatively intended surgery subsequently develop recurrence and die from the disease. MicroRNAs (miRNAs) are aberrantly expressed in cancers and appear to have both diagnostic and prognostic significance. In this study, we identified novel miRNAs associated with recurrence of CRC, and their possible mechanism of action. TaqMan(®) Human MicroRNA Array Set v2.0 was used to profile the expression of 667 miRNAs in 14 normal colon mucosas and 46 microsatellite stable CRC tumors. Four miRNAs (miR-362-3p, miR-570, miR-148 a* and miR-944) were expressed at a higher level in tumors from patients with no recurrence (p<0.015), compared with tumors from patients with recurrence. A significant association with increased disease free survival was confirmed for miR-362-3p in a second independent cohort of 43 CRC patients, using single TaqMan(®) microRNA assays. In vitro functional analysis showed that over-expression of miR-362-3p in colon cancer cell lines reduced cell viability, and proliferation mainly due to cell cycle arrest. E2F1, USF2 and PTPN1 were identified as potential miR-362-3p targets by mRNA profiling of HCT116 cells over-expressing miR-362-3p. Subsequently, these genes were confirmed as direct targets by Luciferase reporter assays and their knockdown in vitro phenocopied the effects of miR-362-3p over-expression. We conclude that miR-362-3p may be a novel prognostic marker in CRC, and hypothesize that the positive effects of augmented miR-362-3p expression may in part be mediated through the targets E2F1, USF2 and PTPN1.

Tumor-specific usage of alternative transcription start sites in colorectal cancer identified by genome-wide exon array analysis.

BACKGROUND: Approximately half of all human genes use alternative transcription start sites (TSSs) to control mRNA levels and broaden the transcriptional output in healthy tissues. Aberrant expression patterns promoting carcinogenesis, however, may arise from alternative promoter usage. RESULTS: By profiling 108 colorectal samples using exon arrays, we identified nine genes (TCF12, OSBPL1A, TRAK1, ANK3, CHEK1, UGP2, LMO7, ACSL5, and SCIN) showing tumor-specific alternative TSS usage in both adenoma and cancer samples relative to normal mucosa. Analysis of independent exon array data sets corroborated these findings. Additionally, we confirmed the observed patterns for selected mRNAs using quantitative real-time reverse-transcription PCR. Interestingly, for some of the genes, the tumor-specific TSS usage was not restricted to colorectal cancer. A comprehensive survey of the nine genes in lung, bladder, liver, prostate, gastric, and brain cancer revealed significantly altered mRNA isoform ratios for CHEK1, OSBPL1A, and TCF12 in a subset of these cancer types.To identify the mechanism responsible for the shift in alternative TSS usage, we antagonized the Wnt-signaling pathway in DLD1 and Ls174T colorectal cancer cell lines, which remarkably led to a shift in the preferred TSS for both OSBPL1A and TRAK1. This indicated a regulatory role of the Wnt pathway in selecting TSS, possibly also involving TP53 and SOX9, as their transcription binding sites were enriched in the promoters of the tumor preferred isoforms together with their mRNA levels being increased in tumor samples. Finally, to evaluate the prognostic impact of the altered TSS usage, immunohistochemistry was used to show deregulation of the total protein levels of both TCF12 and OSBPL1A, corresponding to the mRNA levels observed. Furthermore, the level of nuclear TCF12 had a significant correlation to progression free survival in a cohort of 248 stage II colorectal cancer samples. CONCLUSIONS: Alternative TSS usage in colorectal adenoma and cancer samples has been shown for nine genes, and OSBPL1A and TRAK1 were found to be regulated in vitro by Wnt signaling. TCF12 protein expression was upregulated in cancer samples and correlated with progression free survival.

Evaluation of two commercial global miRNA expression profiling platforms for detection of less abundant miRNAs.

BACKGROUND: microRNAs (miRNA) are short, endogenous transcripts that negatively regulate the expression of specific mRNA targets. miRNAs are found both in tissues and body fluids such as plasma. A major perspective for the use of miRNAs in the clinical setting is as diagnostic plasma markers for neoplasia. While miRNAs are abundant in tissues, they are often scarce in plasma. For quantification of miRNA in plasma it is therefore of importance to use a platform with high sensitivity and linear performance in the low concentration range. This motivated us to evaluate the performance of three commonly used commercial miRNA quantification platforms: GeneChip miRNA 2.0 Array, miRCURY Ready-to-Use PCR, Human panel I+II V1.M, and TaqMan Human MicroRNA Array v3.0. RESULTS: Using synthetic miRNA samples and plasma RNA samples spiked with different ratios of 174 synthetic miRNAs we assessed the performance characteristics reproducibility, recovery, specificity, sensitivity and linearity. It was found that while the qRT-PCR based platforms were sufficiently sensitive to reproducibly detect miRNAs at the abundance levels found in human plasma, the array based platform was not. At high miRNA levels both qRT-PCR based platforms performed well in terms of specificity, reproducibility and recovery. At low miRNA levels, as in plasma, the miRCURY platform showed better sensitivity and linearity than the TaqMan platform. CONCLUSION: For profiling clinical samples with low miRNA abundance, such as plasma samples, the miRCURY platform with its better sensitivity and linearity would probably be superior.

Alternative splicing of SLC39A14 in colorectal cancer is regulated by the Wnt pathway.

Alternative splicing is a crucial step in the generation of protein diversity and its misregulation is observed in many human cancer types. By analyzing 143 colorectal samples using exon arrays, SLC39A14, a divalent cation transporter, was identified as being aberrantly spliced in tumor samples. SLC39A14 contains two mutually exclusive exons 4A and 4B and the exon 4A/4B ratio was significantly altered in adenomas (p = 3.6 × 10(-10)) and cancers (p = 9.4 × 10(-11)), independent of microsatellite stability status. The findings were validated in independent exon array data sets and by quantitative real-time reverse-transcription PCR (qRT-PCR). Aberrant Wnt signaling is a hallmark of colorectal tumorigenesis and is characterized by nuclear ß-catenin. Experimental inactivation of Wnt signaling in DLD1 and Ls174T cells by knockdown of ß-catenin or overexpression of dominant negative TCFs (TCF1 and TCF4) altered the 4A/4B ratio, indicating that SLC39A14 splicing is regulated by the Wnt pathway. An altered 4A/4B ratio was also observed in gastric and lung cancer where Wnt signaling is also known to be aberrantly activated. The splicing factor SRSF1 and its regulator, the kinase SRPK1, were found to be deregulated upon Wnt inactivation in colorectal carcinoma cells. SRPK1 was also found up-regulated in both adenoma samples (p = 1.5 × 10(-5)) and cancer samples (p = 5 × 10(-4)). In silico splicing factor binding analysis predicted SRSF1 to bind predominantly to the cancer associated exon 4B, hence, it was hypothesized that SRPK1 activates SRSF1 through phosphorylation, followed by SRSF1 binding to exon 4B and regulation of SLC39A14 splicing. Indeed, siRNA-mediated knockdown of SRPK1 and SRSF1 in DLD1 and SW480 colorectal cancer cells led to a change in the 4A/4B isoform ratio, supporting a role of these factors in the regulation of SLC39A14 splicing. In conclusion, alternative splicing of SLC39A14 was identified in colorectal tumors and found to be regulated by the Wnt pathway, most likely through regulation of SRPK1 and SRSF1.