Associations between RET tagSNPs and their haplotypes and susceptibility, clinical severity, and thyroid function in patients with differentiated thyroid cancer

Abstract

Background

It is unclear whether common genetic variants of the RET proto-oncogene contribute to disease susceptibility, clinical severity, and thyroid function in differentiated thyroid cancer (DTC).

Methods

A total of 300 DTC patients and 252 healthy controls were enrolled in this study. Seven RET tagging single nucleotide polymorphisms were genotyped using the KASPar platform.

Results

Subgroup analysis showed that concomitant thyroid benign diseases were less likely to occur in DTC subjects with the rs1799939 AG or AG plus AA genotypes (odds ratio (OR) = 1.93 and 1.88, P = 0.009 and 0.011, respectively). A rare haplotype, CGGATAA, was associated statistically with a reduced risk of DTC (OR = 0.18, P = 0.001). Concerning the aggressive features of DTC, higher level of N stage was more likely to occur in subjects carrying the wild-type genotypes at rs1800860 site (for dominant model: OR = 0.48, P = 0.008). Another rare haplotype, CAAGCGT, conferred increased risk for the occurrence of distant metastasis (OR = 7.57, P = 0.009). Notably, higher thyroid stimulating hormone levels and lower parathyroid hormone levels were found in patients with rs2075912, rs2565200, and rs2742240 heterozygotes and rare homozygotes; similar results were observed between PTH levels and rs1800858.

Conclusion

This study provided useful information on RET variants that should be subjected to further study.

About the Authors:

Caiyun He

Contributed equally to this work with: Caiyun He, Jiangjun Ma, Yongle Jiang, Xuan Su

Roles Formal analysis, Funding acquisition, Resources, Software, Writing – original draft

¶‡ First authorship.

Affiliation: Department of Molecular Diagnostics, Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Jiangjun Ma

Contributed equally to this work with: Caiyun He, Jiangjun Ma, Yongle Jiang, Xuan Su

Roles Formal analysis, Software, Writing – original draft

¶‡ First authorship.

Affiliation: Department of Molecular Diagnostics, Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Yongle Jiang

Contributed equally to this work with: Caiyun He, Jiangjun Ma, Yongle Jiang, Xuan Su

Roles Data curation, Validation

¶‡ First authorship.

Affiliation: Sun Yat-Sen University, Guangzhou, China

Xuan Su

Contributed equally to this work with: Caiyun He, Jiangjun Ma, Yongle Jiang, Xuan Su

Roles Data curation, Formal analysis, Investigation, Validation, Writing – original draft

¶‡ First authorship.

Affiliation: Department of Head and Neck, Sun Yat-sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Xiao Zhang

Roles Data curation, Formal analysis

Affiliation: Department of Molecular Diagnostics, Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Weichao Chen

Roles Data curation, Investigation

Affiliation: Department of Head and Neck, Sun Yat-sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Zulu Ye

Roles Data curation

Affiliation: Department of Molecular Diagnostics, Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

Tiancheng Deng

Roles Data curation, Formal analysis

Affiliation: Guangzhou University of Chinese Medicine, Guangzhou, China

Wenze Deng

Roles Data curation, Formal analysis

Affiliation: Guangzhou University of Chinese Medicine, Guangzhou, China

Ankui Yang

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Head and Neck, Sun Yat-sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center for Cancer Medicine, Guangzhou, China

ORCID http://orcid.org/0000-0002-8534-3883

Introduction

Thyroid cancer, which is the most common malignancy in the endocrine system, is highly curable[1]. However, the increased incidence of differentiated thyroid carcinoma (DTC) is a serious problem[2]. Recent research has focused on the pathogenesis of this malignancy, with genetic factors being suggested to play an crucial role in the carcinogenesis of DTC[3]. Increasing efforts are being devoted to identifying susceptibility loci for this disease and its clinical phenotype. As the most common variants in the population, single nucleotide polymorphisms (SNPs) have the potential to be diagnostic, monitoring, and prognostic biomarkers for DTC.

The RET (rearranged during transfection) proto-oncogene has attracted much interest in DTC because aberrant activation of RET by somatic rearrangements or point mutations is a unique genetic event in DTC, which also plays a significant role in thyroid carcinogenesis[4, 5]. The protein encoded by RET functions as a multicomponent complex in which RET binds to a ligand belonging to the glial cell line–derived neurotrophic factor family[6]. In our previous work, we confirmed the association of RET rearrangements with DTC susceptibility and clinical phenotype in the Chinese population[7, 8]. To date, however, a systematic study of association between thecommon SNP of RET and susceptibility to DTC and its clinical phenotype has not been performed.

