Identifying Recurrences Among Non-Metastatic Colorectal Cancer Patients Using National Health Data Registries: Validation and Optimization of a Registry-Based Algorithm in a Modern Danish Cohort.

Purpose: Colorectal cancer (CRC) recurrence is not routinely recorded in Danish health data registries. Here, we aimed to revalidate a registry-based algorithm to identify recurrences in a contemporary cohort and to investigate the accuracy of estimating the time to recurrence (TTR). Patients and Methods: We ascertained data on 1129 patients operated for UICC TNM stage I-III CRC during 2012-2017 registered in the CRC biobank at the Department of Molecular Medicine, Aarhus University Hospital, Denmark. Individual-level data were linked with data from the Danish Colorectal Cancer Group database, Danish Cancer Registry, Danish National Registry of Patients, and Danish Pathology Registry. The algorithm identified recurrence based on diagnosis codes of local recurrence or metastases, the receipt of chemotherapy, or a pathological tissue assessment code of recurrence more than 180 days after CRC surgery. A subgroup was selected for validation of the algorithm using medical record reviews as a reference standard. Results: We found a 3-year cumulative recurrence rate of 20% (95% CI: 17-22%). Manual medical record review identified 80 recurrences in the validation cohort of 522 patients. The algorithm detected recurrence with 94% sensitivity (75/80; 95% CI: 86-98%) and 98% specificity (431/442; 95% CI: 96-99%). The positive and negative predictive values of the algorithm were 87% (95% CI: 78-93%) and 99% (95% CI: 97-100%), respectively. The median difference in TTR (TTRMedical_chart-TTRalgorithm) was -8 days (IQR: -21 to +3 days). Restricting the algorithm to chemotherapy codes from oncology departments increased the positive predictive value from 87% to 94% without changing the negative predictive value (99%). Conclusion: The algorithm detected recurrence and TTR with high precision in this contemporary cohort. Restriction to chemotherapy codes from oncology departments using department classifications improves the algorithm. The algorithm is suitable for use in future observational studies.

MSIMEP: Predicting microsatellite instability from microarray DNA methylation tumor profiles.

Deficiency in DNA MMR activity results in tumors with a hypermutator phenotype, termed microsatellite instability (MSI). Beyond its utility in Lynch syndrome screening algorithms, today MSI has gained importance as predictive biomarker for various anti-PD-1 therapies across many different tumor types. Over the past years, many computational methods have emerged to infer MSI using either DNA- or RNA-based approaches. Considering this together with the fact that MSI-high tumors frequently exhibit a hypermethylated phenotype, herein we developed and validated MSIMEP, a computational tool for predicting MSI status from microarray DNA methylation tumor profiles of colorectal cancer samples. We demonstrated that MSIMEP optimized and reduced models have high performance in predicting MSI in different colorectal cancer cohorts. Moreover, we tested its consistency in other tumor types with high prevalence of MSI such as gastric and endometrial cancers. Finally, we demonstrated better performance of both MSIMEP models vis-à-vis a MLH1 promoter methylation-based one in colorectal cancer.

Transposon-activated POU5F1B promotes colorectal cancer growth and metastasis.

The treatment of colorectal cancer (CRC) is an unmet medical need in absence of early diagnosis. Here, upon characterizing cancer-specific transposable element-driven transpochimeric gene transcripts (TcGTs) produced by this tumor in the SYSCOL cohort, we find that expression of the hominid-restricted retrogene POU5F1B through aberrant activation of a primate-specific endogenous retroviral promoter is a strong negative prognostic biomarker. Correlating this observation, we demonstrate that POU5F1B fosters the proliferation and metastatic potential of CRC cells. We further determine that POU5F1B, in spite of its phylogenetic relationship with the POU5F1/OCT4 transcription factor, is a membrane-enriched protein that associates with protein kinases and known targets or interactors as well as with cytoskeleton-related molecules, and induces intracellular signaling events and the release of trans-acting factors involved in cell growth and cell adhesion. As POU5F1B is an apparently non-essential gene only lowly expressed in normal tissues, and as POU5F1B-containing TcGTs are detected in other tumors besides CRC, our data provide interesting leads for the development of cancer therapies.

