The FOXE1 locus is a major genetic determinant for radiation-related thyroid carcinoma in Chernobyl

Polymorphisms of DNA Repair Genes in Thyroid Cancer

The incidence of thyroid cancer, one of the most common forms of endocrine cancer, is increasing rapidly worldwide in developed and developing countries. Various risk factors can increase susceptibility to thyroid cancer, but particular emphasis is put on the role of DNA repair genes, which have a significant impact on genome stability. Polymorphisms of these genes can increase the risk of developing thyroid cancer by affecting their function. In this article, we present a concise review on the most common polymorphisms of selected DNA repair genes that may influence the risk of thyroid cancer. We point out significant differences in the frequency of these polymorphisms between various populations and their potential relationship with susceptibility to the disease. A more complete understanding of these differences may lead to the development of effective prevention strategies and targeted therapies for thyroid cancer. Simultaneously, there is a need for further research on the role of polymorphisms of previously uninvestigated DNA repair genes in the context of thyroid cancer, which may contribute to filling the knowledge gaps on this subject.

Chromosomal localization of mutated genes in non-syndromic familial thyroid cancer

Familial non-medullary thyroid carcinoma (FNMTC) is a type of thyroid cancer characterized by genetic susceptibility, representing approximately 5% of all non-medullary thyroid carcinomas. While some cases of FNMTC are associated with familial multi-organ tumor predisposition syndromes, the majority occur independently. The genetic mechanisms underlying non-syndromic FNMTC remain unclear. Initial studies utilized SNP linkage analysis to identify susceptibility loci, including the 1q21 locus, 2q21 locus, and 4q32 locus, among others. Subsequent research employed more advanced techniques such as Genome-wide Association Study and Whole Exome Sequencing, leading to the discovery of genes such as IMMP2L, GALNTL4, WDR11-AS1, DUOX2, NOP53, MAP2K5, and others. But FNMTC exhibits strong genetic heterogeneity, with each family having its own pathogenic genes. This is the first article to provide a chromosomal landscape map of susceptibility genes associated with non-syndromic FNMTC and analyze their potential associations. It also presents a detailed summary of variant loci, characteristics, research methodologies, and validation results from different countries.

Radiation-Related Thyroid Cancer

Radiation is an environmental factor that elevates the risk of developing thyroid cancer. Actual and possible scenarios of exposures to external and internal radiation are multiple and diverse. This article reviews radiation doses to the thyroid and corresponding cancer risks due to planned, existing, and emergency exposure situations, and medical, public, and occupational categories of exposures. Any exposure scenario may deliver a range of doses to the thyroid, and the risk for cancer is addressed along with modifying factors. The consequences of the Chornobyl and Fukushima nuclear power plant accidents are described, summarizing the information on thyroid cancer epidemiology, treatment, and prognosis, clinicopathological characteristics, and genetic alterations. The Chornobyl thyroid cancers have evolved in time: becoming less aggressive and driver shifting from fusions to point mutations. A comparison of thyroid cancers from the 2 areas reveals numerous differences that cumulatively suggest the low probability of the radiogenic nature of thyroid cancers in Fukushima. In view of continuing usage of different sources of radiation in various settings, the possible ways of reducing thyroid cancer risk from exposures are considered. For external exposures, reasonable measures are generally in line with the As Low As Reasonably Achievable principle, while for internal irradiation from radioactive iodine, thyroid blocking with stable iodine may be recommended in addition to other measures in case of anticipated exposures from a nuclear reactor accident. Finally, the perspectives of studies of radiation effects on the thyroid are discussed from the epidemiological, basic science, and clinical points of view.

Association of Polygenic Score With Tumor Molecular Subtypes in Papillary Thyroid Carcinoma

CONTEXT

Genome-wide association studies have identified germline variants associated with elevated PTC risk. It is also known that somatic driver mutations contribute to PTC development and as such PTCs can be further categorized into different molecular subtypes based on their somatic alterations. However, it remains unknown whether identified germline variants predictive of PTC risk are associated with specific molecular subtypes.

OBJECTIVE

The primary goal of the present study is to determine whether germline genetic risk, as assessed using a polygenic score (PGS) is associated with molecular subtypes of papillary thyroid carcinoma (PTC), defined based on tumor driver mutation status.

METHODS

This study was carried out using data from The Cancer Genome Atlas (TCGA) thyroid cancer study. A previously validated 10-single-nucleotide variation PGS for PTC derived from genome-wide association study hits was calculated to ascertain germline genetic risk. The primary molecular subtypes of interest were defined by tumor driver mutation status (BRAFV600E-mutated vs RAS-mutated vs "other"). We also explored associations between PGS and molecular subtypes defined by messenger RNA (mRNA) expression, microRNA expression, and DNA methylation patterns. Polytomous logistic regression analysis was used to assess the association between PGS and PTC molecular subtype with and without adjustment for clinical variables. Odds ratios (ORs) with their 95% CIs were estimated.

RESULTS

A total of 359 patients were included in the study. PGS was significantly associated specific tumor molecular subtypes defined by tumor driver mutation status. Increasing germline risk was associated with having a higher odd of BRAFV600E-mutated PTC compared to PTCs without driver mutations in the "other" category. No significant difference was detected in terms of PGS tumor categorization in the RAS subtype compared to BRAFV600E. In exploratory analyses, PGS was also associated with mRNA-, microRNA-, and DNA methylation-defined molecular subtypes, as defined by the TCGA PTC study.

CONCLUSION

PGS has molecular subtype-specific associations in PTC, which has implications for their use in risk prediction.

Pathogenesis of cancers derived from thyroid follicular cells

The genomic simplicity of differentiated cancers derived from thyroid follicular cells offers unique insights into how oncogenic drivers impact tumour phenotype. Essentially, the main oncoproteins in thyroid cancer activate nodes in the receptor tyrosine kinase-RAS-BRAF pathway, which constitutively induces MAPK signalling to varying degrees consistent with their specific biochemical mechanisms of action. The magnitude of the flux through the MAPK signalling pathway determines key elements of thyroid cancer biology, including differentiation state, invasive properties and the cellular composition of the tumour microenvironment. Progression of disease results from genomic lesions that drive immortalization, disrupt chromatin accessibility and cause cell cycle checkpoint dysfunction, in conjunction with a tumour microenvironment characterized by progressive immunosuppression. This Review charts the genomic trajectories of these common endocrine tumours, while connecting them to the biological states that they confer.

Genetic factors for differentiated thyroid cancer in French Polynesia: new candidate loci

Background

Populations of French Polynesia (FP), where France performed atmospheric tests between 1966 and 1974, experience a high incidence of differentiated thyroid cancer (DTC). However, up to now, no sufficiently large study of DTC genetic factors in this population has been performed to reach definitive conclusion. This research aimed to analyze the genetic factors of DTC risk among the native FP populations.

Methods

We analyzed more than 300 000 single nucleotide polymorphisms (SNPs) genotyped in 283 DTC cases and 418 matched controls born in FP, most being younger than 15 years old at the time of the first nuclear tests. We analyzed the genetic profile of our cohort to identify population subgroups. We then completed a genome-wide analysis study on the whole population.

Results

We identified a specific genetic structure in the FP population reflecting admixture from Asian and European populations. We identified three regions associated with increased DTC risk at 6q24.3, 10p12.2, and 17q21.32. The lead SNPs at these loci showed respective p-values of 1.66 × 10^-7, 2.39 × 10^-7, and 7.19 × 10^-7 and corresponding odds ratios of 2.02, 1.89, and 2.37.

Conclusion

Our study results suggest a role of the loci 6q24.3, 10p12.2 and 17q21.32 in DTC risk. However, a whole genome sequencing approach would be better suited to characterize these factors than genotyping with microarray chip designed for the Caucasian population. Moreover, the functional impact of these three new loci needs to be further explored and validated.

Genome-wide association studies for thyroid physiology and diseases

Genetic factors are involved in the etiology of most diseases, but prior to 2000, the methods for identifying such factors were very limited. Genome-wide association study (GWAS), developed in the 2000s, is an analytical method that can be applied to most diseases, including endocrine disorders. GWAS has provided a wealth of information on disease risks and the molecular pathogenesis of many human diseases. This review summarizes key findings from GWAS for thyroid physiology and diseases, and illustrates how GWAS is a powerful research tool to elucidate the molecular mechanisms of the diseases.

