Discovery of common variants associated with low TSH levels and thyroid cancer risk

Assessing thyroid cancer risk using polygenic risk scores

Genome-wide association studies (GWASs) have identified at least 10 single-nucleotide polymorphisms (SNPs) associated with papillary thyroid cancer (PTC) risk. Most of these SNPs are common variants with small to moderate effect sizes. Here we assessed the combined genetic effects of these variants on PTC risk by using summarized GWAS results to build polygenic risk score (PRS) models in three PTC study groups from Ohio (1,544 patients and 1,593 controls), Iceland (723 patients and 129,556 controls), and the United Kingdom (534 patients and 407,945 controls). A PRS based on the 10 established PTC SNPs showed a stronger predictive power compared with the clinical factors model, with a minimum increase of area under the receiver-operating curve of 5.4 percentage points (P ≤ 1.0 × 10^-9). Adding an extended PRS based on 592,475 common variants did not significantly improve the prediction power compared with the 10-SNP model, suggesting that most of the remaining undiscovered genetic risk in thyroid cancer is due to rare, moderate- to high-penetrance variants rather than to common low-penetrance variants. Based on the 10-SNP PRS, individuals in the top decile group of PRSs have a close to sevenfold greater risk (95% CI, 5.4-8.8) compared with the bottom decile group. In conclusion, PRSs based on a small number of common germline variants emphasize the importance of heritable low-penetrance markers in PTC.

Differentiated thyroid cancer and Hashimoto thyroiditis: Utility of the Afirma gene expression classifier

BACKGROUND AND OBJECTIVES

The Afirma gene expression classifier (AGEC) has not been tested or validated in a high-risk group, such as patients with Hashimoto's thyroiditis (HT). We hypothesized that AGEC would perform worse in patients with HT.

METHODS

A retrospective review of patient charts in a single academic institution who underwent thyroidectomy between 2012 and 2017 was conducted. Patients with HT who underwent AGEC were identified to calculate sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV).

RESULTS

We identified 69 patients with HT and atypia of undetermined significance/follicular lesion of undetermined significance (AUS/FLUS) on cytology who underwent AGEC analysis. The mean age of AGEC cohort was 50 years (range, 26-77 years) with 90% female. The median nodule size was 1.9 cm (interquartile range [IQR], 1.2-2.7 cm). Of the 69 patients, 62 showed suspicious AGEC of which 26 showed TC on surgical pathology. Of the seven benign AGEC, two showed TC on surgical pathology. The sensitivity, specificity, PPV, and NPV were 93%, 12%, 42%, and 71%, respectively. Of the entire AGEC cohort, 17 (43%) showed multicentric disease.

CONCLUSIONS

We observed a lower NPV for AGEC to rule out thyroid cancer in patients with HT, which reduces the utility of the test for this population.

FOXE1 regulates migration and invasion in thyroid cancer cells and targets ZEB1

FOXE1 is a thyroid-specific transcription factor essential for thyroid gland development and maintenance of the differentiated state. Interestingly, a strong association has been recently described between FOXE1 expression and susceptibility to thyroid cancer, but little is known about the mechanisms underlying FOXE1-induced thyroid tumorigenesis. Here, we used a panel of human thyroid cancer-derived cell lines covering the spectrum of thyroid cancer phenotypes to examine FOXE1 expression and to test for correlations between FOXE1 expression, the allele frequency of two SNPs and a length polymorphism in or near the FOXE1 locus associated with cancer susceptibility, and the migration ability of thyroid cancer cell lines. Results showed that FOXE1 expression correlated with differentiation status according to histological sub-type, but not with SNP genotype or cell migration ability. However, loss-and-gain-of-function experiments revealed that FOXE1 modulates cell migration, suggesting a role in epithelial-to-mesenchymal transition (EMT). Our previous genome-wide expression analysis identified Zeb1, a major EMT inducer, as a putative Foxe1 target gene. Indeed, gene silencing of FOXE1 decreased ZEB1 expression, whereas its overexpression increased ZEB1 transcriptional activity. FOXE1 was found to directly interact with the ZEB1 promoter. Lastly, ZEB1 silencing decreased the ability of thyroid tumoral cells to migrate and invade, pointing to its importance in thyroid tumor mestastases. In conclusion, we have identified ZEB1 as a bona fide target of FOXE1 in thyroid cancer cells, which provides new insights into the role of FOXE1 in regulating cell migration and invasion in thyroid cancer.

Variants in LRRC34 reveal distinct mechanisms for predisposition to papillary thyroid carcinoma

BACKGROUND

Papillary thyroid carcinoma (PTC) demonstrates high heritability and a low somatic mutation burden relative to other cancers. Therefore, the genetic risk predisposing to PTC is likely due to a combination of low penetrance variants. A recent genome-wide association study revealed the association of PTC with a missense variant, rs6793295, at 3q26 in a gene called Leucine Repeat Rich Containing 34 (LRRC34).

METHODS

We report the mechanisms of PTC risk at 3q26 using a combination of overexpression, mass spectroscopy, knockdown, transcriptome profiling, migration assays and genetic analysis.

RESULTS

We observed differential binding of wild-type and missense LRRC34 to RANBP1. Overexpression of missense LRRC34 reduced RanGTP levels and increased apoptosis. We also identified a second linkage disequilibrium (LD) block upstream of LRRC34 containing regulatory variants with allele-specific expression. Transcriptome profiling of LRRC34 knockdown cells showed changes in genes involved with cellular movement. LRRC34 knockdown reduced the migration of thyroid cancer cell lines. Lastly, we assessed the relative contribution of PTC risk from each locus using haplotype analysis.

CONCLUSIONS

Our study demonstrates two separate mechanisms, one in G protein signalling and the other in transcriptional control, dictating PTC risk at 3q26 using both biochemical and genetic techniques.

Does the TT Variant of the rs966423 Polymorphism in DIRC3 Affect the Stage and Clinical Course of Papillary Thyroid Cancer?

Thyroid cancer (TC) is the most common cancer of the endocrine system. Most new diagnoses are of low-grade papillary thyroid cancer (PTC), suggesting that PTC may be over-diagnosed. However, the incidence of advanced-stage PTC has increased in recent years. It is therefore very important to identify prognostic factors for advanced PTC. Somatic mutation of the BRAF gene at V600E, or the coexistence of the BRAF V600E mutation and mutations in the TERT promoter are associated with more aggressive disease. It would also be valuable to identify genetic risk factors affecting PTC prognosis. We therefore evaluated the impact of the rs966423 polymorphism in the DIRC3 gene, including its relationship with unfavorable histopathological and clinical features and mortality, in differentiated thyroid cancer (DTC). The study included 1466 patients diagnosed with DTC from one center. There was no significant association between the DIRC3 genotype at rs966423 (CC, CT, or TT) and any histopathological or clinic factor examined, including initial response to therapy, response at follow-up, or overall mortality, in DTC patients.

