Genome-wide association study on differentiated thyroid cancer

Causal association of sleep traits with the risk of thyroid cancer: A mendelian randomization study

BACKGROUND

This study was to explore the causal associations of sleep traits including sleep duration, snoring, chronotype, sleep disorders, getting up in the morning, sleeplessness/insomnia and nap during day with the risk of thyroid cancer based on Mendelian randomization (MR) analysis.

METHOD

Summary single nucleotide polymorphism (SNP)-phenotype association data were obtained from published genome-wide association studies (GWASs) using the FinnGen and UK Biobank databases. A series of screening processes were performed to select qualified SNPs strongly related to exposure. We applied the inverse variance weighted (IVW), the Mendelian Randomization robust adjusted profile score (MR-RAPS), the Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO), and the Weighted Median to estimate the causal links between sleep traits and the risk of thyroid cancer. Odds ratio (OR) and 95% confidence interval (CI) were calculated.

RESULTS

The IVW results showed that getting up in the morning (OR = 0.055, 95%CI: 0.004-0.741) and napping during day (OR = 0.031, 95%CI: 0.002-0.462) were associated with decreased risk of thyroid cancer in the Italian population. A 1.30-h decrease of sleep duration was associated with 7.307-fold of thyroid cancer risk in the Finnish population (OR = 7.307, 95%CI: 1.642-32.519). Cronotype could decrease the risk of thyroid cancer in the Finnish population (OR = 0.282, 95%CI: 0.085-0.939). Sleep disorders increased the risk of thyroid cancer in the Finnish population (OR = 2.298, 95%CI: 1.194-4.422). The combined results revealed that sleep duration was correlated with increased risk of thyroid cancer (OR = 5.600, 95%CI: 1.458-21.486).

CONCLUSION

Decreased sleep duration was associated with increased risk of thyroid cancer, which indicated the importance of adequate sleep for the prevention of thyroid cancer.

Causal relationship between gut microbiota and differentiated thyroid cancer: a two-sample Mendelian randomization study

Background

The gut microbiota has been significantly associated with differentiated thyroid cancer (DTC). However, the causal relationship between the gut microbiota and DTC remains unexplored.

Methods

Genome-wide association study (GWAS) summary databases were utilized to select exposures and outcomes. The Mendelian randomization (MR) method was employed to investigate the causal relationship between the gut microbiota and DTC. A sensitivity analysis was performed to assess the reliability of the findings.

Results

Four bacterial traits were associated with the risk of DTC: Class Mollicutes [odds ratio (OR) = 10.953, 95% confidence interval (95% CI): 2.333-51.428, p = 0.002], Phylum Tenericutes (OR = 10.953, 95% CI: 2.333-51.428, p = 0.002), Genus Eggerthella (OR = 3.219, 95% CI: 1.033-10.024, p = 0.044), and Order Rhodospirillales (OR = 2.829, 95% CI: 1.096-7.299, p = 0.032). The large 95% CI range for the Class Mollicutes and the Phylum Tenericutes may be attributed to the small sample size. Additionally, four other bacterial traits were negatively associated with DTC: Genus Eubacterium fissicatena group (OR = 0.381, 95% CI: 0.148-0.979, p = 0.045), Genus Lachnospiraceae UCG008 (OR = 0.317, 95% CI: 0.125-0.801, p = 0.015), Genus Christensenellaceae R-7 group (OR = 0.134, 95% CI: 0.020-0.886, p = 0.037), and Genus Escherichia Shigella (OR = 0.170, 95% CI: 0.037-0.769, p = 0.021).

Conclusion

These findings contribute to our understanding of the pathological mechanisms underlying DTC and provide novel insights for the clinical treatment of DTC.

Thyroid cancer and cardiovascular diseases: a Mendelian randomization study

Background

Multiple observational studies have shown associations between thyroid cancer (TC) and cardiovascular diseases (CVDs). However, the results were inconsistent, and the potential causal genetic relationship remains unclear.

Methods

The genetic instruments of TC and CVDs were derived from data obtained through genome-wide association studies (GWAS). We performed the two-sample Mendelian randomization(MR) methods to investigate the causality of TC on CVDs. Summary-level statistics for CVDs, including heart failure (HF), atrial fibrillation (AF), coronary artery disease (CAD), myocardial infarction (MI), ischemic stroke (IS) and venous thromboembolism (VTE). The primary method employed in this MR analysis was the Inverse Variance Weighted (IVW) approach, and four additional algorithms were used: MR-Egger, weighted median, simple mode, and weighted mode. Additionally, we assessed the reliability of the causal relationship through pleiotropy, heterogeneity and leave-one-out sensitivity analysis.

Results

In this MR analysis, we only detected causality of genetically predicted TC on HF (IVW method, odds ratio (OR) = 1.00134, 95% confidence interval (CI): 1.00023-1.00244, p = 0.017). However, There were no causal associations of TC with CAD, MI, AF, IS, and VTE.

Conclusion

Our results confirmed the causal association between TC and HF. It is crucial to closely monitor the incidence of HF in TC patients and give comprehensive clinical intervention based on conventional treatment.

Association of major depression, schizophrenia and bipolar disorder with thyroid cancer: a bidirectional two-sample mendelian randomized study

BACKGROUND

Major depressive disease (MDD), schizophrenia (SCZ), and bipolar disorder (BD) are common psychiatric disorders, and their relationship with thyroid cancer has been of great interest. This study aimed to investigate the potential causal effects of MDD, SCZ, BD, and thyroid cancer.

METHODS

We used publicly available summary statistics from large-scale genome-wide association studies to select genetic variant loci associated with MDD, SCZ, BD, and thyroid cancer as instrumental variables (IVs), which were quality controlled and clustered. Additionally, we used three Mendelian randomization (MR) methods, inverse variance weighted (IVW), MR-Egger regression and weighted median estimator (WME) methods, to estimate the bidirectional causal relationship between psychiatric disorders and thyroid cancer. In addition, we performed heterogeneity and multivariate tests to verify the validity of the IVs.

RESULTS

