Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2016; 7: 11743.
Published online 2016 June 15. doi:  10.1038/ncomms11743
PMCID: PMC4912618

A multigene mutation classification of 468 colorectal cancers reveals a prognostic role for APC


Colorectal cancer (CRC) is a highly heterogeneous disease, for which prognosis has been relegated to clinicopathologic staging for decades. There is a need to stratify subpopulations of CRC on a molecular basis to better predict outcome and assign therapies. Here we report targeted exome-sequencing of 1,321 cancer-related genes on 468 tumour specimens, which identified a subset of 17 genes that best classify CRC, with APC playing a central role in predicting overall survival. APC may assume 0, 1 or 2 truncating mutations, each with a striking differential impact on survival. Tumours lacking any APC mutation carry a worse prognosis than single APC mutation tumours; however, two APC mutation tumours with mutant KRAS and TP53 confer the poorest survival among all the subgroups examined. Our study demonstrates a prognostic role for APC and suggests that sequencing of APC may have clinical utility in the routine staging and potential therapeutic assignment for CRC.

Colorectal cancer (CRC) is a highly heterogeneous disease with diverse genetic and clinical features that influence therapeutic outcomes1,2. Clinicopathologic staging has been the mainstay of CRC prognosis since Cuthbert Dukes first described it in 1932 (ref. 3). With the emergence of nextgen sequencing technologies, it is now possible to extensively and rapidly evaluate a large number of genes and samples to inform its diagnosis, prognosis and potential response to therapy. Thus, a robust and detailed molecular disease classification—beyond the few currently measured molecular parameters (MSI, KRAS, NRAS and BRAF)—is urgently needed to better understand and treat this disease4,5,6. One major challenge is to identify the multivariable molecular factors that optimally predict distant metastasis and outcome. Recently, The Cancer Genome Atlas (TCGA) project has profiled and characterized the landscape of CRC mutations7. While the TCGA project robustly characterized individual mutational events in 224 samples, their interaction and impact on metastasis and survival was not fully evaluated.

Adenomatous polyposis coli (APC) has long been considered the ‘gatekeeper' gene for the vast majority of CRCs, and, for this reason, its discriminating functionality in predicting prognosis may have been overlooked in the past. APC is a tumour-suppressor gene whose protein product functions as an antagonist of the WNT signalling pathway by binding and regulating the β-catenin protein2. APC may be also involved in other cellular processes including cell migration and adhesion, transcriptional activation, apoptosis and DNA repair8,9. APC inactivation is thought to be an early event in the development of CRC, and may play a pivotal role in the initiation of the adenoma–carcinoma pathway2,10. In recent years, various studies on molecular classification of CRC have been reported5,6,7,11,12,13,14. Although APC is the most frequently mutated, known driver gene in CRC2, it has not generally been included as a factor in clinical prognostic classification and is not generally included in standard CRC-sequencing panels.

In this study, we performed TGS for 1,321 cancer-related genes and microsatellite instability (MSI) analysis on an independent set of 468 clinically characterized, sporadic, colorectal tumours. We present a multigene interactive classification analysis for colorectal carcinoma, which groups individual mutations based on their natural associations with each other. Clinical covariates, including metastasis and survival, are also provided to help further define the prognostic potential of the classification system, with APC mutations identified to play a central role. Subsequent clinical associations result in a prognostic, five-class classification that demonstrates a previously unknown, prognostic role of the APC mutation status in association with APC ‘partnering mutations'.


Identification of 17 significantly mutated genes

A targeted gene sequencing (TGS) analysis was performed on 468 colorectal tumour samples across 1,321 genes known or highly suspected to be associated with human cancer (Supplementary Table 1). To evaluate the MSI status, the Bethesda panel analysis15 was performed on all 468 CRC cases, yielding 61 (13%) MSI-H samples (called MSI henceforth). We found that 59/61 (97%) of MSI tumours had both TGFBR2 and ACVR2A insertion/deletion (indel) mutations in homopolymer regions of the targeted genes. Meanwhile, only one of three tumours with TGFBR2, but not ACVR2A, and one of six tumours with the reverse profile were MSI. Thus, the co-occurrence of mutations of these two genes nearly perfectly characterizes the MSI status. In support of this notion, ACVR2A and TGFBR2 were also reported by TCGA as the top two representative genes in hypermutated tumours, 77% of which were MSI7.

Among the 1,321 targeted genes, a robust regression analysis versus amino-acid length identified 22 genes with elevated non-silent mutation rates at a z-score of [gt-or-equal, slanted]3 (Table 1 and Supplementary Data 1), either overall, or for MSI or MSS (microsatellite stable) tumours. Our method has independently confirmed the significantly mutated genes identified by the TCGA group in non-hypermutated CRC tumours. We have re-identified all genes seen at [gt-or-equal, slanted]7% frequency in that data set. Five of our identified genes, however, had statistically significant (P<0.0001, Fisher Exact) mutation rates that were much lower in TCGA samples (Supplementary Table 2). As the TCGA mutations were determined using paired tumour/normal data, which is the more precise method for determining somatic mutations, we restricted our attention to the remaining 17 genes that had similar mutation frequencies in the TCGA data set (KRAS, TP53, APC, SMAD4, FBXW7, BRAF, TCF7L2, PIK3CA, GNAS, CBX4, ADAMTS18, TAF1L, FAM123B, CSMD3, ITGB4, LRP1B and SYNE1, given in decreasing order of mutation rate). We next calculated ranking of correlation of these genes with the four key driver genes identified in our study (BRAF, APC, KRAS and TP53), MSI status, distant metastasis status, EMT (epithelial–mesenchymal transition) and RAS/MAPK gene expression signature scores16,17 (Table 2). The genes are sorted based on increasing MSI correlation percentile, which is strongly inversely correlated with distant metastasis (Spearman r=−0.84). APC, TP53 and KRAS (and SMAD4 to a lesser extent) are shown to be among the most strongly negatively associated with MSI but positively associated with metastasis (Table 2 and Fig. 1). When limited to only MSS patients (N=407), mutations in these three genes are strongly associated with metastasis (P=0.0051, logistic regression χ2), with KRAS having the strongest individual association (P=0.0010, χ2 test). Notably, our large cohort confirms the inverse relationship of APC mutations with MSI previously reported by others7,18,19,20,21. We selected these frequently mutated APC, TP53 and KRAS genes to characterize this MSS group. TCF7L2 and FAM123B are WNT signalling genes and have much weaker tendencies, but tend to associate with the genes listed above. By contrast, BRAF, ITGB4, CBX4, CSMD3, SYNE1, FBXW7 and TAF1L are strongly correlated with MSI and non-metastatic leaning genes. The remaining four genes, including PIK3CA, lie in between these contrasting groups (Table 2 and Fig. 1).

