Colorectal cancer (CRC) is the third most prevalent cancer type worldwide with a mortality rate of approximately 50%. Elevated cell-surface expression of truncated carbohydrate structures such as Tn antigen (GalNAcα-Ser/Thr) is frequently observed during tumor progression. We have previously demonstrated that the C-type lectin macrophage galactose-type lectin (MGL), expressed by human antigen presenting cells, can distinguish healthy tissue from CRC through its specific recognition of Tn antigen. Both MGL binding and oncogenic BRAF mutations have been implicated in establishing an immunosuppressive microenvironment. Here we aimed to evaluate whether MGL ligand expression has prognostic value and whether this was correlated to BRAF(V600E) mutation status. Using a cohort of 386 colon cancer patients we demonstrate that high MGL binding to stage III tumors is associated with poor disease-free survival, independent of microsatellite instability or adjuvant chemotherapy. In vitro studies using CRC cell lines showed an association between MGL ligand expression and the presence of BRAF(V600E). Administration of specific BRAF(V600E) inhibitors resulted in decreased expression of MGL-binding glycans. Moreover, a positive correlation between induction of BRAF(V600E) and MGL binding to epithelial cells of the gastrointestinal tract was found in vivo using an inducible BRAF(V600E) mouse model. We conclude that the BRAF(V600E) mutation induces MGL ligand expression, thereby providing a direct link between oncogenic transformation and aberrant expression of immunosuppressive glycans. The strong prognostic value of MGL ligands in stage III colon cancer patients, i.e. when tumor cells disseminate to lymph nodes, further supports the putative immune evasive role of MGL ligands in metastatic disease.

Cancer is caused by somatic DNA alterations such as gene point mutations, DNA copy number aberrations (CNA) and structural variants (SVs). Genome-wide analyses of SVs in large sample series with well-documented clinical information are still scarce. Consequently, the impact of SVs on carcinogenesis and patient outcome remains poorly understood. This study aimed to perform a systematic analysis of genes that are affected by CNA-associated chromosomal breaks in colorectal cancer (CRC) and to determine the clinical relevance of recurrent breakpoint genes.

Colorectal cancer develops in a multi-step manner from normal epithelium, through a pre-malignant lesion (so-called adenoma), into a malignant lesion (carcinoma), which invades surrounding tissues and eventually can spread systemically (metastasis). It is estimated that only about 5% of adenomas do progress to a carcinoma.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

High-throughput molecular profiling techniques are routinely generating vast amounts of data for translational medicine studies. Secure access controlled systems are needed to manage, store, transfer and distribute these data due to its personally identifiable nature. The European Genome-phenome Archive (EGA) was created to facilitate access and management to long-term archival of bio-molecular data. Each data provider is responsible for ensuring a Data Access Committee is in place to grant access to data stored in the EGA. Moreover, the transfer of data during upload and download is encrypted. ELIXIR, a European research infrastructure for life-science data, initiated a project (2016 Human Data Implementation Study) to understand and document the ELIXIR requirements for secure management of controlled-access data. As part of this project, a full ecosystem was designed to connect archived raw experimental molecular profiling data with interpreted data and the computational workflows, using the CTMM Translational Research IT (CTMM-TraIT) infrastructure http://www.ctmm-trait.nl as an example. Here we present the first outcomes of this project, a framework to enable the download of EGA data to a Galaxy server in a secure way. Galaxy provides an intuitive user interface for molecular biologists and bioinformaticians to run and design data analysis workflows. More specifically, we developed a tool -- ega_download_streamer - that can download data securely from EGA into a Galaxy server, which can subsequently be further processed. This tool will allow a user within the browser to run an entire analysis containing sensitive data from EGA, and to make this analysis available for other researchers in a reproducible manner, as shown with a proof of concept study.  The tool ega_download_streamer is available in the Galaxy tool shed: https://toolshed.g2.bx.psu.edu/view/yhoogstrate/ega_download_streamer.