We used two-sample bidirectional MR analysis to determine whether there was a positive causal association between MDD and thyroid cancer risk. The results of the IVW analysis (OR = 3.956 95% CI = 1.177-13.299; P = 0.026) and the WME method (OR = 5.563 95% CI = 0.998-31.008; P = 0.050) confirmed that MDD may increase the risk of thyroid cancer. Additionally, our study revealed a correlation between genetic susceptibility to SCZ and thyroid cancer (OR = 1.532 95% CI = 1.123-2.088; P = 0.007). The results of the WME method analysis based on the median estimate (OR = 1.599 95% CI = 1.014-2.521; P = 0.043) also suggested that SCZ may increase the risk of thyroid cancer. Furthermore, our study did not find a causal relationship between BD and thyroid cancer incidence. In addition, the results of reverse MR analysis showed no significant causal relationships between thyroid cancer and MDD, SCZ, or BD (P > 0.05), ruling out the possibility of reverse causality.

CONCLUSIONS

This MR method analysis provides new evidence that MDD and SCZ may be positively associated with thyroid cancer risk while also revealing a correlation between BD and thyroid cancer. These results may have important implications for public health policy and clinical practice. Future studies will help elucidate the biological mechanisms of these associations and potential confounders.

Association of congenital iodine deficiency syndrome and differentiated thyroid cancer: a Mendelian randomization study

BACKGROUND

The effect of iodine deficiency, especially during the fetal period, on thyroid cancer risk remains unclear. The evidence from observational studies is controversial because of the inevitable confounding factors. We studied the causal effect of congenital iodine deficiency on differentiated thyroid cancer (DTC) based on Mendelian randomization (MR).

METHODS

Two-Sample MR analysis was performed using data from published genome-wide association studies, including congenital iodine deficiency syndrome (CIDS) (353 cases, 187,684 controls) and DTC (649 cases, 431 controls) data.

RESULTS

There was a causal relationship between CIDS and DTC (P < 0.05), with CIDS increasing the DTC risk by 37.4% (OR = 1.374, 95%CI = 1.110-1.700). Heterogeneity tests and tests of multiple validities indicated that the results were solid and reliable (all P < 0.05).

CONCLUSIONS

Fetal iodine deficiency increases the risk of DTC, so future clinical studies should focus on the effect of iodine supplementation during pregnancy to reduce the risk of thyroid cancer in the offspring.

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.

Thyroid Cancer Polygenic Risk Score Improves Classification of Thyroid Nodules as Benign or Malignant

CONTEXT

Thyroid nodule ultrasound-based risk stratification schemas rely on the presence of high-risk sonographic features. However, some malignant thyroid nodules have benign appearance on thyroid ultrasound. New methods for thyroid nodule risk assessment are needed.

OBJECTIVE

We investigated polygenic risk score (PRS) accounting for inherited thyroid cancer risk combined with ultrasound-based analysis for improved thyroid nodule risk assessment.

METHODS

The convolutional neural network classifier was trained on thyroid ultrasound still images and cine clips from 621 thyroid nodules. Phenome-wide association study (PheWAS) and PRS PheWAS were used to optimize PRS for distinguishing benign and malignant nodules. PRS was evaluated in 73 346 participants in the Colorado Center for Personalized Medicine Biobank.

RESULTS

When the deep learning model output was combined with thyroid cancer PRS and genetic ancestry estimates, the area under the receiver operating characteristic curve (AUROC) of the benign vs malignant thyroid nodule classifier increased from 0.83 to 0.89 (DeLong, P value = .007). The combined deep learning and genetic classifier achieved a clinically relevant sensitivity of 0.95, 95% CI [0.88-0.99], specificity of 0.63 [0.55-0.70], and positive and negative predictive values of 0.47 [0.41-0.58] and 0.97 [0.92-0.99], respectively. AUROC improvement was consistent in European ancestry-stratified analysis (0.83 and 0.87 for deep learning and deep learning combined with PRS classifiers, respectively). Elevated PRS was associated with a greater risk of thyroid cancer structural disease recurrence (ordinal logistic regression, P value = .002).

CONCLUSION

Augmenting ultrasound-based risk assessment with PRS improves diagnostic accuracy.

Exploring the Potential Association between COVID-19 and Thyroid Cancer: A Mendelian Randomization Study

In order to address the ongoing debate surrounding the potential link between COVID-19 and thyroid cancer, our study was specifically designed to investigate the association between these two factors. We acquired summary data from a genome-wide association study (GWAS) concerning COVID-19 susceptibility and severity of COVID-19 in the European population, with a focus on their relationship with thyroid cancer. We applied three distinct methodologies to evaluate the causality between COVID-19 and thyroid cancer, employing Mendelian randomization (MR)-Egger, weighted median (WM), and inverse variance-weighted (IVW) approaches. Furthermore, we utilized a variety of techniques to assess pleiotropy and heterogeneity, including the MR-Egger intercept, MR-pleiotropy residual sum and outlier method (PRESSO), and Cochran's Q test. The MR analysis revealed associations between the susceptibility of COVID-19 and thyroid cancer (IVW odds ratio [OR]: 2.826, 95% confidence interval [CI]: [0.842, 9.483], P = 0.093) as well as between the risk of COVID-19 hospitalization and thyroid cancer (IVW OR: 1.630, 95% CI: [1.050, 2.529], P = 0.029). However, the relationship between COVID-19 and the occurrence of severe thyroid cancer cases was less evident (IVW OR: 1.061, 95% CI: [0.575, 1.956], P = 0.850). Our sensitivity analyses did not reveal any signs of horizontal pleiotropy or heterogeneity. Our MR study provided compelling evidence supporting a causal connection between the risk of COVID-19 hospitalization and thyroid cancer. Nevertheless, the MR results derived from genetic data do not support a causal link between susceptibility to COVID-19 and the risk of thyroid cancer or between very severe cases of COVID-19 and the risk of thyroid cancer. These findings have significant implications for further investigations into the impact of COVID-19 on health and the etiology of thyroid cancer.

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.

Assessment of causal association between differentiated thyroid cancer and disordered serum lipid profile: a Mendelian randomization study

Background

Research has shown that the disordered serum lipid profile may be associated with the risk of differentiated thyroid cancer (DTC). Whether this association reflect causal effect is still unclear. The aim of this study was to evaluate the causality of circulating lipoprotein lipids on DTC.

Methods

Mendelian randomization (MR) analysis was conducted to evaluate the relationship between the circulating lipoprotein lipids and DTC risk using single-nucleotide polymorphisms (SNPs) from a genome-wide association (GWA) study containing a high-incidence Italian population of 690 cases samples with DTC and 497 controls.

Results

