Thyroid transcription factor FOXE1 interacts with ETS factor ELK1 to co-regulate TERT

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

InTERTwined: how TERT promoter mutations impact BRAFV600E-driven thyroid cancers

Thyroid cancers are often initiated by the acquisition of a BRAFV600E mutation. BRAFV600E-driven thyroid tumors display a wide range of behaviors, from the slow-growing papillary carcinomas to the highly aggressive anaplastic. Mutations in the promoter of TERT (telomerase reverse transcriptase) gene were discovered a decade ago and identified as prevalent events in thyroid cancers. Multiple studies showed that TERT promoter mutations, particularly when co-occurring with BRAFV600E, are markers of poor prognosis across thyroid cancer subtypes, and can be implemented for routine clinical stratification. Mechanistically, TERT promoter mutations reactivate telomerase expression via the differential recruitment of transcriptional complexes. Re-expression of TERT impacts tumor biology, plausibly via both the well-known function of telomerase maintaining telomeres and by affecting other cancer-relevant processes.

Multiple ETS Factors Participate in the Transcriptional Control of TERT Mutant Promoter in Thyroid Cancers

Simple Summary

Mutations in the promoter region of the telomerase reverse transcriptase (TERT) gene are enriched in patients with advanced thyroid tumors. Their consequence is the reactivation of TERT expression through mechanisms that involve specific transcription factors and other signaling inputs. Here, we show that, contrary to what it has been shown in other tumor types, multiple factors of the ETS family are able to control TERT transcription and they do so in both the presence and absence of promoter mutations. We also show that TERT expression is more dependent on the MAPK signaling pathway in thyroid cells without TERT promoter mutations. Our work points to an intricate, still not fully characterized, regulatory network of TERT transcription in thyroid tumors. We caution against the assumption of mechanisms identified in other cancer lineages being identical in TERT-mutant thyroid specimens.

Abstract

Hotspot mutations in the TERT (telomerase reverse transcriptase) gene are key determinants of thyroid cancer progression. TERT promoter mutations (TPM) create de novo consensus binding sites for the ETS (“E26 transformation specific”) family of transcription factors. In this study, we systematically knocked down each of the 20 ETS factors expressed in thyroid tumors and screened their effects on TERT expression in seven thyroid cancer cell lines with defined TPM status. We observed that, unlike in other TPM-carrying cancers such as glioblastomas, ETS factor GABPA does not unambiguously regulate transcription from the TERT mutant promoter in thyroid specimens. In fact, multiple members of the ETS family impact TERT expression, and they typically do so in a mutation-independent manner. In addition, we observe that partial inhibition of MAPK, a central pathway in thyroid cancer transformation, is more effective at suppressing TERT transcription in the absence of TPMs. Taken together, our results show a more complex scenario of TERT regulation in thyroid cancers compared with other lineages and suggest that compensatory mechanisms by ETS and other regulators likely exist and advocate for the need for a more comprehensive understanding of the mechanisms of TERT deregulation in thyroid tumors before eventually exploring TPM-specific therapeutic strategies.

A novel circular RNA circPPFIA1 promotes laryngeal squamous cell carcinoma progression through sponging miR-340-3p and regulating ELK1 expression

Abnormal expression of circular RNA (circRNA) is closely related to the occurrence and development of many cancers. By screening the expression of circRNA, we identified a novel circRNA termed as has_circ_0023326 in laryngeal squamous cell carcinoma (LSCC). We verified the expression of circPPFIA1 and found that it was upregulated in LSCC tissues compared to the adjacent normal tissues. Functional studies were carried out to detect the effect of circPPFIA1 expression on the phenotype of LSCC cells. These results suggest that circPPFIA1 knockdown can suppress the proliferation, migration, and invasion of LSCC cells, while circPPFIA1 overexpression can promote these processes. Mechanistically, miR-340-3p was predicted to be the target miRNA sponged by circPPFIA1 as confirmed through the luciferase assay and rescue experiments. In addition, miR-340-3p was found to target ELK1 and inhibit its expression. Taken together, circPPFIA1 promotes the progression of LSCC via the miR-340-3p/ELK1 signaling axis, which may serve as a novel prognostic or therapeutic target for LSCC.