Secretory leukocyte protease inhibitor (SLPI), a pleiotropic protein expressed by healthy intestinal epithelial cells, functions as an inhibitor of NF-κB and neutrophil proteases and exerts antimicrobial activity. We previously showed SLPI suppresses intestinal epithelial chemokine production in response to microbial contact. Increased SLPI expression was recently detected in various types of carcinoma. In addition, accumulating evidence indicates SLPI expression is favorable for tumor cells. In view of these findings and the abundance of SLPI in the colonic epithelium, we hypothesized SLPI promotes colorectal cancer (CRC) growth and metastasis. Here, we aimed to establish whether SLPI expression in CRC is related to clinical outcome. Using a cohort of 507 patients with CRC who underwent resection of liver metastases, we show that high SLPI protein expression in both liver metastases and primary CRC is associated with significantly shorter overall survival after resection of liver metastases. The prognostic value of SLPI in CRC patients with liver metastases implies a role for SLPI in the formation of metastasis of human CRC. Based on the immune regulatory functions of SLPI, we anticipate that expression of SLPI provides tumors with a mechanism to evade infiltration by immune cells.

Development of cancer is driven by somatic alterations, including numerical and structural chromosomal aberrations. Currently, several computational methods are available and are widely applied to detect numerical copy number aberrations (CNAs) of chromosomal segments in tumor genomes. However, there is lack of computational methods that systematically detect structural chromosomal aberrations by virtue of the genomic location of CNA-associated chromosomal breaks and identify genes that appear non-randomly affected by chromosomal breakpoints across (large) series of tumor samples. 'GeneBreak' is developed to systematically identify genes recurrently affected by the genomic location of chromosomal CNA-associated breaks by a genome-wide approach, which can be applied to DNA copy number data obtained by array-Comparative Genomic Hybridization (CGH) or by (low-pass) whole genome sequencing (WGS). First, 'GeneBreak' collects the genomic locations of chromosomal CNA-associated breaks that were previously pinpointed by the segmentation algorithm that was applied to obtain CNA profiles. Next, a tailored annotation approach for breakpoint-to-gene mapping is implemented. Finally, dedicated cohort-based statistics is incorporated with correction for covariates that influence the probability to be a breakpoint gene. In addition, multiple testing correction is integrated to reveal recurrent breakpoint events. This easy-to-use algorithm, 'GeneBreak', is implemented in R ( www.cran.r-project.org) and is available from Bioconductor ( www.bioconductor.org/packages/release/bioc/html/GeneBreak.html).

Recurrence rates after resection of colorectal cancer liver metastases (CRLM) are high and correlate with worse survival. Postoperative circulating tumour DNA (ctDNA) is a promising prognostic biomarker. Focusing on patients with resected CRLM, this study aimed to evaluate the association between the detection of postoperative ctDNA, pathologic response and recurrence-free survival (RFS).

Gain of chromosome arm 13q is one of the most prevalent DNA copy number alterations associated with colorectal adenoma-to-carcinoma progression. The oncogenic miR-17-92 cluster, located at 13q, was found to be overexpressed in colorectal cancer and in adenomas harboring 13q gain. However, to what extent overexpression of this group of microRNAs actually drives progression to cancer remains to be resolved. Therefore, we aimed to clarify the role of miR-17-92 cluster in the progression from colorectal adenoma to carcinoma. The miR-17-92 cluster was overexpressed in human colorectal adenoma organoids without 13q gain and downstream effects on mRNA expression were investigated, along with functional consequences in vitro and in vivo. Comparison of mRNA sequencing results of organoids overexpressing miR-17-92 and cultures transduced with control vector revealed a miR-17-92 expression signature. This signature appeared to be enriched in an independent series of colorectal cancers and adenomas with 13q gain, confirming that miR-17-92 expression is associated with malignant progression. However, tumor-associated characteristics, such as increased proliferation rate, were not observed in miR-17-92 overexpressing adenoma organoids in vitro. In addition, subcutaneous injection of these organoids in immunodeficient mice was insufficient to cause tumor outgrowth. In conclusion, this study showed that miR-17-92 expression contributes to 13q gain-associated adenoma-to-carcinoma progression, however, this is insufficient to cause malignancy.