Univariate and multivariate mendelian randomization analysis demonstrated that 'total cholesterol', 'HDL cholesterol', 'apolipoprotein B' and 'ratio of apolipoprotein B to apolipoprotein A1' were correlated with DTC. According to sensitivity analysis, our results were reliable. Furthermore, multivariate analysis revealed that there is no causative association between DTC and any of the many cause factors when they interact with one another, suggesting that there was a deep interaction between the four factors, which could affect each other. Finally, the mechanism of the related effects each other as well as the target genes with significant SNP regulatory effects in DTC was explored by conducting functional enrichment analysis and constructing the regulatory networks.

Conclusions

We obtained four exposure factors (total cholesterol, HDL cholesterol, apolipoprotein B and ratio of apolipoprotein B to apolipoprotein A1) closely related to DTC, which laid a theoretical foundation for the treatment of DTC.

Causal analysis of the gut microbiota in differentiated thyroid carcinoma: a two-sample Mendelian randomization study

Objective: Numerous studies have highlighted an association between the gut microbiota (GM) and thyroid tumors. Employing Mendelian randomization methodology, we seek to elucidate the causal link between the gut microbiota and thyroid neoplasms. Methods: We procured data from the Mibiogen database encompassing 211 distinct gut microbiota taxa, alongside extensive genome-wide association studies (GWAS) summary data for differentiated thyroid carcinoma (DTC). Our principal analytical approach involved the application of the Inverse-Variance Weighted method (IVW) within the framework of Mendelian randomization. Simultaneously, we conducted sensitivity analyses to assess result heterogeneity, horizontal pleiotropy, and outcome stability. Results: IVW analysis revealed a dual role of the GM in thyroid carcinoma. The phylum Actinobacteria (OR, 0.249 [95% CI, 0.121-0.515]; p < 0.001) was associated with a decreased risk of DTC. Conversely, the genus Ruminiclostridium9 (OR, 11.276 [95% CI, 4.406-28.860]; p < 0.001), class Mollicutes (OR, 5.902 [95% CI, 1.768-19.699]; p = 0.004), genus RuminococcaceaeUCG004 (OR, 3.831 [95% CI, 1.516-9.683]; p = 0.005), genus Paraprevotella (OR, 3.536 [95% CI, 1.330-9.401]; p = 0.011), and phylum Tenericutes (OR, 5.902 [95% CI, 1.768-19.699]; p = 0.004) were associated with an increased risk of DTC. Conclusion: Our findings underscore that the presence of genus Ruminiclostridium9, class Mollicutes, genus RuminococcaceaeUCG004, genus Paraprevotella, and phylum Tenericutes is associated with an elevated risk of DTC, whereas the presence of the phylum Actinobacteria is linked to a decreased risk. These discoveries enhance our comprehension of the relationship between the GM and DTC.

Unbalanced bidirectional causal association between thyroid cancer and ER-positive breast cancer: should we recommend screening for thyroid cancer in breast cancer patients?

BACKGROUND

The association between breast cancer (BC) and thyroid cancer (TC) has been studied in several epidemiological studies. However, the underlying causal relationship between them is not yet clear.

METHODS

The data from the latest large-sample genome-wide association studies (GWAS) of BC and TC were searched in the public GWAS database. The BC GWAS data included estrogen receptor (ER)-positive and negative subgroups. Two-way two-sample Mendelian Randomization (MR) was used to explore the potential causal relationship between BC and TC. Inverse variance weighting (IVW) and the MR-Egger method were used to combine the estimation of each single nucleotide variation (previous single nucleotide polymorphism). BC was taken as the result, and the effect of TC exposure was analyzed. Then, the effect of BC exposure on the result of TC was analyzed.

RESULTS

Both IVW and MR-Egger results indicated that gene-driven thyroid cancer does not cause estrogen receptor-positive breast cancer and is a protective factor (β = -1.203, SE = 4.663*10-4, P = 0.010). However, gene-driven estrogen receptor-positive breast cancer can lead to the development of thyroid cancer (β = 0.516, SE = 0.220, P = 0.019).

CONCLUSION

From the perspective of gene drive, people with TC are less likely to have ER-positive BC. In contrast, people with ER-positive BC are more likely to have TC. Therefore, it is recommended that patients with BC be screened regularly for TC.

Causal associations between schizophrenia and cancers risk: a Mendelian randomization study

Background

Previous observational studies have reported inconsistent findings regarding the incidence of cancer in patients with schizophrenia compared to the general population. The causal relationship between schizophrenia and cancer remains unclear and requires further investigation.

Objective

To investigate the causal relationship between schizophrenia and cancer.

Methods

In this study, a two-sample Mendelian randomization (MR) analysis was conducted using publicly available genome-wide association studies to determine the causal relationship. The effect estimates were calculated using the random-effects inverse-variance-weighted method.

Results

We determined a causal relationship between genetic predisposition to schizophrenia and cancer, with schizophrenia increasing lung cancer (odds ratio (OR) = 1.0007; 95% confidence interval (CI), 1.0001-1.0013; p = 0.0192), thyroid cancer (OR = 1.5482; CI, 1.1112-2.1569; p =0.0098),colorectal cancer (OR = 1.0009; CI, 1.0001-1.0018; p = 0.0344), ovarian cancer (OR = 1.0770; CI, 1.0352-1.1203; p = 0.0002), breast cancer (OR = 1.0011; CI, 1.0001- 1.0022; p =0.0352) and reduced the risk of malignant neoplasm of the stomach (OR = 0.8502; CI, 0.7230-0.9998; p = 0.0496).

Conclusions

This study conducted a two-sample MR analysis and discovered a positive causal relationship between schizophrenia and breast, ovarian, thyroid, lung, and colorectal cancers. On the other hand, an inverse causal relationship was found between schizophrenia and malignant neoplasm of the stomach.

Role of ST6GAL1 in Thyroid Cancers: Insights from Tissue Analysis and Genomic Datasets

