Genome-wide association study and meta-analysis in Northern European populations replicate multiple colorectal cancer risk loci.

Genome-wide association studies have been successful in elucidating the genetic basis of colorectal cancer (CRC), but there remains unexplained variability in genetic risk. To identify new risk variants and to confirm reported associations, we conducted a genome-wide association study in 1,701 CRC cases and 14,082 cancer-free controls from the Finnish population. A total of 9,068,015 genetic variants were imputed and tested, and 30 promising variants were studied in additional 11,647 cases and 12,356 controls of European ancestry. The previously reported association between the single-nucleotide polymorphism (SNP) rs992157 (2q35) and CRC was independently replicated (p = 2.08 × 10-4 ; OR, 1.14; 95% CI, 1.06-1.23), and it was genome-wide significant in combined analysis (p = 1.50 × 10-9 ; OR, 1.12; 95% CI, 1.08-1.16). Variants at 2q35, 6p21.2, 8q23.3, 8q24.21, 10q22.3, 10q24.2, 11q13.4, 11q23.1, 14q22.2, 15q13.3, 18q21.1, 20p12.3 and 20q13.33 were associated with CRC in the Finnish population (false discovery rate < 0.1), but new risk loci were not found. These results replicate the effects of multiple loci on the risk of CRC and identify shared risk alleles between the Finnish population isolate and outbred populations.

Variation at 2q35 (PNKD and TMBIM1) influences colorectal cancer risk and identifies a pleiotropic effect with inflammatory bowel disease.

