Abstract

Objective. In this study, we aimed to investigate the associations of three single-nucleotide polymorphisms (SNPs) on TITF1/TITF2 (rs944289, rs965513, and rs1443434) with susceptibility to papillary thyroid carcinoma (PTC) and with nodular goiter (NG) in northern Chinese Han populations. Methods. We performed a case-control study comprising 861 PTC patients, 562 NG patients, and 896 normal controls (NCs). One TITF1 SNP (rs944289) and two TITF2 SNPs (rs965513 and rs1443434) were genotyped. Departures from Hardy–Weinberg equilibrium (HWE) in the control group were evaluated using chi-square test. Associations of the SNPs with PTC and with NG were assessed by unconditional logistic regression using the online SNPStats program. Bonferroni correction was performed for multiple tests in genotype analyses. Data analysis was performed by SPSS24.0 unless otherwise specified. Results. For rs944289, T allele was associated with increased risks for both PTC (OR = 1.23, 95% CI: 1.08–1.41, ) and NG (OR = 1.28, 95% CI: 1.10–1.50, ), and TT genotype significantly increased NG risk (recessive model, OR = 1.60, 95% CI: 1.22–2.10, ). For rs965513, no association was observed after Bonferroni correction. For rs1443434, G allele was associated with increased PTC risk (OR = 1.33, 95% CI: 1.10–1.61, ). Moreover, PTC risk increased with the number of total risk alleles of the three SNPs (OR = 1.25, 95% CI: 1.13–1.37, ). After stratified by gender, the risk effect of rs944289 T allele on PTC was only observed in females (OR = 1.29, 95% CI: 1.10–1.50, ). Individuals carrying rs944289-rs965513-rs1443434 haplotypes T-G-G and T-G-T had increased risks of PTC (OR = 1.82, 95% CI: 1.25–2.64, ) and NG (OR = 1.28, 95% CI: 1.06–1.54, ), respectively. Conclusions. There are associations of rs944289 and rs1443434 polymorphisms with PTC risk and association of rs944289 polymorphism with NG risk. Haplotypes T-G-G and T-G-T are risk haplotypes of PTC and NG, respectively.

1. Introduction

Thyroid cancer, as one of the most common tumor of endocrine system, is the fifth leading cancer for the estimated new cancer cases in females [1, 2]. Among all the histological subtypes of thyroid cancer, papillary thyroid carcinoma (PTC) is the most common type, accounting for about 75–85% of all cases with thyroid cancer.

Genetic factors are essential for the pathogenesis of thyroid cancer [3]. Thyroid cancer has the highest heritability (53%) among 15 common cancers, with very low shared environmental effects (1%) [4]. A joint study from five Nordic countries revealed that lifetime cumulative risk of female relatives of PTC patients (2%) is higher than that of males (1%); the more PTC patients a family has, the more risk of thyroid cancer their relatives has; and cumulative risk of thyroid cancer increases threefold in comparison with general population lifetime [5]. Moreover, a study from the Korea National Health and Nutrition Examination Surveys showed that factors related to thyroid function, such as thyroid-stimulating hormone (TSH) and free thyroxine (fT4), also present high heritability: the narrow-sense heritability for TSH is 54%, and that for free thyroxine (fT4) is 56%; the narrow-sense heritability is higher in females than that in males for TSH; and both TSH and fT4 show negative genetic correlations in females after adjusted for environmental effects [6]. In a twin study including 69 monozygotic and 45 dizygotic adult twin pairs, the environment-adjusted heritability of the thickness of thyroid isthmus and thyroid volume (left and right lobes) on thyroid cancer was 50% and 68–79%, respectively [7]. These evidences reflect the high heritability of thyroid cancer.