Association of Thyroid Function Genetic Predictors With Atrial Fibrillation: A Phenome-Wide Association Study and Inverse-Variance Weighted Average Meta-analysis

Importance

Thyroid hormone levels are tightly regulated through feedback inhibition by thyrotropin, produced by the pituitary gland. Hyperthyroidism is overwhelmingly due to thyroid disorders and is well recognized to contribute to a wide spectrum of cardiovascular morbidity, particularly the increasingly common arrhythmia atrial fibrillation (AF).

Objective

To determine the association between genetically determined thyrotropin levels and AF.

Design, Setting, and Participants

This phenome-wide association study scanned 1318 phenotypes associated with a polygenic predictor of thyrotropin levels identified by a previously published genome-wide association study that included participants of European ancestry. North American individuals of European ancestry with longitudinal electronic health records were analyzed from May 2008 to November 2016. Analysis began March 2018.

Main Outcomes and Measures

Clinical diagnoses associated with a polygenic predictor of thyrotropin levels.

Exposures

Genetically determined thyrotropin levels.

Results

Of 37 154 individuals, 19 330 (52%) were men. The thyrotropin polygenic predictor was positively associated with hypothyroidism (odds ratio [OR], 1.10; 95% CI, 1.07-1.14; P = 5 × 10-11) and inversely associated with diagnoses related to hyperthyroidism (OR, 0.64; 95% CI, 0.54-0.74; P = 2 × 10-8 for toxic multinodular goiter). Among nonthyroid associations, the top association was AF/flutter (OR, 0.93; 95% CI, 0.9-0.95; P = 9 × 10-7). When the analyses were repeated excluding 9801 individuals with any diagnoses of a thyroid-related disease, the AF association persisted (OR, 0.91; 95% CI, 0.88-0.95; P = 2.9 × 10-6). To replicate this association, we conducted an inverse-variance weighted average meta-analysis using AF single-nucleotide variant weights from a genome-wide association study of 17 931 AF cases and 115 142 controls. As in the discovery analyses, each SD increase in predicted thyrotropin was associated with a decreased risk of AF (OR, 0.86; 95% CI, 0.79-0.93; P = 4.7 × 10-4). In a set of AF cases (n = 745) and controls (n = 1680) older than 55 years, directly measured thyrotropin levels that fell within the normal range were inversely associated with AF risk (OR, 0.91; 95% CI, 0.83-0.99; P = .04).

Conclusions and Relevance

This study suggests a role for genetically determined variation in thyroid function within a physiologically accepted normal range as a risk factor for AF. The clinical decision to treat subclinical thyroid disease should incorporate the risk for AF as antithyroid medications to treat hyperthyroidism may reduce AF risk and thyroid hormone replacement for hypothyroidism may increase AF risk.

Assessment of the Relationship Between Genetic Determinants of Thyroid Function and Atrial Fibrillation: A Mendelian Randomization Study

Importance

Increased free thyroxine (FT4) and decreased thyrotropin are associated with increased risk of atrial fibrillation (AF) in observational studies, but direct involvement is unclear.

Objective

To evaluate the potential direct involvement of thyroid traits on AF.

Design, Setting, and Participants

Study-level mendelian randomization (MR) included 11 studies, and summary-level MR included 55 114 AF cases and 482 295 referents, all of European ancestry.

Exposures

Genomewide significant variants were used as instruments for standardized FT4 and thyrotropin levels within the reference range, standardized triiodothyronine (FT3):FT4 ratio, hypothyroidism, standardized thyroid peroxidase antibody levels, and hyperthyroidism. Mendelian randomization used genetic risk scores in study-level analysis or individual single-nucleotide polymorphisms in 2-sample MR for the summary-level data.

Main Outcomes and Measures

Prevalent and incident AF.

Results

The study-level analysis included 7679 individuals with AF and 49 233 referents (mean age [standard error], 62 [3] years; 15 859 men [29.7%]). In study-level random-effects meta-analysis, the pooled hazard ratio of FT4 levels (nanograms per deciliter) for incident AF was 1.55 (95% CI, 1.09-2.20; P = .02; I2 = 76%) and the pooled odds ratio (OR) for prevalent AF was 2.80 (95% CI, 1.41-5.54; P = .003; I2 = 64%) in multivariable-adjusted analyses. The FT4 genetic risk score was associated with an increase in FT4 by 0.082 SD (standard error, 0.007; P < .001) but not with incident AF (risk ratio, 0.84; 95% CI, 0.62-1.14; P = .27) or prevalent AF (OR, 1.32; 95% CI, 0.64-2.73; P = .46). Similarly, in summary-level inverse-variance weighted random-effects MR, gene-based FT4 within the reference range was not associated with AF (OR, 1.01; 95% CI, 0.89-1.14; P = .88). However, gene-based increased FT3:FT4 ratio, increased thyrotropin within the reference range, and hypothyroidism were associated with AF with inverse-variance weighted random-effects OR of 1.33 (95% CI, 1.08-1.63; P = .006), 0.88 (95% CI, 0.84-0.92; P < .001), and 0.94 (95% CI, 0.90-0.99; P = .009), respectively, and robust to tests of horizontal pleiotropy. However, the subset of hypothyroidism single-nucleotide polymorphisms involved in autoimmunity and thyroid peroxidase antibodies levels were not associated with AF. Gene-based hyperthyroidism was associated with AF with MR-Egger OR of 1.31 (95% CI, 1.05-1.63; P = .02) with evidence of horizontal pleiotropy (P = .045).

Conclusions and Relevance

Genetically increased FT3:FT4 ratio and hyperthyroidism, but not FT4 within the reference range, were associated with increased AF, and increased thyrotropin within the reference range and hypothyroidism were associated with decreased AF, supporting a pathway involving the pituitary-thyroid-cardiac axis.

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.

Phytoestrogens and Thyroid Cancer Risk: A Population-Based Case-Control Study in Connecticut

BACKGROUND

Very few previous studies have examined the relationship between thyroid cancer risk and intake of phytoestrogens (PE); furthermore, these studies have reached inconsistent results.

METHODS

We analyzed data from a population-based case-control study in Connecticut from 2010 to 2011, including 387 histologically confirmed thyroid cancer cases and 433 population-based controls, with compound data available concerning specific PEs. Multivariate unconditional logistic regression models were used to estimate the associations between specific PEs and the risk of thyroid cancer, adjusting for potential confounders.