Adipocyte signaling affects thyroid-specific gene expression via down-regulation of TTF-2/FOXE1

Obesity affects thyroid gland function. Hypothyroidism, thyroid nodules, goiter, and thyroid cancer are more frequent in patients with higher BMI values. Although these data are supported by many clinical and epidemiological studies, our knowledge is very scarce at the molecular level. In this study, we present the first experimental evidence that adipocyte signaling downregulates the expression of thyroid-specific transcription factor 2 (TTF-2/FoxE1). It plays a crucial role in thyroid development and thyroid homeostasis and it is strictly connected to thyroid cancer as well. We provide in vivo and in vitro evidence that inhibition of TTF-2/FoxE1 gene expression is mediated by adipocyte signaling.

Foxe1 Deletion in the Adult Mouse Is Associated With Increased Thyroidal Mast Cells and Hypothyroidism

CONTEXT

Foxe1 is a key thyroid developmental transcription factor. Germline deletion results in athyreosis and congenital hypothyroidism. Some data suggest an ongoing role for maintaining thyroid differentiation.

OBJECTIVE

We created a mouse model to directly examine the role of Foxe1 in the adult thyroid.

METHODS

A model of tamoxifen-inducible Cre-mediated ubiquitous deletion of Foxe1 was generated in mice of C57BL/6J background (Foxe1flox/flox/Cre-TAM). Tamoxifen or vehicle was administered to Foxe1flox/flox/Cre mice aged 6-8 weeks. Blood was collected at 4, 12, and 20 weeks, and tissues after 12 or 20 weeks for molecular and histological analyses. Plasma total thyroxine (T4), triiodothyronine, and thyrotropin (TSH) were measured. Transcriptomics was performed using microarray or RNA-seq and validated by reverse transcription quantitative polymerase chain reaction.

RESULTS

Foxe1 was decreased by approximately 80% in Foxe1flox/flox/Cre-TAM mice and confirmed by immunohistochemistry. Foxe1 deletion was associated with abnormal follicular architecture and smaller follicle size at 12 and 20 weeks. Plasma TSH was elevated in Foxe1flox/flox/Cre-TAM mice as early as 4 weeks and T4 was lower in pooled samples from 12 and 20 weeks. Foxe1 deletion was also associated with an increase in thyroidal mast cells. Transcriptomic analyses found decreased Tpo and Tg and upregulated mast cell markers Mcpt4 and Ctsg in Foxe1flox/flox/Cre-TAM mice.

CONCLUSION

Foxe1 deletion in adult mice was associated with disruption in thyroid follicular architecture accompanied by biochemical hypothyroidism, confirming its role in maintenance of thyroid differentiation. An unanticipated finding was an increase in thyroidal mast cells. These data suggest a possible explanation for previous human genetic studies associating alleles in/near FOXE1 with hypothyroidism and/or autoimmune thyroiditis.

Genetic signature of differentiated thyroid carcinoma susceptibility: a machine learning approach

To identify a peculiar genetic combination predisposing to differentiated thyroid carcinoma (DTC), we selected a set of single nucleotide polymorphisms (SNPs) associated with DTC risk, considering polygenic risk score (PRS), Bayesian statistics and a machine learning (ML) classifier to describe cases and controls in three different datasets. Dataset 1 (649 DTC, 431 controls) has been previously genotyped in a genome-wide association study (GWAS) on Italian DTC. Dataset 2 (234 DTC, 101 controls) and dataset 3 (404 DTC, 392 controls) were genotyped. Associations of 171 SNPs reported to predispose to DTC in candidate studies were extracted from the GWAS of dataset 1, followed by replication of SNPs associated with DTC risk (P < 0.05) in dataset 2. The reliability of the identified SNPs was confirmed by PRS and Bayesian statistics after merging the three datasets. SNPs were used to describe the case/control state of individuals by ML classifier. Starting from 171 SNPs associated with DTC, 15 were positive in both datasets 1 and 2. Using these markers, PRS revealed that individuals in the fifth quintile had a seven-fold increased risk of DTC than those in the first. Bayesian inference confirmed that the selected 15 SNPs differentiate cases from controls. Results were corroborated by ML, finding a maximum AUC of about 0.7. A restricted selection of only 15 DTC-associated SNPs is able to describe the inner genetic structure of Italian individuals, and ML allows a fair prediction of case or control status based solely on the individual genetic background.

Investigating the Association Between rs2439302 Polymorphism and Thyroid Cancer: A Systematic Review and Meta-Analysis

Background and Aims

The extent of surgical treatment for most patients with thyroid cancer (TC) remains controversial and varies widely. As an emerging technology, genetic testing facilitates tumor typing and disease progression monitoring and is expected to influence the choice of surgical approach for patients with TC. Recent genome-wide association studies (GWASs) have identified that rs2439302 (8p12) variants near NRG1 are associated with TC risk; however, the results remain inconclusive. Therefore, we aimed to perform a meta-analysis to clarify the association between rs2439302 variants and the risk of TC.

Methods

We search eligible studies using Pubmed, Scopus, Embase, Web of Science, and Cochrane library by July 2021. We analyzed the pooled OR and the corresponding 95% confidence interval (95% CI) of the included studies and then conducted subgroup analysis according to the ethnicity. We also performed a sensitivity analysis to validate the findings.

Results

This meta-analysis finally included 7 studies involving 6,090 cases and 14,461 controls. Results showed that the G allele of the rs2439302 polymorphism was a significant risk factor of TC in Allele (G/C), Dominant (GG+GC/CC), Recessive (GG/GC+CC), Homozygote (GG/CC), Heterozygote (GC/CC) models, with pooled ORs of 1.38 (95%CI, 1.31–1.45), 1.51 (95%CI, 1.41–1.62), 1.52 (95%CI, 1.40–1.66), 1.90 (95%CI, 1.71–2.10), and 1.40 (95%CI, 1.30–1.51), respectively. The subgroup analysis showed that rs2439302 polymorphism was associated with higher TC risk in different ethnicities with OR > 1. The sensitivity analysis exhibited that the results were stable by omitting any included studies.

Conclusions

The study revealed that rs2439302 variants were associated with higher TC risk and may have a major influence on the choice of operative approach for patients with TC.

Genetic susceptibility to hereditary non-medullary thyroid cancer

Non-medullary thyroid cancer (NMTC) is the most common type of thyroid cancer. With the increasing incidence of NMTC in recent years, the familial form of the disease has also become more common than previously reported, accounting for 5-15% of NMTC cases. Familial NMTC is further classified as non-syndromic and the less common syndromic FNMTC. Although syndromic NMTC has well-known genetic risk factors, the gene(s) responsible for the vast majority of non-syndromic FNMTC cases are yet to be identified. To date, several candidate genes have been identified as susceptibility genes in hereditary NMTC. This review summarizes genetic predisposition to non-medullary thyroid cancer and expands on the role of genetic variants in thyroid cancer tumorigenesis and the level of penetrance of NMTC-susceptibility genes.

Age-related long-term response in rat thyroid tissue and plasma after internal low dose exposure to 131I

131I is used clinically for therapy, and may be released during nuclear accidents. After the Chernobyl accident papillary thyroid carcinoma incidence increased in children, but not adults. The aims of this study were to compare 131I irradiation-dependent differences in RNA and protein expression in the thyroid and plasma of young and adult rats, and identify potential age-dependent biomarkers for 131I exposure. Twelve young (5 weeks) and twelve adult Sprague Dawley rats (17 weeks) were i.v. injected with 50 kBq 131I (absorbed dose to thyroid = 0.1 Gy), and sixteen unexposed age-matched rats were used as controls. The rats were killed 3-9 months after administration. Microarray analysis was performed using RNA from thyroid samples, while LC-MS/MS analysis was performed on proteins extracted from thyroid tissue and plasma. Canonical pathways, biological functions and upstream regulators were analysed for the identified transcripts and proteins. Distinct age-dependent differences in gene and protein expression were observed. Novel biomarkers for thyroid 131I exposure were identified: (PTH), age-dependent dose response (CA1, FTL1, PVALB (youngsters) and HSPB6 (adults)), thyroid function (Vegfb (adults)). Further validation using clinical samples are needed to explore the role of the identified biomarkers.

Identification of NID1 as a novel candidate susceptibility gene for familial non-medullary thyroid carcinoma using whole-exome sequencing

