We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
A new panel of SNPs to assess thyroid carcinoma risk: a pilot study in a Brazilian admixture population
BMC Medical Geneticsvolume 18, Article number: 140 (2017)
Thyroid cancer is a common malignant disease of the endocrine system with increasing incidence rates over the last few decades. In this study, we sought to analyze the possible association of 45 single nucleotide polymorphisms (SNPs) with thyroid cancer in a population from Rio Grande do Norte, Brazil.
Based on histological analysis by a pathologist, 80 normal thyroid specimens of tissue adjacent to thyroid tumors were obtained from the biobank at the Laboratory of Pathology of Liga Norte Riograndense Contra o Câncer, Natal, RN. Patient samples were then genotyped using the MassARRAY platform (Sequenon, Inc) followed by statistical analysis employing the SNPassoc package in R program. The genotypic frequencies of all 45 SNPs obtained from the International HapMap Project database and based on data from the ancestral populations of European and African origin were used to compose the control study group.
In our study, the following 9 SNPs showed significant differences in their frequency when comparing the study and control groups: rs3744962, rs258107, rs1461855, rs4075022, rs9943744, rs4075570, rs2356508, rs17485896, and rs2651339. Furthermore, the SNPs rs374492 C/T and rs258107 C/T were associated with a relative risk for thyroid carcinoma of 3.78 (p = 6.27 × 10e−5) and 2.91 (p = 8.27 × 10e−5), respectively, after Bonferroni’s correction for multiple comparisons.
These nine polymorphisms could be potential biomarkers of predisposition to thyroid carcinoma in the population from Rio Grande do Norte. However, complementary studies including a control group with samples obtained from healthy subjects in Rio Grande do Norte state, should be conducted to confirm these results.
Thyroid cancer (TC) is the most common endocrine malignancy, and its incidence rate has been increasing noticeably for decades . Hence, if the observed trends are maintained, TC will replace colorectal cancer as the fourth leading cancer diagnosis in the United States by 2030 . Such an increase is likely due to improved diagnosis, and largely or completely reflects the over diagnosis of indolent disease . According to the Brazilian National Cancer Institute, about 6.960 new cases of TC were expected in Brazil in 2016 . The state of Rio Grande do Norte, located in the northeast region of Brazil, presents the fourth highest estimated TC incidence rate in the country, which corresponds to 130 cases per 100,000 women and 30 cases per 100,000 men .
The thyroid gland is mainly composed of two cell types: follicular cells and parafollicular cells. Based on histological and clinical parameters, follicular cell-derived carcinomas, which correspond to more than 90% of thyroid carcinomas, are typically divided into well differentiated, poorly differentiated, and undifferentiated (or anaplastic) carcinomas. Well differentiated thyroid carcinomas (WDTC) are further subdivided histologically as papillary thyroid carcinoma and follicular thyroid carcinoma, which correspond to 80% and 10% of thyroid carcinomas, respectively . Both these carcinomas can progress to poorly differentiated carcinoma or completely lose differentiation, giving rise to anaplastic carcinoma . On the other hand, parafollicular cell-derived medullary carcinoma accounts for about 5% of thyroid carcinomas [5, 7].
The recent progress in understanding the molecular pathogenesis of TC has paved the way for developing more effective treatment strategies. This has mainly resulted from identification of numerous genetic and epigenetic changes including mutations, gene copy-number gain, gene translocation, and aberrant gene methylation, resulting in alterations of the signaling pathways involved in regulating cell proliferation and survival, such as the MAPK and PI3K–AKT pathways among others, which are reshaping TC medicine [8, 9]. In addition, it is widely accepted that genetic predisposition to WDTC is expected to have common low-penetrance and rare moderate-penetrance genetic variants interacting with each other and with the environment determining individual susceptibility .
Recent high-throughput genotyping (GWAS) studies sought to identify single nucleotide polymorphisms (SNPs) associated with increased risk of thyroid tumorigenesis and to better explain the role of genetic variations in predisposition to TC [10,11,12,13,14,15,16,17,18,19,20,21,22]. However, some studies have demonstrated contradictory results, probably due to genetic diversity both within and among populations. This observation raises the question of whether a new panel of SNPs could be used as potential predisposition markers for TC.
Therefore, this study aimed to investigate whether the selected SNPs could be associated with predisposition to TC in patients diagnosed with this neoplasia in the state of Rio Grande do Norte. Considering that the population of this region is characterized by a marked ethnic mixture, this study may significantly contribute to elucidating the molecular basis underlying both predisposition to TC and the effect of interbreed populations on SNP-based association studies.
In this retrospective study, 80 tumor formalin-fixed paraffin-embedded (FFPE) samples of tissue adjacent to thyroid tumors were selected from patients diagnosed with thyroid cancer at the Laboratory of Pathology, Hospital Liga Norte Riograndense Contra o Câncer (Natal, Brazil) between 2003 and 2011. Among these 80 cases, 77 had papillary thyroid carcinoma and 3 had follicular thyroid carcinoma, as confirmed based on histological analysis by a pathologist. Ten 20-μm sections were cut from blocks containing non-tumor tissue using a Rotary Microtome YD-335 (ANCAP, São Paulo, Brazil). Samples were placed in 1.5 mL microtubes at room temperature until processing. For each patient sample, a unique High Profile 818 disposable blade (Leica Microsystems GmbH, Wetzlar, Germany) was used to avoid contamination of genomic DNA.
The control group was obtained from the International HapMap project and consisted of genotype data from 180 healthy individuals: 90 with European ancestry and 90 with African ancestry (HapMap Genome Browse, release#24 Phase 1 & 2 - full dataset).
This research was approved by the Ethics Committee in Research from Hospital Liga Norte Riograndense Contra o Câncer (protocol number 558788).
FFPE tissue samples were deparaffinized with 2.4 mL of xylene (2 × 5 min at 14,000 rpm) and washed with 2.4 mL of 100% ethanol (2 × 5 min at 14,000 rpm). DNA was then isolated using the QIAamp® DNA FFPE Tissue Kit (QIAGEN, Valencia, EUA), according to the manufacturer’s instruction. DNA quantification was performed using the NanoDrop ND-1000 spectrophotometer. DNA samples were diluted to a final concentration of 20 ng/μL.
A panel of 45 SNPs was selected through previous studies performed at the Genotyping National Center (CEGEN) in Santiago de Compostela – Spain (www.usc.es/cegen/), coordinated by Dr. Angel Carracedo, Professor of the Faculty of Medicine, University of Santiago de Compostela (USC) .
SNP genotyping of the case group was performed using the MassARRAY SNP genotyping system (Agena Bioscience San Diego, EUA) according to manufacturer’s instructions at the National Genotyping Center (CEGEN), on a panel of 45 SNP assays. The primers for amplification and extension were designed using the Extend Primer Assay Design software v4. Sequenom iPLEX GOLD chemistry was used for locus-specific amplification, followed by a single-base primer extension reaction, which generated products of different masses that were quantitatively analyzed using MALDI-TOF mass spectrometry. The resulting data were analyzed using TyperAnalyzer software v 4, followed by manual inspection of the spectra by trained personnel . All assays were performed in 384-well plates, including negative controls and a trio of Coriell samples (Na10830, Na10831, and Na12147) for quality control.
Statistical tests were performed using the SNPassoc function in the R software v2.14.1 statistical package (R Development Core Team, 2006) . For each SNP, associations were assessed by applying logistic regression to estimate the odds ratios (ORs) with 95% confidence intervals and P values. A P value <0.05 was considered statistically significant.
Eighty TC patients from Brazil were genotyped for 45 SNPs, to identify potential molecular markers of predisposition to this neoplasia. Among all the patients studied, 77 were diagnosed with papillary thyroid carcinoma and 3 with follicular thyroid carcinoma. The age of one patient was not determined, and the staging, location, and size of tumor were also not determined in some patients. Among the 80 patients, 61 lymph node metastases were not evaluated. These data are shown in Table 1.
The allelic frequencies of 45 SNPs were analyzed and the minimum cutoff value for missing genotypes was defined as 10%. The allelic composition, allele major frequency and the p value for the Hardy-Weinberg equilibrium (HWE) of each SNP are shown in Table 2. Of these 45 SNPs, 11 were excluded from the study because they were not in the HWE (p < 0.05). The allelic frequencies from the Brazilian population were compared with those available in the HapMap database for 180 individuals of European and African ancestry, which constituted the control group.
The association analysis between patients and control groups for SNPs that showed a p-value greater than 0.05 for Hardy-Weinberg equilibrium is shown in Table 3 and Additional file 1: Table S1, which classify SNPs according to the inheritance model: dominant (AA/Ab + bb); recessive (AA + Ab/bb); over dominant (Ab/AA + bb) and co-dominant (AA/Ab/bb). As observed, 7 SNPs showed significant expression in the co-dominant model (rs1461855, rs2356508, rs258107, rs3744962, rs4075022, rs4075570, and rs9943744), 5 in the dominant model (rs1461855, rs258107, rs3744962, rs4075022, and rs9943744), 4 in the recessive model (rs1461855, rs2651339, rs4075570, and rs9943744) and 6 in the over dominant model (rs1461855, rs17485896, rs2356508, rs258107, rs3744962, and rs4075022). Most SNPs presented statistical significance under more than one model of inheritance. Furthermore, the SNPs rs374492 C/T and rs258107 C/T were associated with a relative risk for TC of 3.78 (p = 6.27 × 10e−5) and 2.91 (p = 8.27 × 10e−5), respectively after Bonferroni’s correction for multiple comparisons.
Figure 1 shows the results for logistic regression analysis of significant SNPs according to the model of inheritance, their respective odds ratios (OR), and p-values. Among the 9 SNPs analyzed (rs1461855, rs17485896, rs2356508, rs258107, rs2651339, rs3744962, rs4075022, rs4075570, and rs9943744), 7 showed OR values greater than 1, meaning that those alleles represent a risk for thyroid cancer development. In particular, the SNP rs3744962 showed the greatest risk with an OR value of 3.78. On the other hand, 2 SNPs (rs1461855 and rs2356508) were associated with a protective effect, with OR values less than 1.
The genetic ancestry profile of Brazilian population samples have been largely investigated in the last decades, mainly using uniparental inherited markers such mitochondrial DNA and Y-chromosome polymorphisms [25,26,27,28,29], as well as autosomal STR markers [30, 31]. In addition, some studies used ancestry informative markers such as SNPs and insertion/deletion polymorphisms have been used to estimate genetic ancestry [32,33,34,35,36].
All these previous studies conclude that the genetic ancestry of the Brazilian population is widely heterogeneous and is characterized by extensive admixture from three different ancestral roots: Amerindians, Europeans, and Africans. Furthermore, the composition of the Brazilian population ancestry follows a clear trend related to historical facts [37, 38]. European colonization of Brazil showed successive migration waves, starting with the arrival of around 500,000 Portuguese men between the years 1500 and 1800, when they met the local Amerindian population. Portuguese-Amerindian admixture began soon after arrival of the first colonizers. Simultaneously, about 4 million Africans were compulsorily introduced into Brazil as slaves until the mid-eighteenth century, providing the second major ethnic contribution to the Brazilian population. Later, other waves of immigration to Brazil occurred, mainly from Italy, Portugal, Spain, Germany, Syria, Lebanon, and Japan, adding even more complexity to the already multi-ethnic highly admixed Brazilian population [37, 38].
Brazil is the fifth largest country in the world and is divided into five geographical regions. In all these regions, the European ancestry is predominant with proportions ranging from 60.6% in the Northeast to 77.7% in the South [34, 39,40,41,42]. The African and Amerindian contributions to the population genetic background vary according to the region. The African ancestry proportion was the second largest in the Northeast (30.3%), followed in decreasing order by the Southeast, South, and North. On the other hand, populations in the North consist of a significant proportion of Amerindian ancestry, while remaining relatively uniform in other regions. In particular, the northeast region of Brazil where the state of Rio Grande do Norte is located, presents the highest African ancestry contribution in the country ranging from 18.6 to 56.8% [34, 39,40,41,42,43].
The heterogeneity and admixture of Brazilians constitutes an important issue to be considered in genetic association studies . Besides, it has important clinical implications for the design and interpretation of clinical trials and genetic counseling. Because of its heterogeneous Amerindian, European, and African ancestral roots, the Brazilian population has been considered an important model for population genetics. In this context, our work aimed to genotype and assess the allele frequency of 45 SNPs to identify potential molecular markers of predisposition to TC in the population of Rio Grande do Norte, compared to the frequency of these alleles in European and African populations.
An extensive review of literature showed that considerable progress towards understanding complex diseases has been made in recent years due to development of high-throughput genotyping technologies . Furthermore, several GWAS studies have identified SNPs associated with many complex diseases or traits. In the particular case of TC, genetic risk loci are under characterization in terms of allelic variants or genes responsible for association with the disease, as well as the affected biological pathway and the main cell type driving its pathology .
In this regard, a growing number of studies have identified several SNPs associated with TC risk [10,11,12,13,14,15,16,17, 19, 20, 22, 46,47,48,49]. In the study by Gudmundsson and collaborators , a GWAS analysis was performed in a population from Iceland and showed a strong association of SNPs rs965513 and rs944289 with papillary and follicular thyroid cancer. Furthermore, rs944289 and rs965513 correspond to gene variants of Thyroid transcription factor 1 and 2, respectively, and both TC risk alleles were associated with low concentrations of thyroid stimulating hormone, and SNP rs965513 was associated with low thyroxin concentration and high triiodothyronine concentration . Later, the same group found that rs966423, rs2439302, and rs116909374 variants were also associated with TC .
The association of rs944289 and rs965513 with TC risk was also shown in genetic studies conducted in Asian populations [14, 50, 51]. Additionally, Wang and collaborators  observed that SNPs rs966423 and rs2439302, previously reported to be associated with TC in European populations, were confirmed as risk factors in the Chinese population.
Similarly, another study has demonstrated that SNPs rs2910164, rs965513, rs1867277, rs6983267, and rs944289 showed a significant association with TC risk in patients from the United Kingdom . Furthermore, according to Jendrzejewski , the SNP rs944289 can predispose to papillary thyroid carcinoma through deregulation of PTCSC3 expression, which acts as a tumor suppressor.
A correlation study between Forkhead box E1 (FOXE1) gene variants rs894673, rs1867277, and rs3758249, and histopathological features of TC, suggested that FOXE1 variations generate a higher risk for poor histopathological features in papillary thyroid carcinoma . In addition, Mancikova et al.,  showed an association between SNPs rs7028661 and rs7037324, located near the FOXE1 locus, and TC risk. Moreover, the rare alleles of three SNPs (rs2997312, rs10788123 and rs1254167) showed suggestive evidence of association with higher risk for the disease. On the other hand, the SNP rs4075570 conferred protection in the series studied.
Finally, Lidral  showed that rs7850258 G allele associated with cleft lip, cleft palate, and hypothyroidism, has significantly greater enhancer activity than the allele associated with thyroid cancer (A); and Wokolorczyk et al.,  showed that the SNP rs6983267 could be a good candidate multi-cancer susceptibility marker, once this SNP is associated with a wide range of cancers affecting the colon, prostate, breast, bladder, larynx, lung, kidney, and the thyroid . These are some examples of the rapidly increasing number of studies focused on identifying genetic variants associated with TC predisposition and development.
An important challenge faced by molecular epidemiological association studies of candidate disease-susceptibility genes is to define variants that are functionally implicated in the disease. In addition to identifying polymorphic genetic variants and their possible association with diseases, it is necessary to understand the functional relevance of these SNPs. This is particularly urgent because the amount of genomic information that is available greatly exceeds the information about the function of variants that are presented in human disease studies . Importantly, most of the recent genetic studies only establish the statistical associations of genetic markers and the disease, without supporting evidence of functional relevance. In context of predisposition to TC, new light is given on the molecular mechanisms of genetic variants through advances in molecular technologies.
As an example, the association of polymorphisms in DNA repair genes XRCC1 (rs25487, rs1799782) and XRCC3 (rs861539) with thyroid cancer risk and progression can be considered. As Yan and collaborators showed, the XRCC1 variant can interact with the XRCC3 variant to significantly increase differentiated thyroid carcinoma (DTC) susceptibility . Similarly, a study to assess the role of polymorphisms in the LEP (rs7799039 and rs2167270) and LEPR (rs1137101 and rs1137100) genes in DTC susceptibility and their effect on leptin levels showed that LEP polymorphisms modify serum leptin concentrations in patients with DTC. Furthermore, the polymorphisms rs7799039 and rs1137101 increase the risk of DTC development, though they do not correlate with tumor aggressiveness . As another example, supporting evidence of functional relevance was presented by Ceolin and collaborators , who evaluated the frequency of RET (REarranged during Transfection) proto-oncogene 3’UTR variants (rs76759170 and rs3026785) in Medullary Thyroid Carcinoma (MTC) patients. In silico analysis indicated that these variants might affect the secondary structure of RET mRNA, suggesting that they might play a role in the posttranscriptional control of RET transcripts.
In our study for molecular markers of predisposition to TC, we assessed the allele and genotype frequency of 45 SNPs, including four of the aforementioned SNPs. We observed that SNPs rs2910164 and rs965513, despite being in HWE, did not show associated risk with TC for the Brazilian population. Furthermore, SNPs rs6983267 and rs944289 were not in Hardy-Weinberg equilibrium in our studied groups. As observed, the frequencies of these particular SNPs differ from those observed for European, American and Asian populations, where some of these markers were associated with an increased risk of developing TC. These results reinforce the concept that the genetic constitution of a population, and the contribution of ancestral roots, is an important and influencing parameter in genetic association studies.
Other studies evaluated potential markers for predisposition to thyroid cancer in larger sample sizes [10, 46]. These studies were performed in European populations, but these kinds of studies were not previously performed in a Brazilian population. Therefore, this pilot study presents a good opportunity to examine the feasibility of our approach, exploring a novel panel of SNPs that could potentially predict TC risk. We have thus first investigated these SNPs in a pilot study with a plan to validate our findings in a larger sample in the future. It is critical to improve both our knowledge of thyroid cancer risk factors and our knowledge of markers that predict aggressive disease in order to reduce disease incidence and unnecessary treatments that result in undesirable side effects with long-term financial and clinical impacts.
Of a total of 35 SNPs that were found in HWE in the Brazilian population, logistic regression analysis of 7 SNPs showed statistical significance when evaluated according to the genetic model of inheritance. For the first time, our results suggest that a new panel constituted of SNPs rs3744962, rs258107, rs4075022, rs9943744, rs4075570, rs17485896, and rs2651339 could represent predisposition to TC development in the population of Rio Grande do Norte. Furthermore, they may be considered suitable molecular markers for early diagnosis of the disease. These SNPs could be useful models to predict risk and for genetic counseling in clinical practice, thus avoiding invasive methods of analysis and often inconclusive diagnoses. In addition, SNPs rs1461855 and rs2356508 might possibly be associated with a protective effect against TC development. Nevertheless, complementary studies with a larger patient population as well as a control group composed of samples from the state of Rio Grande do Norte should be conducted to confirm these results.
Differentiated Thyroid Carcinoma
Forkhead box E1
Allele Minor Frequency
Medullary Thyroid Carcinoma
REarranged during Transfection
Well differentiated thyroid carcinomas
Jin S, et al. Signaling pathways in thyroid cancer and their therapeutic implications. J Clin Med Res. 2016;8(4):284–96.
Rahib L, et al. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. 2014;74(11):2913–21.
La Vecchia C, et al. Thyroid cancer mortality and incidence: a global overview. Int J Cancer. 2015;136(9):2187–95.
INCA, I.N.d.C.J.A.G.d.S. 2016: http://www.inca.gov.br/estimativa/2016.
Pallante P, et al. Deregulation of microRNA expression in thyroid neoplasias. Nat Rev Endocrinol. 2014;10(2):88–101.
Leonardi, G. C. et al. microRNAs and thyroid cancer : Biological and clinical significance. Int J Mol Med, 2012. p. 991–999.
Yuan ZM, Yang ZL, Zheng Q. Deregulation of microRNA expression in thyroid tumors. J Zhejiang Univ Sci B. 2014;15(3):212–24.
Xing M. Molecular pathogenesis and mechanisms of thyroid cancer. Nat Rev Cancer. 2013;13(3):184–99.
Croce CM. Oncogenes and cancer. N Engl J Med. 2008;358(5):502–11.
Mancikova V, et al. Thyroid cancer GWAS identifies 10q26.12 and 6q14.1 as novel susceptibility loci and reveals genetic heterogeneity among populations. Int J Cancer. 2015;137(8):1870–8.
Gudmundsson J, et al. Discovery of common variants associated with low TSH levels and thyroid cancer risk. Nat Genet. 2012;44(3):319–22.
Jones AM, 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(3):158–63.
Figlioli G, et al. Novel genome-wide association study-based candidate loci for differentiated thyroid cancer risk. J Clin Endocrinol Metab. 2014;99(10):E2084–92.
Gao Y, et al. Replication and meta-analysis of common gene mutations in TTF1 and TTF2 with papillary thyroid cancer. Medicine (Baltimore). 2015;94(36):e1246.
Lidral AC, et al. A single nucleotide polymorphism associated with isolated cleft lip and palate, thyroid cancer and hypothyroidism alters the activity of an oral epithelium and thyroid enhancer near FOXE1. Hum Mol Genet. 2015;24(14):3895–907.
Marcello MA, et al. Polymorphism in LEP and LEPR may modify Leptin levels and represent risk factors for thyroid cancer. Int J Endocrinol. 2015;2015:173218.
Somuncu E, et al. The investigation of foxe1 variations in papillary thyroid carcinoma. Int J Clin Exp Pathol. 2015;8(10):13458–64.
Wang Y, et al. Association between genetic polymorphisms in the promoter regions of Let-7 and risk of papillary thyroid carcinoma: a case-control study. Medicine (Baltimore). 2015;94(43):e1879.
Ceolin L, et al. Effect of 3'UTR RET variants on RET mRNA secondary structure and disease presentation in Medullary thyroid carcinoma. PLoS One. 2016;11(2):e0147840.
Figlioli G, et al. A comprehensive meta-analysis of case-control association studies to evaluate polymorphisms associated with the risk of differentiated thyroid carcinoma. Cancer Epidemiol Biomark Prev. 2016;
Liyanarachchi S, et al. Cumulative risk impact of five genetic variants associated with papillary thyroid carcinoma. Thyroid. 2013;23(12):1532–40.
Yan L, et al. Association studies between XRCC1, XRCC2, XRCC3 polymorphisms and differentiated thyroid carcinoma. Cell Physiol Biochem. 2016;38(3):1075–84.
Gabriel, S., L. Ziaugra, and D. Tabbaa, SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Protoc Hum Genet, 2009. Chapter 2: Unit 2. 12.
Gonzalez JR, et al. SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007;23(5):644–5.
Abe-Sandes K, Silva WA Jr, Zago MA. Heterogeneity of the Y chromosome in afro-Brazilian populations. Hum Biol. 2004;76(1):77–86.
Alves-Silva J, et al. The ancestry of Brazilian mtDNA lineages. Am J Hum Genet. 2000;67(2):444–61.
Bortolini MC, et al. African-derived south American populations: a history of symmetrical and asymmetrical matings according to sex revealed by bi- and uni-parental genetic markers. Am J Hum Biol. 1999;11(4):551–63.
Carvalho-Silva DR, et al. The phylogeography of Brazilian Y-chromosome lineages. Am J Hum Genet. 2001;68(1):281–6.
Silva WA, et al. MtDNA haplogroup analysis of black Brazilian and sub-Saharan populations: implications for the Atlantic slave trade. Hum Biol. 2006;78(1):29–41.
Callegari-Jacques SM, et al. Historical genetics: spatiotemporal analysis of the formation of the Brazilian population. Am J Hum Biol. 2003;15(6):824–34.
Pimenta JR, et al. Color and genomic ancestry in Brazilians: a study with forensic microsatellites. Hum Hered. 2006;62(4):190–5.
Blanton RE, et al. Genetic ancestry and income are associated with dengue hemorrhagic fever in a highly admixed population. Eur J Hum Genet. 2008;16(6):762–5.
Estrela RC, et al. Distribution of ABCB1 polymorphisms among Brazilians: impact of population admixture. Pharmacogenomics. 2008;9(3):267–76.
Lins TC, et al. Genetic composition of Brazilian population samples based on a set of twenty-eight ancestry informative SNPs. Am J Hum Biol. 2010;22(2):187–92.
Luizon MR, et al. Ancestry informative markers in Amerindians from Brazilian Amazon. Am J Hum Biol. 2008;20(1):86–90.
Suarez-Kurtz G, et al. Self-reported skin color, genomic ancestry and the distribution of GST polymorphisms. Pharmacogenet Genomics. 2007;17(9):765–71.
Aguiar VR, et al. Updated Brazilian STR allele frequency data using over 100,000 individuals: an analysis of CSF1PO, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D21S11, FGA, Penta D, Penta E, TH01, TPOX and vWA loci. Forensic Sci Int Genet. 2012;6(4):504–9.
Cordeiro Q, et al. A review of psychiatric genetics research in the Brazilian population. Rev Bras Psiquiatr. 2009;31(2):154–62.
Saloum de Neves Manta F, et al. Revisiting the genetic ancestry of Brazilians using autosomal AIM-Indels. PLoS One. 2013;8(9):e75145.
Pena SD, et al. The genomic ancestry of individuals from different geographical regions of Brazil is more uniform than expected. PLoS One. 2011;6(2):e17063.
Moura RR, et al. Meta-analysis of Brazilian genetic admixture and comparison with other Latin America countries. Am J Hum Biol. 2015;27(5):674–80.
Santos HC, et al. A minimum set of ancestry informative markers for determining admixture proportions in a mixed American population: the Brazilian set. Eur J Hum Genet. 2016;24(5):725–31.
Gontijo CC, et al. Brazilian quilombos: a repository of Amerindian alleles. Am J Hum Biol. 2014;26(2):142–50.
Tak YG, Farnham PJ. Making sense of GWAS: using epigenomics and genome engineering to understand the functional relevance of SNPs in non-coding regions of the human genome. Epigenetics Chromatin. 2015;8:57.
Ray-Jones H, et al. One SNP at a time: moving beyond GWAS in psoriasis. J Invest Dermatol. 2016;136(3):567–73.
Gudmundsson J, et al. Common variants on 9q22.33 and 14q13.3 predispose to thyroid cancer in European populations. Nat Genet. 2009;41(4):460–4.
Wokolorczyk D, et al. A range of cancers is associated with the rs6983267 marker on chromosome 8. Cancer Res. 2008;68(23):9982–6.
Wang YL, et al. 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(10):689–95.
Wei WJ, et al. Clinical significance of papillary thyroid cancer risk loci identified by genome-wide association studies. Cancer Gene Ther. 2015;208(3):68–75.
Ai L, et al. Associations between rs965513/rs944289 and papillary thyroid carcinoma risk: a meta-analysis. Endocrine. 2014;47(2):428–34.
Ai L, et al. Are the SNPs of NKX2-1 associated with papillary thyroid carcinoma in the Han population of northern China? Front Med. 2014;8(1):113–7.
Jendrzejewski J, et al. The polymorphism rs944289 predisposes to papillary thyroid carcinoma through a large intergenic noncoding RNA gene of tumor suppressor type. Proc Natl Acad Sci U S A. 2012;109(22):8646–51.
Rebbeck TR, Spitz M, Wu X. Assessing the function of genetic variants in candidate gene association studies. Nat Rev Genet. 2004;5(8):589–97.
The authors thanks Soledad Otero, Juan Ansede and Rosanna Abal for all the assistance during the processing and genotyping of the samples.
This study had financial support from FAPERN (EDITAL PPPIV 005/2011) and CNPq for PIBIC/IC fellowship.
Availability of data and materials
-The datasets generated and/or analyzed during the current study are available in the HapMap Genome Browser, release#24 Phase 1 & 2 full dataset repository, [http://www.internationalgenome.org/1000-genomes-browsers].
-The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
As this is a retrospective study where biopsy samples were obtained from 2003 to 2011, the Ethics Committee in Research from Hospital Liga Norte Riograndense Contra o Câncer approved the conduction of the research without informed consent of the participants for this particular case (Protocol number 558788). Nevertheless, the authors protected the privacy and maintain the confidentiality.
Consent for publication
Competing interests. The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The logistic regression analyses between case-control and SNPs polymorphism. (DOCX 129 kb)