Skip to main content

Exome sequencing characterizes the somatic mutation spectrum of early serrated lesions in a patient with serrated polyposis syndrome (SPS)



Serrated or Hyperplastic Polyposis Syndrome (SPS, HPS) is a yet poorly defined colorectal cancer (CRC) predisposition characterised by the occurrence of multiple and/or large serrated polyps throughout the colon. A serrated polyp-CRC sequence (serrated pathway) of CRC formation has been postulated, however, to date only few molecular signatures of serrated neoplasia (BRAF, KRAS, RNF43 mutations, CpG Island Methylation, MSI) have been described in a subset of SPS patients and neither the etiology of the syndrome nor the distinct genetic alterations during tumorigenesis have been identified.


To identify somatic point mutations in potential novel candidate genes of SPS-associated lesions and the involved pathways we performed exome sequencing of eleven early serrated polyps obtained from a 41 year-old female patient with clinically confirmed SPS. For data filtering and analysis, standard pipelines were used. Somatic mutations were identified by comparison with leukocyte DNA and were validated by Sanger sequencing.


The BRAF p.V600E or KRAS p.G12D mutation was identified in six polyps (~50%) and not found in polyps from the distal colon. In addition, we found seven unique rare somatic alterations of seven different genes in four serrated tumours, all of which are missense variants. The variant in ABI3BP and CATSPERB are predicted to be deleterious. No established cancer gene or candidate genes related to serrated tumorigenesis were affected.


Somatic mutations seem to be rare events in early hyperplastic and serrated lesions of SPS patients. Neither frequently affected genes nor enrichment of specific pathways were observed. Thus, other alterations such as non-coding variants or epigenetic changes might be the major driving force of tumour progression in SPS.

Peer Review reports


Serrated colorectal polyps represent a heterogeneous group of lesions that includes hyperplastic polyps (HPs), sessile serrated polyps/adenomas (SSP/A) with or without dysplasia, traditional serrated adenomas (TSA) and mixed polyps with a combination of different characteristics [1, 2].

HPs belong to the most common colorectal polyp types and were traditionally regarded as harmless, non-neoplastic lesions, however, there is growing evidence that at least some HPs have a malignant potential [3, 4]. According to this model, a small percentage of HPs –often large and right-sided – progress to other serrated polyps, in particular SSP/A, and then evolve to colorectal cancer (CRC). Commonly, HPs are smaller than 5 mm in diameter and tend to appear more frequently in the distal colon and rectum whereas SSA/Ps are found more often in the right than the left colon and are larger than 5 mm in diameter [5].

Nowadays, it is assumed that serrated polyps are part of an alternative (serrated) pathway of colorectal carcinogenesis. Similar to the classical adenoma-adenocarcinoma sequence [6], a serrated polyp-carcinoma sequence (serrated neoplasia pathway) was postulated where serrated polyps replace the traditional adenoma as the precursor lesion to CRC [7]. It is assumed, that serrated adenocarcinomas, which represent approximately 10% of sporadic CRC, arise via the serrated polyp-carcinoma sequence.

While the step-wise progression of premalignant lesions to carcinomas in the classical adenoma-carcinoma sequence is characterised by chromosomal instability (CIN) and APC or KRAS mutations, the molecular profiles described in serrated neoplasias include specific BRAF and KRAS mutations, microsatellite instability (MSI), and a CpG island methylator phenotype (CIMP) [8,9,10,11]. Nonetheless, although the serrated pathway is regarded to be the second most important pathway leading to CRC, the molecular steps of tumorigenesis are still largely unknown. Previous studies identified two potential RNA markers for SSP/As (ANXA10, CLDN1) by microarray analyses [12, 13] and somatic mutations in established cancer pathways [14]. Very recently, a unique signature of differentially expressed genes was found in early SSP/As by a comprehensive RNA-Seq analysis that discriminates between SSP/As and HPs [15]. Interestingly, an almost complete gene overlap between sporadic and syndromic SSP/As was observed.

Hyperplastic polyposis syndrome (HPS), now mostly designated as serrated polyposis syndrome (SPS), is a yet clinically defined colorectal polyposis characterised by the presence of multiple and/or large serrated polyps throughout the large intestine. In this way, SPS can be regarded as a model disease of the serrated neoplasia pathway [11]. The empirical WHO diagnostic criteria require 1) at least 5 serrated polyps proximal to the sigmoid colon, with two or more >10 mm, or 2) any number of serrated polyps proximal to the sigmoid colon with an affected 1° relative, or 3) >20 serrated polyps of any size throughout the colon [1].

If left untreated, affected persons and their relatives have an increased lifetime risk for the development of CRC. Therefore, SPS patients undergo endoscopic surveillance with removal of polyps or surgical colonic resection. The prognosis of these early-manifesting lesions can be decisively improved through the timely detection of the tumour by established disease specific surveillance programs, however, the optimal extent and frequency of surveillance is unknown so far. Interestingly, CRC even occurred in patients despite all SSP/As were excised [16].

In contrast to many other inherited gastrointestinal polyposis syndromes, the etiology of SPS remains unclear. Although SPS was initially considered to be non-inherited, familial clustering and a high risk (up to 50%) of CRC in first degree relatives of SPS patients has been described [17]. Based on these reports, a yet unidentified underlying genetic defect seems to play a significant role in SPS, however, the possible modes of inheritance are still unknown. A previous study proposed germline mutations in oncogene-induced senescence pathways as predisposing factor [18]. In addition, both truncating somatic and germline mutations of RNF43 were found in SSA and a few patients with (autosomal dominant) SPS, respectively [18, 19].