Thyroid cancer is the predominant endocrine-related malignancy. ST6 β-galactoside α2,6-sialyltransferase 1 (ST6GAL1) has been studied in various types of cancers; however, the expression and function of ST6GAL1 in thyroid cancer has not been investigated so far. Previously, we conducted two genome-wide association studies and have identified the association of the ST6GAL1 gene with plasma thyroglobulin (Tg) levels. Since Tg levels are altered in thyroid pathologies, in the current study, we wanted to evaluate the expression of ST6GAL1 in thyroid cancer tissues. We performed an immunohistochemical analysis using human thyroid tissue from 89 patients and analyzed ST6GAL1 protein expression in papillary thyroid cancer (including follicular variant and microcarcinoma) and follicular thyroid cancer in comparison to normal thyroid tissue. Additionally, ST6GAL1 mRNA levels from The Cancer Genome Atlas (TCGA, n = 572) and the Genotype-Tissue Expression (GTEx) project (n = 279) were examined. The immunohistochemical analysis revealed higher ST6GAL1 protein expression in all thyroid tumors compared to normal thyroid tissue. TCGA data revealed increased ST6GAL1 mRNA levels in both primary and metastatic tumors versus controls. Notably, the follicular variant of papillary thyroid cancer exhibited significantly higher ST6GAL1 mRNA levels than classic papillary thyroid cancer. High ST6GAL1 mRNA levels significantly correlated with lymph node metastasis status, clinical stage, and reduced survival rate. ST6GAL1 emerges as a potential cancer-associated glycosyltransferase in thyroid malignancies, offering valuable insights into its diagnostic and prognostic significance.

Genetic link between primary sclerosing cholangitis and thyroid dysfunction: a bidirectional two-sample Mendelian randomization study

Background

Observational studies have demonstrated an association between primary sclerosing cholangitis (PSC) and thyroid dysfunction (TD). However, the causal relationship between PSC and TD remains uncertain. The purpose of this study is to investigate the causal associations and specific direction between these two conditions. Gaining insight into the potential causal relationship between PSC and TD is valuable for elucidating the pathogenesis of PSC and for devising innovative approaches for the prevention and treatment of PSC and its associated complications.

Methods

We conducted a bidirectional two-sample Mendelian randomization (MR) analysis to investigate the causal association between PSC and TD, such as autoimmune thyroid disease (AITD), thyroid cancer (TC), thyroid stimulating hormone (TSH), thyrotropin-releasing hormone (TRH), among others. PSC was the exposure variable, while TD was the outcome variable. To identify suitable instrumental variables (IVs), we utilized genome-wide association study (GWAS) datasets to select potential candidate single-nucleotide polymorphisms (SNPs). The primary statistical approach employed was the inverse-variance weighted (IVW) method, which was complemented by a series of sensitivity analyses to assess the robustness of the results by estimating heterogeneity and pleiotropy.

Results

We found that the causal associations between genetically predicted PSC and Graves' disease (GD), hyperthyroidism (IVW OR=1.230, 95%CI: 1.089-1.389, P=0.001; IVW OR=1.001, 95%CI: 1.000-1.002, P=0.000) were statistically significant. The reverse MR analysis indicated that genetic susceptibility to hyperthyroidism (P=0.000) and hypothyroidism (p=0.028) might be the risk of PSC. There was no statistically significant causal association observed between PSC and other TD (IVW P>0.05), with the exception of GD, hyperthyroidism, and hypothyroidism as determined through bidirectional two-sample analysis. To ensure the reliability of our findings, additional sensitivity analyses were conducted, including the leave-one-out (LOO) test, heterogeneity test, and pleiotropic test.

Conclusion

In this study, we conducted an investigation into the causal association between PSC and TD. Our findings indicate that PSC significantly elevates the susceptibility to GD and hyperthyroidism from a statistical perspective. These results shed light on the etiology of PSC and have implications for the management of patients with PSC.

Causal association between inflammatory bowel disease and 32 site-specific extracolonic cancers: a Mendelian randomization study

BACKGROUND

The risk of extracolonic cancer is increased in inflammatory bowel disease (IBD) patients, but it is not clear whether there is a causal relationship. We aimed to systematically estimate the causal relationship between IBD and extracolonic cancers.

METHODS

Independent genetic variants strongly associated with IBD were extracted as instrumental variables from genome-wide association study (GWAS) conducted by the International IBD Genetics Consortium including 12,882 IBD patients, 5956 Crohn's disease (CD) patients, and 6968 ulcerative colitis (UC) patients. Three sources of cancer GWAS were selected as outcome data. Two-sample Mendelian randomization (MR) analysis was conducted to assess the causal effects of IBD on 32 extracolonic cancers. The meta-analysis was applied to assess the combined causal effect with multiple MR results.

RESULTS

IBD, CD, and UC have potential causal associations with oral cavity cancer (IBD: OR = 1.180, 95% CI: 1.059 to 1.316, P = 0.003; CD: OR = 1.112, 95% CI: 1.008 to 1.227, P = 0.034; UC: OR = 1.158, 95% CI: 1.041 to 1.288, P = 0.007). Meta-analysis showed a significant positive causal relationship between IBD and breast cancer (OR = 1.059; 95% CI: 1.033 to 1.086; P < 0.0001) as well as a potential causal relationship between CD and breast cancer (OR = 1.029; 95% CI: 1.002 to 1.055; P = 0.032) based on combining multiple MR results.

CONCLUSIONS

This comprehensive MR analysis suggested that genetically predicted IBD, as well as its subtypes, may be a risk factor in the development of oral cavity and breast cancer.

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Causal associations of Sjögren's syndrome with cancers: a two-sample Mendelian randomization study

BACKGROUND

Several observational studies have explored the associations between Sjögren's syndrome (SS) and certain cancers. Nevertheless, the causal relationships remain unclear. Mendelian randomization (MR) method was used to investigate the causality between SS and different types of cancers.

METHODS

We conducted the two-sample Mendelian randomization with the public genome-wide association studies (GWASs) summary statistics in European population to evaluate the causality between SS and nine types of cancers. The sample size varies from 1080 to 372,373. The inverse variance weighted (IVW) method was used to estimate the causal effects. A Bonferroni-corrected threshold of P < 0.0031 was considered significant, and P value between 0.0031 and 0.05 was considered to be suggestive of an association. Sensitivity analysis was performed to validate the causality. Moreover, additional analysis was used to assess the associations between SS and well-accepted risk factors of cancers.

RESULTS