RESULTS

An elevated risk of thyroid cancer was associated with moderate to high levels of coumestrol intake [OR = 2.48, 95% confidence interval (CI), 1.39-4.43 for 40-80 μg/day; OR = 2.41, 95% CI, 1.32-4.40 for 80-130 μg/day; and OR = 2.38, 95% CI, 1.26-4.50 for >200 μg/day compared with <40 μg/day], and the main elevation in risk appeared among microcarcinomas (≤1 cm). A decreased risk of papillary macrocarcinomas (>1 cm; OR = 0.26, 95% CI, 0.08-0.85 for 1,860-3,110 μg/day compared with <760 μg/day) was associated with moderate genistein intake among women.

CONCLUSIONS

Our study suggests that high coumestrol intake increases the risk of thyroid cancer, especially microcarcinomas, whereas moderate amounts of genistein intake appear to be protective for females with thyroid macrocarcinomas.

IMPACT

The study highlights the importance of distinguishing between microcarcinomas and macrocarcinomas in future research on the etiology of thyroid cancer.

Genetic Variants Implicate Dual Oxidase-2 in Familial and Sporadic Nonmedullary Thyroid Cancer

Highly penetrant hereditary thyroid cancer manifests as familial nonmedullary thyroid cancer (FNMTC), whereas low-penetrance hereditary thyroid cancer manifests as sporadic disease and is associated with common polymorphisms, including rs965513[A]. Whole-exome sequencing of an FNMTC kindred identified a novel Y1203H germline dual oxidase-2 (DUOX2) mutation. DUOX2_Y1203H is enzymatically active, with increased production of reactive oxygen species. Furthermore, patients with sporadic thyroid cancer homozygous for rs965513[A] demonstrated higher DUOX2 expression than heterozygous rs965513[A/G] or homozygous rs965513[A]-negative patients. These data suggest that dysregulated hydrogen peroxide metabolism is a common mechanism by which high- and low-penetrance genetic factors increase thyroid cancer risk. SIGNIFICANCE: This study provides novel insights into the genetic and molecular mechanisms underlying familial and sporadic thyroid cancers.

High Rates of Underlying Thyroid Cancer in Patients Undergoing Thyroidectomy for Hyperthyroidism

BACKGROUND

The rate of thyroid cancer in patients with hyperthyroidism is reported to be rare, and patients with toxic thyroid nodules do not routinely undergo fine-needle aspiration (FNA) to evaluate for malignancy. However, higher rates of malignancy in hyperthyroid patients may exist than previously reported. This study examines the rate of malignancy in patients with hyperthyroidism who have undergone thyroidectomy.

METHODS

A retrospective review of prospectively collected data of 138 patients with hyperthyroidism who underwent thyroidectomy at a single institution was performed. Patients were divided into three groups: Graves' disease (n = 80), toxic multinodular goiter (n = 46), and toxic solitary nodule (n = 12). Patients with previous thyroid surgery were excluded from the study. All patients had biochemical confirmation of hyperthyroidism with thyroid-stimulating hormone <0.1 mIU/L and clinical diagnosis by a referring physician.

RESULTS

Of 138 patients, 22% (31/138) were found to have malignancy on final pathology. The breakdown of malignancy by hyperthyroid condition was as follows: 16% in Graves' disease, 24% in toxic multinodular goiter patients, and 50% in toxic solitary nodule patients.

CONCLUSIONS

There is a clinically significant rate of malignancy seen in patients who undergo thyroidectomy for hyperthyroidism. Patients with distinct thyroid nodules in the presence of hyperthyroidism may have the highest rates of malignancy and should undergo appropriate workup with ultrasound and FNA to exclude underlying malignancy. In cases with suspicious ultrasound features and/or FNA cytopathology, surgical treatment should be considered as initial management.

FOXE1 inhibits cell proliferation, migration and invasion of papillary thyroid cancer by regulating PDGFA

PURPOSE

Forkhead box E1 (FOXE1) plays an important role in the development, proliferation and differentiation of thyroid cells. However, the biological functions of FOXE1 in papillary thyroid cancer (PTC) remain unclear.

MATERIALS AND METHODS

In this study, the level of FOXE1 expression was examined in human PTC tissues and cells. Then, the high expression of FOXE1 was specifically silenced by RNA interference in vitro. Subsequently, FOXE1-shRNA was transfected into PTC cells (TPC-1 and K1). The effects on cell proliferation, migration and invasion were evaluated. In addition, FOXE1 targets were screened by cDNA microarray assays. The correlation between the expression of target gene platelet-derived growth factor A (PDGFA) and clinicopathological features of PTC patients was analysed.

RESULTS

FOXE1 is highly expressed in PTC tissues and PTC cell lines. The silencing of FOXE1 significantly promotes PTC cell proliferation, migration and invasion in vitro. The cDNA microarray analyses show that PDGFA is a critical downstream target gene of FOXE1 in PTC cells. It was also observed that PDGFA is negatively regulated by FOXE1 in PTC. The clinical data indicate that the low expression level of PDGFA is correlated with the small size of PTC.

CONCLUSION

Collectively, the results indicate for the first time that high expression of FOXE1 may function as a tumour suppressor in the early stage of PTC and restrain the proliferation, migration and invasion of PTC by negatively regulating PDGFA expression. Thus, FOXE1 could serve as a prognostic biomarker for PTC.

Identification of Rare Variants Predisposing to Thyroid Cancer

Background: Familial non-medullary thyroid cancer (NMTC) accounts for a relatively small proportion of thyroid cancer cases, but it displays strong genetic predisposition. So far, only a few NMTC susceptible genes and low-penetrance variants contributing to NMTC have been described. This study aimed to identify rare germline variants that may predispose individuals to NMTC by sequencing a cohort of 17 NMTC families. Methods: Whole-genome sequencing and genome-wide linkage analysis were performed in 17 NMTC families. MendelScan and BasePlayer were applied to screen germline variants followed by customized filtering. The remaining candidate variants were subsequently validated by Sanger sequencing. A panel of 277 known cancer predisposition genes was also screened in these families. Results: A total of 41 rare coding candidate variants in 40 genes identified by whole-genome sequencing are reported, including 24 missense, five frameshift, five splice change, and seven nonsense variants. Sanger sequencing confirmed all 41 rare variants and proved their co-segregation with NMTC in the extended pedigrees. In silico functional analysis of the candidate genes using Ingenuity Pathway Analysis showed that cancer was the top category of "Diseases and Disorders." Additionally, a targeted search displayed six variants in known cancer predisposition genes, including one frameshift variant and five missense variants. Conclusions: The data identify rare germline variants that may play important roles in NMTC predisposition. It is proposed that in future research including functional characterization, these variants and genes be considered primary candidates for thyroid cancer predisposition.

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.

