
Cancer Epidemiology Biomarkers & Prevention 17, 1499, June 1, 2008. doi: 10.1158/1055-9965.EPI-07-0235
© 2008 American Association for Cancer Research
Strong Association of Chromosome 1p12 Loci with Thyroid Cancer Susceptibility
Aida Baida1,2,
Mounaim Akdi1,
Eddy González-Flores1,
Pere Galofré3,
Ricard Marcos1,2 and
Antonia Velázquez1,2
1 Grup de Mutagènesi, Unitat de Genètica, Departament de Genètica i de Microbiologia, Edifici Cn, Universitat Autònoma de Barcelona, Bellaterra, Spain; 2 CIBER Epidemiologia y Salud Pública, ISCIII, Madrid, Spain; and 3 Servei de Medicina Nuclear, Hospital Josep Trueta, Girona, Spain
Requests for reprints: Antonia Velázquez, Departament de Genètica i de Microbiologia, Edifici Cn, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain. Phone: 34-93-5813111; Fax: 34-93-5812387. E-mail: antonia.velazquez{at}uab.es
 |
Abstract
|
|---|
Several genes directly related to thyroid cancer development have been described; nevertheless, the genetic pathways of this tumorigenesis process are unknown. Together with environmental factors, susceptibility genes could have an important role in thyroid cancer. Our previous studies suggest that the chromosome 1p12-13 is related to thyroid cancer incidence. Here, we extend the analysis with a case-control association study in a Spanish population. Thus, six single-nucleotide polymorphisms were genotyped, covering 2.4 Mb of the 1p12-13 region. A statistically significant association between thyroid cancer incidence and the rs2145418 and rs4658973 polymorphisms was found (P < 0.0001). No association was detected for the other four polymorphisms studied. The rs2145418 marker showed a significant odds ratio of 5.0 [95% confidence interval (95% CI), 2.85-8.83] and 9.2 (95% CI, 4.50-21.6) for heterozygous and homozygous G-variant alleles, respectively. For rs4658973, the odds ratios were 0.40 (95% CI, 0.26-0.62) and 0.07 (95% CI, 0.03-0.18) for heterozygous and homozygous G-variant alleles, respectively. These markers map into the 1p12 region, and no linkage disequilibrium was found between them, indicating an independent relation of these polymorphisms with thyroid cancer susceptibility. Our data provide evidence of a strong association of the chromosome 1p12 with thyroid cancer risk, and it is the first study describing susceptibility loci for thyroid cancer in this region. (Cancer Epidemiol Biomarkers Prev 2008;17(6):1499–504)
 |
Introduction
|
|---|
Thyroid cancer is the most frequent endocrine cancer with familial and sporadic forms, the latter being the most common type of thyroid cancer (1). Genetic factors directly related to thyroid tumorigenesis have been reported (2, 3), but the genetic pathways of such processes remain to be established (3, 4). Together with the genetic effects, the role of environmental factors in thyroid cancer development is well known, ionizing radiation being the best recognized risk factor for this pathology (5, 6). Nevertheless, the etiology of thyroid cancer is still unknown. In general, it is believed that susceptibility genetic factors are crucial in the genomics of cancer development; thus, association studies aiming to identify susceptibility genes for cancer are taking important relevance. This type of study has recently emerged for thyroid cancer (7-12), although no specific genetic factor for thyroid cancer susceptibility has yet been described. In a previous study, we have reported that the chromosome 1p12-13is related to thyroid cancer susceptibility (13). Several studies also indicate that this region undergoes many aberrations in different types of cancer (14-17), which would suggest a certain role in cancer incidence. Here, we report the localization of susceptibility factors to thyroid cancer in this region. We have carried out a case-control study in a Spanish population analyzing six single-nucleotide polymorphisms covering 2.4 Mb of the chromosome 1p12-13 region to localize susceptibility factors. We found that, in the studied population, the single-nucleotide polymorphisms rs2145418 and rs4658973, which are separated by 377 kb, have a strong association (P < 0.0001) with thyroid cancer.
 |
