Expression Profiling of a Human Thyroid Cell Line Stably Expressing the BRAFV600E Mutation

BRAFV600E Mutation Enhances Estrogen-Induced Metastatic Potential of Thyroid Cancer by Regulating the Expression of Estrogen Receptors

BACKGRUOUND

Cross-talk between mitogen-activated protein kinase and estrogen has been reported; however, the role of BRAFV600E in the estrogen responsiveness of thyroid cancer is unknown. We elucidated the effect of BRAFV600E on the estrogen-induced increase in metastatic potential in thyroid cancer.

METHODS

Using a pair of cell lines, human thyroid cell lines which harbor wild type BRAF gene (Nthy/WT) and Nthy/BRAFV600E (Nthy/V600E), the expression of estrogen receptors (ERs) and estrogen-induced metastatic phenotypes were evaluated. Susceptibility to ERα- and ERβ-selective agents was evaluated to confirm differential ER expression. ESR expression was analyzed according to BRAFV600E status and age (≤50 years vs. >50 years) using The Cancer Genome Atlas (TCGA) data.

RESULTS

Estradiol increased the ERα/ERβ expression ratio in Nthy/V600E, whereas the decreased ERα/ERβ expression ratio was found in Nthy/WT. BRAFV600E-mutated cell lines showed a higher E2-induced increase in metastatic potential, including migration, invasion, and anchorage-independent growth compared with Nthy/WT. An ERα antagonist significantly inhibited migration in Nthy/V600E cells, whereas an ERβ agonist was more effective in Nthy/WT. In the BRAFV600E group, ESR1/ESR2 ratio was significantly higher in younger age group (≤50 years) compared with older age group (>50 years) by TCGA data analysis.

CONCLUSION

Our data show that BRAFV600E mutation plays a crucial role in the estrogen responsiveness of thyroid cancer by regulating ER expression. Therefore, BRAFV600E might be used as a biomarker when deciding future hormone therapies based on estrogen signaling in thyroid cancer patients.

Induction of BRAF^V600E Mutation in Thyroid Cells Leads to Frequent Hypermethylation

Objectives.

The BRAF^V600E mutation is a major driver mutation in papillary thyroid cancer. The aim of this study was to elucidate the correlation between DNA methylation and gene expression changes induced by the BRAF^V600E mutation in thyroid cells.

Methods.

We used Nthy/BRAF cell lines generated by transfection of Nthy/ori cells with the wild-type BRAF gene (Nthy/WT cells) and the V600E mutant-type BRAF gene (Nthy/V600E cells). We performed gene expression microarray and DNA methylation array analyses for Nthy/WT and Nthy/V600E cells. Two types of array data were integrated to identify inverse correlations between methylation and gene expression. The results were verified in silico using data from The Cancer Genome Atlas (TCGA) and in vivo through pyrosequencing and quantitative real-time polymerase chain reaction (qRT-PCR).

Results.

In the Nthy/V600E cells, 199,821 probes were significantly hypermethylated, and 697 genes showed a “hypermethylation-downregulation” pattern in Nthy/V600E. Tumor suppressor genes and apoptosis-related genes were included. In total, 66,446 probes were significantly hypomethylated, and 227 genes showed a “hypomethylation-upregulation” pattern in Nthy/V600E cells. Protooncogenes and developmental protein-coding genes were included. In the TCGA analysis, 491/697 (70.44%) genes showed a hypermethylation-downregulation pattern, and 153/227 (67.40%) genes showed a hypomethylation-upregulation pattern. Ten selected genes showed a similar methylation-gene expression pattern in pyrosequencing and qRT-PCR.

Conclusion.

Induction of the BRAF^V600E mutation in thyroid cells led to frequent hypermethylation. Anticancer genes, such as those involved in tumor suppression or apoptosis, were downregulated by upstream hypermethylation, whereas carcinogenic genes, such as protooncogenes, were upregulated by hypomethylation. Our results suggest that the BRAF^V600E mutation in thyroid cells modulates DNA methylation and results in cancer-related gene expression.

