Single Nucleotide Polymorphisms in the Wnt and BMP Pathways and Colorectal Cancer Risk in a Spanish Cohort

Abstract

Background

Colorectal cancer (CRC) is considered a complex disease, and thus the majority of the genetic susceptibility is thought to lie in the form of low-penetrance variants following a polygenic model of inheritance. Candidate-gene studies have so far been one of the basic approaches taken to identify these susceptibility variants. The consistent involvement of some signaling routes in carcinogenesis provided support for pathway-based studies as a natural strategy to select genes that could potentially harbour new susceptibility loci.

Methodology/Principal Findings

We selected two main carcinogenesis-related pathways: Wnt and BMP, in order to screen the implicated genes for new risk variants. We then conducted a case-control association study in 933 CRC cases and 969 controls based on coding and regulatory SNPs. We also included rs4444235 and rs9929218, which did not fulfill our selection criteria but belonged to two genes in the BMP pathway and had consistently been linked to CRC in previous studies. Neither allelic, nor genotypic or haplotypic analyses showed any signs of association between the 37 screened variants and CRC risk. Adjustments for sex and age, and stratified analysis between sporadic and control groups did not yield any positive results either.

Conclusions/Significance

Despite the relevance of both pathways in the pathogenesis of the disease, and the fact that this is indeed the first study that considers these pathways as a candidate-gene selection approach, our study does not present any evidence of the presence of low-penetrance variants for the selected markers in any of the considered genes in our cohort.

About the Authors:

Ceres Fernández-Rozadilla

Affiliation: Galician Public Foundation of Genomic Medicine (FPGMX), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Genomics Medicine Group, Hospital Clínico, Santiago de Compostela, University of Santiago de Compostela, Galicia, Spain

Luisa de Castro

Affiliation: Gastroenterology Department, Hospital Meixoeiro, Vigo, Galicia, Spain

Juan Clofent

Affiliation: Gastroenterology Department, Hospital La Fe, Valencia, Spain

Alejandro Brea-Fernández

Affiliation: Galician Public Foundation of Genomic Medicine (FPGMX), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Genomics Medicine Group, Hospital Clínico, Santiago de Compostela, University of Santiago de Compostela, Galicia, Spain

Xavier Bessa

Affiliation: Gastroenterology Department, Hospital del Mar, Institut Municipal d'Investigació Médica (IMIM), Pompeu Fabra University, Barcelona, Catalonia, Spain

Anna Abulí

Affiliation: Gastroenterology Department, Hospital del Mar, Institut Municipal d'Investigació Médica (IMIM), Pompeu Fabra University, Barcelona, Catalonia, Spain

Montserrat Andreu

Affiliation: Gastroenterology Department, Hospital del Mar, Institut Municipal d'Investigació Médica (IMIM), Pompeu Fabra University, Barcelona, Catalonia, Spain

Rodrigo Jover

Affiliation: Unidad de Gastroenterología, Hospital General Universitario de Alicante, Alicante, Spain

Rosa Xicola

Affiliation: Section of Digestive Diseases and Nutrition, University of Illinois at Chicago, Chicago, Illinois, United States of America

Xavier Llor

Affiliation: Section of Digestive Diseases and Nutrition, University of Illinois at Chicago, Chicago, Illinois, United States of America

Antoni Castells

Affiliation: Department of Gastroenterology, Hospital Clínic, CIBERehd, IDIBAPS, University of Barcelona, Barcelona, Catalonia, Spain

Sergi Castellví-Bel

Affiliation: Department of Gastroenterology, Hospital Clínic, CIBERehd, IDIBAPS, University of Barcelona, Barcelona, Catalonia, Spain

Angel Carracedo

Affiliation: Galician Public Foundation of Genomic Medicine (FPGMX), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Genomics Medicine Group, Hospital Clínico, Santiago de Compostela, University of Santiago de Compostela, Galicia, Spain

Clara Ruiz-Ponte

* E-mail: [email protected]

Affiliation: Galician Public Foundation of Genomic Medicine (FPGMX), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Genomics Medicine Group, Hospital Clínico, Santiago de Compostela, University of Santiago de Compostela, Galicia, Spain

for the Gastrointestinal Oncology Group of the Spanish Gastroenterological Association

¶Members of the Gastrointestinal Oncology Group of the Spanish Gastroenterological Association are listed inNote S1

Introduction