Figure 1
Multidimensional scaling plot for the correlation of highly mutated genes.
Table 1
Genes with elevated non-silent mutation rates.
Table 2
Percentile ranking of Spearman correlation among 1,082 genes.

Four driver genes had pairwise and significant correlations

As the single gene mutation status does not capture the full complexity of a tumour, we classified our samples via mutation combinations in those genes determined to be commonly mutated, but varying across samples. Among the 17 genes, the most striking negative Spearman correlations were BRAF and KRAS, and BRAF and APC (both with P<0.0001), while the highest positive association (Spearman r=0.26, P<0.0001) was between APC and KRAS (Table 2). The only remaining positively associated gene with APC with Fisher P<0.01 was TP53. Given the strength of these associations and their high mutation frequencies, these four genes were used to comprise the basis of a 10-group classification as shown in Supplementary Fig. 1A. Notably, our data and TCGA data have similar distribution patterns across multigene mutation groups in the 10-group analysis. In our cohort, the distribution patterns of MSI and MSS tumours differ remarkably (Supplementary Fig. 1B), supporting the notion that these CRC types undergo tumorigenesis through very different genetic mechanisms, and, accordingly, have different clinical outcomes22,23.

BRAF and APC/KRAS/TP53 partnering defines MSI and metastasis

Demographic features of the 10 groups are shown in Supplementary Table 3. While the 10-group analysis lays out the mutational patterns in detail, some trending becomes clearer by counting the number of mutated partnering genes (APC, KRAS and TP53). Eighty-four per cent of BRAF-only tumours are MSI compared with only 50% when a TP53 mutation is also present. Of these 45 BRAF-grouped tumours (B or BP), only 9 (20%) were metastatic. In contrast, among the 199 tumours with both APC and TP53 mutations, only 1% (2) are MSI, while 43% have metastatic disease.

The observed counts of the four APC-mutated groups (A, AK, AP and AKP) show that APC usually co-occurs with either KRAS or TP53 mutations, or both (Supplementary Table 3); of 312 such tumours, only 43 (14%) are in Group A (APC only), suggesting that APC mutations need to partner with one or more additional driver mutations to advance to CRCs. Interestingly, the triple-mutant AKP (APC/KRAS/TP53) group had the highest metastatic rate, with 22% developing metastatic disease from a stage 1–3 primary, suggesting a possible evolution from either AK or AP mutants by adding the third partner mutation. Primary site location might also be a deciding factor here (Supplementary Table 3). Notably, in the caecum and ascending right colon, AKPs outnumber APs and AKs, whereas in the other locations the number of APs always exceeded that of AKPs, especially in the rectosigmoid and sigmoid locations.

While APC mutations were highly associated with both KRAS (Fisher P<0.0001) and TP53 (Fisher P=0.0042) mutations, the latter are marginally negatively associated with each other (Fisher P=0.066), as has been reported by other groups4,20. This negative association is confined to the APC-mutated group (Fisher P=0.0005; Supplementary Table 4). These observations (APC is less frequently mutated alone, more commonly mutated with KRAS, TP53 or both) suggest that the APC mutation occurs early in carcinogenesis, a result corroborated by the high frequency of APC mutations in adenomas2,10. However, it is noteworthy that, although APCKRASTP53 is a widely perceived sequential pathway of the key driver mutation events commonly occurring in CRC development1,10,24, our data provide further clues in the development of CRC involving multiple pathways, as 156 tumours lack an APC mutation.

Analysis of APC mutational status

We further analysed APC mutational status in association with MSI, allelic loss and other non-APC WNT pathway mutations (Table 3 and Supplementary Table 5). Among APC-mutated tumours, 38% (118/312) displayed two or more truncating mutations. Only five tumours had more than two APC mutations. These data provide a compelling case for the significance of the second mutation, supporting the bi-allelic nature of APC because of its deviation from a Poisson distribution (χ2 goodness of fit P=2.6 × 10−11). We then divided the APC gene into four zones based on their truncating mutation rates (Table 3). Zones 1 (codons 1–218) and 4 (codons 1,588–2,843) had very few such mutations (N=3 for zone 1; N=2 for zone 4, both in one POLE tumour that also had two other APC-truncating mutations). Zone 3 (codons 1,263–1,587) had a high mutation rate (259 mutations over 325 amino acids), and this zone has been deemed the ‘mutation cluster region'25 because of the high rate of germline mutations in familial adenomatous polyposis (FAP)-associated colon cancer. Zone 2 (codons 219–1,262) had a moderately high rate of APC mutations (N=181 over 1,044 amino acids). With striking statistical significance (χ2 P<1 × 10−15), tumours with two APC mutations had one mutation each from zones 1–3. This is in accord with the double mutation profile established by others26,27, more specifically for FAP colon cancers. Zone 3 mutations would yield 1–3 (out of 14) functional 20-amino-acid β-catenin-binding sites per molecule and zone 1–2 mutations would remove all β-catenin-binding sites. This further supports the notion that ‘just right' WNT/beta-catenin signalling may be necessary for advanced progression of CRC21,26,27. Interestingly, most rare exceptions to this APC double-hit pattern were either in variant isoforms or were hypermutated POLE tumours7 (Table 3).

Table 3
APC truncation mutation profiles.

On the other hand, it was also observed that almost 30% of tumours harboured only one APC mutation without inferred allelic loss, most of which were MSS tumours (Table 3), with a substantial number also lacking any other WNT pathway mutation from our examined list. Analysis of the published TCGA data set (including somatic mutations and copy-number (CN) alterations) shows a similar result: 35% (74/209) of tumours harbour a one-hit mutation in APC without CN loss (Supplementary Table 6), which actually exceeds the number of samples with two APC hits (mutations and/or CN loss): 34% (71/209). Moreover, an independent mutation and a loss of heterozygosity (LOH) analysis of 62 additional CRC tumours were carried out on the selected drivers. While APC, TP53, KRAS/NRAS and BRAF showed similar mutation patterns, the rate of APC allelic loss appeared to be low and was observed only in tumours with APC 0- and 1-truncating mutation (Supplementary Data 2).

Immunohistochemical analysis of β-catenin nuclear staining