The levels of serum thyroid stimulating hormone (TSH), free triiodothyronine (fT3) and thyroxine (fT4), thyroglobulin (TG), and anti-thyroid peroxidase antibody (ATPO) are routine testing in the patients with thyroid cancer. These factors determine the balance of thyroid hormone homeostasis. Whether the common SNP in RET oncogene affects the thyroid function in DTC patients remains unclear and requires further study.

The RET gene is located at 10q11.2 and comprises 20 exons and 21 introns. To date, hundreds of SNPs in this gene have been revealed. However, information connecting RET SNPs with susceptibility to DTC is very limited. To assess comprehensively the role of RET genetic polymorphism, we adopted a candidate gene association study strategy with seven selected potentially functional tagging SNPs (tagSNPs) in RET[9]. We examined their individual and combined effects on the risk of DTC and their relationship with clinical severity and thyroid function.

Material and methods

Study subjects

This protocol was approved by the Ethics Committee of Sun Yat-sen University Cancer Center, Guangdong, China. Written informed consent was obtained from all patients at their first visit. For the genetic association study, a total of 552 subjects, comprising 300 DTC cases and 252 healthy control subjects, were included. The blood collection took place when the patients first visited this center. All the included subjects were recruited from Sun Yat-sen University Cancer Center in Guangzhou of Guangdong Province, China, between July 2011 and June 2016. We excluded two DTC patients from further analysis because of a past history of breast cancer or lung cancer. The healthy subjects comprised individuals with normal manifestation of the thyroid gland and thyroid function under ultrasound and serum examination. Healthy subjects who had a history of nodular goiter, Hashimoto's thyroiditis, or other benign thyroid diseases, and who had a history of other malignancies, were excluded.

Data collection

Medical records and surgical pathology reports were reviewed by Caiyun He, Yongle Jiang and Xuan Su. Information of demographic parameters, pathological characteristics of the tumors, and thyroid function parameters were collected. Slides were reviewed independently by two pathologists to confirm the diagnosis of DTC and concomitant thyroid benign diseases, such as nodular goiter and Hashimoto's thyroiditis. Information on the primary tumor size, stage grouping, extrathyroid extension, and metastasis were assessed based on the National Comprehensive Cancer Network (NCCN Guidelines, Version 1, 2016) on thyroid cancer recommendations (https://www.nccn.org/). The clinicopathological features of T, N and M stages and extrathyroid extension are well-known parameters indicating the clinical severity of thyroid cancer. The measurement and distribution of the clinicopathological features of thyroid cancer patients were showed in S1 Table.

TagSNP selection

