Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Ultrasensitive Detection of Circulating LINE-1 ORF1p as a Specific Multicancer Biomarker

Improved biomarkers are needed for early cancer detection, risk stratification, treatment selection, and monitoring treatment response. Although proteins can be useful blood-based biomarkers, many have limited sensitivity or specificity for these applications. Long INterspersed Element-1 (LINE-1) open reading frame 1 protein (ORF1p) is a transposable element protein overexpressed in carcinomas and high-risk precursors during carcinogenesis with negligible expression in normal tissues, suggesting ORF1p could be a highly specific cancer biomarker. To explore ORF1p as a blood-based biomarker, we engineered ultrasensitive digital immunoassays that detect mid-attomolar (10-17 mol/L) ORF1p concentrations in plasma across multiple cancers with high specificity. Plasma ORF1p shows promise for early detection of ovarian cancer, improves diagnostic performance in a multianalyte panel, provides early therapeutic response monitoring in gastroesophageal cancers, and is prognostic for overall survival in gastroesophageal and colorectal cancers. Together, these observations nominate ORF1p as a multicancer biomarker with potential utility for disease detection and monitoring.

SIGNIFICANCE

The LINE-1 ORF1p transposon protein is pervasively expressed in many cancers and is a highly specific biomarker of multiple common, lethal carcinomas and their high-risk precursors in tissue and blood. Ultrasensitive ORF1p assays from as little as 25 μL plasma are novel, rapid, cost-effective tools in cancer detection and monitoring. See related commentary by Doucet and Cristofari, p. 2502. This article is featured in Selected Articles from This Issue, p. 2489.

Inherited mutations affecting the SRCAP complex are central in moderate-penetrance predisposition to uterine leiomyomas

Uterine leiomyomas (ULs) are benign smooth muscle tumors that are common in premenopausal women. Somatic alterations in MED12, HMGA2, FH, genes encoding subunits of the SRCAP complex, and genes involved in Cullin 3-RING E3 ligase neddylation are mutually exclusive UL drivers. Established predisposition genes explain only partially the estimated heritability of leiomyomas. Here, we examined loss-of-function variants across 18,899 genes in a cohort of 233,614 White European women, revealing variants in four genes encoding SRCAP complex subunits (YEATS4, ZNHIT1, DMAP1, and ACTL6A) with a significant association to ULs, and YEATS4 and ZNHIT1 strikingly rank first and second, respectively. Positive mutation status was also associated with younger age at diagnosis and hysterectomy. Moderate-penetrance UL risk was largely attributed to rare non-synonymous mutations affecting the SRCAP complex. To examine this disease phenotype more closely, we set out to identify inherited mutations affecting the SRCAP complex in our in-house sample collection of Finnish individuals with ULs (n = 860). We detected one individual with an ACTL6A splice-site mutation, two individuals with a YEATS4 missense mutation, and four individuals with DMAP1 mutations: one splice-site, one nonsense, and two missense variants. These individuals had large and/or multiple ULs, were often diagnosed at an early age, and many had family history of ULs. When a somatic second hit was found, ACTL6A and DMAP1 were silenced in tumors by somatic mutation and YEATS4 by promoter hypermethylation. Decreased H2A.Z staining was observed in the tumors, providing further evidence for the pathogenic nature of the germline mutations. Our results establish inactivation of genes encoding SRCAP complex subunits as a central contributor to moderate-penetrance UL predisposition.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

No evidence of EMAST in whole genome sequencing data from 248 colorectal cancers

Microsatellite instability (MSI) is caused by defective DNA mismatch repair (MMR), and manifests as accumulation of small insertions and deletions (indels) in short tandem repeats of the genome. Another form of repeat instability, elevated microsatellite alterations at selected tetranucleotide repeats (EMAST), has been suggested to occur in 50% to 60% of colorectal cancer (CRC), of which approximately one quarter are accounted for by MSI. Unlike for MSI, the criteria for defining EMAST is not consensual. EMAST CRCs have been suggested to form a distinct subset of CRCs that has been linked to a higher tumor stage, chronic inflammation, and poor prognosis. EMAST CRCs not exhibiting MSI have been proposed to show instability of di- and trinucleotide repeats in addition to tetranucleotide repeats, but lack instability of mononucleotide repeats. However, previous studies on EMAST have been based on targeted analysis of small sets of marker repeats, often in relatively few samples. To gain insight into tetranucleotide instability on a genome-wide level, we utilized whole genome sequencing data from 227 microsatellite stable (MSS) CRCs, 18 MSI CRCs, 3 POLE-mutated CRCs, and their corresponding normal samples. As expected, we observed tetranucleotide instability in all MSI CRCs, accompanied by instability of mono-, di-, and trinucleotide repeats. Among MSS CRCs, some tumors displayed more microsatellite mutations than others as a continuum, and no distinct subset of tumors with the previously proposed molecular characters of EMAST could be observed. Our results suggest that tetranucleotide repeat mutations in non-MSI CRCs represent stochastic mutation events rather than define a distinct CRC subclass.