To determine whether the WNT/β-catenin pathway might be activated in tumours with one APC-truncating mutation, we performed an immunohistochemical (IHC) analysis of 52 FFPE (formalin-fixed, paraffin-embedded) tumours (selected from 468 CRCs) with a β-catenin antibody. Results show that β-catenin nuclear staining in APC single-mutation tumours was significantly stronger than APCwt tumours (Welch t-test P=0.014, Fig. 2a,b), indicating increased β-catenin activities in these APC-mutated tumours. For example, a ‘zero' nuclear staining score was seen in six APCwt tumours (out of 14 cases), but in only one single-mutation tumour (out of 17 cases; Supplementary Table 7). Significantly increased nuclear staining was not demonstrated in tumours with either APC multiple mutations (n=12, Welch t-test P=0.37) or single mutation plus inferred allelic loss (n=9, Welch t-test P=0.17), possibly because of small sample sizes (Fig. 2b).

Figure 2
Immunohistochemical staining of β-catenin in 52 colorectal cancers.

Activation of the WNT/β-catenin pathway in APC subgroups

We further investigated WNT/β-catenin activation in 458 CRCs (10 samples without suitable microarray data were excluded), which were divided into four APC groups: (1) APCwt (n=151), (2) 1 truncating mutation without inferred allelic loss (n=124), (3) 2+ truncating mutations (n=124) and (4) 1 truncating mutation plus inferred allelic loss (n=59). A complete set of 64 β-catenin (upregulated)-targeted genes identified in CRC cell lines by Herbst et al.28 was used as a measure of WNT pathway activation. Results show that the WNT targets were significantly activated in APC single-mutation tumours without allelic loss (Welch t-test P<0.0001), and in tumours with APC 2+ mutations (Welch t-test P=0.0004) or with single mutation plus allelic loss (Welch t-test P=0.011; Supplementary Fig. 2A). In addition, no statistically significant difference was seen between zone 1–2 and 3–4 mutations in tumours containing APC-one-truncating mutation either without or with inferred allelic loss (Supplementary Fig. 2B). These data support the notion that the WNT/β-catenin pathway was significantly activated in the APC-mutated tumours, which appeared to be zone-independent.

APC expression and methylation and other WNT pathway genes

We also investigated a possible role of methylation of the APC gene as a ‘second hit' using the TCGA CRCs7, of which 209 samples had mutation, CNA and RNAseq expression data. Here we identified ‘hits' in APC by counting the number of truncating mutations and CN deletions in each sample. We then compared the number of hits to mRNA expression (RNASeq) of the genes of interest. We first examined APC mRNA expression and observed a general trend of decreasing APC expression, with increasing number of APC hits (Supplementary Fig. 3). Differences between 0, 1 and 0, 2 APC hits were significant, but only the difference between 0 and 2 hits was significant after Bonferroni–Holm correction. This difference could be related to mRNA instability induced by truncation. We then examined APC methylation profiles and the relationship to APC expression in each APC mutation group. Promoter-specific probe cg15020645 (Supplementary Fig. 4 and Methods) showed the bimodal pattern often observed in promoter methylation; however, no significant differences were observed between APC-hit groups using the Wilcoxon rank-sum test (Supplementary Fig. 5A). Although the methylation pattern was different in probe cg01240931 (Methods), no significant differences were observed between APC-hit groups (Supplementary Fig. 5B). Indeed, neither probe showed any difference between one and two APC hits, strongly suggesting that APC promoter methylation is not a second hit in APC one-hit APC samples. In addition, we observed no significant correlation between APC promoter methylation and APC gene expression in either APC one-hit or all samples (Supplementary Fig. 6).

In addition to APC methylation, we also investigated a possible role of non-APC-mutated WNT pathway gene mutations as a ‘second hit' in one-hit- APC mutation tumours using the TCGA CRCs that were available with whole-exome-sequencing data. Of 67 WNT pathway genes examined, 63 had mutations and were analysed (Methods). We found no significant correlation between the number of mutated non-APC WNT pathway genes and the expression of the 64 β-catenin-targeted genes (Pearson r=−0.03, P=0.67, Supplementary Fig. 7), suggesting that these non-APC WNT gene mutations are not responsible for WNT pathway activation. We then assessed the frequency of these non-APC WNT gene mutations versus the number of APC hits, and found that more mutations were associated with zero-hit APC (Supplementary Fig. 8). In addition, no differences were seen between one- and two-hit APC groups. Taken together, mutations in WNT pathway genes other than APC did not appear to be a ‘second hit' in the majority of one-hit APC samples.

Prognostic role of APC and the five-class classification

For prognostic analysis of 468 CRCs, we first divided patients by their MSI status for clinical evaluation because of the strong negative association between the MSI status and distant metastasis, and their sharp contrast in mutation counts (median=169 for MSI and 33 for MSS tumours). Among MSI tumours, normally considered to carry a good prognosis overall22,23, Kaplan–Meier analysis shows that BRAF (V600E)-mutated tumours have marginally significantly poorer overall survival (log-rank P=0.090; Fig. 3a).

Figure 3
Kaplan–Meier survival analysis.

In MSS tumours, overall survival was roughly equivalent among the APC, APC/KRAS and APC/TP53 groups, but lower in the APC/KRAS/TP53 (AKP) group, which was nearly equivalent to that seen among APCwt patients (log-rank P=0.0090; Fig. 3b). Notably, our association analyses with two signature scores indicate that, while the AKP triple mutants had higher RAS/MAPK activity16, APCwt tumours were associated with a higher EMT score17,29 (Supplementary Fig. 9), suggesting that AKP and APCwt might undergo distinct mechanisms of tumorigenesis and thereby differ in drug responsiveness. In support of this notion, our previous study showed that EMT was inversely associated with RAS activity17. Moreover, APC two-mutation- and zero-mutation tumours showed worse survival compared with APC one mutation tumours (log-rank P=0.018; Fig. 3c). Since both AKP triply mutated and other APC-mutated tumours had either one or two APC mutations, we further re-grouped these MSS tumours into five classes: Class 0 (n=111): ‘APCwt'; Class 1 (n=135): ‘APC(1), APC(1)/KRAS, APC(1)/TP53'; Class 2 (n=76): ‘APC(2), APC(2)/KRAS, APC(2)/TP53'; Class 3 (n=45): ‘APC(1)/KRAS/TP53' (APC(1)KP for short); and Class 4 (n=40): ‘APC(2)/KRAS/TP53' (APC(2)KP for short). Here APC(1) represents one APC mutation and APC(2) represents two APC mutations. This five-class classification incorporates both APC mutation count and partnering mutation count, and results in more significantly resolved survival analysis (log-rank P=0.0020; Fig. 3d). Compared with Class 1 (reference) tumours, while APCwt (Class 0) tumours still carried substantially worse survival (hazard ratio (HR)=1.94, χ2 P=0.0023), the AKP triply mutated tumours with two APC mutations (Class 4) had the worst survival (HR=2.48, χ2 P=0.0011). However, the AKP tumours with one APC mutation (Class 3) and other APC two-mutation (Class 2) tumours were no longer significantly associated with worse survival (HR=1.48, χ2 P=0.17 and HR=1.11, χ2 P=0.67, respectively) in this refined model, suggesting that either AKP or two APC mutations alone were separately insufficient for prediction of worse survival.