The genetics underlying non-syndromic familial non-medullary thyroid carcinoma (FNMTC) is still poorly understood. To identify susceptibility genes for FNMTC, we performed whole-exome sequencing (WES) in a Brazilian family affected by papillary thyroid carcinoma (PTC) in three consecutive generations. WES was performed in four affected and two unaffected family members. Manual inspection in over 100 previously reported susceptibility genes for FNMTC showed that no variants in known genes co-segregated with disease phenotype in this family. Novel candidate genes were investigated using PhenoDB and filtered using Genome Aggregation (gnomAD) and Online Archive of Brazilian Mutations (ABraOM) population databases. The missense variant p.Ile657Met in the NID1 gene was the only variant that co-segregated with the disease, while absent in unaffected family members and controls. The allele frequency for this variant was <0.0001 in the gnomAD and ABbraOM databases. In silico analysis predicted the variant to be deleterious or likely damaging to the protein function. Somatic mutations in NID1 gene were found in nearly 500 cases of different cancer subtypes in the intOGen platform. Immunohistochemistry analysis showed NID1 expression in PTC cells, while it was absent in normal thyroid tissue. Our findings were corroborated using data from the TCGA cohort. Moreover, higher expression of NID1 was associated with higher likelihood of relapse after treatment and N1b disease in PTCs from the TCGA cohort. Although replication studies are needed to better understand the role of this variant in the FNMTC susceptibility, the NID1 variant (c.1971T>G) identified in this study fulfills several criteria that suggest it as a new FNMTC predisposing gene.

Identification of Genetic Predispositions Related to Ionizing Radiation in Primary Human Skin Fibroblasts From Survivors of Childhood and Second Primary Cancer as Well as Cancer-Free Controls: Protocol for the Nested Case-Control Study KiKme

BACKGROUND

Therapy for a first primary neoplasm (FPN) in childhood with high doses of ionizing radiation is an established risk factor for second primary neoplasms (SPN). An association between exposure to low doses and childhood cancer is also suggested; however, results are inconsistent. As only subgroups of children with FPNs develop SPNs, an interaction between radiation, genetic, and other risk factors is presumed to influence cancer development.

OBJECTIVE

Therefore, the population-based, nested case-control study KiKme aims to identify differences in genetic predisposition and radiation response between childhood cancer survivors with and without SPNs as well as cancer-free controls.

METHODS

We conducted a population-based, nested case-control study KiKme. Besides questionnaire information, skin biopsies and saliva samples are available. By measuring individual reactions to different exposures to radiation (eg, 0.05 and 2 Gray) in normal somatic cells of the same person, our design enables us to create several exposure scenarios for the same person simultaneously and measure several different molecular markers (eg, DNA, messenger RNA, long noncoding RNA, copy number variation).

RESULTS

Since 2013, 101 of 247 invited SPN patients, 340 of 1729 invited FPN patients, and 150 of 246 invited cancer-free controls were recruited and matched by age and sex. Childhood cancer patients were additionally matched by tumor morphology, year of diagnosis, and age at diagnosis. Participants reported on lifestyle, socioeconomical, and anthropometric factors, as well as on medical radiation history, health, and family history of diseases (n=556). Primary human fibroblasts from skin biopsies of the participants were cultivated (n=499) and cryopreserved (n=3886). DNA was extracted from fibroblasts (n=488) and saliva (n=510).

CONCLUSIONS

This molecular-epidemiological study is the first to combine observational epidemiological research with standardized experimental components in primary human skin fibroblasts to identify genetic predispositions related to ionizing radiation in childhood and SPNs. In the future, fibroblasts of the participants will be used for standardized irradiation experiments, which will inform analysis of the case-control study and vice versa. Differences between participants will be identified using several molecular markers. With its innovative combination of experimental and observational components, this new study will provide valuable data to forward research on radiation-related risk factors in childhood cancer and SPNs.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/32395.

Thyroid dose estimates for the genome-wide association study of thyroid cancer in persons exposed in Belarus to 131I after the Chernobyl accident

The Chernobyl accident on 26 April 1986 led to a sharp increase in thyroid cancer (TC) incidence in the individuals exposed to radiation in childhood. The major risk factor for TC was exposure to Iodine-131 (131I). Here, we estimated the thyroid doses due to 131I intake for 2041 participants of the genome-wide association study of TC in Belarusian people exposed to radioactive fallout from the Chernobyl accident. The following parameter-values specially developed in this study were used to estimate individual thyroid doses: (i) scaling factors for adjustment of the model-based doses, (ii) age and gender diet to characterize 131I intake, and (iii) area-, age- and gender-specific S-values for the thyroid gland per 131I decay in the thyroid. The most reliable doses were calculated for 103 people with measured 131I thyroid activity (the arithmetic mean of 1.2 Gy, median 0.52 Gy), and 275 individuals with detailed residential history and dietary data (the arithmetic mean of 0.41 Gy, median 0.24 Gy). The arithmetic mean of thyroid doses among all study participants was 0.23 Gy (median 0.082 Gy); the highest individual dose was 9.0 Gy. Special attention was paid to the reliability and validity of the obtained estimates, in particular for the individuals without 131I thyroid activity measurements and individual data on residential history and diet, by comparing those with the doses from other post-Chernobyl epidemiological studies. Overall, the doses estimated in the current study were in reasonable agreement with previously reported thyroid doses. These doses will be used in the genome-wide association study of TC in people exposed in Belarus to 131I after the Chernobyl accident.

Polygenic Risk Score Improves Risk Stratification and Prediction of Subsequent Thyroid Cancer after Childhood Cancer

BACKGROUND

Subsequent thyroid cancer (STC) is one of the most common malignancies in childhood cancer survivors. We aimed to evaluate the polygenic contributions to STC risk and potential utility in improving risk prediction.

METHODS

A polygenic risk score (PRS) was calculated from 12 independent SNPs associated with thyroid cancer risk in the general population. Associations between PRS and STC risk were evaluated among survivors from St. Jude Lifetime Cohort (SJLIFE) and were replicated in survivors from Childhood Cancer Survivor Study (CCSS). A risk prediction model integrating the PRS and clinical factors, initially developed in SJLIFE, and its performance were validated in CCSS.

RESULTS

Among 2,370 SJLIFE survivors with a median follow-up of 28.8 [interquartile range (IQR) = 21.9-36.1] years, 65 (2.7%) developed STC. Among them, the standardized PRS was associated with an increased rate of STC [relative rate (RR) = 1.57; 95% confidence interval (CI) = 1.24-1.98; P < 0.001]. Similar associations were replicated in 6,416 CCSS survivors, among whom 121 (1.9%) developed STC during median follow-up of 28.9 (IQR = 22.6-34.6) years (RR = 1.52; 95% CI = 1.25-1.83; P < 0.001). A risk prediction model integrating the PRS with clinical factors showed better performance than the model considering only clinical factors in SJLIFE (P = 0.004, AUC = 83.2% vs. 82.1%, at age 40), which was further validated in CCSS (P = 0.010, AUC = 72.9% vs. 70.6%).

CONCLUSIONS

Integration of the PRS with clinical factors provided a statistically significant improvement in risk prediction of STC, although the magnitude of improvement was modest.

IMPACT

PRS improves risk stratification and prediction of STC, suggesting its potential utility for optimizing screening strategies in survivorship care.

FOXE1-Dependent Regulation of Macrophage Chemotaxis by Thyroid Cells In Vitro and In Vivo

Forkhead box E1 (FOXE1) is a lineage-restricted transcription factor involved in thyroid cancer susceptibility. Cancer-associated polymorphisms map in regulatory regions, thus affecting the extent of gene expression. We have recently shown that genetic reduction of FOXE1 dosage modifies multiple thyroid cancer phenotypes. To identify relevant effectors playing roles in thyroid cancer development, here we analyse FOXE1-induced transcriptional alterations in thyroid cells that do not express endogenous FOXE1. Expression of FOXE1 elicits cell migration, while transcriptome analysis reveals that several immune cells-related categories are highly enriched in differentially expressed genes, including several upregulated chemokines involved in macrophage recruitment. Accordingly, FOXE1-expressing cells induce chemotaxis of co-cultured monocytes. We then asked if FOXE1 was able to regulate macrophage infiltration in thyroid cancers in vivo by using a mouse model of cancer, either wild type or with only one functional FOXE1 allele. Expression of the same set of chemokines directly correlates with FOXE1 dosage, and pro-tumourigenic M2 macrophage infiltration is decreased in tumours with reduced FOXE1. These data establish a novel link between FOXE1 and macrophages recruitment in the thyroid cancer microenvironment, highlighting an unsuspected function of this gene in the crosstalk between neoplastic and immune cells that shape tumour development and progression.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians

Differentiated thyroid carcinoma (DTC) incidence is characterized by wide ethnic and geographic variations, with high incidence rates observed in Oceanian populations. Genome-wide association studies (GWAS) identified mainly four DTC susceptibility loci at 9q22.33, 14q13.3, 2q35 and 8p12. Here we performed fine-mapping of the 2q35 and 8p12 loci in the population of the EPITHYR consortium that includes Europeans, Melanesians and Polynesians to identify likely causal variants for DTC risk. We conducted a colocalization analysis using eQTLs data to determine the SNPs with the highest probability of causality. At 2q35, we highlighted rs16857609 located in DIRC3. This SNP has a high probability of causality in the three populations, and a significant association in Europeans (OR = 1.4, p = 1.9 x 10-10). It is also associated with expression of DIRC3 and of the nearby gene IGFBP5 in thyroid tumour cells. At 8p12, we identified rs7844425 which was significantly associated with DTC in Europeans (OR = 1.32, p = 7.6 x 10-8) and rs2439304, which was highlighted by the colocalization analysis but only moderately associated with DTC in our dataset (OR = 1.2, p = 0.001). These SNPs are linked to the expression of NRG1 in thyroid tissue. Hence, our study identified novel variants at 2q35 and 8p12 to be prioritized for further functional studies.

Consequences of atmospheric contamination by radioiodine: the Chernobyl and Fukushima accidents

PURPOSE

After the accidents of nuclear power plants at Chernobyl and at Fukushima, huge amounts of radioactive iodine were released into the atmosphere.

METHODS

We reviewed data on the health consequences of these accidents with a focus on thyroid consequences.

RESULTS

Among the 2 million children who were living in highly contaminated regions in Belarus, Ukraine and Russia, 7000 cases of thyroid cancer had occurred in 2005. This is the most significant radiation-induced consequence of the Chernobyl accident. The increased incidence of thyroid cancer observed in adult population who lived in these highly contaminated regions is at least in major part related to screening and it is not possible to individualize among these thyroid cancers those that are potentially caused by radiation exposure. For populations who lived outside these regions at the time of the accident, there is no detectable consequence of the radiation exposure on the thyroid gland. Among children who lived nearby the Fukushima power plant in 2011, there is currently no evidence of an increased incidence of thyroid cancer. Ultrasonography screening in these individuals detected a number of thyroid cancers that are probably not related to the accident. Because thyroid cancer is frequent, studies have been carried out to distinguish radiation-induced from their sporadic counterparts, and genomic signatures might be helpful.

CONCLUSIONS

The consequences of the Chernobyl accident clearly demonstrate that populations living nearby a nuclear power plant should be protected in case of accident by sheltering, food restrictions and prophylaxis of thyroid irradiation by potassium iodine administration, if the predicted estimated dose to the thyroid gland of children might be >50 mGy. These countermeasures should be applied in priority to children, adolescents and pregnant women; they are safe and effective.

FOXE1 Gene Dosage Affects Thyroid Cancer Histology and Differentiation In Vivo

The transcription factor Forkhead box E1 (FOXE1) is a key player in thyroid development and function and has been identified by genome-wide association studies as a susceptibility gene for papillary thyroid cancer. Several cancer-associated polymorphisms fall into gene regulatory regions and are likely to affect FOXE1 expression levels. However, the possibility that changes in FOXE1 expression modulate thyroid cancer development has not been investigated. Here, we describe the effects of FOXE1 gene dosage reduction on cancer phenotype in vivo. Mice heterozygous for FOXE1 null allele (FOXE1^+/−) were crossed with a BRAF^V600E-inducible cancer model to develop thyroid cancer in either a FOXE1^+/+ or FOXE1^+/− genetic background. In FOXE1^+/+ mice, cancer histological features are quite similar to that of human high-grade papillary thyroid carcinomas, while cancers developed with reduced FOXE1 gene dosage maintain morphological features resembling less malignant thyroid cancers, showing reduced proliferation index and increased apoptosis as well. Such cancers, however, appear severely undifferentiated, indicating that FOXE1 levels affect thyroid differentiation during neoplastic transformation. These results show that FOXE1 dosage exerts pleiotropic effects on thyroid cancer phenotype by affecting histology and regulating key markers of tumor differentiation and progression, thus suggesting the possibility that FOXE1 could behave as lineage-specific oncogene in follicular cell-derived thyroid cancer.

Genetic Mutations and Variants in the Susceptibility of Familial Non-Medullary Thyroid Cancer

Thyroid cancer is the most frequent endocrine malignancy with the majority of cases derived from thyroid follicular cells and caused by sporadic mutations. However, when at least two or more first degree relatives present thyroid cancer, it is classified as familial non-medullary thyroid cancer (FNMTC) that may comprise 3-9% of all thyroid cancer. In this context, 5% of FNMTC are related to hereditary syndromes such as Cowden and Werner Syndromes, displaying specific genetic predisposition factors. On the other hand, the other 95% of cases are classified as non-syndromic FNMTC. Over the last 20 years, several candidate genes emerged in different studies of families worldwide. Nevertheless, the identification of a prevalent polymorphism or germinative mutation has not progressed in FNMTC. In this work, an overview of genetic alteration related to syndromic and non-syndromic FNMTC is presented.

Network-Based Analysis Reveals Association of FOXE1 Gene Polymorphisms in Thyroid Cancer Patients; A Case-Control Study in Southeast of Iran

Thyroid cancer (TC) is the mainly frequent endocrine cancer by different incidence rate in worldwide. However, early prediction of  this cancer is still challenging due to the unclear pathogenicity. In this study with the aid of systems biology approach, performed a holistic study on GSE65144 dataset containing anaplastic thyroid carcinoma tissues. Co-expression network analysis by WGCNA suggested that highly preserved turquoise module with 1,480 genes was significantly correlated to TC. Most of the top 54 hub-genes of this module are functionality correlated to thyroid hormone generation (GO:0006590). Of these 54 hub-genes, FOXE1 has been reported previously to contain mutation asosiated to TC and chosen for experimental validation step. To this end, we conducted a case-control study including 81 TC patients and 165 controls individuals to evaluate the effects of FOXE1 functional polymorphisms (rs1867277) on the development of TC in Sistan and Balouchestan province of Iran. The polymorphisms of FOXE1 gene (rs1867277) assessed by tetra-ARMS PCR technique. Homozygous (GG) and (AA) variant of rs1867277 polymorphism were detected in 26 (32.1%) and 15 (18.5 %) of TC patients, and 66 (40.0%), and 15 (9.1%) in controls, respectively (p-value= 0.03, OR= 2.53). The A allele frequency was 70 (43.2%) in TC patients and 114 (34.5%) in controls (p-value= 0.06, OR= 1.44). Overall, our results suggested that FOXE1 gene could be used as a prognostic marker in TC and also provides information related to FOXE1 functional polymorphisms (rs1867277) in Sistan and Balouchestan province of Iran.
.

Primary cell cultures for the personalized therapy in aggressive thyroid cancer of follicular origin

Thyroid cancer (TC) is the most prevalent endocrine malignancy. More than 90 % of TC is represented by differentiated TC (DTC) arising from the follicular thyroid cells. DTC includes papillary TC (PTC), follicular TC (FTC), and Hürthle cell TC. Anaplastic TC (ATC) accounts for 1% of TC, and it represents 15-40 % of TC death. Current treatment strategies are not completely effective against aggressive DTC or ATC, and mortality is one of the most important challenges. Recently, progresses have been obtained in the understanding of the molecular/genetic basis of TC progression, and new drugs have been introduced [i.e. tyrosine kinase inhibitors (TKIs)], able to block the oncogenic or signaling kinases, associated with cellular growth. Thyroid cell lines, obtained from tumoral cells and chosen for high proliferation in vitro, have been used as preclinical models. Actually, these cells lose the characteristic features of the primary tumor, because they adapt to in vitro growth conditions. For these reasons, the use of these cell lines has important limitations, and more recently human primary cell cultures have been established as monolayer cultures, and investigated for their biological behavior. Moreover, in the past, primary TC cells could be collected only through surgical biopsies, while recently human primary cell cultures can be established also from samples of fine-needle aspiration citology from aggressive dedifferentiated DTC or ATC. Testing in vitro different TKIs in each patient can help to develop new personalized treatments, without using ineffective drugs. In conclusion, personalized medicine and precise oncology, which consider both patients and their disease features, represent the future of the treatment approach, and further progress is needed in this direction.