Transcriptomic and proteomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum.

BACKGROUND: Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH. METHODS: Multiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping was performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors). RESULTS: Overall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high transcriptional ITH and were classified discordantly. Subtyping of proximal MSS tumors were discordant for 50% of the tumors, this ITH was related to differences in the microenvironment. Subtyping of distal MSS tumors were less discordant, here the ITH was more cancer-cell related. The subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05). Unsupervised hierarchical clustering of the protein data identified large ITH at protein level; as the multiregional biopsies clustered together for only 9 out of 29 tumors. CONCLUSION: Our transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.

Novel DNA methylation biomarkers show high sensitivity and specificity for blood-based detection of colorectal cancer-a clinical biomarker discovery and validation study.

BACKGROUND: Early detection plays an essential role to reduce colorectal cancer (CRC) mortality. While current screening methods suffer from poor compliance, liquid biopsy-based strategies for cancer detection is rapidly gaining promise. Here, we describe the development of TriMeth, a minimal-invasive blood-based test for detection of early-stage colorectal cancer. The test is based on assessment of three tumour-specific DNA methylation markers in circulating cell-free DNA. RESULTS: A thorough multi-step biomarker discovery study based on DNA methylation profiles of more than 5000 tumours and blood cell populations identified CRC-specific DNA methylation markers. The DNA methylation patterns of biomarker candidates were validated by bisulfite sequencing and methylation-specific droplet digital PCR in CRC tumour tissue and peripheral blood leucocytes. The three best performing markers were first applied to plasma from 113 primarily early-stage CRC patients and 87 age- and gender-matched colonoscopy-verified controls. Based on this, the test scoring algorithm was locked, and then TriMeth was validated in an independent cohort comprising 143 CRC patients and 91 controls. Three DNA methylation markers, C9orf50, KCNQ5, and CLIP4, were identified, each capable of discriminating plasma from colorectal cancer patients and healthy individuals (areas under the curve 0.86, 0.91, and 0.88). When combined in the TriMeth test, an average sensitivity of 85% (218/256) was observed (stage I: 80% (33/41), stage II: 85% (121/143), stage III: 89% (49/55), and stage IV: 88% (15/17)) at 99% (176/178) specificity in two independent plasma cohorts. CONCLUSION: TriMeth enables detection of early-stage colorectal cancer with high sensitivity and specificity. The reported results underline the potential utility of DNA methylation-based detection of circulating tumour DNA in the clinical management of colorectal cancer.

Contribution of allelic imbalance to colorectal cancer.

Point mutations in cancer have been extensively studied but chromosomal gains and losses have been more challenging to interpret due to their unspecific nature. Here we examine high-resolution allelic imbalance (AI) landscape in 1699 colorectal cancers, 256 of which have been whole-genome sequenced (WGSed). The imbalances pinpoint 38 genes as plausible AI targets based on previous knowledge. Unbiased CRISPR-Cas9 knockout and activation screens identified in total 79 genes within AI peaks regulating cell growth. Genetic and functional data implicate loss of TP53 as a sufficient driver of AI. The WGS highlights an influence of copy number aberrations on the rate of detected somatic point mutations. Importantly, the data reveal several associations between AI target genes, suggesting a role for a network of lineage-determining transcription factors in colorectal tumorigenesis. Overall, the results unravel the contribution of AI in colorectal cancer and provide a plausible explanation why so few genes are commonly affected by point mutations in cancers.

Molecular-Subtype-Specific Biomarkers Improve Prediction of Prognosis in Colorectal Cancer.