The five-class classification for all patients (n=468) is given in Table 4. To biologically analyse this differential behaviour, expression comparison between the five classes was also performed on 399 MSS tumours (from 458 CRCs with complete data) using the EMT17 and RAS/MAPK16 signatures as well as the 64 β-catenin-targeted genes28 (Fig. 4a–c). Notably, despite their different impacts on survival, no significant differences were seen between Class 3 and Class 4 in the EMT scores, RAS scores or the expression of the WNT/β-catenin targets (Fig. 4a–c). We also assessed APC mRNA expression in the five classes and found that Class 4, which had the worst survival, also had the highest mRNA expression, which was significantly above the less lethal classes (Class 1, Welch t-test P<0.0001, Class 2, P=0.011 and Class 3, P=0.0022) but was also marginally significantly higher than the Class 0, which also had poor survival (Welch t-test P=0.062; Fig. 4d). Notably, APC mRNA expression of the tumours that had two APC mutations was significantly higher than that of their respective one APC mutation counterparts (that is, Class 2 versus Class 1 and Class 4 versus Class 3); however, Class 4 tumours also had stronger RAS pathway activation and higher expression of WNT targets than Class 2 tumours (Fig. 4b,c).

Figure 4
Expression (log2) comparison of MSS tumours in five classes (n=399).
Table 4
Five classes grouped by the APC mutation status with partnering mutations.

Cox model predicators of survival and CMS classification

Cox modelling was used to determine the effect of the four key genes on overall survival among MSS tumours, with only an adjustment for older age (Table 5; Model 1). An APC mutation was found to be a moderate good risk factor, with an HR of 0.65 (χ2 P=0.021), while tumours with KRAS and TP53 mutations were in moderate poor risk, with HRs in the 1.4–1.6 range; a BRAF mutation, although rare (only 18/407 tumours, 4.4%), had a dramatic HR of 2.63 (in agreement with several recent reports indicating BRAF(V600E) as a strong negative prognostic marker of CRC30,31,32,33,34,35), while the HR was 1.05 per year over the age of 70. Adding in terms for the AKP group and for double APC mutation tumours (Table 5; Model 2), while allowing individual mutation effects to drop out, suggests that the increased hazard for KRAS- or TP53-mutated tumours may be largely borne by those patients with both. Including the presence of metastatic disease (Table 5; Model 3) resulted in a better fit (70 point rise in the LR χ2) and a 4.51 HR for distant metastasis. Notably, adding metastasis led to a significant HR for two-hit APC mutation tumours (HR=1.58, χ2 P=0.027, Model 3) but decreased the HR of AKP from 1.79 (χ2 P=0.005, Model 2) to 1.43 (χ2 P=0.090, Model 3) probably because of the very strong correlation of AKP with distant metastasis. Furthermore, addition of APC(2)KP (Class 4) that replaced AKP and two-hit APC mutations showed a substantially significant HR (HR=2.10, χ2 P=0.003; Table 5; Model 4), supporting Class 4 tumours as a worst-risk class as shown earlier in the Kaplan-Meier (KM) survival analysis (Fig. 3d).

Table 5
Cox model predictors of overall survival of MSS tumours (N=407).

Recently, an international consortium36 has coalesced six independent (gene expression) classification systems of CRCs into four consensus molecular subtypes (CMS1–4) and a fifth unclassified group (CMS_NA). We applied this system to classify 458 CRCs. The Random Forest classification (RF) and the single sample predictor (SSP) classification were used, both of which captured the mutation features of the common drivers (APC, KRAS, TP53 and BRAF) similar to those reported36. A combined RF–SSP classification is shown in Table 6. In terms of significantly higher or lower observations than expectations, CMS1 and CMS2 showed opposite patterns for a majority of categories listed. In relationship with the five-class classification, Class 0 and Class 1 showed significantly higher or lower observations in two CMS subtypes (CMS1 and CMS2) with opposite directions (Supplementary Table 8).

Table 6
CMS classification of 458 CRCs.

Furthermore, Cox modelling of 399 MSS tumours indicated that, compared with CMS2, worse survival was significantly associated with CMS1 (HR=3.48, χ2 P<0.0001) and CMS_NA (HR=1.82, χ2 P=0.02), and was marginally associated with CMS3 (HR=1.76, χ2 P=0.06) but not with CMS4 (HR=1.34, χ2 P=0.22; Supplementary Table 9; Model 3). However, only the CMS1 tumours remained significantly associated with worse survival in the multivariable model upon adding APC and BRAF plus AKP and APC two-mutation (HR=2.39, χ2 P=0.012, Model 4) or plus APC(2)KP (HR=2.55, χ2 P=0.007, Model 5).


Our robust regression analysis of mutations versus amino-acid length of 468 CRC tumours objectively identified a subset of 17 genes that are commonly mutated in CRC. The striking pairwise, statistically significant, correlations of APC, KRAS, TP53 and BRAF reflect a positive or negative interaction between these common drivers in CRC. These observations have shed light on a crosstalk between multiple pathways deregulated by these driver mutations, with a new focus on the role of APC.

Heretofore, APC was presumed to be an important ‘initiator' gene for the majority of CRCs2,10; thus, its newly discovered prognostic role was likely underestimated. Only with the evaluation of APC in the context of its partnering mutations (with KRAS and TP53) and its bi-allelic mutation status does a prognostic role emerge for APC, with implications for tumour evolution. Survival analysis indicates that two APC-truncating mutations, in the presence of mutant KRAS and TP53 (AKP), carry a substantially worse prognosis than single truncating mutations, but surprisingly are equivalent to highly lethal tumours lacking any APC mutation. These data suggest that wild-type APC tumours are driven by an ‘alternate', non-WNT lethal pathway. It is noteworthy that no prognostic role for APC mutations was previously reported in a study of 107 CRC patients37. The identification of the new prognostic role of APC may be attributed to a more detailed mutational analysis and to a larger sample size.

While tumours with two different APC mutations may fit the bi-allelic theory, tumours with one mutation are of particular interest because they have improved overall survival and statistically significant different biological profiles from both APCwt and two-hit APC tumours. The WNT/β-catenin pathway was activated in these one-hit APC tumours as measured by IHC nuclear staining of a β-catenin antibody and expression of the 64 β-catenin-targeted genes28. Intriguingly, the β-catenin nuclear staining scores and the expression of the 64 targeted genes appear to be less significant in tumours with two-hit APC mutation or one-hit mutation plus inferred allelic loss (Fig. 2b and Supplementary Fig. 2A), suggesting a possible ‘feedback' regulation of WNT activation via non-APC mechanisms in these tumours. This is consistent with the notion of ‘just right' WNT/beta-catenin signalling21,26,27.