Retrotransposon insertions can initiate colorectal cancer and are associated with poor survival

Genomic instability pathways in colorectal cancer (CRC) have been extensively studied, but the role of retrotransposition in colorectal carcinogenesis remains poorly understood. Although retrotransposons are usually repressed, they become active in several human cancers, in particular those of the gastrointestinal tract. Here we characterize retrotransposon insertions in 202 colorectal tumor whole genomes and investigate their associations with molecular and clinical characteristics. We find highly variable retrotransposon activity among tumors and identify recurrent insertions in 15 known cancer genes. In approximately 1% of the cases we identify insertions in APC, likely to be tumor-initiating events. Insertions are positively associated with the CpG island methylator phenotype and the genomic fraction of allelic imbalance. Clinically, high number of insertions is independently associated with poor disease-specific survival.

Retrotransposons are usually dormant in healthy tissue, but become activated during malignancy. Here, in colorectal cancer, Cajuso et al. show that retrotransposon activity associates with clinical features of the disease.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

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.

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.

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.

Abstract 1440: Germline loss-of-function alleles in Finnish colorectal cancer patients

Colorectal cancer, a major cause of cancer deaths worldwide, is often preventable by screening and surveillance, especially in patients at increased risk; however, its heritability is incompletely understood, and studies in isolated populations such as Finland may be valuable for gene discovery. We analyzed whole-genome or whole-exome sequencing data in 322 Finnish colorectal cancer cases of unknown etiology and 928 population-matched controls. To enrich for susceptibility alleles, cases were selected on the basis of early-onset colorectal cancer (age at diagnosis &lt;50 years; 50/322, 16%), synchronous or metachronous colorectal cancer (30/322, 9.3%), familial colorectal cancer (134/322, 42%), personal history of noncolorectal cancer (93/322, 29%), or family history of at least two noncolorectal cancers (124/322, 39%). To increase the prior probability of pathogenicity of variants, the study was restricted to rare and low-frequency protein-coding loss-of-function variants (nonsense, frameshift, and essential splice site variants; allele frequency &lt;3%). Single-variant and gene-level associations were tested, and genes with loss-of-function variants exclusively in cases were pinpointed. An independent set comprising 1,257 cases and 2,638 controls were genotyped using SNP arrays, and genotype imputation provided additional data on 225 (8.5%) of the 2,654 initially discovered loss-of-function variants. Principal component analysis was applied to correct for population structure. In association analyses, with a statistical power of at least 80% for a range of moderate-to-high effect sizes and low allele frequencies (e.g., odds ratio of 3 and allele frequency of 0.5%) at p &lt;0.15, a total of 35 candidate colorectal cancer predisposition genes were identified, although none of them reached exome-wide significance. There were 21 genes with rare loss-of-function alleles in at least two cases but none of the sequenced controls. Whole-genome sequencing data from paired normal and tumor tissues were available in 243 cases, and potential second hits (wild-type allele losses or somatic point mutations) were found in 15 genes. These results may, through further analysis and replication, contribute to our knowledge of the genetic architecture and pathobiology of colorectal cancer susceptibility.

Citation Format: Tomas Tanskanen, Ulrika Hänninen, Tatiana Cajuso, Johanna Kondelin, Riku Katainen, Jiri Hamberg, Niko Välimäki, Esa Pitkänen, Kimmo Palin, Lauri A. Aaltonen. Germline loss-of-function alleles in Finnish colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1440. doi:10.1158/1538-7445.AM2017-1440

Abstract 4381: The mobile genome of colorectal cancer: Characterization of retrotransposon insertions in 202 colorectal cancer whole genomes

The aim of this study is to characterize the landscape of retrotransposon insertions in 202 colorectal cancer (CRC) patients to elucidate their role in colorectal tumorigenesis. Retrotransposons are genetic sequences that can copy themselves into an mRNA intermediate and insert elsewhere in the genome. Germline retrotranspositions can contribute to genetic variation however, high retrotransposon activity can also lead to genetic instability. Albeit retrotransposons are usually repressed in normal adult tissues, they become highly active in several cancer types, being most active in epithelial cancers such as CRC. Moreover, retrotransposon insertions have been difficult to detect with previous methodological approaches, which has resulted in very few genome-wide studies.

Whole genome sequencing (WGS) was performed in 202 tumors and their corresponding normal tissue on the Illumina HiSeq 2000 platform; with 100 base-pairs paired-end reads. Each sample was sequenced to a minimum of 40x median coverage. To detect retrotransposon insertions we utilized TraFiC (Transposome Finder in Cancer). To identify retrotransposon-mediated transductions that were not identifiable using TraFiC, we utilized DELLY and the European database of L1-HS retrotransposon insertions in humans (euL1db). To characterize somatic insertions and decrease the rate of false positives, all insertions present in any of the normal genomes were filtered. After careful visual inspection of the paired-end read data, we estimated 80% of true somatic calls.