Systematic evaluation of cancer-specific genetic risk score for 11 types of cancer in The Cancer Genome Atlas and Electronic Medical Records and Genomics cohorts

BACKGROUND

Genetic risk score (GRS) is an odds ratio (OR)-weighted and population-standardized method for measuring cumulative effect of multiple risk-associated single nucleotide polymorphisms (SNPs). We hypothesize that GRS is a valid tool for risk assessment of most common cancers.

METHODS

Utilizing genotype and phenotype data from The Cancer Genome Atlas (TCGA) and Electronic Medical Records and Genomics (eMERGE), we tested 11 cancer-specific GRSs (bladder, breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, prostate, renal, and thyroid cancer) for association with the respective cancer type. Cancer-specific GRSs were calculated, for the first time in these cohorts, based on previously published risk-associated SNPs using the Caucasian subjects in these two cohorts.

RESULTS

Mean cancer-specific GRS in the population controls of eMERGE approximated the expected value of 1.00 (between 0.98 and 1.02) for all 11 types of cancer. Mean cancer-specific GRS was consistently higher in respective cancer patients than controls for all 11 types of cancer (P < 0.05). When subjects were categorized into low-, average-, and high-risk groups based on cancer-specific GRS (<0.5, 0.5-1.5, and >1.5, respectively), significant dose-response associations of higher cancer-specific GRS with higher OR of respective type of cancer were found for nine types of cancer (P_-trend  < 0.05). More than 64% subjects in the population controls of eMERGE can be classified as high risk for at least one type of these cancers.

CONCLUSION

Validity of GRS for predicting cancer risk is demonstrated for most types of cancer. If confirmed in larger studies, cancer-specific GRS may have the potential for developing personalized cancer screening strategy.

Fasting serum glucose, thyroid-stimulating hormone, and thyroid hormones and risk of papillary thyroid cancer: A case-control study

BACKGROUND

This study was to investigate the association of fasting serum glucose (FSG), thyroid-stimulating hormone (TSH), and thyroid hormones with papillary thyroid cancer (PTC).

METHODS

A total of 649 participants were included in this case-control study. The associations of FSG, TSH, free triiodothyronine (FT3) and free thyroxine (FT4) with PTC were estimated using an unconditional logistic regression.

RESULTS

Compared with the lowest quintile of TSH levels, odds ratios (ORs) and 95% confidence intervals (CIs) for association between PTC risk and highest quintile of TSH levels were 1.67 (95% CI, 0.99-2.83). However, this risk correlation was more significant in PTC cases with ≤1.0 cm tumor size (adjusted OR, 1.95; 95% CI, 1.08-3.54; adjusted P-trend, 0.05). The PTC risk was also inversely associated with the serum FT3 level in all participants (adjusted P-trend, 0.001), but positively associated with the serum FT4 (adjusted P-trend, 0.001) and FSG (adjusted P-trend, 0.01) levels. Among the participants without diabetes, the individuals with high FSG levels and abnormal TSH concentration had an increased PTC risk (adjusted OR, 3.38; 95% CI, 1.78-6.42).

CONCLUSION

The current study provides evidence for the association of FSG, TSH, and thyroid hormones (FT3 and FT4) with PTC risk. However, larger relative studies are needed.

Association between hTERT Polymorphisms and Female Papillary Thyroid Carcinoma

BACKGROUND

A previous genome-wide association study showed that hTERT rs10069690 and rs2736100 polymorphisms were associated with thyroid cancer risk.

OBJECTIVE

This study further investigated the association between increased risk and clinicopathologic characteristics for Papillary Thyroid Carcinoma (PTC) and hTERT polymorphisms rs10069690 or rs2736100 in a Chinese female population.

METHODS

The hTERT genotypes of 276 PTC patients and 345 healthy subjects were determined with regard to SNPs rs10069690 and rs2736100. The association between these SNPs and the risk of PTC and clinicopathologic characteristics was investigated by logistic regression.

RESULTS

We found a significant difference between PTC and rs10069690 (Odds Ratio (OR) = 1.515; P = 0.005), but not between PTC and rs2736100. When the analysis was limited to females, rs10069690 and rs2736100 were both associated with increased risk for PTC in female individuals (OR = 1.647, P = 0.007; OR = 1.339, P = 0.041, respectively). Further haplotype analysis revealed a stimulative effect of haplotypes TC and CA of TERT rs10069690-rs2736100, which increased risk for PTC in female individuals (OR = 1.579, P = 0.014; OR = 0.726, P = 0.025, respectively). Furthermore, the heterozygote A/C of rs2736100 showed significant difference for age (OR = 0.514, P = 0.047).

CONCLUSION

Our finding suggests that hTERT polymorphisms rs10069690 and rs2736100 are associated with increased risk for PTC in Chinese female population and rs2736100 may be related to age. Consistent with US20170360914 and US20170232075, they are expected to be a potential molecular target for anti-cancer therapy.

Fine mapping of 14q13 reveals novel variants associated with different histological subtypes of papillary thyroid carcinoma

The first two genome wide association studies (GWAS) of papillary thyroid carcinoma (PTC) detected five variants associated with PTC. Two of them (rs944289 and rs116909374) are located at 14q13 making that locus an important target of research into the genetic predisposition to PTC. We aimed at uncovering other variants at 14q13 associated with PTC independently from the GWAS variants. We performed next generation sequencing of the 14q13 region and analyzed the allele frequencies of single nucleotide polymorphisms (SNPs) in n = 90 PTC cases vs. n = 379 EUR controls from the 1,000 Genome Project. The variants associated with PTC were validated in an Ohio cohort of n = 1,216 PTC cases and n = 1,416 controls. Next, we analyzed the association between SNPs and expression of nearby genes and clinical parameters. We showed that rs368187 was associated with PTC (OR = 1.31, p = 2.20 × 10^-6 ). Rs1632250, Rs1863347 and rs1755787 showed association with classical PTC (cPTC) (n = 891; OR = 1.24, 2.22 × 10^-3 , OR = 1.31, p = 2.15 × 10^-4 and OR = 1.24, p = 2.06 × 10^-3 , respectively) while variant rs28397092 showed association with follicular variant (n = 243; OR = 1.51, p = 1.36 × 10^-3 ). Rs1863347 was associated with suppression of PTCSC3 in unaffected thyroid tissue (p = 0.026). Rs1632250, rs1863347 and rs1755787 showed association with multifocality (OR = 1.85, p = 0.001, OR = 1.98, p = 0.001 and OR = 1.76, p = 0.003 respectively) and N stage (OR = 1.79, p = 0.014, OR = 1.73, p = 0.023 and OR = 1.81, p = 0.013, respectively) in microPTC (n = 328) while rs368187 was associated with M stage (OR = 0.56, p = 0.034) in cPTC. Our results disclose multiple variants associated with PTC and clinical features in the 14q13 superlocus. We suggest that translational genotype/phenotype studies should take into account not only somatic mutations but also germline variants.

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.

Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves’ disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets.