After correcting the heterogeneity and horizontal pleiotropy, the results indicated that patients with SS were significantly associated with an increased risk of lymphomas (odds ratio [OR] = 1.0010, 95% confidence interval [CI]: 1.0005-1.0015, P = 0.0002) and reduced risks of prostate cancer (OR = 0.9972, 95% CI: 0.9960-0.9985, P = 2.45 × 10-5) and endometrial cancer (OR = 0.9414, 95% CI: 0.9158-0.9676, P = 1.65 × 10-5). Suggestive associations were found in liver and bile duct cancer (OR = 0.9999, 95% CI: 0.9997-1.0000, P = 0.0291) and cancer of urinary tract (OR = 0.9996, 95% CI: 0.9992-1.0000, P = 0.0281). No causal effect of SS on other cancer types was detected. Additional MR analysis indicated that causal effects between SS and cancers were not mediated by the well-accepted risk factors of cancers. No evidence of the causal relationship was observed for cancers on SS.

CONCLUSIONS

SS had significant causal relationships with lymphomas, prostate cancer, and endometrial cancer, and suggestive evidence of association was found in liver and bile duct cancer and cancer of urinary tract, indicating that SS may play a vital role in the incidence of these malignancies.

Immp2l Mutation Induces Mitochondrial Membrane Depolarization and Complex III Activity Suppression after Middle Cerebral Artery Occlusion in Mice

OBJECTIVE

We previously reported that mutations in inner mitochondrial membrane peptidase 2-like (Immp2l) increase infarct volume, enhance superoxide production, and suppress mitochondrial respiration after transient cerebral focal ischemia and reperfusion injury. The present study investigated the impact of heterozygous Immp2l mutation on mitochondria function after ischemia and reperfusion injury in mice.

METHODS

Mice were subjected to middle cerebral artery occlusion for 1 h followed by 0, 1, 5, and 24 h of reperfusion. The effects of Immp2l^+/- on mitochondrial membrane potential, mitochondrial respiratory complex III activity, caspase-3, and apoptosis-inducing factor (AIF) translocation were examined.

RESULTS

Immp2l^+/- increased ischemic brain damage and the number of TUNEL-positive cells compared with wild-type mice. Immp2l^+/- led to mitochondrial damage, mitochondrial membrane potential depolarization, mitochondrial respiratory complex III activity suppression, caspase-3 activation, and AIF nuclear translocation.

CONCLUSION

The adverse impact of Immp2l^+/- on the brain after ischemia and reperfusion might be related to mitochondrial damage that involves depolarization of the mitochondrial membrane potential, inhibition of the mitochondrial respiratory complex III, and activation of mitochondria-mediated cell death pathways. These results suggest that patients with stroke carrying Immp2l^+/- might have worse and more severe infarcts, followed by a worse prognosis than those without Immp2l mutations.

Metachronous papillary thyroid carcinoma and hereditary breast carcinoma with BRCA1 mutation in a Vietnamese woman: A case report

INTRODUCTION AND IMPORTANCE

Papillary thyroid carcinoma is the most common type of malignancy in endocrine tumor. Women with breast cancer have an increased risk of developing thyroid cancer. Especially, patients with BRCA1 germline variants, which belong to the DNA DSB repair system, there may be a genetic susceptibility to thyroid cancer.

CASE PRESENTATION

This study investigated a 47-year-old Vietnamese female patient with BRCA1 mutation, namely NM_007294.3 (BRCA1): c.4998insA (p. Tyr1666Terfs), who developed PTC after one year of breast cancer treatment.

CLINICAL DISCUSSION

Factors that are thought to increase the risk of thyroid cancer after breast cancer treatment include treatment methods, hormonal factors, genetic susceptibility, and others. A hypothesis is breast cancer with BRCA1 mutation increases the risk of thyroid cancer.

CONCLUSION

Understanding the relationship between breast and thyroid cancer helps clinicians for well management of both diseases and also contributes to the development of new preventive strategies, diagnostics and therapeutics as so as a new research direction.

Assessing the causality between thyroid and breast neoplasms: A bidirectional Mendelian randomization study

Aim

This study aimed to evaluate the association between thyroid neoplasms (TN) and the risk of developing breast neoplasms (BN) by assessing data on single nucleotide polymorphisms (SNPs) obtained from the Deutsches Krebsforschungszentrum (DKFZ) and Breast Cancer Association (BCAC).

Methods

Data on SNPs associated with TN and BN were obtained from DKFZ and BCAC, respectively. Secondary data analysis of all pooled data from genome-wide association studies (GWAS) was performed to identify the genetic loci closely associated with TN or BN as instrumental variables (IVs). To evaluate the causal relationship between TN and BN, a bidirectional Mendelian randomization (MR) analysis was performed using MR Egger regression, weighted median, inverse variance weighted (IVW) random effects model, simple mode, weighted mode, maximum likelihood, penalized weighted median, IVW radial, IVW fixed effects, and robust adjusted profile scores (RAPS) method.

Results

The MR in this study demonstrated a modest reverse causal relationship between TN and BN but a significant positive causal relationship between BN and TN.

Conclusions

The MR of this study provided genetic evidence suggesting an association between BN and TN; however, further research is warranted to explore the potential mechanism of interaction between these two malignancies. Moreover, general breast screening should be performed in individuals with TN, but TN screening should be reinforced in individuals with BN.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

Background

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

Case presentation

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

Conclusions

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

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.

Cancers and COVID-19 Risk: A Mendelian Randomization Study

Simple Summary

During the COVID-19 pandemic, cancer patients are regarded as a highly vulnerable population. Given the unavoidable bias and unmeasured confounders in observational studies, the causal effects of cancers on COVID-19 outcomes are largely unknown. In the study, we tried to evaluate the causal effects of cancers on COVID-19 outcomes using the Mendelian randomization (MR) approach. No strong evidence was observed to support a causal role of cancer in COVID-19 development. Previous observational correlations between cancers and COVID-19 outcomes were likely confounded. Large and well-conducted epidemiological studies are required to determine whether cancers causally contribute to increased risk of COVID-19.

Abstract

