DTC susceptibility loci evidenced in family studies on non-syndromic NMTC (in chronological order of discovery).
Abstract
Differentiated thyroid carcinoma (DTC) represents more than 90% of all thyroid cancer histological types. Its incidence has increased at a faster rate than most other malignancies during the last three decades and varies considerably around the world. The familial form of the disease has also become more common than previously reported, accounting for 5−15% of DTC cases. The main established risk factor of thyroid cancer is exposure to ionizing radiation, particularly if occurred during childhood. Thyroid cancer (including DTC) is also characterized by having one of the highest familial risks of any cancer supporting heritable predisposition. In spite of such a high familial risk, linkage analysis in non-syndromic DTC families (i.e. families where DTC is the primary cancer) performed two decades ago mapped several susceptibility loci but did not lead to the identification of high-penetrance causal germline variants. More recently, genome-wide association studies based on population case–control studies identified a limited number of DTC-associated loci and suggested that multiple low penetrance genes are involved in predisposition to DTC. This chapter reviews known genetic factors predisposing to DTC as well as approaches used to map them in various populations, and opens up on alternative strategies that could help to understand DTC tumorigenesis.
Keywords
- differentiated thyroid cancer
- familial risk
- case–control study
- genetic predisposition
- genome-wide association study
1. Introduction
Almost 95% of patients with malignant thyroid tumors have non-medullary thyroid cancer (NMTC), which originates from follicular cells. Most NMTC are well-differentiated and include papillary thyroid carcinoma (PTC) (80−90%) and follicular thyroid carcinoma (FTC) (10−15%) while poorly differentiated thyroid carcinoma and anaplastic carcinoma are rare. C-cell-derived medullary thyroid carcinoma (MTC) are diagnosed in approximately 5% of patients [1]. Altogether differentiated thyroid carcinoma (DTC) represents more than 90% of all thyroid cancer histological types and the most frequent malignancy of the endocrine system (Figure 1) [2].
DTC incidence varies considerably around the world. High incidence was reported in some Pacific islands such as Hawaii, New Caledonia, and French Polynesia [3]. Ethnic differences in incidence within the same country have been noted in Hawaii and New Caledonia with higher rates among Filipinos and Melanesians than in other ethnic groups [4]. The causes underlying these wide geographic and ethnic variations are still poorly understood. If this variation in incidence was attributed to screening practices, it was suggested that environment and inherited genetic risk factors may also play an important role [5]. Clinical awareness of potential risk factors, such as inherited genetic variants allows for earlier recognition of more vulnerable populations, earlier detection, proper treatment, and improved outcomes for patients and their families, justifying current efforts to identify and understand the causal factors and mechanisms of DTC so that effective interventions can be implemented.
Familial associations are often quantified in terms of familial relative risk (FRR). The FRR denotes the risk of disease when a family member is affected compared to the risk level in the general population. Specific types of familial relationships, like first-degree relatives, parent–child, or siblings are generally examined. Remarkably, thyroid cancer displays one of the strongest FRRs among cancers. Large case–control studies conducted in populations from Utah and Sweden showed that the FRR of thyroid cancer for first-degree relatives of probands was 8.5 and 12.4, respectively [6, 7]. Data from studies focusing on DTC conducted in Sweden, Iceland, and Norway showed that the standard incidence risk (SIR) of DTC was between 4.1 and 7.8 for male relatives and between 1.9 and 4.9 for female relatives of the proband [8, 9, 10]. The SIR for PTC was calculated in the Norwegian study as 5.8 and 4.1 for male and female relatives, respectively [10]. The family structures of non-medullary thyroid carcinoma (NMTC) patients in Taiwan were also studied. The prevalence of NMTC in the general population and in first-degree relatives of NMTC patients was 0.16% and 0.64%, respectively. This corresponds to a 5.5-fold increased risk of NMTC for first-degree family members [11].
Like many cancers, thyroid cancers may arise from mutations that may or may not be heritable. They can occur due to any mistake during DNA replication during cell division or may be induced due to the effect of carcinogens on DNA like ionizing radiation. Most thyroid cancers are the result of the accumulation of somatic mutations in the cancer genome, either driver mutations of oncogenesis or passenger mutations [12]. They are not present constitutionally in the individual but only in part of thyroid cells. By contrast, constitutional (germline) variants may predispose to cancer susceptibility and are present in affected individuals in all the body’s nucleic cells, as well as the cancer genome, and may therefore be heritable. Nonetheless, over 90% of all thyroid cancers are sporadic,
2. Familial non-medullary thyroid cancers
Familial NMTC (FNMTC) is clinically defined as the presence of the disease in two or more first-degree relatives of the patient. It encompasses a heterogeneous group of diseases, including diverse syndromic-associated tumors with a preponderance of non-thyroidal tumors and non-syndromic tumors with a preponderance of NMTC. Hereditary cancer syndromes associated with FNMTC account for 5% of all familial cases and include familial adenomatous polyposis (FAP) and its variant Gardner syndrome, Cowden syndrome, Carney complex, Werner syndrome, and DICER1 syndrome. Other syndromes with less established links to the development of NMTC include McCune-Albright syndrome, Ataxia-telangiectasia, Li-Fraumeni syndrome, and Peutz-Jeghers syndrome [16, 17, 18, 19, 20], also no epidemiological studies confirmed a significantly increased risk of NMTC for patients affected by these latter syndromes and for their relatives. Non-syndromic-associated conditions encompass pure familial PTC (fPTC) with or without oxyphilia, fPTC with papillary renal cell carcinoma, and fPTC with multinodular goiter [21, 22, 23]. The clinical characteristics of FNMTC are controversial. Several studies found an earlier age of onset, higher incidence of multifocality and lymph node metastasis, and a more aggressive outcome with more frequent relapses compared with sporadic disease [24], while other studies showed no significant increase in risk of recurrence or disease-related mortality in FNMTC cases compared to sporadic cases [25, 26]. Additionally, the second generation of parent-offspring FNMTC cases presents the disease at a younger age with more severe symptoms, indicating the presence of genetic anticipation [27].
Understanding the genetic basis of a heterogeneous disease such as FNMTC and the identification of biomarkers of disease aggressiveness can help to better stratify risk, allowing predictive screening of at-risk family members, improved surveillance guidance, and clinical management plan.
2.1 Genetic variants associated with risk of syndromic-associated disorders
Germline mutation (or “pathogenic variants”) accounting for syndromic FNMTC are highly penetrant and actionable, meaning that targeted genet testing is recommended when the clinician recognizes the clinical phenotype of the syndrome. Clinical characteristics and genes involved in the predisposition to syndromic FNMTC are summarized hereafter.
2.2 Genetic variants associated with risk of non-syndromic-associated disorders
Initial efforts to identify DTC susceptibility genes were conducted in the late 90s − early 2000s by conducting genome-wide linkage analysis in multigenerational families with multiple affected members, usually with attempt to replicate best hits in an independent set of smaller families. Some candidate genes within the mapped regions have been subsequently screened. To date, seven loci involved in FNMTC susceptibility have been mapped (1q21, 2q21, 8p23.1-p22, 8q24, 12p14, 14q32, 19p13.2), where the causal genes remain to be identified or confirmed in independent family sets. With the introduction of massive-parallel sequencing technologies in diagnostic and research laboratories in the 2010s, some of these regions have been more extremely screened highlighting new candidates (
Locus (name) | Family/cases in the discovery study | Candidate gene | Replication study | Ref. | Approach |
---|---|---|---|---|---|
14q32 ( | 1 French Canadian family, 18 MNG cases, 2 cases also with PTC | Investigation of 37 NMTC families (88 DTC cases) indicates that only a small proportion of FNMTC is attributable to | [23, 44] | Linkage analysis [44]; sequencing of | |
19p13.2 ( | 1 French family, 9 PTC cases (atypical carcinomas and adenomas with cell oxyphilia); 1 family, 3 cases | No linkage in subsequent study involving 56 NMTC families [45]; WES in the original family identified | [45, 46, 47] | Linkage analysis [45, 47]; WES + functional assay + sequencing of | |
1q21 ( | 1 family, 5 PTC cases, 2 family members with papillary renal neoplasia | No replication study. | [22] | Linkage analysis | |
2q21 ( | 1 Tasmanian family, 8 PTC cases | Linkage confirmed in an independent set of 80 families. | [48] | Linkage analysis | |
8p23.1-p22 ( | 1 Portuguese family, 5 PTC cases, 11 members with benign thyroid lesions | 17 candidates excluded. | No replication study. | [49] | Linkage analysis + gene expression profiling + sequencing + LOH analysis. |
8q24 | 1 family with PTC and melanoma | Linkage confirmed using 25 additional PTC families (86 cases and 13 obligate carriers) | [50] | Linkage analysis + target sequencing+ gene expression analysis. | |
12q14 | 38 families, 108 PTC cases | Association study on tag SNPs spanning the locus performed in 2 cases-control series from Ohio and Poland; biochemical assays showed that missense variants p.Gln149His and p.Arg617Cys impair the ability of SRGAP1 to inactivate CDC42. | [51] | Linkage analysis + association study + sequencing of | |
4q32 | 1 US family, 11 PTC cases and 2 cases with anaplastic thyroid carcinoma | (putative enhancer) | Not replicated in 38 NMTC families; rare single nucleotide variant absent in 2676 sporadic cases and 2470 controls from Poland and Ohio. | [52] | Linkage analysis + functional study |
1 family, 6 PTC cases | Target sequencing of | [53, 54] | WES + investigation of the expression of HABP2 in thyroid tissue samples + functional assay on variant p.Gly534Glu [53]; target sequencing of | ||
16p13.3 | 1 family, 6 PTC cases | Identified variant c.1037C>T (p.Ser346Phe) not found in 138 familial PTC cases; association study involving 1170 sporadic PTC cases and 1404 controls confirmed association with PTC. | [55] | WES, then genotyping of identified variant in 138 other PTC; association study in sporadic PTC and unrelated controls; RNA-Seq to assess effect on efficiency of RNA splicing machinery. | |
15q23 | 34 Chinese families, 77 cases | Variants c.961G>A and c.1100T>C identified in 2 families of the original Chinese study [56] not found in 33 Italian FNMTC families [57]. | [56, 57] | WES + functional study [56]; targeted sequencing to search for the 2 previously described variants only [57] | |
19q13.33 | 1 family with 5 NMTC cases | Variant rs78530808 (NOP53 p.Asp31His) identified in 3 of 44 additional families and absent in unaffected spouses; Functional studies showed oncogenic function but high frequency in the general population (MAF 1.8%) suggests a low-penetrant variant, possibly a modifier. | [58] | WES + targeted sequencing of candidate variants in familial cases and controls + functional assays. | |
22q12.1 | 5 NMTC families, 23 cases | Variants enrichment in MAPK/ERK and PI3K/AKT signaling pathways | No replication study. | [59] | WGS |
14q12 | 1 family with 10 members affected with PTC and/or melanoma + 23 NMTC families, 34 cases | No replication study. | [60, 61] | WGS in the key family with PTC and melanoma + target sequencing of | |
7q31.33 | 1 NMTC family, 8 cases | Another germline likely pathogenic variant in Lack of likely pathogenic | [59, 64] | WGS + functional analysis, including quantification of relative telomere length in variant carriers and noncarriers |
3. Genetic factors associated with sporadic DTC
3.1 Findings from association studies
Sporadic DTC, which represents the majority of all NMTC, is considered as a complex disease caused by multiple environmental and genetic risk factors. Common polymorphisms in low-penetrance genes or altering their expression are hypothesized to play an important role in sporadic DTC. The numerous association studies conducted in the general population used a case–control design that compares the distribution of the polymorphisms in a population of affected versus unaffected individuals. Two types of association studies can be performed: the candidate gene approach that focuses on a limited number of polymorphisms located in genes select based on
The candidate gene approach usually focuses on functional or tags polymorphisms located in specific genes selected based on their potential functional impact on the gene product or on its potential interaction with known environmental/lifestyle risk factors for the disease. Candidate genes selected for analysis in association studies on DTC are mostly involved in the following biological pathways: DNA repair, cell cycle regulation, and apoptosis, xenobiotic metabolism, thyroid function, MAPK pathway, immune response and inflammation, and obesity [65, 66, 67].
Figlioli et al. [67] proposed an exhaustive review of polymorphisms investigated in association studies published before September 2013. They reviewed 100 original articles and five meta-analyses and reported 91 significant SNPs (over 316 analyzed) from 127 genes. They also conducted a meta-analysis on 46 SNPs, which were reported by at least two studies, and reported 13 significant SNPs, of which six are located in the coding sequence of candidate genes:
With the completion of the human genome sequencing in 2003, GWAS involving hundreds of thousands to millions of SNPs across the human genome became more and more common. GWAS are usually composed with one discovery phase that aims to analyze a large number of variants, followed by a replication phase consisting of validation of the most significant variants in an independent sample. Since 2009, seven GWAS were published on NMTC risk, of which six were conducted in individuals of European ancestry [68, 69, 70, 71, 72, 73] and one was conducted in individuals of Asian ancestry from Korea [74] (Table 2). The latest GWAS conducted by Truong et al. [73] also included a small discovery sample of individuals of Oceanian ancestry but with no replication set. Among the GWAS conducted in the European ancestry population, one study focused on radiation-related PTC, with cases recruited in Belarus and aged 0−18 years old at the time of the Chernobyl accident [69]. All these GWAS were based on a relatively low number of cases compared to other cancers GWAS and the largest study included 3,001 cases is a meta-analysis of several studies with no replication phase [72].
Reference | Populations | Histology | Number of cases/controls in discovery phase | Number of cases/controls in replication phase |
---|---|---|---|---|
[68] | Iceland | NMTC | 192 /37,196 | |
Spain | NMTC | 89 /1,343 | ||
USA (Columbus) | PTC | 294/384 | ||
[69] | Belarus (Gomel) and Russia | PTC | 401/620 | 259/648 |
Age <18 years old at the time of Chernobyl accident | ||||
[70] | Italy | PTC | 701/499 | |
Italy | DTC | 1213/989 | ||
Italy | DTC | 326/730 | ||
Poland | DTC | 468/470 | ||
UK | DTC | 509/1,118 | ||
Spain | DTC | 443/420 | ||
[71] | Spain | DTC | 398/502 | |
Italy | DTC | 541/532 | ||
Spain (Galicia) | DTC | 240/531 | ||
Spain (Catalonia) | DTC | 354/408 | ||
[74] | Korea | DTC | 470/8,279 | 615/605 |
[72] | Iceland | NMTC | 1,003/278,991 | |
Netherlands | NMTC | 85/4,956 | ||
USA (Columbus) | PTC | 1,580/1,628 | ||
USA (Houston) | PTC | 250/363 | ||
Spain (Zaragosa) | NMTC | 83/1,612 | ||
[73] | European descents from France, French Polynesia, New Caledonia, Belarus, Cuba | DTC | 1,554/1,973 | |
Oceanian from France Polynesia, New Caledonia | DTC | 301/348 | ||
USA (Columbus) | PTC | 1,580/1,628 | ||
Italy | DTC | 649/431 |
Tables 3 and 4 summarize the significant and suggestive loci highlighted by these seven GWAS. All these variants are located in intronic or intergenic regions, except rs6793295 which is a missense variant in
Locus | Nearest gene(s) | Location | Lead SNPs | EA/OA | EAF | OR | p-value | Ref. | Ancestry | Remarks | Replicated |
---|---|---|---|---|---|---|---|---|---|---|---|
9q22.33 | FOXE1, | PTCSC2 intron | rs965513 | A/G | 0.34 | 1.75 | 1.7 x 10–27 | [68] | Eur | [75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89] | |
PTCSC2 | 0.34 | 1.65 | 4.8 x 10–12 | [69] | Eur | radiation-related DTC | |||||
_ | 1.78 | 2.7 x 10–10 | [70] | Eur | |||||||
0.33 | 1.65 | 2.7 x 10–23 | [71] | Eur | |||||||
PTCSC2 intron | rs1588635 | A/C | 0.40 | 1.69 | 2.0 x 10–58 | [72] | Eur | r2 = 0.99 with rs965513 in EUR | |||
0.39 | 1.64 | 2.1 x 10–21 | [73] | Eur | |||||||
Intergenic | rs72753537 | C/T | 0.07 | 1.41 | 7.7 x 10–6 | [74] | Asian | r2 = 0.001 with rs965513 in EAS | |||
14q13.3 | NKX2-1, | PTCSC3 promoter | rs944289 | T/C | 0.57 | 1.37 | 2.0 x 10–9 | [68] | Eur | [75, 79, 80, 82, 83, 84, 87, 88, 89, 92, 93, 94] | |
PTCSC3 | T/C | 0.56 | 1.24 | 1.5 x 10–5 | [71] | Eur | |||||
T/C | 0.46 | 1.25 | 1.4 x 10–6 | [74] | Asian | ||||||
Intergenic | rs116909374 | T/C | 0.03 | 1.81 | 1.1 x 10–16 | [72] | Eur | r2 = 0.006 with rs944289 in EUR | [79, 82, 88, 93] | ||
0.06 | 2.33 | 1.6 x 10–10 | [73] | Eur | |||||||
Intergenic | rs368187 | G/C | 0.58 | 1.39 | 5.1 x 10–23 | [72] | Eur | r2 = 0.70 with rs944289 in EUR | [93] | ||
0.63 | 1.47 | 3.8 x 10–13 | [73] | Eur | r2<0.01 with rs116909374 in EUR | ||||||
Intron | rs34081947 | T/C | 0.41 | 1.27 | 1.2 x 10–7 | [74] | Asian | r2 = 0.61 with rs944289 in EAS r2 = 0.90 with rs368187 in EAS | |||
2q35 | DIRC3 | Intron | rs6759952 | T/C | 0.38 | 1.25 | 6.4 x 10–10 | [70] | Eur | [108] | |
Noncoding transcript | rs11693806 | C/G | 0.32 | 1.43 | 1.5 x 10–24 | [72] | Eur | r2 = 0.45 with rs6759952 in EUR | |||
Noncoding transcript | rs3821098 | T/C | 0.28 | 1.39 | 1.1 x 10–10 | [73] | Eur | r2 = 0.98 with rs11693806 in EUR | |||
Intron | rs12990503 | G/C | 0.63 | 1.34 | 3.5 x 10–9 | [74] | Asian | r2 = 0.97 with rs11693806 in EAS | |||
8p12 | NRG1 | Intron | rs2466076 | G/T | 0.48 | 1.32 | 1.5 x 10–17 | [72] | Eur | r2 = 0.94 with rs2439302 in EUR | |
Intron | rs142450470 | T/- | 0.35 | 1.33 | 1.5 x 10–7 | [73] | Eur | r2 = 0.45 with rs2439302 in EUR | |||
Significant for PTC (p = 9.4 x 10−9) | |||||||||||
Intron | rs2439302 | G/C | 0.21 | 1.37 | 1.4 x 10–9 | [74] | Asian | r2 = 0.95 with rs2466076 in EAS r2 = 0.39 with rs142450470 in EAS | [73, 84, 88, 89, 92, 94, 98] | ||
Intron | rs6996585 | G/A | 0.23 | 1.39 | 1.1 x 10–10 | [74] | Asian | r2 = 0.67 with rs2439302 in EAS | |||
Intron | rs12542743 | C/T | 0.25 | 1.36 | 4.6 x 10–10 | [74] | Asian | r2 = 0.37 with rs2439302 in EAS |
Locus | Nearest gene(s) | Location | Lead SNPs | EA/OA | EAF | OR | p-value | Ref. | Ancestry | Remarks | Replicated |
---|---|---|---|---|---|---|---|---|---|---|---|
1q42.2 | Intron | rs12129938 | A/G | 0.79 | 1.32 | 4.0 × 10–11 | [72] | Eur | [109] | ||
Intron | rs4649295 | T/C | 0.82 | 1.43 | 6.0 × 10–8 | [74] | Asian | r2 = 0.67 with rs12129938 in EAS | [109] | ||
10q.24.33 | Intergenic | rs7902587 | T/C | 0.11 | 1.41 | 5.4 × 10–11 | [72] | Eur | |||
5q22.1 | Intron | rs73227498 | A/T | 0.87 | 1.37 | 3.0 × 10–10 | [72] | Eur | |||
15q22.33 | Intron | rs2289261 | C/G | 0.68 | 1.23 | 3.1 × 10–9 | [72] | Eur | |||
Intron | rs56062135 | T/C | 0.25 | 1.24 | 4.9 × 10–9 | [72] | Eur | ||||
16q23.2 | Intergenic | rs16950982 | G/A | 0.37 | 1.22 | 4.7 × 10–9 | [73] | Eur | |||
12q14.3 | Intron | rs11175834 | T/C | 0.15 | 1.37 | 4.3 × 10–8 | [74] | Asian | |||
1p13.3 | Intron | rs4915076 | T/C | 0.70 | 1.33 | 8.5 × 10–8 | [74] | Asian | [110] | ||
3q26.2 | Exon, p.Ser249Gly | rs6793295 | T/C | 0.76 | 1.23 | 2.7 × 10–8 | [72] | Eur | |||
5p15.33 | Intron | rs10069690 | T/C | 0.27 | 1.20 | 3.2 × 10–7 | [72] | Eur | [111, 112] | ||
Intron | rs7726159 | A/C | 0.39 | 1.17 | 5 × 10–6 | [73] | Eur | r2 = 0.44 with rs10069690 in EUR | |||
19p12 | Intron | rs7260863 | T/C | 0.20 | 1.22 | 8.7 × 10–7 | [73] | Eur | |||
4q21.1 | Intergenic | rs1874564 | G/A | 0.69 | 1.31 | 2.0 × 10–7 | [74] | Asian | |||
3p14.2 | Intergenic | rs9858271 | G/A | 0.43 | 1.26 | 6.8 × 10–7 | [74] | Asian | |||
1p31.3 | Intron | rs334729 | C/T | 0.05 | 1.43 | 7.6 × 10–7 | [73] | Eur | |||
7q31.1 | Intron | rs7800391 | T/C | 0.34 | 1.25 | 5.7 × 10–6 | [70] | Eur | Specific to Italian pop | ||
Intron | rs10238549 | C/T | 0.63 | 1.27 | 4.1 × 10–6 | [70] | Eur | Specific to Italian pop | |||
3q25.32 | Intron | rs7617304 | A/G | 0.25 | 1.25 | 4.6 × 10–5 | [70] | Eur | Specific to Italian pop | ||
10q26.12 | Intron | rs1254167 | C/G | 0.08 | 1.38 | 5.9 × 10–5 | [71] | Eur | |||
Intron | rs10788123 | T/C | 0.20 | 1.69 | 3.2 × 10–5 | [71] | Eur | ||||
9q34.3 | Intergenic | rs10781500 | C/T | 0.60 | 1.23 | 3.5 × 10–5 | [70] | Eur | Specific to Italian pop | ||
19p13.2 | Intron | rs7248104 | A/G | 0.36 | 1.22 | 2.0 × 10–5 | [74] | Asian | |||
6q14.1 | Intergenic | rs4075570 | G/A | 0.35 | 0.82 | 2.0 × 10–4 | [71] | Eur | [113] |
3.1.1 Locus 9p22.33
At 9p22.33, the most robust association reported in GWAS was for rs96551 or rs1588635, a highly correlated proxy (r2 = 0.99 in the European population from the 1000 Genomes Project). Rs965513 has been associated with radiation-related DTC as well as with sporadic DTC, and it was subsequently consistently replicated in several different populations of European ancestry [75, 76, 77, 78, 79, 80, 81, 82], of Asian ancestry [81, 83, 84, 85] as well as in admixed populations from Oceania [79, 86], Cuba [87], Colombia [88], or Kazakhstan [89].
In 2009, a study focusing on the role of
3.1.2 Locus 14q13.3
Rs944289 was the first variant identified by GWAS at 14q13.3 [68]. This association was replicated in subsequent studies conducted in diverse populations [75, 79, 80, 82, 83, 84, 87, 88, 89, 92, 93, 94]. However, the most recent GWAS, conducted in several European populations [72, 73] and in one Asian population from Korea [74], analyzed a higher number of SNPs (using genotyped and imputed SNPs) and reported the strongest signal for respectively rs368187 and rs34081947 (r2 = 0.98 in Europeans and r2 = 0.90 in East Asians from the 1000 Genomes Project), which are in moderate LD with rs944289 (r2<0.70 in Europeans or East Asians from the 1000 Genomes Project). The recent European GWAS also reported an independent signal at rs116909374 and rs368187. Rs116909374 replicated only in studies on European ancestry populations [79, 82, 88, 93] as this SNP is monomorphic or very rare in Asian populations. Interestingly, rs368187 is in high LD with rs34081947, which was highlighted by the GWAS conducted in the Korean population (r2 = 0.98 in East Asians from the 1000 Genome Project) (Table 3).
In 2012, Jendrzejewski et al. [95] described a novel lincRNA named
3.1.3 Locus 2q35
The association between DTC and variants at 2q35 was first highlighted in 2009 in a GWAS that investigated genetic factors associated with thyroid stimulating hormone (TSH) levels in blood in 27,758 Icelandic individuals. The authors further investigated the role of the top SNPs associated with circulating TSH levels in DTC susceptibility of which rs966423 at 2q35 [82]. Another GWAS on DTC risk, conducted in 2013 in an Italian population [70], reported rs6759952 as the lead SNP at 2q35, which is moderately correlated to rs966423 (r2 = 0.69 in Europeans from the 1000 Genomes project). The most recent GWAS from 2017 [72, 73, 74] reported three new SNPs (rs11693806, rs3821098, and rs12990503) with the strongest signal at 2q35, which are highly correlated to each other but only moderately with the two previous reported SNPs (Table 2).
Finally, the recent
3.1.4 Locus 8p12
At 8p12, rs2439302 was first highlighted in the Icelandic GWAS on TSH levels, and it was found to be associated with DTC risk in subsequent analyses [82]. The association between DTC and rs2439302 was then replicated in different populations [84, 88, 89, 92, 94, 98]. The recent GWAS [72, 73, 74] replicated this SNP and also reported several other variants within the gene
In 2018, He et al. [99] reported that the risk allele [G] of rs2439302 was associated with the expression of multiple NGR1 isoforms in normal thyroid tissue. They also suggested that multiple enhancer variants exist at this locus that may have a combinatory effect on the expression of
3.1.5 Other loci
Among the other loci reported by GWAS, only SNPs at four loci (1q42.2, 5p15.33, 1p13.3, 6q14.1) were replicated in independent studies on DTC (Table 4). Interestingly, some of the variants reported in Table 4 were also previously associated with other diseases or traits. For instance, the SNPs at 5p15.33 (
3.2 Polygenic risk scores
Based on findings from GWAS, polygenic risk scores (PRS), which are calculated by computing the sum of risk alleles of identified susceptibility SNPs weighted by the effect size estimate from the GWAS, were proposed to predict DTC risk. Several studies evaluated a DTC PRS in different populations [114, 115, 116, 117, 118]. The most recent studies used PRS including 10 to 12 SNPs reported by the meta-analysis of GWAS [72] and estimated odds ratios per standard deviation of PRS from 1.55 to 2.31. Liyanarachchi et al. [118] estimated that about 8% of the genetic predisposition to PTC could be accounted for by 10 SNPs (rs12129938, rs11693806, rs6793295, rs73227498, rs2466076, rs1588635, rs7902587, rs368187, rs116909374, rs2289261). They also estimated that individuals of European ancestry in the highest decile of PRS had a 6.9 higher risk to develop PTC than individuals in the lowest decile group. A recent study reported that the PRS improved significantly predictive scores based on clinical factors in the prediction of subsequent thyroid cancers in childhood cancer survivors of European ancestry [119]. Future studies should investigate the combined effect of PRS and exposure to lifestyle and environmental factors in order to enhance individualized DTC risk prediction. There is also a need to extend the DTC PRS to other ethnic groups.
4. Conclusion
Despite the solid evidence for heritability of thyroid cancer, only a handful of variants have been significantly associated with an increased risk of DTC representing the most common form of thyroid cancer. The high heritability of the disease is likely due to the contribution of rare high-penetrance variants in some cases and the combination of common low-penetrance variants in others, as well as the influence of common shared environmental factors in DTC-prone families or in specific groups of the general population. So far, efforts to identify DTC predisposing genes outside of syndromic FNMTC led to the identification of mainly low-to-moderate penetrance genes, and routine genetic testing for these genes is not recommended. Further large studies to characterize their penetrance and function and to identify new DTC-associated loci or alternative hereditary mechanisms such as epigenetic modifications are required to improve our understanding of DTC tumorigenesis. Ultimately, risk prediction models integrating family history of DTC, PRS, and some modifiable risk factors (obesity, exposure to ionizing radiations from medical diagnostic procedure, etc.) may help stratify individuals according to their risk of developing DTC, which can be useful for elaborating screening policies. Moreover, inherited genetic factors can also impact the final outcome of the disease such as histological subtypes, localization of metastases, or molecular profiling of the tumor, and their characterization can help to predict effectiveness of the initial treatment [120].
References
- 1.
Chrisoulidou A, Boudina M, Tzemailas A, Doumala E, Iliadou PK, Patakiouta F, et al. Histological subtype is the most important determinant of survival in metastatic papillary thyroid cancer. Thyroid Research. 2011; 4 (1):12. DOI: 10.1186/1756-6614-4-12 - 2.
Jemal A, Siegel R, Xu J, Ward E. Cancer statistics, 2010. CA: A Cancer Journal for Clinicians. 2010; 60 (5):277-300. DOI: 10.3322/caac.20073 - 3.
Bray F, Ferlay J, Soerjomataram I, Siegel R, Torre L, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: A Cancer Journal for Clinicians. 2018; 68 (6):394-424 - 4.
Truong T, Rougier Y, Dubourdieu D, Guihenneuc-Jouyaux C, Orsi L, Hémon D, et al. Time trends and geographic variations for thyroid cancer in New Caledonia, a very high incidence area (1985-1999). European Journal of Cancer Prevention. 2007; 16 (1995):62-70 - 5.
Ito Y, Nikiforov YE, Schlumberger M, Vigneri R. Increasing incidence of thyroid cancer: Controversies explored. Nature Reviews. Endocrinology. 2013; 9 (3):178-184. DOI: 10.1038/nrendo.2012.257 - 6.
Goldgar DE, Easton DF, Cannon-Albright LA, Skolnick MH. Systematic population-based assessment of cancer risk in first-degree relatives of cancer probands. Journal of the National Cancer Institute. 1994; 86 (21):1600-1608 - 7.
Dong C, Hemminki K. Modification of cancer risks in offspring by sibling and parental cancers from 2,112,616 nuclear families. International Journal of Cancer. 2001; 92 (1):144-150 - 8.
Hemminki K, Dong C. Familial relationships in thyroid cancer by histo-pathological type. International Journal of Cancer. 2000; 85 (2):201-205 - 9.
Hrafnkelsson J, Tulinius H, Jonasson JG, Sigvaldason H. Familial non-medullary thyroid cancer in Iceland. Journal of Medical Genetics. 2001; 38 (3):189-191. DOI: 10.1136/jmg.38.3.189 - 10.
Frich L, Glattre E, Akslen LA. Familial occurrence of nonmedullary thyroid cancer: A population-based study of 5673 first-degree relatives of thyroid cancer patients from Norway. Cancer Epidemiology, Biomarkers & Prevention. 2001; 10 (2):113-117 - 11.
Lin HT, Liu FC, Lin SF, Kuo CF, Chen YY, Yu HP. Familial aggregation and heritability of nonmedullary thyroid cancer in an asian population: A nationwide cohort study. The Journal of Clinical Endocrinology and Metabolism. 2020; 105 (7). DOI: 10.1210/clinem/dgaa191 - 12.
Xing M. Molecular pathogenesis and mechanisms of thyroid cancer. Nature Reviews. Cancer. 2013; 13 (3):184-199. DOI: 10.1038/nrc3431 - 13.
Romei C, Ciampi R, Elisei R. A comprehensive overview of the role of the RET proto-oncogene in thyroid carcinoma. Nature Reviews. Endocrinology. 2016; 12 (4):192-202. DOI: 10.1038/nrendo.2016.11 - 14.
Hedayati M, Zarif Yeganeh M, Sheikholeslami S, Afsari F. Diversity of mutations in the RET proto-oncogene and its oncogenic mechanism in medullary thyroid cancer. Critical Reviews in Clinical Laboratory Sciences. 2016; 53 (4):217-227. DOI: 10.3109/10408363.2015.1129529 - 15.
Accardo G, Conzo G, Esposito D, Gambardella C, Mazzella M, Castaldo F, et al. Genetics of medullary thyroid cancer: An overview. International Journal of Surgery. 2017; 41 (Suppl. 1):S2-S6. DOI: 10.1016/j.ijsu.2017.02.064 - 16.
Hincza K, Kowalik A, Kowalska A. Current knowledge of germline genetic risk factors for the development of non-medullary thyroid cancer. Genes (Basel). 2019; 10 (7):482. DOI: 10.3390/genes10070482 - 17.
Schultz KAP, Rednam SP, Kamihara J, Doros L, Achatz MI, Wasserman JD, et al. PTEN, DICER1, FH, and their associated tumor susceptibility syndromes: Clinical features, genetics, and surveillance recommendations in childhood. Clinical Cancer Research. 2017; 23 (12):e76-e82. DOI: 10.1158/1078-0432.CCR-17-0629 - 18.
Son EJ, Nose V. Familial follicular cell-derived thyroid carcinoma. Frontiers in Endocrinology. 2012; 3 :61. DOI: 10.3389/fendo.2012.00061 - 19.
Legrand MA, Raverot G, Nicolino M, Chapurlat R. GNAS mutated thyroid carcinoma in a patient with Mc Cune Albright syndrome. Bone Reports. 2020; 13 :100299. DOI: 10.1016/j.bonr.2020.100299 - 20.
Lacka K, Maciejewski A, Stawny B, Lacki JK. Follicular thyroid cancer in a patient with pendred syndrome. Annales d’Endocrinologie. 2021; 82 (6):622-624. DOI: 10.1016/j.ando.2021.05.001 - 21.
Malchoff CD, Malchoff DM. Familial nonmedullary thyroid carcinoma. Cancer Control. 2006; 13 (2):106-110 - 22.
Malchoff CD, Sarfarazi M, Tendler B, Forouhar F, Whalen G, Joshi V, et al. Papillary thyroid carcinoma associated with papillary renal neoplasia: Genetic linkage analysis of a distinct heritable tumor syndrome. The Journal of Clinical Endocrinology and Metabolism. 2000; 85 (5):1758-1764 - 23.
Rio Frio T, Bahubeshi A, Kanellopoulou C, Hamel N, Niedziela M, Sabbaghian N, et al. DICER1 mutations in familial multinodular goiter with and without ovarian Sertoli-Leydig cell tumors. JAMA. 2011; 305 (1):68-77. DOI: 10.1001/jama.2010.1910 - 24.
Moses W, Weng J, Kebebew E. Prevalence, clinicopathologic features, and somatic genetic mutation profile in familial versus sporadic nonmedullary thyroid cancer. Thyroid. 2011; 21 (4):367-371. DOI: 10.1089/thy.2010.0256 - 25.
Pitoia F, Cross G, Salvai ME, Abelleira E, Niepomniszcze H. Patients with familial non-medullary thyroid cancer have an outcome similar to that of patients with sporadic papillary thyroid tumors. Arquivos Brasileiros de Endocrinologia & Metabologia. 2011; 55 (3):219-223. DOI: 10.1590/s0004-27302011000300007 - 26.
Capezzone M, Fralassi N, Secchi C, Cantara S, Brilli L, Pilli T, et al. Long-term clinical outcome in familial and sporadic papillary thyroid carcinoma. The European Thyroid Journal. 2020; 9 (4):213-220. DOI: 10.1159/000506955 - 27.
Park YJ, Ahn HY, Choi HS, Kim KW, Park DJ, Cho BY. The long-term outcomes of the second generation of familial nonmedullary thyroid carcinoma are more aggressive than sporadic cases. Thyroid. 2012; 22 (4):356-362. DOI: 10.1089/thy.2011.0163 - 28.
Dinarvand P, Davaro EP, Doan JV, Ising ME, Evans NR, Phillips NJ, et al. Familial adenomatous polyposis syndrome: An update and review of extraintestinal manifestations. Archives of Pathology & Laboratory Medicine. 2019; 143 (11):1382-1398. DOI: 10.5858/arpa.2018-0570-RA - 29.
Harach HR, Williams GT, Williams ED. Familial adenomatous polyposis associated thyroid carcinoma: A distinct type of follicular cell neoplasm. Histopathology. 1994; 25 (6):549-561. DOI: 10.1111/j.1365-2559.1994.tb01374.x - 30.
Soravia C, Sugg SL, Berk T, Mitri A, Cheng H, Gallinger S, et al. Familial adenomatous polyposis-associated thyroid cancer: A clinical, pathological, and molecular genetics study. The American Journal of Pathology. 1999; 154 (1):127-135. DOI: 10.1016/S0002-9440(10)65259-5 - 31.
Cameselle-Teijeiro J, Chan JK. Cribriform-morular variant of papillary carcinoma: A distinctive variant representing the sporadic counterpart of familial adenomatous polyposis-associated thyroid carcinoma? Modern Pathology. 1999; 12 (4):400-411 - 32.
Herraiz M, Barbesino G, Faquin W, Chan-Smutko G, Patel D, Shannon KM, et al. Prevalence of thyroid cancer in familial adenomatous polyposis syndrome and the role of screening ultrasound examinations. Clinical Gastroenterology and Hepatology. 2007; 5 (3):367-373. DOI: 10.1016/j.cgh.2006.10.019 - 33.
Crippa S, Saletti P, Barizzi J, Mazzucchelli L. The clinical management in familial adenomatous polyposis deserves continuous monitoring for thyroid carcinoma. BML Case Reports. 2012; 2012 :1-3. DOI: 10.1136/bcr-2012-007046 - 34.
Tan MH, Mester JL, Ngeow J, Rybicki LA, Orloff MS, Eng C. Lifetime cancer risks in individuals with germline PTEN mutations. Clinical Cancer Research. 2012; 18 (2):400-407. DOI: 10.1158/1078-0432.CCR-11-2283 - 35.
Ngeow J, Mester J, Rybicki LA, Ni Y, Milas M, Eng C. Incidence and clinical characteristics of thyroid cancer in prospective series of individuals with Cowden and Cowden-like syndrome characterized by germline PTEN, SDH, or KLLN alterations. The Journal of Clinical Endocrinology and Metabolism. 2011; 96 (12):E2063-E2071. DOI: 10.1210/jc.2011-1616 - 36.
Bertherat J, Horvath A, Groussin L, Grabar S, Boikos S, Cazabat L, et al. Mutations in regulatory subunit type 1A of cyclic adenosine 5’-monophosphate-dependent protein kinase (PRKAR1A): Phenotype analysis in 353 patients and 80 different genotypes. The Journal of Clinical Endocrinology and Metabolism. 2009; 94 (6):2085-2091. DOI: 10.1210/jc.2008-2333 - 37.
Bossis I, Stratakis CA. Minireview: PRKAR1A: Normal and abnormal functions. Endocrinology. 2004; 145 (12):5452-5458. DOI: 10.1210/en.2004-0900 - 38.
Nikiforov YE. RET/PTC rearrangement in thyroid tumors. Endocrine Pathology. 2002; 13 (1):3-16. DOI: 10.1385/ep:13:1:03 - 39.
Yu CE, Oshima J, Hisama FM, Matthews S, Trask BJ, Schellenberg GD. A YAC, P1, and cosmid contig and 17 new polymorphic markers for the Werner syndrome region at 8p12-p21. Genomics. 1996; 35 (3):431-440. DOI: 10.1006/geno.1996.0382 - 40.
Lauper JM, Krause A, Vaughan TL, Monnat RJ Jr. Spectrum and risk of neoplasia in Werner syndrome: A systematic review. PLoS One. 2013; 8 (4):e59709. DOI: 10.1371/journal.pone.0059709 - 41.
Slade I, Bacchelli C, Davies H, Murray A, Abbaszadeh F, Hanks S, et al. DICER1 syndrome: Clarifying the diagnosis, clinical features and management implications of a pleiotropic tumour predisposition syndrome. Journal of Medical Genetics. 2011; 48 (4):273-278. DOI: 10.1136/jmg.2010.083790 - 42.
Fuziwara CS, Kimura ET. MicroRNAs in thyroid development, function and tumorigenesis. Molecular and Cellular Endocrinology. 2017; 456 :44-50. DOI: 10.1016/j.mce.2016.12.017 - 43.
Khan NE, Bauer AJ, KAP S, Doros L, Decastro RM, Ling A, et al. Quantification of thyroid cancer and multinodular goiter risk in the DICER1 syndrome: A family-based cohort study. The Journal of Clinical Endocrinology and Metabolism. 2017; 102 (5):1614-1622. DOI: 10.1210/jc.2016-2954 - 44.
Bignell GR, Canzian F, Shayeghi M, Stark M, Shugart YY, Biggs P, et al. Familial nontoxic multinodular thyroid goiter locus maps to chromosome 14q but does not account for familial nonmedullary thyroid cancer. American Journal of Human Genetics. 1997; 61 (5):1123-1130 - 45.
Lesueur F, Stark M, Tocco T, Ayadi H, Delisle MJ, Goldgar DE, et al. Genetic heterogeneity in familial nonmedullary thyroid carcinoma: Exclusion of linkage to RET, MNG1, and TCO in 56 families. NMTC consortium. The Journal of Clinical Endocrinology and Metabolism. 1999; 84 (6):2157-2162 - 46.
Diquigiovanni C, Bergamini C, Evangelisti C, Isidori F, Vettori A, Tiso N, et al. Mutant MYO1F alters the mitochondrial network and induces tumor proliferation in thyroid cancer. International Journal of Cancer. 2018; 143 (7):1706-1719. DOI: 10.1002/ijc.31548 - 47.
Canzian F, Amati P, Harach HR, Kraimps JL, Lesueur F, Barbier J, et al. A gene predisposing to familial thyroid tumors with cell oxyphilia maps to chromosome 19p13.2. American Journal of Human Genetics. 1998; 63 (6):1743-1748 - 48.
McKay JD, Lesueur F, Jonard L, Pastore A, Williamson J, Hoffman L, et al. Localization of a susceptibility gene for familial nonmedullary thyroid carcinoma to chromosome 2q21. American Journal of Human Genetics. 2001; 69 (2):440-446 - 49.
Cavaco BM, Batista PF, Sobrinho LG, Leite V. Mapping a new familial thyroid epithelial neoplasia susceptibility locus to chromosome 8p23.1-p22 by high-density single-nucleotide polymorphism genome-wide linkage analysis. The Journal of Clinical Endocrinology and Metabolism. 2008; 93 (11):4426-4430 - 50.
He H, Nagy R, Liyanarachchi S, Jiao H, Li W, Suster S, et al. A susceptibility locus for papillary thyroid carcinoma on chromosome 8q24. Cancer Research. 2009; 69 (2):625-631 - 51.
He H, Bronisz A, Liyanarachchi S, Nagy R, Li W, Huang Y, et al. SRGAP1 is a candidate gene for papillary thyroid carcinoma susceptibility. The Journal of Clinical Endocrinology and Metabolism. 2013; 98 (5):E973-E980. DOI: 10.1210/jc.2012-3823 - 52.
He H, Li W, Wu D, Nagy R, Liyanarachchi S, Akagi K, et al. Ultra-rare mutation in long-range enhancer predisposes to thyroid carcinoma with high penetrance. PLoS One. 2013; 8 (5):e61920. DOI: 10.1371/journal.pone.0061920 - 53.
Gara SK, Jia L, Merino MJ, Agarwal SK, Zhang L, Cam M, et al. Germline HABP2 mutation causing familial nonmedullary thyroid cancer. The New England Journal of Medicine. 2015; 373 (5):448-455. DOI: 10.1056/NEJMoa1502449 - 54.
Tomsic J, He H, de la Chapelle A. HABP2 mutation and nonmedullary thyroid cancer. The New England Journal of Medicine. 2015; 373 (21):2086. DOI: 10.1056/NEJMc1511631 - 55.
Tomsic J, He H, Akagi K, Liyanarachchi S, Pan Q , Bertani B, et al. A germline mutation in SRRM2, a splicing factor gene, is implicated in papillary thyroid carcinoma predisposition. Scientific Reports. 2015; 5 :10566. DOI: 10.1038/srep10566 - 56.
Ye F, Gao H, Xiao L, Zuo Z, Liu Y, Zhao Q , et al. Whole exome and target sequencing identifies MAP2K5 as novel susceptibility gene for familial non-medullary thyroid carcinoma. International Journal of Cancer. 2019; 144 (6):1321-1330. DOI: 10.1002/ijc.31825 - 57.
Cirello V, Colombo C, Persani L, Fugazzola L. Absence of the MAP2K5 germline variants c.G961A and c.T1100C in a wide series of familial nonmedullary thyroid carcinoma Italian families. International Journal of Cancer. 2019; 145 (2):600. DOI: 10.1002/ijc.32244 - 58.
Orois A, Gara SK, Mora M, Halperin I, Martinez S, Alfayate R, et al. NOP53 as a candidate modifier locus for familial non-medullary thyroid cancer. Genes (Basel). 2019; 10 (11). DOI: 10.3390/genes10110899 - 59.
Srivastava A, Kumar A, Giangiobbe S, Bonora E, Hemminki K, Forsti A, et al. Whole genome sequencing of familial non-medullary thyroid cancer identifies germline alterations in MAPK/ERK and PI3K/AKT signaling pathways. Biomolecules. 2019; 9 (10). DOI: 10.3390/biom9100605 - 60.
He H, Li W, Comiskey DF, Liyanarachchi S, Nieminen TT, Wang Y, et al. A truncating germline mutation of TINF2 in individuals with thyroid cancer or melanoma results in longer telomeres. Thyroid. 2020; 30 (2):204-213. DOI: 10.1089/thy.2019.0156 - 61.
Wang Y, Liyanarachchi S, Miller KE, Nieminen TT, Comiskey DF Jr, Li W, et al. Identification of rare variants predisposing to thyroid cancer. Thyroid. 2019; 29 (7):946-955. DOI: 10.1089/thy.2018.0736 - 62.
Wilson TL, Hattangady N, Lerario AM, Williams C, Koeppe E, Quinonez S, et al. A new POT1 germline mutation-expanding the spectrum of POT1-associated cancers. Familial Cancer. 2017; 16 (4):561-566. DOI: 10.1007/s10689-017-9984-y - 63.
Orois A, Badenas C, Reverter JL, Lopez V, Potrony M, Mora M, et al. Lack of mutations in POT1 gene in selected families with familial non-medullary thyroid cancer. Hormone Therapy for Cancer. 2020; 11 (2):111-116. DOI: 10.1007/s12672-020-00383-5 - 64.
Srivastava A, Miao B, Skopelitou D, Kumar V, Kumar A, Paramasivam N, et al. A germline mutation in the POT1 gene is a candidate for familial non-medullary thyroid cancer. Cancers (Basel). 2020; 12 (6):1441. DOI: 10.3390/cancers12061441 - 65.
Adjadj E, Schlumberger M, Vathaire F De. Germ-line DNA polymorphisms and susceptibility to differentiated thyroid cancer. The Lancet Oncology 2009;10(2):181-190. DOI: 10.1016/S1470-2045(09)70020-8 - 66.
Landa I, Robledo M. Association studies in thyroid cancer susceptibility: Are we on the right track? Journal of Molecular Endocrinology. 2011; 47 :R43-R58. DOI: 10.1530/JME-11-0005 - 67.
Figlioli G, Elisei R, Romei C, Melaiu O, Cipollini M, Bambi F, et al. A comprehensive meta-analysis of case – control association studies to evaluate polymorphisms associated with the risk of differentiated thyroid carcinoma. Cancer Causes & Control. 2016; 25 (4):700-713. DOI: 10.1158/1055-9965.EPI-15-0652 - 68.
Gudmundsson J, Sulem P, Gudbjartsson DF, Jonasson JG, Sigurdsson A, Bergthorsson JT, et al. Common variants on 9q22.33 and 14q13.3 predispose to thyroid cancer in European populations. Nature Genetics. 2009; 41 (4):460-464. DOI: 10.1038/ng.339 - 69.
Takahashi M, Saenko VA, Rogounovitch TI, Kawaguchi T, Drozd VM, Takigawa-Imamura H, et al. The FOXE1 locus is a major genetic determinant for radiation-related thyroid carcinoma in Chernobyl. Human Molecular Genetics. 2010; 19 (12):2516-2523. DOI: 10.1093/hmg/ddq123 - 70.
Köhler A, Chen B, Gemignani F, Elisei R, Romei C, Figlioli G, et al. Genome-wide association study on differentiated thyroid cancer. The Journal of Clinical Endocrinology and Metabolism. 2013; 98 (10):E1674-E1681. DOI: 10.1210/jc.2013-1941 - 71.
Mancikova V, Cruz R, Inglada-Pérez L, Fernández-Rozadilla C, Landa I, Cameselle-Teijeiro J, et al. Thyroid cancer GWAS identifies 10q26.12 and 6q14.1 as novel susceptibility loci and reveals genetic heterogeneity among populations. International Journal of Cancer. 2015; 137 (8):1870-1878. DOI: 10.1002/ijc.29557 - 72.
Gudmundsson J, Thorleifsson G, Sigurdsson JK, Stefansdottir L, Jonasson JG, Gudjonsson SA, et al. A genome-wide association study yields five novel thyroid cancer risk loci. Nature Communications. 2017; 8 :14517. DOI: 10.1038/ncomms14517 - 73.
Truong T, Lesueur F, Sugier P, Guibon J, Xhaard C, Karimi M, et al. Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium. International Journal of Cancer. 2021; 148 (12):2935-2946. DOI: 10.1002/ijc.33488 - 74.
Son HY, Hwangbo Y, Yoo SK, Im SW, Yang SD, Kwak SJ, et al. Genome-wide association and expression quantitative trait loci studies identify multiple susceptibility loci for thyroid cancer. Nature Communications. 2017; 8 :15966. DOI: 10.1038/ncomms15966 - 75.
Jones AM, Howarth KM, Martin L, Gorman M, Mihai R, Moss L, 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. Journal of Medical Genetics. 2012; 49, 63 (3):158. DOI: 10.1136/jmedgenet-2011-100586 - 76.
Tomaz RA, Sousa I, Silva JG, Santos C, Teixeira MR, Leite V, et al. FOXE1 polymorphisms are associated with familial and sporadic nonmedullary thyroid cancer susceptibility. Clinical Endocrinology (Oxf). 2012;77(6):926-933. DOI: 10.1111/j.1365-2265.2012.04505.x - 77.
Penna-Martinez M, Epp F, Kahles H, Ramos-Lopez E, Hinsch N, Hansmann ML, et al. FOXE1 association with differentiated thyroid cancer and its progression. Thyroid. 2014; 24 (5):845-851 - 78.
Jendrzejewski J, Liyanarachchi S, Nagy R, Senter L, Wakely PE, Thomas A, et al. Papillary thyroid carcinoma: Association between germline DNA variant markers and clinical parameters. Thyroid. 2016; 26 (9):1276-1284 - 79.
Tcheandjieu C, Lesueur F, Sanchez M, Baron-Dubourdieu D, Guizard AV, Mulot C, et al. Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 9q22.33 and 14q13.3 detects novel candidate functional SNPs in Europeans from metropolitan France and Melanesians from New Caledonia. International Journal of Cancer. 2016; 139 (3):617-627 - 80.
Kula D, Kalemba M, Puch Z, Polaska J, Swierniak M, Rusinek D, et al. Age at diagnosis and gender modify the risk of 9q22 and 14q13 polymorphisms for papillary thyroid carcinoma. Endokrynologia Polska. 2017; 68 (3):283-289 - 81.
Nikitski AV, Rogounovitch TI, Bychkov A, Takahashi M, Yoshiura KI, Mitsutake N, et al. Genotype analyses in the Japanese and belarusian populations reveal independent effects of rs965513 and rs1867277 but do not support the role of FOXE1 polyalanine tract length in conferring risk for papillary thyroid carcinoma. Thyroid. 2017; 27 (2):224-235 - 82.
Gudmundsson J, Sulem P, Gudbjartsson DF, Jonasson JG, Masson G, He H, et al. Discovery of common variants associated with low TSH levels and thyroid cancer risk. Nature Genetics. 2012; 44 (3):319-322 - 83.
Matsuse M, Takahashi M, Mitsutake N, Nishihara E, Hirokawa M, Kawaguchi T, et al. The FOXE1 and NKX2-1 loci are associated with susceptibility to papillary thyroid carcinoma in the Japanese population. Journal of Medical Genetics. 2011; 48 (9):645-648. DOI: 10.1136/jmedgenet-2011-100063 - 84.
Wang YL, Feng SH, Guo SC, Wei WJ, Li DS, Wang Y, 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. Journal of Medical Genetics. 2013; 50 (10):689-695. DOI: 10.1136/jmedgenet-2013-101687 - 85.
Guo S, Wang YL, Li Y, Jin L, Xiong M, Ji QH, et al. Significant SNPs have limited prediction ability for thyroid cancer. Cancer Medicine. 2014; 3 (3):731-735. DOI: 10.1002/cam4.211 - 86.
Maillard S, Damiola F, Clero E, Pertesi M, Robinot N, Rachédi F, 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 (4):e0123700. DOI: 10.1371/JOURNAL.PONE.0123700 - 87.
Pereda CM, Lesueur F, Pertesi M, Robinot N, Lence-Anta JJ, Turcios S, et al. Common variants at the 9q22.33, 14q13.3 and ATM loci, and risk of differentiated thyroid cancer in the Cuban population. BMC Genetics. 2015; 16 :22. DOI: 10.1186/s12863-015-0180-5 - 88.
Estrada-Florez AP, Bohórquez ME, Sahasrabudhe R, Prieto R, Lott P, Duque CS, et al. Clinical features of Hispanic thyroid cancer cases and the role of known genetic variants on disease risk. Medicine (Baltimore). 2016; 95 (32):e4148. DOI: 10.1097/MD.0000000000004148 - 89.
Mussazhanova Z, Rogounovitch TI, Saenko VA, Krykpayeva A, Espenbetova M, Azizov B, et al. The contribution of genetic variants to the risk of papillary thyroid carcinoma in the kazakh population: Study of common single nucleotide polymorphisms and their clinicopathological correlations. Frontiers in Endocrinology (Lausanne). 2020; 11 :543500. DOI: 10.3389/fendo.2020.543500 - 90.
Landa I, Ruiz-Llorente S, Montero-Conde C, Inglada-Perez L, Schiavi F, Leskela S, et al. The variant rs1867277 in FOXE1 gene confers thyroid cancer susceptibility through the recruitment of USF1/USF2 transcription factors. PLoS Genetics. 2009; 5 (9):e1000637 - 91.
He H, Li W, Liyanarachchi S, Jendrzejewski J, Srinivas M, Davuluri RV, et al. Genetic predisposition to papillary thyroid carcinoma: Involvement of FOXE1, TSHR, and a novel lincRNA Gene, PTCSC2. The Journal of Clinical Endocrinology and Metabolism. 2015; 100 (1):E164-E172 - 92.
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 Genetics. 2015; 208 (3):68-75 - 93.
Jendrzejewski J, Liyanarachchi S, Eiterman A, Thomas A, He H, Nagy R, et al. Fine mapping of 14q13 reveals novel variants associated with different histological subtypes of papillary thyroid carcinoma. International Journal of Cancer. 2019; 144 (3):503-512. DOI: 10.1002/ijc.31933 - 94.
Rogounovitch TI, Bychkov A, Takahashi M, Mitsutake N, Nakashima M, Nikitski AV, et al. The common genetic variant rs944289 on chromosome 14q13.3 associates with risk of both malignant and benign thyroid tumors in the Japanese population. Thyroid. 2015; 25 (3):333-340 - 95.
Jendrzejewski J, He H, Radomska HS, Li W, Tomsic J, Liyanarachchi S, et al. The polymorphism rs944289 predisposes to papillary thyroid carcinoma through a large intergenic noncoding RNA gene of tumor suppressor type. Proceedings of the National Academy of Sciences of the United States of America. 2012; 109 (22):8646-8651. DOI: 10.1073/PNAS.1205654109 - 96.
Guibon J, Sugier PE, Kulkarni O, Karimi M, Bacq-Daian D, Besse C, et al. Fine–mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians. Oncotarget. 2021; 12 (5):493-506. DOI: 10.18632/ONCOTARGET.27888 - 97.
Michailidou K, Beesley J, Lindstrom S, Canisius S, Dennis J, Lush MJ, et al. Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer. Nature Genetics. 2015; 47 (4):373-380. DOI: 10.1038/ng.3242 - 98.
Liyanarachchi S, Wojcicka A, Li W, Czetwertynska M, Stachlewska E, Nagy R, et al. Cumulative risk impact of five genetic variants associated with papillary thyroid carcinoma. Thyroid. 2013; 23 (12):1532-1540. DOI: 10.1089/thy.2013.0102 - 99.
He H, Li W, Liyanarachchi S, Wang Y, Yu L, Genutis LK, et al. The role of NRG1 in the predisposition to papillary thyroid carcinoma. The Journal of Clinical Endocrinology and Metabolism. 2018; 103 (4):1369-1379. DOI: 10.1210/JC.2017-01798 - 100.
Pooley KA, Bojesen SE, Weischer M, Nielsen SF, Thompson D, Amin Al Olama A, et al. A genome-wide association scan (GWAS) for mean telomere length within the COGS project: Identified loci show little association with hormone-related cancer risk. Human Molecular Genetics. 2013; 22 (24):5056-5064. DOI: 10.1093/HMG/DDT355 - 101.
Li Y, Xiang C, Shen N, Deng L, Luo X, Yuan P, et al. Functional polymorphisms on chromosome 5p15.33 disturb telomere biology and confer the risk of non-small cell lung cancer in Chinese population. Molecular Carcinogenesis. 2019; 58 (6):913-921. DOI: 10.1002/MC.22980 - 102.
Rode L, Nordestgaard BG, Bojesen SE. Long telomeres and cancer risk among 95 568 individuals from the general population. International Journal of Epidemiology. 2016; 45 (5):1634-1643. DOI: 10.1093/IJE/DYW179 - 103.
Terry KL, Tworoger SS, Vitonis AF, Wong J, Titus-Ernstoff L, De Vivo I, et al. Telomere length and genetic variation in telomere maintenance genes in relation to ovarian cancer risk. Cancer Epidemiology, Biomarkers & Prevention. 2012; 21 (3):504-512 - 104.
Phelan CM, Kuchenbaecker KB, Tyrer JP, Kar SP, Lawrenson K, Winham SJ, et al. Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer. Nature Genetics. 2017; 49 (5):680-691. DOI: 10.1038/NG.3826 - 105.
Turner AW, Martinuk A, Silva A, Lau P, Nikpay M, Eriksson P, et al. Functional analysis of a novel genome-wide association study signal in SMAD3 that confers protection from coronary artery disease. Arteriosclerosis, Thrombosis, and Vascular Biology. 2016; 36 (5):972-983. DOI: 10.1161/ATVBAHA.116.307294 - 106.
Wu M, Assassi S, Salazar GA, Pedroza C, Gorlova OY, Chen WV, et al. Genetic susceptibility loci of idiopathic interstitial pneumonia do not represent risk for systemic sclerosis: A case control study in Caucasian patients. Arthritis Research & Therapy. 2016; 18 (1). DOI: 10.1186/S13075-016-0923-3 - 107.
Rosenthal EA, Crosslin DR, Gordon AS, Carrell DS, Stanaway IB, Larson EB, et al. Association between triglycerides, known risk SNVs and conserved rare variation in SLC25A40 in a multi-ancestry cohort. BMC Medical Genomics. 2021; 14 (1). DOI: 10.1186/S12920-020-00854-2 - 108.
Song SS, Huang S, Park S. Association of polygenetic risk scores related to cell differentiation and inflammation with thyroid cancer risk and genetic interaction with dietary intake. Cancers (Basel). 2021; 13 (7):1510. DOI: 10.3390/cancers13071510 - 109.
Hao R, Han P, Zhang L, Bi Y, Yan J, Li H, et al. Genetic polymorphisms in the PCNXL2 gene are risk factors for thyroid cancer in the Chinese population. Future Oncology. 2021; 17 (34):4677-4686. DOI: 10.2217/FON-2021-0748 - 110.
Ren W, Xu C, Wang S, Li H, Dai H, Yang F, et al. The effect of VAV3 polymorphisms on thyroid cancer. Endocrine. 2022; 75 (1):178-184. DOI: 10.1007/s12020-021-02827-6 - 111.
Liu Y, Shi F. Associations of telomerase reverse transcriptase rs10069690 and rs2736100 polymorphisms with papillary thyroid carcinoma. European Journal of Cancer Prevention. 2020:259-265. DOI: 10.1097/CEJ.0000000000000536 - 112.
He G, Song T, Zhang Y, Chen X, Xiong W, Chen H, et al. TERT rs10069690 polymorphism and cancers risk: A meta-analysis. Molecular Genetics & Genomic Medicine. 2019; 7 (10):e00903. DOI: 10.1002/MGG3.903 - 113.
dos Santos ICC, Genre J, Marques D, da Silva AMG, dos Santos JC, de Araújo JNG, et al. A new panel of SNPs to assess thyroid carcinoma risk: A pilot study in a Brazilian admixture population. BMC Medical Genetics. 2017; 18 (1):140. DOI: 10.1186/S12881-017-0502-8 - 114.
Figlioli G, Chen B, Elisei R, Romei C, Campo C, Cipollini M, et al. Novel genetic variants in differentiated thyroid cancer and assessment of the cumulative risk. Scientific Reports. 2015; 5 :8922 - 115.
Graff RE, Cavazos TB, Thai KK, Kachuri L, Rashkin SR, Hoffman JD, et al. Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts. Nature Communications. 2021; 12 (1). DOI: 10.1038/s41467-021-21288-z - 116.
Park J, Choi J yeob, Choi J, Chung S, Song N, Park SK, et al. Gene-environment interactions relevant to estrogen and risk of breast cancer: Can gene-environment interactions be detected only among candidate SNPs from genome-wide association studies? Cancers (Basel). 2021; 13 (10):2370 - 117.
Wang L, Desai H, Verma SS, Le A, Hausler R, Verma A, et al. Performance of polygenic risk scores for cancer prediction in a racially diverse academic biobank. Genetics in Medicine. 2022; 24 (3):601-609. DOI: 10.1016/j.gim.2021.10.015 - 118.
Liyanarachchi S, Gudmundsson J, Ferkingstad E, He H, Jonasson JG, Tragante V, et al. Assessing thyroid cancer risk using polygenic risk scores. Proceedings of the National Academy of Sciences of the United States of America. 2020; 117 (11):5997-6002. DOI: 10.1073/pnas.1919976117 - 119.
Song N, Liu Q , Wilson CL, Sapkota Y, Ehrhardt MJ, Gibson TM, et al. Polygenic risk score improves risk stratification and prediction of subsequent thyroid cancer after childhood cancer. Cancer Epidemiology, Biomarkers & Prevention. 2021; 30 (11):2096-2104. DOI: 10.1158/1055-9965.EPI-21-0448 - 120.
Jendrzejewski JP, Sworczak K, Comiskey DF, de la Chapelle A. Clinical implications of GWAS variants associated with differentiated thyroid cancer. Endokrynologia Polska. 2019; 70 (5):423-429. DOI: 10.5603/EP.a2019.0027