Thyroid dysfunction is a common public health problem and associated with cardiovascular co-morbidities. Here, the authors carry out genome-wide meta-analysis for thyroid hormone (TH) levels, hyper- and hypothyroidism and identify SLC17A4 as a TH transporter and AADAT as a TH metabolizing enzyme.

Thyroid cancers of follicular origin in a genomic light: in-depth overview of common and unique molecular marker candidates

In recent years, thyroid malignances have become more prevalent, especially among women. The most common sporadic types of thyroid tumors of follicular origin include papillary, follicular and anaplastic thyroid carcinomas. Although modern diagnosis methods enable the identification of tumors of small diameter, tumor subtype differentiation, which is imperative for the correct choice of treatment, is still troublesome. This review discusses the recent advances in the field of molecular marker identification via next-generation sequencing and microarrays. The potential use of these biomarkers to distinguish among the most commonly occurring sporadic thyroid cancers is presented and compared. Geographical heterogeneity might be a differentiator, although not necessarily a limiting factor, in biomarker selection. The available data advocate for a subset of mutations common for the three subtypes as well as mutations that are unique for a particular tumor subtype. Tumor heterogeneity, a known issue occurring within solid malignancies, is also discussed where applicable. Public databases with datasets derived from high-throughput experiments are a valuable source of information that aid biomarker research in general, including the identification of molecular hallmarks of thyroid cancer.

Association of Polygenic Risk Scores for Multiple Cancers in a Phenome-wide Study: Results from The Michigan Genomics Initiative

Health systems are stewards of patient electronic health record (EHR) data with extraordinarily rich depth and breadth, reflecting thousands of diagnoses and exposures. Measures of genomic variation integrated with EHRs offer a potential strategy to accurately stratify patients for risk profiling and discover new relationships between diagnoses and genomes. The objective of this study was to evaluate whether polygenic risk scores (PRS) for common cancers are associated with multiple phenotypes in a phenome-wide association study (PheWAS) conducted in 28,260 unrelated, genotyped patients of recent European ancestry who consented to participate in the Michigan Genomics Initiative, a longitudinal biorepository effort within Michigan Medicine. PRS for 12 cancer traits were calculated using summary statistics from the NHGRI-EBI catalog. A total of 1,711 synthetic case-control studies was used for PheWAS analyses. There were 13,490 (47.7%) patients with at least one cancer diagnosis in this study sample. PRS exhibited strong association for several cancer traits they were designed for, including female breast cancer, prostate cancer, melanoma, basal cell carcinoma, squamous cell carcinoma, and thyroid cancer. Phenome-wide significant associations were observed between PRS and many non-cancer diagnoses. To differentiate PRS associations driven by the primary trait from associations arising through shared genetic risk profiles, the idea of "exclusion PRS PheWAS" was introduced. Further analysis of temporal order of the diagnoses improved our understanding of these secondary associations. This comprehensive PheWAS used PRS instead of a single variant.

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.

Association of Established Thyroid-stimulating Hormone and Free Thyroxine Genetic Variants with Hashimoto's Thyroiditis

Hashimoto's thyroiditis (HT), the most frequent autoimmune thyroid disease (AITD), is characterized by chronic inflammation of the thyroid gland that usually results in hypothyroidism. Thyroid-stimulating hormone (TSH) and free thyroxine (FT4) levels are used as clinical determinants of thyroid function. The main aim of this study was to explore the association of established TSH and FT4 genetic variants with HT. We performed a case-control analysis using 23 genetic markers in 200 HT patients and 304 controls. Additionally, we tested the association of selected variants with several thyroid-related quantitative traits in HT cases only. Two genetic variants showed nominal association with HT: rs11935941 near NR3C2 gene (p = 0.0034, OR = 0.57, 95% CI = 0.39-0.83) and rs1537424 near MBIP gene (p = 0.0169, OR = 0.72, 95% CI = 0.55-0.94). Additionally, three SNPs showed nominal association with thyroglobulin antibody (TgAb) levels: rs4804416 in INSR gene (p = 0.0073, β = -0.51), rs6435953 near IGFBP5 gene (p = 0.0081, β = 0.75), and rs1537424 near MBIP gene (p = 0.0117, β = 0.49). GLIS3 genetic variant rs10974423 showed nominal association with thyroid peroxidase antibody (TPOAb) levels (p = 0.0465, β = -0.56) and NRG1 genetic variant rs7825175 was nominally associated with thyroid gland volume (p = 0.0272, β = -0.18). All detected loci were previously related to thyroid function or pathology. Findings from our study suggest biological relevance of NR3C2 and MBIP with HT, although these loci require additional confirmation in a larger replication study.

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.

The Role of NRG1 in the Predisposition to Papillary Thyroid Carcinoma

Context

Previous genome-wide association studies have shown that single-nucleotide polymorphism (SNP) rs2439302 in chromosome 8p12 is significantly associated with papillary thyroid carcinoma (PTC) risk and dysregulated NRG1 expression. The underlying mechanisms remain to be discovered.

Objective

To evaluate the expression of NRG1 isoforms, candidate functional variants, and potential genes downstream of NRG1 in thyroid tissue.

Methods

Quantitative reverse transcription polymerase chain reaction was applied for gene expression analysis. SNaPshot assay, haplotype, and computer analyses were performed to evaluate candidate functional variants. Other functional assays [chromatin immunoprecipitation (ChIP) assay, luciferase assay, small interfering RNA knockdown, and RNA sequencing] were performed.

Results

Three NRG1 isoforms (NM_004495, NM_013958, and NM_001160008) tested were highly expressed in thyroid tissue. The expression levels of the three isoforms were significantly correlated with the genotypes of rs2439302. A DNA block of ~32 kb containing the risk G allele of rs2439302 was revealed, harboring multiple candidate functional variants. ChIP assay for active chromatin markers indicated at least nine regions in the DNA block showing strong H3Kme1 and H3K27Ac signals in thyroid tissue. Luciferase reporter assays revealed differential allelic activities associated with seven SNPs. Knocking down NRG1 in primary thyroid cells revealed downstream or interacting genes related to NRG1.