Observational studies have shown increased COVID-19 risk among cancer patients, but the causality has not been proven yet. Mendelian randomization analysis can use the genetic variants, independently of confounders, to obtain causal estimates which are considerably less confounded. We aimed to investigate the causal associations of cancers with COVID-19 outcomes using the MR analysis. The inverse-variance weighted (IVW) method was employed as the primary analysis. Sensitivity analyses and multivariable MR analyses were conducted. Notably, IVW analysis of univariable MR revealed that overall cancer and twelve site-specific cancers had no causal association with COVID-19 severity, hospitalization or susceptibility. The corresponding p-values for the casual associations were all statistically insignificant: overall cancer (p = 0.34; p = 0.42; p = 0.69), lung cancer (p = 0.60; p = 0.37; p = 0.96), breast cancer (p = 0.43; p = 0.74; p = 0.43), endometrial cancer (p = 0.79; p = 0.24; p = 0.83), prostate cancer (p = 0.54; p = 0.17; p = 0.58), thyroid cancer (p = 0.70; p = 0.80; p = 0.28), ovarian cancer (p = 0.62; p = 0.96; p = 0.93), melanoma (p = 0.79; p = 0.45; p = 0.82), small bowel cancer (p = 0.09; p = 0.08; p = 0.19), colorectal cancer (p = 0.85; p = 0.79; p = 0.30), oropharyngeal cancer (p = 0.31; not applicable, NA; p = 0.80), lymphoma (p = 0.51; NA; p = 0.37) and cervical cancer (p = 0.25; p = 0.32; p = 0.68). Sensitivity analyses and multivariable MR analyses yielded similar results. In conclusion, cancers might have no causal effect on increasing COVID-19 risk. Further large-scale population studies are needed to validate our findings.

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.

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.

Association between dried fruit intake and pan-cancers incidence risk: A two-sample Mendelian randomization study

Background

Observational studies have revealed that dried fruit intake may be associated with cancer incidence; however, confounding factors make the results prone to be disturbed. Therefore, we conducted a two-sample Mendelian randomization (MR) study to explore the causal relationship between dried fruit intake and 11 site-specific cancers.

Materials and methods

Forty-three single nucleoside polymers (SNPs) with robust genome-wide association study (GWAS) evidence, strongly correlated with dried fruit intake, were used as instrumental variables (IVs) in this study. The summary-level genetic datasets of site-specific cancers were obtained from the Oncoarray oral cavity and oropharyngeal cancer consortium, International Lung Cancer Consortium, Breast Cancer Association Consortium (BCAC), Ovarian Cancer Association Consortium, PanScan1, and GWAS of other scholars. We analyzed the causality between dried fruit intake and 11 site-specific cancers using the inverse-variance-weighted (IVW) and weighted median (WM) methods. For the results of the MR analysis, Cochran's Q test was used to check for heterogeneity, and multiplicative random effects were used to evaluate the heterogeneity further. Gene pleiotropy was tested using MR-Egger regression and MR-PRESSO methods. In addition, the main results of this study were validated by using the summary statistical data from the FinnGen and UK Biobank databases, and adjusted body mass index (BMI), years of education, fresh fruit intake, and vitamin C using multivariable MR analysis to ensure the stability of the research results.

Results

The evidence from IVW analyses showed that each increase of dried fruit intake by one standard deviation was statistically significantly associated with 82.68% decrease of oral cavity/pharyngeal cancer incidence risk (P = 0.0131), 67.01% decrease of lung cancer incidence risk (P = 0.0011), 77% decrease of squamous cell lung cancer incidence risk (P = 0.0026), 53.07% decrease of breast cancer incidence risk (P = 4.62 × 10^-5), 39.72% decrease of ovarian cancer incidence risk (P = 0.0183), 97.26% decrease of pancreatic cancer incidence risk (P = 0.0280), 0.53% decrease of cervical cancer incidence risk (P = 0.0482); however, there was no significant effect on lung adenocarcinoma (P = 0.4343), endometrial cancer (P = 0.8742), thyroid cancer (P = 0.6352), prostate cancer (P = 0.5354), bladder cancer (P = 0.8996), and brain cancer (P = 0.8164). In the validation part of the study results, the causal relationship between dried fruit intake and lung cancer (P = 0.0043), squamous cell lung cancer (P = 0.0136), and breast cancer (P = 0.0192) was determined. After adjusting for the potential impact of confounders, the causal relationship between dried fruit intake and lung cancer (P = 0.0034), squamous cell lung cancer (P = 0.046), and breast cancer (P = 0.0001) remained. The sensitivity analysis showed that our results were stable and reliable.

Conclusion

The intake of dried fruits may have a protective effect against some site-specific cancers. Therefore, health education and a reasonable adjustment of dietary proportions may help in the primary prevention of cancer.

Susceptibility Genes and Chromosomal Regions Associated With Non-Syndromic Familial Non-Medullary Thyroid Carcinoma: Some Pathogenetic and Diagnostic Keys

Thyroid cancer is the malignant tumor that is increasing most rapidly in the world, mainly at the expense of sporadic papillary thyroid carcinoma. The somatic alterations involved in the pathogenesis of sporadic follicular cell derived tumors are well recognized, while the predisposing alterations implicated in hereditary follicular tumors are less well known. Since the genetic background of syndromic familial non-medullary carcinoma has been well established, here we review the pathogenesis of non-syndromic familial non-medullary carcinoma emphasizing those aspects that may be useful in clinical and pathological diagnosis. Non-syndromic familial non-medullary carcinoma has a complex and heterogeneous genetic basis involving several genes and loci with a monogenic or polygenic inheritance model. Most cases are papillary thyroid carcinoma (classic and follicular variant), usually accompanied by benign thyroid nodules (follicular thyroid adenoma and/or multinodular goiter). The possible diagnostic and prognostic usefulness of the changes in the expression and/or translocation of various proteins secondary to several mutations reported in this setting requires further confirmation. Given that non-syndromic familial non-medullary carcinoma and sporadic non-medullary thyroid carcinoma share the same morphology and somatic mutations, the same targeted therapies could be used at present, if necessary, until more specific targeted treatments become available.

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.

Pro64His (rs4644) Polymorphism Within Galectin-3 Is a Risk Factor of Differentiated Thyroid Carcinoma and Affects the Transcriptome of Thyrocytes Engineered via CRISPR/Cas9 System