Colorectal cancer (CRC) is characterized by major inter-tumor diversity that complicates the prediction of disease and treatment outcomes. Recent efforts help resolve this by sub-classification of CRC into natural molecular subtypes; however, this strategy is not yet able to provide clinicians with improved tools for decision making. We here present an extended framework for CRC stratification that specifically aims to improve patient prognostication. Using transcriptional profiles from 1,100 CRCs, including >300 previously unpublished samples, we identify cancer cell and tumor archetypes and suggest the tumor microenvironment as a major prognostic determinant that can be influenced by the microbiome. Notably, our subtyping strategy allowed identification of archetype-specific prognostic biomarkers that provided information beyond and independent of UICC-TNM staging, MSI status, and consensus molecular subtyping. The results illustrate that our extended subtyping framework, combining subtyping and subtype-specific biomarkers, could contribute to improved patient prognostication and may form a strong basis for future studies.

Comparative analysis of 12 different kits for bisulfite conversion of circulating cell-free DNA.

Blood circulating cell-free DNA (cfDNA) is becoming popular in the search of promising predictive and prognostic biomarkers. Among these biomarkers, cfDNA methylation markers have especially gained considerable attention. A significant challenge in the utilization of cfDNA methylation markers is the limited amount of cfDNA available for analyses; reportedly, bisulfite conversion (BSC) reduce cfDNA amounts even further. Nevertheless, few efforts have focused on ensuring high cfDNA conversion efficiency and recovery after BSC. To compare cfDNA recovery of different BSC methods, we compared 12 different commercially available BSC kits. We tested whether DNA recovery was affected by the molecular weight and/or quantity of input DNA. We also tested BSC efficiency for each kit. We found that recovery varied for DNA fragments of different lengths: certain kits recovered short fragments better than others, and only 3 kits recovered DNA fragments of <100 bp well. In contrast, DNA input amount did not seem to affect DNA recovery: for quantities spanning between 820 and ∼25,000 genome equivalents per BSC, a linear relation was found between input and recovery amount. Overall, mean recovery ranged between 9 and 32%, with BSC efficiency of 97-99.9%. When plasma cfDNA was used as input for BSC, recovery varied from 22% for the poorest and 66% for the best performing kits, while conversion efficiency ranged from 96 to 100% among different kits. In conclusion, clear performance differences exist between commercially available BSC kits, both in terms of DNA recovery and conversion efficiency. The choice of BSC kit can substantially impact the amount of converted cfDNA available for downstream analysis, which is critical in a cfDNA methylation marker setting.

The proto-oncogene TWIST1 is regulated by microRNAs.

Upregulation of the proto-oncogene Twist1 is highly correlated with acquired drug resistance and poor prognosis in human cancers. Altered expression of this multifunctional transcription factor is also associated with inherited skeletal malformations. The mammalian Twist1 3'UTRs are highly conserved and contain a number of potential regulatory elements including miRNA target sites. We analyzed the translational regulation of TWIST1 using luciferase reporter assays in a variety of cell lines. Among several miRNAs tested, miR-145a-5p, miR-151-5p and a combination of miR-145a-5p + miR-151-5p and miR-151-5p + miR-337-3p were able to significantly repress Twist1 translation. This phenomena was confirmed with both exogenous and endogenous miRNAs and was dependent on the presence of the predicted target sites in the 3'UTR. Furthermore, the repression was sensitive to LNA-modified miRNA antagonists and resulted in decreased migratory potential of murine embryonic fibroblast cells. Understanding the in vivo mechanisms of this oncogene's regulation might open up a possibility for therapeutic interference by gene specific cancer therapies.

Functional selection of shRNA loops from randomized retroviral libraries.

Gene silencing by RNA interference (RNAi) can be achieved by the ectopic expression of tailored short hairpin RNAs (shRNAs) which after export to the cytoplasm are processed by Dicer and incorporated into the RNA induced silencing complex (RISC). Design rules for shRNAs have been the focus of several studies, but only a few reports have turned the attention to the sequence of the loop-region. In this work we selected high-functional and low-functional shRNA loops from retroviral hairpin-loop-libraries in an RNAi reporter assay. The procedure revealed a very significant and stem sequence-dependent effect of the loop on shRNA function and although neither strong consensus loop sequence nor structural motifs could be identified, a preferred loop sequence (5'-UGUGCUU-3') was found to support robust knock down with little stem sequence dependency. These findings will serve as a guide for designing shRNAs with improved knock down capacity.