Conclusions

Our data suggest a role for transcriptional regulation of NRG1 in the predisposition to PTC.

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.

The role of SMAD3 in the genetic predisposition to papillary thyroid carcinoma

PURPOSE

To identify and characterize the functional variants, regulatory gene networks, and potential binding targets of SMAD3 in the 15q22 thyroid cancer risk locus.

METHODS

We performed linkage disequilibrium (LD) and haplotype analyses to fine map the 15q22 locus. Luciferase reporter assays were applied to evaluate the regulatory effects of the candidate variants. Knockdown by small interfering RNA, microarray analysis, chromatin immunoprecipitation (ChIP) and quantitative real-time polymerase chain reaction assays were performed to reveal the regulatory gene network and identify its binding targets.

RESULTS

We report a 25.6-kb haplotype within SMAD3 containing numerous single-nucleotide polymorphisms (SNPs) in high LD. SNPs rs17293632 and rs4562997 were identified as functional variants of SMAD3 by luciferase assays within the LD region. These variants regulate SMAD3 transcription in an allele-specific manner through enhancer elements in introns of SMAD3. Knockdown of SMAD3 in thyroid cancer cell lines revealed its regulatory gene network including two upregulated genes, SPRY4 and SPRY4-IT1. Sequence analysis and ChIP assays validated the actual binding of SMAD3 protein to multiple SMAD binding element sites in the region upstream of SPRY4.

CONCLUSION

Our data provide a functional annotation of the 15q22 thyroid cancer risk locus.

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.

A new panel of SNPs to assess thyroid carcinoma risk: a pilot study in a Brazilian admixture population

BACKGROUND

Thyroid cancer is a common malignant disease of the endocrine system with increasing incidence rates over the last few decades. In this study, we sought to analyze the possible association of 45 single nucleotide polymorphisms (SNPs) with thyroid cancer in a population from Rio Grande do Norte, Brazil.

METHODS

Based on histological analysis by a pathologist, 80 normal thyroid specimens of tissue adjacent to thyroid tumors were obtained from the biobank at the Laboratory of Pathology of Liga Norte Riograndense Contra o Câncer, Natal, RN. Patient samples were then genotyped using the MassARRAY platform (Sequenon, Inc) followed by statistical analysis employing the SNPassoc package in R program. The genotypic frequencies of all 45 SNPs obtained from the International HapMap Project database and based on data from the ancestral populations of European and African origin were used to compose the control study group.

RESULTS

In our study, the following 9 SNPs showed significant differences in their frequency when comparing the study and control groups: rs3744962, rs258107, rs1461855, rs4075022, rs9943744, rs4075570, rs2356508, rs17485896, and rs2651339. Furthermore, the SNPs rs374492 C/T and rs258107 C/T were associated with a relative risk for thyroid carcinoma of 3.78 (p = 6.27 × 10e^-5) and 2.91 (p = 8.27 × 10e^-5), respectively, after Bonferroni's correction for multiple comparisons.

CONCLUSIONS

These nine polymorphisms could be potential biomarkers of predisposition to thyroid carcinoma in the population from Rio Grande do Norte. However, complementary studies including a control group with samples obtained from healthy subjects in Rio Grande do Norte state, should be conducted to confirm these results.

Pathological processes and therapeutic advances in radioiodide refractory thyroid cancer

While in most patients with non-medullary thyroid cancer (TC), disease remission is achieved by thyroidectomy and ablation of tumor remnants by radioactive iodide (RAI), a substantial subgroup of patients with metastatic disease present tumor lesions that have acquired RAI resistance as a result of dedifferentiation. Although oncogenic mutations in BRAF, TERT promoter and TP53 are associated with an increased propensity for induction of dedifferentiation, the role of genetic and epigenetic aberrations and their effects on important intracellular signaling pathways is not yet fully elucidated. Also immune, metabolic, stemness and microRNA pathways have emerged as important determinants of TC dedifferentiation and RAI resistance. These signaling pathways have major clinical implications since their targeting could inhibit TC progression and could enable redifferentiation to restore RAI sensitivity. In this review, we discuss the current insights into the pathological processes conferring dedifferentiation and RAI resistance in TC and elaborate on novel advances in diagnostics and therapy to improve the clinical outcome of RAI-refractory TC patients.

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.

Genome-Wide Association Studies of Cancer in Diverse Populations

Genome-wide association studies (GWAS) of cancer have identified more than 700 risk loci, of which approximately 80% were first discovered in European ancestry populations, approximately 15% in East Asians, 3% in multiethnic scans, and less than 1% in African and Latin American populations. These percentages closely mirror the distribution of samples included in the discovery phase of cancer GWAS to date (84% European, 11% East Asian, 4% African, and 1% Latin American ancestry). GWAS in non-European ancestry populations have provided insight into ancestry-specific variation in cancer and have pointed to regions of susceptibility that are of particular importance in certain populations. Uncovering and characterizing cancer risk loci in diverse populations is critical for understanding underlying biological mechanisms and developing future genetic risk prediction models in non-European ancestry populations. New GWAS and continued collaborations will be required to eliminate population inequalities in the number of studies, sample sizes, and variant content on GWAS arrays, and to better align genetic research in cancer to the global distribution of race/ethnicity Cancer Epidemiol Biomarkers Prev; 27(4); 405-17. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in 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.

Thyroid-Stimulating Hormone, Thyroid Hormones, and Risk of Papillary Thyroid Cancer: A Nested Case-Control Study

Background: The effects of thyroid-stimulating hormone (TSH) and thyroid hormones on the development of human papillary thyroid cancer (PTC) remain poorly understood.Methods: The study population consisted of 741 (341 women, 400 men) histologically confirmed PTC cases and 741 matched controls with prediagnostic serum samples stored in the Department of Defense Serum Repository. Concentrations of TSH, total T3, total T4, and free T4 were measured in serum samples. Conditional logistic regression models were used to calculate ORs and 95% confidence intervals (CI).Results: The median time between blood draw and PTC diagnosis was 1,454 days. Compared with the middle tertile of TSH levels within the normal range, serum TSH levels below the normal range were associated with an elevated risk of PTC among women (OR, 3.74; 95% CI, 1.53-9.19) but not men. TSH levels above the normal range were associated with an increased risk of PTC among men (OR, 1.96; 95% CI, 1.04-3.66) but not women. The risk of PTC decreased with increasing TSH levels within the normal range among both men and women (P_trend = 0.0005 and 0.041, respectively).Conclusions: We found a significantly increased risk of PTC associated with TSH levels below the normal range among women and with TSH levels above the normal range among men. An inverse association between PTC and TSH levels within the normal range was observed among both men and women.Impact: These results could have significant clinical implications for physicians who are managing patients with abnormal thyroid functions and those with thyroidectomy. Cancer Epidemiol Biomarkers Prev; 26(8); 1209-18. ©2017 AACR.