Colorectal cancer (CRC) is one of the main forms of cancer, being the second most frequent neoplasm in both sexes and one of the most important morbidity causes in the western world [1]. The genetic contribution to CRC has been estimated to be around 35% by extensive twin studies [2]. However, highly penetrant variants, that cause mendelian predisposition syndromes, account only for, at most, 5% of the disease cases [3]. The remaining genetic susceptibility is thought to follow a polygenic model, with an interplay of multiple low-penetrance allelic variants appearing in high frequency in the general population, and each conferring a modest effect on disease risk [4], [5].

Candidate-gene studies have been one of the most commonly used tools in the screening for new variants affecting CRC risk. Gene selection in these studies is mainly based on the functional implications of a possible association, and thus genes selected have either been chosen because of the previous presence of other high/low risk alleles [6], or their participation in a pathway implicated in the pathogenesis of the disease [7]. Candidate-gene studies can be performed by either direct approaches, where the variants genotyped are presumed to be the underlying cause of the disease because of their location (variants in exonic or regulatory regions), or by indirect approaches, where tag SNPs take advantage of the linkage disequilibrium properties of the human genome to try and screen the most of the variability in a given gene.

This latter approach has also allowed, together with the development of high-throughput technologies, the implementation of new hypothesis-free approaches (in opposition with hypothesis-based candidate-gene approaches), covering the majority of the genome (genome-wide association studies or GWAS). This implementation has successfully led to the identification of some new susceptibility loci [8]–[14], including rs4444235 and rs9929218, that fall within reach of two genes belonging to the BMP pathway. Nevertheless, these have been found to predict only a small proportion of the disease susceptibility, with the remaining yet to be discovered [15].

We hence aimed to find such susceptibility variants through a candidate-gene approach screening a selected number of variants within two cellular pathways that have consistently been linked to CRC tumorogenesis: the Wnt and the BMP signaling pathways [16], [17].

The Wnt pathway contains genes that have for long been known to be responsible of some hereditary CRC syndromes, such as APC and familial adenomatous polyposis [18]. Moreover, somatic alterations in APC are found in almost 80% of the sporadic colorectal cancers, and Wnt signaling activation is involved in the best part of sporadic colorectal carcinomas [19]. On the other hand, the BMP pathway acts as positive regulator of some of the Wnt proteins [17], and the tumor suppressive role of this signaling pathway in the pathogenesis of CRC and other cancers is well established [20], [21]. Besides, mutations in two of its genes, SMAD4 and BMPR1A, are responsible for juvenile polyposis syndrome, another hereditary CRC condition [22]. Considering all this information, we thought it would be interesting to screen some of the genetic variability within these pathways for any evidence of new CRC related variants that could explain at least part of the missing heritability. Our approach was mainly functional, for only SNPs within exonic or cis-regulatory sequences (5′ and 3′ unstranslated regions) were selected to analyse their relationship with CRC susceptibility.

Results and Discussion

Following our pathway-based candidate-gene selection method, we performed our study in a total of 45 SNPs that were in either exonic or regulatory regions, in an overall of 21 genes from both the Wnt and BMP pathways. Details of SNP features and association values for the 37 SNPs that successfully passed quality control criteria are shown on Table 1. None of the screened SNPs were significantly associated with an altered risk of CRC, considering odds-ratios and related p values for allelic and genotypic tests (trend, dominant and recessive). Logistic regression for age and sex adjustment was performed, although it did not improve p value results. Haplotype analysis results were consistent in both Unphased and Haploview, and did not show any signs of positive associations either for any of the 8 genes for which this analysis was performed (AXIN1, HDAC9, BMP4, DACT1, CDH3, CDH1, BTRC, and APC), (Figure S1). Stratification analysis comparing sporadic and familial cases was also implemented, but it did not provide any evidence of differences in susceptibilities between the groups that could be a sign of any specific associations within either of the groups (Table 2).

[Figure omitted. See PDF.]

Table 1. Description of the 37 SNPs that passed quality control criteria and their associated p values.

https://doi.org/10.1371/journal.pone.0012673.t001

[Figure omitted. See PDF.]

Table 2. Association values for stratified analysis in familial and sporadic CRC groups.

https://doi.org/10.1371/journal.pone.0012673.t002

Thus, our strategy has not managed to detect any new susceptibility loci for CRC risk.

Pathway-based expectations have proved to be quite discouraging in the literature as well, for strong candidate pathways, such as DNA-repair ones, surprisingly failed too in identifying any new risk variants [7], [23]–[24]. In addition to this, most of the genetic variants that have been found to be associated with disease are located in intergenic regions, with potential functions that are yet unknown.