To interpret our observation on WNT activation in one-hit APC mutation tumours using the classic two-hit model of APC mutation in FAP, we considered several possible mechanisms as a ‘second hit'. First, allelic loss, which has been well documented to be interdependent with a single APC mutation20,38,39, and somatic recombination (Robertsonian translocation)40, might account for the ‘second hit' in 59 of 187 one-hit APC mutation samples. This is predicted by examining alternate read percentages (those mutations with percentages [gt-or-equal, slanted]50% suggest such loss of heterozygosity and/or recombination); however, this methodology could not distinguish these two possible mechanisms.

Second, a secondary WNT pathway gene mutation might contribute to tumorigenicity in remaining single APC-mutated tumours. An analysis of the TCGA CRCs, whose whole-exome-sequencing data allowed us to examine WNT pathway genes more comprehensively, did identify a substantial number of APCwt (~53%), one-hit (~33%) or two-hit (~30%) APC tumours, which had one or more non-APC WNT pathway gene mutations. These ‘other' WNT mutations, however, did not appear to significantly contribute to WNT/β-catenin activation. It is also noteworthy that the majority of other WNT pathway alterations were missense mutations that are predicted to be neutral or to have low functional impact. In addition, a great majority of the ‘other' WNT pathway genes are known to positively contribute to WNT pathway activation. Thus, inactivating (truncating) mutations in the WNT pathway genes might more likely suppress WNT activity rather than activate it in the majority of CRCs.

Third, methylation has been thought to play an important role in CRC. This notion is supported by the fact that the CMS1 tumours (that had a significant lower rate of APC mutations) show a hypermethylation profile36, potentially initiating tumorigenesis when APC mutation is absent. Hypermethylation of the APC promoter in CRCs has been reported previously41,42,43,44. Our analysis on the TCGA CRCs7 also identified a significant number of methylation events associated with the APC promoter, but, importantly, no statistically significant difference was observed between zero-, one- and two-hit APC tumours. Moreover, no significant correlation was obtained between APC methylation and mRNA expression. Thus, our data do not support methylation of the APC promoter as a ‘second hit' in one-hit APC tumours. This agrees with several recently published studies showing that promoter methylation in APC plays an insignificant role in substituting for truncating mutations44,45. Collectively, we have identified a significant number of tumours harbouring one-hit APC mutations, many of which retained WNT activation that could not be explained by a few possible ‘second-hit' mechanisms.

A refined survival analysis based on the APC mutation status and partnering mutated genes (KRAS, TP53) led us to develop a prognostic, five-class genetic classification of CRC. When all patients are considered, Class 0 tumours are significantly associated with MSI that is also strongly correlated with BRAF mutation. Notably, a recent report on a BRAF(V600E) mouse model suggests that BRAF may be an early driver gene for MSI-like tumours46. Class 0 MSS are predicted to be mesenchymal-like and are positively correlated with a ‘composite' expression signature score that predicts poorer outcome in stage 1–3 patients29. Notably, most of the CRC classification studies report that MSI and EMT are the only two subtypes that have been consistently shown to have a statistically significant, prognostic role5,6,7,11,12,13. For Class 4 MSS tumours, an explanation for their worse survival is not evident, as no significant difference is seen between this class and its one APC mutation counterpart Class 3 in terms of EMT, RAS and WNT/β-catenin expression signatures. However, we found that APC mRNA expression of Class 4 is significantly higher than that for Classes 1–3 and marginally significantly higher than Class 0.

It has been previously reported that APC expression may be associated with another function of APC by binding DNA polymerase β (Pol-β) and flap endonuclease 1 (Fen-1) to inhibit their DNA repair activities9,47,48,49,50,51. Since the DNA-repair-inhibitory domain of APC was located upstream of mutation cluster region (MCR)9,49, most two-hit APC mutation tumours could retain an intact DNA-repair-inhibitory domain on the allele(s) with a mutation at or downstream of the APC ‘mutation cluster region' (identified here as zones 3 and 4). Thus, a possible linkage of Class 4 tumours' worse survival with APC's DNA-repair-inhibitory function needs to be further investigated.

Using the recently reported CMS classification system36, we found that CMS1, instead of CMS4 MSS tumours, were significantly associated with worse survival. This difference in CRC prognostic properties between our analysis and reported results36 might result from the highly heterogeneous nature of CRCs.

In conclusion, we have developed a prognostic, five-class, multigene mutation classification system for CRC, with APC playing a central role. We believe that routine clinical APC mutation assessment, in addition to other known classifiers BRAF, KRAS and TP53, is useful and expedient in predicting outcomes and may ultimately help improve the resolution of appropriate therapy by identifying high- versus low-risk subpopulations.


Tumour specimen and DNA extraction

A cohort of the 468 colorectal adenocarcinoma patients for whom primary or metastatic tissue (367 primary lesions from stage 1–4 patients and 101 metastatic lesions) was obtained between October 2006 and September 2010 and subjected to TGS was included in the study. Primary and metastatic samples were both included on the basis of our previous work, demonstrating a high degree of mutation overlap between matched primary/met samples52. In all cases, tissue and clinical data were collected from patients under the approval of the University of South Florida institutional review board, and informed consent was obtained from the participating patients. All the tumours were collected from curative survival resections and snap-frozen in liquid nitrogen within 15–20 min of extirpation. Tumours then underwent a macrodissection quality-control process to ensure that >80% tumour was present in the specimen that underwent sequence analysis. Normal tissue, necrotic tissue and excessive stromal tissues were dissected away from the specimen under frozen-section control. DNA was then extracted from the specimen for subsequent TGS, comprising 1,321 cancer-associated genes (Supplementary Table 1) that were selected by a joint committee (Merck & Moffitt Cancer Center).

Independent targeted sequencing cohort

Overall, 62 primary colon adenocarcinomas (including all histologic subtypes) were profiled by Molecular Health, utilizing a comprehensive 613 gene panel (The Woodlands, TX). The assay utilizes custom hybrid-capture, full exonic sequencing on the Illumina 2500 instrument (Illumina Inc., San Diego, CA). Proprietary analytic pipelines were used to map and call variants (single-nucleotide variants, frameshift mutations and CN alterations).