Association between forkhead box E1 polymorphisms and risk of non-syndromic cleft lip with or without cleft palate: A meta-analysis

OBJECTIVE

The purpose of the present work was to investigate the association between forkhead box E1 (FOXE1) and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Relevant studies were searched in several professional databases up to 31 July 2019. The pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated using a fixed-effect model or a random-effect model to analyse the relationship between FOXE1 polymorphisms and NSCL/P.

RESULTS

A total of four single nucleotide polymorphisms (SNPs), including rs3758249, rs4460498, rs1443434 and rs10217225, were analysed. The overall findings showed that FOXE1 rs4460498 was statistically associated with NSCL/P (including cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO)). Genotypes CC and CT of rs4460498 were significantly more closely correlated with NSCL/P (including CL/P and CPO) than genotype TT (NSCL/P: TT vs CC, OR = 0.630, P = .000; TT vs TC + CC, OR = 0.775, P = .020; CL/P: TT vs CC, OR = 0.664, P = .000; TT vs TC + CC, OR = 0.738, P = .006. CPO: TT vs CC, OR = 0.761, P = .027; TT vs TC + CC, OR = 0.792, P = .045). For rs10217225, only the TT genotype might have contributed to the elevated risk of CL/P (TT vs CC OR = 2.236, P = .000). The other FOXE1 polymorphisms were not associated with NSCLP, CL/P or CPO.

CONCLUSION

The meta-analysis provided confirmation that the polymorphism of FOXE1 rs10217225 was correlated with an increased risk of CL/P, and the polymorphism of FOXE1 rs4460498 was a protective factor for NSCL/P, including CLP and CPO.

Machine learning on genome-wide association studies to predict the risk of radiation-associated contralateral breast cancer in the WECARE Study

The purpose of this study was to identify germline single nucleotide polymorphisms (SNPs) that optimally predict radiation-associated contralateral breast cancer (RCBC) and to provide new biological insights into the carcinogenic process. Fifty-two women with contralateral breast cancer and 153 women with unilateral breast cancer were identified within the Women’s Environmental Cancer and Radiation Epidemiology (WECARE) Study who were at increased risk of RCBC because they were ≤ 40 years of age at first diagnosis of breast cancer and received a scatter radiation dose > 1 Gy to the contralateral breast. A previously reported algorithm, preconditioned random forest regression, was applied to predict the risk of developing RCBC. The resulting model produced an area under the curve (AUC) of 0.62 (p = 0.04) on hold-out validation data. The biological analysis identified the cyclic AMP-mediated signaling and Ephrin-A as significant biological correlates, which were previously shown to influence cell survival after radiation in an ATM-dependent manner. The key connected genes and proteins that are identified in this analysis were previously identified as relevant to breast cancer, radiation response, or both. In summary, machine learning/bioinformatics methods applied to genome-wide genotyping data have great potential to reveal plausible biological correlates associated with the risk of RCBC.

Association of rs944289, rs965513, and rs1443434 in TITF1/TITF2 with Risks of Papillary Thyroid Carcinoma and with Nodular Goiter in Northern Chinese Han Populations

OBJECTIVE

In this study, we aimed to investigate the associations of three single-nucleotide polymorphisms (SNPs) on TITF1/TITF2 (rs944289, rs965513, and rs1443434) with susceptibility to papillary thyroid carcinoma (PTC) and with nodular goiter (NG) in northern Chinese Han populations.

METHODS

We performed a case-control study comprising 861 PTC patients, 562 NG patients, and 896 normal controls (NCs). One TITF1 SNP (rs944289) and two TITF2 SNPs (rs965513 and rs1443434) were genotyped. Departures from Hardy-Weinberg equilibrium (HWE) in the control group were evaluated using chi-square test. Associations of the SNPs with PTC and with NG were assessed by unconditional logistic regression using the online SNPStats program. Bonferroni correction was performed for multiple tests in genotype analyses. Data analysis was performed by SPSS24.0 unless otherwise specified.

RESULTS