MicroRNA expression profiling of carcinoma in situ cells of the testis.

Testicular germ cell tumours, seminoma (SE) and non-seminoma (NS), of young adult men develop from a precursor cell, carcinoma in situ (CIS), which resembles foetal gonocytes and retains embryonic pluripotency. We used microarrays to analyse microRNA (miRNA) expression in 12 human testis samples with CIS cells and compared it with miRNA expression profiles of normal adult testis, testis with Sertoli-cell-only that lacks germ cells, testis tumours (SE and embryonal carcinoma (EC), an undifferentiated component of NS) and foetal male and female gonads. Principal components analysis revealed distinct miRNA expression profiles characteristic for each of the different tissue types. We identified several miRNAs that were unique to testis with CIS cells, foetal gonads and testis tumours. These included miRNAs from the hsa-miR-371-373 and -302-367 clusters that have previously been reported in germ cell tumours and three miRNAs (hsa-miR-96, -141 and -200c) that were also expressed in human epididymis. We found several miRNAs that were upregulated in testis tumours: hsa-miR-9, -105 and -182-183-96 clusters were highly expressed in SE, while the hsa-miR-515-526 cluster was high in EC. We conclude that the miRNA expression profile changes during testis development and that the miRNA profile of adult testis with CIS cells shares characteristic similarities with the expression in foetal gonocytes.

In vivo screening of modified siRNAs for non-specific antiviral effect in a small fish model: number and localization in the strands are important.

Small interfering RNAs (siRNAs) are promising new active compounds in gene medicine but the induction of non-specific immune responses following their delivery continues to be a serious problem. With the purpose of avoiding such effects chemically modified siRNAs are tested in screening assay but often only examining the expression of specific immunologically relevant genes in selected cell populations typically blood cells from treated animals or humans. Assays using a relevant physiological state in biological models as read-out are not common. Here we use a fish model where the innate antiviral effect of siRNAs is functionally monitored as reduced mortality in challenge studies involving an interferon sensitive virus. Modifications with locked nucleic acid (LNA), altritol nucleic acid (ANA) and hexitol nucleic acid (HNA) reduced the antiviral protection in this model indicative of altered immunogenicity. For LNA modified siRNAs, the number and localization of modifications in the single strands was found to be important and a correlation between antiviral protection and the thermal stability of siRNAs was found. The previously published sisiRNA will in some sequences, but not all, increase the antiviral effect of siRNAs. The applied fish model represents a potent tool for conducting fast but statistically and scientifically relevant evaluations of chemically optimized siRNAs with respect to non-specific antiviral effects in vivo.

Re-inspection of small RNA sequence datasets reveals several novel human miRNA genes.

BACKGROUND: miRNAs are key players in gene expression regulation. To fully understand the complex nature of cellular differentiation or initiation and progression of disease, it is important to assess the expression patterns of as many miRNAs as possible. Thereby, identifying novel miRNAs is an essential prerequisite to make possible a comprehensive and coherent understanding of cellular biology. METHODOLOGY/PRINCIPAL FINDINGS: Based on two extensive, but previously published, small RNA sequence datasets from human embryonic stem cells and human embroid bodies, respectively [1], we identified 112 novel miRNA-like structures and were able to validate miRNA processing in 12 out of 17 investigated cases. Several miRNA candidates were furthermore substantiated by including additional available small RNA datasets, thereby demonstrating the power of combining datasets to identify miRNAs that otherwise may be assigned as experimental noise. CONCLUSIONS/SIGNIFICANCE: Our analysis highlights that existing datasets are not yet exhaustedly studied and continuous re-analysis of the available data is important to uncover all features of small RNA sequencing.