Prevalence of colorectal cancer (CRC) in the British Bangladeshi population (BAN) is low compared to British Caucasians (CAU). Genetic background may influence mutations and disease features.

We describe a statistical method for the characterization of genomic aberrations in single nucleotide polymorphism microarray data acquired from cancer genomes. Our approach allows us to model the joint effect of polyploidy, normal DNA contamination and intra-tumour heterogeneity within a single unified Bayesian framework. We demonstrate the efficacy of our method on numerous datasets including laboratory generated mixtures of normal-cancer cell lines and real primary tumours.

Aberrant activation of the SRC family kinase hematopoietic cell kinase (HCK) triggers hematological malignancies as a tumor cell-intrinsic oncogene. Here we find that high HCK levels correlate with reduced survival of colorectal cancer patients. Likewise, increased Hck activity in mice promotes the growth of endogenous colonic malignancies and of human colorectal cancer cell xenografts. Furthermore, tumor-associated macrophages of the corresponding tumors show a pronounced alternatively activated endotype, which occurs independently of mature lymphocytes or of Stat6-dependent Th2 cytokine signaling. Accordingly, pharmacological inhibition or genetic reduction of Hck activity suppresses alternative activation of tumor-associated macrophages and the growth of colon cancer xenografts. Thus, Hck may serve as a promising therapeutic target for solid malignancies.

Molecular indicators of colorectal cancer prognosis have been assessed in several studies, but most analyses have been restricted to a handful of markers. We aimed to identify prognostic biomarkers for colorectal cancer by sequencing panels of multiple driver genes.

In order to identify new genes with differential expression in early intestinal tumours, we performed mRNA (messenger ribonucleic acid) expression profiling of 16 human and 63 mouse adenomas. All individuals had germline APC mutations to ensure that tumorigenesis was driven by 'second hits' at APC. Using stringent filtering to identify changes consistent between humans and mice, we identified 60 genes up-regulated and 151 down-regulated in tumours. For 22 selected genes--including known Wnt targets--expression differences were confirmed by qRT-PCR (quantitative reverse transcription polymerase chain reaction). Most, but not all, differences were also present in colorectal carcinomas. In situ analysis showed a complex picture. Expression of up-regulated genes in adenomas was usually uniform/diffuse (e.g. ITGA6) or prominent in the tumour core (e.g. LGR5); in normal tissue, these genes were expressed at crypt bases or the transit amplifying zone. Down-regulated genes were often undetectable in adenomas, but in normal tissue were expressed in mesenchyme (e.g. GREM1/2) or differentiated cells towards crypt tops (e.g. SGK1). In silico analysis of TCF4-binding motifs showed that some of our genes were probably direct Wnt targets. Previous studies, mostly focused on human tumours, showed partial overlap with our 'expression signature', but 37 genes were unique to our study, including TACSTD2, SEMA3F, HOXA9 and IER3 (up-regulated), and TAGLN, GREM1, GREM2, MAB21L2 and RARRES2 (down-regulated). Combined analysis of our and published human data identified additional genes differentially expressed in adenomas, including decreased BMPs (bone morphogenetic proteins) and increased BUB1/BUB1B. Several of the newly identified, differentially expressed genes represent potential diagnostic or therapeutic targets for intestinal tumours.

In genome-wide association studies (GWASs) of colorectal cancer, we have identified two genomic regions in which pairs of tagging-single nucleotide polymorphisms (tagSNPs) are associated with disease; these comprise chromosomes 1q41 (rs6691170, rs6687758) and 12q13.13 (rs7163702, rs11169552). We investigated these regions further, aiming to determine whether they contain more than one independent association signal and/or to identify the SNPs most strongly associated with disease. Genotyping of additional sample sets at the original tagSNPs showed that, for both regions, the two tagSNPs were unlikely to identify a single haplotype on which the functional variation lay. Conversely, one of the pair of SNPs did not fully capture the association signal in each region. We therefore undertook more detailed analyses, using imputation, logistic regression, genealogical analysis using the GENECLUSTER program and haplotype analysis. In the 1q41 region, the SNP rs11118883 emerged as a strong candidate based on all these analyses, sufficient to account for the signals at both rs6691170 and rs6687758. rs11118883 lies within a region with strong evidence of transcriptional regulatory activity and has been associated with expression of PDGFRB mRNA. For 12q13.13, a complex situation was found: SNP rs7972465 showed stronger association than either rs11169552 or rs7136702, and GENECLUSTER found no good evidence for a two-SNP model. However, logistic regression and haplotype analyses supported a two-SNP model, in which a signal at the SNP rs706793 was added to that at rs11169552. Post-GWAS fine-mapping studies are challenging, but the use of multiple tools can assist in identifying candidate functional variants in at least some cases.