Background: Galectin-3 (LGALS3) is an important glycoprotein involved in the malignant transformation of thyrocytes acting in the extracellular matrix, cytoplasm, and nucleus where it regulates TTF-1 and TCF4 transcription factors. Within LGALS3 gene, a common single-nucleotide polymorphism (SNP) (c.191C>A, p.Pro64His; rs4644) encoding for the variant Proline to Histidine at codon 64 has been extensively studied. However, data on rs4644 in the context of thyroid cancer are lacking. Thus, the aim of the present work was to evaluate the role of the rs4644 SNP as risk factor for differentiated thyroid cancer (DTC) and to determine the effect on the transcriptome in thyrocytes. Methods: A case/control association study in 1223 controls and 1142 unrelated consecutive DTC patients was carried out to evaluate the association between rs4644-P64H and the risk of DTC. We used the nonmalignant cell line Nthy-Ori (rs4644-C/A) and the CRISPR/Cas9 technique to generate isogenic cells carrying either the rs4644-A/A or rs4644-C/C homozygosis. Then, the transcriptome of the derivative and unmodified parental cells was analyzed by RNA-seq. Genes differentially expressed were validated by quantitative reverse transcription PCR and further tested in the parental Nthy-Ori cells after LGALS3 gene silencing, to investigate whether the expression of target genes was dependent on galectin-3 levels. Results: rs4644 AA genotype was associated with a reduced risk of DTC (compared with CC, OR_adj = 0.66; 95% confidence interval = 0.46-0.93; P_ass = 0.02). We found that rs4644 affects galectin-3 as a transcriptional coregulator. Among 34 genes affected by rs4644, HES1, HSPA6, SPC24, and NHS were of particular interest since their expression was rs4644-dependent (CC>AA for the first and AA>CC for the others), also in 574 thyroid tissues of Genotype-Tissue Expression (GTEx) biobank. Moreover, the expression of these genes was regulated by LGALS3-silencing. Using the proximity ligation assay in Nthy-Ori cells, we found that the TTF-1 interaction was genotype dependent. Conclusions: Our data show that in thyroid, rs4644 is a trans-expression quantitative trait locus that can modify the transcriptional expression of downstream genes, through the modulation of TTF-1.

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.

Mitochondrial ATP-Dependent Proteases-Biological Function and Potential Anti-Cancer Targets

Cells must eliminate excess or damaged proteins to maintain protein homeostasis. To ensure protein homeostasis in the cytoplasm, cells rely on the ubiquitin-proteasome system and autophagy. In the mitochondria, protein homeostasis is regulated by mitochondria proteases, including four core ATP-dependent proteases, m-AAA, i-AAA, LonP, and ClpXP, located in the mitochondrial membrane and matrix. This review will discuss the function of mitochondrial proteases, with a focus on ClpXP as a novel therapeutic target for the treatment of malignancy. ClpXP maintains the integrity of the mitochondrial respiratory chain and regulates metabolism by degrading damaged and misfolded mitochondrial proteins. Inhibiting ClpXP genetically or chemically impairs oxidative phosphorylation and is toxic to malignant cells with high ClpXP expression. Likewise, hyperactivating the protease leads to increased degradation of ClpXP substrates and kills cancer cells. Thus, targeting ClpXP through inhibition or hyperactivation may be novel approaches for patients with malignancy.

Association of Polygenetic Risk Scores Related to Cell Differentiation and Inflammation with Thyroid Cancer Risk and Genetic Interaction with Dietary Intake

Simple Summary

Global thyroid cancer incidence is increasing, especially in women. Genetic and environmental factors mutually contribute to its incidence. We aimed to identify genetic variants to influence thyroid cancer risk and determine their interactions with lifestyles in a large city hospital-based cohort (495 thyroid cancer patients and 56,439 control). The best polygenetic model included DIRC3_rs6759952, GAP43_rs13059137, NRG1_rs7834206, PROM1_rs72616195, LRP1B_rs1369535, and LOC100507065_rs11175834, tumorigenesis and cancer cell differentiation-related genes. Their high polygenetic risk scores (PRS) increased thyroid cancer risk by 3.90-fold compared to low-PRS. Thyroid cancer risk was elevated in females, high white blood cell counts, and high energy, low alcohol, and high seaweed intakes by 4.21, 4.03, 7.00, 4.11, and 4.02-fold, respectively. These factors interacted with PRS: the women with high-PRS elevated thyroid cancer risk much among women with high daily energy, seaweeds, and alcohol intake. These results could be applied to personalized nutrition plans to reduce thyroid cancer risk.

Abstract

The incidence of thyroid cancer continues to increase steadily, and this increasing incidence cannot be attributed solely to the overdiagnosis of microcarcinoma or technical advancements in detection methods and may also depend on environmental and genetic factors. However, the impacts and interactions of genetic and environmental factors remain controversial, and they may differ in Eastern and Western countries. The study’s purpose was to identify single nucleotide polymorphisms of genes related to cell differentiation and inflammation to influence thyroid cancer incidence and determine interactions with lifestyles in a large city hospital-based cohort. Genetic variants were selected by genome-wide association study with thyroid cancer participants (case; n = 495) and controls without cancers (n = 56,439). SNPs having gene–gene interactions were selected by generalized multifactor dimensionality reduction. Polygenic risk scores (PRSs) were generated by summing the number of selected SNP risk alleles. PRSs of the best model included 6 SNPs, that is, DIRC3_rs6759952, GAP43_rs13059137, NRG1_rs7834206, PROM1_rs72616195, LRP1B_rs1369535, and LOC100507065_rs11175834. Participants with a high-PRS had a higher thyroid cancer risk by 3.9-fold than those with a low-PRS. The following variables were related to an increased thyroid cancer risk; female (OR = 4.21), high white blood cell count (OR = 4.03), and high energy (OR = 7.00), low alcohol (OR = 4.11), and high seaweed (OR = 4.02) intakes. These variables also interacted with PRS to influence thyroid cancer risk. Meat/noodle diet patterns interacted with PRSs to increase thyroid cancer risk (p = 0.0023). In conclusion, women with a high-PRS associated with cell differentiation and inflammation were at an elevated thyroid cancer risk. Daily energy, seaweeds, and alcohol intake interacted with PRS for thyroid cancer risk. These results could be applied to personalized nutrition plans to reduce the risk of thyroid cancer.

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.

Inherited Thyroid Tumors With Oncocytic Change