(1) Background: This study aimed to develop a machine learning model based on radiomics of pretreatment magnetic resonance imaging (MRI) 3D T2W contrast sequence scans combined with clinical parameters (CP) to predict neoadjuvant chemoradiotherapy (nCRT) response in patients with locally advanced rectal carcinoma (LARC). The study also assessed the impact of radiomics dimensionality on predictive performance. (2) Methods: Seventy-five patients were prospectively enrolled with clinicopathologically confirmed LARC and nCRT before surgery. Tumor properties were assessed by calculating 2141 radiomics features. Least absolute shrinkage selection operator (LASSO) and multivariate regression were used for feature selection. (3) Results: Two predictive models were constructed, one starting from 72 CP and 107 radiomics features, and the other from 72 CP and 1862 radiomics features. The models revealed moderately advantageous impact of increased dimensionality, with their predictive respective AUCs of 0.86 and 0.90 in the entire cohort and 0.84 within validation folds. Both models outperformed the CP-only model (AUC = 0.80) which served as the benchmark for predictive performance without radiomics. (4) Conclusions: Predictive models developed in this study combining pretreatment MRI radiomics and clinicopathological features may potentially provide a routine clinical predictor of chemoradiotherapy responders, enabling clinicians to personalize treatment strategies for rectal carcinoma.

Genetic predisposition controlled by susceptibility quantitative trait loci (QTLs) contributes to a large proportion of common cancers. Studies of genetics of cancer susceptibility, however, did not address systematically the relationship between susceptibility to cancers in different organs. We present five sets of data on genetic architecture of colon and lung cancer susceptibility in mice, humans and rats. They collectively show that the majority of genes for colon and lung cancer susceptibility are linked pair-wise and are likely identical or related. Four CcS/Dem recombinant congenic strains, each differing from strain BALB/cHeA by a different small random subset of ±12.5% of genes received from strain STS/A, suggestively show either extreme susceptibility or extreme resistance for both colon and lung tumors, which is unlikely if the two tumors were controlled by independent susceptibility genes. Indeed, susceptibility to lung cancer (Sluc) loci underlying the extreme susceptibility or resistance of such CcS/Dem strains, mapped in 226 (CcS-10 x CcS-19)F2 mice, co-localize with susceptibility to colon cancer (Scc) loci. Analysis of additional Sluc loci that were mapped in OcB/Dem strains and Scc loci in CcS/Dem strains, respectively, shows their widespread pair-wise co-localization (P  =  0.0036). Finally, the majority of published human and rat colon cancer susceptibility genes map to chromosomal regions homologous to mouse Sluc loci. 12/12 mouse Scc loci, 9/11 human and 5/7 rat colon cancer susceptibility loci are close to a Sluc locus or its homologous site, forming 21 clusters of lung and colon cancer susceptibility genes from one, two or three species. Our data shows that cancer susceptibility QTLs can have much broader biological effects than presently appreciated. It also demonstrates the power of mouse genetics to predict human susceptibility genes. Comparison of molecular mechanisms of susceptibility genes that are organ-specific and those with trans-organ effects can provide a new dimension in understanding individual cancer susceptibility.

Colorectal adenomas form a biologically and clinically distinct intermediate stage in development of colorectal cancer (CRC) from normal colon epithelium. Only 5% of adenomas progress into adenocarcinomas, indicating that malignant transformation requires other biological alterations than those involved in adenoma formation. The present study aimed to explore which cancer-related biological processes are affected during colorectal adenoma-to-carcinoma progression and to identify key genes within these pathways that can serve as tumor markers for malignant transformation. The activity of 12 cancer-related biological processes was compared between 37 colorectal adenomas and 31 adenocarcinomas, using the pathway analysis tool Gene Set Enrichment Analysis. Expression of six gene sets was significantly increased in CRCs compared to adenomas, representing chromosomal instability, proliferation, differentiation, invasion, stroma activation, and angiogenesis. In addition, 18 key genes were identified for these processes based on their significantly increased expression levels. For AURKA and PDGFRB, increased mRNA expression levels were verified at the protein level by immunohistochemical analysis of a series of adenomas and CRCs. This study revealed cancer-related biological processes whose activities are increased during malignant transformation and identified key genes which may be used as tumor markers to improve molecular characterization of colorectal tumors.