Materials and Methods
|
|---|
Subjects and DNA Isolation
In this study, two populations of unrelated Spanish subjects were genotyped. Blood samples were collected from 202 healthy individuals (117 women and 85 men; mean age, 41.2 ± 13.2 years) and from 227 thyroid cancer patients (167 women and 60 men; mean age, 43.5 ± 14.6 years) from the Nuclear Medicine Service at the University Hospital Vall d'Hebron (Barcelona). Tumors of these patients were classified as papillary (171), follicular (39), or Hürthle cell (6) carcinomas. This information was not available for 11 individuals at the moment of the study, and these were considered as unclassified. All individuals gave informed consent, and clearance from the ethical committee of our institutions was obtained.
DNA isolation was done using a standard phenol-chloroform method and 30 to 100 µL of Tris-EDTA [10 mmol/L Tris; 0.2 mmol/L EDTA (pH 7.5)].
Genotyping
Case and control populations were genotyped by analyzing six single-nucleotide polymorphisms markers localized in the 1p12-13 region: rs4659200, rs3765945, rs2145418, rs7515409, rs4658973, and rs1241. The selection of the single-nucleotide polymorphisms was based on the information available in the public databases, on allele frequency (common variants; minor allele frequency, >0.2), and on the possibility of PCR-RFLP analysis. The six single-nucleotide polymorphisms were chosen to cover a region of 2.4 Mb of the chromosome 1p12-13, and none of the selected single-nucleotide polymorphisms was considered haplotype tagged single-nucleotide polymorphisms by the HapMap database. Relative map position of these markers is shown in Fig. 1
.