The characteristic molecular profile of serrated neoplasia mentioned above was described in a subset of SPS-associated polyps. Probably, the current definition of SPS subsumes several clinical subtypes with different cancer risks and prognosis, however, in the absence of valid biomarkers no distinction between those subtypes is possible yet.

As the etiology of the syndrome and the distinct genetic alterations during tumorigenesis are still poorly understood, this study aimed to characterise the spectrum of somatic mutations in protein-coding genes present in early serrated lesions from different colonic parts of a single SPS patient.


Patient / data collection

The female patient presented at 38 years of age due to unclear abdominal trouble. A colonoscopy resulted in the diagnosis of a colorectal polyposis. Initially, several flat, broad-based polyps resembling a dense lawn were found within the caecum in addition to at least 10 more, partly flat, partly pedunculated polyps with a diameter of up to 20 mm. Four more colonoscopies were performed; altogether more than 45 polyps covering the whole colon were identified and partly removed. Around 20 polyps were histologically examined, the vast majority was classified as hyperplastic polyp and a few as serrated polyp, all without severe dysplasia. Neither typical adenomatous polyps nor carcinomas were found. A gastroduodenoscopy showed normal results. The WHO diagnostic criteria for SPS were met by the patient. A MUTYH-associated polyposis (MAP) was excluded by screening of the MUTYH gene as described elsewhere [20].

At the age of 41 years, a subtotal colectomy with ileoanal anastomosis was performed. The removed colon was examined carefully by an experienced pathologist from the Institute of Pathology in Cologne, the reference pathology of the German HNPCC consortium. Overall, 21 polyps up to 18 mm were seen, removed, and immediately fresh-frozen; all of them were classified as hyperplastic or sessile serrated polyps later on.

Beside of this, the patient had an inconspicuous medical history without important pre-existing disease. The family history was unremarkable, in particular regarding gastrointestinal polyps, extraintestinal tumours, or cancer, however, the father and his three half siblings had no colonoscopy.

The study was approved by the local ethics review board (Medical Faculty of the University of Bonn ethics review board no. 208/12), and a written informed consent was obtained from the patient.

DNA extraction from blood and tissue samples

Leukocyte-derived genomic DNA was isolated from peripheral blood with standard protocols as described elsewhere [21]. DNA from macrodissected, fresh-frozen serrated polyp tissue was extracted using the proteinase K treatment followed phenol extraction as described elsewhere [22].

Targeted genetic profiling of polyps

The targeted mutation screening was performed with the commercial kit TruSight®Tumor 15 Sequencing Panel (Illumina, San Diego) which includes 15 established cancer-associated genes (AKT1, BRAF, EGFR, ERBB2, FOXL2, GNA11, GNAQ, KIT, KRAS, MET, NRAS, PDGFRA, PIK3CA, RET, TP53). DNA concentration was determined with the Quanti Fluor dsDNA Systems Assay Kit (Promega, Mannheim, Germany) on the Quantus Fluorometer (Promega). Twenty nanogram DNA was used for library preparation. Target enrichment and high-throughput sequencing was performed according to the manufacturers protocol. All samples were sequenced on an Illumina MiSeq sequencer (Illumina, San Diego, CA). Pooled libraries (5 ng/μl) were spiked with 1% PhiX DNA (Illumina) and paired-end sequenced was performed with the “MiSeq reagent Kit V3 (600-cycles)” (Illumina). FastQ files generated by the MiSeq Reporter were used as data output. All variants identified were confirmed by high resolution melting (HRM) analysis and by Sanger sequencing as previously reported [23].

Microsatellite analysis

Microsatellite analysis was performed on matched tumour and normal DNA samples. This involved use of the National Cancer Institute (NCI) reference marker panel for the evaluation of microsatellite instability (MSI) in colorectal cancer. This panel consists of two mononucleotide (BAT25, BAT26), and three dinucleotide (D2S123, D5S346 and D17S250) repeats [24]. Tumour DNA was extracted from macrodissected tumour tissue. Normal DNA was extracted from peripheral blood leukocytes. Tumours were scored as highly unstable (MSI-H) if two or more of these five markers exhibited additional alleles, and as stable (MSS) if none of the five markers showed instability.

Whole exome sequencing (WES)

For whole-exome sequencing, 1 μg of DNA was fragmented with sonication technology (Bioruptor, Diagenode). The fragments were end-repaired and adaptor-ligated, including incorporation of sample index barcodes. After size selection, a pool of all 6 libraries was subjected to an enrichment process with the SeqCap EZ Human Exome Library version 3.0 kit (Roche NimbleGen). The final libraries were sequenced on an Illumina HiSeq 2000 sequencing instrument (Illumina, San Diego, U.S.) with a paired-end 2 × 100 bp protocol.

Exome analysis and filtering

Primary data were filtered using the Illumina Realtime Analysis (RTA) software version 1.8. Subsequently, the reads were mapped to the human genome reference build GRCh37 using the BWA alignment algorithm [25]. GATK version 1.6 [26] was used to mark duplicated reads, perform local realignment around short insertions and deletions, recalibrate the base quality scores and call SNVs and short Indels together with SAMtools version 0.1.18 and Dindel version 1.01 [27]. On average, 104.7 million reads were uniquely mapped. 80.9% of all targeted regions were covered >30 times, with a mean coverage of 89. Scripts developed in-house at the Cologne Center for Genomics were applied to detect protein changes, affected donor and acceptor splice sites, and overlaps with known variants. Acceptor and donor splice site mutations were analyzed using a Maximum Entropy model [28].