To identify new risk loci for colorectal cancer (CRC), we conducted a meta-analysis of seven genome-wide association studies (GWAS) with independent replication, totalling 13 656 CRC cases and 21 667 controls of European ancestry. The combined analysis identified a new risk association for CRC at 2q35 marked by rs992157 (P = 3.15 × 10-8, odds ratio = 1.10, 95% confidence interval = 1.06-1.13), which is intronic to PNKD (paroxysmal non-kinesigenic dyskinesia) and TMBIM1 (transmembrane BAX inhibitor motif containing 1). Intriguingly this susceptibility single-nucleotide polymorphism (SNP) is in strong linkage disequilibrium (r2 = 0.90, D' = 0.96) with the previously discovered GWAS SNP rs2382817 for inflammatory bowel disease (IBD). Following on from this observation we examined for pleiotropy, or shared genetic susceptibility, between CRC and the 200 established IBD risk loci, identifying an additional 11 significant associations (false discovery rate [FDR]) < 0.05). Our findings provide further insight into the biological basis of inherited genetic susceptibility to CRC, and identify risk factors that may influence the development of both CRC and IBD.

Mendelian randomisation implicates hyperlipidaemia as a risk factor for colorectal cancer.

While elevated blood cholesterol has been associated with an increased risk of colorectal cancer (CRC) in observational studies, causality is uncertain. Here we apply a Mendelian randomisation (MR) analysis to examine the potential causal relationship between lipid traits and CRC risk. We used single nucleotide polymorphisms (SNPs) associated with blood levels of total cholesterol (TC), triglyceride (TG), low-density lipoprotein (LDL), and high-density lipoprotein (HDL) as instrumental variables (IV). We calculated MR estimates for each risk factor with CRC using SNP-CRC associations from 9,254 cases and 18,386 controls. Genetically predicted higher TC was associated with an elevated risk of CRC (odds ratios (OR) per unit SD increase = 1.46, 95% confidence interval [CI]: 1.20-1.79, p = 1.68 × 10-4 ). The pooled ORs for LDL, HDL, and TG were 1.05 (95% CI: 0.92-1.18, p = 0.49), 0.94 (95% CI: 0.84-1.05, p = 0.27), and 0.98 (95% CI: 0.85-1.12, p = 0.75) respectively. A genetic risk score for 3-hydoxy-3-methylglutaryl-coenzyme A reductase (HMGCR) to mimic the effects of statin therapy was associated with a reduced CRC risk (OR = 0.69, 95% CI: 0.49-0.99, p = 0.046). This study supports a causal relationship between higher levels of TC with CRC risk, and a further rationale for implementing public health strategies to reduce the prevalence of hyperlipidaemia.

Mendelian randomisation analysis strongly implicates adiposity with risk of developing colorectal cancer.

BACKGROUND: Observational studies have associated adiposity with an increased risk of colorectal cancer (CRC). However, such studies do not establish a causal relationship. To minimise bias from confounding we performed a Mendelian randomisation (MR) analysis to examine the relationship between adiposity and CRC. METHODS: We used SNPs associated with adult body mass index (BMI), waist-hip ratio (WHR), childhood obesity and birth weight as instrumental variables in a MR analysis of 9254 CRC cases and 18 386 controls. RESULTS: In the MR analysis, the odds ratios (ORs) of CRC risk per unit increase in BMI, WHR and childhood obesity were 1.23 (95% CI: 1.02-1.49, P=0.033), 1.59 (95% CI: 1.08-2.34, P=0.019) and 1.07 (95% CI: 1.03-1.13, P=0.018), respectively. There was no evidence for association between birth weight and CRC (OR=1.22, 95% CI: 0.89-1.67, P=0.22). Combining these data with a concurrent MR-based analysis for BMI and WHR with CRC risk (totalling to 18 190 cases, 27 617 controls) provided increased support, ORs for BMI and WHR were 1.26 (95% CI: 1.10-1.44, P=7.7 × 10(-4)) and 1.40 (95% CI: 1.14-1.72, P=1.2 × 10(-3)), respectively. CONCLUSIONS: These data provide further evidence for a strong causal relationship between adiposity and the risk of developing CRC highlighting the urgent need for prevention and treatment of adiposity.

Eleven candidate susceptibility genes for common familial colorectal cancer.

Hereditary factors are presumed to play a role in one third of colorectal cancer (CRC) cases. However, in the majority of familial CRC cases the genetic basis of predisposition remains unexplained. This is particularly true for families with few affected individuals. To identify susceptibility genes for this common phenotype, we examined familial cases derived from a consecutive series of 1514 Finnish CRC patients. Ninety-six familial CRC patients with no previous diagnosis of a hereditary CRC syndrome were included in the analysis. Eighty-six patients had one affected first-degree relative, and ten patients had two or more. Exome sequencing was utilized to search for genes harboring putative loss-of-function variants, because such alterations are likely candidates for disease-causing mutations. Eleven genes with rare truncating variants in two or three familial CRC cases were identified: UACA, SFXN4, TWSG1, PSPH, NUDT7, ZNF490, PRSS37, CCDC18, PRADC1, MRPL3, and AKR1C4. Loss of heterozygosity was examined in all respective cancer samples, and was detected in seven occasions involving four of the candidate genes. In all seven occasions the wild-type allele was lost (P = 0.0078) providing additional evidence that these eleven genes are likely to include true culprits. The study provides a set of candidate predisposition genes which may explain a subset of common familial CRC. Additional genetic validation in other populations is required to provide firm evidence for causality, as well as to characterize the natural history of the respective phenotypes.

Exome sequencing reveals frequent inactivating mutations in ARID1A, ARID1B, ARID2 and ARID4A in microsatellite unstable colorectal cancer.

ARID1A has been identified as a novel tumor suppressor gene in ovarian cancer and subsequently in various other tumor types. ARID1A belongs to the ARID domain containing gene family, which comprises of 15 genes involved, for example, in transcriptional regulation, proliferation and chromatin remodeling. In this study, we used exome sequencing data to analyze the mutation frequency of all the ARID domain containing genes in 25 microsatellite unstable (MSI) colorectal cancers (CRCs) as a first systematic effort to characterize the mutation pattern of the whole ARID gene family. Genes which fulfilled the selection criteria in this discovery set (mutations in at least 4/25 [16%] samples, including at least one nonsense or splice site mutation) were chosen for further analysis in an independent validation set of 21 MSI CRCs. We found that in addition to ARID1A, which was mutated in 39% of the tumors (18/46), also ARID1B (13%, 6/46), ARID2 (13%, 6/46) and ARID4A (20%, 9/46) were frequently mutated. In all these genes, the mutations were distributed along the entire length of the gene, thus distinguishing them from typical MSI target genes previously described. Our results indicate that in addition to ARID1A, other members of the ARID gene family may play a role in MSI CRC.

Identification of candidate oncogenes in human colorectal cancers with microsatellite instability.

Microsatellite instability can be found in approximately 15% of all colorectal cancers. To detect new oncogenes we sequenced the exomes of 25 colorectal tumors and respective healthy colon tissue. Potential mutation hot spots were confirmed in 15 genes; ADAR, DCAF12L2, GLT1D1, ITGA7, MAP1B, MRGPRX4, PSRC1, RANBP2, RPS6KL1, SNCAIP, TCEAL6, TUBB6, WBP5, VEGFB, and ZBTB2; these were validated in 86 tumors with microsatellite instability. ZBTB2, RANBP2, and PSRC1 also were found to contain hot spot mutations in the validation set. The form of ZBTB2 associated with colorectal cancer increased cell proliferation. The mutation hot spots might be used to develop personalized tumor profiling and therapy.

Exome sequencing in diagnostic evaluation of colorectal cancer predisposition in young patients.

OBJECTIVE: Early-onset colorectal cancer (CRC), defined here as age of onset less than 40 years, develops frequently in genetically predisposed individuals. Next-generation sequencing is an increasingly available option in the diagnostic workup of suspected hereditary susceptibility, but little is known about the practical feasibility and additional diagnostic yield of the technology in this patient group. MATERIALS AND METHODS: We analyzed 38 young CRC patients derived from a set of 1514 CRC cases. All 38 tumors had been tested in our laboratory for microsatellite instability (MSI), and Sanger sequencing had been used to screen for MLH1 and MSH2 mutations in MSI cases. Also, gastrointestinal polyposis had been diagnosed clinically and molecularly. Family histories were acquired from national registries. If inherited syndromes had not been diagnosed in routine diagnostic efforts (n = 23), normal tissue DNA was analyzed for mutations in a comprehensive set of high-penetrance genes (MLH1, MSH2, MSH6, PMS2, APC, MUTYH, SMAD4, BMPR1A, LKB1/STK11, and PTEN) by exome sequencing. RESULTS: CRC predisposition syndromes were confirmed in 42% (16/38) of early-onset CRC patients. Hereditary nonpolyposis colorectal cancer was diagnosed in 12 (32%) patients, familial adenomatous polyposis in three (7.9%), and juvenile polyposis in one (2.6%) patient. Exome sequencing revealed one additional MLH1 mutation. Over half of the patients had advanced cancers (Dukes C or D, 61%, 23/38). The majority of nonsyndromic patients had unaffected first-degree relatives and microsatellite-stable tumors. CONCLUSIONS: Microsatellite instability positivity or gastrointestinal polyposis characterized all patients with unambiguous highly penetrant germline mutations. In our series, exome sequencing produced little added value in diagnosing the underlying predisposition conditions.

A Finnish founder mutation in RAD51D: analysis in breast, ovarian, prostate, and colorectal cancer.

BACKGROUND: RAD51D and RAD54L are involved in homologous recombination, and rare mutations in RAD51D were recently found in breast-ovarian cancer families. This study investigated RAD51D and RAD54L for mutations in breast and ovarian cancer patients in the Finnish population. METHODS: The study sequenced the RAD51D and RAD54L genes in 95 breast and/or ovarian cancer families and genotyped the identified mutation in an additional 2200 breast and 553 ovarian cancer patients and 2102 population controls. To investigate the role of the mutation in other common cancers, 1094 prostate and 980 colorectal cancer patients were genotyped. RESULTS: In the screening of RAD51D, one deleterious founder mutation c.576+1G>A was identified in two breast-ovarian cancer families. No mutations were found in RAD54L. Altogether, the c.576+1G>A mutation was detected in 5/707 patients with a personal or family history of ovarian cancer (OR 9.16, 95% CI 1.07 to 78.56; p=0.024), with the highest frequency among breast-ovarian cancer families (3/105 vs 1/1287 controls, OR 37.82, 95% CI 3.90 to 366.91; p=0.0016), but no elevated frequency among breast cancer patients/families (2/2105, p=1). The mutation was not found among prostate or colorectal cancer patients. CONCLUSIONS: The results of this study on familial and unselected breast, ovarian, colorectal, and prostate cancer patients suggest that RAD51D is primarily a moderate penetrance susceptibility gene for ovarian cancer, with clinical significance for the carriers.

Screening of Finnish RAD51C founder mutations in prostate and colorectal cancer patients.

BACKGROUND: Rare, heterozygous germline mutations in the RAD51C gene have been found in breast and ovarian cancer families. In the Finnish population, we have identified two founder mutations in RAD51C that increase the risk of ovarian cancer but not breast cancer in the absence of ovarian cancer. Risk for other cancers has not been studied. METHODS: To study the role of RAD51C mutations in other common cancer types, we genotyped the Finnish RAD51C founder mutations c.837 + 1G > A and c.93delG in 1083 prostate cancer patients and 802 colorectal cancer patients using TaqMan Real-Time PCR. RESULTS: No RAD51C mutations c.837 + 1G > A or c.93delG were detected among the prostate or colorectal cancer patients. CONCLUSIONS: The results suggest that the RAD51C mutations do not predispose to prostate or colorectal cancer.

High frequency of TTK mutations in microsatellite-unstable colorectal cancer and evaluation of their effect on spindle assembly checkpoint.

Frameshift mutations frequently accumulate in microsatellite-unstable colorectal cancers (MSI CRCs) typically leading to downregulation of the target genes due to nonsense-mediated messenger RNA decay. However, frameshift mutations that occur in the 3' end of the coding regions can escape decay, which has largely been ignored in previous works. In this study, we characterized nonsense-mediated decay-escaping frameshift mutations in MSI CRC in an unbiased, genome wide manner. Combining bioinformatic search with expression profiling, we identified genes that were predicted to escape decay after a deletion in a microsatellite repeat. These repeats, located in 258 genes, were initially sequenced in 30 MSI CRC samples. The mitotic checkpoint kinase TTK was found to harbor decay-escaping heterozygous mutations in exon 22 in 59% (105/179) of MSI CRCs, which is notably more than previously reported. Additional novel deletions were found in exon 5, raising the mutation frequency to 66%. The exon 22 of TTK contains an A(9)-G(4)-A(7) locus, in which the most common mutation was a mononucleotide deletion in the A(9) (c.2560delA). When compared with identical non-coding repeats, TTK was found to be mutated significantly more often than expected without selective advantage. Since TTK inhibition is known to induce override of the mitotic spindle assembly checkpoint (SAC), we challenged mutated cancer cells with the microtubule-stabilizing drug paclitaxel. No evidence of checkpoint weakening was observed. As a conclusion, heterozygous TTK mutations occur at a high frequency in MSI CRCs. Unexpectedly, the plausible selective advantage in tumourigenesis does not appear to be related to SAC.

Somatic mutations and germline sequence variants in patients with familial colorectal cancer.

It is estimated that up to 35% of colorectal cancers (CRC) can be explained by hereditary factors. However, genes predisposing to highly penetrant CRC syndromes account for only a small fraction of all cases. Thus, most CRCs still remain molecularly unexplained. A recent systematic sequencing study on well-annotated human protein coding genes identified 280 somatically mutated candidate cancer genes (CAN genes) in breast and colorectal cancer. It is estimated that 8% of all reported cancer genes show both somatic and germline mutations. Therefore, the identified CAN genes serve as a distinct set of candidates for being involved in hereditary susceptibility. The aim of this study was to evaluate the role of colorectal CAN genes in familial CRC. Samples from 45 familial CRCs without known cancer predisposing mutations were screened for somatic and germline variants in 15 top-ranked CAN genes. Six of the genes were found to be somatically mutated in our tumor series. We identified 22 nonsynonymous somatic mutations of which the majority was of missense type. In germline, three novel nonsynonymous variants were identified in the following genes: CSMD3, EPHB6 and C10orf137, and none of the variants were present in 890 population-matched healthy controls. It is possible that the identified germline variants modulate predisposition to CRC. Functional validation and larger sample sets, however, will be required to clarify the role of the identified germline variants in CRC susceptibility.

Comprehensive evaluation of coding region point mutations in microsatellite-unstable colorectal cancer.

Microsatellite instability (MSI) leads to accumulation of an excessive number of mutations in the genome, mostly small insertions and deletions. MSI colorectal cancers (CRCs), however, also contain more point mutations than microsatellite-stable (MSS) tumors, yet they have not been as comprehensively studied. To identify candidate driver genes affected by point mutations in MSI CRC, we ranked genes based on mutation significance while correcting for replication timing and gene expression utilizing an algorithm, MutSigCV Somatic point mutation data from the exome kit-targeted area from 24 exome-sequenced sporadic MSI CRCs and respective normals, and 12 whole-genome-sequenced sporadic MSI CRCs and respective normals were utilized. The top 73 genes were validated in 93 additional MSI CRCs. The MutSigCV ranking identified several well-established MSI CRC driver genes and provided additional evidence for previously proposed CRC candidate genes as well as shortlisted genes that have to our knowledge not been linked to CRC before. Two genes, SMARCB1 and STK38L, were also functionally scrutinized, providing evidence of a tumorigenic role, for SMARCB1 mutations in particular.

Comprehensive Evaluation of Protein Coding Mononucleotide Microsatellites in Microsatellite-Unstable Colorectal Cancer.

Approximately 15% of colorectal cancers exhibit microsatellite instability (MSI), which leads to accumulation of large numbers of small insertions and deletions (indels). Genes that provide growth advantage to cells via loss-of-function mutations in microsatellites are called MSI target genes. Several criteria to define these genes have been suggested, one of them being simple mutation frequency. Microsatellite mutation rate, however, depends on the length and nucleotide context of the microsatellite. Therefore, assessing the general impact of mismatch repair deficiency on the likelihood of mutation events is paramount when following this approach. To identify MSI target genes, we developed a statistical model for the somatic background indel mutation rate of microsatellites to assess mutation significance. Exome sequencing data of 24 MSI colorectal cancers revealed indels at 54 million mononucleotide microsatellites of three or more nucleotides in length. The top 105 microsatellites from 71 genes were further analyzed in 93 additional MSI colorectal cancers. Mutation significance and estimated clonality of mutations determined the most likely MSI target genes to be the aminoadipate-semialdehyde dehydrogenase AASDH and the solute transporter SLC9A8 Our findings offer a systematic profiling of the somatic background mutation rate in protein-coding mononucleotide microsatellites, allowing a full cataloging of the true targets of MSI in colorectal cancer. Cancer Res; 77(15); 4078-88. ©2017 AACR.

Pro-inflammatory fatty acid profile and colorectal cancer risk: A Mendelian randomisation analysis.

BACKGROUND: While dietary fat has been established as a risk factor for colorectal cancer (CRC), associations between fatty acids (FAs) and CRC have been inconsistent. Using Mendelian randomisation (MR), we sought to evaluate associations between polyunsaturated (PUFA), monounsaturated (MUFA) and saturated FAs (SFAs) and CRC risk. METHODS: We analysed genotype data on 9254 CRC cases and 18,386 controls of European ancestry. Externally weighted polygenic risk scores were generated and used to evaluate associations with CRC per one standard deviation increase in genetically defined plasma FA levels. RESULTS: Risk reduction was observed for oleic and palmitoleic MUFAs (OROA = 0.77, 95% CI: 0.65-0.92, P = 3.9 × 10-3; ORPOA = 0.36, 95% CI: 0.15-0.84, P = 0.018). PUFAs linoleic and arachidonic acid had negative and positive associations with CRC respectively (ORLA = 0.95, 95% CI: 0.93-0.98, P = 3.7 × 10-4; ORAA = 1.05, 95% CI: 1.02-1.07, P = 1.7 × 10-4). The SFA stearic acid was associated with increased CRC risk (ORSA = 1.17, 95% CI: 1.01-1.35, P = 0.041). CONCLUSION: Results from our analysis are broadly consistent with a pro-inflammatory FA profile having a detrimental effect in terms of CRC risk.

Systematic search for rare variants in Finnish early-onset colorectal cancer patients.

The heritability of colorectal cancer (CRC) is incompletely understood, and the contribution of undiscovered rare variants may be important. In search of rare disease-causing variants, we exome sequenced 22 CRC patients who were diagnosed before the age of 40 years. Exome sequencing data from 95 familial CRC patients were available as a validation set. Cases with known CRC syndromes were excluded. All patients were from Finland, a country known for its genetically homogenous population. We searched for rare nonsynonymous variants with allele frequencies below 0.1% in 3,374 Finnish and 58,112 non-Finnish controls. In addition, homozygous and compound heterozygous variants were studied. No genes with rare loss-of-function variants were present in more than one early-onset CRC patient. Three genes (ADAMTS4, CYTL1, and SYNE1) harbored rare loss-of-function variants in both early-onset and familial CRC cases. Five genes with homozygous variants in early-onset CRC cases were found (MCTP2, ARHGAP12, ATM, DONSON, and ROS1), including one gene (MCTP2) with a homozygous splice site variant. All discovered homozygous variants were exclusive to one early-onset CRC case. Independent replication is required to associate the discovered variants with CRC. These findings, together with a lack of family history in 19 of 22 (86%) early-onset patients, suggest genetic heterogeneity in unexplained early-onset CRC patients, thus emphasizing the requirement for large sample sizes and careful study designs to elucidate the role of rare variants in CRC susceptibility.

3'-UTR poly(T/U) repeat of EWSR1 is altered in microsatellite unstable colorectal cancer with nearly perfect sensitivity.

Approximately 15% of colorectal cancers exhibit instability of short nucleotide repeat regions, microsatellites. These tumors display a unique clinicopathologic profile and the microsatellite instability status is increasingly used to guide clinical management as it is known to predict better prognosis as well as resistance to certain chemotherapeutics. A panel of five repeats determined by the National Cancer Institute, the Bethesda panel, is currently the standard for determining the microsatellite instability status in colorectal cancer. Recently, a quasimonomorphic mononucleotide repeat 16T/U at the 3' untranslated region of the Ewing sarcoma breakpoint region 1 gene was reported to show perfect sensitivity and specificity in detecting mismatch repair deficient colorectal, endometrial, and gastric cancers in two independent populations. To confirm this finding, we replicated the analysis in 213 microsatellite unstable colorectal cancers from two independent populations, 148 microsatellite stable colorectal cancers, and the respective normal samples by PCR and fragment analysis. The repeat showed nearly perfect sensitivity for microsatellite unstable colorectal cancer as it was altered in 212 of the 213 microsatellite unstable (99.5%) and none of the microsatellite stable colorectal tumors. This repeat thus represents the first potential single marker for detecting microsatellite instability.

CTCF/cohesin-binding sites are frequently mutated in cancer.

Cohesin is present in almost all active enhancer regions, where it is associated with transcription factors. Cohesin frequently colocalizes with CTCF (CCCTC-binding factor), affecting genomic stability, expression and epigenetic homeostasis. Cohesin subunits are mutated in cancer, but CTCF/cohesin-binding sites (CBSs) in DNA have not been examined for mutations. Here we report frequent mutations at CBSs in cancers displaying a mutational signature where mutations in A•T base pairs predominate. Integration of whole-genome sequencing data from 213 colorectal cancer (CRC) samples and chromatin immunoprecipitation sequencing (ChIP-exo) data identified frequent point mutations at CBSs. In contrast, CRCs showing an ultramutator phenotype caused by defects in the exonuclease domain of DNA polymerase ɛ (POLE) displayed significantly fewer mutations at and adjacent to CBSs. Analysis of public data showed that multiple cancer types accumulate CBS mutations. CBSs are a major mutational hotspot in the noncoding cancer genome.

Characterization of the colorectal cancer-associated enhancer MYC-335 at 8q24: the role of rs67491583.

Recent genome-wide association studies have identified multiple regions at 8q24 that confer susceptibility to many cancers. In our previous work, we showed that the colorectal cancer (CRC) risk variant rs6983267 at 8q24 resides within a TCF4 binding site at the MYC-335 enhancer, with the risk allele G having a stronger binding capacity and Wnt responsiveness. Here, we searched for other potential functional variants within MYC-335. Genetic variation within MYC-335 was determined in samples from individuals of European, African, and Asian descent, with emphasis on variants in putative transcription factor binding sites. A 2-bp GA deletion rs67491583 was found to affect a growth factor independent (GFI) binding site and was present only in individuals with African ancestry. Chromatin immunoprecipitation performed in heterozygous cells showed that the GA deletion had an ability to reduce binding of the transcriptional repressors GFI1 and GFI1b. Screening of 1,027 African American colorectal cancer cases and 1,773 healthy controls did not reveal evidence for association (odds ratio: 1.17, 95% confidence interval: 0.97-1.41, P = 0.095). In this study, rs67491583 was identified as another functional variant in the CRC-associated enhancer MYC-335, but further studies are needed to establish the role of rs67491583 in the colorectal cancer predisposition of African Americans.