For rs944289, T allele was associated with increased risks for both PTC (OR = 1.23, 95% CI: 1.08-1.41, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50, P=0.002) and NG (OR = 1.28, 95% CI: 1.10-1.50.

CONCLUSIONS

There are associations of rs944289 and rs1443434 polymorphisms with PTC risk and association of rs944289 polymorphism with NG risk. Haplotypes T-G-G and T-G-T are risk haplotypes of PTC and NG, respectively.

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

OBJECTIVE

Risk for developing papillary thyroid carcinoma (PTC), the most common endocrine malignancy, is thought to be mediated by lifestyle, environmental exposures and genetic factors. Recent progress in the genome-wide association studies of thyroid cancer leads to the identification of several genetic variants conferring risk to this malignancy across different ethnicities. We set out to elucidate the impact of selected single nucleotide polymorphisms (SNPs) on PTC risk and to evaluate clinicopathological correlations of these genetic variants in the Kazakh population for the first time.

METHODS

Eight SNPs were genotyped in 485 patients with PTC and 1,008 healthy control Kazakh subjects. The association analysis and multivariable modeling of PTC risk by the genetic factors, supplemented with rigorous statistical validation, were performed.

RESULT

Five of the eight SNPs: rs965513 (FOXE1/PTCSC2, P = 1.3E-16), rs1867277 (FOXE1 5'UTR, P = 7.5E-06), rs2439302 (NRG1 intron 1, P = 4.0E-05), rs944289 (PTCSC3/NKX2-1, P = 4.5E-06) and rs10136427 (BATF upstream, P = 9.8E-03) were significantly associated with PTC. rs966423 (DIRC3, P = 0.07) showed a suggestive association. rs7267944 (DHX35) was associated with PTC risk in males (P = 0.02), rs1867277 (FOXE1) conferred the higher risk in subjects older than 55 years (P = 7.0E-05), and rs6983267 (POU5F1B/CCAT2) was associated with pT3-T4 tumors (P = 0.01). The contribution of genetic component (unidirectional independent effects of rs965513, rs944289, rs2439302 and rs10136427 adjusted for age and sex) to PTC risk in the analyzed series was estimated to be 30-40%.

CONCLUSION

Genetic factors analyzed in the present work display significant association signals with PTC either on the whole group analysis or in particular clinicopathological groups and account for about one-third of the risk for PTC in the Kazakh population.

Radiation-Induced Thyroid Cancers: Overview of Molecular Signatures

Enormous amounts of childhood thyroid cancers, mostly childhood papillary thyroid carcinomas (PTCs), after the Chernobyl nuclear power plant accident have revealed a mutual relationship between the radiation exposure and thyroid cancer development. While the internal exposure to radioactive ¹³¹I is involved in the childhood thyroid cancers after the Chernobyl accident, people exposed to the external radiation, such as atomic-bomb (A-bomb) survivors, and the patients who received radiation therapy, have also been epidemiologically demonstrated to develop thyroid cancers. In order to elucidate the mechanisms of radiation-induced carcinogenesis, studies have aimed at defining the molecular changes associated with the thyroid cancer development. Here, we overview the literatures towards the identification of oncogenic alterations, particularly gene rearrangements, and discuss the existence of radiation signatures associated with radiation-induced thyroid cancers.

Targeted next-generation sequencing in papillary thyroid carcinoma patients looking for germline variants predisposing to the disease

PURPOSESS

The purpose of this study was using next-generation sequencing technique to explore the potential association between germline variants of 14 targeted genes and papillary thyroid carcinoma (PTC) predisposition as well as disease progression.

METHODS

In all, 516 subjects were enrolled in this study including 416 PTC patients and 100 healthy controls. PTC patients were divided into distant metastasis group and non-distant metastasis group. Patients in distant metastasis group were further divided into radioiodine-refractory PTC (RR-PTC) and non-RR-PTC depending on their response to radioiodine therapy. Genomic DNA was extracted from peripheral blood sample and MiSeq Benchtop Sequencer was used for sequencing.

RESULTS

We found rs11246050 in NLRP6 (dominant model, OR/95% CI: 2.028/1.091-3.769, p = 0.025), rs2286742 and rs3740530 in HABP2 (recessive model, OR/95% CI: 9.644/1.307-71.16, p = 0.026 and 3.989/1.413-11.26, p = 0.009), rs2736098 in TERT (recessive model, OR/95% CI: 2.322/1.028-5.242. p = 0.042) and rs62054619 in GAS8-AS1 (recessive model, OR/95% CI: 2.219/1.067-4.617, p = 0.033) were associated with the risk of PTC. rs1137282 in KRAS (dominant model, OR/95% CI: 0.5430/0.3192-0.9236, p = 0.024), rs1347591 and rs4461062 in NUP93 (dominant model, OR/95% CI: 0.6121/0.4128-0.9076, p = 0.015 and 0.6156/0.4157-0.9117, p = 0.015) were associated with low risk of distant metastatic disease in PTC patients. rs33954691 in TERT was associated with the risk of RR-PTC under dominant model (OR/95% CI: 3.161/1.596-6.262).

CONCLUSIONS

Germline variants of related genes could be associated with the susceptibility of PTC as well as disease progression (distant metastasis and radioiodine-refractory status). However, these results must be further verified and the potential biological functions of these germline variants in the pathogenesis of PTC remain to be determined in future studies.

Genome-Wide Association Study Reveals Distinct Genetic Susceptibility of Thyroid Nodules From Thyroid Cancer

CONTEXT

Thyroid nodules are very common, and 7% to 15% of them are diagnosed as thyroid cancer. However, the inherited genetic risk factors for thyroid nodules and their associations with thyroid cancer remain unknown.

OBJECTIVE

To identify the genetic variants associated with susceptibility to thyroid nodules in comparison with thyroid cancer.

DESIGN AND SETTING

We performed a three-stage genome-wide association study for thyroid nodules. The discovery stage involved a genome-wide scan of 811 subjects with thyroid nodules and 691 subjects with a normal thyroid from a population-based cohort. Replication studies were conducted in an additional 1981 cases and 3100 controls from the participants of a health checkup. We also performed expression quantitative trait loci analysis of public data.

RESULTS

The most robust association was observed in TRPM3 (rs4745021) in the joint analysis (OR, 1.26; P = 6.12 × 10-8) and meta-analysis (OR, 1.28; P = 2.11 × 10-8). Signals at MBIP/NKX2-1 were replicated but did not reach genome-wide significance in the joint analysis (rs2415317, P = 4.62 × 10-5; rs944289, P = 8.68 × 10-5). The expression quantitative trait loci analysis showed that TRPM3 expression was associated with the rs4745021 genotype in thyroid tissues.

CONCLUSIONS

To the best of our knowledge, we have performed the first genome-wide association study of thyroid nodules and identified a susceptibility locus associated with thyroid nodules, suggesting that thyroid nodules have a genetic predisposition distinct from that of thyroid cancer.

Incidence of Second Malignancy in Patients with Papillary Thyroid Cancer from Surveillance, Epidemiology, and End Results 13 Dataset

Increased risk of second primary malignancy (SPM) in papillary thyroid cancer (PTC) has been reported. Here, we present the most updated incidence rates of second primary malignancy from original diagnosis of PTC by using the data from the Surveillance, Epidemiology, and End Results. In this cohort, 3,200 patients developed SPM, a substantially higher number than in the reference population of 2,749 with observed to expected ratio (O/E) of 1.16 (95% CI; 1.12–1.21). Bone and joint cancer had the highest O/E ratio of 4.26 (95% confidence interval [CI] 2.33–7.15) followed by salivary gland (O/E 4.15; 95% CI 2.76–6.0) and acute lymphocytic leukemia (O/E 3.98; 95% CI 2.12–6.8). Mean age at the diagnosis of SPM was 64.4 years old. Interestingly, incidence of colorectal cancer was lower in thyroid cancer survivors compared to general population (large intestine O/E 0.3; 95% CI 0.06–0.88, rectum O/E 0.6; 95% CI 0.41–0.85); however, this was not observed in patients who underwent radiation therapy. The incidence of SPM at all sites was higher during 2000–2012 compared to 1992–1999 (O/E 1.24 versus 1.10). Surprisingly, patients with micropapillary cancer had higher incidence of SPM than counterparts with a larger tumor in radiation group (O/E of 1.40 versus 1.15). O/E of all cancers were higher in males compared to females with O/E of 1.41 versus 1.17 during the period of 2000–2012. Diagnosis of PTC before age 50, especially at age 30–34, was associated with higher incidence of overall SPM (age 30–34; O/E 1.43; 95% CI; 1.19–1.71). Efficient monitoring strategies that include age at the time of thyroid cancer diagnosis, exposure to radiation, gender, and genetic susceptibility may successfully detect SPM earlier in the disease course. This is especially important given the excellent prognosis of the initial thyroid cancer itself.

Genome-Wide Association Studies of Autoimmune Thyroid Diseases, Thyroid Function, and Thyroid Cancer

Thyroid diseases, including autoimmune thyroid diseases and thyroid cancer, are known to have high heritability. Family and twin studies have indicated that genetics plays a major role in the development of thyroid diseases. Thyroid function, represented by thyroid stimulating hormone (TSH) and free thyroxine (T4), is also known to be partly genetically determined. Before the era of genome-wide association studies (GWAS), the ability to identify genes responsible for susceptibility to thyroid disease was limited. Over the past decade, GWAS have been used to identify genes involved in many complex diseases, including various phenotypes of the thyroid gland. In GWAS of autoimmune thyroid diseases, many susceptibility loci associated with autoimmunity (human leukocyte antigen [HLA], protein tyrosine phosphatase, non-receptor type 22 [PTPN22], cytotoxic T-lymphocyte associated protein 4 [CTLA4], and interleukin 2 receptor subunit alpha [IL2RA]) or thyroid-specific genes (thyroid stimulating hormone receptor [TSHR] and forkhead box E1 [FOXE1]) have been identified. Regarding thyroid function, many susceptibility loci for levels of TSH and free T4 have been identified through genome-wide analyses. In GWAS of differentiated thyroid cancer, associations at FOXE1, MAP3K12 binding inhibitory protein 1 (MBIP)-NK2 homeobox 1 (NKX2-1), disrupted in renal carcinoma 3 (DIRC3), neuregulin 1 (NRG1), and pecanex-like 2 (PCNXL2) have been commonly identified in people of European and Korean ancestry, and many other susceptibility loci have been found in specific populations. Through GWAS of various thyroid-related phenotypes, many susceptibility loci have been found, providing insights into the pathogenesis of thyroid diseases and disease co-clustering within families and individuals.

Genetic Polymorphism Predisposing to Differentiated Thyroid Cancer: A Review of Major Findings of the Genome-Wide Association Studies

Thyroid cancer has one of the highest hereditary component among human malignancies as seen in medical epidemiology investigations, suggesting the potential meaningfulness of genetic studies. Here we review researches into genetic variations that influence the chance of developing non-familial differentiated thyroid cancer (DTC), focusing on the major findings of the genome-wide association studies (GWASs) of common single-nucleotide polymorphisms (SNPs). To date, eight GWAS have been performed, and the association of a number of SNPs have been reproduced in dozens of replication investigations across different ethnicities, including Korea and Japan. Despite the cumulative effect of the strongest SNPs demonstrates gradual increase in the risk for cancer and their association signals are statistically quite significant, the overall prediction ability for DTC appears to be very limited. Thus, genotyping of common SNPs only would be insufficient for evidence-based counseling in clinical setting at present. Further studies to include less significant and rare SNPs, non-SNP genetic information, gene-gene interactions, ethnicity, non-genetic and environmental factors, and development of more advanced computational algorithms are warranted to approach to personalized disease risk prediction and prognostication.

Role of biomarkers in predicting the occurrence of thyroid neoplasms in radiation-exposed children

With increasing numbers of childhood cancer survivors who were treated with radiation, there is a need to evaluate potential biomarkers that could signal an increased risk of developing thyroid cancer. We aimed to examine the relationships between thyrotropin and thyroglobulin levels and the risk of developing thyroid nodules and cancer in a cohort of radiation-exposed children. 764 subjects who were irradiated in the neck area as children were examined and followed for up to 25 years. All subjects underwent a clinical examination, measurements of thyrotropin, thyroglobulin levels and thyroid imaging. At baseline, 216 subjects had thyroid nodules and 548 did not. Of those with nodules, 176 underwent surgery with 55 confirmed thyroid cancers. During the follow-up, 147 subjects developed thyroid nodules including 22 with thyroid cancer. Thyroglobulin levels were higher in subjects with prevalent thyroid nodules (26.1 ng/mL vs 9.37 ng/mL; P < 0.001) and in those who had an initial normal examination but later developed thyroid nodules (11.2 ng/mL vs 8.87 ng/mL; P = 0.017). There was no relationship between baseline thyrotropin levels and the prevalent presence or absence of thyroid nodules, whether a prevalent neoplasm was benign or malignant, subsequent development of thyroid nodules during follow-up or whether an incident nodule was benign or malignant. In conclusion, in radiation-exposed children, higher thyroglobulin levels indicated an increased risk of developing thyroid nodules but did not differentiate between benign and malignant neoplasms. There was no association between the baseline TSH level and the risk of developing thyroid nodules or cancer.

Thyroid transcription factor FOXE1 interacts with ETS factor ELK1 to co-regulate TERT

BACKGROUND

Although FOXE1 was initially recognized for its role in thyroid organogenesis, more recently a strong association has been identified between the FOXE1 locus and thyroid cancer. The role of FOXE1 in adult thyroid, and in particular regarding cancer risk, has not been well established. We hypothesised that discovering key FOXE1 transcriptional partners would in turn identify regulatory pathways relevant to its role in oncogenesis.

RESULTS

In a transcription factor-binding array, ELK1 was identified to bind FOXE1. We confirmed this physical association in heterologously transfected cells by IP and mammalian two-hybrid assays. In thyroid tissue, endogenous FOXE1 was shown to bind ELK1, and using ChIP assays these factors bound thyroid-relevant gene promoters TPO and TERT in close proximity to each other. Using a combination of electromobility shift assays, TERT promoter assays and siRNA-silencing, we found that FOXE1 positively regulated TERT expression in a manner dependent upon its association with ELK1. Treating heterologously transfected thyroid cells with MEK inhibitor U0126 inhibited FOXE1-ELK1 interaction, and reduced TERT and TPO promoter activity.

METHODOLOGY

We investigated FOXE1 interactions within in vitro thyroid cell models and human thyroid tissue using a combination of immunoprecipitation (IP), chromatin IP (ChIP) and gene reporter assays.

CONCLUSIONS

FOXE1 interacts with ELK1 on thyroid relevant gene promoters, establishing a new regulatory pathway for its role in adult thyroid function. Co-regulation of TERT suggests a mechanism by which allelic variants in/near FOXE1 are associated with thyroid cancer risk.

Mismatch repair single nucleotide polymorphisms and thyroid cancer susceptibility

Thyroid cancer (TC) is the most common endocrine malignancy and its incidence continues to rise worldwide. Ionizing radiation exposure is the best established etiological factor. Heritability is high; however, despite valuable contribution from recent genome-wide association studies, the current understanding of genetic susceptibility to TC remains limited. Several studies suggest that altered function or expression of the DNA mismatch repair (MMR) system may contribute to TC pathogenesis. Therefore, the present study aimed to evaluate the potential role of a panel of MMR single nucleotide polymorphisms (SNPs) on the individual susceptibility to well-differentiated TC (DTC). A case-control study was performed involving 106 DTC patients and 212 age- and gender-matched controls, who were all Caucasian Portuguese. Six SNPs present in distinct MMR genes (MLH1 rs1799977, MSH3 rs26279, MSH4 rs5745325, PMS1 rs5742933, MLH3 rs175080 and MSH6 rs1042821) were genotyped through TaqMan^® assays and genotype-associated risk estimates were calculated. An increased risk was observed in MSH6 rs1042821 variant homozygotes [adjusted odds ratio (OR)=3.42, 95% CI: 1.04-11.24, P=0.04, under the co-dominant model; adjusted OR=3.84, 95% CI: 1.18-12.44, P=0.03, under the recessive model]. The association was especially evident for the follicular histotype and female sex. The association was also apparent when MSH6 was analysed in combination with other MMR SNPs such as MSH3 rs26279. Interestingly, two other SNP combinations, both containing the MSH6 heterozygous genotype, were associated with a risk reduction, suggesting a protective effect for these genotype combinations. These data support the idea that MMR SNPs such as MSH6 rs1042821, alone or in combination, may contribute to DTC susceptibility. This is coherent with the limited evidence available. Nevertheless, further studies are needed to validate these findings and to establish the usefulness of these SNPs as genetic susceptibility biomarkers for DTC so that, in the near future, cancer prevention policies may be optimized under a personalized medicine perspective.

Long noncoding RNAs: emerging players in thyroid cancer pathogenesis

Thyroid cancer continues to be the most common malignancy of endocrine glands. The incidence of thyroid cancer has risen significantly over the past 4 decades and has emerged as a major health issue. In recent years, significant progress has been achieved in our understanding of the molecular mechanisms of thyroid carcinogenesis, resulting in significant diagnostic, prognostic and therapeutic implications; yet, it has not reached a satisfactory level. Identifying novel molecular therapeutic targets and molecules for diagnosis and prognosis is expected to advance the overall management of this common malignancy. Long noncoding RNAs (lncRNAs) are implicated in the regulation of various key cellular genes involved in cell differentiation, proliferation, cell cycle, apoptosis, migration and invasion mainly through modulation of gene expression. Recent studies have established that lncRNAs are deregulated in thyroid cancer. In this review, we discuss extensively the tumor-suppressive (for example, LINC00271, MEG3, NAMA, PTCSC1/2/3, etc.) and oncogenic (for example, ANRIL, FAL1, H19, PVT1, etc.) roles of various lncRNAs and their possible disease associations implicated in thyroid carcinogenesis. We briefly summarize the strategies and mechanisms of lncRNA-targeting agents. We also describe the potential role of lncRNAs as prospective novel therapeutic targets, and diagnostic and prognostic markers in thyroid cancer.

rs965513 polymorphism as a common risk marker is associated with papillary thyroid cancer

Papillary thyroid cancer (PTC) is the most common type of thyroid cancer. With the rapid development of genome-wide association studies (GWAS), many genome variants associated with susceptibility to PTC have been identified, including the single nucleotide polymorphism rs965513 (9q22.33) near FOXE1. To evaluate the association between rs965513 and PTC in different ethnicities and countries, we conducted a meta-analysis using relatively large-scale samples from 23 studies (N = 163,136; 20,736 cases and 142,400 controls) by searching the PubMed and Google Scholar databases. Significant heterogeneity caused by different populations among the selected studies was observed. The A allele of rs965513 polymorphism was shown to be highly associated with risk of thyroid cancer, with odds ratios of 1.58 (95% CI 1.32–1.90) in all populations, 1.65 (95% CI 1.31–2.07)) in Caucasian populations and 1.49 in Asian populations. Compared to the dominant and recessive models, we observed the highest odds ratio (OR = 2.80, 95% CI 2.12–3.69) in the homozygous model. These results revealed that the rs965513 polymorphism is a risk factor for thyroid cancer

Lessons from Fukushima: Latest Findings of Thyroid Cancer After the Fukushima Nuclear Power Plant Accident

The increase in risk for late-onset thyroid cancer due to radiation exposure is a potential health effect after a nuclear power plant accident mainly due to the release of radioiodine in fallout. The risk is particularly elevated in those exposed during infancy and adolescence. To estimate the possibility and extent of thyroid cancer occurrence after exposure, it is of utmost importance to collect and analyze epidemiological information providing the basis for evaluation of radiation risk, and to consider radiobiology and molecular genetics. In this regard, the dose-response of cancer risk, temporal changes in the rates of thyroid cancer, its histopathological types and subtypes, and frequency of underlying genetic abnormalities are important. At present, however, it is difficult or impossible to distinguish radiation-induced thyroid cancer from spontaneous/sporadic thyroid cancer because molecular radiation signatures, biomarkers of radiation exposure, or genetic factors specific to radiation-induced cancer have not yet been identified. The large-scale ultrasound screening in Fukushima Prefecture of Japan demonstrated a high detection rate of thyroid cancer in young individuals, revealing 116 and 71 cases in the first and second rounds, respectively, among the same cohort of approximately 300,000 subjects. These findings raised concerns among residents and the public that it might be due to putative exposure to radiation from the accident at Fukushima Daiichi Nuclear Power Plant. This review summarizes evaluations by international organizations and reviews scientific publications by the authors and others on childhood thyroid cancer, especially those relevant to radiation, including basic studies on molecular mechanisms of thyroid carcinogenesis. Clinical details are also provided on surgical cases in Fukushima Prefecture, and the effect of thyroid ultrasound screening is discussed. Correct understanding of issues relating to radiation and the thyroid are essential for interpretation of thyroid cancer in Fukushima.

Current Scope and Challenges in Phenome-Wide Association Studies

PURPOSE OF REVIEW

Over many decades, researchers have been designing studies to investigate the relationship between genotypes and phenotypes to gain an understanding about the effect of genetics on disease. Recently, a high-throughput approach called phenome-wide associations studies (PheWAS) have been extensively used to identify associations between genetic variants and many diseases and traits simultaneously. In this review, we describe the value of PheWAS along with methodological issues and challenges in interpretation for current applications of PheWAS.

RECENT FINDINGS

PheWAS have uncovered a paradigm to identify new associations for genetic loci across many diseases. The application of PheWAS have been effective with phenotype data from electronic health records, epidemiological studies, and clinical trials data.

SUMMARY

The key strength of a PheWAS is to identify the association of one or more genetic variants with multiple phenotypes, which can showcase interconnections among the phenotypes due to shared genetic associations. While the PheWAS approach appears promising, there are a number of challenges that need to be addressed to provide additional robustness to PheWAS findings.

Genome-wide association and expression quantitative trait loci studies identify multiple susceptibility loci for thyroid cancer

Thyroid cancer is the most common cancer in Korea. Several susceptibility loci of differentiated thyroid cancer (DTC) were identified by previous genome-wide association studies (GWASs) in Europeans only. Here we conducted a GWAS and a replication study in Koreans using a total of 1,085 DTC cases and 8,884 controls, and validated these results using expression quantitative trait loci (eQTL) analysis and clinical phenotypes. The most robust associations were observed in the NRG1 gene (rs6996585, P=1.08 × 10-10) and this SNP was also associated with NRG1 expression in thyroid tissues. In addition, we confirmed three previously reported loci (FOXE1, NKX2-1 and DIRC3) and identified seven novel susceptibility loci (VAV3, PCNXL2, INSR, MRSB3, FHIT, SEPT11 and SLC24A6) associated with DTC. Furthermore, we identified specific variants of DTC that have different effects according to cancer type or ethnicity. Our findings provide deeper insight into the genetic contribution to thyroid cancer in different populations.

Investigation of DNA repair-related SNPs underlying susceptibility to papillary thyroid carcinoma reveals MGMT as a novel candidate gene in Belarusian children exposed to radiation

BACKGROUND

Genetic factors may influence an individual's sensitivity to ionising radiation and therefore modify his/her risk of developing papillary thyroid carcinoma (PTC). Previously, we reported that common single nucleotide polymorphisms (SNPs) within the DNA damage recognition gene ATM contribute to PTC risk in Belarusian children exposed to fallout from the Chernobyl power plant accident. Here we explored in the same population the contribution of a panel of DNA repair-related SNPs in genes acting downstream of ATM.

METHODS

The association of 141 SNPs located in 43 DNA repair genes was examined in 75 PTC cases and 254 controls from the Gomel region in Belarus. All subjects were younger than 15 years at the time of the Chernobyl accident. Conditional logistic regressions accounting for radiation dose were performed with PLINK using the additive allelic inheritance model, and a linkage disequilibrium (LD)-based Bonferroni correction was used for correction for multiple testing.

RESULTS

The intronic SNP rs2296675 in MGMT was associated with an increased PTC risk [per minor allele odds ratio (OR) 2.54 95% CI 1.50, 4.30, P per allele = 0.0006, P corr.= 0.05], and gene-wide association testing highlighted a possible role for ERCC5 (P Gene = 0.01) and PCNA (P Gene = 0.05) in addition to MGMT (P Gene = 0.008).

CONCLUSIONS

These findings indicate that several genes acting in distinct DNA repair mechanisms contribute to PTC risk. Further investigation is needed to decipher the functional properties of the methyltransferase encoded by MGMT and to understand how alteration of such functions may lead to the development of the most common type of thyroid cancer.

Inherited variants in genes somatically mutated in thyroid cancer

Background

Tumour suppressor genes when mutated in the germline cause various cancers, but they can also be somatically mutated in sporadic tumours. We hypothesized that there may also be cancer-related germline variants in the genes commonly mutated in sporadic well-differentiated thyroid cancer (WDTC).

Methods

We performed a two-stage case-control association study with a total of 2214 cases and 2108 healthy controls from an Italian population. By genotyping 34 single nucleotide polymorphisms (SNPs), we covered a total of 59 missense SNPs and SNPs located in the 5' and 3' untranslated regions (UTRs) of 10 different genes.

Results

The Italian1 series showed a suggestive association for 8 SNPs, from which three were replicated in the Italian2 series. The meta-analysis revealed a study-wide significant association for rs459552 (OR: 0.84, 95%CI: 0.75–0.94) and rs1800900 (OR: 1.15, 95%CI: 1.05–1.27), located in the APC and GNAS genes, respectively. The APC rs459552 is a missense SNP, located in a conserved amino acid position, but without any functional consequences. The GNAS rs1800900 is located at a conserved 5'UTR and according to the experimental ENCODE data it may affect promoter and histone marks in different cell types.

Conclusions

The results of this study yield new insights on WDTC, showing that inherited variants in the APC and GNAS genes can play a role in the etiology of thyroid cancer. Further studies are necessary to better understand the role of the identified SNPs in the development of WDTC and to functionally validate our in silico predictions.

A genome-wide association study yields five novel thyroid cancer risk loci

The great majority of thyroid cancers are of the non-medullary type. Here we report findings from a genome-wide association study of non-medullary thyroid cancer, including in total 3,001 patients and 287,550 controls from five study groups of European descent. Our results yield five novel loci (all with Pcombined<3 × 10-8): 1q42.2 (rs12129938 in PCNXL2), 3q26.2 (rs6793295 a missense mutation in LRCC34 near TERC), 5q22.1 (rs73227498 between NREP and EPB41L4A), 10q24.33 (rs7902587 near OBFC1), and two independently associated variants at 15q22.33 (rs2289261 and rs56062135; both in SMAD3). We also confirm recently published association results from a Chinese study of a variant on 5p15.33 (rs2736100 near the TERT gene) and present a stronger association result for a moderately correlated variant (rs10069690; OR=1.20, P=3.2 × 10-7) based on our study of individuals of European ancestry. In combination, these results raise several opportunities for future studies of the pathogenesis of thyroid cancer.

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

BACKGROUND

Several functional single-nucleotide polymorphisms (SNPs) at the FOXE1 locus on chromosome 9q22.33 have been associated with the risk for papillary thyroid carcinoma (PTC). This study set out to elucidate whether their effects are independent, using genotyping results in populations of Asian and European descent.

METHODS

SNPs rs965513 and rs1867277 and a polymorphic region determining the length of the FOXE1 polyalanine (poly-Ala) tract were genotyped in 501 patients with PTC and 748 healthy individuals from Japan, and in 660 patients and 820 population controls from Belarus. Functional analysis of transactivation activities of FOXE1 isoforms with varying number of alanine repeats was performed by a Dual-Luciferase^® Assay.

RESULTS

All three polymorphisms were significantly associated with PTC in both populations on univariate analysis. However, conditional analysis revealed independent effects of rs965513 and rs1867277 SNPs but not of the FOXE1 poly-Ala polymorphism. The independent effect of the lead rs965513 SNP was observed in both populations, while that of rs1867277 was only identified in the Japanese population, in which linkage disequilibrium between the three polymorphisms is markedly weaker. Despite the strong decrease in transcriptional activity with increasing FOXE1 poly-Ala tract length, no difference in transactivation potential of the FOXE1 poly-Ala isoforms could be seen after adjustment for the minimal promoter activity in the reporter vectors. Plasmids encoding FOXE1 isoforms of increasing poly-Ala tract length were also found to produce less FOXE1 protein after cell transfection.

CONCLUSIONS

The functional variants rs965513 and rs1867277 independently contribute to genetic predisposition to PTC, while a contributing role of the FOXE1 poly-Ala polymorphism could not be confirmed.