Variant filtering was performed with the VARBANK graphical user interface ( ), comparison with public variant databases and subsequent manual check to select relevant, likely pathogenic somatic mutations affecting protein coding genes. Variants were included if they 1) are truncating variants (nonsense mutations, frameshift deletions/insertions, mutations located at exon-flanking, highly conserved intronic splice sites), or apparent missense mutations predicted to be pathogenic by in-silico prediction tools (see below); 2) are located at a region with a minimal total coverage (read depth) of at least 10×; 3) had a variant allele frequency of 35% in tumour DNA; and 4) were not detected in corresponding matching normal tissue (leukocyte DNA).

Subsequently, all common variants which are reported with a minor allele frequency (MAF) of >1% in the germline according to population-based databases (NCBI dbSNP, ExAC, 1000Genomes) were removed since those variants are very likely to represent rare polymorphisms or low-penetrant variants rather than pathogenic somatic driver mutations of tumorigenesis. In addition, filtering was performed against an in-house database containing variants from 511 exomes from epilepsy patients in order to exclude pipeline-related artifacts or population-specific rare polymorphisms. The remaining variants were further checked by visual inspection using the VARBANK read browser to exclude obvious false positive variants / artifacts.

Afterwards, 123 genes which had been proposed as driver candidate genes of the serrated pathway due to genetic alterations or epigenetic silencing in serrated lesions, including established causative genes or candidate genes for (hereditary) colorectal tumours (Additional file 1: Table S1) [19, 29,30,31,32,33,34,35,36] were analyzed for the presence of any rare variant.

Sanger sequencing

All rare variants in promising candidate genes were validated by Sanger sequencing of tumour DNA, assuming that relevant somatic mutations are present in an allelic fraction of ≥10%. DNA sequences were obtained from the UCSC Genome Browser (hg19). Primers were designed using Primer3 V.0.4.0 ( ) (Additional file 1: Table S2). Stored genomic DNA was extracted from leukocytes and tumour DNA was used to amplify the coding regions and adjacent intronic sequences of the respective genes. PCR products were purified using the QIA quick PCR purification kit (Qiagen, Hilden, Germany), and sequenced on an ABI 3500xl Genetic Analyzer (Life Technologies) using the BigDye terminator kit version 1.1 (Life Technologies). The cDNA bases were numbered according to the gene reference sequence in GenBank, where 1 corresponds to the A of the ATG translation initiation codon.

In-silico analysis

Splicing efficiencies of the normal and mutant sequences were calculated using the splice prediction program NNSPLICE 0.9 from BDGP (the Berkeley Drosophila Genome Project). Potentially deleterious effects of putative missense variants were predicted using the variant-based in-silico tools CADD (cutoff on deleteriousness: values >15) [37], PolyPhen-2, MutationTaster, and SIFT, as well as the functional gene constraint z- and %ExAC_RVIS-scores [38, 39]. For all rare interesting candidate variants, the allele frequencies in controls were again checked using the gnomAD [38], COSMIC, and TCGA databases.


Eleven representative polyps from all parts of the colon had been selected which cover a spectrum from very small (4 mm) to medium-sized (13 mm) lesions. All were classified as hyperplastic or sessile serrated polyps (Table 1). All examined polyps demonstrated microsatellite stability (MSS) (data not shown).

Table 1 Results of histopathology and targeted molecular profiling of the eleven polyps

Prior to further filtering and validation steps of the exome sequencing data, we looked for variants in 123 genes including RNF43 and other driver candidate genes of the serrated pathway and established causative genes or candidate genes for (hereditary) colorectal tumours. Thereby we identified the two common hotspot mutations BRAF V600E (c.1799 T > A;pVal600Glu; NM_004333.4) and KRAS G12D (c.35G > A;p.Gly12Asp; NM_033360.3) (8–17% of mutated reads) in 6/11 tumour samples (Table 1). No further putative pathogenic variant could be identified in any of the remaining genes although the respective genomic regions had a minimum of 30× coverage in all polyp-derived DNA samples tested except for sample T69 which showed a lower coverage in general, and very few single exons of some genes with a coverage of <25× while RNF43 was covered higher than 100× in all samples and exons.

Targeted genetic profiling of all polyps for 15 established cancer-associated genes (TruSight®Tumor15, Illumina) confirmed the BRAF V600E and KRAS G12D mutations in 5 (46%) and one (9%) polyps, respectively. No mutation was identified in the remaining 13 cancer-associated genes. All mutations could be further validated by Sanger sequencing (Additional file 2: Figure S1).

Fig. 1

Sanger sequencing confirmed seven somatic missense mutations (upper row: tumour DNA; lower row: leukocyte DNA of the patient; for COL8A1, DNAI1 and VGLL2 reverse sequences are shown)

After applying the above mentioned filter criteria, the vast majority of variants was excluded. All of the remaining 25 unique variants were single base pair substitutions, the predominant mutation type (75%) were predicted missense variants caused by G > T transversions (Additional file 1: Table S2). However, all but one (CATSPERB:c.2309C > A; p.Pro770His) of the G > T transversion could not be validated by Sanger sequencing, pointing to technical artifacts; in almost all of them the read depth was low (10-19×). Corresponding to the NGS results, seven somatic variants affecting seven genes were confirmed (both directions); all of them are predicted to be missense variants and all showed a low fraction of the variant allele indicating tissue heterogeneity / mosaicism in the tumour sample (Table 2, Fig. 1). All affected genes are expressed in colorectal tissue.

Table 2 Details of seven predicted missense variants which were validated by Sanger sequencing