Comparative Proteomic Profiling of Ectosomes Derived from Thyroid Carcinoma and Normal Thyroid Cells Uncovers Multiple Proteins with Functional Implications in Cancer

Proteins carried by tumor-derived ectosomes play an important role in cancer progression, and are considered promising diagnostic markers. In the present study, a shotgun nanoLC–MS/MS proteomic approach was applied to profile and compare the protein content of ectosomes released in vitro by normal human thyroid follicular epithelial Nthy-ori 3-1 cells and human anaplastic thyroid carcinoma (TC) 8305C cells. Additionally, the pro-migratory and pro-proliferative effects of Nthy-ori 3-1- and 8305C-derived ectosomes exerted on the recipient cells were assessed in wound closure and Alamar Blue assays. A total of 919 proteins were identified in all replicates of 8305C-derived ectosomes, while Nthy-ori 3-1-derived ectosomes contained a significantly lower number of 420 identified proteins. Qualitative analysis revealed 568 proteins present uniquely in 8305C-derived ectosomes, suggesting their applicability in TC diagnosis and management. In addition, 8305C-derived ectosomes were able to increase the proliferation and motility rates of the recipient cells, likely due to the ectosomal transfer of the identified cancer-promoting molecules. Our description of ectosome protein content and its related functions provides the first insight into the role of ectosomes in TC development and progression. The results also indicate the applicability of some of these ectosomal proteins for further investigation regarding their potential as circulating TC biomarkers.

Lactate Dehydrogenase A as a Potential New Biomarker for Thyroid Cancer

BACKGROUND

Several cancers show increased levels of lactate dehydrogenase A (LDHA), which are associated with cancer progression. However, it remains unclear whether LDHA levels are associated with papillary thyroid cancer (PTC) aggressiveness or with the presence of the PTC prognostic marker, the BRAFV600E mutation. This study aimed to evaluate the potential of LDHA as a PTC prognostic marker.

METHODS

LDHA expression was examined in 83 PTC tissue specimens by immunohistochemistry. Human thyroid cell lines were genetically manipulated to overexpress BRAFV600E or were treated with a BRAF-specific short hairpin RNA (shBRAF), whose effects on LDHA expression were evaluated by Western blotting. Data from 465 PTC patients were obtained from The Cancer Genome Atlas (TCGA) database and analyzed to validate the in vitro results.

RESULTS

LDHA was aberrantly overexpressed in PTC. Intense immunostaining for LDHA was observed in PTC specimens carrying mutated BRAF, whereas the intensity was less in wild-type BRAF samples. Overexpression of BRAFV600E resulted in LDHA upregulation, whereas treatment with shBRAF downregulated LDHA in human thyroid cell lines. Furthermore, LDHA mRNA expression was significantly elevated and associated with BRAFV600E expression in thyroid cancer tissues from TCGA database. Additionally, LDHA overexpression was found to be correlated with aggressive clinical features of PTC, such as lymph node metastases and advanced tumor stages.

CONCLUSION

LDHA overexpression is associated with the BRAFV600E mutation and an aggressive PTC behavior. Therefore, LDHA may serve as a biomarker and therapeutic target in PTC.

PTPN11 Knockdown Prevents Changes in the Expression of Genes Controlling Cell Cycle, Chemotherapy Resistance, and Oncogene-Induced Senescence in Human Thyroid Cells Overexpressing BRAF V600E Oncogenic Protein