The mean age of the cohort is 59.9 years, with a range of 37–89 years. The mean age of patients with zero, one and two APC mutations is 59.0, 58.0 and 63.4 years, respectively. Results are shown in Supplementary Data 2.

Sequence and data analysis

A median of 15,530,823 total reads (1.40 Gbases) was generated for each sample across 1,321 genes covering 3.8 MB, which were targeted using SureSelect custom designs (Agilent Technologies Inc., Santa Clara, CA) using GAIIx sequencing technology (Illumina Inc.) by BGI (Shenzhen, China). The Burrows-Wheeler Aligner (BWA53) was used to align sequences to human reference hg19. The Genome Analysis ToolKit (GATK)54 was used for insertion/deletion realignment, quality score recalibration and variant identification. The median of average depth of coverages across samples was 140.54 × . Depth was high across the targeted region: the median percentage of targeted bases covered by [gt-or-equal, slanted]10 and [gt-or-equal, slanted]20 high-quality bases were 94.02% and 90.06%, respectively (Supplementary Fig. 10). ANNOVAR55 was used to annotate mutations. Matched normal samples were generally not available for comparison; somatic mutations were enriched by removing variants found at >1% in the 1000 Genomes project data set. In addition, 264 non-matching normal tissue samples were sequenced with the colon tumour samples, and were used to filter out normal single nucleotide variants (SNVs) and potential artefacts. Finally, data from 523 normal tissue or blood samples from the TCGA breast data set were downloaded from dbGAP and analysed using the same BWA/GATK pipeline as our sequenced samples. A total of 21,179 non-silent variants were identified as being present in normal tissues and were used for filtering suspected germline variants. Mutation counts were then based on filtered data. Mutations were manually examined with samtools tview56 to identify and remove suspected alignment and other artefacts. A detailed description of 22 significantly mutated genes is given in Supplementary Data 1.

Statistical methods

Robust regression was performed using the ROBUSTREG procedure (SAS 9.1, Cary, NC), with the MM approach. A value of 3 or more was used to define outlier genes with elevated mutation rates. Associations between the presence and absence of mutations in different genes, and comparison of mutation frequencies between the Moffitt 468 data set and the TCGA data for CRC was made using the Fisher exact test. Survival curves were constructed using the Kaplan–Meier method, and the log-rank test was used to assess statistical significance. The living patients were followed for a median of 34 months. Four Cox proportional hazards models were developed to explore the prognostic significance of the four key driver genes. A piecewise linear adjustment, called Age [gt-or-equal, slanted]70, which is 0 for patients <70 and linearly increasing for those 70 or more, was used to adjust for an increasing likelihood of death based on age. Model 2 adds in the effect of the AKP group (defined below) and presence of two or more APC truncation mutations, while model 3 shows the impact of metastatic status on the model. Both models 2 and 3 used backward elimination, with a significance-level-to-stay of 0.10 to remove nonsignificant factors. In Model 4, AKP and APC 2 mutations were replaced by APC(2)KP. A similar analysis was also carried out with CMS subtypes.

IHC analysis of β-catenin staining

Overall, 52 FFPE samples were selected from 468 CRCs for IHC analysis using a β-catenin antibody. Briefly, the FFPE slides were stained using a Ventana Discovery XT automated system (Ventana Medical Systems, Tucson) as per the manufacturer's protocol with proprietary reagents. The mouse monoclonal antibody for β-catenin (#760-4242, Ventana, 50 tests per vial at predilute concentration) was used by a Ventana anti-mouse secondary antibody. IHC staining was analysed by a pathologist, and Allred scores were given (Supplementary Table 7).

Expression of β-catenin-targeted genes

A set of 64 ‘consensus' β-catenin (upregulated) genes were adopted from a recent study of Herbst et al.28 A mean Log2 expression was calculated for each set of genes in APC subgroups. P-values for comparison were obtained by two-tailed, unequal variance Welch t-test.

EMT and RAS signature scores

EMT and RAS signature scores were previously developed16,17 and calculated in 468 samples similarly as described in our recent report29.

Analyses of TCGA CRCs

To determine the potential contribution of APC CN deletion, truncation, methylation and mutations of other WNT pathway genes, we have examined data from an earlier CRC publication by the TCGA consortium7. Level 3 data were downloaded The mRNA expression was expressed by the mean of (log2 rpkm+0.00748 floor) of RNASeq data. For promoter methylation analysis, we first identified six probes at the 5′ end of the APC gene (Supplementary Fig. 4); five of these probes behaved similarly to each other and the methylation ratio beta values were highly correlated (range of pairwise Pearson r=0.964–0.978). We selected cg15020645 as a representative of this ‘promoter' group. Probe cg01240931 was located further inside the gene, and it showed a more modest correlation with the representative ‘promoter' group probe of 0.52; therefore, we examined this probe as well. We then plotted the distribution of methylation beta values in each of the APC ‘hit' groups. Notably, one or more mutations were observed in TCGA samples for 62 non-APC WNT pathway genes (out of total 65 genes adapted from Nusse's the Wnt homepage In addition, 10 missense and 2 truncating mutations were also observed in RNF43, but no mutation was seen in its closely related ZNRF3, both of which were recently reported to play an important role in negatively regulating WNT pathway57,58,59,60. Thus, the number of mutations for these 63 genes was applied for correlation or comparison analysis. P-values were obtained from the exact Wilcoxon rank-sum test.

Additional information

Accession codes: The targeted exome sequencing data have been deposited in the database of Genotypes and Phenotypes (dbGaP) under accession code phs001111.v1.p1.

How to cite this article: Schell, M. J. et al. A multigene mutation classification of 468 colorectal cancers reveals a prognostic role for APC. Nat. Commun. 7:11743 doi: 10.1038/ncomms11743 (2016).

Supplementary Material

Supplementary Information:

Supplementary Figures 1-10, Supplementary Tables 1-9 and Supplementary References.

Supplementary Data 1:

Targeted exome sequencing data of 22 significantly mutated genes in 468 CRCs.xls

Supplementary Data 2:

Targeted gene sequencing and LOH data of driver genes in 62 CRCs.xlsx


This work was supported by NIH grant U01CA157960 (to T.J.Y.). We also acknowledge Total Cancer Care, which is enabled, in part, by the support of the DeBartolo Family, and we thank the patients who provided data and tissue to the Total Cancer Care Consortium. Our study also received assistance from the Biostatistics, Cancer Informatics, Collaborative Data Services and Tissue Core Facilities at the H. Lee Moffitt Cancer Center and Research Institute, an NCI-designated Comprehensive Cancer Center, supported under NIH grant P30-CA76292. We thank Jiannong Li (Moffitt Cancer Center, Tampa, FL) for assistance in statistical analysis and Sandra Indiviglio (Molecular Health, The Woodlands) for assistance in mutation analysis.