Each variant was present in a single polyp only. Four mutations occurred in the same polyp (T71), the other three were identified each in another polyp (T69, T70, and T72). No known cancer-associated gene was among the seven candidates. Data mining according to gene function and pathways demonstrated that the majority is involved in processes relevant for tumorigenesis such as cell differentiation and proliferation (Table 2). Most of the variants are neither reported as somatic, tumour-related variant in TCGA or COSMIC, nor as germline variant in large control cohorts (ExAC database) so far. Just the variants in DNAI1 and CALD1 have been reported in ExAC very rarely (heterozygote frequency < 0.01%); the DNAI1 variant was found in a hepatocellular carcinoma (COSMIC) [40]. The ABI3BP and CATSPERB variants are predicted to be damaging by all three in-silico tools (Table 2).


In recent years, a characteristic molecular pattern was identified in both sporadic and SPS-associated serrated lesions, namely specific activating somatic mutations in two members of the mitogen-activated protein kinase (MAPK) cascade (BRAF and KRAS), a CpG Island Methylator Phenotype (CIMP), microsatellite instability (MSI) (reviewed in [7]), and a unique expression signature [15]. The frequency of specific alterations strongly depends on the number, location (right versus left colon) and histologic subtype of polyps examined. In addition, truncating mutations of the tumor suppressor gene RNF43, a transmembrane E3 ubiquitin protein ligase which acts as a Wnt signaling inhibitor by targeting the Frizzled receptor for degradation, was identified in sporadic SSP/A and TSA and as germline mutation in a small subset of SPS families [18, 19, 41].

Our results of exome sequencing, targeted genetic profiling and microsatellite analysis are consistent with published literature: the BRAF missense mutation V600E was found in 67–88% of HPs and 61–83% of SSP/A; the most common KRAS missense mutations (G12D, G12 V, G13D) in 6–17% of HP and in 7–25% of SSP/A. A CIMP, which might result in MLH1 or MGMT (O6-methylguanine-DNA methyltransferase) silencing, was observed in 41–73% of HPs and 44–77% of SSP/A, respectively [30]. In contrast, MSI was rarely examined in serrated polyps and then usually not found in HPs and SSP/A except for a Japanese study where it was identified in 36% of SSP/A [42, 43]. However, in serrated adenocarcinomas, the BRAF mutation, MSI and CIMP are all detected in more than 80% of samples [30] suggesting that high levels of MSI are more relevant in advanced stages of serrated tumorigenesis. Further genetic alterations which have been described in more advanced tumours include targets of MSI (BAX, TGFßR2, IGFR2, p53, often present late during adenoma/carcinoma progression), CDKN2A/p16, MINT genes, 3p, and 15q (CRAC1) [10, 33, 44]. Consistently, it was shown that the number of methylated genes increased significantly in the order of aberrant crypt foci (ACF) to SSA/P to serrated cancer [36].

However, beside of these few marker alterations, the molecular steps and predisposing genetic factors of sporadic and SPS-associated serrated lesions are largely unknown yet. To describe and characterise the mutation spectrum and pattern of serrated lesions, we have chosen to select a number of representative polyps present in a single 41 years old female patient with clinically confirmed SPS. This approach using a single patient is not sufficient to identify potential predisposing genetic factors of SPS (germline variants), however, it offers the opportunity to investigate a large number of polyps of different stages and colonic regions with the same constitutional genetic background. The recent identification of an SSP/A specific expression profile indicates that the underlying mechanisms are operating in both sporadic and syndromic SSP/As, so that it is likely that the results of SPS-associated polyps might be also representative for the sporadic counterpart [15].

The synchronus or metachronous occurrence of dozens of serrated polyps strongly argues for an underlying genetic basis, although it remains unclear so far, whether the predisposing genetic factors mainly act in a monogenic fashion, or contribute as low or moderately penetrant variants to a more complex, oligo/polygenic trait. Thus, it can be hypothesized that hyperplastic and serrated polyps originating from an SPS patient, share genetic features or involved pathways of tumorigenesis to a certain degree.

The aim of the present study was the systematic and exome-wide identification of point mutations (single base pair substitutions, small indels) in protein-coding genes as potential novel oncogenic drivers of serrated tumorigenesis rather than a comprehensive molecular approach covering all levels of possible alterations including methylation pattern and gross genomic alterations. By targeted and exome sequencing of eleven early hyperplastic / serrated polyps from all parts of the colon, we found a somatic mutational profile characteristic for SPS patients with the specific oncogenic BRAF and KRAS mutations present in more than half of the examined polyps which is comparable to other studies [45, 46]. It is known that the frequency of detected BRAF and KRAS mutations observed in SPS patients depends on the number and the location of polyps examined [45, 47].

Usually, SSP/As are larger than HPs and localized predominantly in the proximal colon [48, 49] which is similar to our findings. Interestingly, the specific BRAF or KRAS mutations were identified in almost all polyps from the proximal but no polyp from the distal colon although the relevant genomic regions had a minimum of 30× and 50× coverage for BRAF and KRAS, respectively, in the DNA derived from mutation negative polyps.

The association of proximal polyps with a BRAF mutation and MLH1 methylation was already reported by others who suggested that this relation increases the risk of progression to malignancy [50].

All BRAF and KRAS mutations identified in the exome sequencing data of the polyps were confirmed by targeted molecular analysis which demonstrates the ability of our approach to identify relevant variants with a fraction of variant reads 3 5% in tumour DNA.