ETS like-1 protein ELK1-induced lncRNA LINC01638 accelerates the progression of papillary thyroid cancer by regulating Axin2 through Wnt/β-catenin signaling pathway

Papillary thyroid carcinoma (PTC) characterized by distant metastasis is a major public health issue among women worldwide. LncRNA LINC01638 is reportedly a critical oncogene in the development of certain cancers. However, the biological function of LINC01638 in PTC is currently unclear. The goal of this study was to identify LINC01638 expression level and its role in PTC progression. The expression of LINC01638 was detected applying qRT-PCR. CCK-8 assay, colony formation assay, immunofluorescence staining and flow cytometric analysis were performed to assess cell proliferation and cell cycle. In addition, cell migration and invasion were examined via wound healing assay, transwell assay and western blot analysis. We found that LINC0163 was upregulated in PTC cells compared with normal thyroid gland epithelial cell line Nthy-ori3-1. ELK1 could act as a transcription factor of LINC01638 and induce LINC01638 expression. LINC01638 silencing inhibited cell proliferation, migration and invasion, and obstructed the progress of TPC-1 cell cycle. LINC0163 silencing activated Axin2 while suppressing the expressions of β-catenin, Cyclin-D1 and c-MYC. Rescue experiment utilizing the transfection of Axin2 overexpression plasmid weakened LINC01638 overexpression-enhanced TPC-1 cell proliferation, metastasis, cell cycle progress and Wnt/β-catenin pathway. These results indicate that LINC0163 regulates PTC progression via inhibition of Wnt/β-catenin and activation of Axin2, which may develop into a novel therapeutic strategy for PTC treatment.

The ETS Inhibitor YK-4-279 Suppresses Thyroid Cancer Progression Independent of TERT Promoter Mutations

Hotspot mutations in the core promoter region of the telomerase reverse transcriptase (TERT) gene have been well established to associate with aggressive clinical characteristics, radioiodine refractory, tumor recurrence, and mortality in thyroid cancer. Several E-twenty-six (ETS) transcription factors were reported to selectively bound to the mutant TERT promoter and activated TERT expression. In this study we aimed to investigate whether TERT promoter mutations confer sensitivity to ETS inhibitor YK-4-279 in thyroid cancer cells and whether this inhibitor could be served as a potential therapeutic agent for thyroid cancer. In vitro assays showed that YK-4-279 treatment sharply suppressed cell viability, colony formation, migration, and invasion, as well as induced cell cycle arrest and apoptosis in a panel of thyroid cancer cells. The cell viability after YK-4-279 treatment was similar between cell lines harboring mutant and wild-type TERT promoters. Furthermore, YK-4-279 treatment reduced both luciferase activity and mRNA expression of TERT independent of TERT promoter mutation status. Data from RNA-seq further revealed that YK-4-279 significantly affected biological processes including DNA replication and cell cycle. Reduced DNA helicase activity and decreased expression of several helicase genes were observed after YK-4-279 treatment. Moreover, YK-4-279 significantly inhibited tumor growth and induced apoptosis in a xenograft mice model. Thus, ETS inhibitor YK-4-279 suppressed TERT expression and conferred anti-tumor activity in a TERT promoter mutation-independent manner, and it could be a potential agent for the treatment of advanced thyroid cancers.

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

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

Network-Based Genetic Profiling Reveals Cellular Pathway Differences Between Follicular Thyroid Carcinoma and Follicular Thyroid Adenoma