A genome-wide association study of thyroid stimulating hormone and free thyroxine in Danish children and adolescents

Background

Hypothyroidism is associated with obesity, and thyroid hormones are involved in the regulation of body composition, including fat mass. Genome-wide association studies (GWAS) in adults have identified 19 and 6 loci associated with plasma concentrations of thyroid stimulating hormone (TSH) and free thyroxine (fT4), respectively.

Objective

This study aimed to identify and characterize genetic variants associated with circulating TSH and fT4 in Danish children and adolescents and to examine whether these variants associate with obesity.

Methods

Genome-wide association analyses of imputed genotype data with fasting plasma concentrations of TSH and fT4 from a population-based sample of Danish children, adolescents, and young adults, and a group of children, adolescents, and young adults with overweight and obesity were performed (N = 1,764, mean age = 12.0 years [range 2.5−24.7]). Replication was performed in additional comparable samples (N = 2,097, mean age = 11.8 years [1.2−22.8]). Meta-analyses, using linear additive fixed-effect models, were performed on the results of the discovery and replication analyses.

Results

No novel loci associated with TSH or fT4 were identified. Four loci previously associated with TSH in adults were confirmed in this study population (PDE10A (rs2983511: β = 0.112SD, p = 4.8 ∙ 10⁻¹⁶), FOXE1 (rs7847663: β = 0.223SD, p = 1.5 ∙ 10⁻²⁰), NR3C2 (rs9968300: β = 0.194SD), p = 2.4 ∙ 10⁻¹¹), VEGFA (rs2396083: β = 0.088SD, p = 2.2 ∙ 10⁻¹⁰)). Effect sizes of variants known to associate with TSH or fT4 in adults showed a similar direction of effect in our cohort of children and adolescents, 11 of which were associated with TSH or fT4 in our study (p<0.0002). None of the TSH or fT4 associated SNPs were associated with obesity in our cohort, indicating no pleiotropic effects of these variants on obesity.

Conclusion

In a group of Danish children and adolescents, four loci previously associated with plasma TSH concentrations in adults, were associated with plasma TSH concentrations in children, suggesting comparable genetic determinants of thyroid function in adults and children.

Effects of RET, NRG1 and NRG3 Polymorphisms in a Chinese Population with Hirschsprung Disease

The RET proto-oncogene was identified as a major locus involved in Hirschsprung disease (HSCR). A genome-wide association study (GWAS) and whole exome sequencing identified NRG1 and NRG3 as additional HSCR susceptibility loci. We investigated the effects of RET (rs2506030 and rs2435357), NRG1 (rs2439302, rs16879552 and rs7835688) and NRG3 (rs10748842, rs10883866 and rs6584400) polymorphisms in a Chinese population with HSCR. We assessed single nucleotide polymorphisms (SNPs) in the RET, NRG1 and NRG3 genes in a cohort of 362 sporadic HSCR patients and 1,448 normal controls using a TaqMan genotyping assay. Significant associations were found between HSCR risk and rs2506030, rs2435357, rs2439302 and rs7835688 (odds ratio [OR] 1.64, P = 1.72E-06; 2.97, P = 5.15E-33; 1.84, P = 9.36E-11; and 1.93, P = 1.88E-12, respectively). Two locus analyses of SNPs indicated increased disease risks of HSCR between NRG1 rs2439302 and RET rs2435357 or rs2506030. RET rs2506030 (GG genotype) and rs2435357 (TT genotype), in combination with NRG1 rs2439302 (GG genotype), were strongly associated with the highest risk of HSCR (OR = 56.53, P = 4.50E-07) compared with the two loci or a single SNP of either RET or NRG1. Our results support the association between genetic variation of RET and NRG1 and susceptibility to HSCR in the Chinese population.

Genetics of thyroid function

Recent studies show that subtle variations in thyroid function, including subclinical thyroid dysfunction, and even variation in thyroid function within the normal range, are associated with morbidity and mortality. It is estimated that 40-65% of the inter-individual variation in serum TSH and FT4 levels is determined by genetic factors. To identify these factors, various linkage and candidate gene studies have been performed in the past, which have identified only a few genes. In the last decade, genome-wide association studies identified many new genes, while recent whole-genome sequencing efforts have also been proven to be effective. In the current review, we provide a systematic overview of these studies, including strengths and limitations. We discuss new techniques which will further clarify the genetic basis of thyroid function in the near future, as well as the potential use of these genetic markers in personalizing the management of thyroid disease patients.

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.

MYH9 binds to lncRNA gene PTCSC2 and regulates FOXE1 in the 9q22 thyroid cancer risk locus

A locus on chromosome 9q22 harbors a SNP (rs965513) firmly associated with risk of papillary thyroid carcinoma (PTC). The locus also comprises the forkhead box E1 (FOXE1) gene, which is implicated in thyroid development, and a long noncoding RNA (lncRNA) gene, papillary thyroid cancer susceptibility candidate 2 (PTCSC2). How these might interact is not known. Here we report that PTCSC2 binds myosin-9 (MYH9). In a bidirectional promoter shared by FOXE1 and PTCSC2, MYH9 inhibits the promoter activity in both directions. This inhibition can be reversed by PTCSC2, which acts as a suppressor. RNA knockdown of FOXE1 in primary thyroid cells profoundly interferes with the p53 pathway. We propose that the interaction between the lncRNA, its binding protein MYH9, and the coding gene FOXE1 underlies the predisposition to PTC triggered by rs965513.

Association between IL-27 gene polymorphisms and risk of papillary thyroid carcinoma

AIM

IL-27 has potent antitumor effects. We aimed to examine the contribution of single nucleotide polymorphisms in IL-27 to the risk of papillary thyroid carcinoma (PTC).

MATERIALS & METHODS

 IL-27 rs153109 and rs17855750 were analyzed in 496 PTC patients and 629 controls, using a polymerase chain reaction-restriction fragment length polymorphism method.

RESULTS

The rs153109 AG and AG/GG genotypes were significantly associated with increased risks for PTC. Significantly increased PTC risk was also associated with rs17855750 GT and GT/GG genotypes. Combined genotypes of rs153109 AG/GG and rs17855750 GT/GG increased the risk of PTC (p < 0.05).