HapMap genotype data of the Chinese Han Beijing (CHB) population (Release 27, Phase I +II +III, http://www.HapMap.org) were extracted within extended gene regions of RET in 2016, which encompassed 2 kb of upstream and downstream flanking sequence. The CHB RET gene has 33 common SNPs. As described in our previous studies, tagSNPs were selected by using the Haploview software[10, 11]. Briefly, tagSNPs were chosen initially based on pairwise linkage disequilibrium information to maximally represent (r2 > 0.8) common SNPs (minor allele frequency (MAF) > 0.05) by Haploview 4.2. We then used bioinformatics search to prioritize the tagSNPs for genotyping, based on their predicted functional effects (http://snpinfo.niehs.nih.gov/). For example, rs1799939 was predicted to be located at a splicing site, suggesting a potential effect on this gene's splicing regulation.

DNA isolation and genotyping of RET tagSNPs

For each subject, 3 mL of fasting blood was collected for DNA isolation. Genomic DNA was isolated from peripheral blood lymphocytes by the routine phenol–chloroform method, as previously described[10, 11]. Each DNA sample was diluted to a working concentration of 50 ng/μL for genotyping. Assay design and SNP genotyping were performed by Gene Company Limited (Shanghai, China) using the KASPar platform according to the manufacturer’s instructions[12]. All samples were randomized on 384-well plates and blinded for disease status. Twenty randomly selected samples were genotyped repeatedly, and the results were 100% concordant.

Statistical analysis

The Hardy-Weinberg Equilibrium (HWE) of the genotype distribution of each SNP was detected by a chi-square test in the control group. The odds ratio (OR) and corresponding 95% confidence interval (CI) were calculated to measure the strength of the association between the genotype and the risk of DTC. The genetic effect of a single tagSNP on disease risk was detected by Logistic regression analysis controlling for sex and age. For the haplotype analysis of the seven RET tagSNPs (rs17028, rs1799939, rs1800858, rs1800860, rs2075912, rs2565200, and rs2742240), we set the other haplotypes pooled together as a reference, and assessed the genetic effect of each haplotype with a frequency of at least more than 0.03 in the healthy controls.

For the association analysis between the genotype and clinicopathological features, we only performed association analysis in DTC patients since the control subjects have normal manifestation of thyroid. Logistic regression analysis adjusting for sex and age was employed. For the analysis of thyroid function, the difference of fT3, fT4, TSH, Anti-TG, TG, PTH and calcitonin levels between two groups were compared by the Mann–Whitney U test because these quantitative variables did not follow a normal distribution. The corresponding variables are shown as the median (with 25–75% quartiles).

All these analyses were performed using SPSS 17.0 software (SPSS, Chicago, IL, USA), except for the haplotype analysis, which was performed using the online software SHEsis[13] (http://analysis.bio-x.cn/myAnalysis.php). All P values were two sided, and a P value < 0.05 was considered statistically significant. The authenticity of this article has been validated by uploading the key raw data onto the Research Data Deposit public platform (www.researchdata.org.cn), with the approval RDD number as RDDB2017000172.

Results

Functional prediction of tagSNPs

Seven tagSNPs in RET with potential functions were included in this study and their predicted functional effects are summarized in Table 1. Among the seven tagSNPs, rs1799939, rs1800858, and rs1800860 are located within exon regions and were predicted to affect exonic splicing enhancer (ESE) or exonic splicing silencer (ESS) binding site activity, or even abolish a protein domain. The other four tagSNPs, rs17028, rs2075912, rs2565200, and rs2742240 are located in 3’ untranslated region (UTR) and were predicted to be located within microRNA (miRNA) binding sites. The distribution of the genotype frequencies of the seven tagSNPs in healthy controls obeyed the Hardy-Weinberg equilibrium test (Table 1, all P > 0.05). Therefore, these seven tagSNPs were included for the association analysis.

[Figure omitted. See PDF.]

Table 1. Functional prediction of seven tagSNPs of the RET gene.

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

Association between RET tagSNPs and susceptibility, clinicopathological features, and thyroid function of DTC

When we considered the influence of the RET tagSNPs on susceptibility to DTC, none of the seven tagSNPs showed an association with risk of DTC (Table 2, all P > 0.05). In the subgroup analysis, rs1799939 AG or AG plus AA genotypes showed an association with increased risk for the DTC cases without concomitant thyroid benign diseases (OR = 1.93 and 1.88, respectively, see S2 Table). No association was found between RET SNPs and the risk of DTC cases with concomitant thyroid benign diseases (see S3 Table).

[Figure omitted. See PDF.]

Table 2. Association between RET tagSNP and susceptibility to thyroid cancer.

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

We further analyzed the relationship between RET tagSNPs and the clinicopathological features of DTC (Table 3). When compared with the CC or CC plus TC genotypes, TT variants at rs17028 was boarderline associated with a reduced risk of developing a higher T stage exceeding T1 (for TT vs. CC: OR = 0.26, 95% CI: 0.08–0.85, P = 0.026; for TT vs. CC and TC: OR = 0.25, 95% CI: 0.08–0.78, P = 0.017), and seemed to associated with a reduced risk for the occurrence of extrathyroid extension (for TT vs. CC: OR = 0.34, 95% CI: 0.11–1.00, P = 0.050; for TT vs. CC and TC: OR = 0.34, 95% CI: 0.12–0.97, P = 0.044). The rs1800860 GA variant had a 0.45-fold reduced risk for the occurrence of lymph node metastasis compared with the common GG genotype (P = 0.005); a similar effect was observed for GA plus AA variants (OR = 0.48, 95% CI: 0.28–0.83, P = 0.008).

[Figure omitted. See PDF.]

Table 3. Association between RET tagSNP and clinicopathological features of thyroid cancer.

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

We also explored whether RET variants affect pre-operative thyroid function in DTC patients (Table 4). Variants in the rs2075912, rs2565200 and rs2742240 loci were associated statistically with increased thyroid stimulating hormone (TSH) levels but decreased parathyroid hormone (PTH) levels. The PTH levels in DTC cases carrying GA genotype at rs1800858 site were also statistically lower than those in the wild-type carriers. In addition, subjects carrying the rare rs17028 TT homozygote genotype showed a borderline association with a higher PTH levels compared with the common CC homozygotes.

[Figure omitted. See PDF.]

Table 4. Association between RET tagSNP and pre-operative thyroid function.

https://doi.org/10.1371/journal.pone.0187968.t004

Association between haplotype of RET and clinicopathological features and risk of DTC

For DTC susceptibility, a rare haplotype of CGGATAA of rs17028-rs1799939-rs1800858-rs1800860-rs2075912-rs2565200-rs2742240 showed a significant association with reduced risk of developing DTC, demonstrating an OR of 0.18 (95% CI: 0.06–0.54, P values for Pearson’s Chi-square test or Fisher’s exact test = 0.001, Table 5).

[Figure omitted. See PDF.]

Table 5. Association between haplotypes of RET tagSNPs and susceptibility to thyroid cancer.

https://doi.org/10.1371/journal.pone.0187968.t005

For the clinicopathological features of DTC, another rare haplotype CAAGCGT of rs17028-rs1799939-rs1800858-rs1800860-rs2075912-rs2565200-rs2742240 conferred increased risk for the occurrence of distant metastasis (OR = 7.57, 95% CI = 1.01–56.56, P values for Pearson’s Chi-square test or Fisher’s exact test = 0.009, Table 6). No associations were observed between haplotypes of the seven RET tagSNPs and T stage, N stage, extrathyroid extension, number of cancer lesions, and occurrence of concomitant thyroid benign diseases in DTC (all P > 0.05, data not shown).

[Figure omitted. See PDF.]

Table 6. Association between haplotypes of RET tagSNPs and the occurrence of distant metastasis in thyroid cancer.

https://doi.org/10.1371/journal.pone.0187968.t006

Discussion

Genetic variants of the proto-oncogene RET have attracted much research attention in recent studies of cancer causation[4, 5]. Our evaluation of a subgroup of DTC patients without concomitant thyroid benign diseases showed that DTC was more likely to occur in subjects with the rs1799939 AG or AG plus AA genotypes compared with the common GG genotype. Concerning the aggressive features of DTC, higher level of N stage was more likely to exist in subjects carrying the common homozygotes of rs1800860. Of note, higher TSH levels and lower PTH levels were found in patients with rs2075912, rs2565200, and rs2742240 heterozygotes and rare homozygotes, and similar results were observed between PTH levels and rs1800858.

We observed that the carriers with rs1799939 variants were more susceptible to DTC without concomitant thyroid benign diseases. rs1799939 is one of the most frequently studied SNPs in RET, and leads to an amino acid change from glycine to serine. In this study, online informatics software also predicted that rs1799939 would affect the splicing of RET transcripts. In another Chinese study, Huang and Yang reported that the rs1799939 AA genotype conferred a 3.76-fold increased risk of thyroid cancer[14]. A positive association was also reported in Khan et al.’s study in an Indian population[15]. However, a negative association for this SNP was found in Santos et al.’s study in a Portuguese Population[16]. Besides an independent effect of a single SNP, we observed an association between a rare haplotype, CGGATAA, (incorporated by rs17028-rs1799939-rs1800858-rs1800860-rs2075912-rs2565200-rs2742240 SNPs) and reduced risk of DTC, compared with other haplotypes, including the A risky allele at the rs1799939 site, which excluded chance as an explanation for the observed associations. The mechanisms by which the RET polymorphisms may confer an increased susceptibility to DTC remain to be determined. A pathogenic role of this SNP in DTC should be further evaluated in vitro.

We also found relationships between two RET tagSNPs and the aggressive features of DTC. The homozygous genotype of rs17028 showed hints of association with a reduced risk of the occurrence of a higher T stage, as well as extrathyroid extension, and the heterozygous genotype of rs1800860 contributed to a reduced risk for the occurrence of lymph node metastasis. rs17028 and rs1800860 have been reported in patients with Hirschsprung disease[17]. In addition, rs1800860 was reported to be associated with the risk of differentiated thyroid cancer[18]. This is an initial study reporting the frequencies of rs17028 and rs1800860 and their relationship with the clinicopathological features in DTC. As for the molecular function, a study focusing on kidney size in newborns suggested that rs1800860 was identified within an exonic splicing enhancer[19]. As predicted in this study, the variant A in the mRNA reduced the affinity for spliceosome proteins, enhanced the likelihood of aberrant mRNA splicing, and diminished the level of the functional transcript in human cells. So far, there is no known function for rs17028. The variant T allele of rs17028, located in the RET 3’ UTR, was predicted to impair the binding affinity of miR-24. He et al.’s study suggested overexpression of miR-24 in DTC[20]. Therefore, we hypothesized that the weakened binding affinity of miR-24 might constrain tumor progression in DTC patients harboring the T allele. Whether this SNP affects RET function in DTC in this manner requires further study.

In contrast to the susceptibility and clinicopathological features, RET variants have rarely been evaluated in thyroid function. Among the clinical parameters of thyroid function, the most frequently affected factors are TSH and PTH, which are critical factors for monitoring patients after surgery. Our previous study found that RET fusion genes correlated with higher TSH levels after surgery. In this study, we found several SNPs that were related to pre-operative thyroid function, including rs2075912, rs2565200, rs2742240, and rs1800858. Strikingly, these four SNPs did not show an association with susceptibility and clinical phenotype of DTC, whilst the DTC risk and clinical phenotype associated rs1799939 and rs1800860 were not related with thyroid function. This phenomenon provoked speculation that RET variants exclusively affect the development of DTC and thyroid function. Among the four SNPs, rs2075912, rs2565200, and rs2742240, are located in the RET 3’ UTR and were predicted to affect the binding of certain miRNAs. This suggested a biological function of the RET variants in the regulation of miRNAs that affect the synthesis of thyroid factors. Although the SNPs were not related to the susceptibility and TNM staging in DTC, it remains unclear whether these RET SNPs could predispose patients the recurrence of DTC, which requires further study.

It is a pity that there are several limitations in our study. In this study, the statistical analyses were performed on seven RET SNPs. And our study population is relatively small, especially for the haplotype analysis for distant metastasis in thyroid cancer; thereby, impact of RET SNPs on disease risk and disease severity might be probably underestimated. View in this manner, we considered a more conservative significance level of 0.05 rather than 0.05/7, as recommended previously[21, 22]. However, we are aware that the statistical results of our analyses should be considered to be exploratory, and hence did not adjust for multiple testing.

In conclusion, DTC without concomitant thyroid benign diseases was more likely to occur in subjects with the rs1799939 AG variant compared with those carrying the common GG genotype. An association was found between the rare haplotype CGGATAA and reduced risk of DTC compared with other haplotypes incorporating the risky A allele at rs1799939 site, which excluded chance as an explanation for the observed associations. A higher level of N stage was more likely to exist in subjects carrying the wild-type genotype at rs1800860. Notably, higher TSH levels and lower PTH levels were found in patients with rs2075912, rs2565200, and rs2742240 heterozygotes and rare homozygotes, and similar results were observed between PTH levels and the GA genotype of rs1800858.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Measurement and distribution of the clinicopathological features of thyroid cancer patients.

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

(DOCX)

S2 Table. Association between RET tagSNP and susceptibility to thyroid cancer in DTC patients without concomitant thyroid benign diseases.

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

(DOCX)

S3 Table. Association between RET tagSNP and susceptibility to thyroid cancer in DTC patients with concomitant diseases.

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

(DOCX)

Citation: He C, Ma J, Jiang Y, Su X, Zhang X, Chen W, et al. (2017) Associations between RET tagSNPs and their haplotypes and susceptibility, clinical severity, and thyroid function in patients with differentiated thyroid cancer. PLoS ONE12(11): e0187968. https://doi.org/10.1371/journal.pone.0187968

References

1. Haugen BR, Alexander EK, Bible KC, Doherty GM, Mandel SJ, Nikiforov YE, et al. 2015 American Thyroid Association Management Guidelines for Adult Patients with Thyroid Nodules and Differentiated Thyroid Cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid: official journal of the American Thyroid Association.2016;26(1):1–133.

2. Vigneri R, Malandrino P, Vigneri P. The changing epidemiology of thyroid cancer: why is incidence increasing? Current opinion in oncology.2015;27(1):1–7. pmid:25310641

3. Nikiforov YE, Nikiforova MN. Molecular genetics and diagnosis of thyroid cancer. Nature reviews Endocrinology.2011;7(10):569–80. pmid:21878896

4. Arighi E, Borrello MG, Sariola H. RET tyrosine kinase signaling in development and cancer. Cytokine & growth factor reviews.2005;16(4–5):441–67.

5. Nikiforov YE. RET/PTC rearrangement in thyroid tumors. Endocrine pathology.2002;13(1):3–16. pmid:12114746

6. Durbec P, Marcos-Gutierrez CV, Kilkenny C, Grigoriou M, Wartiowaara K, Suvanto P, et al. GDNF signalling through the Ret receptor tyrosine kinase. Nature.1996;381(6585):789–93. pmid:8657282

7. Su X, Li Z, He C, Chen W, Fu X, Yang A. Radiation exposure, young age, and female gender are associated with high prevalence of RET/PTC1 and RET/PTC3 in papillary thyroid cancer: a meta-analysis. Oncotarget.2016;7(13):16716–30. pmid:26918339

8. Su X, He C, Ma J, Tang T, Zhang X, Ye Z, et al. RET/PTC Rearrangements Are Associated with Elevated Postoperative TSH Levels and Multifocal Lesions in Papillary Thyroid Cancer without Concomitant Thyroid Benign Disease. PloS one.2016;11(11):e0165596. pmid:27802347

9. Tabor HK, Risch NJ, Myers RM. Candidate-gene approaches for studying complex genetic traits: practical considerations. Nature reviews Genetics.2002;3(5):391–7. pmid:11988764

10. He C, Tu H, Sun L, Xu Q, Li P, Gong Y, et al. Helicobacter pylori-related host gene polymorphisms associated with susceptibility of gastric carcinogenesis: a two-stage case-control study in Chinese. Carcinogenesis.2013;34(7):1450–7. pmid:23455381

11. He CY, Sun LP, Xu Q, Liu JW, Jiang JY, Dong NN, et al. PGC TagSNP and its interaction with H. pylori and relation with gene expression in susceptibility to gastric carcinogenesis. PloS one.2014;9(12):e115955. pmid:25551587

12. Wang B, Xu Q, Yang HW, Sun LP, Yuan Y. The association of six polymorphisms of five genes involved in three steps of nucleotide excision repair pathways with hepatocellular cancer risk. Oncotarget.2016;7(15):20357–67. pmid:26967386

13. Shi YY, He L. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell research.2005;15(2):97–8. pmid:15740637

14. Huang RX, Yang F. RET polymorphisms might be the risk factors for thyroid cancer. International journal of clinical and experimental pathology.2015;8(5):5793–7. pmid:26191299

15. Khan MS, Pandith AA, Iqbal M, Naykoo NA, Khan SH, Rather TA, et al. Possible Impact of RET Polymorphism and Its Haplotypic Association Modulates the Susceptibility to Thyroid Cancer. Journal of cellular biochemistry.2015;116(8):1712–8. pmid:25736215

16. Santos M, Azevedo T, Martins T, Rodrigues FJ, Lemos MC. Association of RET genetic polymorphisms and haplotypes with papillary thyroid carcinoma in the Portuguese population: a case-control study. PloS one.2014;9(10):e109822. pmid:25330015

17. Pan ZW, Luo CF, Liu ZJ, Li JC. RET 3'UTR polymorphisms and its protective role in Hirschsprung disease in southeastern Chinese. Journal of pediatric surgery.2012;47(9):1699–705. pmid:22974609

18. Ho T, Li G, Zhao C, Wei Q, Sturgis EM. RET polymorphisms and haplotypes and risk of differentiated thyroid cancer. The Laryngoscope.2005;115(6):1035–41. pmid:15933516

19. Zhang Z, Quinlan J, Hoy W, Hughson MD, Lemire M, Hudson T, et al. A common RET variant is associated with reduced newborn kidney size and function. Journal of the American Society of Nephrology: JASN.2008;19(10):2027–34. pmid:18820179

20. He H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S, et al. The role of microRNA genes in papillary thyroid carcinoma. Proceedings of the National Academy of Sciences of the United States of America.2005;102(52):19075–80. pmid:16365291

21. Bender R, Lange S. Adjusting for multiple testing—when and how? J Clin Epidemiol.2001;54(4):343–9. pmid:11297884

22. Labovitz S. Criteria for Selecting a Significance Level: A Note on the Sacredness of .05. The American Sociologist 1968;3(3): 220–2.