View larger version (6K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 1. Map position of the six studied single-nucleotide polymorphisms at chromosome 1p12-13 and relative to the BAT-40 marker. Single-nucleotide polymorphisms associated with thyroid cancer susceptibility are indicated in bold.
|
|
Genotypes were generated by RFLP after PCR amplification. The rs7515409 and rs4658973 markers were genotyped on 202 control and 227 patient samples. From all originally collected samples, a reduced number was available in the subsequent analysis of the rs4659200, rs3765945, rs2145418, and rs1241 single-nucleotide polymorphisms; therefore, 136 control and 201 patient samples were only used to genotype these markers. Taking into consideration the variant allele frequencies (range, 0.24-0.48) of these polymorphisms, the sample size of the study was sufficient to detect odds ratios of 2.2 and 1.9, respectively, with
= 0.05 and a power of 95% (G-Power software; ref. 18).
Primers used to amplify each single-nucleotide polymorphism, as well as the PCR product lengths, the restriction enzymes used, and the RFLP fragments are shown in Table 1
.
View this table:
[in this window]
[in a new window]
|
Table 1. PCR amplification primers and restriction enzymes used for genotyping six single-nucleotide polymorphisms of the 1p12 region
|
|
PCRs were carried out in a final volume of 25 µL using 100 ng of DNA in 1x PCR buffer (10 mmol/L of Tris-HCl, 50 mmol/L of KCl, and 0.1% Triton X-100), 2.5 mmol/L of MgCl2, 0.2 mmol/L of each deoxynucleotide triphosphate, 0.1 µmol/L of each primer, and 0.75 U of Taq DNA polymerase (Promega). Amplifications were as follows: an initial denaturing step at 94°C for 4 min, 30 cycles at 94°C for 30 s, 54°C for 40 s, 72°C for 1 min, and a final step at 72°C for 4 min. PCR products were then digested using the restriction enzymes AflIII (New England Biolabs), SapI (New England Biolabs), AcsI (Roche), BsrI (New England Biolabs), TaqI (Roche), and NlaIII (New England Biolabs) for rs4659200, rs3765945, rs2145418, rs7515409, rs4658973, and rs1241, respectively. The digestion products were resolved on 2.5% agarose gels. Genotyping errors were discarded by successful duplicate analysis of 10 samples for each single-nucleotide polymorphism. In addition, an alternative genotyping method was also done in 117 control and 114 patient samples using the MassArray technique (Sequenom, Inc.) to genotype the rs2145418 and rs4658973 markers. RFLP and MassArray gave an identical genotype profile, in which sample errors were not detected.
Statistical Analysis
Allele frequencies for the six single-nucleotide polymorphisms were calculated from the genotype frequencies and compared by the
2 test or the Armitage's trend test, when Hardy-Weinberg equilibrium failed (that is, rs2145418 and rs1241). Hardy-Weinberg equilibrium was tested in the control population. The degree of pairwise linkage disequilibrium between markers was expressed as D' and r2 values. Although both coefficients represent the most common measures of linkage disequilibrium, the use of r2 is recommended in association studies (19).
Association between single-nucleotide polymorphisms and thyroid cancer susceptibility was calculated by the odds ratio and 95% confidence intervals (95% CI) using logistic regression analyses adjusted for age, sex, and alcohol and smoking habits.
We assumed a model of multiplicative allelic effect (codominant model). All analyses were done with a 5% level of significance using the SNP.stats software.4
 |
Results
|
|---|
We previously reported that the chromosome 1p12-13 was associated with thyroid cancer susceptibility in a case-control study using the microsatellite BAT-40 marker (13). In the present study, we have conducted a fine mapping of the susceptibility loci in this region by extending the analysis to an average distance of 2.4 Mb. We have used six additional polymorphic markers to examine for any association with an increased risk for thyroid cancer. As shown in Fig. 1, rs4659200 and rs3765945 are surrounding our marker of reference, BAT-40 (proximal extreme of the 2.4-Mb region); rs2145418, rs7515409, and rs4658973 map to the middle of the analyzed region; and rs1241 is located on the distal extreme of this region.
Because of the different frequency of females between our control (58%) and patient (74%) populations, for each single-nucleotide polymorphism, the genotype distribution about gender was analyzed before doing the association analysis, but no significant difference was found (data not shown). Consequently, the different proportion of females between control and patient groups does not influence our association analysis.
Single-nucleotide polymorphisms, nucleotide changes, and allele frequencies in cases and controls are shown in Table 2
. In the control population, the observed genotype frequencies did not differ from those expected from the Hardy-Weinberg equilibrium for the rs4659200, rs3765945, rs7515409, and rs4658973 polymorphisms. However, for rs1241 and rs2145418, significant deviations from Hardy-Weinberg equilibrium were observed (P = 0.007 and P = 0.005, respectively). In these cases, genotyping errors were discarded as being the cause of the Hardy-Weinberg deviations because nonerrors were detected by duplicate genotype analysis (see Materials and Methods). Appropriate correction for these deviations was then used when comparing allele frequencies between cases and controls for these two polymorphisms (Armitage's trend test), as shown in Table 2. On the other hand, we checked the pairwise linkage disequilibrium coefficients D' and r2, and we did not find strong linkage disequilibrium between any combination of markers (Table 3
), indicating an independent transmission pattern of inheritance. Even those markers with near position showed weak association, for example, rs4659200 and rs3765945 that are 48.5 kb apart and present a weak association (D' = 0.35, r2 = 0.04).
View this table:
[in this window]
[in a new window]
|
Table 3. Pairwise linkage disequilibrium coefficients D' and r2 between the six studied single-nucleotide polymorphisms
|
|
As shown in Table 2, significant differences between control and patient groups were found with respect to the allele frequencies of the rs2145418 and rs4658973 single-nucleotide polymorphisms. The G variant of rs2145418 is overrepresented in patients (P < 0.0001), whereas the G variant of rs4658973 is underrepresented in this cohort (P < 0.0001). Genotype distribution of the six single-nucleotide polymorphisms in the control and patient cohorts, and their odds ratios are shown in Table 4
. We want to point out that no difference was found between the crude odds ratio and the adjusted odds ratio for age, sex, alcohol, and tobacco, which discard these features acting as confounding factors in this study. Therefore, the crude values of the odds ratio were used to assess the association between thyroid cancer and polymorphic variants. Significant odds ratios for the rs2145418 and rs4658973 polymorphisms were found (P < 0.0001 for each marker; see Table 4), which indicates a strong association between thyroid cancer and these single-nucleotide polymorphisms. No association was detected for the other four polymorphisms studied.
View this table:
[in this window]
[in a new window]
|
Table 4. Genotype distribution and association analysis of six single-nucleotide polymorphisms at the 1p12-13 region
|
|
Association between thyroid cancer and the G-variant allele of the rs2145418 was shown by the significant odds ratios found, with values of 5.0 for heterozygous carriers (95% CI, 2.85-8.83; P < 0.0001) and of 9.2 for the homozygous (95% CI, 4.50-21.6; P < 0.0001). Therefore, carrying only one G allele is sufficient to increase 5-fold the risk for thyroid cancer. For the rs4658973 marker, the T allele is associated with thyroid cancer, which provides a second marker for thyroid cancer susceptibility. In this case, the T variant is the frequent allele in the control population (T, 0.55; G, 0.45), being even more frequent in patients (T, 0.73; G, 0.26). We used the homozygous carriers of the frequent allele as the reference genotype to calculate the odds ratio (see Table 4); therefore, in the case of the rs4658973 marker, the odds ratios for the GT and GG genotypes take the inverse value that is inferior than 1, that is, 0.40 (95% CI, 0.26-0.62; P < 0.0001) and 0.07 (95% CI, 0.03-0.18; P < 0.0001), respectively. The odds ratio of the T allele heterozygous carriers is 5.7 times higher than the odds ratio of the homozygous GG, indicating that the individuals with at least one risk allele (T variant) for this locus have increased risk for thyroid cancer. The rs2145418 and rs4658973 markers show weak association in the control population (D' = 0.05, r2 = 0.022), suggesting that the effect in thyroid cancer susceptibility is not due to linkage disequilibrium between both markers. Thus, each of these single-nucleotide polymorphisms would be associated to different unknown susceptibility genes that are randomly transmitted. On the other hand, because the interaction between genotypes at different loci might influence phenotype outcomes, we analyzed the interaction of genotypes at the rs2145418 and rs4658973 loci. However, nonstatistically significant interaction was found (P = 0.089). In addition, because clinicopathologic characteristics might influence the risk for thyroid cancer of the different genotypes of the rs2145418 and rs4658973 polymorphisms in the patients group, we investigated this effect, but no differences with respect to sex, age of cancer diagnosis, and type of thyroid cancer were found.
 |
Discussion
|
|---|
Studies on genetic susceptibility for thyroid cancer could have important relevance to understand the etiology of this pathology. Thus, association studies aimed to identify genes for thyroid cancer risk have been recently published (9-12), although no specific genes for thyroid cancer susceptibility have yet been described. Our previous study, using the BAT-40 marker in a case-control study, indicated that the chromosome 1p12-13 is related to thyroid cancer susceptibility (13). Therefore, in the present study, we aimed to map the susceptibility loci on this region. In this context, we have scanned a region of
2.4 Mb surrounding the BAT-40 marker by genotyping six single-nucleotide polymorphisms in a Spanish population. This study has allowed us to identify in our population two independent genetic susceptibility loci on chromosome 1p12.
Association analysis using a case-control design showed a strong association between the rs2145418 marker and thyroid cancer incidence, and similar results were found with the rs4658973 marker. Thus, the variant allele for the rs2145418 marker and the overrepresented allele for the rs4658973 marker in the patients group show an estimation of 5.0- and 5.7-fold increase of risk to thyroid cancer in heterozygous carries, respectively. This risk is double in homozygous carriers, which is in accordance with a codominant model. These polymorphisms are 377 kb distant, they show weak association, and no genotype interaction was found between them, indicating that their relation with thyroid cancer susceptibility is independent. At present, no risk factors for thyroid cancer or other type of cancer have been reported to map in this region, although alterations in the1p12-13 region have been described in different tumors (14-17).
The single-nucleotide polymorphisms analyzed in this study were selected as genetic markers for susceptibility loci. These single-nucleotide polymorphisms reside in noncoding sequences; therefore, we would not expect rs4658973 and rs2145418 to be the cause of their association with thyroid cancer. Nevertheless, the disease-associated rs4658973 marker maps into a region containing three overlapping genes (Fig. 1). The WDR3 (member of the conserved family of WD-repeat genes), the SPAG1 (sperm-associated antigen 17 gene), and the GDAP2 (ganglioside-induced differentiation-associated protein 2 gene). Specifically, rs4658973 is located within the 24 introns of the WDR3 gene, which overlaps the UTR 3'of the SPAG17 gene. In addition, the WDR3 and GDAP2 genes run in opposite directions, and both promoter regions are overlapped. Thereby, any of these genes could be related to thyroid cancer susceptibility. In particular, WDR3, as a member of the WDR-repeat proteins, could be involved in several cellular processes, including cell cycle progression and signal transduction (20, 21). An irregular expression of some of the WDR-repeat proteins has been found in certain types of cancer (22-25). On the other hand, GDAP2 belongs to the GDAP1-10 protein family involved in signal transduction pathways (26). Then, we cannot rule out the possible effect of the WDR3 and/or GDAP2 genes in the risk for thyroid cancer, although more investigations are needed to define the thyroid cancer susceptibility causal gene in this region.
On the other hand, there are neither defined nor predicted genes in the rs2145418 region. The closer gene region to rs2145418 is the region containing the rs4658973 marker that we have previously mentioned, which is 377 kb upstream (see Fig. 1). The other end of rs2145418, 550 kb downstream, maps the TBX15 gene involved in developmental processes (27, 28). Hence, nonevidences exist for hypothetical functional genes that could explain the strong association of the rs2145418 marker with thyroid cancer found in our study. Alternatively, the rs2145418 single-nucleotide polymorphism could reside in a regulatory sequence that might alter the expression of causal genes mapping near or a certain distance from this marker.
In conclusion and as we previously suggested, our results show a strong association between the chromosome 1p12 and thyroid cancer risk in a Spanish population. Our studies are the first describing susceptibility loci for thyroid cancer in this region and provide justification for further genetic analysis of the chromosome 1p12 region in other populations, as well as the investigation of the causal genes for thyroid cancer susceptibility in this region of chromosome 1.
 |
Acknowledgments
|
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank the subjects who participated in this study and the members of the Nuclear Medicine Service, Hospital Vall d'Hebron (Barcelona), for providing patient blood samples; Sara Gutiérrez and Alba Hernández for the support in collecting and isolating DNA samples; Lourdes Alarcon for the support in the PCR analysis; and Victor Moreno for the support in the statistical analysis.
 |
Footnotes
|
|---|
Grant support: Spanish Ministry of Education and Science (project SAF2007-63338), the Generalitat de Catalunya (CIRIT, 2005SGR-00136), and predoctoral fellowships from the Universitat Autònoma de Barcelona (A. Baida and M. Akdi) and the Caroline Foundation (E. González-Flores).
4 http://bioinfo.iconcologia.net/snpstats/ 
Received 3/15/07;
revised 2/21/08;
accepted 3/11/08.
 |
References
|
|---|
- Sherman SI. Thyroid carcinoma. Lancet 2003;361:501–11.[CrossRef][Medline]
- Hunt JL. Molecular mutations in thyroid carcinogenesis. Am J Clin Pathol 2002;118:S116–27.[Medline]
- Kondo T, Ezzat S, Asa Sl. Pathogenetic mechanisms in thyroid follicular-cell neoplasia. Nat Rev Cancer 2006;6:292–306.[CrossRef][Medline]
- Vecchio G, Santoro M. Oncogenes and thyroid cancer. Clin Chem Lab Med 2000;38:113–6.[CrossRef][Medline]
- Kazakov VS, Demidchik EP, Astakhova LN. Thyroid cancer after Chernobyl. Nature 1992;359:21.[Medline]
- Inskip PD. Thyroid cancer after radiotherapy for childhood cancer. Med Pediatr Oncol 2001;36:568–73.[CrossRef][Medline]
- Hernández A, Cespedes W, Xamena N, et al. Glutathione S-transferase polymorphisms in thyroid cancer patients. Cancer Lett 2003;190:37–44.[CrossRef][Medline]
- Cybulski C, Gorski B, Huzarski T, et al. CHEK2 is a multiorgan cancer susceptibility gene. Am J Hum Genet 2004;75:1131–5.[CrossRef][Medline]
- Iuliano R, Le Pera I, Cristofaro C, et al. The tyrosine phosphatase PTPRJ/DEP-1 genotype affects thyroid carcinogenesis. Oncogene 2004;23:8432–8.[CrossRef][Medline]
- Matakidou A, Hamel N, Polpat S, et al. Risk of non-medullary thyroid cancer influenced by polymorphic variation in the thyroglobulin gene. Carcinogenesis 2004;25:369–73.[Abstract/Free Full Text]
- Silva SN, Gil OM, Oliveira VC, et al. Association of polymorphisms in ERCC2 gene with non-familial thyroid cancer risk. Cancer Epidemiol Biomarkers Prev 2005;4:2407–12.
- Lonn S, Bhatti P, Alexander BH, et al. Papillary thyroid cancer and polymorphic variants in TSHR- and RET-related genes: a nested case-control study within a cohort of U.S. radiologic technologists. Cancer Epidemiol Biomarkers Prev 2007;6:174–7.
- Baida A, Farrington SM, Galofré P, et al. Thyroid cancer susceptibility and THRA1 and BAT-40 repeats polymorphisms. Cancer Epidemiol Biomarkers Prev 2005;14:638–42.[Abstract/Free Full Text]
- Rudolph B, Harbott J, Lampert F. Fragile sites and neuroblastoma: fragile site at 1p13.1 and other points on lymphocyte chromosomes from patients and family members. Cancer Genet Cytogenet 1988;31:83–94.[Medline]
- Bieche I, Champeme MH, Matifas F, et al. Demonstration of two regions involved in chromosome 1p deletion in breast tumors. Bull Cancer 1994;81:108–13.[Medline]
- Mitelman F, Mertens F, Johansson B. A breakpoint map of recurrent chromosomal rearrangements in human neoplasia. Nat Genet 1997;15:417–74.[CrossRef][Medline]
- Zhang J, Glatfelter AA, Taetle R, et al. Frequent alterations of evolutionarily conserved regions of chromosome 1 in human malignant melanoma. Cancer Genet Cytogenet 1999;111:119–23.[Medline]
- Faul F, Erdfelder E, Lang A-G, Buchner A. G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods 2007;39:175–91.[Medline]
- Ardlie KG, Kruglyak L, Seielstad M. Patterns of linkage disequilibrium in the human genome. Nat Rev Genet 2002;3:299–309.[CrossRef][Medline]
- Van der Voorn L, Ploegh HL. The WD-40 repeat. FEBS Lett 1992;307:131–134.[CrossRef][Medline]
- Smith TF, Gaitatzes C, Saxena K, et al. The WD repeat: a common architecture for diverse functions. Trends Biochem Sci 1999;245:181–5.
- Chernova OB, Hunyadi A, Malaj E, et al. A novel member of the WD-repeat gene family, WDR11, maps to the 10q26 region and is disrupted by a chromosome translocation in human glioblastoma cells. Oncogene 2001;20:5378–92.[CrossRef][Medline]
- Honore B, Baandrup U, Nielsen S, et al. Endonuclein is a cell cycle regulated WD-repeat protein that is up-regulated in adenocarcinoma of the pancreas. Oncogene 2002;21:1123–9.[Medline]
- Proikas-Cezanne T, Waddell S, Gaugel A, et al. WIPI-1
(WIPI49), a member of the novel 7-bladed WIPI protein family, is aberrantly expressed in human cancer and is linked to starvation-induced autophagy. Oncogene 2004;23:9314–25.[CrossRef][Medline] - Silva FP, Hamamoto R, Nakamura Y, et al. WDRPUH, a novel WD-repeat-containing protein, is highly expressed in human hepatocellular carcinoma and involved in cell proliferation. Neoplasia 2005;7:348–55.[CrossRef][Medline]
- Liu H, Nakagawa T, Kanematsu T, et al. Isolation of 10 differentially expressed cDNAs in differentiated Neuro2a cells induced through controlled expression of the GD3 synthase gene. J Neurochem 1999;72:1781–90.[CrossRef][Medline]
- Begemann G, Gibert Y, Meyer A, et al. Cloning of zebrafish T-box genes tbx15 and tbx18 and their expression during embryonic development. Mech Dev 2002;114:137–41.[CrossRef][Medline]
- Singh MK, Petry M, Haenig B, et al. The T-box transcription factor Tbx15 is required for skeletal development. Mech Dev 2005;122:131–44.[CrossRef][Medline]