Familial non-medullary thyroid carcinoma (FNMTC) corresponds to 5-10% of all follicular cell-derived carcinoma (FCDTC). Oncocytic thyroid tumors have an increased incidence in the familial context in comparison with sporadic FCDTC, encompassing benign and malignant tumors in the same family presenting with some extent of cell oxyphilia. This has triggered the interest of our and other groups to clarify the oncocytic change, looking for genetic markers that could explain the emergence of this phenotype in thyroid benign and malignant lesions, focusing on familial aggregation. Despite some advances regarding the identification of the gene associated with retinoic and interferon-induced mortality 19 (GRIM-19), as one of the key candidate genes affected in the "Tumor with Cell Oxyphilia" (TCO) locus, most of the mutations follow a pattern of "private mutations", almost exclusive to one family. Moreover, no causative genetic alterations were identified so far in most families. The incomplete penetrance of the disease, the diverse benign and malignant phenotypes in the affected familial members and the variable syndromic associations create an additional layer of complexity for studying the genetic alterations in oncocytic tumors. In the present review, we summarized the available evidence supporting genomic-based mechanisms for the oncocytic change, particularly in the context of FNMTC. We have also addressed the challenges and gaps in the aforementioned mechanisms, as well as molecular clues that can explain, at least partially, the phenotype of oncocytic tumors and the respective clinico-pathological behavior. Finally, we pointed to areas of further investigation in the field of oncocytic (F)NMTC with translational potential in terms of therapy.

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.

Characterizing the function of EPB41L4A in the predisposition to papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common histotype of thyroid carcinoma. The heritability of PTC is high compared to other cancers, but its underlying causes are unknown. A recent genome-wide association study revealed the association of a variant at the 5q22 locus, rs73227498, with PTC predisposition. We report that rs17134155, a variant in high linkage disequilibrium with rs73227498, is located in an enhancer region downstream of coding transcripts of EPB41L4A. Rs17134155 showed significant enhancer activity in luciferase assays, and haplotypes containing the protective allele of this variant conferred a significantly lower risk of PTC. While the index SNP, rs73227498, acted as a significant cis-eQTL for expression of EPB41L4A, rs17134155 was a significant cis-sQTL for the alternative splicing of a non-coding transcript of EPB41L4A, called EPB41L4A-203. We also performed knockdown of EPB41L4A followed by microarray analysis. Some of the top differentially-expressed genes were represented among regulators of the WNT/β-catenin signaling pathway. Our results indicate that an enhancer region at 5q22 regulates the expression and splicing of EPB41L4A transcripts. We also provide evidence that EPB41L4A expression is involved in regulating growth and differentiation pathways, suggesting that decreased expression of EPB41L4A is a mechanism in the predisposition to PTC.

IL-1RN gene polymorphisms reduces thyroid cancer risk in Chinese Han population

The purpose of this study was to investigate the effect of IL-1RN polymorphisms on thyroid cancer (TC) risk in Han population. Genotypes of four single nucleotide polymorphisms (SNPs) (rs17042888, rs928940, rs3181052, and rs452204) were analyzed by Agena MassARRAY. Meanwhile, we used the logistic regression to calculate the odds ratios (ORs) and 95% confidence intervals (CIs), and significant differences were evaluated by t test and χ² test. Findings found that allele "G" of rs452204 and rs3181052 in interleukin-1 receptor antagonist (IL-1RN) reduced the risk of TC. (OR = 0.72, 95% CI = 0.55-0.94, p = .017; OR = 0.73, 95% CI = 0.56-0.94, p = .017, respectively). Hierarchical analysis indicated that three SNPs (rs17042888, rs3181052, and rs452204) significantly reduced the risk of TC among females or individuals older than 48 years (p < .05). Our findings indicate that IL-1RN polymorphisms may contribute to a protective role against TC risk.

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.

Linkage disequilibrium and haplotype analysis of Src and Yes1 genes in thyroid cancer

Purpose: This study was planned to examine the effects of Src and Yes1 single nucleotide polymorphism (SNPs) on the risk of thyroid cancer in 499 patients and 500 controls. Materials & methods: Three SNPs of Src gene and three SNPs of Yes1 gene were analyzed using Tetra-primer ARMS-PCR followed by sequencing. Results: rs121913314 of Src gene genotype TT showed 32-fold increased risk of thyroid cancer and rs2305994 of Yes1 genotypes TT and CT showed 2.7-fold and 16-fold increased risk in thyroid cancer (p < 0.0001). Haplotype analysis revealed that CATGCC, CATGCT, CATGTC, CATGTT, TATGCC and TATGTTA haplotypes are associated with thyroid cancer risk. Conclusion: Results showed that genotypes and allele distribution of Src and Yes1 genes are significantly linked with increased risk of thyroid cancer.

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.

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.

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.

VDR gene FokI polymorphism as a poor prognostic factor for papillary thyroid cancer

This is the first study to investigate the effect of vitamin D receptor ( VDR) gene single-nucleotide polymorphisms on the clinicopathologic features of papillary thyroid cancer in Turkey. A total of 165 patients with papillary thyroid cancer and 172 controls were included in this case-control study. VDR gene single-nucleotide polymorphisms FokI (rs2228570), BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) were evaluated using reverse-transcription polymerase chain reaction. VDR gene polymorphisms BsmI, ApaI, and TaqI did not differ between the papillary thyroid cancer group and control group (p > 0.05, each). BsmI, ApaI, and TaqI were not associated with papillary thyroid cancer risk. The VDR gene FokI CT/TT genotype was associated with an increased papillary thyroid cancer risk (CT vs CC: odds ratio = 1.71, 95% confidence interval = 1.15-2.76, p = 0.028; TT vs CC: odds ratio = 2.44, 95% confidence interval = 1.29-4.62, p = 0.005; CT/TT vs CC: odds ratio = 1.88, 95% confidence interval = 1.20-2.96, p = 0.006; CT/CC vs TT: odds ratio = 1.80, 95% confidence interval = 1.05-3.20, p = 0.041). VDR gene polymorphisms were not in linkage disequilibrium. The FokI TT genotype was associated with having T3 and T4, stage III/IV, extra-thyroidal invasion. The FokI CT/TT or TT genotype was associated with developing N1 status, multifocality, tumor size ≥10 mm, and treatment with radioiodine therapy. Persistence/recurrence did not differ between the FokI genotypes. Carriers of the FokI T allele were at an increased risk of more advanced tumor-node-metastasis stage, greater tumor size, multifocality, and extra-thyroidal invasion of papillary thyroid cancer compared with the CC genotype. VDR gene FokI T allele and TT genotype correlated with aggressiveness of papillary thyroid cancer; thus, FokI could be useful as a poor prognostic factor to assess the high risk of papillary thyroid cancer.

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.

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.

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.