The presence of tumor immune infiltrating cells (TILs), particularly CD8(+) T-cells, is a robust predictor of outcome in patients with colorectal cancer (CRC). We revisited TIL abundance specifically in patients with microsatellite stable (MSS) CRC without evidence of lymph node or metastatic spread. Examination of the density of CD8(+) T-cells in primary tumors in the context of other pro-oncogenic markers was performed to investigate potential regulators of TILs. Two independent cohorts of patients with MSS T2-4N0M0 CRC, enriched for cases with atypical relapse, were investigated. We quantified CD8(+) and CD45RO(+) -TILs, inflammatory markers, NFkBp65, pStat3, Cyclo-oxygenase-2 (COX2) and GRP78 as well as transcription factors (TF), β-catenin and MYB. High CD8(+) TILs correlated with a better relapse-free survival in both cohorts (p = 0.002) with MYB and its target gene, GRP78 being higher in the relapse group (p = 0.001); no difference in pSTAT3 and p65 was observed. A mouse CRC (CT26) model was employed to evaluate the effect of MYB on GRP78 expression as well as T-cell infiltration. MYB over-expressing in CT26 cells increased GRP78 expression and the analysis of tumor-draining lymph nodes adjacent to tumors showed reduced T-cell activation. Furthermore, MYB over-expression reduced the efficacy of anti-PD-1 to modulate CT26 tumor growth. This high MYB and GRP78 show a reciprocal relationship with CD8(+) TILs which may be useful refining the prediction of patient outcome. These data reveal a new immunomodulatory function for MYB suggesting a basis for further development of anti-GRP78 and/or anti-MYB therapies.

The B-type Raf kinase (BRAF) V600E mutation is a well-established biomarker for poor prognosis in metastatic colorectal cancer (mCRC) and is a highly attractive drug target. A barrier to the development of new therapies targeting BRAF V600E in mCRC is the low prevalence of mutations (approximately 10 %) and the current need for access to sequencing-based technologies which are not routinely available outside of large cancer centres. Availability of a standardised immunohistochemistry (IHC) test, more suited to routine pathology practice, would provide much broader access to patient identification. We sought to evaluate the accuracy and clinical utility of a recently developed BRAF V600E IHC method as a prognostic biomarker in a large cohort of community-based CRC patients. Archival tumour samples from 505 patients with stage I-IV CRC were immunohistochemically tested with two antibodies, pBR1 for total BRAF and VE1 for BRAF V600E. Cases were assessed by two blinded pathologists, and results were compared to BRAF V600E mutation status determined using DNA sequencing. Discordant cases were retested with a BRAF V600E SNaPshot assay. BRAF mutation status was correlated with overall survival (OS) in stage IV CRC. By DNA sequencing and IHC, 505 and 477 patients were respectively evaluable. Out of 477 patients, 56 (11. 7 %) had BRAF V600E mutations detected by sequencing and 63 (13.2 %) by IHC. Using DNA sequencing results as the reference, sensitivity and specificity for IHC were 98.2 % (55/56) and 98.1 % (413/421), respectively. IHC had a positive predictive value (PPV) of 87.3 % (55/63) and a negative predictive value (NPV) of 99.8 % (413/414). Compared to DNA sequencing plus retesting of available discordant cases by SNaPshot assay, IHC using the VE1 antibody had a 100 % sensitivity (59/59), specificity (416/416), NPV (416/416) and PPV (59/59). Stage IV CRC patients with BRAF V600E protein detected by IHC exhibited a significantly shorter overall survival (hazard ratio = 2.20, 95 % CI 1.26-3.83, p = 0.005), consistent with other published series. Immunohistochemistry using the BRAF V600E VE1 antibody is an accurate diagnostic assay in CRC. The test provides a simple, clinically applicable method of testing for the BRAF V600E mutation in routine practice.

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.