Molecular mechanisms (gene mutation, metastasis suppressors abnormally expression, epigenetic modification, microRNA regulation, and signaling pathway alteration) are involved in thyroid cancer pathogenesis [8, 9]. The BRAFV600E mutation and RAS mutation has been identified as the most prevalent oncogenic event in thyroid cancer [9]. Novel risk loci on different genes have further been identified and validated in different populations, providing a new dimension to this extensive line of thyroid cancer research [10, 11]. TITF1 and TITF2 are located on chromosome 14q13.3 and 9q22.33, respectively, encoding thyroid-specific transcription factors (TITF1 and TITF2) [12, 13]. TITF1 is a homeodomain-containing nuclear transcription factor, being crucial in thyroid development and differentiation by regulating thyroid-specific genes thyroglobulin (Tg), thyroperoxidase (TPO), and thyrotropin receptor (TSHR) [14]. TITF2 contains a forkhead domain and encodes a phosphorylated protein. Mutations in TITF2 lead to the Bamforth syndrome, penetrating as thyroid agenesis [15]. These indicate that important loci mutations on TITF1 and TITF2 may modify the expression of these two transcription factors and influence their function on thyroid gland. Recently, different single-nucleotide polymorphisms (SNPs) in TITF1 (also called NKX2-1 and TTF1)/TITF2 (also called FOXE1 and TTF2) are found in PTC, suggesting that rs944289, rs965513, and rs1443434 may be implicated in PTC risk [2, 13]. However, these studies are performed in populations with different ethnic background, conferring inconsistent results [1618].

In this study, we evaluated the association of rs944289, rs965513, and rs1443434 in TITF1/TITF2 with PTC risks and with nodular goiter (NG) in northern Chinese Han populations.

2. Materials and Methods

2.1. Patients and Samples

A total of 861 patients with PTC, 562 patients with NG, and 896 normal controls (NCs) were included in this study. The patients with PTC and the patients with NG were recruited from January 2012 to December 2014 in the China-Japan Union Hospital of Jilin University, China. The inclusion criteria of PTC and NG patients included the followings: (1) northern Chinese Han populations; (2) pathologically diagnosed based on the Revised American Thyroid Association Management Guidelines for patients with thyroid nodules and differentiated thyroid cancer [19]; (3) histologic type of PTC was pure papillary carcinoma other than follicular variant of PTC (fvPTC). The NCs were selected from a large-scale community-based cross-sectional survey of chronic disease and risk factors among adults in Jilin province conducted in July 2012. All the NCs were northern Chinese Han populations in community with the following exclusion criteria: (1) with metabolic disease (hypertension, diabetes, and dyslipidemia), (2) with any thyroid disease (thyroiditis, benign and malignant thyroid tumors), (3) with other malignancy, and (4) with mental disorders. Demographic and clinical information of the patients and controls were obtained through inquiry or consulting medical records. Peripheral blood was provided by each participant for genotyping. The study was approved by the ethics committee of the School of Public Health, Jilin University, and written informed consents were obtained from the patients and controls in the study.

2.2. Selection of SNPs

We selected rs944289 in TITF1, and rs965513 and rs1443434 in TITF2 in this study based on previous genome-wide association studies (GWAS) and confirmatory studies [18, 2022].

2.3. DNA Extraction and SNP Genotyping

Genomic DNA extraction and SNP genotyping were performed as in our previously published study [23]. Briefly, 5 ml of peripheral blood obtained from each patient and control was stored in EDTA-containing tubes. Genomic DNA was extracted using a commercial DNA extraction kit (ClotBlood DNA kit, Cwbio, Beijing, China) according to the manufacturer’s instructions. Primers for polymerase chain reaction (PCR) were designed by Assay Designer 3.1 (Supplementary Table S1). Genotypes were detected by matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF-MS) using the MassARRAY system (Sequenom, San, Diego, CA, USA). The genotyping rate of rs944289, rs965513, and rs1443434 was 99.09%, 97.15%, and 95.34, respectively.

2.4. Statistical Analysis