Still, in light of the recent discoveries that followed up the analysis of genome-wide data, both Wnt and BMP have earned a renewed fame. The susceptibility locus found on 8q24 (rs6983267) has been linked to an enhanced Wnt signaling through its interaction with TCF4 [25], [26], and a meta-analysis conducted on a series of GWAS data succeeded in associating two variants in the BMP4 and CDH1 gene regions with the disease (rs4444235 and rs9929218, respectively)[8].

Even though this is actually the first association study that considers the pathways as a whole for gene selection, some of the genes included in our analysis (i.e APC, CCND1, CDH1 and TCF7) had already been screened for risk alleles [6], [27]–[30]. It is quite remarkable that there has been a growing debate over some of these loci, specially the p.V1822D variant in APC (rs459552). This missense change is widely documented in the literature, with some studies defending it as neutral (this study and others)[31], and some conferring its minor allele a protective effect [6], [28]. Lack of appropriate study power, resultant from insufficient number of samples has been a major problem in many of these studies and thus most of them have not provided very convincing results [32].

Although our study had over 80% power to detect OR as low as 1.21 with minor allele frequencies of 0.30 (57% of our SNPs), and 1.24 for MAFs down to 0.2 (78% of the SNPs), assuming a log-additive model and α = 0.05, we were unable to detect any positive associations suggesting the presence of any new CRC susceptibility variants. Nevertheless, it is quite remarkable that, albeit our failure to replicate the associations for the BMP4 and CDH1 SNPs, this is the first study that investigates any of the so-called 10 new GWAS-discovered susceptibility loci in a Southern-European population.

Despite our negative results, we must consider that we did not whatsoever comprehensively cover all possible low-penetrance variants within the selected genes. This is mainly due to the fact that our strategy was purely functional, selecting the variants that were a priori good candidates to be directly associated with the disease. This indeed may constitute a limitation in the study, for most of the genetic variation within the loci was not investigated. Thus, we believe further efforts should be made to screen a wider variety of loci within these pathways, specially considering the previous positive associations described so far for both Wnt and BMP-related genes.

Pondering the potential odds ratios of the variants described so far (1.11, CI 1.08–1.15 and 0.91, CI 0.89–0.94 for rs4444235 and rs9929218, respectively), we assume larger cohorts may be required to detect such subtle effects. On the other hand, when considering candidate-gene approaches, it would also be useful to meta-analyse previous studies and pull the information across of them altogether in the search of evidences of potential new pathways linked to the pathogenesis of the disease.

Materials and Methods

Study populations

Subjects were 933 CRC patients and 969 controls that belonged to the EPICOLON project, a prospective, multicentre, population-based epidemiology survey studying the incidence and features of familial and sporadic CRC in the Spanish population [33]. Cases were selected across 11 hospitals in Spain as all patients with a de-novo histologically confirmed diagnosis of colorectal adenocarcinoma and who attended 11 community hospitals across Spain between November 2006 and December 2007. Patients in whom CRC developed in the context of familial adenomatous polyposis or inflammatory bowel disease, and cases where patients or family refused to participate in the study were excluded. Demographic, clinical and tumour-related characteristics of probands, as well as a detailed family history were obtained using a pre-established questionnaire, and registered in a single database. Of these, 592 (63%) were male and 341 (37%) female. Median age for cases was 73 (range 26–95), whereas mean was 71(SD±10.7). Hospital-based controls were recruited together with cases and were confirmed to have no cancer or prior history of neoplasm, and no family history of CRC. All controls were randomly selected and matched with cases for sex and age (±5 years) in a 1∶1 ratio. Both cases and controls were of European ancestry and from Spain.

Ethics statement

The study was approved by the “Comité Ético de Investigación Clínica de Galicia”, and each of the institutional review boards of the hospitals where samples were collected (“Ethics Committee of the Hospital Clínic-Barcelona”, “Ethics Committee of the Hospital del Mar-Barcelona”, “Ethics Committee of the Hospital German Trias i Pujol-Barcelona”, “Ethics Committee of the Hospital Sant Pau-Barcelona,” “Ethics Committee of the Hospital Universitari Arnau de Vilanova-Lleida”, “Ethics Committee of the Hospital General-Alicante”, “Ethics Committee of the Hospital de Donosti”, “Ethics Committee of the Hospital General de Asturias-Oviedo”, “Ethics Committee of the Hospital Clinico-Zaragoza”, “Ethics Committee of the Hospital de Calahorra-La Rioja”, “Ethics Committee of the Hospital Meixoeiro-Vigo”). All samples were obtained with written informed consent reviewed by the ethical board of the corresponding hospital.

