Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 9q22.33 and 14q13.3 detects novel candidate functional SNPs in Europeans from metropolitan France and Melanesians from New Caledonia

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

CONTEXT

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Effect of Cigarette Smoking on Thyroid Cancer: Meta-Analysis

BACKGROUND

Although smoking is generally carcinogenic, its effect on thyroid cancers is still subject to controversy. The purpose of this study was to summarize the role of smoking in relation to thyroid cancer occurrence.

METHODS

We performed a meta-analysis of 24 eligible studies: 21 case-control studies and three prospective cohort studies. The summary odds ratio (OR) and 95% confidence interval (CI) of all studies were acquired based on random effect model. Further subgroup analyses were conducted according to gender, histological type of thyroid cancer, and smoking status of patients for the case-control studies.

RESULTS

The summary effect size indicated a negative association of smoking for thyroid cancer (OR, 0.798; 95% CI, 0.681 to 0.935). From the subgroup analyses for the case-control studies, reduced risk of thyroid cancer was observed in both men (OR, 0.734; 95% CI, 0.553 to 0.974) and women (OR, 0.792; 95% CI, 0.700 to 0.897). The protective effect of smoking was observed in studies in which thyroid cancer was limited to differentiated thyroid cancers (DTCs) (OR, 0.798; 95% CI, 0.706 to 0.902).

CONCLUSION

Our results suggests that smoking may have a protective effect on thyroid cancer, especially on DTCs. Further studies with larger sample sizes should be conducted in elucidating the dose and time dependent effect of smoking on thyroid cancer with specific focus on the types of thyroid cancers.

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.

Dietary Inflammatory Index and Differentiated Thyroid Carcinoma Risk: A Population-Based Case-Control Study in New Caledonia

Despite research efforts, current knowledge of the etiology of thyroid carcinoma remains limited. To explore the potential role of diet-induced inflammation, we examined the association between differentiated thyroid cancer risk and the energy-adjusted Dietary Inflammatory Index (E-DII) in a population-based case-control study conducted in New Caledonia, a Pacific archipelago with one of the highest recorded thyroid cancer incidence rates in the world. The E-DII was computed from food frequency questionnaire information on usual dietary intake. Logistic regression analyses were performed on data from 324 histologically confirmed cases of papillary or follicular carcinoma, diagnosed from 1993 to 1999, and 402 controls. Positive associations between E-DII and thyroid cancer risk were observed (comparing extreme tertiles, odds ratio = 1.67, 95% confidence interval: 1.08, 2.58; P for trend = 0.002), with stronger associations found for larger carcinomas (P for trend = 0.0005). Stratified analyses showed an increased risk of thyroid cancer associated with the E-DII among Southern province residents (P for trend = 0.003), Melanesian women (P for trend = 0.02), obese participants (P for trend = 0.006), and ever-smokers (P for trend = 0.0005). Our results suggest that a proinflammatory diet-especially when concomitant with other inflammation-inducing conditions or habits (e.g., obesity, smoking)-is associated with increased risk of thyroid carcinoma.

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

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

Association between SNPs in Long Non-coding RNAs and the Risk of Female Breast Cancer in a Chinese Population

Long non-coding RNAs (LncRNAs) have been reported to be involved in tumorigenesis and tumor progression. Single nucleotide polymorphisms (SNPs) in the lncRNAs also play a vital role in carcinogenesis. The aim of this study was to assess the relationships between the four selected tagSNPs (rs944289, rs3787016, rs1456315, rs7463708) in the lncRNAs and the risk of female breast cancer in a Chinese population. A case-control study was carried out involving in a total of 439 breast cancer patients and 439 age-matched healthy controls. The genotyping was performed with Sequenom MassARRAY and the expression of estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor-2 (HER-2) in tumor tissues was measured by the immunohistochemistry (IHC) assay. We found that rs3787016 TT genotype (adjusted odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.09-2.41, P = 0.018) was associated with an increased risk of female breast cancer, especially among the patients with premenopausal status (adjusted OR = 2.55, 95% CI = 1.30-4.97, P = 0.006). Moreover, a statistically significant increased risk of the rs3787016 TT genotype was observed among the patients with advanced tumor stage (Ⅲ and Ⅳ), poor histological grade (G3-G4), positive lymph node involvement, positive expression of ER and PR and negative expression of HER-2; rs7463708 GT and GT/GG genotype were associated with decreased risk of breast cancer in the subgroup of patients with postmenopausal status (GT versus (vs.) TT: adjusted OR = 0.67, 95% CI = 0.46-0.99, P = 0.043; GT/GG vs. TT: adjusted OR = 0.68, 95% CI = 0.47-0.98, P = 0.041) and tumor late-stage (GT vs. TT: adjusted OR = 0.65, 95% CI = 0.43-0.97, P = 0.037; GT/GG vs. TT: adjusted OR = 0.65, 95% CI = 0.44-0.96, P = 0.029). In short, rs3787016 TT genotype was associated with increased breast cancer risk and clinicopathologic features of the tumor, especially among premenopausal women.

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

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

Interaction Landscape of Inherited Polymorphisms with Somatic Events in Cancer

Recent studies have characterized the extensive somatic alterations that arise during cancer. However, the somatic evolution of a tumor may be significantly affected by inherited polymorphisms carried in the germline. Here, we analyze genomic data for 5,954 tumors to reveal and systematically validate 412 genetic interactions between germline polymorphisms and major somatic events, including tumor formation in specific tissues and alteration of specific cancer genes. Among germline-somatic interactions, we found germline variants in RBFOX1 that increased incidence of SF3B1 somatic mutation by 8-fold via functional alterations in RNA splicing. Similarly, 19p13.3 variants were associated with a 4-fold increased likelihood of somatic mutations in PTEN. In support of this association, we found that PTEN knockdown sensitizes the MTOR pathway to high expression of the 19p13.3 gene GNA11 Finally, we observed that stratifying patients by germline polymorphisms exposed distinct somatic mutation landscapes, implicating new cancer genes. This study creates a validated resource of inherited variants that govern where and how cancer develops, opening avenues for prevention research.Significance: This study systematically identifies germline variants that directly affect tumor evolution, either by dramatically increasing alteration frequency of specific cancer genes or by influencing the site where a tumor develops. Cancer Discovery; 7(4); 410-23. ©2017 AACR.See related commentary by Geeleher and Huang, p. 354This article is highlighted in the In This Issue feature, p. 339.