Subsequently, all variants called in the exome sequencing data were filtered in multiple stringent steps to select for rare, non-polymorphic variants. No relevant variants were detected in 123 potential candidate genes for serrated tumorigenesis published in recent years including RNF43, several targets of CIMP, and MSI [19, 33]. Finally, 25 somatic variants, all of which were single base pair substitutions, were identified. Of these variants, 18 could not be validated by Sanger sequencing although the ratio of variant alleles in exome data was clearly over the detection threshold of Sanger sequencing. Interestingly, all of the non-validated variants were G > T transversions and in almost all of them the read depth was low, pointing to technical artifacts which is in line with previous observations in HiSeq exome data [51, 52].

Altogether, seven rare variants could be validated by Sanger sequencing in the examined polyps. They affect seven different genes and all are predicted to result in missense changes. Of note, five of the genes are related to functions relevant to tumour development such as proliferation, cell differentiation, or cell invasion. However, all variants occurred each in a single polyp only and none of the involved genes seem to belong to a shared pathway, i.e. there was no evidence for recurrently mutated genes or functionally related gene groups. Among the promising candidate genes, only the two variants of ABI3BP and CATSPERB are predicted to be deleterious. Thus, proof of the causal relevance of these variants remains challenging.

ABI3BP (also known as TARSH, NESH, NESHBP) is a positive regulator of cell-substrate adhesion and involved in an extracellular matrix organization. It plays an important role in proliferation and replicative senescence and may serve as a trigger of tumour development [53]. ABI3BP prevents chromosome instability in a p53-independent manner, resulting in suppression of carcinogenesis [54]. Down-regulation of ABI3BP expression has been found in thyroid tumorigenesis [55] and lung cancer cell lines [56], and thus, deleterious mutations of ABI3AP might contribute the tumour development.

CATSPERB is involved in cell differentiation. The CATSPER gene family is known to be relevant for spermatogenesis [57] and CATSPER SNPs have been reported to be associated with bone mineral density in premenopausal women [58], but there is no study which relates CATSPER variants to cancer so far.

Our results are consistent with the present knowledge regarding the pathology of early serrated polyps. There are only few published studies that aimed to identify genetic changes in early serrated polyps. Apart from RNF43, no convincing novel alterations were found across studies so far. In recent years, a number of candidate driver genes for SSP/A development and progression have been identified by expression analysis, however, the replication of these findings and the evaluation of the clinical and biological relevance has yet to be demonstrated. Our findings do not indicate that point mutations of these candidate genes are the underlying mechanism leading to aberrant expression. Recently, Sakai et al. performed a targeted mutation screening of 126 candidate cancer driver genes in 25 SSA/Ps and found somatic mutations in several cancer pathways [14]. However, although the SSA/Ps were larger (27 mm on average) compared to those examined in this study, they could not identify any mutation in 2, no mutation besides of BRAF in 9, and only one additional mutation besides of BRAF in 7 SSA/Ps. Therefore, the negative results of the present study compared to the study of Sakai et al. might indicate considerable genetic heterogeneity of somatic events in serrated polyps. In addition, the difference might be caused by the smaller size of polyps in the present study since less somatic events are expected to be present in early stages of tumorigenesis.

We cannot rule out that variants might have been overlooked due to a very low allelic fraction (tissue heterogeneity) or incomplete coverage. Moreover, we did not cover the whole spectrum of potential genetic and epigenetic alterations in serrated lesions since the analysis of variants outside the coding part of the genome (e.g. promoters, deep intronic regions, intergenic regions), aberrant methylation pattern, and structural rearrangements was beyond the scope of this study.


To the best of our knowledge, this is the first study which performed exome sequencing in a number of serrated polyps from a single patient to identify potential novel drivers of serrated tumorigenesis. Our data indicate that somatic mutations beyond the well-known driver mutations in the established (BRAF/KRAS) and novel (RNF43) genes seem to be rare events in early BRAF/KRAS-related serrated lesions of SPS patients. No frequently affected genes and no enrichment of specific pathways have been observed. In addition, none of the current candidate genes known to be affected by epigenetic and/or expression changes harboured a relevant point mutation. Thus, other alterations such as specific epigenetic changes, in particular hypermethylation of tumour suppressor genes, or mutations in regulatory regions, which result in aberrant expression, might be the major driving force of tumour progression in sporadic and syndromic serrated polyps as proposed by other investigators [15, 32, 45]. This would be consistent with the recent observation that aberrant crypt foci (ACF) as the earliest precancerous lesions of the serrated pathway, arise by a BRAF mutation and methylation of a few genes and develop into SSP/As through accumulated methylation of a limited number of additional genes [15, 36].


1° :

The first degree


Aberrant crypt foci


Burrow-wheeler aligner


CpG island methylator phenotype


Chromosomal instability


Colorectal cancer


Deoxyribonucleic acid


Genome Analysis Toolkit


Hereditary non polyposis colorectal cancer


Hyperplastic polyposis syndrome


Hyperplastic polyps


Minor allele frequency


MUTYH-associated polyposis


Mitogen-activated protein kinase




Microsatellite instability


Microsatellite stability


National Cancer Institute


Next generation sequencing


Sequence alignment/map tools


Single nucleotide variant


Serrated polyposis syndrome


Sessile serrated polyps/adenomas


Traditional serrated adenomas


Whole exome sequencing