MVN and AL are employees of Merck, Sharpe and Dohme. D.M.G was an employee of Merck, Sharpe and Dohme when this work was carried out. GAB-W is an employee of Molecular Health. The remaining authors declare no competing financial interests.

Author contributions M.J.S., M.Y., J.K.T., F.Y.L., A.M., D.C., A.N.A.M., M.V.N., B.Y., A.L., G.A.B.-W., D.M.G. and T.J.Y. designed and performed experiments, and analysed data; M.J.S., M.Y. and T.J.Y. wrote the manuscript; J.K.T., F.Y.L., A.M., D.C., A.N.A.M., M.V.N., B.Y., A.L., G.A.B.-W. and D.G helped in revising the manuscript; M.J.S. and T.J.Y. conceived the hypothesis and co-supervised the project.


  • Markowitz S. D. & Bertagnolli M. M. Molecular origins of cancer: molecular basis of colorectal cancer. N. Engl. J. Med. 361, 2449–2460 (2009). [PMC free article] [PubMed]
  • Fearon E. R. Molecular genetics of colorectal cancer. Annu. Rev. Pathol. 6, 479–507 (2011). [PubMed]
  • Dukes C. Histological grading of rectal cancer: (section of pathology). Proc. R Soc. Med. 30, 371–376 (1937). [PMC free article] [PubMed]
  • Domingo E. et al. . Use of multivariate analysis to suggest a new molecular classification of colorectal cancer. J. Pathol. 229, 441–448 (2013). [PMC free article] [PubMed]
  • De Sousa E. M. F. et al. . Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions. Nat. Med. 19, 614–618 (2013). [PubMed]
  • Sadanandam A. et al. . A colorectal cancer classification system that associates cellular phenotype and responses to therapy. Nat. Med. 19, 619–625 (2013). [PMC free article] [PubMed]
  • Cancer Genome Atlas, N. Comprehensive molecular characterization of human colon and rectal cancer. Nature 487, 330–337 (2012). [PMC free article] [PubMed]
  • Brocardo M. & Henderson B. R. APC shuttling to the membrane, nucleus and beyond. Trends Cell Biol. 18, 587–596 (2008). [PubMed]
  • Jaiswal A. S. & Narayan S. A novel function of adenomatous polyposis coli (APC) in regulating DNA repair. Cancer Lett. 271, 272–280 (2008). [PMC free article] [PubMed]
  • Fearon E. R. & Vogelstein B. A genetic model for colorectal tumorigenesis. Cell 61, 759–767 (1990). [PubMed]
  • Budinska E. et al. . Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer. J. Pathol. 231, 63–76 (2013). [PMC free article] [PubMed]
  • Marisa L. et al. . Gene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic value. PLoS Med. 10, e1001453 (2013). [PMC free article] [PubMed]
  • Roepman P. et al. . Colorectal cancer intrinsic subtypes predict chemotherapy benefit, deficient mismatch repair and epithelial-to-mesenchymal transition. Int. J. Cancer 134, 552–562 (2014). [PMC free article] [PubMed]
  • De Roock W. et al. . Effects of KRAS, BRAF, NRAS, and PIK3CA mutations on the efficacy of cetuximab plus chemotherapy in chemotherapy-refractory metastatic colorectal cancer: a retrospective consortium analysis. Lancet Oncol. 11, 753–762 (2010). [PubMed]
  • Rodriguez-Bigas M. A. et al. . A National Cancer Institute Workshop on Hereditary Nonpolyposis Colorectal Cancer Syndrome: meeting highlights and Bethesda guidelines. J. Natl Cancer Inst. 89, 1758–1762 (1997). [PubMed]
  • Dry J. R. et al. . Transcriptional pathway signatures predict MEK addiction and response to selumetinib (AZD6244). Cancer Res. 70, 2264–2273 (2010). [PMC free article] [PubMed]
  • Loboda A. et al. . EMT is the dominant program in human colon cancer. BMC Med. Genomics 4, 9 (2011). [PMC free article] [PubMed]
  • Jass J. R. et al. . APC mutation and tumour budding in colorectal cancer. J. Clin. Pathol. 56, 69–73 (2003). [PMC free article] [PubMed]
  • Salahshor S. et al. . Colorectal cancer with and without microsatellite instability involves different genes. Genes Chromosomes Cancer 26, 247–252 (1999). [PubMed]
  • Samowitz W. S. et al. . APC mutations and other genetic and epigenetic changes in colon cancer. Mol. Cancer Res. 5, 165–170 (2007). [PubMed]
  • Christie M. et al. . Different APC genotypes in proximal and distal sporadic colorectal cancers suggest distinct WNT/beta-catenin signalling thresholds for tumourigenesis. Oncogene 32, 4675–4682 (2013). [PMC free article] [PubMed]
  • Vilar E. & Gruber S. B. Microsatellite instability in colorectal cancer-the stable evidence. Nat. Rev. 7, 153–162 (2010). [PMC free article] [PubMed]
  • Boland C. R. & Goel A. Microsatellite instability in colorectal cancer. Gastroenterology 138, 2073–2087 e2073 (2010). [PMC free article] [PubMed]
  • Jones S. et al. . Comparative lesion sequencing provides insights into tumor evolution. Proc. Natl Acad. Sci. USA 105, 4283–4288 (2008). [PubMed]
  • Klaus A. & Birchmeier W. Wnt signalling and its impact on development and cancer. Nat. Rev. Cancer 8, 387–398 (2008). [PubMed]
  • Fodde R., Smits R. & Clevers H. APC, signal transduction and genetic instability in colorectal cancer. Nat. Rev. Cancer 1, 55–67 (2001). [PubMed]
  • Albuquerque C. et al. . The 'just-right' signaling model: APC somatic mutations are selected based on a specific level of activation of the beta-catenin signaling cascade. Hum. Mol. Genet. 11, 1549–1560 (2002). [PubMed]
  • Herbst A. et al. . Comprehensive analysis of beta-catenin target genes in colorectal carcinoma cell lines with deregulated Wnt/beta-catenin signaling. BMC Genomics 15, 74 (2014). [PMC free article] [PubMed]
  • Schell M. J. et al. . A composite gene expression signature optimizes prediction of colorectal cancer metastasis and outcome. Clin. Cancer Res. 22, 734–745 (2016). [PMC free article] [PubMed]
  • Thiel A. & Ristimaki A. Toward a molecular classification of colorectal cancer: the role of BRAF. Front. Oncol. 3, 281 (2013). [PMC free article] [PubMed]
  • Samowitz W. S. et al. . Poor survival associated with the BRAF V600E mutation in microsatellite-stable colon cancers. Cancer Res. 65, 6063–6069 (2005). [PubMed]
  • Farina-Sarasqueta A. et al. . The BRAF V600E mutation is an independent prognostic factor for survival in stage II and stage III colon cancer patients. Ann. Oncol. 21, 2396–2402 (2010). [PubMed]
  • Roth A. D. et al. . Prognostic role of KRAS and BRAF in stage II and III resected colon cancer: results of the translational study on the PETACC-3, EORTC 40993, SAKK 60-00 trial. J. Clin. Oncol. 28, 466–474 (2010). [PubMed]
  • Ogino S. et al. . Predictive and prognostic roles of BRAF mutation in stage III colon cancer: results from intergroup trial CALGB 89803. Clin. Cancer Res. 18, 890–900 (2012). [PMC free article] [PubMed]
  • Douillard J. Y. et al. . Panitumumab-FOLFOX4 treatment and RAS mutations in colorectal cancer. N. Engl. J. Med. 369, 1023–1034 (2013). [PubMed]
  • Guinney J. et al. . The consensus molecular subtypes of colorectal cancer. Nat. Med. 21, 1350–1356 (2015). [PMC free article] [PubMed]
  • Conlin A., Smith G., Carey F. A., Wolf C. R. & Steele R. J. The prognostic significance of K-ras, p53, and APC mutations in colorectal carcinoma. Gut 54, 1283–1286 (2005). [PMC free article] [PubMed]
  • Rowan A. J. et al. . APC mutations in sporadic colorectal tumors: a mutational ‘hotspot' and interdependence of the ‘two hits'. Proc. Natl Acad. Sci. USA 97, 3352–3357 (2000). [PubMed]
  • Crabtree M. et al. . Refining the relation between 'first hits' and 'second hits' at the APC locus: the 'loose fit' model and evidence for differences in somatic mutation spectra among patients. Oncogene 22, 4257–4265 (2003). [PubMed]
  • Therman E., Susman B. & Denniston C. The nonrandom participation of human acrocentric chromosomes in Robertsonian translocations. Ann. Hum. Genet. 53, 49–65 (1989). [PubMed]
  • Chen J. et al. . Molecular analysis of APC promoter methylation and protein expression in colorectal cancer metastasis. Carcinogenesis 26, 37–43 (2005). [PubMed]
  • Hiltunen M. O. et al. . Hypermethylation of the APC (adenomatous polyposis coli) gene promoter region in human colorectal carcinoma. Int. J. Cancer 70, 644–648 (1997). [PubMed]
  • Sakamoto Y., Kitazawa R., Maeda S. & Kitazawa S. Methylation of CpG loci in 5'-flanking region alters steady-state expression of adenomatous polyposis coli gene in colon cancer cell lines. J. Cell Biochem. 80, 415–423 (2001). [PubMed]
  • Segditsas S. et al. . Promoter hypermethylation leads to decreased APC mRNA expression in familial polyposis and sporadic colorectal tumours, but does not substitute for truncating mutations. Exp. Mol. Pathol. 85, 201–206 (2008). [PubMed]
  • Guo Y. et al. . Association of aberrant DNA methylation in Apc(min/+) mice with the epithelial-mesenchymal transition and Wnt/beta-catenin pathways: genome-wide analysis using MeDIP-seq. Cell Biosci. 5, 24 (2015). [PMC free article] [PubMed]
  • Rad R. et al. . A genetic progression model of Braf(V600E)-induced intestinal tumorigenesis reveals targets for therapeutic intervention. Cancer Cell 24, 15–29 (2013). [PMC free article] [PubMed]
  • Jaiswal A. S., Armas M. L., Izumi T., Strauss P. R. & Narayan S. Adenomatous polyposis coli interacts with flap endonuclease 1 to block its nuclear entry and function. Neoplasia 14, 495–508 (2012). [PMC free article] [PubMed]
  • Jaiswal A. S., Balusu R., Armas M. L., Kundu C. N. & Narayan S. Mechanism of adenomatous polyposis coli (APC)-mediated blockage of long-patch base excision repair. Biochemistry 45, 15903–15914 (2006). [PMC free article] [PubMed]
  • Jaiswal A. S. & Narayan S. Assembly of the base excision repair complex on abasic DNA and role of adenomatous polyposis coli on its functional activity. Biochemistry 50, 1901–1909 (2011). [PMC free article] [PubMed]
  • Jaiswal A. S. et al. . Adenomatous polyposis coli-mediated accumulation of abasic DNA lesions lead to cigarette smoke condensate-induced neoplastic transformation of normal breast epithelial cells. Neoplasia 15, 454–460 (2013). [PMC free article] [PubMed]
  • Narayan S., Jaiswal A. S. & Balusu R. Tumor suppressor APC blocks DNA polymerase beta-dependent strand displacement synthesis during long patch but not short patch base excision repair and increases sensitivity to methylmethane sulfonate. J. Biol. Chem. 280, 6942–6949 (2005). [PubMed]
  • Kim R. et al. . Co-evolution of somatic variation in primary and metastatic colorectal cancer may expand biopsy indications in the molecular era. PLoS ONE 10, e0126670 (2015). [PMC free article] [PubMed]
  • Li H. & Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009). [PMC free article] [PubMed]
  • DePristo M. A. et al. . A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498 (2011). [PMC free article] [PubMed]
  • Wang K., Li M. & Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 38, e164 (2010). [PMC free article] [PubMed]
  • Li H. et al. . The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009). [PMC free article] [PubMed]
  • Giannakis M. et al. . RNF43 is frequently mutated in colorectal and endometrial cancers. Nat. Genet. 46, 1264–1266 (2014). [PMC free article] [PubMed]
  • Koo B. K. et al. . Tumour suppressor RNF43 is a stem-cell E3 ligase that induces endocytosis of Wnt receptors. Nature 488, 665–669 (2012). [PubMed]
  • Loregger A. et al. . The E3 ligase RNF43 inhibits Wnt signaling downstream of mutated beta-catenin by sequestering TCF4 to the nuclear membrane. Sci. Signal. 8, ra90 (2015). [PubMed]
  • Zebisch M. & Jones E. Y. ZNRF3/RNF43--A direct linkage of extracellular recognition and E3 ligase activity to modulate cell surface signalling. Prog. Biophys. Mol. Biol. 118, 112–118 (2015). [PubMed]
  • Curia M. C. et al. . Transcript dosage effect in familial adenomatous polyposis: model offered by two kindreds with exon 9 APC gene mutations. Hum. Mutat. 11, 197–201 (1998). [PubMed]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group