DNA extraction

DNA was obtained from frozen peripheral blood; extraction was performed in a CHEMAGEN robot (Chemagen Biopolymer-Technologie AG, Baesweiler, Germany) in accordance with the manufacturer's instructions, at the Galician Public Fundation of Genomic Medicine in Santiago de Compostela. Cases and controls were extracted in mixed batches to avoid any kind of bias.

Candidate-gene selection

Both Wnt and BMP pathways were initially selected after the findings of Nishanian et al. [34], who demonstrated the interaction between these two pathways. Both pathways were thoroughly investigated through the Cancer Genome Anatomy Project site [35], but we failed to find any information regarding the BMP pathway in either this or other web browsers. For that reason, Wnt genes were selected by browsing the pathway through Biocarta [36], whereas BMP genes had to be strictly selected from previous literature [17], [34]. Forty-one genes were finally selected to be included in the analysis.

SNP selection and genotyping

SNP selection criteria only considered functional markers with minor allele frequencies above 0.05 and at least two independent validation criteria as established in dbSNP [37]. This included all exonic variants selected with Pupasuite [38] and gene-regulatory regions in cis (5′or 3′ UTR ends), as defined by the FESD web browser [39]. 5′UTR variants were only included when they complied to the abovementioned criteria and were presumed to be in the potential binding site of a known transctiptional binding factor. 3′ UTR variants were included because of their potential relationship with miRNA binding regions [40]. Because some of the selected genes had no SNPs of such these kinds in any of the three browsers at the time of SNP selection, they ultimately had to be dropped out of the study. Finally, 43 SNPs were chosen within 21 genes to be screened as potential direct modifiers of CRC susceptibility (Table 3).

[Figure omitted. See PDF.]

Table 3. Description of all genes selected from both pathways and SNPs screened within each of them.

https://doi.org/10.1371/journal.pone.0012673.t003

rs4444235 and rs9929218 are two variants lying in the near-by and intronic regions of BMP4 and CDH1, respectively, that have been recently reported to be associated with the disease [8]. Considering that the SNPs that we had chosen within these two genes were not good taggers for these two variants (r-squared values were 0.6 for the SNPs in BMP4, and 0.02 for those in CHD1) (Figure 1), we decided to include them in our study as well, although they did not fulfill our selection criteria, making the total number of interrogated SNPs rise to 45.

[Figure omitted. See PDF.]

Figure 1. Linkage disequilibrium blocks for the BMP4 and CDH1 genes.

R-squared relationships between SNP pairs: A. rs4444235-rs17563 in BMP4 and B. rs9929218-rs1801552 in CDH1.

https://doi.org/10.1371/journal.pone.0012673.g001

Genotyping was performed with the MassARRAY (Sequenom Inc., San Diego, USA) technology at the Santiago de Compostela node of the Spanish Genotyping Center. Calling of genotypes was done with Sequenom Typer v4.0 software using all the data from the study simultaneously.

Statistical analyses

Quality control was performed, first by excluding both SNPs and samples with genotype success rates below 95%, with the help of the Genotyping Data Filter (GDF) [41]. Genotypic distributions for all SNPs in controls were consistent with Hardy-Weinberg equilibrium as assessed using a X2 test (1df). All p-values obtained were ≥0.05, thereby excluding the possibility of genotyping artifacts (data not shown). Population stratification was assessed with Structure v2.2 [42]. Briefly, the posibility of different scenarios was tested assuming a different number of underlying populations (k ranging from 1 to 4), allowing for a large number of iterations (25 K in the burn-in period followed by 500 K repetitions). The mean log likelihood was estimated for the data for a given k (referred to as L(K)) in each run. We as well performed multiple runs for each value of k computing the overall mean L(K) and its standard deviation. All results seemed to be concordant with the original assumption of a single existing population. Moreover, additional procedures for better confounding variable visualization were undertaken by means of a Principal Component Analysis (PCA) using the EIGENSOFT tool smartpca [43], although number of markers was very low. No differences were found of population stratification between cases and controls for either STRUCTURE or the first 10 components of the PCA analysis (Figure S2). After quality control 1746 samples (854 cases and 892 controls) and 37 SNPs remained for further analyses.