World Health Organization


  1. 1.

    Snover DC, Batts KP. Serrated colorectal Neoplasia. Surg Pathol Clin. 2010;3(2):207–40.

    Article  PubMed  Google Scholar 

  2. 2.

    Delker DA, McGettigan BM, Kanth P, Pop S, Neklason DW, Bronner MP, Burt RW, Hagedorn CH. RNA sequencing of sessile serrated colon polyps identifies differentially expressed genes and immunohistochemical markers. PLoS One. 2014;9(2):e88367.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Kim MJ, Lee EJ, Suh JP, Chun SM, Jang SJ, Kim DS, Lee DH, Lee SH, Youk EG. Traditional serrated adenoma of the colorectum: clinicopathologic implications and endoscopic findings of the precursor lesions. Am J Clin Pathol. 2013;140(6):898–911.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Jass JR, Whitehall VL, Young J, Leggett BA. Emerging concepts in colorectal neoplasia. Gastroenterology. 2002;123(3):862–76.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Leonard DF, Dozois EJ, Smyrk TC, Suwanthanma W, Baron TH Sr, Cima RR, Larson DW. Endoscopic and surgical management of serrated colonic polyps. Br J Surg. 2011;98(12):1685–94.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Fearon ER, Vogelstein B. A genetic model for colorectal tumorigenesis. Cell. 1990;61(5):759–67.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Yamane L, Scapulatempo-Neto C, Reis RM, Guimaraes DP. Serrated pathway in colorectal carcinogenesis. World J Gastroenterol. 2014;20(10):2634–40.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Leggett B, Whitehall V. Role of the serrated pathway in colorectal cancer pathogenesis. Gastroenterology. 2010;138(6):2088–100.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Jass JR, Young J, Leggett BA. Hyperplastic polyps and DNA microsatellite unstable cancers of the colorectum. Histopathology. 2000;37(4):295–301.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Kim H, An JY, Noh SH, Shin SK, Lee YC, Kim H. High microsatellite instability predicts good prognosis in intestinal-type gastric cancers. J Gastroenterol Hepatol. 2011;26(3):585–92.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Rosty C, Young JP, Walsh MD, Clendenning M, Walters RJ, Pearson S, Pavluk E, Nagler B, Pakenas D, Jass JR, et al. Colorectal carcinomas with KRAS mutation are associated with distinctive morphological and molecular features. Mod Pathol. 2013;26(6):825–34.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Gonzalo DH, Lai KK, Shadrach B, Goldblum JR, Bennett AE, Downs-Kelly E, Liu X, Henricks W, Patil DT, Carver P, et al. Gene expression profiling of serrated polyps identifies annexin A10 as a marker of a sessile serrated adenoma/polyp. J Pathol. 2013;230(4):420–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Caruso M, Fung KY, Moore J, Brierley GV, Cosgrove LJ, Thomas M, Cheetham G, Brook E, Fraser LM, Tin T, et al. Claudin-1 expression is elevated in colorectal cancer precursor lesions harboring the BRAF V600E mutation. Transl Oncol. 2014;7(4):456–63.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Sakai E, Fukuyo M, Ohata K, Matsusaka K, Doi N, Mano Y, Takane K, Abe H, Yagi K, Matsuhashi N, et al. Genetic and epigenetic aberrations occurring in colorectal tumors associated with serrated pathway. Int J Cancer. 2016;138(7):1634–44.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Kanth P, Bronner MP, Boucher KM, Burt RW, Neklason DW, Hagedorn CH, Delker DA. Gene signature in sessile serrated polyps identifies colon cancer subtype. Cancer Prev Res (Phila). 2016;9(6):456–65.

    CAS  Article  Google Scholar 

  16. 16.

    IJspeert JE, Rana SA, Atkinson NS, van Herwaarden YJ, Bastiaansen BA, van Leerdam ME, Sanduleanu S, Bisseling TM, Spaander MC, Clark SK, et al. Clinical risk factors of colorectal cancer in patients with serrated polyposis syndrome: a multicentre cohort analysis. Gut. 2017;66(2):278–84.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Guarinos C, Sanchez-Fortun C, Rodriguez-Soler M, Alenda C, Paya A, Jover R. Serrated polyposis syndrome: molecular, pathological and clinical aspects. World J Gastroenterol. 2012;18(20):2452–61.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Gala MK, Mizukami Y, Le LP, Moriichi K, Austin T, Yamamoto M, Lauwers GY, Bardeesy N, Chung DC. Germline mutations in oncogene-induced senescence pathways are associated with multiple sessile serrated adenomas. Gastroenterology. 2014;146(2):520–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Yan HH, Lai JC, Ho SL, Leung WK, Law WL, Lee JF, Chan AK, Tsui WY, Chan AS, Lee BC, et al. RNF43 germline and somatic mutation in serrated neoplasia pathway and its association with BRAF mutation. Gut. 2017;66(9):1645–56.

    Article  PubMed  Google Scholar 

  20. 20.

    Aretz S, Uhlhaas S, Goergens H, Siberg K, Vogel M, Pagenstecher C, Mangold E, Caspari R, Propping P, Friedl W. MUTYH-associated polyposis: 70 of 71 patients with biallelic mutations present with an attenuated or atypical phenotype. Int J Cancer. 2006;119(4):807–14.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Pitera R, Pitera JE, Mufti GJ, Salisbury JR. Modification of standard proteinase K/phenol method for extraction of DNA from small tumour biopsies. Pathol Res Pract. 1993;189(8):882–7.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Ney JT, Froehner S, Roesler A, Buettner R, Merkelbach-Bruse S. High-resolution melting analysis as a sensitive prescreening diagnostic tool to detect KRAS, BRAF, PIK3CA, and AKT1 mutations in formalin-fixed, paraffin-embedded tissues. Arch Pathol Lab Med. 2012;136(9):983–92.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Boland CR, Thibodeau SN, Hamilton SR, Sidransky D, Eshleman JR, Burt RW, Meltzer SJ, Rodriguez-Bigas MA, Fodde R, Ranzani GN, et al. A National Cancer Institute workshop on microsatellite instability for cancer detection and familial predisposition: development of international criteria for the determination of microsatellite instability in colorectal cancer. Cancer Res. 1998;58(22):5248–57.

    CAS  PubMed  Google Scholar 

  25. 25.

    Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Albers CA, Lunter G, MacArthur DG, McVean G, Ouwehand WH, Durbin R. Dindel: accurate indel calls from short-read data. Genome Res. 2011;21(6):961–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Yeo G, Burge CB. Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol. 2004;11(2–3):377–94.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Kambara T, Simms LA, Whitehall VL, Spring KJ, Wynter CV, Walsh MD, Barker MA, Arnold S, McGivern A, Matsubara N, et al. BRAF mutation is associated with DNA methylation in serrated polyps and cancers of the colorectum. Gut. 2004;53(8):1137–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    O’Brien MJ, Yang S, Mack C, Xu H, Huang CS, Mulcahy E, Amorosino M, Farraye FA. Comparison of microsatellite instability, CpG island methylation phenotype, BRAF and KRAS status in serrated polyps and traditional adenomas indicates separate pathways to distinct colorectal carcinoma end points. Am J Surg Pathol. 2006;30(12):1491–501.

    Article  PubMed  Google Scholar 

  31. 31.

    Sjoblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N, et al. The consensus coding sequences of human breast and colorectal cancers. Science. 2006;314(5797):268–74.

    Article  PubMed  Google Scholar 

  32. 32.

    Esteller M. Epigenetics in cancer. N Engl J Med. 2008;358(11):1148–59.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Grady WM, Carethers JM. Genomic and epigenetic instability in colorectal cancer pathogenesis. Gastroenterology. 2008;135(4):1079–99.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Kaji E, Uraoka T, Kato J, Hiraoka S, Suzuki H, Akita M, Saito S, Tanaka T, Ohara N, Yamamoto K. Externalization of saw-tooth architecture in small serrated polyps implies the presence of methylation of IGFBP7. Dig Dis Sci. 2012;57(5):1261–70.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Beggs AD, Jones A, Shepherd N, Arnaout A, Finlayson C, Abulafi AM, Morton DG, Matthews GM, Hodgson SV, Tomlinson IP. Loss of expression and promoter methylation of SLIT2 are associated with sessile serrated adenoma formation. PLoS Genet. 2013;9(5):e1003488.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Inoue M, Takahashi S, Soeda H, Shimodaira H, Watanabe M, Miura K, Sasaki I, Kato S, Ishioka C. Gene-expression profiles correlate with the efficacy of anti-EGFR therapy and chemotherapy for colorectal cancer. Int J Clin Oncol. 2015;20(6):1147–55.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46(3):310–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB et al: Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536(7616):285-291.

  39. 39.

    Petrovski S, Wang Q, Heinzen EL, Allen AS, Goldstein DB. Genic intolerance to functional variation and the interpretation of personal genomes. PLoS Genet. 2013;9(8):e1003709.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A, et al. COSMIC: mining complete cancer genomes in the catalogue of somatic mutations in cancer. Nucleic Acids Res. 2011;39(Database issue):D945–50.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Bond CE, McKeone DM, Kalimutho M, Bettington ML, Pearson SA, Dumenil TD, Wockner LF, Burge M, Leggett BA, Whitehall VL. RNF43 and ZNRF3 are commonly altered in serrated pathway colorectal tumorigenesis. Oncotarget. 2016;7(43):70589–600.

    PubMed  PubMed Central  Google Scholar 

  42. 42.

    Guarinos C, Sanchez-Fortun C, Rodriguez-Soler M, Perez-Carbonell L, Egoavil C, Juarez M, Serradesanferm A, Bujanda L, Fernandez-Banares F, Cubiella J, et al. Clinical subtypes and molecular characteristics of serrated polyposis syndrome. Clin Gastroenterol Hepatol. 2013;11(6):705–11. quiz e746

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Konishi K, Yamochi T, Makino R, Kaneko K, Yamamoto T, Nozawa H, Katagiri A, Ito H, Nakayama K, Ota H, et al. Molecular differences between sporadic serrated and conventional colorectal adenomas. Clin Cancer Res. 2004;10(9):3082–90.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Fujita K, Yamamoto H, Matsumoto T, Hirahashi M, Gushima M, Kishimoto J, Nishiyama K, Taguchi T, Yao T, Oda Y. Sessile serrated adenoma with early neoplastic progression: a clinicopathologic and molecular study. Am J Surg Pathol. 2011;35(2):295–304.

    Article  PubMed  Google Scholar 

  45. 45.

    Kim KM, Lee EJ, Ha S, Kang SY, Jang KT, Park CK, Kim JY, Kim YH, Chang DK, Odze RD. Molecular features of colorectal hyperplastic polyps and sessile serrated adenoma/polyps from Korea. Am J Surg Pathol. 2011;35(9):1274–86.

    Article  PubMed  Google Scholar 

  46. 46.

    Pai RK, Hartman DJ, Gonzalo DH, Lai KK, Downs-Kelly E, Goldblum JR, Liu X, Patil DT, Bennett AE, Plesec TP, et al. Serrated lesions of the appendix frequently harbor KRAS mutations and not BRAF mutations indicating a distinctly different serrated neoplastic pathway in the appendix. Hum Pathol. 2014;45(2):227–35.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Kwon HJ, Cho NY, Chang MS, Kim YS, Kang GH. Intermediate serrated polyp as an intermediate lesion of hyperplastic polyp and sessile serrated polyp/adenoma in terms of morphological and molecular features. Hum Pathol. 2014;45(8):1759–65.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Schramm C, Kaiser M, Drebber U, Gruenewald I, Franklin J, Kuetting F, Bowe A, Hoffmann V, Gatzke S, Toex U, et al. Factors associated with reclassification of hyperplastic polyps after pathological reassessment from screening and surveillance colonoscopies. Int J Color Dis. 2016;31(2):319–25.

    Article  Google Scholar 

  49. 49.

    Okamoto K, Kitamura S, Kimura T, Nakagawa T, Sogabe M, Miyamoto H, Muguruma N, Takayama T. Clinicopathological characteristics of serrated polyps as precursors to colorectal cancer: current status and management. J Gastroenterol Hepatol. 2017;32(2):358–67.

    Article  PubMed  Google Scholar 

  50. 50.

    Fernando WC, Miranda MS, Worthley DL, Togashi K, Watters DJ, Leggett BA, Spring KJ. The CIMP phenotype in BRAF mutant serrated polyps from a prospective colonoscopy patient cohort. Gastroenterol Res Pract. 2014;2014:374926.

    Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Do H, Dobrovic A. Dramatic reduction of sequence artefacts from DNA isolated from formalin-fixed cancer biopsies by treatment with uracil- DNA glycosylase. Oncotarget. 2012;3(5):546–58.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Costello M, Pugh TJ, Fennell TJ, Stewart C, Lichtenstein L, Meldrim JC, Fostel JL, Friedrich DC, Perrin D, Dionne D, et al. Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation. Nucleic Acids Res. 2013;41(6):e67.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Uekawa N, Terauchi K, Nishikimi A, Shimada J, Maruyama M. Expression of TARSH gene in MEFs senescence and its potential implication in human lung cancer. Biochem Biophys Res Commun. 2005;329(3):1031–8.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Wakoh T, Uekawa N, Terauchi K, Sugimoto M, Ishigami A, Shimada J, Maruyama M. Implication of p53-dependent cellular senescence related gene, TARSH in tumor suppression. Biochem Biophys Res Commun. 2009;380(4):807–12.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Latini FR, Hemerly JP, Oler G, Riggins GJ, Cerutti JM. Re-expression of ABI3-binding protein suppresses thyroid tumor growth by promoting senescence and inhibiting invasion. Endocr Relat Cancer. 2008;15(3):787–99.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Terauchi K, Shimada J, Uekawa N, Yaoi T, Maruyama M, Fushiki S. Cancer-associated loss of TARSH gene expression in human primary lung cancer. J Cancer Res Clin Oncol. 2006;132(1):28–34.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Liu J, Xia J, Cho KH, Clapham DE, Ren D. CatSperbeta, a novel transmembrane protein in the CatSper channel complex. J Biol Chem. 2007;282(26):18945–52.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Koller DL, Ichikawa S, Lai D, Padgett LR, Doheny KF, Pugh E, Paschall J, Hui SL, Edenberg HJ, Xuei X, et al. Genome-wide association study of bone mineral density in premenopausal European-American women and replication in African-American women. J Clin Endocrinol Metab. 2010;95(4):1802–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank the patient for participating in the study.


