Quantitative assessment of common genetic variants on FOXE1 and differentiated thyroid cancer risk

Transcription Factors: The Fulcrum Between Cell Development and Carcinogenesis

Higher eukaryotic development is a complex and tightly regulated process, whereby transcription factors (TFs) play a key role in controlling the gene regulatory networks. Dysregulation of these regulatory networks has also been associated with carcinogenesis. Transcription factors are key enablers of cancer stemness, which support the maintenance and function of cancer stem cells that are believed to act as seeds for cancer initiation, progression and metastasis, and treatment resistance. One key area of research is to understand how these factors interact and collaborate to define cellular fate during embryogenesis as well as during tumor development. This review focuses on understanding the role of TFs in cell development and cancer. The molecular mechanisms of cell fate decision are of key importance in efforts towards developing better protocols for directed differentiation of cells in research and medicine. We also discuss the dysregulation of TFs and their role in cancer progression and metastasis, exploring TF networks as direct or indirect targets for therapeutic intervention, as well as specific TFs' potential as biomarkers for predicting and monitoring treatment responses.

The Dominant Role of Forkhead Box Proteins in Cancer

Forkhead box (FOX) proteins are multifaceted transcription factors that are significantly implicated in cancer, with various critical roles in biological processes. Herein, we provide an overview of several key members of the FOXA, FOXC, FOXM1, FOXO and FOXP subfamilies. Important pathophysiological processes of FOX transcription factors at multiple levels in a context-dependent manner are discussed. We also specifically summarize some major aspects of FOX transcription factors in association with cancer research such as drug resistance, tumor growth, genomic alterations or drivers of initiation. Finally, we suggest that targeting FOX proteins may be a potential therapeutic strategy to combat 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.

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

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

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

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

Association of FOXE1 polyalanine repeat region with thyroid cancer is dependent on tumour size

OBJECTIVE

Polymorphisms in the thyroid transcription factor forkhead factor E1 (FOXE1) gene have been implicated in the genetic susceptibility to differentiated thyroid cancer, but little is known about their effect on tumour characteristics. The objective of this study was to determine the contribution of the FOXE1 polyalanine repeat region to the susceptibility to thyroid cancer and to its clinical characteristics.

DESIGN, PATIENTS AND MEASUREMENTS

A total of 500 patients with sporadic thyroid cancer (440 papillary and 60 follicular thyroid carcinoma) and 502 healthy controls were included in this case-control association study. The number of FOXE1 alanine repeats in each subject was determined by PCR and multiplex fragment analysis by capillary electrophoresis. FOXE1 genotype and allele frequencies among groups were compared by logistic regression and adjusted for sex and age at diagnosis. Data were analysed according to cancer subtype, tumour size and the presence of lymph node or distant metastasis.

RESULTS

FOXE1 alleles with 16 or more alanine repeats were more frequent in patients with tumour size > 1 cm compared to tumour size ≤ 1 cm (adjusted OR 1·44; 95% CI 1·05-1·88; P = 0·019). Genotypes containing at least one allele with 16 or more alanine repeats were associated with larger tumour size (adjusted OR 1·71; 95% CI 1·15-2·57; P = 0·009). No significant differences were observed between cancer subtypes or the presence/absence of metastasis.

CONCLUSIONS

FOXE1 polyalanine repeat polymorphisms are associated with thyroid cancer, but only for tumours larger than 1 cm, suggesting a role in disease progression.

Genotype Analyses in the Japanese and Belarusian Populations Reveal Independent Effects of rs965513 and rs1867277 but Do Not Support the Role of FOXE1 Polyalanine Tract Length in Conferring Risk for Papillary Thyroid Carcinoma

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Spatiotemporal Co-existence of Female Thyroid and Breast Cancers in Hangzhou, China

Thyroid and breast cancers (TC, BC) are common female malignant tumors worldwide. Studies suggest that TC patients have a higher BC risk, and vice versa. However, it has not been investigated quantitatively if there is an association between the space-time TC and BC incidence distributions at the population level. This work aims to answer this question. 5358 TC and 8784 BC (female) cases were diagnosed in Hangzhou (China, 2008-2012). Pearson and Spearman rank correlation coefficients of the TC and BC incidences were high, and their patterns were geographically similar. The spatiotemporal co-existence of TC and BC distributions was investigated using the integrative disease predictability (IDP) criterion: if TC-BC association is part of the disease mapping knowledge bases, it should yield improved space-time incidence predictions. Improved TC (BC) incidence predictions were generated when integrating both TC and BC data than when using only TC (BC) data. IDP consistently demonstrated the spatiotemporal co-existence of TC and BC distributions throughout Hangzhou (2008-2012), which means that when the population experiences high incidences of one kind of cancer attention should be paid to the other kind of cancer too. The strength of TC-BC association was measured by the IDP coefficients and incidence prediction accuracy.

Multiple functional variants in long-range enhancer elements contribute to the risk of SNP rs965513 in thyroid cancer