The MAPK (RAS/BRAF/MEK/ERK) signaling pathway is a kinase cascade involved in the regulation of cell proliferation, differentiation, and survival in response to external stimuli. The V600E mutation in the BRAF gene has been detected in various tumors, resulting in a 500-fold increase in BRAF kinase activity. However, monotherapy with selective BRAF V600E inhibitors often leads to reactivation of MAPK signaling cascade and emergence of drug resistance. Therefore, new targets are being developed for the inhibition of components of the aberrantly activated cascade. It was recently discovered that resistance to BRAF V600E inhibitors may be associated with the activity of the tyrosine phosphatase SHP-2 encoded by the PTPN11 gene. In this paper, we analyzed transcriptional effects of PTPN11 gene knockdown and selective suppression of BRAF V600E in a model of thyroid follicular epithelium. We found that the siRNA-mediated knockdown of PTPN11 after vemurafenib treatment prevented an increase in the expression CCNA1 and NOTCH4 genes involved in the formation of drug resistance of tumors. On the other hand, downregulation of PTPN11 expression blocked the transcriptional activation of genes (p21, p15, p16, RB1, and IGFBP7) involved in cell cycle regulation and oncogene-induced senescence in response to BRAF V600E expression. Therefore, it can be assumed that SHP-2 participates not only in emergence of drug resistance in cancer cells, but also in oncogene-induced cell senescence.

BRAFV600E Transduction of an SV40-Immortalized Normal Human Thyroid Cell Line Induces Dedifferentiated Thyroid Carcinogenesis in a Mouse Xenograft Model

Background: Despite active studies of the clinical importance of BRAF^V600E, suitable research models to investigate the role of this mutation in the etiopathogenesis of human thyroid cancers are limited. Thus, we generated cell lines by transducing the simian virus (SV)-40 immortalized human thyroid cell line Nthy-ori 3-1 (Nthy) with lentiviral vectors expressing either BRAF^WT (Nthy/WT) or BRAF^V600E. Nthy/WT and Nthy/V600E cells were then xenografted into mice to evaluate the carcinogenic role of BRAF^V600E. Methods: Each cell line was subcutaneously injected into NOD.Cg-Prkdc^scid Il2rg^tm1Wjl/SzJ mice, and a pathological analysis was performed. The effects of the mutation were further verified by using a BRAF^V600E-selective inhibitor (PLX-4032, vemurafenib). The transcriptome was analyzed by RNA sequencing and compared with data from The Cancer Cell Line Encyclopedia and Gene Expression Omnibus. Results: While Nthy/WT was not tumorigenic in vivo, Nthy/V600E formed tumors reaching 2784.343 mm³ in 4 weeks, on average. A pathological analysis indicated that Nthy/V600E tumors were dedifferentiated thyroid cancer. We found metastases in the lung, liver, and relevant lymph nodes. A transcriptomic analysis revealed 5512 differentially expressed genes (DEGs) between the mutant and wild-type cell lines, and more DEGs were shared with anaplastic thyroid cancer than with papillary thyroid cancer. BRAF^V600E activated the cell cycle mainly by regulating G1/S phases. PLX-4032 treatment significantly inhibited tumor growth and metastasis. Conclusions: Our data show that BRAF^V600E plays a pivotal role in the carcinogenic transformation of an SV40-transfected immortalized normal human thyroid cell line. This xenograft model is expected to contribute to studies of the etiopathogenesis and treatment of highly malignant thyroid cancers.

Caveolin-1 Regulates CCL5 and PPARγ Expression in Nthy-ori 3-1 Cells: Possible Involvement of Caveolin-1 and CCL5 in the Pathogenesis of Hashimoto's Thyroiditis

BACKGROUND

Hashimoto's thyroiditis (HT) is characterized by lymphocytic infiltration of the thyroid parenchyma, which ultimately leads to tissue destruction and loss of function. Caveolin-1 (Cav-1) is an essential structural constituent of lipid rafts in the plasma membrane of cells and is reported to be significantly reduced in thyrocytes from HT patients. However, the mechanism of Cav-1 involvement in HT pathogenesis is still largely unclear.

METHODS