CONCLUSION

These findings showed that IL-27 rs153109 and rs17855750 might be related to the tumorigenesis of PTC.

DIRC3 and near NABP1 genetic polymorphisms are associated laryngeal squamous cell carcinoma patient survival

Laryngeal squamous cell carcinoma (LSCC) is one of the most common and aggressive malignancies of the upper digestive tract. The present study is a retrospective analysis of data from a prospective longitudinal study. A total of 170 male LSCC patients (average age, 60.75±10.082) at the First Affiliated Hospital of Xi'an Jiaotong University School of Medicine were recruited between January 2002 and April 2013 for this study. We assessed correlations between patient characteristics and survival, and sequenced genomic DNA from patient peripheral blood samples. We found that the single nucleotide polymorphisms (SNPs), rs11903757, with closest proximity to NABP1 and SDPR, and rs966423 in DIRC3, were associated with survival in LSCC patients. Median follow-up was 38 months (range 3-122) and median survival time was 48 months. LSCC patients with total laryngectomy, poor differentiation, T3-T4 stage, N1-N2 stage or III-IV TNM stage had reduced survival. This is the first study to demonstrate that the rs11903757 GT (HR=2.036; 95% CI, 1.071-3.872; p=0.030) and rs966423 TT (HR=11.677; 95% CI, 3.901-34.950; p=0.000) genotypes predict poor patient outcome. These polymorphisms may serve as useful clinical markers to predict patient survival, and to guide individual patient therapeutic decisions.

Thyroid Function Characteristics and Determinants: The Rotterdam Study

BACKGROUND

Information on determinants and change of thyroid function over time is sparse and conflicting but crucial for clinical interpretation and research. Therefore, our aim was to systematically investigate determinants of thyroid-stimulating hormone (TSH), free thyroxine (FT4) (as markers of thyroid function), their mutual relation (as marker of thyroid function set point) and changes in thyroid function over time.

METHODS

We included 9402 participants from the Rotterdam Study not taking thyroid medication and with available thyroid function measurements. Repeated measurements (6.5-year interval) were available for 1225 participants. The association of selected determinants with TSH, FT4, and their mutual relation (reflecting thyroid function set point) was estimated using linear regression models using restricted cubic splines with three knots. The factors investigated were age, sex, body mass index (BMI), tobacco smoking, alcohol use, thyroperoxidase antibodies (TPOAb), and common genetic factors.

RESULTS

Most influential determinants of TSH were age, smoking, genetic determinants, and TPOAb levels (p < 0.001). For FT4, most influential determinants were age, BMI, sex, genetic determinants and TPOAb levels (p < 0.001). Older age, female sex, and increased TPOAb levels were associated with a stronger relation between TSH and FT4. TSH levels did not change over time, irrespective of age. FT4 levels increased over time, most prominently in those older than 65 years of age (mean increase of 4.5 pmol/L).

CONCLUSIONS

The main factors that influence the relationship between thyroid hormone and molar concentrations of TSH in our population-based cohort study are age, smoking, BMI, TPOAb levels, and common genetic variants. The set point that determines TSH secretion as it relates to negative thyroid hormone feedback is modified by age, sex and TPOAb positivity. FT4 levels increase over time, with a more pronounced increase in the elderly, while TSH values seem stable over time. Our results question the current notion of an increase of TSH with increasing age.

Radiation Exposure and Thyroid Cancer Risk After the Fukushima Nuclear Power Plant Accident in Comparison with the Chernobyl Accident

The actual implementation of the epidemiological study on human health risk from low dose and low-dose rate radiation exposure and the comprehensive long-term radiation health effects survey are important especially after radiological and nuclear accidents because of public fear and concern about the long-term health effects of low-dose radiation exposure have increased considerably. Since the Great East Japan earthquake and the Fukushima Daiichi Nuclear Power Plant accident in Japan, Fukushima Prefecture has started the Fukushima Health Management Survey Project for the purpose of long-term health care administration and medical early diagnosis/treatment for the prefectural residents. Especially on a basis of the lessons learned from the Chernobyl accident, both thyroid examination and mental health care are critically important irrespective of the level of radiation exposure. There are considerable differences between Chernobyl and Fukushima regarding radiation dose to the public, and it is very difficult to estimate retrospectively internal exposure dose from the short-lived radioactive iodines. Therefore, the necessity of thyroid ultrasound examination in Fukushima and the intermediate results of this survey targeting children will be reviewed and discussed in order to avoid any misunderstanding or misinterpretation of the high detection rate of childhood thyroid cancer.

Papillary Thyroid Carcinoma: Association Between Germline DNA Variant Markers and Clinical Parameters

BACKGROUND

Papillary thyroid cancer (PTC) is reported to be highly heritable in epidemiological studies. Genome-wide association studies (GWAS) have uncovered several variants associated with PTC predisposition. It remains unknown whether these variants might contribute to better clinical stratification of PTC patients.

METHODS

In order to assess the usefulness of germline genetic analyses in the management of PTC patients, the genotypes of five variants (rs965513, rs944289, rs116909374, rs2439302, and rs966423) were determined in 1216 PTC patients and 1416 controls. Additionally, the expression of seven genes located close to GWAS variants (PTCSC3, MBIP, NKX2-1, FOXE1, DIRC3, PTCSC2, and NRG1) were measured in 73 PTC paired tumor/normal tissues, respectively. Next, the association was analyzed between the genotypes of the germline variants and the levels of gene expression with clinical/pathological features such as age, sex, TNM staging, multifocality status, extrathyroidal expansion, and MACIS score.

RESULTS

The risk allele of rs965513 was associated with larger tumor size (p = 0.025) and extrathyroidal expansion (odd ratio [OR] = 1.29, p = 0.045). The variant rs2439302 showed association with lymph node metastasis (OR = 1.24, p = 0.016), and multifocality status of the tumor (OR = 1.24, p = 0.012). The expression of MBIP was associated with T stage (p = 0.010). MBIP and PTCSC3 displayed lower expression in PTC tissue in males than in females (p = 0.025 and p = 0.036, respectively). NKX2-1 displayed lower expression in patients with N1 stage (p = 0.040).

CONCLUSIONS

The studied germline risk alleles predisposing to PTC were associated with a more aggressive course of the disease reflected by larger tumor diameter, higher multifocality rate, and more advanced N stage at the time of diagnosis. These results show that germline variants not only predispose to PTC but also might impact its clinical course. However, these associations were only moderate, and further large multi-ethnic studies are required to evaluate the usefulness of these germline variants in the clinical stratification of PTC patients.