Upfront cytoreductive surgery with HIPEC (CRS-HIPEC) is the standard treatment for isolated resectable colorectal peritoneal metastases (PM) in the Netherlands. This study investigates whether addition of perioperative systemic therapy to CRS-HIPEC improves oncological outcomes.

Non-invasive approaches for cell-free DNA (cfDNA) assessment provide an opportunity for cancer detection and intervention. Here, we use a machine learning model for detecting tumor-derived cfDNA through genome-wide analyses of cfDNA fragmentation in a prospective study of 365 individuals at risk for lung cancer. We validate the cancer detection model using an independent cohort of 385 non-cancer individuals and 46 lung cancer patients. Combining fragmentation features, clinical risk factors, and CEA levels, followed by CT imaging, detected 94% of patients with cancer across stages and subtypes, including 91% of stage I/II and 96% of stage III/IV, at 80% specificity. Genome-wide fragmentation profiles across ~13,000 ASCL1 transcription factor binding sites distinguished individuals with small cell lung cancer from those with non-small cell lung cancer with high accuracy (AUC = 0.98). A higher fragmentation score represented an independent prognostic indicator of survival. This approach provides a facile avenue for non-invasive detection of lung cancer.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) have been identified in bacteria, archaea and mitochondria of plants, but not in eukaryotes. Here, we report the discovery of 12,572 putative CRISPRs randomly distributed across the human chromosomes, which we termed hCRISPRs. By using available transcriptome datasets, we demonstrate that hCRISPRs are distinctively expressed as small non-coding RNAs (sncRNAs) in cell lines and human tissues. Moreover, expression patterns thereof enabled us to distinguish normal from malignant tissues. In prostate cancer, we confirmed the differential hCRISPR expression between normal adjacent and malignant primary prostate tissue by RT-qPCR and demonstrate that the SHERLOCK and DETECTR dipstick tools are suitable to detect these sncRNAs. We anticipate that the discovery of CRISPRs in the human genome can be further exploited for diagnostic purposes in cancer and other medical conditions, which certainly will lead to the development of point-of-care tests based on the differential expression of the hCRISPRs.

Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non-poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints.

Macrophages are versatile cells, which phenotype is profoundly influenced by their environment. Pro-inflammatory classically activated or M1 macrophages, and anti-inflammatory alternatively-activated or M2 macrophages represent two extremes of a continuum of functional states. Consequently, macrophages that are present in tumors can exert tumor-promoting and tumor-suppressing activity, depending on the tumor milieu. In this study we investigated how human monocytes-the precursors of macrophages-are influenced by carcinoma cells of different origin. We demonstrate that monocytes, stimulated with breast cancer supernatant, showed increased expression of interleukin (IL)-10, IL-8 and chemokines CCL17 and CCL22, which are associated with an alternatively-activated phenotype. By contrast, monocytes that were cultured in supernatants of colon cancer cells produced more pro-inflammatory cytokines (e.g., IL-12 and TNFα) and reactive oxygen species. Secretome analysis revealed differential secretion of proteins by colon and breast cancer cell lines, of which the proteoglycan versican was exclusively secreted by colon carcinoma cell lines. Reducing active versican by blocking with monoclonal antibodies or shRNA diminished pro-inflammatory cytokine production by monocytes. Thus, colon carcinoma cells polarize monocytes toward a more classically-activated anti-tumorigenic phenotype, whereas breast carcinomas predispose monocytes toward an alternatively activated phenotype. Interestingly, presence of macrophages in breast or colon carcinomas correlates with poor or good prognosis in patients, respectively. The observed discrepancy in macrophage activation by either colon or breast carcinoma cells may therefore explain the dichotomy between patient prognosis and macrophage presence in these different tumors. Designing new therapies, directing development of monocytes toward M1 activated tumor macrophages in cancer patients, may have great clinical benefits.