Association tests were performed by chi-squared tests for every single SNP and haplotypes where possible with both Haploview v4.0 [44] and Unphased [45]. In short, LD patterns across genes for which more than one SNP was genotyped were checked in Haploview and tested for association using Unphased (to check in any of the haplotypes was associated) and Haploview (to see which of the haplotypes was associated). Genotypic association tests, logistic regression analysis for sex and age adjustment, and stratified analysis between sporadic and familial groups were estimated with PLINK v1.03 [46]. OR and 95% confidence intervals were calculated for each statistic, and to address the issue of multiple-testing, permutation tests and the Bonferroni correction were used. Study power was estimated with CATS software [47].

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Haplotype structure and analysis for the 8 genes for which more than one SNP was genotyped. The table shows association values for each SNP generated by Haploview.

https://doi.org/10.1371/journal.pone.0012673.s001

(3.40 MB TIF)

Figure S2.

Principal component analysis plot for the first vs. second component, comparing our case and control populations.

https://doi.org/10.1371/journal.pone.0012673.s002

(0.96 MB TIF)

Note S1.

https://doi.org/10.1371/journal.pone.0012673.s003

(0.03 MB DOC)

Acknowledgments

We are sincerely grateful to all patients participating in this study who were recruited in 11 Spanish hospitals as part of the EPICOLON project. We thank Maria Magdalena-Castro (MMC), Olga Lortes (OL) and Eva Fernández (EF) for their excellent technical assistance.

Author Contributions

Conceived and designed the experiments: CRP. Performed the experiments: CFR. Analyzed the data: CFR ABF AA RMX SCB A. Carracedo. Contributed reagents/materials/analysis tools: LdC JC XB MA RJ XL A. Castells. Wrote the paper: CFR. Revised the manuscript and approved the final version: A. Castells SCB A. Carracedo.

Citation: Fernández-Rozadilla C, de Castro L, Clofent J, Brea-Fernández A, Bessa X, Abulí A, et al. (2010) Single Nucleotide Polymorphisms in the Wnt and BMP Pathways and Colorectal Cancer Risk in a Spanish Cohort. PLoS ONE5(9): e12673. https://doi.org/10.1371/journal.pone.0012673

References

1. Ferlay J, Bray F, Pisani P, Parkin DM (2004) GLOBOCAN 2002: Cancer Incidence, Mortality and Prevalence Worldwide IARC Cancer Base No. 5. version 2.0.J. FerlayF. BrayP. PisaniDM Parkin2004GLOBOCAN 2002: Cancer Incidence, Mortality and Prevalence Worldwide IARC Cancer Base No. 5. version 2.0.

2. Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, et al. (2000) Environmental and heritable factors in the causation of cancer–analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med 343: 78–85.P. LichtensteinNV HolmPK VerkasaloA. IliadouJ. Kaprio2000Environmental and heritable factors in the causation of cancer–analyses of cohorts of twins from Sweden, Denmark, and Finland.N Engl J Med3437885

3. de la Chapelle A (2004) Genetic predisposition to colorectal cancer. Nat Rev Cancer 4: 769–80.A. de la Chapelle2004Genetic predisposition to colorectal cancer.Nat Rev Cancer476980

4. Castells A, Castellví-Bel S, Balaguer F (2009) Concepts in familial colorectal cancer: where do we stand and what is the future? Gastroenterology 137: 404–409.A. CastellsS. Castellví-BelF. Balaguer2009Concepts in familial colorectal cancer: where do we stand and what is the future?Gastroenterology137404409

5. Houlston RS, Peto J (2004) The search for low-penetrance cancer susceptibility alleles. Oncogene 23: 6471–6.RS HoulstonJ. Peto2004The search for low-penetrance cancer susceptibility alleles.Oncogene2364716

6. Chen SP, Tsai ST, Jao SW, Huang YL, Chao YC, et al. (2006) Single nucleotide polymorphisms of the APC gene and colorectal cancer risk: a case-control study in Taiwan. BMC Cancer 6: 83.SP ChenST TsaiSW JaoYL HuangYC Chao2006Single nucleotide polymorphisms of the APC gene and colorectal cancer risk: a case-control study in Taiwan.BMC Cancer683

7. Naccarati A, Pardini B, Hemminki K, Vodicka P (2007) Sporadic colorectal cancer and individual susceptibility: a review of the association studies investigating the role of DNA repair genetic polymorphisms. Mutat Res 635: 118–45.A. NaccaratiB. PardiniK. HemminkiP. Vodicka2007Sporadic colorectal cancer and individual susceptibility: a review of the association studies investigating the role of DNA repair genetic polymorphisms.Mutat Res63511845