This work was supported by the German Cancer Aid (Deutsche Krebshilfe e.V. Bonn, Grant number 108421). This funding source had no involvement in the study design; the collection, analysis, or interpretation of data; the writing of the report; or the decision to submit the manuscript for publication.

Availability of data and materials

Not applicable

Databases / web-resources

CADD (Combined Annotation Dependent Depletion):

COSMIC (Catalogue of somatic mutations in cancer):

ExAC Browser (Exome Aggregation Consortium):

gnomAD Browser (Genome Aggregation Database):

IGV (Integrative Genomics Viewer):

Primer3 v.0.4.0:

RVIS (Residual Variation Intolerance Score):




SIFT tool:

TCGA (The Cancer Genome Atlas):

Author information




SH was involved in protocol optimization, data analysis, and the manuscript preparation. SA designed and supervised the study, discussed results, and wrote and revised the manuscript. IS, SP, JA, were involved in providing and checking clinical data, data analysis, and protocol optimization. JK, MK, and MO performed histological examinations and somatic targeted mutation analysis. DD, HT, PT were involved in NGS data analysis. PK, RH, CS, GK, MO, JN, RB contributed to sample and clinical data collection and manuscript revision. PH was involved in initial study design, supervison, and manuscript revision. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Stefan Aretz.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the local ethics review board (Medical Faculty of the University of Bonn ethics review board no. 208/12), and a written informed consent was obtained from the patient.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Established causative or candidate genes for (hereditary) colorectal tumours and published genes related to serrated polyps (n = 74). Table S2. Primers used for validating 25 variants by Sanger sequencing. (PDF 239 kb)

Additional file 2: Figure S1.

Sanger sequencing confirmed the KRAS G12D mutation in tumour T69 and the BRAF V600E mutation in tumour T71 and T76. (PDF 72 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Horpaopan, S., Kirfel, J., Peters, S. et al. Exome sequencing characterizes the somatic mutation spectrum of early serrated lesions in a patient with serrated polyposis syndrome (SPS). Hered Cancer Clin Pract 15, 22 (2017).

Download citation


  • Familial colorectal cancer (CRC)
  • Exome sequencing
  • Hyperplastic polyposis syndrome (HPS)
  • Serrated polyposis syndrome (SPS)
  • Serrated pathway
  • Somatic mutations