The N-Myc Downstream-Regulated Gene 4 (NDRG4), a prominent biomarker for colorectal cancer (CRC), is specifically expressed by enteric neurons. Considering that nerves are important members of the tumor microenvironment, we here establish different Ndrg4 knockout (Ndrg4-/- ) CRC models and an indirect co-culture of primary enteric nervous system (ENS) cells and intestinal organoids to identify whether the ENS, via NDRG4, affects intestinal tumorigenesis. Linking immunostainings and gastrointestinal motility (GI) assays, we show that the absence of Ndrg4 does not trigger any functional or morphological GI abnormalities. However, combining in vivo, in vitro, and quantitative proteomics data, we uncover that Ndrg4 knockdown is associated with enlarged intestinal adenoma development and that organoid growth is boosted by the Ndrg4-/- ENS cell secretome, which is enriched for Nidogen-1 (Nid1) and Fibulin-2 (Fbln2). Moreover, NID1 and FBLN2 are expressed in enteric neurons, enhance migration capacities of CRC cells, and are enriched in human CRC secretomes. Hence, we provide evidence that the ENS, via loss of Ndrg4, is involved in colorectal pathogenesis and that ENS-derived Nidogen-1 and Fibulin-2 enhance colorectal carcinogenesis.

Removal of colorectal adenomas is an effective strategy to reduce colorectal cancer (CRC) mortality rates. However, as only a minority of adenomas progress to cancer, such strategies may lead to overtreatment. The present study aimed to characterize adenomas by in-depth molecular profiling, to obtain insights into altered biology associated with the colorectal adenoma-to-carcinoma progression. We obtained low-coverage whole genome sequencing, RNA sequencing and tandem mass spectrometry data for 30 CRCs, 30 adenomas and 18 normal adjacent colon samples. These data were used for DNA copy number aberrations profiling, differential expression, gene set enrichment and gene-dosage effect analysis. Protein expression was independently validated by immunohistochemistry on tissue microarrays and in patient-derived colorectal adenoma organoids. Stroma percentage was determined by digital image analysis of tissue sections. Twenty-four out of 30 adenomas could be unambiguously classified as high risk (n = 9) or low risk (n = 15) of progressing to cancer, based on DNA copy number profiles. Biological processes more prevalent in high-risk than low-risk adenomas were related to proliferation, tumor microenvironment and Notch, Wnt, PI3K/AKT/mTOR and Hedgehog signaling, while metabolic processes and protein secretion were enriched in low-risk adenomas. DNA copy number driven gene-dosage effect in high-risk adenomas and cancers was observed for POFUT1, RPRD1B and EIF6. Increased POFUT1 expression in high-risk adenomas was validated in tissue samples and organoids. High POFUT1 expression was also associated with Notch signaling enrichment and with decreased goblet cells differentiation. In-depth molecular characterization of colorectal adenomas revealed POFUT1 and Notch signaling as potential drivers of tumor progression.

The core spliceosomal Sm proteins were recently proposed as cancer-selective lethal targets in non-small cell lung cancer (NSCLC). In contrast, the loss of the commonly mutated cancer target SF3B1 appeared to be toxic to non-malignant cells as well. In the current study, the transcriptomes of A549 NSCLC cells, in which SF3B1 or SNRPD3 was silenced, were compared using RNA sequencing. The skipping of exon 4 of the proteasomal subunit beta type-3 (PSMB3) mRNA, resulting in a shorter PSMB3-S variant, occurred only after silencing SNRPD3. This observation was extended to the other six Sm genes. Remarkably, the alternative splicing of PSMB3 mRNA upon Sm gene silencing was not observed in non-malignant IMR-90 lung fibroblasts. Furthermore, PSMB3 was found to be overexpressed in NSCLC clinical samples and PSMB3 expression correlated with Sm gene expression. Moreover, a high PSMB3 expression corresponds to worse survival in patients with lung adenocarcinomas. Finally, silencing the canonical full-length PSMB3-L, but not the shorter PSMB3-S variant, was cytotoxic and was accompanied by a decrease in proteasomal activity. Together, silencing Sm genes, but not SF3B1, causes a cytotoxic alternative splicing switch in the PSMB3 mRNA in NSCLC cells only.