8. Houlston RS, Webb E, Broderick P, Pittman AM, Di Bernardo MC, et al. (2008) Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer. Nat Genet 40: 1426–35.RS HoulstonE. WebbP. BroderickAM PittmanMC Di Bernardo2008Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer.Nat Genet40142635

9. Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, et al. (2007) A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nat Genet 39: 984–8.I. TomlinsonE. WebbL. Carvajal-CarmonaP. BroderickZ. Kemp2007A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21.Nat Genet399848

10. Zanke BW, Greenwood CM, Rangrej J, Kustra R, Tenesa A, et al. (2007) Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24. Nat Genet 39: 989–94.BW ZankeCM GreenwoodJ. RangrejR. KustraA. Tenesa2007Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24.Nat Genet3998994

11. Broderick P, Carvajal-Carmona L, Pittman AM, Webb E, Howarth K, et al. (2007) A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk. Nat Genet 39: 1315–7.P. BroderickL. Carvajal-CarmonaAM PittmanE. WebbK. Howarth2007A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk.Nat Genet3913157

12. Jaeger E, Webb E, Howarth K, Carvajal-Carmona L, Rowan A, et al. (2008) Common genetic variants at the CRAC1 (HMPS) locus on chromosome 15q13.3 influence colorectal cancer risk. Nat Genet 40: 26–8.E. JaegerE. WebbK. HowarthL. Carvajal-CarmonaA. Rowan2008Common genetic variants at the CRAC1 (HMPS) locus on chromosome 15q13.3 influence colorectal cancer risk.Nat Genet40268

13. Tomlinson IP, Webb E, Carvajal-Carmona L, Broderick P, Howarth K, et al. (2008) A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3. Nat Genet 40(5): 623–30.IP TomlinsonE. WebbL. Carvajal-CarmonaP. BroderickK. Howarth2008A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3.Nat Genet40562330

14. Tenesa A, Farrington SM, Prendergast JG, Porteous ME, Walker M, et al. (2008) Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nat Genet 40(5): 631–7.A. TenesaSM FarringtonJG PrendergastME PorteousM. Walker2008Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21.Nat Genet4056317

15. Cazier JB, Tomlinson I (2010) General lessons from large-scale studies to identify human cancer predisposition genes. J Pathol 220(2): 255–62.JB CazierI. Tomlinson2010General lessons from large-scale studies to identify human cancer predisposition genes.J Pathol220225562

16. Segditsas S, Tomlinson I (2006) Colorectal cancer and genetic alterations in the Wnt pathway. Oncogene 25: 7531–7.S. SegditsasI. Tomlinson2006Colorectal cancer and genetic alterations in the Wnt pathway.Oncogene2575317

17. Huang HC, Klein PS (2004) Interactions between BMP and Wnt signaling pathways in mammalian cancers. Cancer Biol Ther 3: 676–8.HC HuangPS Klein2004Interactions between BMP and Wnt signaling pathways in mammalian cancers.Cancer Biol Ther36768

18. Half E, Bercovich D, Rozen P (2009) Familial adenomatous polyposis. Orphanet J Rare Dis 4: 22.E. HalfD. BercovichP. Rozen2009Familial adenomatous polyposis.Orphanet J Rare Dis422

19. Rowan AJ, Lamlum H, Ilyas M, Wheeler J, Straub J, et al. (2000) APC mutations in sporadic colorectal tumors: A mutational “hotspot” and interdependence of the “two hits”. Proc Natl Acad Sci U S A 97: 3352–7.AJ RowanH. LamlumM. IlyasJ. WheelerJ. Straub2000APC mutations in sporadic colorectal tumors: A mutational “hotspot” and interdependence of the “two hits”.Proc Natl Acad Sci U S A9733527

20. Deng H, Ravikumar TS, Yang WL (2009) Overexpression of bone morphogenetic protein 4 enhances the invasiveness of Smad4-deficient human colorectal cancer cells. Cancer Lett 281: 220–31.H. DengTS RavikumarWL Yang2009Overexpression of bone morphogenetic protein 4 enhances the invasiveness of Smad4-deficient human colorectal cancer cells.Cancer Lett28122031

21. Deng H, Makizumi R, Ravikumar TS, Dong H, Yang W, et al. (2007) Bone morphogenetic protein-4 is overexpressed in colonic adenocarcinomas and promotes migration and invasion of HCT116 cells. Exp Cell Res 313: 1033–44.H. DengR. MakizumiTS RavikumarH. DongW. Yang2007Bone morphogenetic protein-4 is overexpressed in colonic adenocarcinomas and promotes migration and invasion of HCT116 cells.Exp Cell Res313103344