Continuous variables were presented with mean and standard deviation, and categorical variables were presented with percentage. Differences between groups were tested on the basis of independent-samples t-test or chi-square test using SPSS Statistics 24.0 (IBM Canada Ltd.). Departure from Hardy–Weinberg equilibrium (HWE) in NCs was calculated on the basis of chi-square test using SPSS Statistics 24.0 (IBM Canada Ltd.). Associations between SNPs and PTC/NG under five models of inheritance (codominant, dominant, recessive, overdominant, and allele model) were performed using online SNPStats program (http://bioinfo.iconcologia.net/SNPStats) [24]. Strength of association was assessed by odds ratio (OR) and corresponding 95% confidence intervals (CIs). The best model of inheritance for each SNP was selected on the basis of the Akaike information criterion (AIC). Bonferroni correction was performed to reduce type I error in multiple testing in view of 15 comparisons between groups (three SNPs × five models); thus, significant value was set at 0.05/15 = 0.003 for genotype analyses. Linkage disequilibrium analysis between the three SNPs and haplotype analysis was analyzed using the online SNPStats program. Power calculations were done using Quanto 1.2.4 software. Two-sided test with value less than 0.05 was considered statistically significant otherwise specified.

3. Results

3.1. Characteristics of the Study Subjects

The age and gender distribution were shown in Table 1. Significant difference of age distribution was observed between NG and NC (). The female-to-male ratio was more than 3.0 in the three groups, and no significant difference of gender distribution was found among the three groups (χ2 = 0.008, ). Clinical characteristics of PTC were shown in Supplementary Table S2. In this study, 15.7% of the patients with PTC had a family history; 34.7% of the patients with PTC were at stage 0-I; and 18.4% of the PTCs had lymph node metastasis.

3.2. Genotype and Allele Distributions of the Three SNPs

Significant differences of genotype and allele distributions among the groups were found for all the three SNPs (all ) (Table 2). For genotype distribution among the groups, chi-square test of linear-by-linear association was also performed, and significant trend associations were found for all the three SNPs (). No departure from HWE was observed in NC group (for rs944289, ; for rs965513, ; and for rs1443434, ). Linkage disequilibrium between the three SNPs was analyzed using the online SNPStats program, and no significant linkage disequilibrium was found () (Supplementary Table S3). To make our results more dependable and avoid false-positive results, we also analyzed the difference of genotype and allele distributions with 1000 genome Han Chinese in Beijing (CHB) as controls, and the results were shown in Supplementary Table S4. Except for rs1443434 (), no significant distribution difference was found for rs944289 and rs965513 ().

3.3. Association between the Three SNPs and PTC/NG

We used codominant, dominant, recessive, overdominant, and allele models to evaluate associations of rs944289, rs965513, and rs1443434 polymorphisms with PTC and with NG (Tables 35). Because of the significant age difference between NG and NC, evaluation of associations between SNPs and NG was adjusted by age. Codominant model was detected as the best model for all the three SNPs in the analysis of associations between the SNPs and PTC risk, and recessive model was detected as the best model for all the three SNPs in the analysis of association between the SNPs and NG risk. For rs944289, T allele was associated with increased risks for both PTC (OR = 1.23, 95% CI: 1.08–1.41, ) and NG (OR = 1.28, 95% CI: 1.10–1.50, ). TT genotype significantly increased PTC risk (codominant model, OR = 1.52, 95% CI: 1.15–1.99, ) and NG risk (recessive model, OR = 1.60, 95% CI: 1.22–2.10, ); however, the risk effect of TT genotype on PTC was not significant after Bonferroni correction. For rs965513, A allele was associated with increased PTC risk (OR = 1.34, 95% CI: 1.08–1.67, ), and AA genotype significantly increased PTC risk under codominant model (OR = 3.81, 95% CI: 1.25–11.63, ). NG risk increased significantly for AA genotype under recessive model (OR = 4.53, 95% CI: 1.41–14.59, ). However, after Bonferroni correction, no association was still significant for rs965513. For rs1443434, G allele was associated with increased PTC risk (OR = 1.33, 95% CI: 1.10–1.61, ). Under the codominant model, GG genotype of rs1443434 significantly increased PTC risk (OR = 2.21, 95% CI: 1.15–4.23, ), but the association was not significant after Bonferroni correction. Powers were calculated under dominant and recessive models using Quanto 1.2.4 software on the basis of the sample size and minor allele frequencies in our study (Supplementary Table S5). We also performed the association analyses with 1000 genome CHB as controls (Supplementary Table S6S8). No genotype or allele on the three SNPs was found significantly associated with PTC (). For rs1443434, G allele was associated with increased risk of NG (OR = 1.87, 95% CI: 1.05–3.31, ), but the association was not significant after Bonferroni correction.

3.4. Cumulative Risk Effect of the Three SNPs on PTC

We performed chi-square test of linear-by-linear association and logistic regression to evaluate the combined effect of the three SNPs on PTC risk (Table 6). Total counts of risk alleles for each subject were calculated, with a range of 0–5. As a result, percentages of PTCs with more risk alleles (three, four, and five) were higher than that of NGs and NCs (), and with the increasing of each risk allele, PTC risk increased 0.25-fold (OR = 1.25, 95% CI: 1.13–1.37, ).

3.5. Haplotype Analysis

Haplotype analysis revealed that individuals carrying rs944289-rs965513-rs1443434 haplotypes T-G-G and T-G-T had increased risks of PTC (OR = 1.82, 95% CI: 1.25–2.64, ) and of NG (OR = 1.28, 95% CI: 1.06–1.54, ) (Supplementary Table S9-S10), respectively.

3.6. Association between the Three SNPs and PTC Risks after Stratified by Gender

Genotype and allele frequencies of the three SNPs in females and males are shown in Supplementary Table S11. After stratified by gender, associations between the three SNPs and PTC risk under codominant, dominant, recessive, overdominant, and allele models are shown in Tables 79. Significant associations with PTC risks were found for rs944289 and rs965513 in females, and for rs1443434 in males. Codominant model was detected as the best model for all the significant comparisons. For rs944289, T allele (OR = 1.29, 95% CI: 1.10–1.50, ) and TT homozygote (codominant model: OR = 1.67, 95% CI: 1.22–2.28, ) significantly increased PTC risk in females, but the risk effect of TT genotype was not significant after Bonferroni correction. For rs965513, A allele (OR = 1.30, 95% CI: 1.01–1.67, ) and AA homozygote (codominant model: OR = 4.37, 95% CI: 0.92–20.65, ) significantly increased PTC risk in females, but the associations were not significant after Bonferroni correction any more. No significant association was found in males for rs965513 and rs944289. For rs1443434, G allele (OR = 1.74, 95% CI: 1.17–2.59, ) and GG homozygote (codominant model: OR = 4.91, 95% CI: 1.03–23.51, ) significantly increased PTC risk in males; however, the associations were not significant after Bonferroni correction. No significant association was found in females for rs1443434.

4. Discussion

In our present study, we evaluated the effect of TITF1 rs944289 and TITF2 rs965513 and rs1443434 polymorphisms on both PTC and NG, and we found that polymorphisms of rs944289 and rs1443434 are associated with PTC and NG, and PTC risk increases with the number of risk alleles of the three SNPs. After stratified by gender, a significant association between rs944289 polymorphism and PTC risk was only observed in females. Moreover, haplotypes T-G-G and T-G-T (rs944289-rs965513-rs1443434) are risk haplotypes of PTC and NG, respectively.

The polymorphism rs944289 predisposes to PTC through PTC susceptibility candidate 3 (PTCSC3) gene. PTCSC3 is a large intergenic noncoding RNA gene [25]. The risk allele of SNP rs944289 suppresses PTCSC3 by destroying a transcription factor-binding site in the promoter of PTCSC3 [26]. Our results showed that rs944289 polymorphism was associated with PTC risk, but no similar result was obtained when we used 1000 genome CHBs as controls, and we speculated that sample size was the main reason to the difference. In 1000 genome online database, only 45 CHBs was genotyped for the three SNPs, whereas in our own case-control study, 896 Han Chinese individuals were included. Moreover, geographical difference between 1000 genome CHBs (Beijing) and our controls (mainly from the Northeast China) may influence on these results. After stratified by gender, a significant association between rs944289 polymorphism and PTC risk was only observed in females, validating that the sex difference is involved in PTC pathogenesis. However, no statistically significant association between rs944289 and the risk of female breast cancer was also observed in a Chinese population [27], suggesting that rs944289 polymorphism may be required for female PTC specially, rather than breast cancer.

Although several meta-analyses have been performed to explore the association between these SNPs and thyroid cancer susceptibility, they still get no consistent result. Study of Chen et al. revealed the rs965513 polymorphism is significantly associated with risk of differentiated thyroid carcinoma (DTC) in Caucasians, but not in East Asians. The rs944289 polymorphism is associated with DTC in both Caucasians and East Asians [28]. Kang et al.’s study showed that the polymorphisms of rs944289 and rs965513 are not associated with thyroid cancer risk, but the rs1443434 polymorphism was significantly associated with thyroid cancer risk in overall population [29]. In Zhuang et al.’s study, rs965513 A allele is significantly associated with increased risk of thyroid cancer in both Caucasians and Asians [30]. In Gao et al.’s study, significant associations are found between the polymorphisms of rs944289 and rs965513 and risk of PTC risk in Caucasians and Asians [31]. In Figlioli et al.’s study, using GWAS and meta-analysis approaches confirms rs965513 as a risk locus of DTC [32].

Our present study focused on associations of the three SNPs (rs944289, rs965513, and rs1443434) and PTC risk in northern Chinese Han populations, whose ethnic background is different from the previous studies [1618]. As a hospital-based case-control designed study, we must admit the potential selection bias which is common in such studies. In addition, sample size of our study (861 PTCs, 562 NGs, and 896 NCs) is relatively small, influencing the statistical power. Besides, environmental factors, such as smoking and drinking, are also remarkable risk factors for thyroid cancer, which were not analyzed in our study. Furthermore, we only selected three SNPs of interest for analysis, not taking other molecular subtyping into account. BRAFV600E mutation, for example, which is an important mutation in thyroid cancer, may influence the genetic structure of the studied population, affecting the risk evaluation. Finally, studies on mechanism between rs944289 and rs1443434 polymorphisms and PTC should be performed to demonstrate the effect. Nevertheless, our study provides possible biomarkers for PTC diagnose and provides clews for further mechanism study, especially in northern Chinese Han populations.

5. Conclusions

The results of this study showed that polymorphisms of rs944289 and rs1443434 are associated with both PTC and NG in Chinese Han populations, and PTC risk increases with the number of total risk alleles of the three SNPs. Moreover, haplotypes T-G-G and T-G-T (rs944289-rs965513-rs1443434) are risk haplotypes of PTC and NG, respectively.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xin Zhang and Yulu Gu contributed equally to this work.

Acknowledgments

The authors would like to thank Lizhe Ai, Guang Yang, Qingqing, Zhang, Hongqin Xu, Chong Wang, and Wen Wen for their assistance in data collection and experiment. This work was supported by the National Natural Science Foundation of China (grant numbers 81602907, 81903389, and 81702606), the funds from Science and Technology Commission of Jilin Province (grant number 17240BGH010135414), and the funds from Commission of Health and Family Planning of Jilin Province (grant number 2017Q037). The National Natural Science Foundation of China plays role in the study design and in the manuscript writing, editing, approval, and decision to publish, and the others play role in sample collection, analysis, and interpretation of data.

Supplementary Materials

Table S1: the PCR amplification primer sequence. Table S2: clinical characteristics of PTCs (n = 861). Table S3: linkage disequilibrium between the three SNPs. Table S4: genotype and allele distributions of the three SNPs in PTCs, NGs, and 1000genome CHBs. Table S5: summary of Quanto results for power calculation. Table S6: association between rs944289 polymorphism and PTC/NG with 1000genome CHB as controls. Table S7: association between rs965513 polymorphism and PTC/NG with 1000genome CHB as controls. Table S8: association between rs1443434 polymorphism and PTC/NG with 1000genome CHB as controls. Table S9: associations between TITF1/TITF2 haplotypes and risk of PTC. Table S10: associations between TITF1/TITF2 haplotypes and risk of NG. Table S11: genotype and allele frequencies of the three SNPs by gender (n, %). (Supplementary Materials)