Molecular mechanisms underlying the pathogenesis and progression of malignant thyroid cancers, such as follicular thyroid carcinomas (FTCs), and how these differ from benign thyroid lesions, are poorly understood. In this study, we employed network-based integrative analyses of FTC and benign follicular thyroid adenoma (FTA) lesion transcriptomes to identify key genes and pathways that differ between them. We first analysed a microarray gene expression dataset (Gene Expression Omnibus GSE82208, n = 52) obtained from FTC and FTA tissues to identify differentially expressed genes (DEGs). Pathway analyses of these DEGs were then performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) resources to identify potentially important pathways, and protein-protein interactions (PPIs) were examined to identify pathway hub genes. Our data analysis identified 598 DEGs, 133 genes with higher and 465 genes with lower expression in FTCs. We identified four significant pathways (one carbon pool by folate, p53 signalling, progesterone-mediated oocyte maturation signalling, and cell cycle pathways) connected to DEGs with high FTC expression; eight pathways were connected to DEGs with lower relative FTC expression. Ten GO groups were significantly connected with FTC-high expression DEGs and 80 with low-FTC expression DEGs. PPI analysis then identified 12 potential hub genes based on degree and betweenness centrality; namely, TOP2A, JUN, EGFR, CDK1, FOS, CDKN3, EZH2, TYMS, PBK, CDH1, UBE2C, and CCNB2. Moreover, transcription factors (TFs) were identified that may underlie gene expression differences observed between FTC and FTA, including FOXC1, GATA2, YY1, FOXL1, E2F1, NFIC, SRF, TFAP2A, HINFP, and CREB1. We also identified microRNA (miRNAs) that may also affect transcript levels of DEGs; these included hsa-mir-335-5p, -26b-5p, -124-3p, -16-5p, -192-5p, -1-3p, -17-5p, -92a-3p, -215-5p, and -20a-5p. Thus, our study identified DEGs, molecular pathways, TFs, and miRNAs that reflect molecular mechanisms that differ between FTC and benign FTA. Given the general similarities of these lesions and common tissue origin, some of these differences may reflect malignant progression potential, and include useful candidate biomarkers for FTC and identifying factors important for FTC pathogenesis.

ETS Factor ETV5 Activates the Mutant Telomerase Reverse Transcriptase Promoter in Thyroid Cancer

Background: Co-occurrence of TERT (telomerase reverse transcriptase) promoter (TERTp) mutations with BRAF/RAS mutations is associated with significantly more aggressive thyroid cancer. TERTp mutations are hypothesized to generate de novo binding sites for ETS transcription factors, which are themselves activated by BRAF/RAS-stimulated MEK-ERK activity. To date, a detailed study of this mechanism has been limited to only a few cancer types, and we hypothesized that ETS factors involved in TERTp activation could vary between different cancers. Methodology: Here we sought to identify ETS factor(s) required for TERTp activation in thyroid cancer, using a combination of in silico analyses of TCGA data, and experimentation using in vitro thyroid cell models analyzed by quantitative reverse transcription-PCR, immunoprecipitation (IP), chromatin IP, and gene reporter assays. Results: We found that ETV5 was abundantly expressed in papillary thyroid cancers from the TCGA data set, and in thyroid cancer cell line models. Furthermore, ETV5 was found to preferentially bind to the -124 bp(T) TERTp allele and stimulate TERT transcription in thyroid cancer cells devoid of GA binding protein transcription factor (GABP) activity. We also found that ETV5 functionally cooperates with the transcription factor FOXE1 to further enhance TERTp activity, a mechanism that may at least partially explain why FOXE1 represents a significant genetic determinant of thyroid cancer risk. Conclusions: ETS factors that activate mutant TERTp vary between cancer types, and here we show for the first time that ETV5 demonstrates mutant allele-specific affinity for TERTp in thyroid cancer, a property that has previously only been attributable to GABP.

Molecular characterization of papillary thyroid carcinoma: a potential three-lncRNA prognostic signature

Purpose

Papillary thyroid carcinoma (PTC), the most frequent type of malignant thyroid tumor, lacks novel and reliable biomarkers of patients’ prognosis. In the current study, we mined The Cancer Genome Atlas (TCGA) to develop lncRNA signature of PTC.

Patients and methods

