Abstract
Fang Wang1,*, Dehui Yan1,*, Xu Ji2,*, Jun Han3, Meijun Chen3, Hong Qiao3, Shaojun Zhang1
1College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, 150081, China
2Department of Otolaryngology, The First Affiliated Hospital of China Medical University, Shenyang, 110001, China
3Department of Endemic Disease, the Second Affiliated Hospital, Harbin Medical University, Harbin, 150086, China
*These authors contributed equally to this work
Correspondence to:
Shaojun Zhang, email: [email protected]
Hong Qiao, email: [email protected]
Keywords: thyroid cancer, genome-wide association studies, FOXE1, rs965513, meta-analysis
Received: November 06, 2015 Accepted: April 16, 2016 Published: May 12, 2016
ABSTRACT
Papillary thyroid cancer (PTC) is the most common type of thyroid cancer. With the rapid development of genome-wide association studies (GWAS), many genome variants associated with susceptibility to PTC have been identified, including the single nucleotide polymorphism rs965513 (9q22.33) near FOXE1. To evaluate the association between rs965513 and PTC in different ethnicities and countries, we conducted a meta-analysis using relatively large-scale samples from 23 studies (N = 163,136; 20,736 cases and 142,400 controls) by searching the PubMed and Google Scholar databases. Significant heterogeneity caused by different populations among the selected studies was observed. The A allele of rs965513 polymorphism was shown to be highly associated with risk of thyroid cancer, with odds ratios of 1.58 (95% CI 1.32–1.90) in all populations, 1.65 (95% CI 1.31–2.07)) in Caucasian populations and 1.49 in Asian populations. Compared to the dominant and recessive models, we observed the highest odds ratio (OR = 2.80, 95% CI 2.12–3.69) in the homozygous model. These results revealed that the rs965513 polymorphism is a risk factor for thyroid cancer
INTRODUCTION
Thyroid cancer (TC) is the most common malignancy in the endocrine system [1] and the fifth leading malignancy in female patients [2], and papillary thyroid cancer (PTC) is the most common type of thyroid cancer. Recent studies have improved our understanding of the pathogenesis of PTC, including the identification of genetic alterations that activate a common effector pathway involving the RET-Ras-BRAF signaling cascade, as well as other unique chromosomal rearrangements [3]. In addition, a large scale genome-wide association study (GWAS) identified many single-nucleotide polymorphisms (SNPs) that are significantly associated with PTC, such as USF1 on chromosome 1, FOXE1 on chromosome 9, ATM on chromosome 11, NKX2-1 on chromosome 14 [4], XRCC1 on chromosome 19, XRCC3 on chromosome 14 [5], and ALMS1 on chromosome 2 [6].
The rs965513 polymorphism located near FOXE1 was first identified as significantly associated with thyroid cancer by Gudmundsson J in 2009. That study showed that two common variants are associated with thyroid cancer, specifically rs965513 on 9q22.33 (OR = 1.75; P = 1.7 × 10−27) and rs944289 on 14q13.3 (OR = 1.37; P = 2.0 × 10−9) [7]. Currently, rs965513 has shown no association with events involved in the progression of PTC, such as invasion and tumor stage [8]. Due to the heterogeneity of susceptibility to cancer, such that no association between rs965513 and PTC has been observed in the US [9], it is important to investigate whether rs965513 is associated with thyroid cancer risk in all ethnicities and populations. Studies following Gudmundsson J have investigated this association in Cuba, France, Iceland, The US, Spain, China, Poland, Russia, The UK, Germany, Belarus, Portugal and Japan and provide the opportunity to evaluate whether rs965513 can be used as a common marker in other populations [4, 7–18].
Meta-analysis has been described as combining and analyzing quantitative evidence from related studies to produce results based on a whole body of research [19]. Thus, accounting for the importance of the variance and the inconsistency of results, we evaluated the genetic heterogeneity of rs965513 polymorphism in multiple populations by searching the PubMed and Google Scholar database and performed a meta-analysis to achieve a higher statistical power.
RESULTS
Literature search
Forty-one articles were selected from the PubMed and Google Scholar databases. Based on the inclusion and exclusion criteria, 13 articles, including 23 independent studies, were included in our analysis. More detailed information about the decision to include or exclude the selected studies can be found in Figure 1. In total, 20,736 cases and 142,400 controls from Cuba, France, Iceland, The US, Spain, China, Poland, Russia, The UK, Germany, Belarus, Portugal and Japan were included in our meta-analysis. Table 1 shows the main characteristics of the included studies: the name of the first author, the year of publication, the population or ethnicity, and the numbers of cases and controls.
Figure 1: The flow chart for identifying relevant studies.
Table 1: The main characteristics of the included studies
Author | Year | Countryn | Ethnicity | Allele G (freq) | Allele A (freq) | ||
|---|---|---|---|---|---|---|---|
Case | Control | Case | Control | ||||
Pereda CM | 2015 | Cuba | Mixed | 257 (0.640) | 317 (0.750) | 149 (0.360) | 107 (0.250) |
Maillard S | 2015 | French | Polynesian | 235 (0.730) | 392 (0.790) | 85 (0.270) | 104 (0.210) |
Gudmundsson J2009 | 2009 | Iceland | Caucasian | 590 (0.510) | 48206 (0.648) | 568 (0.490) | 26186 (0.352) |
Gudmundsson J2009 | 2009 | US | Caucasian | 312 (0.529) | 516 (0.671) | 276 (0.471) | 252 (0.329) |
Gudmundsson J2009 | 2009 | Spain | Caucasian | 98 (0.556) | 1768 (0.658) | 80 (0.444) | 918 (0.342) |
Wang YL | 2013 | China | Asian | 1489 (0.880) | 1847 (0.920) | 201 (0.120) | 163 (0.080) |
Liyanarachchi S | 2013 | US | Caucasian | 707 (0.490) | 1269 (0.680) | 733 (0.510) | 613 (0.320) |
Liyanarachchi S | 2013 | Poland | Caucasian | 1880 (0.540) | 2471 (0.650) | 1606 (0.460) | 1351 (0.350) |
Wei WJ | 2015 | China | Asian | 246 (0.893) | 1230 (0.876) | 30 (0.107) | 170 (0.124) |
Takahashi M | 2010 | Russia | Caucasian | 189 (0.487) | 238 (0.666) | 199 (0.513) | 120 (0.334) |
Takahashi M | 2010 | Russia | Caucasian | 224 (0.524) | 581 (0.648) | 204 (0.476) | 315 (0.352) |
Takahashi M | 2010 | Russia | Caucasian | 279 (0.538) | 820 (0.633) | 239 (0.462) | 476 (0.367) |
Jones AM | 2012 | UK | Caucasian | 768 (0.511) | 8225 (0.672) | 734 (0.489) | 4015 (0.328) |
Penna-Martinez | 2014 | Germany | Caucasian | 263 (0.549) | 348 (0.644) | 223 (0.451) | 192 (0.356) |
Damiola | 2014 | Byelorussia | Caucasian | 78 (0.600) | 303 (0.584) | 56 (0.400) | 153 (0.416) |
Tomaz (FNMTC) | 2012 | Portugal | Caucasian | 53 (0.442) | 168 (0.646) | 67(0.558) | 92(0.354) |
Tomaz (NMTC) | 2012 | Portugal | Caucasian | 63 (0.394) | 168 (0.646) | 97 (0.606) | 92 (0.354) |
Tomaz (all) | 2012 | Portugal | Caucasian | 116 (0.414) | 168 (0.646) | 164 (0.586) | 92 (0.354) |
Denny | 2011 | US | Caucasian | 1881 (0.714) | 6589 (0.652) | 753 (0.286) | 3517 (0.348) |
Denny | 2011 | US | Caucasian | 1881 (0.714) | 7376 (0.655) | 753 (0.286) | 3888 (0.345) |
Matsuse | 2011 | Japan | Asian | 872 (0.910) | 5213 (0.943) | 86 (0.090) | 315 (0.057) |
Matsuse | 2011 | Japan | Asian | 692 (0.908) | 5213 (0.943) | 70 (0.092) | 315 (0.057) |
Matsuse | 2011 | Japan | Asian | 175 (0.921) | 5213 (0.943) | 15 (0.079) | 315 (0.057) |
Mixed in Ethnicity was the mixture of European and African.
Heterogeneity test
The genetic heterogeneity of the rs965513 polymorphism was evaluated based on the additive, dominant, recessive and homozygous models and the data from the selected studies (Table 2). Significant heterogeneity was observed among these studies. In the additive model (A vs. G) and the dominant model (AA + AG vs. GG), extreme heterogeneity was observed among the 23 selected studies (additive model: P < 0.0001 and I2 = 95.4%; dominant model: P < 0.0001 and I2 = 76.7%). The recessive model (AA vs. AG + GG) and the homozygous model (AA vs. GG) showed large heterogeneity among the 23 selected studies (recessive model: P = 0.0061 and I2 = 62.6%; homozygous model: P = 0.0016 and I2 = 67.8%). There was no significant heterogeneity observed in Asian populations. However, we found significant heterogeneity in Caucasian populations.
Table 2: The result of heterogeneity test
Population | Risk model | Q | P | I2 | 95% CI |
|---|---|---|---|---|---|
All | Additive (A vs G) | 480.41 | < 0.0001 | 95.4% | 94.1%–96.4% |
Recessive (AA vs AG+GG) | 21.41 | 0.0061 | 62.6% | 23.1%–81.9% | |
Dominant (AA +AG vs GG) | 34.36 | < 0.0001 | 76.7% | 55.5%–87.8% | |
Homozygous (AA vs GG) | 24.85 | 0.0016 | 67.8% | 35.2%–84% | |
Asian | Additive (A vs G) | 7.55 | 0.1091 | 47% | 0%–80.6% |
Recessive (AA vs AG+GG) | 1.54 | 0.214 | 35.2% | — | |
Dominant (AA +AG vs GG) | 4.41 | 0.0358 | 77.3% | — | |
Homozygous (AA vs GG) | 1.69 | 0.1932 | 40.9% | — | |
Caucasian | Additive (A vs G) | 467.45 | < 0.0001 | 96.8% | 95.8%–97.5% |
Recessive (AA vs AG+GG) | 18.19 | 0.0027 | 72.5% | 36.7%–88.1% | |
Dominant (AA +AG vs GG) | 19.68 | 0.0014 | 74.6% | 42.3–88.8% | |
Homozygous (AA vs GG) | 19.83 | 0.0013 | 74.8% | 42.8%–88.9% |
— represent the missing data because that there were few studies in Asian population 95% CI of I2 did not can be estimate.
We used meta-regression to further investigate the potential sources of heterogeneity, including publication year, country, ethnicity, sample size, bias of sample size, type of control population and SNP genotyping techniques. We found that publication year, sample size and bias of sample size were significant, and these factors may be the sources of the observed heterogeneity (Supplementary Table 1). The absence of significant effects of country, ethnicity and type of control population implied that these factors cannot be cause of the heterogeneity among the studies.
Meta-analysis
We performed a meta-analysis to calculate the overall ORs using the random effect model in all populations and in Caucasian populations and using the fixed effect model in Asian populations based on heterogeneity analysis. The risk of thyroid cancer associated with the A allele was 1.58-fold that of the G allele (Figure 2A, OR = 1.58, 95% CI 1.32–1.90). Moreover, we found that Caucasian populations had higher risk than Asian populations (Figure 2B and 2C, OR = 1.65 vs. 1.49), and this result could be strengthened by including more studies. Because 9 of the 23 selected studies provided the number of rs965513 genotypes or provided sufficient data to calculate the number of rs965513 genotypes, meta-analyses using dominant, recessive and homozygous models were conducted among these 9 studies. The association between rs965513 and thyroid cancer was also significant in the dominant model (Figure 3A, AA + AG vs. GG, OR = 1.78, 95% CI 1.48– 2.15), the recessive model (Figure 3B, AA vs. AG + GG, OR = 2.10, 95% CI 1.66–2.64), and the homozygous model (Figure 3C, AA vs. GG, OR = 2.80, 95% CI 2.12–3.69). In addition, the genotype AA was most strongly associated with risk of thyroid cancer in Caucasian populations (Figure 3D–3F). Only two studies provided the number of rs965513 genotypes, so the risks of genotypes in Asian populations could not be estimated.
Figure 2: The forest plot for A vs G of the rs965513 polymorphism. (A) All populations. (B) Caucasian populations. (C) Asian populations.
Figure 3: The forest plot for genotypes of the rs965513 polymorphism. (A) Dominant model in all populations. (B) Recessive model in all populations. (C) Homozygous model in all populations. (D) Dominant model in Caucasian populations. (E) Recessive model in Caucasian populations. (F) Homozygous model in Caucasian populations.
Sensitivity and publication bias analysis
We performed a one-way sensitivity analysis to evaluate the robustness of the results of this meta-analysis. The pooled ORs from different populations were not influenced by removal of one study under four genetic models (Supplementary Tables 2 and 3), suggesting that the results of this meta-analysis are stable. Begg’s and Egger’s tests were performed to evaluate publication bias. Although the funnel plots of the additive model were asymmetrical inverted funnels (Figure 4A), the results of both Begg’s test and Egger’s test were not significant (Table 3). In addition, the funnel plots of the dominant, recessive and homozygous models are symmetrical inverted funnels (Figure 4B–4D), which suggests no significant publication bias. The above results suggest that the conclusions of our meta-analysis are credible and stable.
Figure 4: The funnel plots for publication bias analysis of the rs965513 polymorphism. (A) Additive model. (B) Dominant model. (C) Recessive model. (D) Homozygous model.
Table 3: The result of Begg and Egger’s tests
Risk model | Egger’ s test | Begg’ s test | ||
|---|---|---|---|---|
T statistic | P value | Z statistic | P value | |
Additive | 2.07 | 0.05086 | –0.87 | 0.38 |
Recessive | –0.85 | 0.42 | –0.28 | 0.36 |
Dominant | –1.44 | 0.19 | –0.56 | 0.04 |
Homozygous | –1.31 | 0.23 | –0.22 | 0.48 |
DISCUSSION
FOXE1, which is also called TTF2 (Thyroid transcription factor 2), is a transcription factor involved in thyroid gland development (thyroid formation, migration and morphogenesis control [20–23]) and in the maintenance of differentiation in the thyroid [24] and which is highly expressed in thyroid follicular cells [25, 56]. According to Goldgar DE and Eng C, multiple low- to moderate-penetrance genes (LPGs) interacting with each other and with the environment may result in thyroid cancer [27, 28]. FOXE1 is a likely LPG in this content because FOXE1 is the center of a regulatory network of transcription factors and cofactors that initiate thyroid differentiation [22]. The influence of FOXE1 on thyroid cancer has been investigated by Rihab Kallel [24]. The llelic, genotypic and phenotypic analyses strongly suggested that the length of the alanine stretch in FOXE1 modulates genetic susceptibility to papillary thyroid cancer. They reported that the 16-Ala allele and homozygous 16/16 genotype showed increased risk of thyroid cancer development. However, subjects with the 14-Ala allele seemed to be protected against the occurrence of this pathology [24]. The presence of rs965513 near FOXE1, located on 9q22, showed a ~1.8-fold odds ratio of PTC risk through genome-wide association. He et al. investigated the molecular mechanism by which rs965513 regulated the expression of FOXE1, generating susceptibility to thyroid cancer [29]. They found rs965513 located in a linkage disequilibrium block ~33 kb including at least three regulatory elements functioning as enhancers. The region overlapped with the promoter region of FOXE1, and the variability of genotypes was associated with differential activity levels of an enhancer, further leading to variations in FOXE1 expression that resulted in altered risk of thyroid cancer.
The significant association between the rs965513 polymorphism and thyroid cancer was first identified by Gudmundsson J in 2009 [7]. Over the following few years, subsequent studies have continued to explore this association and have reported both consistent and conflicting results. Geng et al. reviewed 10 studies and observed that the A allele of rs965513 had a 1.31-fold risk of thyroid cancer [30]. However, Kang et al. evaluated the association of 12 SNPs in FOXE1 and PTC, and reported that rs965513 showed no association with PTC [31]. It is important to assess the genetic architecture of the rs965513 polymorphism across different ethnicities and populations. There were several similar meta-analyses that have assessed the risk of rs965513 in PTC [31–35]. However, we used a larger sample size of 13 articles corresponding to 23 studies (20,736 cases vs. 142,400 controls) found in PubMed and Google Scholar to reevaluate this association, which will help to accurately assess the risk of rs965513 in PTC. Kang J. et al. included 8 studies corresponding to 2,085 cases and 10,341 controls [31]. Zhuang Y. et al. included 13 studies involving 8491 cases, 103,218 control and 629 family members [32]. Ai L et al. included 6 studies and a total of 52,363 individuals (5,193 cases vs. 47,170 controls) [33]. Zhu et al. included 14 studies and a total of 9828 subjects [34]. Gao et al. included 16 studies and 8119 cases vs. 66,936 controls [35]. Moreover, we evaluated the risk of thyroid cancer with rs965513 under several genetic models including the additive, recessive, dominant and homozygous models, to assess the increased levels of PTC risk under different genotypes. Zhuang et al. evaluated the allelic, dominant and recessive models [32]. Ai L et al. performed meta-analysis using additive model [33]. Zhu et al. performed meta-analysis using the additive, heterozygous and homozygous models [34]. Lastly, we considered differential associations between rs965513 and PTC in different populations and performed stratified meta-analysis separately in Asian and Caucasian populations. Ai L et al. performed meta-analysis in a mixed population [33].
In our meta-analysis, the genetic heterogeneity of rs965513 among the selected studies was evaluated, and significant heterogeneity was observed in the additive, dominant, recessive and homozygous models. The genetic heterogeneity may be caused by differences in the publication year, sample size and bias of sample size through meta-regression. Because more studies corresponded to Caucasian populations, significant heterogeneity was observed in Caucasian populations but not in Asian populations. Next, meta-analyses were conducted separately for Caucasian, Asian and all populations. Our results showed that the A allele of rs965513 had a 1.58-fold risk of thyroid cancer in all populations, a 1.65-fold risk (95% CI 1.31–2.07)) in Caucasian populations and a 1.49-fold risk in Asian populations. Compared to the dominant and recessive models, the homozygous model showed the highest odds ratio (OR = 2.80, 95% CI 2.12–3.69) in all populations and in the Caucasian populations. The one-way sensitivity analysis suggested that the results of this meta-analysis were stable. We found that the Caucasian populations had higher risk than Asian populations, which further supports the previous findings.
In addition to the rs965513 polymorphism, two other variants of FOXE1, rs1867277 and rs71369530, are also significantly associated with thyroid cancer [4]. The variant rs1867277 is located within the 5′ untranslated region (UTR) and is involved in the allele-specific transcriptional regulation of FOXE1 through recruitment of the transcription factors USF1/USF2 [4]. Jones et al. reconstructed haplotypes at these two loci (rs965513 and rs1867277) and estimated the ORs associated with having each of the three possible risk haplotypes compared to the non-risk haplotype. They reported that carrying the haplotype with both risk alleles significantly increases the risk of thyroid cancer, while carrying haplotypes with a single risk allele at either rs965513 or rs1867277 somewhat increases the risk [14]. The variant rs71369530 is a poly-alanine expansion in the FOXE1 coding region. Martyn’s study revealed that the poly-alanine expansion of FOXE1 (rs71369530) is significantly associated with PTC in Caucasian subjects (OR = 2.23, 95% CI 1.42–3.50) [36].
We conducted meta-analyses using additive, dominant, and recessive models at the same time for more powerful results. However, our study has certain limitations. We attempted to obtain exact genotype numbers from all studies used in our analysis for the dominant and recessive models, but only some of the genotype numbers were available [4, 8, 10–16]. Future studies can supplement our results using dominant and recessive models.
MATERIALS AND METHODS
Literature search
In stage 1, we searched PubMed and Google Scholar to select all possible studies with key words including ‘rs965513’ and ‘thyroid cancer’ or ‘FOXE1’ and ‘thyroid cancer’. The literature search was updated on April 7, 2015. Then, in stage 2, we used Google Scholar (http://scholar.google.com/) to query the articles citing these studies and all of the references therein as identified using PubMed. We selected only published articles written in English.
Inclusion criteria
The studies that were selected met the following criteria: (1) the study must use a case-control design; (2) the study evaluated the association between rs965513 polymorphism and thyroid cancer; (3) the study provided the number of rs965513 genotypes or (4) provided sufficient data to calculate the number of rs965513 genotypes; and (5) the study provided an OR with 95% CI as well as the P value or (6) the study provided sufficient data to calculate the OR and 95% CI.
Data extraction
We extracted the following information from every study: (1) the name of the first author; (2) the year of publication; (3) the population or ethnicity; (4) the numbers of cases and controls; (5) the genotype number of rs965513 polymorphisms in cases and controls; (6) the number of rs965513 genotypes or (7) information needed to calculate the numbers of rs965513 genotypes; (8) the OR with 95% CI or (9) information needed to calculate the OR and 95% CI. All related calculations were completed using R software.
Genetic model
The rs965513 polymorphism includes the two alleles G and A, of which A is the minor allele. A is assumed to be the high-risk allele and G the low-risk allele. We selected the additive, dominant, recessive and homozygous genetic model for further meta-analysis. The additive model can be described as the A allele versus the G allele [37].
Heterogeneity test
Cochran’s Q test was used to evaluate the genetic heterogeneity among the included studies, which approximately follows a χ2 distribution with k-1 degrees of freedom (k stands for the number of studies for analysis). I2 = (Q-(k-1)) ÷ Q × 100%, ranging from 0 to 100% [38], was also used. I2 is a measure of heterogeneity and is a statistic that indicates the percentage of variance in a meta-analysis that is attributable to study heterogeneity [39]. Low, moderate, large and extreme heterogeneity corresponded to 0–25%, 25–50%, 50–75% and 75–100%, respectively [39]. The significance levels for heterogeneity are defined as P < 0.01 and I2 > 50%. Meta-regression was used to investigate the potential sources of heterogeneity and was performed using the metafor package in R software.
Meta-analysis
For case in which there is no significant heterogeneity among the included studies, the pooled OR was calculated using the fixed effect model; otherwise, the OR was calculated using random-effect model. The Z test was used to determine the significance of ORs. All statistical tests for heterogeneity and meta-analysis were performed using the meta package in R software.
Publication bias analyses
Funnel plots were used to evaluate the potential publication bias [40]. Begg’s and Egger’s tests were used to evaluate the asymmetry of the funnel plot [40]. It is assumed that in the absence of publication bias, the largest study will be plotted near the average, and smaller studies will be spread evenly on both sides of the average, creating a roughly funne-shaped distribution. Deviation from this shape can indicate publication bias. The publication bias analyses were performed using the metafor package in R software.
ACKNOWLEDGMENTS AND FUNDING
We are thankful for the support of the National Natural Science Foundation of China [grant numbers 31401075, 61402139] and the Natural Scientific Research Fund of Heilongjiang Province [grant number QC2013C018].
CONFLICTS OF INTEREST
The authors declare no competing financial interests.
Authors’ contributions
SZ and HQ conceived the hypothesis and revised the manuscript. DY and XJ searched the literature and extracted the data. FW performed the meta-analysis and wrote the manuscript. JH and MC filtered the literature.
REFERENCES
1. DeLellis RA. Pathology and genetics of thyroid carcinoma. J Surg Oncol. 2006; 94:662–9.
2. Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013; 63:11–30.
3. Kondo T, Ezzat S, Asa SL. Pathogenetic mechanisms in thyroid follicular-cell neoplasia. Nat Rev Cancer. 2006; 6:292–306.
4. Pereda CM, Lesueur F, Pertesi M, Robinot N, Lence-Anta JJ, Turcios S, Velasco M, Chappe M, Infante I, Bustillo M, Garcia A, Clero E, Xhaard C, et al. Common variants at the 9q22.33, 14q13.3 and ATM loci, and risk of differentiated thyroid cancer in the Cuban population. BMC Genet. 2015; 16:22.
5. Wang X, Zhang K, Liu X, Liu B, Wang Z. Association between XRCC1 and XRCC3 gene polymorphisms and risk of thyroid cancer. Int J Clin Exp Pathol. 2015; 8:3160–7.
6. Papadakis M, Meyer A, Schuster F, Weyerbrock N, Corinth C, Dotzenrath C. Follicular variant of papillary thyroid cancer in Alstrom syndrome. Fam Cancer. 2015; 14:599–602.
7. Gudmundsson J, Sulem P, Gudbjartsson DF, Jonasson JG, Sigurdsson A, Bergthorsson JT, He H, Blondal T, Geller F, Jakobsdottir M, Magnusdottir DN, Matthiasdottir S, Stacey SN, et al. Common variants on 9q22.33 and 14q13.3 predispose to thyroid cancer in European populations. Nat Genet. 2009; 41:460–4.
8. Wei WJ, Lu ZW, Wang Y, Zhu YX, Wang YL, Ji QH. Clinical significance of papillary thyroid cancer risk loci identified by genome-wide association studies. Cancer Genet. 2015; 208:68–75.
9. Denny JC, Crawford DC, Ritchie MD, Bielinski SJ, Basford MA, Bradford Y, Chai HS, Bastarache L, Zuvich R, Peissig P, Carrell D, Ramirez AH, Pathak J, et al. Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies. Am J Hum Genet. 2011; 89:529–42.
10. Maillard S, Damiola F, Clero E, Pertesi M, Robinot N, Rachedi F, Boissin JL, Sebbag J, Shan L, Bost-Bezeaud F, Petitdidier P, Doyon F, Xhaard C, et al. Common variants at 9q22.33, 14q13.3, and ATM loci, and risk of differentiated thyroid cancer in the French Polynesian population. PLoS One. 2015; 10:e0123700.
11. Wang YL, Feng SH, Guo SC, Wei WJ, Li DS, Wang Y, Wang X, Wang ZY, Ma YY, Jin L, Ji QH, Wang JC. Confirmation of papillary thyroid cancer susceptibility loci identified by genome-wide association studies of chromosomes 14q13, 9q22, 2q35 and 8p12 in a Chinese population. J Med Genet. 2013; 50:689–95.
12. Liyanarachchi S, Wojcicka A, Li W, Czetwertynska M, Stachlewska E, Nagy R, Hoag K, Wen B, Ploski R, Ringel MD, Kozlowicz-Gudzinska I, Gierlikowski W, Jazdzewski K, et al. Cumulative risk impact of five genetic variants associated with papillary thyroid carcinoma. Thyroid. 2013; 23:1532–40.
13. Takahashi M, Saenko VA, Rogounovitch TI, Kawaguchi T, Drozd VM, Takigawa-Imamura H, Akulevich NM, Ratanajaraya C, Mitsutake N, Takamura N, Danilova LI, Lushchik ML, Demidchik YE, et al. The FOXE1 locus is a major genetic determinant for radiation-related thyroid carcinoma in Chernobyl. Hum Mol Genet. 2010; 19:2516–23.
14. Jones AM, Howarth KM, Martin L, Gorman M, Mihai R, Moss L, Auton A, Lemon C, Mehanna H, Mohan H, Clarke SE, Wadsley J, Macias E, et al. Thyroid cancer susceptibility polymorphisms: confirmation of loci on chromosomes 9q22 and 14q13, validation of a recessive 8q24 locus and failure to replicate a locus on 5q24. J Med Genet. 2012; 49:158–63.
15. Penna-Martinez M, Epp F, Kahles H, Ramos-Lopez E, Hinsch N, Hansmann ML, Selkinski I, Grunwald F, Holzer K, Bechstein WO, Zeuzem S, Vorlander C, Badenhoop K. FOXE1 association with differentiated thyroid cancer and its progression. Thyroid. 2014; 24:845–51.
16. Damiola F, Byrnes G, Moissonnier M, Pertesi M, Deltour I, Fillon A, Le Calvez-Kelm F, Tenet V, McKay-Chopin S, McKay JD, Malakhova I, Masyakin V, Cardis E, et al. Contribution of ATM and FOXE1 (TTF2) to risk of papillary thyroid carcinoma in Belarusian children exposed to radiation. Int J Cancer. 2014; 134:1659–68.
17. Tomaz RA, Sousa I, Silva JG, Santos C, Teixeira MR, Leite V, Cavaco BM. FOXE1 polymorphisms are associated with familial and sporadic nonmedullary thyroid cancer susceptibility. Clin Endocrinol (Oxf). 2012; 77:926–33.
18. Matsuse M, Takahashi M, Mitsutake N, Nishihara E, Hirokawa M, Kawaguchi T, Rogounovitch T, Saenko V, Bychkov A, Suzuki K, Matsuo K, Tajima K, Miyauchi A, et al. The FOXE1 and NKX2–1 loci are associated with susceptibility to papillary thyroid carcinoma in the Japanese population. J Med Genet. 2011; 48:645–8.
19. Riley RD, Lambert PC, Abo-Zaid G. Meta-analysis of individual participant data: rationale, conduct, and reporting. BMJ. 2010; 340:c221.
20. Dathan N, Parlato R, Rosica A, De Felice M, Di Lauro R. Distribution of the titf2/foxe1 gene product is consistent with an important role in the development of foregut endoderm, palate, and hair. Dev Dyn. 2002; 224:450–6.
21. De Felice M, Ovitt C, Biffali E, Rodriguez-Mallon A, Arra C, Anastassiadis K, Macchia PE, Mattei MG, Mariano A, Scholer H, Macchia V, Di Lauro R. A mouse model for hereditary thyroid dysgenesis and cleft palate. Nat Genet. 1998; 19:395–8.
22. Parlato R, Rosica A, Rodriguez-Mallon A, Affuso A, Postiglione MP, Arra C, Mansouri A, Kimura S, Di Lauro R, De Felice M. An integrated regulatory network controlling survival and migration in thyroid organogenesis. Dev Biol. 2004; 276:464–75.
23. Zannini M, Avantaggiato V, Biffali E, Arnone MI, Sato K, Pischetola M, Taylor BA, Phillips SJ, Simeone A, Di Lauro R. TTF-2, a new forkhead protein, shows a temporal expression in the developing thyroid which is consistent with a role in controlling the onset of differentiation. EMBO J. 1997; 16:3185–97.
24. Kallel R, Belguith-Maalej S, Akdi A, Mnif M, Charfeddine I, Galofre P, Ghorbel A, Abid M, Marcos R, Ayadi H, Velazquez A, Hadj Kacem H. Genetic investigation of FOXE1 polyalanine tract in thyroid diseases: new insight on the role of FOXE1 in thyroid carcinoma. Cancer Biomark. 2010; 8:43–51.
25. Civitareale D, Lonigro R, Sinclair AJ, Di Lauro R. A thyroid-specific nuclear protein essential for tissue-specific expression of the thyroglobulin promoter. EMBO J. 1989; 8:2537–42.
26. Francis-Lang H, Zannini M, De Felice M, Berlingieri MT, Fusco A, Di Lauro R. Multiple mechanisms of interference between transformation and differentiation in thyroid cells. Mol Cell Biol. 1992; 12:5793–800.
27. Goldgar DE, Easton DF, Cannon-Albright LA, Skolnick MH. Systematic population-based assessment of cancer risk in first-degree relatives of cancer probands. J Natl Cancer Inst. 1994; 86:1600–8.
28. Eng C. Familial papillary thyroid cancer—many syndromes, too many genes? J Clin Endocrinol Metab. 2000; 85:1755–7.
29. He H, Li W, Liyanarachchi S, Srinivas M, Wang Y, Akagi K, Wang Y, Wu D, Wang Q, Jin V, Symer DE, Shen R, Phay J, et al. Multiple functional variants in long-range enhancer elements contribute to the risk of SNP rs965513 in thyroid cancer. Proc Natl Acad Sci U S A. 2015; 112:6128–33.
30. Geng P, Ou J, Li J, Liao Y, Wang N, Xie G, Sa R, Liu C, Xiang L, Liang H. TITF1 and TITF2 loci variants indicate significant associations with thyroid cancer. Endocrine. 2015; 50:598–607.
31. Kang J, Deng XZ, Fan YB, Wu B. Relationships of FOXE1 and ATM genetic polymorphisms with papillary thyroid carcinoma risk: a meta-analysis. Tumour Biol. 2014; 35:7085–96.
32. Zhuang Y, Wu W, Liu H, Shen W. Common genetic variants on FOXE1 contributes to thyroid cancer susceptibility: evidence based on 16 studies. Tumour Biol. 2014; 35:6159–66.
33. Ai L, Liu X, Yao Y, Yu Y, Sun H, Yu Q. Associations between rs965513/rs944289 and papillary thyroid carcinoma risk: a meta-analysis. Endocrine. 2014; 47:428–34.
34. Zhu H, Xi Q, Liu L, Wang J, Gu M. Quantitative assessment of common genetic variants on FOXE1 and differentiated thyroid cancer risk. PLoS One. 2014; 9:e87332.
35. Gao Y, Chen F, Niu S, Lin S, Li S. Replication and Meta-Analysis of Common Gene Mutations in TTF1 and TTF2 with Papillary Thyroid Cancer. Medicine (Baltimore). 2015; 94:e1246.
36. Bullock M, Duncan EL, O’Neill C, Tacon L, Sywak M, Sidhu S, Delbridge L, Learoyd D, Robinson BG, Ludgate M, Clifton-Bligh RJ. Association of FOXE1 polyalanine repeat region with papillary thyroid cancer. J Clin Endocrinol Metab. 2012; 97:E1814–9.
37. Lewis CM, Knight J. Introduction to genetic association studies. Cold Spring Harb Protoc. 2012; 2012:297–306.
38. Liu G, Zhang S, Cai Z, Ma G, Zhang L, Jiang Y, Feng R, Liao M, Chen Z, Zhao B, Li K. PICALM gene rs3851179 polymorphism contributes to Alzheimer’s disease in an Asian population. Neuromolecular Med. 2013; 15:384–8.
39. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003; 327:557–60.
40. Jiang Y, Zhang R, Zheng J, Liu P, Tang G, Lv H, Zhang L, Shang Z, Zhan Y, Lv W, Shi M. Meta-analysis of 125 rheumatoid arthritis-related single nucleotide polymorphisms studied in the past two decades. PLoS One. 2012; 7:e51571.