22. Chen HM, Fang JY (2009) Genetics of the hamartomatous polyposis syndromes: a molecular review. Int J Colorectal Dis 24(8): 865–74.HM ChenJY Fang2009Genetics of the hamartomatous polyposis syndromes: a molecular review.Int J Colorectal Dis24886574

23. Tao H, Shinmura K, Suzuki M, Kono S, Mibu R, et al. (2008) Association between genetic polymorphisms of the base excision repair gene MUTYH and increased colorectal cancer risk in a Japanese population. Cancer Sci 99: 355–60.H. TaoK. ShinmuraM. SuzukiS. KonoR. Mibu2008Association between genetic polymorphisms of the base excision repair gene MUTYH and increased colorectal cancer risk in a Japanese population.Cancer Sci9935560

24. Schafmayer C, Buch S, Egberts JH, Franke A, Brosch M, et al. (2007) Genetic investigation of DNA-repair pathway genes PMS2, MLH1, MSH2, MSH6, MUTYH, OGG1 and MTH1 in sporadic colon cancer. Int J Cancer 121: 555–8.C. SchafmayerS. BuchJH EgbertsA. FrankeM. Brosch2007Genetic investigation of DNA-repair pathway genes PMS2, MLH1, MSH2, MSH6, MUTYH, OGG1 and MTH1 in sporadic colon cancer.Int J Cancer1215558

25. Pomerantz MM, Ahmadiyeh N, Jia L, Herman P, Verzi MP, et al. (2009) The 8q24 cancer risk variant rs6983267 shows long-range interaction with MYC in colorectal cancer. Nat Genet 41: 882–4.MM PomerantzN. AhmadiyehL. JiaP. HermanMP Verzi2009The 8q24 cancer risk variant rs6983267 shows long-range interaction with MYC in colorectal cancer.Nat Genet418824

26. Tuupanen S, Turunen M, Lehtonen R, Hallikas O, Vanharanta S, et al. (2009) The common colorectal cancer predisposition SNP rs6983267 at chromosome 8q24 confers potential to enhanced Wnt signaling. Nat Genet 41: 885–90.S. TuupanenM. TurunenR. LehtonenO. HallikasS. Vanharanta2009The common colorectal cancer predisposition SNP rs6983267 at chromosome 8q24 confers potential to enhanced Wnt signaling.Nat Genet4188590

27. Picelli S, Zajac P, Zhou XL, Edler D, Lenander C, et al. (2010) Common variants in human CRC genes as low-risk alleles. Eur J Cancer 46(6): 1041–8.S. PicelliP. ZajacXL ZhouD. EdlerC. Lenander2010Common variants in human CRC genes as low-risk alleles.Eur J Cancer46610418

28. Slattery ML, Samowitz W, Ballard L, Schaffer D, Leppert M, et al. (2001) A molecular variant of the APC gene at codon 1822: its association with diet, lifestyle, and risk of colon cancer. Cancer Res 61: 1000–4.ML SlatteryW. SamowitzL. BallardD. SchafferM. Leppert2001A molecular variant of the APC gene at codon 1822: its association with diet, lifestyle, and risk of colon cancer.Cancer Res6110004

29. Liu B, Zhang Y, Jin M, Ni Q, Liang X, et al. (2010) Association of selected polymorphisms of CCND1, p21, and caspase8 with colorectal cancer risk. Mol Carcinog 49: 75–84.B. LiuY. ZhangM. JinQ. NiX. Liang2010Association of selected polymorphisms of CCND1, p21, and caspase8 with colorectal cancer risk.Mol Carcinog497584

30. Hazra A, Fuchs CS, Chan AT, Giovannucci EL, Hunter DJ (2008) Association of the TCF7L2 polymorphism with colorectal cancer and adenoma risk. Cancer Causes Control 19: 975–80.A. HazraCS FuchsAT ChanEL GiovannucciDJ Hunter2008Association of the TCF7L2 polymorphism with colorectal cancer and adenoma risk.Cancer Causes Control1997580

31. Ruiz-Ponte C, Vega A, Conde R, Barros F, Carracedo A (2001) The Asp1822Val variant of the APC gene is a common polymorphism without clinical implications. J Med Genet 38: E33.C. Ruiz-PonteA. VegaR. CondeF. BarrosA. Carracedo2001The Asp1822Val variant of the APC gene is a common polymorphism without clinical implications.J Med Genet38E33