Cav-1 expression in thyroid tissues from HT patients and euthyroid nodular goiter tissues was detected by immunohistochemistry staining. Cav-1 knockdown and overexpression were constructed by lentiviral transfection in the human thyroid follicular epithelial cell (TFC) line of Nthy-ori 3-1. The mRNA expression levels of chemokines in TFCs were determined by quantitative real-time PCR (qPCR). Cav-1 and peroxisome proliferator-activated receptor gamma (PPARγ) levels were analysed by qPCR and Western blot analysis. The migration ability of peripheral blood mononuclear cells (PBMCs) was detected by the Transwell assay.

RESULTS

In this study, Cav-1 and PPARγ expression was reduced in the thyroid tissues from HT patients. In vitro experiments showed that the expressions of chemokine (C-C motif) ligand 5 (CCL5) and migration of PBMCs were markedly increased, while the level of PPARγ was significantly decreased after the lentivirus-mediated knockdown of Cav-1 in Nthy-ori 3-1 cells. Interestingly, pioglitazone, a PPARγ agonist, not only upregulated PPARγ and Cav-1 proteins significantly, but also effectively reversed the Cav-1-knockdown-induced upregulation of CCL5 in Nthy-ori 3-1 cells and reduced the infiltration of lymphocytes.

CONCLUSION

The inhibition of Cav-1 upregulated the CCL5 expression and downregulated the PPARγ expression in TFC while pioglitazone, a PPARγ agonist, reversed the detrimental consequence. This outcome might be a potential target for the treatment of lymphocyte infiltration into the thyroid gland and HT development.

Constitutive Expression of NRAS with Q61R Driver Mutation Activates Processes of Epithelial-Mesenchymal Transition and Leads to Substantial Transcriptome Change of Nthy-ori 3-1 Thyroid Epithelial Cells

The Q61R mutation of the NRAS gene is one of the most frequent driver mutations of thyroid cancer. Tumors with this mutation are characterized by invasion into blood vessels and formation of distant metastases. To study the role of this mutation in the growth of thyroid cancer, we developed a model system on the basis of thyroid epithelial cell line Nthy-ori 3-1 transduced by a lentiviral vector containing the NRAS gene with the Q61R mutation. It was found that the expression of NRAS(Q61R) in thyroid epithelial cells has a profound influence on groups of genes involved in the formation of intercellular contacts, as well as in processes of epithelial-mesenchymal transition and cell invasion. The alteration in the expression of these genes affects the phenotype of the model cells, which acquire traits of mesenchymal cells and demonstrate increased ability for survival and growth without attachment to the substrate. The key regulators of these processes are transcription factors belonging to families SNAIL, ZEB, and TWIST, and in different types of tumors the contribution of each individual factor can vary greatly. In our model system, phenotype change correlates with an increase in the expression of SNAIL2 and TWIST2 factors, which indicates their possible role in regulating invasive growth of thyroid cancer with the mutation of NRAS(Q61R).

Interaction of BRAF-induced ETS factors with mutant TERT promoter in papillary thyroid cancer

Synergistic effects of BRAFV600E and TERT promoter mutations on the poor clinical outcomes in papillary thyroid cancer (PTC) have been demonstrated. The potential mechanism of this phenomenon has been proposed: MAPK pathway activation by the BRAFV600E mutation may upregulate E-twenty six (ETS) transcription factors, increasing TERT expression by binding to the ETS-binding site generated by the TERT promoter mutation; however, it has not yet been fully proven. This article provides transcriptomic insights into the interaction between BRAFV600E and TERT promoter mutations mediated by ETS factors in PTC. RNA sequencing data on 266 PTCs from The Cancer Genome Atlas and 65 PTCs from our institute were analyzed for gene expression changes and related molecular pathways, and the results of transcriptomic analyses were validated by in vitro experiments. TERT mRNA expression was increased by the coexistence of BRAFV600E and TERT promoter mutations (fold change, 16.17; q-value = 7.35 × 10-12 vs no mutation). In the ETS family of transcription factors, ETV1, ETV4 and ETV5 were upregulated by the BRAFV600E/MAPK pathway activation. These BRAFV600E-induced ETS factors selectively bound to the mutant TERT promoter. The molecular pathways activated by BRAFV600E were further augmented by adding the TERT promoter mutation, and the pathways related to immune responses or adhesion molecules were upregulated by TERT expression. The mechanism of the synergistic effect between BRAFV600E and TERT promoter mutations on cancer invasiveness and progression in PTC may be explained by increased TERT expression, which may result from the BRAF-induced upregulation of several ETS transcription factors.

