- Open Access
Somatic targeted mutation profiling of colorectal cancer precursor lesions
BMC Medical Genomics volume 15, Article number: 143 (2022)
Most colorectal cancers (CRC) arise from precursor lesions. This study aimed to characterize the mutation profile of colorectal cancer precursor lesions in a Brazilian population.
In total, 90 formalin-fixed paraffin-embedded colorectal precursor lesions, including 67 adenomas, 7 sessile serrated lesions, and 16 hyperplastic polyps, were analyzed by next-generation sequencing using a panel of 50 oncogenes and tumor suppressor genes. The genetic ancestry of the patients was estimated.
Somatic driver mutations were identified in 66.7% of cases, including alterations in APC (32.2%), TP53 (20.0%), KRAS (18.9%), BRAF (13.3%) and EGFR (7.8%). Adenomas displayed a higher number of mutations, mainly in APC, compared to serrated polyps (73.1% vs. 47.8%, p = 0.026). Advanced adenomas had a significantly higher frequency of mutation in KRAS and a high overall mutation rate than early adenomas (92.9% vs. 59%, p = 0.006). A high degree of ancestry admixture was observed in the population studied, with a predominance of European components (mean of 73%) followed by African (mean of 11.3%). No association between genetic ancestry and type of lesions was found. The mutation profile of Brazilian colorectal precursor lesions exhibits alteration in APC, KRAS, TP53, and BRAF at different frequencies according to lesion type.
These results bestow the knowledge of CRC's biologic history and support the potential of these biomarkers for precursor lesions detection in CRC screening of the Brazilian population.
Colorectal cancer (CRC) is the third most incident cancer worldwide, resulting in 915,880 deaths in 2020 . In Brazil, CRC ranks second in incidence in men and women , and a continuous increase in both incidence and mortality is expected in the coming years [3,4,5]. The development of colorectal cancer is a multi-stage evolution process that occurs through a progressive accumulation of molecular alterations in the colon epithelium cells, which can be transformed into pre-malignant lesions and cancer [6, 7].
Although viewed as a single disease, from the molecular and morphological point of view, CRC is a heterogeneous disease that is believed to arise mainly from two different types of precursor lesions: adenoma and serrated polyps [7, 8]. In the classic sequence of adenoma to carcinoma progression model, the development of colorectal cancer originates from aberrant crypts, progressing to early adenoma, advanced adenoma, and subsequently carcinoma . Serrated polyps have recently been recognized as important precursor lesions and account for approximately 15%-30% of all cases of CRC [8,9,10,11]. According to the updated World Health Organization (WHO) classification, serrated polyps include hyperplastic polyp (HP), sessile serrated lesion (SSL), and traditional serrated adenoma (TSA) . Among these lesions, the sessile serrated lesion and traditional serrated adenoma are more likely to evolve into cancer. Despite being the most frequent type of lesion, hyperplastic polyps are considered to have no malignant potential [13, 14]. Moreover, patients with proximal serrated polyps, particularly those larger than 10 mm, are associated with an increased risk of developing CRC .
The molecular mechanisms underlying the progression through the canonical pathway frequently comprise somatic mutations in oncogenes such as KRAS and tumor suppressor genes, such as APC, TP53, and SMAD4 . In addition, mutations in BRAF have an important play in the serrated pathway [7, 15]. Recently, our group performed a mutational portrait of Brazilian CRC patients and observed a similar molecular portrait than reported worldwide .
Screening for CRC can reduce incidence and mortality by detecting and removing precursor lesions . This strategy is possible due to the long period of progression from a precursor lesion (adenoma) to cancer, which ranges from 7 to 10 years . Colonoscopy is the most suitable and reliable diagnostic tool for CRC screening, yet it is not feasible for large-scale due to its risks and high cost. In organized population-based screening, fecal occult blood tests are preferred , with the FIT (fecal immunochemical test) widely used. Still, FIT is not perfect; its sensitivity for CRC ranges from 60 to 80% and only 20–30% for advanced adenoma . To overcome these limitations, DNA-based analysis of body fluids–liquid biopsy–such as blood or feces, can increase the accuracy of FIT. Two commercialized assays are FDA approved for CRC screening: Epi ProColon 2.0, analyzing SEPT9 DNA methylation in blood; Cologuard test®, stool-based that analyses NDRG4 and BMP3 methylation and KRAS mutation [20, 21].
Therefore, identifying genetic alterations in precursor lesions can lead to molecular-based strategies, improving the sensitivity, specificity, and impact of colorectal screening and surveillance programs. Nevertheless, few studies evaluated the mutation status in South America CRC precursor lesions [22, 23]. Results of the first two years (first round) of our colorectal cancer screening program (from Barretos Cancer Hospital program) in Brazil showed a successful implementation with a high test return participation rate, colonoscopy completion, and detection lesion rates. In addition, participant’s blood and FIT (fecal immunochemical test) tests has been stored in the Barretos Cancer Hospital biobank to allow future biomarker studies and consequently improve lesions detection rates .
Herein, we performed the mutation analysis by next-generation sequencing (NGS) of 50 oncogenes and tumor suppressor genes in colorectal cancer precursor lesions and also evaluated the genetic ancestry composition of the Brazilian samples included in the study. Somatic nucleotide variants were identified in all types of precursor lesions, most of them more prevalent in the adenoma group. Moreover, the spectrum of mutated genes was different between adenomas and serrated polyps. These results extend our knowledge of the molecular biological features of precursor lesions and the natural history of colorectal cancer.
This retrospective study analyzed a total of 90 formalin-fixed paraffin-embedded (FFPE) samples removed from 87 patients during diagnostic colonoscopy from 2014 to 2016 at Barretos Cancer Hospital . Samples included 67 adenomas (39 early and 28 advanced adenomas), 7 sessile serrated lesions and 16 hyperplastic polyps. The 87 patients were between 49 and 88 years of age. Subjects with a personal history of familial adenomatous polyposis (FAP) or Lynch syndrome were excluded [26, 27].
The Institutional Research Board of the Barretos Cancer Hospital approved the study on Feb 4, 2016 (number ID: 1074/2016). Informed consent was waived due to the retrospective nature of this study. The study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki.
Lesions were classified according to Paris classification  and histological analysis using WHO criteria . Villous structures in > 25% of adenoma was required for tubulovillous adenoma. If > 75% of the adenoma has a villous architecture, it was diagnosed as villous adenoma. Adenomas were classified as advanced if > 1 cm in diameter or presented more than 25% of villous structures in histology or high-grade intraepithelial neoplasia. Table 1 summarizes the histopathological data of the 90 samples.
DNA was isolated from FFPE tissue as previously reported . Hematoxylin and eosin stained slides were reviewed by a pathologist and the contents of precursor lesions samples were more than 50%. The lesion area was delimited and macrodissected from six 10-μm-thick unstained tissue sections. Xylene and ethanol (100%, 70%, 50%) were used for paraffin removal.
FFPE genomic DNA was extracted using the QiaAmp DNA micro Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. DNA concentration was determined using Qubit™ dsDNA HS assay kit (Thermo Fisher Scientific, Eugene, Oregon, USA) on the Qubit 2.0 Fluorometer (Thermo Fisher Scientific).
Library preparation and sequencing
Next-generation sequencing library preparation was conducted to amplify 10 ng of genomic DNA using AmpliSeq™ Cancer HotSpot Panel v2 panel kit (Thermo Fisher Scientific) and was performed with Ion Ampliseq™ Kit for Chef DL8 in the Ion Chef instrument. This panel includes primer for regions of the following 50 cancer driver genes—SMARCB1, RB1, TP53, ERBB4, FBXW7, BRAF, KIT, GNAS, HRAS, EGFR, PDGFRA, PIK3CA, CDKN2A, ERBB2, ABL1, JAK2, KRAS, NRAS, NOTCH1, ATM, FGFR1, STK11, PTPN11, APC, SMAD4, PTEN, SMO, CTNNB1, RET, IDH2, SRC, EZH2, VHL, MPL, NPM1, FLT3, FGFR3, CDH1, KDR, HNF1A, MLH1, ALK, IDH1, GNAQ, AKT1, JAK3, FGFR2, GNA11, MET, CSF1R. Pooled libraries were submitted to emulsion PCR, enrichment of beads containing the template and chip loading in the Ion Chef instrument using Ion PGM™ Hi-Q™ View Chef Kit according to the manufacturer’s instructions (Thermo Fisher Scientific, Waltham, MA). The final library was sequenced in an Ion 318 v2 chip on Ion Torrent PGM using Ion PGM™ Hi-Q™ View supplies.
Sequencing data were processed in the Ion PGM™ Torrent Server and generated reads were aligned to the reference genome (hg19) using TMAP (Torrent Mapping Alignment Program) in the Torrent Suite™ Software (ThermoFisher). Variant calling and annotation were performed using the Ion Reporter™ Software (version 5.10).
Only variants with sequencing depth of at least 200 × and variant allele frequency (VAF) > 10% were retained. Intronic and synonymous variants were filtered out, as well variants with frequency higher than 1% in population database (ExAC), available in The Cancer Genome Interpreter tool (CGI). This platform was also employed to verify the status of driver variants and only known variants or predicted as driver in colorectal cancer-related genes were retained. In addition, the remaining variants with frequency higher than 1% in the Brazilian genomic variants database (ABraOM) were excluded from further analysis.
Samples harboring the BRAF V600E mutation (n = 3) were selected for the qualitative validation of the variant with the real-time PCR assay cobas 4800 BRAF V600 Mutation Test (Roche Molecular Diagnostics). The test was performed using cobas 4800 System according to the manufacturer’s instructions.
Genetic ancestry analysis
DNA isolation from peripheral blood samples was performed using the QIAmp DNA Blood Mini Kit (QIAGEN, Hilden, Germany) following the manufacturer’s instructions and by Biobank Barretos procedures . DNA concentration was determined using NanoDrop™ Spectrophotometer (Thermo Scientific).
The genetic ancestry of 81 patients with colorectal precursor lesions with available blood was determined by 46 autosomal ancestry informative markers (AIMs), which consist of insertion-deletion polymorphisms (INDELs) as described . A multiplex PCR was followed by a fragment analysis performed by ABI 3500xL Genetic Analyzer (Applied Biosystems). The analysis of genotypes was performed with GeneMapper Software v4.1(Applied Biosystems).
Genetic data of the Human Genome Diversity Project Center d’Etude du Polymorphisme Humain (HGDP-CEPH)  was used for the classification of the proportion of ancestry for each of the main populations: African, European, East Asian and, Native American using the Structure Software v2.3.4 [32, 33].
Descriptive statistics were expressed in number, percentage, mean and standard deviation. The Chi-square test and Fisher's exact test were used to compare the mutation rates between the different pathological and clinical features of colorectal lesions. The association of the genetic ancestry component (AFR-African, EUR-European, EAS-Asian, NAM-Native American) with precursor lesions was performed using the Kruskal–Wallis test. P values were adjusted by Bonferroni correction for multiple comparisons analysis.
All statistical analyses were performed using SPSS software (v.21) and R software (v. 3.6.1).
Summary of identified driver variants
The mean amplicon sequence coverage of 200 × used in AmpliSeq Cancer HotSpot v2 panel was 99.27%. The average coverage per amplicon, number of mapped reads, on target percent and mean depth per sample are shown in Additional file 1: Fig. S1 and Additional file 4: Table S1.
Among the 90 lesions included in this study, 60 (66.7%) showed at least one driver variant. Overall, a total of 124 somatic non-synonymous driver variants were identified in 16 genes. Sequencing of colorectal lesions showed a mean of 1.4 (range of 1–8) driver mutations per lesion. The following mutation frequencies were identified: 33.3% (30/90) of lesions showed no driver mutations, 34.4% (31/90) of lesions carried one mutation per case, 14.4% (13/90) had two, 5.6% (5/90) had three mutations and 11 samples (12.2%) had four or more mutations.
The spectrum of identified mutations included 92 missense, 20 nonsense, 10 frameshift and 2 splice site variants. APC alterations were the most common in our samples (32.2%), followed by TP53 (20.0%), KRAS (18.9%), BRAF (13.3%) and EGFR (7.8%) (Fig. 1, Additional file 2: Fig. S2 and Additional file 5: Table S2 and Additional file 6: Table S3). Mutations in three selected BRAF V600E mutated cases were validated using cobas® 4800 BRAF V600 Mutation Test.
Variants in colorectal lesions exhibit differences based on lesion histological type
We further analyzed the variants according to the lesion histological type and localization (Fig. 1). The frequency of driver mutations in adenomas was slightly higher (73.1% of cases with at least one variant) when compared to serrated polyps (47.8%, p = 0.117) (Table 2). In addition, a mean of 1.6 driver mutations per case were identified in adenomas and 0.7 in serrated polyps (p = 0.01).
Number of variants were significantly higher in advanced adenomas than early adenomas with 2.2 vs.1.2 variants per lesion (p = 0.03), respectively, and were detected in 92.9% and 59.0% of the lesions, respectively (p = 0.006). Slightly differences were also observed between serrated polyps, although not significant, with more variants detected in sessile serrated lesions when compared to hyperplastic polyps (85.7% vs 31.3%, p = 0.081). The mean numbers of mutations were 1.4 in SSL and 0.4 in hyperplastic polyps (p = 0.096) (Additional file 3: Fig. S3).
Differences in frequency of variants in the 16 genes between adenoma, sessile serrated lesions and hyperplastic polyps are shown in Table 3. The differences observed were higher frequency of APC mutations in adenomas when compared to serrated sessile lesion and hyperplastic polyp, while variants in BRAF were more prevalent in sessile serrated lesions (Table 3). When comparing early and advanced adenomas, a higher frequency of mutations in KRAS was observed in the late stage of the lesion (p = 0.001, Fig. 1a, Additional file 7: Table S4). Although we observed a higher frequency in BRAF when comparing hyperplastic polyps and SSL (Fig. 1b, Additional file 7: Table S4), this difference was not significant when adjusted analysis were performed.
Analysis comparing mutational status in adenomas according to the location of lesions revealed that GNAS (p = 0.003) was significantly mutated in rectal adenomas (Additional file 8: Table S5). In serrated polyps, none of the genes were significant altered according to the location (Additional file 9: Table S6).
Molecular pathways associated to colorectal cancer in precursor lesions
Several genes associated with important signaling pathways in cancer were identified as frequently mutated (Fig. 2). The Wnt pathway (APC, SMAD4 and, CTNNB1 genes) was the signaling pathway with the highest frequency of alterations (37.8%). The highest frequency of alteration in this pathway was found among samples of adenoma group when compared to serrated polyps (49.3% vs 4.3%, p < 0.001), with increased frequency of mutations in advanced adenomas, followed by early adenomas, SSLs and hyperplastic polyps (p = 0.001, Fig. 2).
KRAS, BRAF, EGFR, HRAS and FGFR3 genes of MAPK pathway were mutated in 33.3% of the cases. BRAF and KRAS genes were mutually exclusive mutated, except for one advanced adenoma. Mutations were more frequently identified among SSLs, followed by advanced adenomas, hyperplastic polyps, and early adenomas (p < 0.001).
Regarding the PI3K-AKT pathway, mutations in the EGFR, PTEN, PIK3CA, and FGFR3 genes were found in 13.3% of the cases. Mutations in genes involved in the p53 pathway (TP53 and ATM) were found in 20% of the lesions, mostly with mutations in the TP53 gene.
Genetic ancestry of patients with colorectal lesions
The genetic ancestry component was obtained for 80 patients (91.9%) included in the study. The results indicated a high admixture and heterogeneity of the ancestry proportion of Brazilian samples, with the mean of ancestral proportions as follow: 73% (SD = 23.6%) for EUR, 11.3% (SD = 13.2%) for AFR, 9.1% (SD = 18.0%) for EAS and 6.6% (SD = 7.9%) for NAM (Fig. 3). We further compared the average genetic ancestry according to the lesion type, but no significant association was found (Additional file 10: Table S7). According to the mutation status, no differences were observed between the patient’s ancestry (Additional file 11: Table S8).
In the present study, we carried out a molecular profile of 50 cancer-related genes in precursor lesions of CRC. Adenomas exhibited mutations in genes already known to be involved in colorectal carcinogenesis, such as APC, KRAS, TP53, and FBXW7. On the other hand, serrated polyps showed low frequency in APC and TP53 genes and a high frequency of BRAF gene mutations.
Our findings corroborate the molecular differences previously reported in these two major distinct pathways of carcinogenesis . According to the classic adenoma-carcinoma progression model, the progressive accumulation of genetic alterations leads to carcinoma development from the normal mucosa . Recent studies have added complexity to this model, demonstrating the presence of molecular heterogeneity in the early stages of the development of colorectal lesions and mutations in several genes considered drivers for CRC [35,36,37,38]. As expected, we found a higher average of driver mutations in advanced adenomas than in early adenomas. The acquisition rate of mutations is increased in adenomas than normal tissue, and the mutational burden in advanced adenomas has been reported to be similar to cancer tissues, even when only driver mutations are analyzed [35, 36, 39].
We also reported a slightly lower frequency of mutations among serrated polyps when compared to adenomas. Few studies addressed this issue [38, 40]. When comparing only SSL and adenomas, these authors show no difference in the frequency of mutations among these groups, likewise our study. Further, we also observed a slightly higher mutation frequency in SSLs than in hyperplastic polyps, which are lesions with lower malignancy potential. Recently, hyperplastic polyps and serrated sessile lesions were associated with the Consensus Molecular Subtype 1 (CMS1), which often has microsatellite instability (MSI) and hypermutation [41, 42].
Additionally, we found a significant difference in the Wnt, MAPK, PI3K-AKT, and p53 signaling pathways between adenomas and serrated polyps. Alterations in the Wnt pathway are an initial event in the adenoma-carcinoma progression, predominantly due to mutations in the APC (40.3% to 80.0%) followed by the CTNNB1 gene (11.9–20.0%) [43,44,45]. Our study found a lower frequency of APC (41.8%) and CTNNB1 (3.3%) mutations, which can be because we did not analyze the whole coding sequence, but the major hotspot regions of both genes. In the serrated polyps pathway, the Wnt signaling is reported to be less targeted [40, 46], following our findings.
Activation of the MAPK pathway is also observed in CRC, with mutations mainly found in KRAS and BRAF oncogenes [47, 48]. We found 33.3% of our samples harboring mutations in this pathway, with mutations in the KRAS gene slightly more frequent in the adenoma group (22.4%) and BRAF predominantly present in the SSL group. In the adenoma group, the KRAS mutation frequency is within the variation observed in other studies (10.7% to 60.0%) [38, 43, 49, 50]. For the Brazilian population, previous reports on the frequency of KRAS mutation in adenomas have reported a lower frequency than we found (13.6%) . This difference could be explained by the higher sensitivity of NGS used in this study compared to Sanger sequencing to detect low-frequency variants [51, 52]. We also observed a higher frequency of mutations in the KRAS gene in advanced adenomas than in early adenomas, similar to previous studies [38, 53], including reports on the Brazilian population .
A high frequency of mutations in the MAPK pathway genes was observed in serrated polyps, mainly due to the activating BRAF gene mutations in SSLs. This is consistent with a previously reported frequency of BRAF mutations in our population . Mutations in the BRAF gene have been consistently related to SSL with a high frequency of samples (8.7–88%) harboring mutations [23, 40, 45, 50, 54]. Interestingly, the main activating mutation BRAF V600E was found only in serrated polyps, as previously reported [23, 49].
Mutations in TP53 are generally observed during the transition from adenoma to carcinoma [35, 55, 56]. Recent studies reported a lower frequency of TP53 mutations in early or low growth rate adenomas and a higher mutation frequency during the progression of early to advanced adenomas [38, 49]. No significant difference was observed between early and advanced adenomas in our data. Nevertheless, our results agree with the Vogelstein model, where TP53 is associated with the adenoma-carcinoma transition. The frequency of mutations in adenomas was lower when compared to the frequency of mutations in CRC cases previously reported in our population (25.4% in adenomas vs. 56.0% in cancer) . Besides, mutations in TP53 in the serrated polyps were found only in SSLs, which was already described [40, 57].
According to previous reports, genes of the PI3K-AKT pathway were also mutated in our samples [47, 58]. This pathway may present mutations in precursor lesions, focusing on advanced adenomas or traditional sessile adenomas , suggesting a role in the late steps of both adenoma-carcinoma and serrated pathways progression. In agreement with these data, we observed a slightly higher frequency of mutations in genes of this pathway in advanced adenomas than in early adenomas. In the advanced adenomas, we found 3.6% of samples harboring mutations in PIK3CA. Mutations in this gene are found in regions of carcinoma in situ . However, it is not an initial event during clonal diversification in carcinogenesis, as observed in studies of clonal evolution in CRC [35, 56, 60]. Also, PIK3CA mutations are found in cancer-associated adenomas (20.0–30.0%) [36, 43], or lower frequency in advanced adenomas (3.2%) , similar to the frequency observed in our study.
The presence of mutations in the GNAS gene is frequent in CRC  and has been reported in advanced adenomas [35, 54, 62]. Although mutations in GNAS in serrated polyps have already been reported, its frequency is not high and is related to more advanced lesions . Corroborating these data, we identified mutations in this gene only in advanced adenomas samples and absent in serrated polyps.
Previous studies have reported that polyps and colorectal cancer are more frequent among African Americans than non-Hispanic Whites [64,65,66]. In the present study, as expected [30, 67, 68], we observed a high heterogeneity of the ancestry proportions in our study population, yet, we did not find any difference between genetic ancestry and the groups of precursor lesions evaluated. This result could be due to the small number of cases within each group analyzed. Nevertheless, this is the first study to analyze the mutation profile of CRC precursor lesions in this high admixture population, contributing to overcoming disparities and reducing inequalities in the knowledge of colorectal genomic studies [69, 70]. Yet, the identification of somatic alteration in a heterogeneous ancestry population may have distinct medical significance across population groups . Despite major findings, our study’s limitations lie in the absence of paired normal tissue samples compared to the profile mutation of lesion samples. To overcome this issue, the variants identified were filtered in databases, such as ABraOM (Brazilian population) and ExAC (international population). Also, the nature of the targeted sequencing, which does not cover the whole coding sequencing of the cancer genes, could underestimate the mutation frequencies. Finally, the absence of critical CRC-related genes, such as TCF7L2 and FAM123B , could limit our results’ interpretation.
In summary, our study reports the mutation profile of colorectal precursor lesions in Brazilian patients for the first time. We observed the highest mutation frequency in known CRC driver genes, including APC, TP53, KRAS, and BRAF, with differences according to the type of lesion analyzed, with a higher rate of mutations in adenomas. Moreover, a higher number of mutations were found in advanced adenomas compared to early adenomas and in SSL compared to hyperplastic polyps. Collectively, these findings support the potential of these biomarkers for precursor lesions detection in CRC screening of the Brazilian population.
Availability of data and materials
The original contributions presented in the study are publicly available. This data can be found here: http://www.ncbi.nlm.nih.gov/bioproject/772782 (BIOPROJECT accession number: PRJNA772782).
Brazilian genomic variants database
Autosomal ancestry informative markers
The Cancer Genome Interpreter
Consensus Molecular Subtype 1
Exome Aggregation Consortium
Familial adenomatous polyposis
Fecal Immunochemical Test
Human Genome Diversity Project Center d’Etude du Polymorphisme Humain
Polymerase Chain Reaction
Sessile serrated lesion
Traditional serrated adenoma
Torrent Mapping Alignment Program
Variant allele frequency
World Health Organization
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
INCA. Estimativa 2020: Incidência de Câncer no Brasil. Ministério da Saúde: Instituto Nacional de Câncer José de alencar Gomes da Silva (INCA); 2019.
Martin FL, Morais CLM, Sakita JY, Uyemura SA, Kannen V. Age-related and gender-related increases in colorectal cancer mortality rates in Brazil between 1979 and 2015: projections for continuing rises in disease. J Gastrointest Cancer. 2020;3:1147.
Souza DLB, Jerez-Roig J, Cabral FJ, De Lima JRF, Rutalira MK, Costa JAG. Colorectal cancer mortality in Brazil: predictions until the year 2025 and cancer control implications. Dis Colon Rectum. 2014;57:1082–9.
Arnold M, Sierra MS, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global patterns and trends in colorectal cancer incidence and mortality. Gut. 2017;66:683–91.
Vogelstein B, Fearon ER, Hamilton SR, Kern SE, Preisinger AC, Leppert M, et al. Genetic alterations during CRC development. N Engl J Med. 1988;319:525–32.
Kuipers EJ, Grady WM, Lieberman D, Seufferlein T, Sung JJ, Boelens PG, et al. Colorectal cancer. Nat Rev Dis Prim. 2015;1:887.
Li D, Liu L, Fevrier HB, Alexeeff SE, Doherty AR, Raju M, et al. Increased risk of colorectal cancer in individuals with a history of serrated polyps. Gastroenterology. 2020;159:502-511.e2.
East JE, Vieth M, Rex DK. Serrated lesions in colorectal cancer screening: detection, resection, pathology and surveillance. Gut. 2015;64:991–1000.
Kedrin D, Gala MK. Genetics of the serrated pathway to colorectal cancer. Clin Transl Gastroenterol. 2015;6:e84–93.
Yamane L, Scapulatempo-Neto C, Reis RM, Guimarães DP. Serrated pathway in colorectal carcinogenesis. World J Gastroenterol. 2014;20:2634–40.
WHO Classification of Tumours Editorial Board. WHO classification of Tumours. Digestive System Tumours. 5th edition. Lyon: World Health Organization; 2019.
Leggett B, Whitehall V. Role of the serrated pathway in colorectal cancer pathogenesis. Gastroenterology. 2010;138:2088–100.
Patai ÁV, Molnár B, Tulassay Z, Sipos F. Serrated pathway: alternative route to colorectal cancer. World J Gastroenterol. 2013;19:607–15.
Noffsinger AE. Serrated polyps and colorectal cancer: new pathway to malignancy. Annu Rev Pathol Mech Dis. 2009;4:343–64.
dos Santos W, Sobanski T, de Carvalho AC, Evangelista AF, Matsushita M, Berardinelli GN, et al. Mutation profiling of cancer drivers in Brazilian colorectal cancer. Sci Rep. 2019;9:1–13.
Winawer S, Zauber A, Ho M, O’Brien M, Gottlieb LSSS, et al. Prevention of colorectal cancer by colonoscopic polypectomy. The National Polyp Study Workgroup. N Engl J Med. 1993;329:1977–81.
Schreuders EH, Ruco A, Rabeneck L, Schoen RE, Sung JJY, Young GP, et al. Colorectal cancer screening: a global overview of existing programmes. Gut. 2015;64:1637–49.
Robertson DJ, Lee JK, Boland CR, Dominitz JA, Giardiello FM, Johnson DA, et al. Recommendations on fecal immunochemical testing to screen for colorectal neoplasia: a consensus statement by the US multi-society task force on colorectal cancer. Gastroenterology. 2017;152:1217-1237.e3.
Imperiale TF, Ransohoff DF, Itzkowitz SH, Levin TR, Lavin P, Lidgard GP, et al. Multitarget stool DNA testing for colorectal-cancer screening. N Engl J Med. 2014;370:1287–97.
DeVos T, Tetzner R, Model F, Weiss G, Schuster M, Distler J, et al. Circulating methylated SEPT9 DNA in plasma is a biomarker for colorectal cancer. Clin Chem. 2009;55:1337–46.
Cruvinel-Carloni A, Yamane L, Scapulatempo-Neto C, Guimarães D, Reis RM. Absence of TERT promoter mutations in colorectal precursor lesions and cancer. Genet Mol Biol. 2018;41:82–4.
Yamane LS, Scapulatempo-Neto C, Alvarenga L, Oliveira CZ, Berardinelli GN, Almodova E, et al. KRAS and BRAF mutations and MSI status in precursor lesions of colorectal cancer detected by colonoscopy. Oncol Rep. 2014;32:1419–26.
Guimarães DP, Mantuan LA, de Oliveira MA, Junior RL, da Costa AM, Rossi S, et al. The performance of colorectal cancer screening in Brazil: the first two years of the implementation program in barretos cancer hospital. Cancer Prev Res. 2021;14:241–52.
Guimarães DP, Fregnani JH, Reis RM, Taveira LN, Scapulatempo-Neto C, Matsushita M, et al. Comparison of a new-generation fecal immunochemical test (FIT) with guaiac fecal occult blood test (gFOBT) in detecting colorectal neoplasia among colonoscopy-referral patients. Anticancer Res. 2019;39:261–9.
de Oliveira JC, Viana DV, Zanardo C, Santos EMM, de Paula AE, Palmero EI, et al. Genotype-phenotype correlation in 99 familial adenomatous polyposis patients: a prospective prevention protocol. Cancer Med. 2019;8:2114–22.
Schneider NB, Pastor T, de Paula AE, Achatz MI, dos Santos ÂR, Vianna FSL, et al. Germline MLH1, MSH2 and MSH6 variants in Brazilian patients with colorectal cancer and clinical features suggestive of Lynch Syndrome. Cancer Med. 2018;7:2078–88.
Participants in the Paris Workshop. The Paris endoscopic classification of superficial neoplastic lesions: esophagus, stomach, and colon. Gastrointest Endosc. 2003;58(6 Supplement):S3-43.
Neuber AC, Tostes CH, Ribeiro AG, Marczynski GT, Komoto TT, Rogeri CD, et al. The biobank of barretos cancer hospital: 14 years of experience in cancer research. Cell Tissue Bank. 2021;9:447.
Durães RO, Berardinelli GN, da Costa AM, Scapulatempo-Neto C, Pereira R, Oliveira MA, et al. Role of genetic ancestry in 1002 Brazilian colorectal cancer patients from barretos cancer hospital. Front Oncol. 2020;10:885.
Pereira R, Phillips C, Pinto N, Santos C, dos Santos SEB, Amorim A, et al. Straightforward inference of ancestry and admixture proportions through ancestry-informative insertion deletion multiplexing. PLoS One. 2012;7:3347.
Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
Fearon E, Vogelstein B. A genetic model for colorectal tumorigenesis. Cell. 1990;61:759–67.
Cross W, Kovac M, Mustonen V, Temko D, Davis H, Baker AM, et al. The evolutionary landscape of colorectal tumorigenesis. Nat Ecol Evol. 2018;2:1661–72.
Saito T, Niida A, Uchi R, Hirata H, Komatsu H, Sakimura S, et al. A temporal shift of the evolutionary principle shaping intratumor heterogeneity in colorectal cancer. Nat Commun. 2018;9:1–11.
Sottoriva A, Kang H, Ma Z, Graham TA, Salomon MP, Zhao J, et al. A Big Bang model of human colorectal tumor growth. Nat Genet. 2015;47:209–16.
Lin SH, Raju GS, Huff C, Ye Y, Gu J, Chen JS, et al. The somatic mutation landscape of premalignant colorectal adenoma. Gut. 2018;67:1299–305.
Lugli N, Dionellis VS, Ordóñez-Morán P, Kamileri I, Sotiriou SK, Huelsken J, et al. Enhanced rate of acquisition of point mutations in mouse intestinal adenomas compared to normal tissue. Cell Rep. 2017;19:2185–92.
Murakami T, Akazawa Y, Yatagai N, Hiromoto T, Sasahara N, Saito T, et al. Molecular characterization of sessile serrated adenoma/polyps with dysplasia/carcinoma based on immunohistochemistry, next-generation sequencing, and microsatellite instability testing: a case series study. Diagn Pathol. 2018;13:1–10.
Chang K, Willis JA, Reumers J, Taggart MW, San Lucas FA, Thirumurthi S, et al. Colorectal premalignancy is associated with consensus molecular subtypes 1 and 2. Ann Oncol. 2018;29:2061–7.
Komor MA, Bosch LJW, Bounova G, Bolijn AS, Delis-van Diemen PM, Rausch C, et al. Consensus molecular subtype classification of colorectal adenomas. J Pathol. 2018;246:266–76.
Druliner BR, Wang P, Bae T, Baheti S, Slettedahl S, Mahoney D, et al. Molecular characterization of colorectal adenomas with and without malignancy reveals distinguishing genome, transcriptome and methylome alterations. Sci Rep. 2018;8:1–10.
Liu Z, Yang C, Li X, Luo W, Roy B, Xiong T, et al. The landscape of somatic mutation in sporadic Chinese colorectal cancer. Oncotarget. 2018;9:27412–22.
van Lanschot MCJ, Carvalho B, Rausch C, Snaebjornsson P, van Engeland M, Kuipers EJ, et al. Molecular profiling of longitudinally observed small colorectal polyps: a cohort study. EBioMedicine. 2019;39:292–300.
Patai ÁV, Barták BK, Péterfia B, Micsik T, Horváth R, Sumánszki C, et al. Comprehensive DNA methylation and mutation analyses reveal a methylation signature in colorectal sessile serrated adenomas. Pathol Oncol Res. 2017;23:589–94.
Bailey MH, Tokheim C, Porta-Pardo E, Sengupta S, Bertrand D, Weerasinghe A, et al. Comprehensive characterization of cancer driver genes and mutations. Cell. 2018;173:371-385.e18.
Guinney J, Dienstmann R, Wang X, De Reyniès A, Schlicker A, Soneson C, et al. The consensus molecular subtypes of colorectal cancer. Nat Med. 2015;21:1350–6.
Sievers CK, Zou LS, Pickhardt PJ, Matkowskyj KA, Albrecht DM, Clipson L, et al. Subclonal diversity arises early even in small colorectal tumours and contributes to differential growth fates. Gut. 2017;66:2132–40.
Yi C, Huang Y, Yu X, Li X, Zheng S, Ding K, et al. Clinicopathologic distribution of KRAS and BRAF mutations in a Chinese population with colorectal cancer precursor lesions. Oncotarget. 2016;7:17265–74.
Rohlin A, Wernersson J, Engwall Y, Wiklund L, Björk J, Nordling M. Parallel sequencing used in detection of mosaic mutations: comparison with four diagnostic DNA screening techniques. Hum Mutat. 2009;30:1012–20.
Jamuar SS, Lam A-TN, Kircher M, D’Gama AM, Wang J, Barry BJ, et al. Somatic mutations in cerebral cortical malformations. N Engl J Med. 2014;371:733–43.
Yan HHN, Lai JCW, Ho SL, Leung WK, Law WL, Lee JFY, et al. RNF43 germline and somatic mutation in serrated neoplasia pathway and its association with BRAF mutation. Gut. 2017;66:1645–56.
Nikolaev SI, Sotiriou SK, Pateras IS, Santoni F, Sougioultzis S, Edgren H, et al. A single-nucleotide substitution mutator phenotype revealed by exome sequencing of human colon adenomas. Cancer Res. 2012;72:6279–89.
Kim R, Schell MJ, Teer JK, Greenawalt DM, Yang M, Yeatman TJ. Co-evolution of somatic variation in primary and metastatic colorectal cancer may expand biopsy indications in the molecular era. PLoS ONE. 2015;10:1–12.
Uchi R, Takahashi Y, Niida A, Shimamura T, Hirata H, Sugimachi K, et al. Integrated multiregional analysis proposing a new model of colorectal cancer evolution. PLoS Genet. 2016;12:1–24.
Konishi K, Yamochi T, Makino R, Kaneko K, Yamamoto T, Nozawa H, et al. Molecular differences between sporadic serrated and conventional colorectal adenomas. Clin Cancer Res. 2004;10:3082–90.
Muzny DM, Bainbridge MN, Chang K, Dinh HH, Drummond JA, Fowler G, et al. Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012;487:330–7.
Yamamoto E, Suzuki H, Yamano HO, Maruyama R, Nojima M, Kamimae S, et al. Molecular dissection of premalignant colorectal lesions reveals early onset of the CpG island methylator phenotype. Am J Pathol. 2012;181:1847–61.
Kim TM, An CH, Rhee JK, Jung SH, Lee SH, Baek IP, et al. Clonal origins and parallel evolution of regionally synchronous colorectal adenoma and carcinoma. Oncotarget. 2015;6:27725–35.
Whitehall VLJ, Rickman C, Bond CE, Ramsnes I, Greco SA, Umapathy A, et al. Oncogenic PIK3CA mutations in colorectal cancers and polyps. Int J Cancer. 2012;131:813–20.
Chang P-Y, Chen J-S, Chang S-C, Wang M-C, Chang N-C, Wen Y-H, et al. Acquired somatic TP53 or PIK3CA mutations are potential predictors of when polyps evolve into colorectal cancer. Oncotarget. 2017;8:72352–62.
Hashimoto T, Yamashita S, Yoshida H, Taniguchi H, Ushijima T, Yamada T, et al. WNT pathway gene mutations are associated with the presence of dysplasia in colorectal sessile serrated adenoma/polyps. Am J Surg Pathol. 2017;41:1188–97.
Hernandez-Suarez G, Sanabria MC, Serrano M, Herran OF, Perez J, Plata JL, et al. Genetic ancestry is associated with colorectal adenomas and adenocarcinomas in Latino populations. Eur J Hum Genet. 2014;22:1208–16.
Lebwohl B, Capiak K, Neugut AI, Kastrinos F. Risk of colorectal adenomas and advanced neoplasia in Hispanic, black and white patients undergoing screening colonoscopy. Aliment Pharmacol Ther. 2012;35:1467–73.
Schroy PC, Coe A, Chen CA, O’Brien MJ, Heeren TC. Prevalence of advanced colorectal neoplasia in white and black patients undergoing screening colonoscopy in a Safety-Net Hospital. Ann Intern Med. 2013;159:13.
Alves-Silva J, da Silva SM, Guimarães PEM, Ferreira ACS, Bandelt HJ, Pena SDJ, et al. The ancestry of Brazilian mtDNA lineages. Am J Hum Genet. 2000;67:444–61.
Giolo SR, Soler JMP, Greenway SC, Almeida MAA, De Andrade M, Seidman JG, et al. Brazilian urban population genetic structure reveals a high degree of admixture. Eur J Hum Genet. 2012;20:111–6.
Saulsberry L, Olopade OI. Precision oncology: directing genomics and pharmacogenomics toward reducing cancer inequities. Cancer Cell. 2021;39:730–3.
Popejoy AB, Fullerton SM. Genomics is failing on diversity. Nature. 2016;538:161–4.
The authors thank the Barretos Cancer Hospital and Biobank for the processing of blood samples. WS is a recipient of a fellowship of the postgraduate program of Barretos Cancer Hospital. MBR is a recipient of a post-doc fellowship of the Brazilian Federal Agency for the Support and Evaluation of Graduate Education (CAPES, Brazil). JP is a recipient of an undergraduate fellowship from Barretos Cancer Hospital. RMR is a recipient of the National Council for Scientific and Technological Development (CNPq, Brazil) Productivity Fellowship.
This study was supported by Public Ministry of Labor Campinas (Research, Prevention and Education of Occupational Cancer, Brazil).
Ethics approval and consent to participate
The Ethics Committee - Barretos Cancer Hospital approved the study on number ID: 1074/2016. The need for Informed consent was waived by the Ethics Committee of Barretos Cancer Hospital due to the retrospective nature of this study.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Boxplots reporting the median coverage of the 207 amplicons sequenced across the 50 genes for all precursor lesion samples. The average depth of all amplicons was 1631.5x per sample (ranging from 0x to 13823x).
Lollipop plots showing the distribution of mutations in 16 genes identified altered in colorectal cancer precursor lesions. The Y-axis represents the number of mutations at each residue. Truncating mutations are represented by black circles and green circles indicated a missense mutation. Plot from cBioPortal (http://www.cBioPortal).
Number of driver mutations detected in the sequencing differs between the four classes of precursor lesions (Kruskal-Wallis test p < 0.001). The boxplot shows the median number of mutations observed across the early adenomas, advanced adenomas, hyperplastic polyps and sessile serrated lesion. Mann-Whitney test was used to determine the statistical significance with Bonferroni multiple comparisons correction: *p < 0.05; **p < 0.001.
Sequencing metrics per sample.
Frequency of detected mutations in colorectal cancer precursor lesions.
Detected mutations in colorectal cancer precursor lesions per patient.
Frequency of the most common altered genes in early and advanced adenoma, and serrated polyps.
Association between mutation status and location for adenoma samples.
Association between mutation status and location for serrated polyp samples.
Ancestry background proportions for all four ethnic groups according to lesion type.
Ancestry background proportions for all four ethnic groups according to mutation status.
About this article
Cite this article
dos Santos, W., dos Reis, M.B., Porto, J. et al. Somatic targeted mutation profiling of colorectal cancer precursor lesions. BMC Med Genomics 15, 143 (2022). https://doi.org/10.1186/s12920-022-01294-w
- Serrated polyps
- Molecular profiling