The [A] allele of SNP rs965513 in 9q22 has been consistently shown to be highly associated with increased papillary thyroid cancer (PTC) risk with an odds ratio of ∼1.8 as determined by genome-wide association studies, yet the molecular mechanisms remain poorly understood. Previously, we noted that the expression of two genes in the region, forkhead box E1 (FOXE1) and PTC susceptibility candidate 2 (PTCSC2), is regulated by rs965513 in unaffected thyroid tissue, but the underlying mechanisms were not elucidated. Here, we fine-mapped the 9q22 region in PTC and controls and detected an ∼33-kb linkage disequilibrium block (containing the lead SNP rs965513) that significantly associates with PTC risk. Chromatin characteristics and regulatory element signatures in this block disclosed at least three regulatory elements functioning as enhancers. These enhancers harbor at least four SNPs (rs7864322, rs12352658, rs7847449, and rs10759944) that serve as functional variants. The variant genotypes are associated with differential enhancer activities and/or transcription factor binding activities. Using the chromosome conformation capture methodology, long-range looping interactions of these elements with the promoter region shared by FOXE1 and PTCSC2 in a human papillary thyroid carcinoma cell line (KTC-1) and unaffected thyroid tissue were found. Our results suggest that multiple variants coinherited with the lead SNP and located in long-range enhancers are involved in the transcriptional regulation of FOXE1 and PTCSC2 expression. These results explain the mechanism by which the risk allele of rs965513 predisposes to thyroid cancer.

Double-negative feedback loop between microRNA-422a and forkhead box (FOX)G1/Q1/E1 regulates hepatocellular carcinoma tumor growth and metastasis

Growing evidence indicates that the aberrant expression of microRNAs (miRNAs) contributes to tumor development; however, the function of miRNAs in human hepatocellular carcinoma (HCC) remains largely undefined. In this study, we report that microRNA-422a (miR-422a) is significantly down-regulated in HCC tumor samples and cell lines compared with normal controls, and its expression level is negatively correlated with pathological grading, recurrence, and metastasis. The restoration of miR-422a expression in HCC tumor cells significantly inhibited cell proliferation and migration in vitro. At the same time, the overexpression of miR-422a in HCC tumor cells significantly inhibits tumor growth and liver metastasis in xenograft tumor models. A mechanistic study identified three genes, forkhead box G1 (FOXG1), FOXQ1, and FOXE1, as miR-422a targets in the regulation of HCC development. We also investigated the function of the three targets themselves in HCC tumorigenesis using RNAi manipulation and demonstrated that the knockdown of these targets led to significant inhibition of tumor cell proliferation and migration both in vitro and in vivo. More interestingly, a potential miR-422a promoter region was identified. Both the promoter activity and miR-422a expression were negatively regulated by the three targets, indicating that a double-negative feedback loop exists between miR-422a and its targets. Moreover, we explored the therapeutic potential of miR-422a in HCC treatment and found that the therapeutic delivery of miR-422a significantly inhibited tumor development in a xenograft tumor model and a diethylnitrosamine-induced primary HCC model.

CONCLUSION

Our findings show the critical roles of miR-422a and its targets--FOXG1, FOXQ1, and FOXE1--in the regulation of HCC development and provide new potential candidates for HCC therapy.

Genetic predisposition to papillary thyroid carcinoma: involvement of FOXE1, TSHR, and a novel lincRNA gene, PTCSC2

CONTEXT

By genome-wide association studies, the risk allele [A] of SNP rs965513 predisposes strongly to papillary thyroid carcinoma (PTC). It is located in a gene-poor region of 9q22, some 60 kb from the FOXE1 gene. The underlying mechanisms remain to be discovered.

OBJECTIVE

Our objective was to identify novel transcripts in the 9q22 locus and correlate gene expression levels with the genotypes of rs965513.

DESIGN

We performed 3' and 5' rapid amplification of cDNA ends and RT-PCR to detect novel transcripts. One novel transcript was forcibly expressed in a cell line followed by gene expression array analysis. We genotyped rs965513 from PTC patients and measured gene expression levels by real-time RT-PCR in unaffected thyroid tissue and matched tumor.

SETTING

This was a laboratory-based study using cells from clinical tissue samples and a cancer cell line.

MAIN OUTCOME MEASURES

We detected previously uncharacterized transcripts and evaluated the gene expression levels and the correlation with the risk allele of rs965513, age, gender, chronic lymphocyte thyroiditis (CLT), and TSH levels.

RESULTS

We found a novel long intergenic noncoding RNA gene and named it papillary thyroid cancer susceptibility candidate 2 (PTCSC2). Transcripts of PTCSC2 are down-regulated in PTC tumors. The risk allele [A] of rs965513 was significantly associated with low expression of unspliced PTCSC2, FOXE1, and TSHR in unaffected thyroid tissue. We also observed a significant association of age and CLT with PTCSC2 unspliced transcript levels. The correlation between the rs965513 genotype and the PTCSC2 unspliced transcript levels remained significant after adjusting for age, gender, and CLT. Forced expression of PTCSC2 in the BCPAP cell line affected the expression of a subset of noncoding and coding transcripts with enrichment of genes functionally involved in cell cycle and cancer.

CONCLUSIONS

Our data suggest a role for PTCSC2, FOXE1, and TSHR in the predisposition to PTC.