Methylation dependent down-regulation of G0S2 leads to suppression of invasion and improved prognosis of IDH1-mutant glioma

Isocitrate dehydrogenase (IDH) mutations are a prognostic factor in diffuse glioma. However, the mechanism by which these mutations improve prognosis are not clear. In a subset of IDH-mutant glioma, remodeling of the methylome results in the glioma-CpG island methylator phenotype (G-CIMP) and transcriptional reorganization. In this study, we focus on G0/G1 switch 2 (G0S2), which is highly downregulated in G-CIMP glioma. We found that G0S2 expression tended to increase as the WHO grade increased, and G0S2 knockdown inhibited glioma invasion. Additionally, we revealed that the overexpression of the DNA demethylase Ten-eleven translocation 2 (TET2) in IDH1-plasmid transfected glioblastoma multiforme (GBM) cells restored G0S2 expression. These results indicate that G0S2 is epigenetically silenced in IDH1-mutant glioma. In addition, the stereotactic delivery of glioma cells with decreased G0S2 expression in the mouse brain resulted in prolonged survival. The Cancer Genome Atlas (TCGA) analysis also indicated that survival is longer in the lower G0S2 expression group than in the higher G0S2 expression group. Moreover, G0S2 expression was higher in recurrent tumor specimens than at the initial diagnosis in the same patient. These results provide one explanation for the improved survival in IDH1-mutant glioma as well as a new epigenetic target for glioma treatment.

Caveolin-1 regulates autophagy activity in thyroid follicular cells and is involved in Hashimoto's thyroiditis disease

Hashimoto's thyroiditis (HT) is considered a T helper-type 1 (Th1) cytokine-dominant autoimmune thyroid disease. Caveolin-1 (Cav-1), a part of the thyroxisome multiprotein complex, is localized at the apical pole of thyrocytes and is indispensable for synthesis of thyroid hormones and modulation of oxidative stress in order to avoid cell damage and apoptosis. Reduced autophagy induces thyroid follicular cells (TFC) apoptosis by activating reactive oxygen species (ROS) in HT patients. Nevertheless, whether Cav-1 has roles in the regulation of autophagy remains largely unclear. In this study, we examined Th1 cytokines and Cav-1 expression in HT thyroid tissues, determined the effects of interleukin-1beta (IL-1β) and interferon-gamma (IFN-γ) on Cav-1 and autophagy activity in TFC, and investigated the association between Cav-1 and autophagy activity in vitro. Our results indicate that higher levels of IL-1β and IFN-γ and lower levels of Cav-1 were expressed in thyroid tissues of HT patients than in those of normal controls. Cav-1 mRNA and protein levels were significantly decreased in TFC exposed to IL-1β and IFN-γ, accompanied by decreased expression of autophagy-related protein LC3B-II. Interestingly, small interfering RNA (siRNA)-mediated Cav-1 knockdown in TFC reduced LC3B-II protein expression. Taken together, these results suggest that lack of Cav-1 expression inhibited autophagy activity in TFC exposed to Th1 cytokines (IL-1β and IFN-γ), which might be a novel pathogenetic mechanism of HT.

DNA microarray-based resonance light scattering assay for multiplexed detection of DNA mutation in papillary thyroid cancer

A DNA microarray-based resonance light scattering assay has been developed for multiplexed detection of papillary thyroid carcinoma (PTC) related genic mutation.