The intersection of PTC lncRNAs was obtained from the TCGA database using integrative computational method. By the univariate and multivariate Cox analysis, key lncRNAs were identified to construct the prognostic model. Then, all patients were divided into the high-risk group and low-risk group to perform the Kaplan–Meier (K–M) survival curves and time-dependent receiver operating characteristic (ROC) curve, estimating the prognostic power of the prognostic model. Functional enrichment analysis was also performed. Finally, we verified the results of the TCGA analysis by the Gene Expression Omnibus (GEO) databases and quantitative real-time PCR (qRT-PCR).

Results

After the comprehensive analysis, a three-lncRNA signature (PRSS3P2, KRTAP5-AS1 and PWAR5) was obtained. Interestingly, patients with low-risk scores tended to gain obviously longer survival time, and the area under the time-dependent ROC curve was 0.739. Furthermore, gene ontology (GO) and pathway analysis revealed the tumorigenic and prognostic function of the three lncRNAs. We also found three potential transcription factors to help understand the mechanisms of the PTC-specific lncRNAs. Finally, the GEO databases and qRT-PCR validation were consistent with our TCGA bioinformatics results.

Conclusion

We built a three-lncRNA signature by mining the TCGA database, which could effectively predict the prognosis of PTC.

Does the Polymorphism in the Length of the Polyalanine Tract of FOXE1 Gene Influence the Risk of Thyroid Dysgenesis Occurrence?

Background. Recent data have suggested that polymorphisms in the length of the polyalanine tract (polyA) of FOXE1 gene may act as a susceptibility factor for thyroid dysgenesis. The main purpose of this study was to investigate the influence of polyA of FOXE1 gene on the risk of thyroid dysgenesis. Method. A case-control study was conducted in a sample of 90 Brazilian patients with thyroid dysgenesis and 131 controls without family history of thyroid disease. Genomic DNA was isolated from peripheral blood samples and the genotype of each individual was determined by automated sequencing. Results. More than 90% of genotypes found in the group of patients with thyroid dysgenesis and in controls subjects were represented by sizes 14 and 16 polymorphisms in the following combinations: 14/14, 14/16, and 16/16. Genotypes 14/16 and 16/16 were more frequent in the control group, while genotype 14/14 was more frequent in the group of patients with thyroid dysgenesis. There was no difference between agenesis group and control group. Genotype 14/14 when compared to genotypes 14/16 and 16/16A showed an association with thyroid dysgenesis. Conclusion. PolyA of FOXE1 gene alters the risk of thyroid dysgenesis, which may explain in part the etiology of this disease.

RNA-sequencing investigation identifies an effective risk score generated by three novel lncRNAs for the survival of papillary thyroid cancer patients

Scholars are striving to apply molecular biology involving long non-coding RNA (lncRNA) in the prognostication of papillary thyroid cancer (PTC). However, the clinical role of lncRNAs in the prognostic setting of PTC is still unclear. Herein, a comprehensive inquiry was performed to screen lncRNA expression profiling with 507 PTC patients from The Cancer Genome Atlas RNA-sequencing datasets. A total of 734 lncRNAs were detected to be aberrantly expressed, among which three novel lncRNAs including AC079630.2, CRNDE and CTD-2171N6.1 were markedly related to the progression and survival of PTC. Furthermore, the aberrant expression of these lncRNAs could be verified by other cohorts from gene expression omnibus (GEO) as detected by microarrays. Next, we established a three-lncRNA signature and divided the PTC patients into two subgroups of high- and low-risk. Interestingly, patients with high-risk tended to gain obviously poorer outcome. Most importantly, this three-lncRNA signature was an independent biomarker to predict the patient survival of PTC. The accurate molecular roles of these three lncRNAs remains unclarified and warrants further characterization, but our current data propose that they might play pivotal roles in PTC tumorigenesis and more importantly, these novel lncRNAs are closely related to patients’ survival. These discoveries will have far-reaching consequences with respect to molecular prediction of PTC.