32. Kemp Z, Thirlwell C, Sieber O, Silver A, Tomlinson I (2004) An update on the genetics of colorectal cancer. Hum Mol Genet 13 Spec No 2: R177–85.Z. KempC. ThirlwellO. SieberA. SilverI. Tomlinson2004An update on the genetics of colorectal cancer.Hum Mol Genet 13 Spec No 2R17785

33. Pinol V, Castells A, Andreu M, Castellvi-Bel S, Alenda C, et al. (2005) Accuracy of revised Bethesda guidelines, microsatellite instability, and immunohistochemistry for the identification of patients with hereditary nonpolyposis colorectal cancer. JAMA 293: 1986–94.V. PinolA. CastellsM. AndreuS. Castellvi-BelC. Alenda2005Accuracy of revised Bethesda guidelines, microsatellite instability, and immunohistochemistry for the identification of patients with hereditary nonpolyposis colorectal cancer.JAMA293198694

34. Nishanian TG, Kim JS, Foxworth A, Waldman T (2004) Suppression of tumorigenesis and activation of Wnt signaling by bone morphogenetic protein 4 in human cancer cells. Cancer Biol Ther 3: 667–75.TG NishanianJS KimA. FoxworthT. Waldman2004Suppression of tumorigenesis and activation of Wnt signaling by bone morphogenetic protein 4 in human cancer cells.Cancer Biol Ther366775

35. The Cancer Genome Anatomy Project webpage (2006) The Cancer Genome Anatomy Project webpage2006http://cgap.nci.nih.gov/. http://cgap.nci.nih.gov/.

36. Biocarta (2006) Biocarta2006http://www.biocarta.com/genes/index.asp. http://www.biocarta.com/genes/index.asp.

37. 126.dbSNP (accessed 2007). http://www.ncbi.nlm.nih.gov/projects/SNP/. 126.dbSNP (accessed 2007). http://www.ncbi.nlm.nih.gov/projects/SNP/.

38. Conde L, Vaquerizas JM, Dopazo H, Arbiza L, Reumers J, et al. (2006) PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes. Nucleic Acids Res 34: W621–5.L. CondeJM VaquerizasH. DopazoL. ArbizaJ. Reumers2006PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes.Nucleic Acids Res34W6215

39. Kang HJ, Choi KO, Kim BD, Kim S, Kim YJ (2005) FESD: a Functional Element SNPs Database in human. Nucleic Acids Res 33: D518–22.HJ KangKO ChoiBD KimS. KimYJ Kim2005FESD: a Functional Element SNPs Database in human.Nucleic Acids Res33D51822

40. Jackson RJ, Standart N (2007) How do microRNAs regulate gene expression? Science. (367).re1 p.RJ JacksonN. Standart2007How do microRNAs regulate gene expression?Science367re1

41. Genotyping Data Filter Software (2009) Genotyping Data Filter Software2009http://bioinformatics.cesga.es/gdf/nav_input.php. http://bioinformatics.cesga.es/gdf/nav_input.php.

42. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–59.JK PritchardM. StephensP. Donnelly2000Inference of population structure using multilocus genotype data.Genetics15594559

43. Patterson N, Price AL, Reich D (2006) Population structure and Eigenanalysis. PLoS Genet 2(12).N. PattersonAL PriceD. Reich2006Population structure and Eigenanalysis.PLoS Genet212

44. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–5.JC BarrettB. FryJ. MallerMJ Daly2005Haploview: analysis and visualization of LD and haplotype maps.Bioinformatics212635

45. Dudbridge F (2003) Pedigree disequilibrium tests for multilocus haplotypes. Genet Epidemiol 25: 115–21.F. Dudbridge2003Pedigree disequilibrium tests for multilocus haplotypes.Genet Epidemiol2511521

46. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a toolset for whole-genome association and population-based linkage-analysis. American Journal of Human Genetics. 81. : 559–75.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a toolset for whole-genome association and population-based linkage-analysis.American Journal of Human Genetics8155975http://pngu.mgh.harvard.edu/~purcell/plink/. http://pngu.mgh.harvard.edu/~purcell/plink/.

47. Skol AD, Scott LJ, Abecasis GR, Boehnke M (2006) Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet 38: 209–13.AD SkolLJ ScottGR AbecasisM. Boehnke2006Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies.Nat Genet3820913