We detected high retrotransposon activity in CRC with remarkable variation among patients. The mean number of insertions per tumor was 25, ranging from 0 to 197 insertions. In addition, we were able to identify few highly active retrotransposons that accounted for 70% of the retrotransposon-mediated transductions. We identified a total of 5065 somatic insertions, 2% of these insertions were located in exons, whereas 45% were in introns. We identified several genes with recurrent insertions, some of these loci being known fragile sites. However, we identified 13 known cancer genes with two or more insertions. Furthermore, we identified two patients with insertions in exon 16 of APC, suggesting that these insertions could be initiating tumorigenic events. To conclude, we found retrotransposon insertions to be a common phenomenon in CRC, contributing to genome instability and colorectal tumorigenesis.

Citation Format: Tatiana Cajuso, Päivi Sulo, Riku Katainen, Ulrika Hänninen, Tomas Tanskanen, Johanna Kondelin, Jiri Hamberg, Niko Välimäki, Kimmo Palin, Outi Kilpivaara, Esa Pitkänen, Lauri A. Aaltonen. The mobile genome of colorectal cancer: Characterization of retrotransposon insertions in 202 colorectal cancer whole genomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4381. doi:10.1158/1538-7445.AM2017-4381

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.

NF-κB Mediates the Expression of TBX15 in Cancer Cells

TBX15 is a T-box transcription factor essential for development, also proposed as a marker in prostate cancer; and, recently, its antiapoptotic function indicates a role in carcinogenesis. Regulation of TBX15 is uncovered. In this study, we investigated the regulation of TBX15 expression in human cancer cells, by analyzing the regulatory function of a 5’-distal conserved region of TBX15. Bisulfite sequencing showed high methylation of the CpG island contained in this region that was not correlated with TBX15 mRNA levels, in the cancer cell lines analyzed; however, after 5-aza-dC treatment of TPC-1 cells an increase of TBX15 expression was observed. We also found a significant response of TBX15 to TNF-α activation of the NF-κB pathway using five cancer cell lines, and similar results were obtained when NF-κB was activated with PMA/ionomycin. Next, by luciferase reporter assays, we identified the TBX15 regulatory region containing two functional NF-κB binding sites with response to NF-κBp65, mapping on the -3302 and -3059 positions of the TBX15 gene. Moreover, a direct interaction of NF-κBp65 with one of the two NF-κB binding sites was indicated by ChIP assays. In summary, we provide novel data showing that NF-κB signaling up-regulates TBX15 expression in cancer cells. Furthermore, the link between TBX15 and NF-κB found in this study may be important to understand cancer and development processes.

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.

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.

Frequent L1 retrotranspositions originating from TTC28 in colorectal cancer

L1 element retrotranspositions have been found to alter expression of genes neighboring the insertion sites, potentially involving them in tumorigenesis and tumor progression. In colorectal cancer (CRC), L1 insertions have been found to target genes with a role in tumorigenesis. Structural changes such as L1 insertions are identifiable by whole genome sequencing (WGS). In this study, we observed frequent somatic L1 retrotranspositions originating from TTC28 using deep coverage WGS data from 92 CRC tumor-normal sample pairs. In two cases the event had targeted NOVA1 gene (p=0.025). In addition, a germline retrotransposition event from TTC28 to GABRA4 was found to be a common polymorphism in the Finnish population. Thus while some events may be tumorigenic, others are likely to be neutral. Our data contradict a recent study where a similar signal in TTC28 was interpreted as a common inactivating translocation. While much work remains to be performed to understand the biological significance of retrotranspositions in cancer, accurate identification of these events is a prerequisite for success.

Identification of 33 candidate oncogenes by screening for base-specific mutations

Background:

Genes with recurrent codon-specific somatic mutations are likely drivers of tumorigenesis and potential therapeutic targets. Hypermutable cancers may represent a sensitive system for generation and selection of oncogenic mutations.

Methods:

We utilised exome-sequencing data on 25 sporadic microsatellite-instable (MSI) colorectal cancers (CRCs) and searched for base-specific somatic mutation hotspots.

Results:

We identified novel mutation hotspots in 33 genes. Fourteen genes displayed mutations in the validation set of 254 MSI CRCs: ANTXR1, MORC2, CEP135, CRYBB1, GALNT9, KRT82, PI15, SLC36A1, CNTF, GLDC, MBTPS1, OR9Q2, R3HDM1 and TTPAL. A database search found examples of the hotspot mutations in multiple cancer types.

Conclusions:

This work reveals a variety of new recurrent candidate oncogene mutations to be further scrutinised as potential therapeutic targets.

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.