Role of DNA Repair Variants and Diagnostic Radiology Exams in Differentiated Thyroid Cancer Risk: A Pooled Analysis of Two Case-Control Studies.

BACKGROUND: Given the increased use and diversity of diagnostic procedures, it is important to understand genetic susceptibility to radiation-induced thyroid cancer. METHODS: On the basis of self-declared diagnostic radiology examination records in addition to existing literature, we estimated the radiation dose delivered to the thyroid gland from diagnostic procedures during childhood and adulthood in two case-control studies conducted in France. A total of 1,071 differentiated thyroid cancer (DTC) cases and 1,188 controls from the combined studies were genotyped using a custom-made Illumina OncoArray DNA chip. We focused our analysis on variants in genes involved in DNA damage response and repair pathways, representing a total of 5,817 SNPs in 571 genes. We estimated the OR per milli-Gray (OR/mGy) of the radiation dose delivered to the thyroid gland using conditional logistic regression. We then used an unconditional logistic regression model to assess the association between DNA repair gene variants and DTC risk. We performed a meta-analysis of the two studies. RESULTS: The OR/mGy was 1.02 (95% confidence interval, 1.00-1.03). We found significant associations between DTC and rs7164173 in CHD2 (P = 5.79 × 10-5), rs6067822 in NFATc2 (P = 9.26 × 10-5), rs1059394 and rs699517 both in ENOSF1/THYS, rs12702628 in RPA3, and an interaction between rs7068306 in MGMT and thyroid radiation doses (P = 3.40 × 10-4). CONCLUSIONS: Our results suggest a role for variants in CDH2, NFATc2, ENOSF1/THYS, RPA3, and MGMT in DTC risk. IMPACT: CDH2, NFATc2, ENOSF1/THYS, and RPA3 have not previously been shown to be associated with DTC risk.

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

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

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

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

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.